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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

An application of small-gap equations in sealing devices

The study of a thin, incompressible Newtonian fluid layer trapped between two almost parallel, sliding surfaces has been actively pursued in the last decades. This subject includes lubrication applications such as slider bearings or the sealing of non-pressurized fluids with rubber rotary shaft seals. In the present work we analyze numerically the flow of lubricant fluid through a micro-gap of sealing devices. The first stage of this study is carried out assuming that a 'small-gap' parameter delta attains an extreme value in the Navier-Stokes equations. The precise meaning of small-gap is achieved by the particular limit delta = 0 which, within the bounds of the hypotheses, predicts transport of lubricant through the sealed area by centrifugal instabilities. Numerical results obtained with the penalty function approximation in the finite element method are presented. In particular, the influence of inflow and outflow boundary conditions, and their impact in the simulated flow, are discussed.

Vionnet, Carlos A.↗
NASA NTRS1999

Open-Mode Debonding Analysis of Curved Sandwich Panels Subjected to Heating and Cryogenic Cooling on Opposite Faces

Increasing use of curved sandwich panels as aerospace structure components makes it vital to fully understand their thermostructural behavior and identify key factors affecting the open-mode debonding failure. Open-mode debonding analysis is performed on a family of curved honeycomb-core sandwich panels with different radii of curvature. The curved sandwich panels are either simply supported or clamped, and are subjected to uniform heating on the convex side and uniform cryogenic cooling on the concave side. The finite-element method was used to study the effects of panel curvature and boundary condition on the open-mode stress (radial tensile stress) and displacement fields in the curved sandwich panels. The critical stress point, where potential debonding failure could initiate, was found to be at the midspan (or outer span) of the inner bonding interface between the sandwich core and face sheet on the concave side, depending on the boundary condition and panel curvature. Open-mode stress increases with increasing panel curvature, reaching a maximum value at certain high curvature, and then decreases slightly as the panel curvature continues to increase and approach that of quarter circle. Changing the boundary condition from simply supported to clamped reduces the magnitudes of open-mode stresses and the associated sandwich core depth stretching.

Ko, William L.↗
NASA NTRS2003

RACE and Calculations of Three-dimensional Distributed Cavity Phase Shifts

The design for RACE, a Rb-clock flight experiment for the ISS, is described. The cold collision shift and multiple launching (juggling) have important implications for the design and the resulting clock accuracy and stability. We present and discuss the double clock design for RACE. This design reduces the noise contributions of the local oscillator and simplifies and enhances an accuracy evaluation of the clock. As we try to push beyond the current accuracies of clocks, new systematic errors become important. The best fountain clocks are using cylindrical TE(sub 011) microwave cavities. We recently pointed out that many atoms pass through a node of the standing wave microwave field in these cavities. Previous studies have shown potentially large frequency shifts for atoms passing through nodes in a TE(sub 013) cavity. The shift occurs because there is a small traveling wave component due to the absorption of the copper cavity walls. The small traveling wave component leads to position dependent phase shifts. To study these effects, we perform Finite Element calculations. Three-dimensional Finite Element calculations require significant computer resources. Here we show that the cylindrical boundary condition can be Fourier decomposed to a short series of two-dimensional problems. This dramatically reduces the time and memory required and we obtain (3D) phase distributions for a variety of cavities. With these results, we will be able to analyze this frequency shift in fountain and future space clocks.

Li, Ruoxin↗
NASA NTRS2006

Strip Yield Model Numerical Application to Different Geometries and Loading Conditions

A new numerical method based on the strip-yield analysis approach was developed for calculating the Crack Tip Opening Displacement (CTOD). This approach can be applied for different crack configurations having infinite and finite geometries, and arbitrary applied loading conditions. The new technique adapts the boundary element / dislocation density method to obtain crack-face opening displacements at any point on a crack, and succeeds by obtaining requisite values as a series of definite integrals, the functional parts of each being evaluated exactly in a closed form.

Hatamleh, Omar↗
NASA NTRS1991

Shear buckling of specially orthotropic plates with centrally located cutouts

There is significant industry demand for a method of analyzing the shear buckling of composite plates with cutouts. This method should be able to easily accommodate frequent changes in model design; inflexibility being the major drawback to current finite element methods. The approach taken here is broken into two problems, prebuckling and buckling. To solve for the prebuckling stresses, complex variable equations are used in conjunction with boundary collocation. The least squares approach is utilized to improve the accuracy of the results. The buckling problem is solved using the Ritz method. A product of the Ritz method is a complicated integral equation which is solved using numerical integration. To date, the aforementioned method of determining the prebuckling stresses was verified against infinite plate theory and finite elements. Recent efforts have focused on reducing run time, improving accuracy and adding more boundary condition choices.

Klang, Eric C.↗
NASA NTRS1988

A comparison of experimental and calculated thin-shell leading-edge buckling due to thermal stresses

High-temperature thin-shell leading-edge buckling test data are analyzed using NASA structural analysis (NASTRAN) as a finite element tool for predicting thermal buckling characteristics. Buckling points are predicted for several combinations of edge boundary conditions. The problem of relating the appropriate plate area to the edge stress distribution and the stress gradient is addressed in terms of analysis assumptions. Local plasticity was found to occur on the specimen analyzed, and this tended to simplify the basic problem since it effectively equalized the stress gradient from loaded edge to loaded edge. The initial loading was found to be difficult to select for the buckling analysis because of the transient nature of thermal stress. Multiple initial model loadings are likely required for complicated thermal stress time histories before a pertinent finite element buckling analysis can be achieved. The basic mode shapes determined from experimentation were correctly identified from computation.

Jenkins, Jerald M.↗
NASA NTRS1990

Boundary element analysis of sound scattered by a moving surface

A solution for the acoustic field scattered from a uniformly moving rigid body in the presence of a harmonic incident source has been obtained using a boundary integral method. A derivation of the Kirchhoff formula given by Farassat and Myers (1988) for moving surfaces forms the basis for the analysis, and the development of a boundary integral method for the solution of scattering problems from moving rigid bodies is described. Finite elements are used in conjunction with the Galerkin method in order to solve the integral equation that results from the Kirchhoff formula when the observer point is placed on the moving body surface. Once appropriate surface field values are known they are inserted back into the formula in order to predict the field scattered off the body. Tests, including the so called superposition method, are carried out in order to validate the technique and to establish some confidence in its accuracy. Application of the superposition method to moving bodies is presented, and results of the two approaches are discussed. Sample calculations of scattering from a simple body are presented to illustrate the effects of variations in relevant parameters.

Myers, M. K.↗
NASA NTRS1985

Probabilistic Finite Element Development

The probabilistic finite element computer program known as Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) is being developed for the analysis of critical structural components for reusable space propulsion systems. First year efforts involve the formulation of the probabilistic analysis strategy and the development of a probabilistic linear analysis code. The ultimate goal of the 3-year program is the development of a finite element code capable of performing nonlinear dynamic analysis of structures having stochastic material properties, geometry, and boundary conditions and subjected to random loading. Three levels of sophistication are envisioned for the stochastic description of the structural problem, namely: (1) homogeneous random variable for stiffness, mass, damping, and external loading; (2) stochastic characterization of variables at the element level, with specified interelement correlations; and (3) stochastic interpolation of variables within a finite element. Two alternative probabilistic analysis methods will be developed, allowing for all three levels of modeling sophistication.

Nagtegaal, J.↗
NASA NTRS1975

Finite elements and fluid dynamics

Difficulties concerning a use of the finite element method in the solution of the nonlinear equations of fluid dynamics are partly related to various 'hidden' instabilities which often arise in fluid calculations. The instabilities are typically due to boundary effects or nonlinearities. It is shown that in certain cases these instabilities can be avoided if certain conservation laws are satisfied, and that the latter are often intimately related to finite elements.

Fix, G.↗
NASA NTRS1993

Massively parallel computation of RCS with finite elements

One of the promising combinations of finite element approaches for scattering problems uses Whitney edge elements, spherical vector wave-absorbing boundary conditions, and bi-conjugate gradient solution for the frequency-domain near field. Each of these approaches may be criticized. Low-order elements require high mesh density, but also result in fast, reliable iterative convergence. Spherical wave-absorbing boundary conditions require additional space to be meshed beyond the most minimal near-space region, but result in fully sparse, symmetric matrices which keep storage and solution times low. Iterative solution is somewhat unpredictable and unfriendly to multiple right-hand sides, yet we find it to be uniformly fast on large problems to date, given the other two approaches. Implementation of these approaches on a distributed memory, message passing machine yields huge dividends, as full scalability to the largest machines appears assured and iterative solution times are well-behaved for large problems. We present times and solutions for computed RCS for a conducting cube and composite permeability/conducting sphere on the Intel ipsc860 with up to 16 processors solving over 200,000 unknowns. We estimate problems of approximately 10 million unknowns, encompassing 1000 cubic wavelengths, may be attempted on a currently available 512 processor machine, but would be exceedingly tedious to prepare. The most severe bottlenecks are due to the slow rate of mesh generation on non-parallel machines and the large transfer time from such a machine to the parallel processor. One solution, in progress, is to create and then distribute a coarse mesh among the processors, followed by systematic refinement within each processor. Elimination of redundant node definitions at the mesh-partition surfaces, snap-to-surface post processing of the resulting mesh for good modelling of curved surfaces, and load-balancing redistribution of new elements after the refinement are auxiliary steps expected to result in a robust low i/o system for very large finite element problems.

Parker, Jay↗
NASA NTRS1986

A new approach to the solution of boundary value problems involving complex configurations

A new approach for solving certain types of boundary value problems about complex configurations is presented. Numerical algorithms from such diverse fields as finite elements, preconditioned Krylov subspace methods, discrete Fourier analysis, and integral equations are combined to take advantage of the memory, speed and architecture of current and emerging supercomputers. Although the approach has application to many branches of computational physics, the present effort is concentrated in areas of Computational Fluid Dynamics (CFD) such as steady nonlinear aerodynamics, time harmonic unsteady aerodynamics, and aeroacoustics. The most significant attribute of the approach is that it can handle truly arbitrary boundary geometries and eliminates the difficult task of generating surface fitted grids.

Rubbert, P. E.↗
NASA NTRS1985

A first-order Green's function approach to supersonic oscillatory flow: A mixed analytic and numeric treatment

A frequency domain Green's Function Method for unsteady supersonic potential flow around complex aircraft configurations is presented. The focus is on the supersonic range wherein the linear potential flow assumption is valid. In this range the effects of the nonlinear terms in the unsteady supersonic compressible velocity potential equation are negligible and therefore these terms will be omitted. The Green's function method is employed in order to convert the potential flow differential equation into an integral one. This integral equation is then discretized, through standard finite element technique, to yield a linear algebraic system of equations relating the unknown potential to its prescribed co-normalwash (boundary condition) on the surface of the aircraft. The arbitrary complex aircraft configuration (e.g., finite-thickness wing, wing-body-tail) is discretized into hyperboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed to vary linearly within each panel. The long range goal is to develop a comprehensive theory for unsteady supersonic potential aerodynamic which is capable of yielding accurate results even in the low supersonic (i.e., high transonic) range.

Freedman, M. I.↗
NASA NTRS2002

Finite Element A Posteriori Error Estimation for Heat Conduction

This research investigates residual-based a posteriori error estimates for finite element approximations of heat conduction in single-layer and multi-layered materials. The finite element approximation, based upon hierarchical modelling combined with p-version finite elements, is described with specific application to a two-dimensional, steady state, heat-conduction problem. Element error indicators are determined by solving an element equation for the error with the element residual as a source, and a global error estimate in the energy norm is computed by collecting the element contributions. Numerical results of the performance of the error estimate are presented by comparisons to the actual error. Two methods are discussed and compared for approximating the element boundary flux. The equilibrated flux method provides more accurate results for estimating the error than the average flux method. The error estimation is applied to multi-layered materials with a modification to the equilibrated flux method to approximate the discontinuous flux along a boundary at the material interfaces. A directional error indicator is developed which distinguishes between the hierarchical modeling error and the finite element error. Numerical results are presented for single-layered materials which show that the directional indicators accurately determine which contribution to the total error dominates.

Lang, Christapher G.↗
NASA NTRS2004

Deformations of Flat Unsymmetric Laminates Subjected to Inplane Loads

The geometrically nonlinear deformation response of initially flat unsymmetric cross-ply laminates subjected to an inplane compressive load and two sets of boundary conditions is studied. Stability of the deformations is considered. At issue is whether or not the plate remains flat with increased compressive loading, and whether it buckles. A semi-infinite unsymmetric cross-ply laminate is used to show the combined effects of geometric nonlinearities and bending-stretch coupling. Finite element results for finite laminates are then presented, and it is shown that to a large degree the boundary conditions control the character of the deformation response. It appears that clamped boundary conditions support buckling behavior, in the classic sense of bifurcation, whereas simply-supported conditions do not.

Hyer, Michael W.↗
NASA NTRS1989

Approximate minimum-time trajectories for 2-link flexible manipulators

Powell's nonlinear programming code, VF02AD, was used to generate approximate minimum-time tip trajectories for 2-link semi-rigid and flexible manipulator movements in the horizontal plane. The manipulator is modeled with an efficient finite-element scheme for an n-link, m-joint system with horizontal-plane bending only. Constraints on the trajectory include boundary conditions on position and energy for a rest-to-rest maneuver, straight-line tracking between boundary positions, and motor torque limits. Trajectory comparisons utilize a change in the link stiffness, EI, to transition from the semi-rigid to flexible case. Results show the level of compliance necessary to excite significant modal behavior. Quiescence of the final configuration is examined with the finite-element model.

Eisler, G. R.↗
DOE OSTI2021

A data-driven approach to modeling cancer cell mechanics during microcirculatory transport

In order to understand the effect of cellular level features on the transport of circulating cancer cells in the microcirculation, there has been an increasing reliance on high-resolution in silico models. Accurate simulation of cancer cells flowing with blood cells requires resolving cellular-scale interactions in 3D, which is a significant computational undertaking warranting a cancer cell model that is both computationally efficient yet sufficiently complex to capture relevant behavior. Given that the characteristics of metastatic spread are known to depend on cancer type, it is crucial to account for mechanistic behavior representative of a specific cancer’s cells. To address this gap, in the present work we develop and validate a means by which an efficient and popular membrane model-based approach can be used to simulate deformable cancer cells and reproduce experimental data from specific cell lines. Here, cells are modeled using the immersed boundary method (IBM) within a lattice Boltzmann method (LBM) fluid solver, and the finite element method (FEM) is used to model cell membrane resistance to deformation. Through detailed comparisons with experiments, we (i) validate this model to represent cancer cells undergoing large deformation, (ii) outline a systematic approach to parameterize different cell lines to optimally fit experimental data over a range of deformations, and (iii) provide new insight into nucleated vs. non-nucleated cell models and their ability to match experiments. While many works have used the membrane-model based method employed here to model generic cancer cells, no quantitative comparisons with experiments exist in the literature for specific cell lines undergoing large deformation. Here, we describe a phenomenological, data-driven approach that can not only yield good agreement for large deformations, but explicitly detail how it can be used to represent different cancer cell lines. This model is readily incorporated into cell-resolved hemodynamic transport simulations, and thus offers significant potential to complement experiments towards providing new insights into various aspects of cancer progression.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Thermal Management System for an Electric Machine With Additively Manufactured Hollow Conductors With Integrated Heat Pipes

Here, this paper discusses steps taken to develop a novel thermal management system for an aircraft propulsion electric machine containing additively manufactured coils integrated with heat pipes aimed at boosting its specific power. Experimental setups are used to size and characterize heat pipes for the application and 3D thermal finite element analysis is used to determine optimum heat transfer coefficient of the convective boundaries. For some of the convective boundaries, fin-based surface area enhancement is required to reach a target combined overall heat transfer coefficient and surface area performance (UA). This enhancement is worked out using a combination of the Engineering Equation Solver tool and 3D thermal FEA. The thermal management system's UA, and by extension its specific power, sensitivity to coolant temperature is explored. Temperature distribution plots of optimized machine components are also presented and discussed. Lastly, additional heat pipe testing is carried out to study its maximum heat transfer capability's sensitivity to condenser coolant temperature and configuration.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗