SEARCH · Search NASA

Results for “Boundary finite element”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Boundary Layer Transition Induced by Distributed Roughness Array

The effects of a finite, spanwise-periodic array of cylindrical roughness elements on boundary layer transition over a NACA 0012 airfoil are investigated at a chord-based Reynolds number of 1.44×105 by using hotwire anemometry and infrared thermography. Both the number and the spanwise spacing of roughness elements in the array are varied in order to study their effect on the wake flow topology. Spanwise interaction between the roughness elements has an effect on the connection and the merging of neighbouring low-speed regions, which results in the formation of merged low-speed blobs (MLSs) that modify the spatial distribution and the amplitudes of the velocity streaks. When the spanwise distance between adjacent roughness elements equals 1.5 times the cylinder diameter, the transition location moves rapidly upstream. In this case, the two neighbouring low-speed regions overlap with each other in the near wake of the roughness, leading to the maximum growth in the velocity streak amplitude and the velocity fluctuations. The number of roughness elements affects the total number of MLSs within the boundary layer. For a single MLS behind a pair of cylinders, the Kelvin-Helmholtz instability dominates the growth of velocity fluctuations around the three-dimensional shear layers. When three cylinders are placed in the array, two MLSs appear in the near wake, which coalesce in to one low-speed blob downstream before the onset of transition, revealing the importance of Kelvin-Helmholtz instability.

Qingqing Ye↗

NASA NTRS1992

Finite element nonlinear panel flutter with arbitrary temperatures in supersonic flow

A finite element frequency domain method for predicting nonlinear flutter response of panels with temperature effects is presented. By using the principle of virtual work, the element nonlinear stiffness formulation for a panel under a combined thermal and aerodynamic loads is derived on the bases of von Karman's large deflection plate theory, the first-order piston theory aerodynamics and the quasi-steady thermal stress theory. The system equations of motion can be mathematically separated into two sets of equations and then solved in sequence. The first set of equations yields the panel thermal-aerodynamic equilibrium and the second set of equations of motion leads to the flutter limit-cycle oscillations. Stability and flutter boundaries can also be obtained from the two sets of system equations. Finite element large amplitude limit-cycle flutter results at different uniform temperatures are obtained for a simply supported square panel and are compared with existing Galerkin/time integration and other finite element solutions. Effects of nonuniform temperature distributions, panel length-to-width ratios, and boundary conditions on flutter responses of rectangular and triangular panels are presented.

Xue, David Y.↗

NASA NTRS1989

Analytical modeling of structure-soil systems for lunar bases

The study of the behavior of granular materials in a reduced gravity environment and under low effective stresses became a subject of great interest in the mid 1960's when NASA's Surveyor missions to the Moon began the first extraterrestrial investigation and it was found that Lunar soils exhibited properties quite unlike those on Earth. This subject gained interest during the years of the Apollo missions and more recently due to NASA's plans for future exploration and colonization of Moon and Mars. It has since been clear that a good understanding of the mechanical properties of granular materials under reduced gravity and at low effective stress levels is of paramount importance for the design and construction of surface and buried structures on these bodies. In order to achieve such an understanding it is desirable to develop a set of constitutive equations that describes the response of such materials as they are subjected to tractions and displacements. This presentation examines issues associated with conducting experiments on highly nonlinear granular materials under high and low effective stresses. The friction and dilatancy properties which affect the behavior of granular soils with low cohesion values are assessed. In order to simulate the highly nonlinear strength and stress-strain behavior of soils at low as well as high effective stresses, a versatile isotropic, pressure sensitive, third stress invariant dependent, cone-cap elasto-plastic constitutive model was proposed. The integration of the constitutive relations is performed via a fully implicit Backward Euler technique known as the Closest Point Projection Method. The model was implemented into a finite element code in order to study nonlinear boundary value problems associated with homogeneous as well as nonhomogeneous deformations at low as well as high effective stresses. The effect of gravity (self-weight) on the stress-strain-strength response of these materials is evaluated. The calibration of the model is performed via three techniques: (1) physical identification, (2) optimized calibration at the constitutive level, and (3) optimized calibration at the finite element level (Inverse Identification). Activities are summarized in graphic and outline form.

Macari-Pasqualino, Jose Emir↗

NASA NTRS1996

Finite Element Analysis of Geodesically Stiffened Cylindrical Composite Shells Using a Layerwise Theory

Layerwise finite element analyses of geodesically stiffened cylindrical shells are presented. The layerwise laminate theory of Reddy (LWTR) is developed and adapted to circular cylindrical shells. The Ritz variational method is used to develop an analytical approach for studying the buckling of simply supported geodesically stiffened shells with discrete stiffeners. This method utilizes a Lagrange multiplier technique to attach the stiffeners to the shell. The development of the layerwise shells couples a one-dimensional finite element through the thickness with a Navier solution that satisfies the boundary conditions. The buckling results from the Ritz discrete analytical method are compared with smeared buckling results and with NASA Testbed finite element results. The development of layerwise shell and beam finite elements is presented and these elements are used to perform the displacement field, stress, and first-ply failure analyses. The layerwise shell elements are used to model the shell skin and the layerwise beam elements are used to model the stiffeners. This arrangement allows the beam stiffeners to be assembled directly into the global stiffness matrix. A series of analytical studies are made to compare the response of geodesically stiffened shells as a function of loading, shell geometry, shell radii, shell laminate thickness, stiffener height, and geometric nonlinearity. Comparisons of the structural response of geodesically stiffened shells, axial and ring stiffened shells, and unstiffened shells are provided. In addition, interlaminar stress results near the stiffener intersection are presented. First-ply failure analyses for geodesically stiffened shells utilizing the Tsai-Wu failure criterion are presented for a few selected cases.

Gerhard, Craig Steven↗

NASA NTRS1972

NASTRAN thermal analyzer: A general purpose finite element heat transfer computer program

The program not only can render temperature distributions in solids subjected to various thermal boundary conditions, including effects of diffuse-gray thermal radiation, but is fully compatible in capacity and in the finite-element model representation with that of its structural counterpart in the NASTRAN system. The development history of the finite-element approach for determining temperatures is summarized. The scope of analysis capability, program structure, features, and limitations are given with the objective of providing NASTRAN users with an overall veiw of the NASTRAN thermal analyzer.

Lee, H.↗

NASA NTRS1995

Adaptive Finite Element Methods for Continuum Damage Modeling

The paper presents an application of adaptive finite element methods to the modeling of low-cycle continuum damage and life prediction of high-temperature components. The major objective is to provide automated and accurate modeling of damaged zones through adaptive mesh refinement and adaptive time-stepping methods. The damage modeling methodology is implemented in an usual way by embedding damage evolution in the transient nonlinear solution of elasto-viscoplastic deformation problems. This nonlinear boundary-value problem is discretized by adaptive finite element methods. The automated h-adaptive mesh refinements are driven by error indicators, based on selected principal variables in the problem (stresses, non-elastic strains, damage, etc.). In the time domain, adaptive time-stepping is used, combined with a predictor-corrector time marching algorithm. The time selection is controlled by required time accuracy. In order to take into account strong temperature dependency of material parameters, the nonlinear structural solution a coupled with thermal analyses (one-way coupling). Several test examples illustrate the importance and benefits of adaptive mesh refinements in accurate prediction of damage levels and failure time.

Min, J. B.↗

NASA NTRS1977

Calculation methods for compressible turbulent boundary layers, 1976

Equations and closure methods for compressible turbulent boundary layers are discussed. Flow phenomena peculiar to calculation of these boundary layers were considered, along with calculations of three dimensional compressible turbulent boundary layers. Procedures for ascertaining nonsimilar two and three dimensional compressible turbulent boundary layers were appended, including finite difference, finite element, and mass-weighted residual methods.

Bushnell, D. M.↗

NASA NTRS1985

A numerical algorithm for endochronic plasticity and comparison with experiment

A numerical algorithm based on the finite element method of analysis of the boundary value problem in a continuum is presented, in the case where the plastic response of the material is given in the context of endochronic plasticity. The relevant constitutive equation is expressed in incremental form and plastic effects are accounted for by the method of an induced pseudo-force in the matrix equations. The results of the analysis are compared with observed values in the case of a plate with two symmetric notches and loaded longitudinally in its own plane. The agreement between theory and experiment is excellent.

Valanis, K. C.↗

NASA NTRS1986

Combined bearing and bypass loading on a graphite/epoxy laminate

A combined experimental and analytical study was conducted to determine the behavior of a graphite/epoxy laminate subjected to combined bearing and bypass loading. Single-fastener quasi-isotropic specimens were loaded at various bearing-bypass ratios until damage was produced at the fastener hole. Damage-onset strengths and damage modes were then analyzed using local hole-boundary stresses calculated by a finite-element analysis. The tension data showed the expected linear interaction for combined bearing and bypass loading with damage developing in the net-section tension mode. However, the compression bearing-bypass strength showed an unexpected interaction involving the bearing mode.

Crews, J. H., Jr.↗

NASA NTRS1987

Calculation of stresses and forces between the skin and stiffener in composite panels

A methodology for computing the stresses at the interface between the skin and stiffener in stiffened composite panels is described. The methodology is based on finite-element analyses and on elasticity solution. The finite-element analyses are standard, while the elasticity solution is based on eigenvalue expansions of the stress function. The eigenvalue expansions are assumed to be valid in the local region where the stiffener flange terminates. The local elasticity solution is coupled to the global finite-element analysis using collocation on the boundary of the local region. Accuracy and convergence of the method are discussed and several examples of its utility are presented.

Hyer, M. W.↗

NASA NTRS1987

On a finite element CFD algorithm for compressible, viscous and turbulent aerodynamic flows

This paper develops and analyses individual construction aspects of an efficient and accurate finite element algorithm for prediction of viscous and turbulent flow fields of impact in aerodynamics. The theoretical construction employs a Taylor weak statement (TWS) for coincident embedding of stability mechanisms within a classic Galerkin finite element formulation of semidiscrete approximation error orthogonalization. A wide variety of the stabilizing mechanisms of independently derived CFD algorithms are contained within the TWS theory. An implicit construction that meets the requirement of efficient convergence to steady state is developed. The theoretical asymptotic error estimates of the TWS finite element algorithm for supersonic and viscous boundary layer flows are verified. Application to a three-dimensional turbulent flow is cited.

Baker, A. J.↗

NASA NTRS1990

Free-vibration analysis of three-dimensional solids by BEM

This paper discusses the calculation of natural frequencies and mode shapes of three-dimensional solids, using the boundary element method. The metod developed is based on the use of particular integrals of the elastic equilibrium equations, and employs only real arithmetic. It represents an extension to three dimensions of work previously carried out for two-dimensional problems by Ahmad and Banerjee (1986) and Nardini and Brebbia (1982). The method has been incorporated, for multiregion analysis, in the BEST3D (Boundary Element Solution Technology, Three-Dimensional) computer program. Results of the boundary element calculations are compared with analytical, finite element, and experimental results.

Wilson, R. B.↗

NASA NTRS1989

Elevated temperature crack growth

Alloy 718 crack growth experiments were conducted to assess the ability of the selected path-independent (P-I) integrals to describe the elevated temperature crack growth behavior. These tests were performed on single edge notch (SEN) specimens under displacement control with multiple extensometers to monitor the specimen and crack mouth opening displacement (CMOD). The displacements in these tests were sufficiently high to induce bulk cyclic inelastic deformation of the specimen. Under these conditions, the linear elastic fracture mechanics (LEFM) parameter K does not correlate the crack growth data. The experimentally measured displacement gradients at the end of specimen gage length were used as the boundary conditions in elastic-plastic finite element method (FEM) analyses. These analyses were performed with a node release approach using CYANIDE, a GEAE FEM code, which included a gap element which is capable of efficiently simulating crack closure. Excellent correlation was obtained between the experimentally measured and predicted variation of stress and CMOD with crack length and the stress-CMOD loops for Alloy 718 tests conducted at 538 C. This confirmed the accuracy of the FEM crack growth simulation approach. The experimentally measured crack growth rate data correlated well the selected P-I integrals. These investigations have produced significant progress in developing P-I integrals as non-linear fracture mechanics parameters. The results suggest that this methodology has the potential of accurately describing elevated temperature crack growth behavior under the combined influence of thermal cycling and bulk elastic-inelastic deformation states.

Kim, K. S.↗

NASA NTRS1992

Thermomechanical buckling and postbuckling of multilayered composite panels

A study is made of the thermomechanical buckling and postbuckling responses of flat unstiffened composite panels. The panels are subjected to combined temperature change and applied edge displacement. The analysis is based on a first-order shear deformation, von Karman type nonlinear plate theory. A mixed formulation is used with the fundamental unknowns consisting of the generalized displacements and the stress resultants of the plate. An efficient multiple-parameter reduction method is used in conjunction with mixed finite element models, for determining the stability boundary and postbuckling response. The reduction method is also used for evaluating the sensitivity coefficients which measure the sensitivity of the buckling and postbuckling responses to variations in the different lamination and material parameters of the panel. Numerical results are presented showing the effects of variations in the laminate stacking sequence, fiber orientation, number of layers and aspect ratio of the panels on their thermomechanical buckling and postbuckling responses and their sensitivity coefficients.

Noor, Ahmed K.↗

NASA NTRS1992

Postbuckling of multilayered composite plates subjected to combined axial and thermal loads

A study is made of the postbuckling response of composite plates subjected to combined axial and thermal loadings. The analysis is based on a first-order shear deformation, von Karman type of nonlinear plate theory. A mixed formulation is used with the fundamental unknowns consisting of the generalized displacements and the stress resultants of the plate. An efficient reduction method is used in conjunction with mixed finite element models for determining the stability boundary and the postbuckling response of the plate. Sensitivity derivatives are evaluated and used to study the sensitivity of the postbuckling response to variations in the different lamination and material parameters of the plate. For quasi-isotropic plates, numerical results are presented showing the effects of variations in the material characteristics and fiber orientation of individual layers on the postbuckling response of the plate.

Noor, Ahmed K.↗

NASA NTRS1993

Composite micromechanical modeling using the boundary element method

The use of the boundary element method for analyzing composite micromechanical behavior is demonstrated. Stress-strain, heat conduction, and thermal expansion analyses are conducted using the boundary element computer code BEST-CMS, and the results obtained are compared to experimental observations, analytical calculations, and finite element analyses. For each of the analysis types, the boundary element results agree reasonably well with the results from the other methodologies, with explainable discrepancies. Overall, the boundary element method shows promise in providing an alternative method to analyze composite micromechanical behavior.

Goldberg, Robert K.↗

NASA NTRS2016

Tangle-Free Finite Element Mesh Motion for Ablation Problems

In numerical simulations involving boundaries that evolve in time, the primary challenge is updating the computational mesh to reflect the physical changes in the domain. In particular, the fundamental objective for any such \mesh motion" scheme is to maintain mesh quality and suppress unphysical geometric anamolies and artifacts. External to a physical process of interest, mesh motion is an added component that determines the specifics of how to move the mesh given certain limited information from the main system. This paper develops a set of boundary conditions designed to eliminate tangling and internal collision within the context of PDE-based mesh motion (linear elasticity). These boundary conditions are developed for two- and three-dimensional meshes. The paper presents detailed algorithms for commonly occuring topological scenarios and explains how to apply them appropriately. Notably, the techniques discussed herein make use of none of the specifics of any particular formulation of mesh motion and thus are more broadly applicable. The two-dimensional algorithms are validated by an extensive verification procedure. Finally, many examples of diverse geometries in both two- and three-dimensions are shown to showcase the capabilities of the tangle-free boundary conditions.

Droba, Justin↗

NASA NTRSDate not supplied

Fluid-Structure Interaction Simulations of the ASPIRE SR01 Supersonic Parachute

High-fidelity computational fluid dynamics (CFD) simulations have so far only played exploratory and supporting roles in the study and qualification of parachutes for planetary entry. The last few years have seen the maturation of coupled computational methods that are capable of modeling the complex fluid-structure interactions between a parachute canopy and the supersonic flow in the wake of an entry vehicle in flight conditions. One of the primary goals of these methods is to predict the peak opening load experienced by the parachute during inflation. The Launch, Ascent, and Vehicle Aerodynamics (LAVA)team is developing efficient, high-fidelity numerical methods to perform such challenging fluid-structure interaction simulations. A loose coupling approach is used to advance the solutions of a Cartesian ghost cell immersed boundary method CFD solver and a finite element computational structural dynamics (CSD) solver in space and time. The coupled solver is employed to simulate the ASPIRE SR01 flight test, where a build-to-print version of the Mars Science Laboratory parachute was inflated in supersonic conditions in the upper terrestrial atmosphere. The simulations conducted in the current paper predict a peak opening load that is within 10% of that from the flight test. Grid convergence with respect to the volume and structural domains is demonstrated, and less than 1% variation in the peak opening load is predicted between all grid resolutions.

ARMD↗