Search NASA⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Finite element two-dimensional panel flutter at high supersonic speeds and elevated temperature

A consistent finite element formulation for flutter of two-dimensional panels including effects of temperature variation along panel length is presented. The von Karman nonlinear strain-displacement relation is used to account for the large-deflections, and the quasi-steady first-order piston theory is employed for aerodynamic loading in the finite-element formulation. The panel flutter under combined thermal/aerodynamic loading can be separated into solving two problems in sequence. They are: (1) postbuckling deflection, stresses and buckling stability boundary, and (2) critical dynamic pressure and linear flutter boundary. The solution procedures in solving these problems are presented in detail. Finite element results of two-dimensional panels with uniform temperature change are compared with classic analytical solutions. Influence of nonuniform temperature distributions on critical buckling temperature, flutter deflection and stresses, and flutter boundary is also given.

Xue, David Y.↗
NASA NTRS1983

Finite-element methods for steady solidification problems

Four Galerkin finite-element methods are tested for solving the free-boundary problem that describes steady solidification. The formulations differ in the solution method used to account for the unknown shape of the melt/solid interface, in the interphase condition (either balance of heat flux or equilibrium of temperature) distinguished for locating the interface, and in the technique used for solving the systems of algebraic equations that result from the finite-element approximations. Methods that use the melting point isotherm to locate the melt/solid interface are found more accurate and efficient than formulations based on the interfacial energy balance. Solution by a Galerkin-Newton algorithm of the free-boundary problem transformed to a fixed domain is most efficient when the field problem in each phase is made nonlinear by including radiation from the melt and solid to the surroundings.

Ettouney, H. M.↗
NASA NTRS2012

Multiscale Modeling of Grain Boundaries in ZrB2: Structure, Energetics, and Thermal Resistance

A combination of ab initio, atomistic and finite element methods (FEM) were used to investigate the structures, energetics and lattice thermal conductance of grain boundaries for the ultra high temperature ceramic ZrB2. Atomic models of idealized boundaries were relaxed using density functional theory. Information about bonding across the interfaces was determined from the electron localization function. The Kapitza conductance of larger scale versions of the boundary models were computed using non-equilibrium molecular dynamics. The interfacial thermal parameters together with single crystal thermal conductivities were used as parameters in microstructural computations. FEM meshes were constructed on top of microstructural images. From these computations, the effective thermal conductivity of the polycrystalline structure was determined.

Lawson, John W.↗
DOE OSTI2021

Simulation of the PBX 9501 Self-Ignited Heavily Confined Cookoff Test

The PBX 9501 self-ignited thermally damaged explosive (SITX) test was carried out by Holmes et al. A 6-inch diameter sphere made of the PBX 9501 was heated inside a steel confinement vessel until it underwent a thermal explosion. This test provided measured temperature data that was used to define a boundary condition to carry out post-experiment finite element (FE) thermal simulations and comparisons of the FE thermal response with the experimental time to thermal ignition and internal temperature measurements at the center of the PBX 9501 sphere. Finite element models of this experiment were constructed using the heat transfer code Aria. A cookoff model for PBX 9501 was used to simulate the thermal decomposition of PBX 9501 and to predict the time to thermal ignition. An inverse heat conduction scheme was designed and applied to estimate a 1D temperature boundary condition.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS2023

Uncertainty Analysis of Slug Calorimeters in the HyMETS Arc-Jet Facility

The objective of this work is to perform an uncertainty analysis of the deduced stagnation heat flux environment on a slug calorimeter for conditions that span the performance envelope of the Hypersonic Materials Environmental Test System arc-jet facility located at NASA Langley Research Center. Analytical solutions are developed for boundary-value problems on the slug element accounting for non-ideal effects, including spatial variation in the slug heat flux, multi-dimensional thermal conduction, and back-face losses, which departs from the state-of-the-art method derived from the American Society of Testing and Materials. Boundary-value problem definitions are informed by preliminary finite element thermal analysis of the slug calorimeter assembly (including both slug and housing) and just the slug element. The analytical solutions are presented in a general sense and in a truncated form from error analysis. Results are shown in optimizing and validating the analytical models against available slug back-face thermal data. The optimization results indicate that the appropriate epistemic uncertainty of the deduced stagnation heat flux on the slug calorimeter is at most±2.5% for both a high-and low-enthalpy test condition. In addition, a numerical approach is used to determine the aleatory (probabilistic) uncertainty component in the slug stagnation heat flux by applying a marching least-squares slope routine through the steady-state portion of the slug back-face thermal response. Results indicate a compromise between the number of samples and the filter frequency of slug back-face thermal data points when evaluating the standard deviation of the deduced stagnation heat flux statistics. When combining the mixed uncertainty, both aleatory and epistemic, the interval of uncertainty in the deduced stagnation heat flux is determined to be up to ±4%, which is at least a 60% reduction from the standard uncertainty used in the state-of-the-art method.

uncertainty↗
DOE OSTI2025

Modeling Electrodeposition in 3D Porous Architectures for Solid-State Li-Metal Batteries

Li-metal storage in three-dimensional (3D) electrodes is considered a potential dendrite-mitigation strategy. The large surface area and high porosity of these electrodes result in reduced local Li-plating current densities. The porous topology provides a scaffold for Li-deposition and stripping, maintaining both mechanical integrity and Li accessibility. The goal of this study is to understand how characteristics, such as geometry and material properties, affect the current distribution and deposition pattern. To this end, we developed a computational method to track material growth driven by electrodeposition within a complex geometry. This method ensures that the finite-element discretization remains conforming to the moving boundary while preserving an adequate mesh quality, and thus maintains solution accuracy. Using this new computational tool, we analyze the conditions under which porous anode architectures effectively expand the surface area of the charge-transfer interface, and self-regulate current density and dendrite growth.

3D electrode architectures↗
DOE OSTI2021

Ribbon Wrinkling Origin and Its Effect on Tape-Wound Core Shape

Improved understanding of advanced manufacturing techniques is needed to fabricate next- generation soft magnetic tape-wound cores. These utilize metal amorphous nano-composites (MANCs) in the form of long thin ribbons. MANCs exhibit attractive soft magnetic properties for use in high frequency and high power density magnetic components. Of particular interest here, the permeability and hysteresis of MANC ribbons can be tuned in a strain annealing process. An unexpected consequence is that wider ribbons develop longitudinal wrinkle patterns. This can substantially reduce the tape-wound core stacking factor. However, when winding a wrinkled ribbon into a toroidal core, the wrinkles can completely flatten, or crinkle. By this phenomenon, the stacking factor can be recovered. To understand the crinkling effect we have developed analytical and finite element analysis (FEA) models. Experiments with known boundary conditions demonstrate the transition from wrinkled to crinkled versus winding curvature. The models utilize an elastic modulus measured via ribbon tensile testing and a ribbon thickness value measurement made by profilometry. The analytical model provides insight into the crinkling effect but cannot explain the length of the transition region from wrinkling to crinkling. With FEA, we successfully predict the transition region and depending on wrinkling amplitude can determine when a ribbon will wind into a useable core. Such mechanics-based insight provides critical feedback for allowable processing parameters in the strain annealing process.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS1985

A theoretical study of non-adiabatic surface effects for a model in the NTF cryogenic wind tunnel

A theoretical analysis was made of the severity and effect of nonadiabatic surface conditions for a model in the NTF cryogenic wind tunnel. The nonadiabatic condition arises from heaters that are used to maintain a constant thermal environment for instrumentation internal to the model. The analysis was made for several axi-symmetric representations of a fuselage cavity, using a finite element heat conduction code. Potential flow and boundary layer codes were used to calculate the convection condition for the exterior surface of the model. The results of the steady state analysis show that it is possible to maintain the surface temperature very near the adiabatic value, with the judicious use of insulating material. Even for the most severe nonadiabatic condition studied, the effects on skin friction drag and displacement thickness were only marginally significant. The thermal analysis also provided an estimate of the power required to maintain a specified cavity temperature.

Macha, J. M.↗
NASA NTRS1985

Importance of anisotropic bending stiffnesses on buckling of symmetrically laminated composite plates loaded in compression

The differential equation governing buckling of symmetrically laminated composite plates loaded in compression is presented in nondimensional form. From this equation, nondimensional material coefficients are obtained and a nondimensional parameter is presented that is used to assess when anisotropic bending stiffnesses can be neglected in a buckling analysis. Results obtained using finite element analyses are presented that show how boundary conditions, aspect ratio, fiber orientation, stacking sequence, and thickness affect the importance of the anisotropic bending stiffnesses.

Nemeth, M. P.↗
NASA NTRS1996

Implementation of Hybrid V-Cycle Multilevel Methods for Mixed Finite Element Systems with Penalty

The goal of this paper is the implementation of hybrid V-cycle hierarchical multilevel methods for the indefinite discrete systems which arise when a mixed finite element approximation is used to solve elliptic boundary value problems. By introducing a penalty parameter, the perturbed indefinite system can be reduced to a symmetric positive definite system containing the small penalty parameter for the velocity unknown alone. We stabilize the hierarchical spatial decomposition approach proposed by Cai, Goldstein, and Pasciak for the reduced system. We demonstrate that the relative condition number of the preconditioner is bounded uniformly with respect to the penalty parameter, the number of levels and possible jumps of the coefficients as long as they occur only across the edges of the coarsest elements.

Lai, Chen-Yao G.↗
NASA NTRS2002

Eigenvalues of Rectangular Waveguide Using FEM With Hybrid Elements

A finite element analysis using hybrid triangular-rectangular elements is developed to estimate eigenvalues of a rectangular waveguide. Use of rectangular vector-edge finite elements in the vicinity of the PEC boundary and triangular elements in the interior region more accurately models the physical nature of the electromagnetic field, and consequently quicken the convergence.

Deshpande, Manohar D.↗
NASA NTRS1992

Multi-Mode Excitation and Data Reduction for Fatigue Crack Characterization in Conducting Plates

Advances in the technique of fatigue crack characterization by resonant modal analysis have been achieved through a new excitation mechanism and data reduction of multiple resonance modes. A non-contacting electromagnetic device is used to apply a time varying Lorentz force to thin conducting sheets. The frequency and direction of the Lorentz force are such that resonance modes are generated in the test sample. By comparing the change in frequency between distinct resonant modes of a sample, detecting and sizing of fatigue cracks are achieved and frequency shifts caused by boundary condition changes can be discriminated against. Finite element modeling has been performed to verify experimental results.

Wincheski, B.↗
NASA NTRS1992

Mathematical modelling of flow distribution in the human cardiovascular system

The paper presents a detailed model of the entire human cardiovascular system which aims to study the changes in flow distribution caused by external stimuli, changes in internal parameters, or other factors. The arterial-venous network is represented by 325 interconnected elastic segments. The mathematical description of each segment is based on equations of hydrodynamics and those of stress/strain relationships in elastic materials. Appropriate input functions provide for the pumping of blood by the heart through the system. The analysis employs the finite-element technique which can accommodate any prescribed boundary conditions. Values of model parameters are from available data on physical and rheological properties of blood and blood vessels. As a representative example, simulation results on changes in flow distribution with changes in the elastic properties of blood vessels are discussed. They indicate that the errors in the calculated overall flow rates are not significant even in the extreme case of arteries and veins behaving as rigid tubes.

NASA Center JSC↗
NASA NTRS1978

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.↗
NASA NTRSDate not supplied

Effect of Boundary Conditions on Process-Induced Stresses in a Plain Weave Unit Cell

Woven polymer matrix composites (PMCs) are leveraged in aerospace applications for their desirable specific properties, yet they are vulnerable to high residual stresses during manufacturing and their complex geometry makes experimental results difficult to observe. Process modeling is needed to characterize the effects of the curing and predict end stress states. Finite element software can be used to model woven architectures, however accurate representation of processing conditions remains a challenge when it comes to selecting boundary conditions. The effect of BCs on process-induced stress within woven PMCs is studied. The commercial Finite Element Analysis (FEA) software Abaqus is coupled with user-written subroutines in a process modeling framework. A two-dimensionally (2D) woven PMC repeating unit cell (RUC) is modeled with TexGen and Abaqus. Virtual curing is imposed on the bulk matrix. The BC study is conducted with Free, Periodic, Flat, and Flat-Free configurations. Results show that the end stress state is sensitive to the boundary condition assumptions. Flat BC results show great agreement with Periodic BCs. Residual stress results from process modeling are then compared with a linear-elastic thermal cooldown analysis in Abaqus. Cooldown results indicate an overestimation in matrix stresses compared with process modeling.

micromechanics↗
NASA NTRS1994

A mathematical model for simulating noise suppression of lined ejectors

A mathematical model containing the essential features embodied in the noise suppression of lined ejectors is presented. Although some simplification of the physics is necessary to render the model mathematically tractable, the current model is the most versatile and technologically advanced at the current time. A system of linearized equations and the boundary conditions governing the sound field are derived starting from the equations of fluid dynamics. A nonreflecting boundary condition is developed. In view of the complex nature of the equations, a parametric study requires the use of numerical techniques and modern computers. A finite element algorithm that solves the differential equations coupled with the boundary condition is then introduced. The numerical method results in a matrix equation with several hundred thousand degrees of freedom that is solved efficiently on a supercomputer. The model is validated by comparing results either with exact solutions or with approximate solutions from other works. In each case, excellent correlations are obtained. The usefulness of the model as an optimization tool and the importance of variable impedance liners as a mechanism for achieving broadband suppression within a lined ejector are demonstrated.

Watson, Willie R.↗