Search NASA⌕ Search

SEARCH · Search NASA

Results for “node nonlinearity”

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

Fully-Coupled Fluid-Structure Interaction Simulations of a Supersonic Parachute

A validated computational fluid-structure interaction method for simulating the complex interaction between the large deformation of very thin, highly deformable structures and compressible flows is extended to consider large-scale problems in supersonic flows using parallel computing. The coupled fluid-structure interaction system is solved in a partitioned, or weakly-coupled, manner. The foundations of the applied fluid-structure interaction method are a higher-order, block-structured Cartesian, sharp immersed boundary method for the compressible Navier-Stokes equations and a computational structural dynamics solver employing a geometrically nonlinear 3-node shell element based on the mixed interpolation of tensorial components formulation. The method is applied to large deformation fluid-structure interaction validation cases before being applied to the inflation of a supersonic parachute in the upper Martian atmosphere where the goal is to demonstrate the capabilities of the solver when considering large-scale problems in supersonic flows.

Boustani, Jonathan↗

A three-dimensional nonlinear Timoshenko beam based on the core-congruential formulation

A three-dimensional, geometrically nonlinear two-node Timoshenkoo beam element based on the total Larangrian description is derived. The element behavior is assumed to be linear elastic, but no restrictions are placed on magnitude of finite rotations. The resulting element has twelve degrees of freedom: six translational components and six rotational-vector components. The formulation uses the Green-Lagrange strains and second Piola-Kirchhoff stresses as energy-conjugate variables and accounts for the bending-stretching and bending-torsional coupling effects without special provisions. The core-congruential formulation (CCF) is used to derived the discrete equations in a staged manner. Core equations involving the internal force vector and tangent stiffness matrix are developed at the particle level. A sequence of matrix transformations carries these equations to beam cross-sections and finally to the element nodal degrees of freedom. The choice of finite rotation measure is made in the next-to-last transformation stage, and the choice of over-the-element interpolation in the last one. The tangent stiffness matrix is found to retain symmetry if the rotational vector is chosen to measure finite rotations. An extensive set of numerical examples is presented to test and validate the present element.

Crivelli, Luis A.↗

Development of computer program NAS3D using Vector processing for geometric nonlinear analysis of structures

An algorithm for vectorized computation of stiffness matrices of an 8 noded isoparametric hexahedron element for geometric nonlinear analysis was developed. This was used in conjunction with the earlier 2-D program GAMNAS to develop the new program NAS3D for geometric nonlinear analysis. A conventional, modified Newton-Raphson process is used for the nonlinear analysis. New schemes for the computation of stiffness and strain energy release rates is presented. The organization the program is explained and some results on four sample problems are given. The study of CPU times showed that savings by a factor of 11 to 13 were achieved when vectorized computation was used for the stiffness instead of the conventional scalar one. Finally, the scheme of inputting data is explained.

Mangalgiri, P. D.↗

Low-Thrust Trajectory Optimization with Simplified SQP Algorithm

The problem of low-thrust trajectory optimization in highly perturbed dynamics is a stressing case for many optimization tools. Highly nonlinear dynamics and continuous thrust are each, separately, non-trivial problems in the field of optimal control, and when combined, the problem is even more difficult. This paper de-scribes a fast, robust method to design a trajectory in the CRTBP (circular restricted three body problem), beginning with no or very little knowledge of the system. The approach is inspired by the SQP (sequential quadratic programming) algorithm, in which a general nonlinear programming problem is solved via a sequence of quadratic problems. A few key simplifications make the algorithm presented fast and robust to initial guess: a quadratic cost function, neglecting the line search step when the solution is known to be far away, judicious use of end-point constraints, and mesh refinement on multiple shooting with fixed-step integration.In comparison to the traditional approach of plugging the problem into a “black-box” NLP solver, the methods shown converge even when given no knowledge of the solution at all. It was found that the only piece of information that the user needs to provide is a rough guess for the time of flight, as the transfer time guess will dictate which set of local solutions the algorithm could converge on. This robustness to initial guess is a compelling feature, as three-body orbit transfers are challenging to design with intuition alone. Of course, if a high-quality initial guess is available, the methods shown are still valid.We have shown that endpoints can be efficiently constrained to lie on 3-body repeating orbits, and that time of flight can be optimized as well. When optimizing the endpoints, we must make a trade between converging quickly on sub-optimal endpoints or converging more slowly on end-points that are arbitrarily close to optimal. It is easy for the mission design engineer to adjust this trade based on the problem at hand.The biggest limitation to the algorithm at this point is that multi-revolution transfers (greater than 2 revolutions) do not work nearly as well. This restriction comes in because the relationship between node 1 and node N becomes increasingly nonlinear as the angular distance grows. Trans-fers with more than about 1.5 complete revolutions generally require the line search to improve convergence. Future work includes: Comparison of this algorithm with other established tools; improvements to how multiple-revolution transfers are handled; parallelization of the Jacobian computation; in-creased efficiency for the line search; and optimization of many more trajectories between a variety of 3-body orbits.

Parrish, Nathan L.↗

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

Optimal nodal point distribution for improved accuracy in computational fluid dynamics

In applying finite-difference techniques to flow field problems, the accuracy attained for a fixed number of node points can be improved using unequally-spaced node points. The distribution of these node points is chosen here by minimizing a measure of local truncation error with respect to the parameters which define a transformation between the computational space of equally-spaced node points and the physical space of unequally-spaced node points. The problem then becomes a nonlinear programming problem. Numerical results are presented for two one-dimensional test problems: the Blasius boundary layer problem and the inviscid Burgers' equation.

Pierson, B. L.↗

Modeling of flow systems for implementation under KATE

The modeling of flow systems is a task currently being investigated at Kennedy Space Center in parallel with the development of the KATE artificial intelligence system used for monitoring diagnosis and control. Various aspects of the modeling issues are focussed on with particular emphasis on a water system scheduled for demonstration within the KATE environment in September of this year. LISP procedures were written to solve the continuity equations for three internal pressure nodes using Newton's method for simultaneous nonlinear equations.

Whitlow, Jonathan E.↗

Detailed analysis and test correlation of a stiffened composite wing panel

New finite element analysis methods are examined by application to a complicated composite wing panel from the V-22 rotorcraft. A detailed FEM model with a relatively coarse mesh of 9-node elements was generated, and linear and nonlinear stress analyses, first-ply failure analyses, and buckling analyses were conducted. At low values of applied load, i.e., up to the design ultimate load of the panel, the linear stress analysis accurately predicted the strains and structural response characteristics of the panel.

Davis, D. D., Jr.↗

Three dimensional inelastic finite element analysis of laminated composites

Formulations of the inelastic response of laminated composites to thermal and mechanical loading are used as the basis for development of the computer NALCOM (Nonlinear Analysis of Laminated Composites) computer program which uses a fully three dimensional isoparametric finite element with 24 nodes and 72 degrees of freedom. An incremental solution is performed with nonlinearities introduced as pseudoloads computed for initial strains. Equilibrium iteration may be performed at every step. Elastic and elastic-plastic response of boron/epoxy and graphite/epoxy graphite/epoxy and problems of curing 0/90 sub s Gr/Ep laminates with and without circular holes are analyzed. Mechanical loading of + or - 45sub s Gr/Ep laminates is modeled and symmetry conditions which exist in angle-ply laminates are discussed. Results are compared to experiments and other analytical models when possible. All models are seen to agree reasonably well with experimetnal results for off-axis tensile coupons. The laminate analyses show the three dimensional effects which are present near holes and free corners.

Griffin, O. H., Jr.↗

Post-Buckling Analysis of Curved Honeycomb Sandwich Panels Containing Interfacial Disbonds

A numerical study on the effect of facesheet-core disbonds on the post-buckling response of curved honeycomb sandwich panels is presented herein. This work was conducted as part of the development of a damage tolerance plan for the next-generation Space Launch System heavy lift launch vehicle payload fairing. As such, the study utilized full-scale fairing barrel segments as the structure of interest. The panels were composed of carbon fiber reinforced polymer facesheets and aluminum honeycomb core. The panels were analyzed numerically using the finite element method incorporating geometric nonlinearity. In a predetermined circular region, facesheet and core nodes were detached to simulate a disbond, between the outer mold line facesheet and honeycomb core, induced via low-speed impact. Surface-to-surface contact in the disbonded region was invoked to prevent interpenetration of the facesheet and core elements and obtain realistic stresses in the core. The diameter of this disbonded region was varied and the effect of the size of the disbond on the post-buckling response was observed. Significant changes in the slope of the edge load-deflection response were used to determine the onset of global buckling and corresponding buckling load. Finally, several studies were conducted to determine the sensitivity of the numerical predictions to refinement in the finite element mesh.

elastic buckling↗

High-impact dynamic-response analysis of nonlinear structures

Program predicts expected deformations and stresses in nonlinear simple geometric structures subjected to high-impact loading. Technique is based on node-wise predictor-corrector approach and requires moderate computer storage and run time for most problems. Program extends to include physical and geometrical nonlinearities.

Gupta, K. K.↗

SHARP simulation of discontinuities in highly convective steady flow

For steady multidimesional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

Simple high-accuracy resolution program for convective modelling of discontinuities

For steady multidimensional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

Application of the method of integral relations to unsteady fluid flow problems with shocks

A mixed method using both the method of integral relations and the finite difference technique is developed for the solution of unsteady flow problems. The integral relations method is based upon a chosen interpolating function dependent only on the time domain. The resulting local semidiscrete finite element equations obtained are assembled into a global form. The spatial derivatives at the nodes are replaced by finite difference operators and the discretized nonlinear algebraic system is solved by an iterative scheme. Solutions are obtained for the one-dimensional gasdynamics equation, the one-dimensional wave equation and Burger's model of turbulence. Agreement with other numerical and analytical solutions is excellent for the gasdynamics problem and satisfactory for the Burger equation in cases of small viscous effects.

Wadia, A. R.↗

Core reactivity estimation in space reactors using recurrent dynamic networks

A recurrent multilayer perceptron network topology is used in the identification of nonlinear dynamic systems from only the input/output measurements. The identification is performed in the discrete time domain, with the learning algorithm being a modified form of the back propagation (BP) rule. The recurrent dynamic network (RDN) developed is applied for the total core reactivity prediction of a spacecraft reactor from only neutronic power level measurements. Results indicate that the RDN can reproduce the nonlinear response of the reactor while keeping the number of nodes roughly equal to the relative order of the system. As accuracy requirements are increased, the number of required nodes also increases, however, the order of the RDN necessary to obtain such results is still in the same order of magnitude as the order of the mathematical model of the system. It is believed that use of the recurrent MLP structure with a variety of different learning algorithms may prove useful in utilizing artificial neural networks for recognition, classification, and prediction of dynamic systems.

Parlos, Alexander G.↗

Nonlinear Response of Composite Panels Under Combined Acoustic Excitation and Aerodynamic Pressure

A finite element formulation is presented for the analysis of large deflection response of composite panels subjected to aerodynamic pressure- at supersonic flow and high acoustic excitation. The first-order shear deformation theory is considered for laminated composite plates, and the von Karman nonlinear strain-displacement relations are employed for the analysis of large deflection panel response. The first-order piston theory aerodynamics and the simulated Gaussian white noise are employed for the aerodynamic and acoustic loads, respectively. The nonlinear equations of motion for an arbitrarily laminated composite panel subjected to a combined aerodynamic and acoustic pressures are formulated first in structure node degrees-of-freedom. The system equations are then transformed and reduced to a set of coupled nonlinear equations in modal coordinates. Modal participation is defined and the in-vacuo modes to be retained in the analysis are based on the modal participation values. Numerical results include root mean square values of maximum deflections, deflection and strain response time histories, probability distributions, and power spectrum densities. Results showed that combined acoustic and aerodynamic loads have to be considered for panel analysis and design at high dynamic pressure values.

Abdel-Motagaly, K.↗

Nonplanar interface morphologies during unidirectional solidification of a binary alloy. II - Three-dimensional computations

A finite difference method is used to obtain three-dimensional steady-state solutions for nonplanar interface morphologies in order to study the situation of equal thermal properties in the crystal and melt with negligible latent heat release. Stable steady-state solutions corresponding to two-dimensional bands and three-dimensional hexagonal nodes, as well as to rectangular interface planiforms, are found using a model of an aluminum-chromium alloy with a distribution coefficient of greater than one. Hexagonal nodes are predicted near the onset of instability, in agreement with weakly nonlinear theory.

Mcfadden, G. B.↗

Analytical and computational models of shells; Proceedings of the Symposium, ASME Winter Annual Meeting, San Francisco, CA, Dec. 10-15, 1989

Topics presented include asymptotic analysis and computation for shells, the edge effects in the Reissner-Mindlin plate theory, nonlinear problems of geometrically exact shell theories, and developments in variational methods for high performance plate and shell elements. Also presented are an assumed strain solid element model for geometrically nonlinear shell analysis, shell finite elements with six degrees of freedom per node, hierarchic plate and shell models based on p-extension, and a simple shell element formulation for large-scale elastoplastic analysis. Also discussed are the assessment of computational models for multilayered composite cylinders, shell models for impact analysis, analysis of shell structures subjected to contact-impacts, and the application of shell theory to cardiac mechanics.

Noor, Ahmed K.↗