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 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Thermodynamics of hydrogen-helium mixtures at high pressure and finite temperature

A technique is reviewed for calculating thermodynamic quantities for mixtures of light elements at high pressure, in the metallic state. Ensemble averages are calculated with Monte Carlo techniques and periodic boundary conditions. Interparticle potentials are assumed to be coulombic, screened by the electrons in dielectric function theory. This method is quantitatively accurate for alloys at pressures above about 10 Mbar. An alloy of equal parts hydrogen and helium by mass appears to remain liquid and mixed for temperatures above about 3000 K, at pressures of about 15 Mbar. The additive volume law is satisfied to within about 10%, but the Gruneisen equation of state gives poor results. A calculation at 1300 K shows evidence of a hydrogen-helium phase separation.

Hubbard, W. B.↗

NASA NTRS1971

Transpiration and film cooling boundary layer computer program. Volume 2: Computer program and user's manual

A finite difference turbulent boundary layer computer program which allows for mass transfer wall cooling and equilibrium chemistry effects is presented. The program is capable of calculating laminar or turbulent boundary layer solutions for an arbitrary ideal gas or an equilibrium hydrogen oxygen system. Either two dimensional or axisymmetric geometric configurations may be considered. The equations are solved, in nondimension-alized physical coordinates, using the implicit Crank-Nicolson technique. The finite difference forms of the conservation of mass, momentum, total enthalpy and elements equations are linearized and uncoupled, thereby generating easily solvable tridiagonal sets of algebraic equations. A detailed description of the computer program, as well as a program user's manual is provided. Detailed descriptions of all boundary layer subroutines are included, as well as a section defining all program symbols of principal importance. Instructions are then given for preparing card input to the program and for interpreting the printed output. Finally, two sample cases are included to illustrate the use of the program.

Gloss, R. J.↗

NASA NTRS2000

Propulsion Aeroelastic Analysis Developed for Flutter and Forced Response

The NASA Glenn Research Center at Lewis Field develops new technologies to increase the fuel efficiency of aircraft engines, improve the safety of engine operation, reduce emissions, and reduce engine noise. With the development of new designs for fans, compressors, and turbines to achieve these goals, the basic aeroelastic requirements are that there should be no flutter (self-excited vibrations) or high resonant blade stresses (due to forced response) in the operating regime. Therefore, an accurate prediction and analysis capability is required to verify the aeroelastic soundness of the designs. Such a three-dimensional viscous propulsion aeroelastic analysis capability has been developed at Glenn with support from the Advanced Subsonic Technology (AST) program. This newly developed aeroelastic analysis capability is based on TURBO, a threedimensional unsteady aerodynamic Reynolds-averaged Navier-Stokes turbomachinery code developed previously under a grant from Glenn. TURBO can model the viscous flow effects that play an important role in certain aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), flutter in the presence of shock and boundary-layer interaction, and forced response due to wakes and shock impingement. In aeroelastic analysis, the structural dynamics representation of the blades is based on normal modes. A finite-element analysis code is used to calculate these in-vacuum vibration modes and the associated natural frequencies. In an aeroelastic analysis using the TURBO code, flutter and forced response are modeled as being uncoupled. To calculate if a blade row will flutter, one prescribes the motion of the blade to be a harmonic vibration in a specified in-vacuum normal mode. An aeroelastic analysis preprocessor is used to generate the displacement field required for the analysis. The work done by aerodynamic forces on the vibrating blade during a cycle of vibration is calculated. If this work is positive, the blade is dynamically unstable, since it will extract energy from the flow, leading to an increase in the blade s oscillation amplitude. The forced-response excitations on a blade row are calculated by modeling the flow through two adjacent blade rows using the TURBO code. The blades are assumed to be rigid. As an option, a single blade row can be modeled with the upstream blade row influence represented by a time-varying disturbance (gust) at the inlet boundary. The unsteady forces on a blade row from such analyses are used in a structural analysis along with the blade structural dynamics characteristics and aerodynamic damping associated with blade vibration to calculate the resulting dynamic stresses on the blade.

Bakhle, Milind A.↗

NASA NTRS1980

An implicit finite-difference code for inviscid and viscous cascade flow

An implicit finite-difference code is developed to solve either inviscid or viscous flow about two-dimensional cascade blade elements. General coordinate transformations are used so that boundaries can coincide with coordinate lines, and an automatic grid generation routine based on elliptic partial differential equations is employed to mesh arbitrary cascade elements. Characteristic combinations of the differential equations are used at inflow and outflow boundaries. Computed results for both inviscid and viscous flow are compared with other existing cascade solutions and experimental data.

Steger, J. L.↗

NASA NTRS1982

Digital system for structural dynamics simulation

State-of-the-art digital hardware and software for the simulation of complex structural dynamic interactions, such as those which occur in rotating structures (engine systems). System were incorporated in a designed to use an array of processors in which the computation for each physical subelement or functional subsystem would be assigned to a single specific processor in the simulator. These node processors are microprogrammed bit-slice microcomputers which function autonomously and can communicate with each other and a central control minicomputer over parallel digital lines. Inter-processor nearest neighbor communications busses pass the constants which represent physical constraints and boundary conditions. The node processors are connected to the six nearest neighbor node processors to simulate the actual physical interface of real substructures. Computer generated finite element mesh and force models can be developed with the aid of the central control minicomputer. The control computer also oversees the animation of a graphics display system, disk-based mass storage along with the individual processing elements.

Krauter, A. I.↗

NASA NTRS1989

Shear buckling of specially orthotropic plates with centrally located cutouts

The shear loading of a composite rectangular plate with a centrally located circular cutout was analyzed in order to predict the buckling load of the plate. The first step in this analysis is to calculare the in-plane stress distribution of the plate prior to buckling. This problem can be solved using finite element methods, but a quicker and more cost efficient method is desired. The method chosen to determine the pre-buckling stresses is that of boundary collocation using complex variables. The real valued force functions are written in terms of two complex valued functions, each of which is a function of a different complex variable. The solution of the generalized biharmonic function is a superposition of functions of these two complex variables. For this analysis, the force functions are represented with a Laurent series. The constants in these series are found using boundary collocation. Two force equations for complex variables are satisfied over the plate boundaries to obtain these constants. The stresses in the plate are then known for any particular locations on the plate given the applied shear on the edges of the plate.

Klang, Eric C.↗

NASA NTRS1991

Interaction of the sea breeze with a river breeze in an area of complex coastal heating

The interaction of the sea-breeze circulation with a river-breeze circulation in an area of complex coastal heating (east coast of Florida) was studied using a 3D finite-element mesoscale model. The model simulations are compared with temperature and wind fields observed on a typical fall day during the Kennedy Space Center Atmospheric Boundary Layer Experiment. The results from numerical experiments designed to isolate the effect of the river breeze indicate that the convergence in the sea-breeze front is suppressed when it passes over the cooler surface of the rivers.

Zhong, Shiyuan↗

NASA NTRS1994

Out of plane gradients in laminated composites

The research under this grant focused on three different, but related, research projects. All three projects were concerned with predicting the interlaminar stresses that develop in laminated composites in the presence of free edges. The first project was dedicated to the development of a new semi-analytical approximate solution for predicting interlaminar stresses. The second project was devoted to a study of the influence of boundary conditions, geometric parameters, and material modelling on interlaminar stresses in stiffened composite panels. The third project was dedicated to the development of a new finite element for efficient analysis of interlaminar stresses in curved, stiffened panels. A summary of each of the projects is given.

Herakovich, Carl T.↗

NASA NTRS2008

Measurement and Prediction of the Thermomechanical Response of Shape Memory Alloy Hybrid Composite Beams

An experimental and numerical investigation into the static and dynamic responses of shape memory alloy hybrid composite (SMAHC) beams is performed to provide quantitative validation of a recently commercialized numerical analysis/design tool for SMAHC structures. The SMAHC beam specimens consist of a composite matrix with embedded pre-strained SMA actuators, which act against the mechanical boundaries of the structure when thermally activated to adaptively stiffen the structure. Numerical results are produced from the numerical model as implemented into the commercial finite element code ABAQUS. A rigorous experimental investigation is undertaken to acquire high fidelity measurements including infrared thermography and projection moire interferometry for full-field temperature and displacement measurements, respectively. High fidelity numerical results are also obtained from the numerical model and include measured parameters, such as geometric imperfection and thermal load. Excellent agreement is achieved between the predicted and measured results of the static and dynamic thermomechanical response, thereby providing quantitative validation of the numerical tool.

Davis, Brian↗

NASA NTRS2005

High-Aperture-Efficiency Horn Antenna

A horn antenna (see Figure 1) has been developed to satisfy requirements specific to its use as an essential component of a high-efficiency Ka-band amplifier: The combination of the horn antenna and an associated microstrip-patch antenna array is required to function as a spatial power divider that feeds 25 monolithic microwave integrated-circuit (MMIC) power amplifiers. The foregoing requirement translates to, among other things, a further requirement that the horn produce a uniform, vertically polarized electromagnetic field in its patches identically so that the MMICs can operate at maximum efficiency. The horn is fed from a square waveguide of 5.9436-mm-square cross section via a transition piece. The horn features cosine-tapered, dielectric-filled longitudinal corrugations in its vertical walls to create a hard boundary condition: This aspect of the horn design causes the field in the horn aperture to be substantially vertically polarized and to be nearly uniform in amplitude and phase. As used here, cosine-tapered signifies that the depth of the corrugations is a cosine function of distance along the horn. Preliminary results of finite-element simulations of performance have shown that by virtue of the cosine taper the impedance response of this horn can be expected to be better than has been achieved previously in a similar horn having linearly tapered dielectric- filled longitudinal corrugations. It is possible to create a hard boundary condition by use of a single dielectric-filled corrugation in each affected wall, but better results can be obtained with more corrugations. Simulations were performed for a one- and a three-corrugation cosine-taper design. For comparison, a simulation was also performed for a linear- taper design (see Figure 2). The three-corrugation design was chosen to minimize the cost of fabrication while still affording acceptably high performance. Future designs using more corrugations per wavelength are expected to provide better field responses and, hence, greater aperture efficiencies.

Pickens, Wesley↗

NASA NTRS2015

Secondary Instability of Stationary Crossflow Vortices in Mach 6 Boundary Layer Over a Circular Cone

Hypersonic boundary layer flows over a circular cone at moderate incidence can support strong crossflow instability. Due to more efficient excitation of stationary crossflow vortices by surface roughness, such boundary layer flows may transition to turbulence via rapid amplification of the high-frequency secondary instabilities of finite amplitude stationary crossflow vortices. The amplification characteristics of these secondary instabilities are investigated for crossflow vortices generated by an azimuthally periodic array of roughness elements over a 7-degree half-angle circular cone in a Mach 6 free stream. Depending on the local amplitude of the stationary crossflow mode, the most unstable secondary disturbances either originate from the second (i.e., Mack) mode instabilities of the unperturbed boundary layer or correspond to genuine secondary instabilities that reduce to stable disturbances at sufficiently small amplitudes of the stationary crossflow vortex. The predicted frequencies of dominant secondary disturbances are similar to those measured during wind tunnel experiments at Purdue University and the Technical University of Braunschweig, Germany.

Boundary layer flow↗

NASA NTRS1972

Application of finite-element method in the computation of temperature with emphasis on radiative exchanges.

Analyses pertaining to the solution of heat transfer problems in combined modes based on the finite-element method are presented. Two elements - a triangular element employing two spatial variables and a multi-faceted bar element employing one spatial variable - with nonlinear radiation on the boundaries are detailed. The radiative effects considered on the diffuse-gray surface elements include: (1) directional radiant fluxes from distant sources, (2) radiative exchanges including surfaces of prescribed temperatures, and (3) radiative exchanges including elements whose temperatures are not known a priori. The nonlinear part of the boundary condition is treated in two different ways: (1) a consistent linearization method, and (2) a direct energy distribution method. Applications of these elements together with solution algorithms to three sample problems exhibit solution capability and obtainable accuracy.

Lee, H.-P.↗

NASA NTRS1975

Duct acoustics and acoustic finite element method

A comprehensive solution method - Acoustic Finite Element Method - has been proposed. In order to generate confidence in the method proposed, the problem of influence of shear flow in a two-dimensional channel on sound propagation from an arbitrary source distribution has been theoretically formulated. From the preliminary investigation of application of acoustic finite element method, it seems that it can provide solutions to problems concerning the transmission of sound in variable area ducts, or in ducts in which the boundary layer depends on the streamwise variable, or in ducts with discontinuities in lining impedance and with end termination impedances, both the inlet and exhaust jet flows being taken into account.

Kapur, A.↗

NASA NTRS1975

A rigorous solution of the Navier-Stokes equations for unsteady viscous flow at high Reynolds numbers around oscillating airfoils

A method based on the Navier-Stokes equations was developed for analyzing the unsteady incompressible viscous flow around oscillating airfoils at high Reynolds numbers. The Navier-Stokes equations have been integrated in their classical Helmholtz vorticity transport equation form, and the instantaneous velocity field at each time step was determined by the solution of Poisson's equation. A refined finite element was utilized to allow for a conformable solution of the stream function and its first space derivatives at the element interfaces. A corresponding set of accurate boundary conditions was applied; thus obtaining a rigorous solution for the velocity field. The details of the computational procedure and examples of computed results describing the unsteady flow characteristics around the airfoil are presented.

Bratanow, T.↗

NASA NTRS1989

Compact finite volume methods for the diffusion equation

The paper describes an approach to treating initial-boundary-value problems by finite volume methods in which the parallel between differential and difference arguments is closely maintained. By using intrinsic geometrical properties of the volume elements, it is possible to describe discrete versions of the div, curl, and grad operators which lead, using summation-by-parts techniques, to familiar energy equations as well as the div curl = 0 and curl grad = 0 identities. For the diffusion equation, these operators describe compact schemes whose convergence is assured by the energy equations and which yield both the potential and the flux vector with second-order accuracy. A simplified potential form is especially useful for obtaining numerical results by multigrid and ADI methods.

Rose, Milton E.↗

NASA NTRS1992

International Congress on Recent Developments in Air- and Structure-Borne Sound and Vibration, 2nd, Auburn Univ., AL, Mar. 4-6, 1992, Proceedings. Vols. 1-3

The papers presented in this volume cover a variety of topics, including sound intensity, structural intensity, modal analysis and synthesis, statistical energy analysis and energy methods, passive and active damping, and boundary element methods. Attention is also given to diagnostics and condition monitoring, material characterization and nondestructive evaluation, active noise and vibration control, sound radiation and scattering, and finite element analysis.

Crocker, Malcolm J.↗

NASA NTRS1996

The Prediction of Ducted Fan Engine Noise Via a Boundary Integral Equation Method

A computationally efficient Boundary Integral Equation Method (BIEM) for the prediction of ducted fan engine noise is presented. The key features of the BIEM are its versatility and the ability to compute rapidly any portion of the sound field without the need to compute the entire field. Governing equations for the BIEM are based on the assumptions that all acoustic processes are linear, generate spinning modes, and occur in a uniform flow field. An exterior boundary value problem (BVP) is defined that describes the scattering of incident sound by an engine duct with arbitrary profile. Boundary conditions on the duct walls are derived that allow for passive noise control treatment. The BVP is recast as a system of hypersingular boundary integral equations for the unknown duct surface quantities. BIEM solution methodology is demonstrated for the scattering of incident sound by a thin cylindrical duct with hard walls. Numerical studies are conducted for various engine parameters and continuous portions of the total pressure field are computed. Radiation and duct propagation results obtained are in agreement with the classical results of spinning mode theory for infinite ducts.

Boundary integral equation↗

NASA NTRS2012

Vibroacoustic Characterization of a New Hybrid Wing-Body Fuselage Concept

A lighter, more robust airframe design is required to withstand the loading inherent to next generation non cylindrical commercial airliners. The Pultruded Rod Stitched Efficient Unitized Structure concept, a highly integrated composite design involving a stitched and co-cured substructure, has been developed to meet such requirements. While this structure has been shown to meet the demanding out-of-plane loading requirements of the flat-sided pressurized cabin design, there are concerns that the stiff co-cured details will result in relatively high acoustic radiation efficiencies at frequencies well below the thin skin acoustic coincidence frequency. To address this concern and establish a set of baseline vibroacoustic characteristics, a representative test panel was fabricated and a suite of tests were conducted that involved measurements of panel vibration and radiated sound power during point force and diffuse acoustic field excitations. Experimental results are shown and compared with Finite Element and Statistical Energy Analysis model predictions through the use of modal and energy correlation techniques among others. The behavior of the structure subject to turbulent boundary layer excitation is also numerically examined.

Allen, Albert R.↗