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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Finite element analysis of low speed viscous and inviscid aerodynamic flows

A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.

Baker, A. J.↗

NASA NTRS1984

Hybrid stress finite elements for large deformations of inelastic solids

A new hybrid stress finite element algorithm, based on a generalization of Fraeijs de Veubeke's complementary energy principle is presented. Analyses of large quasistatic deformation of inelastic solids (hypoelastic, plastic, viscoplastic) are within its capability. Principle variables in the formulation are the nominal stress rate and spin. A brief account is given of the boundary value problem in these variables, and the 'equivalent' variational principle. The finite element equation, along with initial positions and stresses, comprise an initial value problem. Factors affecting the choice of time integration schemes are discussed. Results found by application of the new algorithm are compared to those obtained by a velocity based finite element algorithm.

Reed, K. W.↗

NASA NTRS2008

Multiscale Modeling of Intergranular Fracture in Aluminum: Constitutive Relation For Interface Debonding

Intergranular fracture is a dominant mode of failure in ultrafine grained materials. In the present study, the atomistic mechanisms of grain-boundary debonding during intergranular fracture in aluminum are modeled using a coupled molecular dynamics finite element simulation. Using a statistical mechanics approach, a cohesive-zone law in the form of a traction-displacement constitutive relationship, characterizing the load transfer across the plane of a growing edge crack, is extracted from atomistic simulations and then recast in a form suitable for inclusion within a continuum finite element model. The cohesive-zone law derived by the presented technique is free of finite size effects and is statistically representative for describing the interfacial debonding of a grain boundary (GB) interface examined at atomic length scales. By incorporating the cohesive-zone law in cohesive-zone finite elements, the debonding of a GB interface can be simulated in a coupled continuum-atomistic model, in which a crack starts in the continuum environment, smoothly penetrates the continuum-atomistic interface, and continues its propagation in the atomistic environment. This study is a step towards relating atomistically derived decohesion laws to macroscopic predictions of fracture and constructing multiscale models for nanocrystalline and ultrafine grained materials.

Yamakov, V.↗

NASA NTRS2006

Shielding Effectiveness in a Two-Dimensional Reverberation Chamber Using Finite-Element Techniques

Reverberation chambers are attaining an increased importance in determination of electromagnetic susceptibility of avionics equipment. Given the nature of the variable boundary condition, the ability of a given source to couple energy into certain modes and the passband characteristic due the chamber Q, the fields are typically characterized by statistical means. The emphasis of this work is to apply finite-element techniques at cutoff to the analysis of a two-dimensional structure to examine the notion of shielding-effectiveness issues in a reverberating environment. Simulated mechanical stirring will be used to obtain the appropriate statistical field distribution. The shielding effectiveness (SE) in a simulated reverberating environment is compared to measurements in a reverberation chamber. A log-normal distribution for the SE is observed with implications for system designers. The work is intended to provide further refinement in the consideration of SE in a complex electromagnetic environment.

Bunting, Charles F.↗

NASA NTRS2013

Interior Noise Predictions in the Preliminary Design of the Large Civil Tiltrotor (LCTR2)

A prediction scheme was established to compute sound pressure levels in the interior of a simplified cabin model of the second generation Large Civil Tiltrotor (LCTR2) during cruise conditions, while being excited by turbulent boundary layer flow over the fuselage, or by tiltrotor blade loading and thickness noise. Finite element models of the cabin structure, interior acoustic space, and acoustically absorbent (poro-elastic) materials in the fuselage were generated and combined into a coupled structural-acoustic model. Fluctuating power spectral densities were computed according to the Efimtsov turbulent boundary layer excitation model. Noise associated with the tiltrotor blades was predicted in the time domain as fluctuating surface pressures and converted to power spectral densities at the fuselage skin finite element nodes. A hybrid finite element (FE) approach was used to compute the low frequency acoustic cabin response over the frequency range 6-141 Hz with a 1 Hz bandwidth, and the Statistical Energy Analysis (SEA) approach was used to predict the interior noise for the 125-8000 Hz one-third octave bands.

Grosveld, Ferdinand W.↗

NASA NTRS1993

Numerical and experimental investigation of the bending response of thin-walled composite cylinders

A numerical and experimental investigation of the bending behavior of six eight-ply graphite-epoxy circular cylinders is presented. Bending is induced by applying a known end-rotation to each end of the cylinders, analogous to a beam in bending. The cylinders have a nominal radius of 6 inches, a length-to-radius ratio of 2 and 5, and a radius-to-thickness ratio of approximately 160. A (+/- 45/0/90)S quasi-isotropic layup and two orthotropic layups, (+/- 45/0 sub 2)S and (+/- 45/90 sub 2)S, are studied. A geometrically nonlinear special-purpose analysis, based on Donnell's nonlinear shell equations, is developed to study the prebuckling responses and gain insight into the effects of non-ideal boundary conditions and initial geometric imperfections. A geometrically nonlinear finite element analysis is utilized to compare with the prebuckling solutions of the special-purpose analysis and to study the buckling and post buckling responses of both geometrically perfect and imperfect cylinders. The imperfect cylinder geometries are represented by an analytical approximation of the measured shape imperfections. Extensive experimental data are obtained from quasi-static tests of the cylinders using a test fixture specifically designed for the present investigation. A description of the test fixture is included. The experimental data are compared to predictions for both perfect and imperfect cylinder geometries. Prebuckling results are presented in the form of displacement and strain profiles. Buckling end-rotations, moments, and strains are reported, and predicted mode shapes are presented. Observed and predicted moment vs. end-rotation relations, deflection patterns, and strain profiles are illustrated for the post buckling responses. It is found that a geometrically nonlinear boundary layer behavior characterizes the prebuckling responses. The boundary layer behavior is sensitive to laminate orthotropy, cylinder geometry, initial geometric imperfections, applied end-rotation, and non-ideal boundary conditions. Buckling end-rotations, strains, and moments are influenced by laminate orthotropy and initial geometric imperfections. Measured buckling results correlate well with predictions for the geometrically imperfect specimens. The postbuckling analyses predict equilibrium paths with a number of scallop-shaped branches that correspond to unique deflection patterns. The observed postbuckling deflection patterns and measured strain profiles show striking similarities to the predictions in some cases. Ultimate failure of the cylinders is attributed to an interlaminar shear failure mode along the nodal lines of the postbuckling deflection patterns.

Fuchs, J. P.↗

NASA NTRS2011

Ascent Heating Thermal Analysis on Spacecraft Adaptor Fairings

When the Crew Exploration Vehicle (CEV) is launched, the spacecraft adaptor (SA) fairings that cover the CEV service module (SM) are exposed to aero heating. Thermal analysis is performed to compute the fairing temperatures and to investigate whether the temperatures are within the material limits for nominal ascent aeroheating case. The ascent heating is analyzed by using computational fluid dynamics (CFD) and engineering codes at Marshall Space Flight Center. The aeroheating environment data used for this work is known as Thermal Environment 3 (TE3) heating data. One of the major concerns is with the SA fairings covering the CEV SM and the SM/crew launch vehicle (CLV) flange interface. The TE3 heating rate is a function of time, wall temperature, and the spatial locations. The implementation of the TE3 heating rate as boundary conditions in the thermal analysis becomes challenging. The ascent heating thermal analysis on SA fairings and SM/CLV flange interface are performed using two commercial software packages: Cullimore & Ring (C&R) Thermal Desktop (TD) 5.1 and MSC Patran 2007r1 b. TD is the pre-and post-processor for SINDA, which is a finite-difference-based solver. In TD, the geometry is built and meshed, the boundary conditions are defined, and then SINDA is used to compute temperatures. MSC Pthermal is a finite-element- based thermal solver. MSC Patran is the pre- and post-processor for Pthermal. Regarding the boundary conditions, the convection, contact resistance, and heat load can be imposed in different ways in both programs. These two software packages are used to build the thermal model for the same analysis to validate each other and show the differences in the modeling details.

Wang, Xiao Yen↗

NASA NTRS2000

Higher Order Modeling In the BEM/FEM Hybrid Formulation

Hybrid formulations using low order curl-conforming bases to represent the total electric field within a finite element region and low order divergence-conforming bases to represent equivalent electric and magnetic currents on the boundary are well known. However, higher-order divergence and curl-conforming bases have been shown to provide significant benefits in convergence rates and accuracy when employed in strictly integral equation and strictly finite element formulations. In this paper, a hybrid electric field formulation employing higher order bases is presented. The paper addresses benefits and issues associated with using higher order divergence-and curl-conforming bases in the hybrid finite element/boundary element electric field formulation. The method of singularity subtraction may be used to compute the self terms of the boundary integral when the bases are of low order. But this method becomes laborious and requires great care when the divergence conforming bases are of higher order. In order to handle these singularities simply and accurately, a generalized Gaussian quadrature method is employed in which the expansion functions account for the singularity. In preliminary tests of the higher order hybrid formulation, the equivalent electric current induced by scattering of a plane wave from a square dielectric cylinder is examined. Accurate results are obtained using only a two-triangle mesh when the current basis is of order 4 or 5. Additional results are presented comparing the error obtained using higher order bases to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained with higher order bases are compared to those obtained with lower order bases for selected sample problems.

Fink, Patrick W.↗

NASA NTRSDate not supplied

Finite Element Modeling for Compression Strength Prediction of Honeycomb Cores with Geometric Imperfections Measured using X-ray CT Imaging

Aluminum honeycomb cores have been used extensively in composite sandwich panels due to having high bending rigidity while maintaining low density. During manufacturing of the aluminum honeycomb, the thin metal walls are susceptible to imperfections that deviate from an ideal honeycomb shape. Computational tools to quantify imperfections and investigate their effects on quasi-static compression and impact response can inform design criteria to construct safer and lighter launch vehicle structures. This paper presents the results of a study on finite element modeling of metallic honeycomb cores (HCC) with geometric imperfections measured using X-ray Computer Tomography (CT), to predict the compression response. Finite element (FE) models are constructed with measured imperfections with appropriate boundary conditions. It is shown that the average behavior of larger 10x11 cell models can be predicted by sampling single cell models from the larger domain. Amplification of imperfections for cells at the center of sampled specimens are mode switched from the larger out-of-plane imperfection to the mode shape of the adjacent cells resulting in strengthening of the center specimens.

aluminum honeycomb↗

NASA NTRS1980

Stress-intensity factors for internal surface cracks in cylindrical pressure vessels

The paper presents stress-intensity factors for a wide range of semi-elliptical surface cracks on the inside of pressurized cylinders. The ratio of crack depth to crack length ranged from 0.2 to 1; the ratio of crack depth to wall thickness ranged from 0.2 to 0.8; and the ratio of wall thickness to vessel radius was 0.1 to to 0.25. The stress-intensity factors were calculated by a three-dimensional finite-element method using singularity elements along the crack front and linear-strain elements elsewhere. An equation for the stress-intensity factors was obtained which applies over a wide range of configuration parameters and was within about 5 percent of the present results. A comparison was also made between the results and other analyses of internal surface cracks in cylinders. The results from a boundary-integral equation method were in agreement and those from another finite-element method were in fair agreement (+ or - 8 percent) with the results.

Newman, J. C., Jr.↗

DOE OSTI2020

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2010

Probabilistic Structural Analysis Program

NASA/NESSUS 6.2c is a general-purpose, probabilistic analysis program that computes probability of failure and probabilistic sensitivity measures of engineered systems. Because NASA/NESSUS uses highly computationally efficient and accurate analysis techniques, probabilistic solutions can be obtained even for extremely large and complex models. Once the probabilistic response is quantified, the results can be used to support risk-informed decisions regarding reliability for safety-critical and one-of-a-kind systems, as well as for maintaining a level of quality while reducing manufacturing costs for larger-quantity products. NASA/NESSUS has been successfully applied to a diverse range of problems in aerospace, gas turbine engines, biomechanics, pipelines, defense, weaponry, and infrastructure. This program combines state-of-the-art probabilistic algorithms with general-purpose structural analysis and lifting methods to compute the probabilistic response and reliability of engineered structures. Uncertainties in load, material properties, geometry, boundary conditions, and initial conditions can be simulated. The structural analysis methods include non-linear finite-element methods, heat-transfer analysis, polymer/ceramic matrix composite analysis, monolithic (conventional metallic) materials life-prediction methodologies, boundary element methods, and user-written subroutines. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. NASA/NESSUS 6.2c is structured in a modular format with 15 elements.

Pai, Shantaram S.↗

NASA NTRS1997

Analysis of Smart Composite Structures Including Debonding

Smart composite structures with distributed sensors and actuators have the capability to actively respond to a changing environment while offering significant weight savings and additional passive controllability through ply tailoring. Piezoelectric sensing and actuation of composite laminates is the most promising concept due to the static and dynamic control capabilities. Essential to the implementation of these smart composites are the development of accurate and efficient modeling techniques and experimental validation. This research addresses each of these important topics. A refined higher order theory is developed to model composite structures with surface bonded or embedded piezoelectric transducers. These transducers are used as both sensors and actuators for closed loop control. The theory accurately captures the transverse shear deformation through the thickness of the smart composite laminate while satisfying stress free boundary conditions on the free surfaces. The theory is extended to include the effect of debonding at the actuator-laminate interface. The developed analytical model is implemented using the finite element method utilizing an induced strain approach for computational efficiency. This allows general laminate geometries and boundary conditions to be analyzed. The state space control equations are developed to allow flexibility in the design of the control system. Circuit concepts are also discussed. Static and dynamic results of smart composite structures, obtained using the higher order theory, are correlated with available analytical data. Comparisons, including debonded laminates, are also made with a general purpose finite element code and available experimental data. Overall, very good agreement is observed. Convergence of the finite element implementation of the higher order theory is shown with exact solutions. Additional results demonstrate the utility of the developed theory to study piezoelectric actuation of composite laminates with pre-existing debonding. Significant changes in the modes shapes and reductions in the control authority result due to partially debonded actuators. An experimental investigation addresses practical issues, such as circuit design and implementation, associated with piezoelectric sensing and actuation of composite laminates. Composite specimens with piezoelectric transducers were designed, constructed and tested to validate the higher order theory. These specimens were tested with various stacking sequences, debonding lengths and gains for both open and closed loop cases. Frequency changes of 15% and damping on the order of more than 20% of critical damping, via closed loop control, was achieved. Correlation with the higher order theory is very good. Debonding is shown to adversely affect the open and closed loop frequencies, damping ratios, settling time and control authority.

Chattopadhyay, Aditi↗

DOE OSTI2020

Construction of a Code Verification Matrix for Heat Conduction With Finite Element Code Applications

When establishing the pedigree of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model. During this process, numerical results on various meshes are systematically compared to a reference analytic solution. The selection of analytic solutions can be a laborious process, as it is difficult to establish adequate code confidence without performing redundant work. In this report we address this issue by applying a physics-based process that establishes a set of reference problems. In this process, code simulation options are categorized and systematically tested, which ensures that gaps in testing are easily identified and addressed. The resulting problems are primarily intended for code verification analysis but may also be useful for comparison to other simulation codes, troubleshooting activities, or training exercises. The process is used to select fifteen code verification problems relevant for the one-dimensional steady-state heat conduction equation. These problems are applicable to a wide variety of simulation tools, but, in this work, a demonstration is performed using the finite element-based nuclear fuel performance code BISON. Convergence to the analytic solution at the theoretical rate is quantified for a selection of the problems, which establishes a baseline pedigree for the code. Not only can this standard set of conduction solutions be used for verification of other codes, but also the physics-based process for selecting problems can be utilized to quantify and expand testing for any simulation tool.

42 ENGINEERING↗

DOE OSTI2020

Ultrasonic waves from radial mode excitation of a piezoelectric disc on the surface of an elastic solid

In this work, we discuss the ultrasonic wave propagation characteristics in an isotropic elastic solid due to the radial mode excitation of a piezoelectric disc actuator bonded to its surface. Finite element simulations using coupled electromechanical modeling is employed to investigate the wave propagation behavior. We find that the radial mode vibrations on the surface of the elastic solid generates all the three types of ultrasonic waves: longitudinal, shear, and surface waves. The waves in the solid are comprised of a central lobe and multiple side lobes based on the frequency of the radial mode excitation. The central lobe is predominantly composed of longitudinal waves and the side lobes are composed of shear waves. While the longitudinal waves have a strong central lobe and weak side lobes that are not fully developed, shear waves consists of fully developed side lobes. In addition, we observe that the longitudinal waves have fewer side lobes inside the solid compared to the shear waves. A semi-analytical approach is presented to explain the above observations. The interface traction boundary condition between the piezoelectric disc and the elastic solid is approximated as a truncated 'Bessel' excitation over the area of the piezoelectric disc in contact with the solid. A Fourier–Bessel series expansion technique is used to investigate the wave propagation behavior from the above traction boundary condition. The results obtained explain the observations from the finite element simulations with regard to the number and distribution of side lobes pertaining to the longitudinal and shear waves generated from the radial modes of the piezoelectric disc. The proposed semi-analytical approach is general and can be applied to any arbitrary axisymmetric excitations on the surface of an elastic solid.

36 MATERIALS SCIENCE↗

NASA NTRS1972

Structural evaluation of Marman V-band coupling and flange with conoseal gasket

Results are described of a development test program directed at evaluating the structural capabilities of the Marman V-band coupling and flange with conoseal gasket. The intended end use was for the 75K NERVA flight engine propellant lines. Of major importance in the structural evaluation was the ability to predict stresses throughout the assembly for a variety of loading conditions. Computer finite element analysis was used to predict these stresses but, for the subject configuration, large uncertainties were introduced in modeling the complex geometry and boundary conditions. The purpose of the structural tests was to obtain actual stresses and deflections for correlation with, and updating of the finite element model. Results of the incomplete test program are inconclusive with respect to determining suitability for use on the NERVA engine.

Oates, J. H.↗

DOE OSTI2024

Computational modeling of phononic pseudocrystal isolators

Methods for the efficient computational prediction of the performance of phononic pseudocrystals (structured materials capable of blocking extraordinary ranges of frequency) in COMSOL and other comparable finite element method codes are set forth. These methods include boundary conditions that make possible halving the size of the computational domain. Also included is an introduction of elastic energy density methods for assessing the extinction of elastic waves within the patterned region.

Swift, Stephen Hales↗

NASA NTRS1979

Stress-concentration factors for finite orthotropic laminates with a circular hole and uniaxial loading

Stresses were calculated for finite-width orthotropic laminates with a circular hole and remote uniaxial loading using a two-dimensional finite element analysis with both uniform stress and uniform displacement boundary conditions. Five different laminates were analyzed. Computed results are presented for selected combinations of hole diameter/sheet-width ratio d/w and length-to-width ratio L/w. For small L/w values, the stress-concentration factors K sub tn were significantly different for the uniform stress and uniform displacement boundary conditions. Typically, for the uniform stress condition, the K sub tn values were much larger than for the infinite strip reference condition; however, for the uniform displacement condition, they were only slightly smaller than for this reference. The results for long strips are also presented as width correction factors. For d/w less or = 0.33, these width correction factors are nearly equal for all five laminates.

Hong, C. S.↗