Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary element method”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Nonconforming mortar element methods: Application to spectral discretizations

Spectral element methods are p-type weighted residual techniques for partial differential equations that combine the generality of finite element methods with the accuracy of spectral methods. Presented here is a new nonconforming discretization which greatly improves the flexibility of the spectral element approach as regards automatic mesh generation and non-propagating local mesh refinement. The method is based on the introduction of an auxiliary mortar trace space, and constitutes a new approach to discretization-driven domain decomposition characterized by a clean decoupling of the local, structure-preserving residual evaluations and the transmission of boundary and continuity conditions. The flexibility of the mortar method is illustrated by several nonconforming adaptive Navier-Stokes calculations in complex geometry.

Maday, Yvon↗
NASA NTRS1991

Electronic neural network for dynamic resource allocation

A VLSI implementable neural network architecture for dynamic assignment is presented. The resource allocation problems involve assigning members of one set (e.g. resources) to those of another (e.g. consumers) such that the global 'cost' of the associations is minimized. The network consists of a matrix of sigmoidal processing elements (neurons), where the rows of the matrix represent resources and columns represent consumers. Unlike previous neural implementations, however, association costs are applied directly to the neurons, reducing connectivity of the network to VLSI-compatible 0 (number of neurons). Each row (and column) has an additional neuron associated with it to independently oversee activations of all the neurons in each row (and each column), providing a programmable 'k-winner-take-all' function. This function simultaneously enforces blocking (excitatory/inhibitory) constraints during convergence to control the number of active elements in each row and column within desired boundary conditions. Simulations show that the network, when implemented in fully parallel VLSI hardware, offers optimal (or near-optimal) solutions within only a fraction of a millisecond, for problems up to 128 resources and 128 consumers, orders of magnitude faster than conventional computing or heuristic search methods.

Thakoor, A. P.↗
NASA NTRS1992

In-flight surface-flow measurements on a subsonic transport high-lift flap system

As part of a multiphased program for subsonic transport high-lift flight research, flight tests were conducted on the Transport Systems Research Vehicle (B737-100 aircraft) at the NASA Langley Research Center, to obtain detailed flow characteristics of the high-lift flap system for correlation with computational and wind-tunnel investigations. Pressure distributions, skin friction, and flow-visualization measurements were made on a triple-slotted flap system for a range of flap deflections, chord Reynolds numbers (10 to 21 million), and Mach numbers (0.16 to 0.36). Experimental test results are given for representative flap settings indicating flow separation on the fore-flap element for the largest flap deflection. Comparisons of the in-flight flow measurements were made with predictions from available viscous multielement computational methods modified with simple-sweep theory. Computational results overpredicted the experimentally measured pressures, particularly in the case involving separation of the fore lap, indicating the need for better modeling of confluent boundary layers and three-dimensional sweep effects.

Yip, Long P.↗
NASA NTRS1998

The Wavelet Element Method: Realization and Additional Features in 2D and 3D - Part 2

The Wavelet Element Method (WEM) provides a construction of multiresolution systems and biorthogonal wavelets on fairly general domains. These are split into subdomains that are mapped to a single reference hypercube. Tensor products of scaling functions and wavelets defined on the unit interval are used on the reference domain. By introducing appropriate matching conditions across the interelement boundaries, a globally continuous biorthogonal wavelet basis on the general domain is obtained. This construction does not uniquely define the basis functions but rather leaves some freedom for fulfilling additional features. In this paper we detail the general construction principle of the WEM to the 1D, 2D and 3D cases. We address additional features such as symmetry, vanishing moments and minimal support of the wavelet functions in each particular dimension. The construction is illustrated by using biorthogonal spline wavelets on the interval.

Canuto, Claudio↗
NASA NTRS2012

Nonlinear Analysis of Bonded Composite Tubular Lap Joints

The present study describes a semi-analytical solution method for predicting the geometrically nonlinear response of a bonded composite tubular single-lap joint subjected to general loading conditions. The transverse shear and normal stresses in the adhesive as well as membrane stress resultants and bending moments in the adherends are determined using this method. The method utilizes the principle of virtual work in conjunction with nonlinear thin-shell theory to model the adherends and a cylindrical shear lag model to represent the kinematics of the thin adhesive layer between the adherends. The kinematic boundary conditions are imposed by employing the Lagrange multiplier method. In the solution procedure, the displacement components for the tubular joint are approximated in terms of non-periodic and periodic B-Spline functions in the longitudinal and circumferential directions, respectively. The approach presented herein represents a rapid-solution alternative to the finite element method. The solution method was validated by comparison against a previously considered tubular single-lap joint. The steep variation of both peeling and shearing stresses near the adhesive edges was successfully captured. The applicability of the present method was also demonstrated by considering tubular bonded lap-joints subjected to pure bending and torsion.

E. Oterkus↗
NASA NTRS1991

Internally mounted thin-liquid-film skin-friction meter - Comparison with floating element method with and without pressure gradient

A new, robust oil film skin friction meter was designed and constructed. This enables skin friction measurements remotely and from within the model, as well as avoiding the need to know the location of the leading edge of the film. The instrument was tested by comparing measurements with those given by a floating element gage in a zero pressure gradient flat plate turbulent boundary layer. Both instruments agreed satisfactorily with the well-known curve for this case. Significant discrepancies between the two instruments were observed in the case of adverse and favorable pressure gradients. The discrepancies were of opposite sign for opposite-sign pressure gradients as is consistent with the error expected from floating-element gages. Additional confidence in the oil film technique is supplied by the good agreement of the behavior of the film profile with predictions from lubrication theory.

Hornung, Hans↗
NASA NTRS1988

A new boundary element formulation for two- and three-dimensional thermoelasticity using particular integrals

New two- and three-dimensional BEM formulations are developed for steady-state and transient uncoupled thermoelasticity. These new procedures differ from previous work in that additional surface of volume integration is not required to incorporate thermal loads in the analysis. Instead, thermal body forces are introduced in the BEM system via particular integrals. The present formulation is implemented in a general-purpose multiregion system, and examples are presented to demonstrate the accuracy and versatility of the method.

Henry, Donald P., Jr.↗
NASA NTRS2016

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗
NASA NTRS1981

Computer program for aerodynamic and blading design of multistage axial-flow compressors

A code for computing the aerodynamic design of a multistage axial-flow compressor and, if desired, the associated blading geometry input for internal flow analysis codes is presented. Compressible flow, which is assumed to be steady and axisymmetric, is the basis for a two-dimensional solution in the meridional plane with viscous effects modeled by pressure loss coefficients and boundary layer blockage. The radial equation of motion and the continuity equation are solved with the streamline curvature method on calculation stations outside the blade rows. The annulus profile, mass flow, pressure ratio, and rotative speed are input. A number of other input parameters specify and control the blade row aerodynamics and geometry. In particular, blade element centerlines and thicknesses can be specified with fourth degree polynomials for two segments. The output includes a detailed aerodynamic solution and, if desired, blading coordinates that can be used for internal flow analysis codes.

Crouse, J. E.↗
NASA NTRS1985

A first-order time-domain Green's function approach to supersonic unsteady flow

A time-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. 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, in space through standard finite-element technique, and in time through finite-difference, 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 is discretized into hyperboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed to vary linearly within each panel. Consistent with the spatial linear (first-order) finite-element approximations, the potential and co-normalwash are assumed to vary linearly in time. The long range goal of our research is to develop a comprehensive theory for unsteady supersonic potential aerodynamics which is capable of yielding accurate results even in the low supersonic (i.e., high transonic) range.

Freedman, M. I.↗
NASA NTRS1989

Space Shuttle Main Engine structural analysis and data reduction/evaluation. Volume 6: Primary nozzle diffuser analysis

The primary nozzle diffuser routes fuel from the main fuel valve on the Space Shuttle Main Engine (SSME) to the nozzle coolant inlet mainfold, main combustion chamber coolant inlet mainfold, chamber coolant valve, and the augmented spark igniters. The diffuser also includes the fuel system purge check valve connection. A static stress analysis was performed on the diffuser because no detailed analysis was done on this part in the past. Structural concerns were in the area of the welds because approximately 10 percent are in areas inaccessible by X-ray testing devices. Flow dynamics and thermodynamics were not included in the analysis load case. Constant internal pressure at maximum SSME power was used instead. A three-dimensional, finite element method was generated using ANSYS version 4.3A on the Lockheed VAX 11/785 computer to perform the stress computations. IDEAS Supertab on a Sun 3/60 computer was used to create the finite element model. Rocketdyne drawing number RS009156 was used for the model interpretation. The flight diffuser is denoted as -101. A description of the model, boundary conditions/load case, material properties, structural analysis/results, and a summary are included for documentation.

Foley, Michael J.↗
NASA NTRS1995

Numerical study to assess sulfur hexafluoride as a medium for testing multielement airfoils

A methodology is described for computing viscous flows of air and sulfur hexafluoride (SF6). The basis is an existing flow solver that calculates turbulent flows in two dimensions on unstructured triangular meshes. The solver has been modified to incorporate the thermodynamic model for SF6 and used to calculate the viscous flow over two multielement airfoils that have been tested in a wind tunnel with air as the test medium. Flows of both air and SF6 at a free-stream Mach number of 0.2 and a Reynolds number of 9 x 10(exp 6) are computed for a range of angles of attack corresponding to the wind-tunnel test. The computations are used to investigate the suitability of SF6 as a test medium in wind tunnels and are a follow-on to previous computations for single-element airfoils. Surface-pressure, lift, and drag coefficients are compared with experimental data. The effects of heavy gas on the details of the flow are investigated based on computed boundary-layer and skin-friction data. In general, the predictions in SF6 vary little from those in air. Within the limitations of the computational method, the results presented are sufficiently encouraging to warrant further experiments.

Bonhaus, Daryl L.↗
NASA NTRS2016

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation, surface-to-surface radiation exchange, and flowfield coupling. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗
NASA NTRS2016

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation and contact interfaces, and example simulations are included. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗
NASA NTRS2002

Advances in space radiation shielding codes

Early space radiation shield code development relied on Monte Carlo methods and made important contributions to the space program. Monte Carlo methods have resorted to restricted one-dimensional problems leading to imperfect representation of appropriate boundary conditions. Even so, intensive computational requirements resulted and shield evaluation was made near the end of the design process. Resolving shielding issues usually had a negative impact on the design. Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary concept to the final design. For the last few decades, we have pursued deterministic solutions of the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard Finite Element Method (FEM) geometry common to engineering design methods. A single ray trace in such geometry requires 14 milliseconds and limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given.

manned↗
NASA NTRS2001

Computational Aeroacoustics by the Space-time CE/SE Method

In recent years, a new numerical methodology for conservation laws-the Space-Time Conservation Element and Solution Element Method (CE/SE), was developed by Dr. Chang of NASA Glenn Research Center and collaborators. In nature, the new method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its rigorous treatment of the fluxes and geometry, it is different from the existing schemes. The CE/SE scheme features: (1) space and time treated on the same footing, the integral equations of conservation laws are solve( for with second order accuracy, (2) high resolution, low dispersion and low dissipation, (3) novel, truly multi-dimensional, simple but effective non-reflecting boundary condition, (4) effortless implementation of computation, no numerical fix or parameter choice is needed, an( (5) robust enough to cover a wide spectrum of compressible flow: from weak linear acoustic waves to strong, discontinuous waves (shocks) appropriate for linear and nonlinear aeroacoustics. Currently, the CE/SE scheme has been developed to such a stage that a 3-13 unstructured CE/SE Navier-Stokes solver is already available. However, in the present paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen as a prototype and is sketched in Section 2. Then applications of the CE/SE scheme to linear, nonlinear aeroacoustics and airframe noise are depicted in Sections 3, 4, and 5 respectively to demonstrate its robustness and capability.

Loh, Ching Y.↗
NASA NTRS1976

An investigation of flap-lag stability of wind turbine rotors in the presence of velocity gradients and helicopter rotors in forward flight

The flap-lag equations of motion of a torsionally rigid and centrally hinged spring-restrained rigid blade are developed using quasi-steady blade-element aerodynamic theory. These equations have periodic coefficients and are applicable to both wind-turbine rotors with velocity gradients and helicopter rotors in forward flight. By solving these equations both by the Floquet-Liapunov method and by an approximate method, the range of applicability of the latter method is established. Flap-lag and pure flapping stability boundaries illustrating the effects of flow and rotor parameters are presented. Finally, some discussion on the techniques for generating the Floquet transition matrix and on the relative merits of the coordinate axes system is presented.

Kaza, K. R. V.↗
NASA NTRS1986

Evaluation and certification of heater assemblies developed for thermal vacuum acceptance testing

Preparation of Rockwell International's Thermal Vacuum Chamber for acceptance testing of a mass produced satellite required the development of unique quartz lamp and hot wire heater assemblies. Testing performed on the basic elements of these heaters is described, as is the final testing done to certify that the heater assemblies meet the thermal requirements for acceptance testing. The methods and procedures of thermal mapping used during the development and final certification of these heater assemblies are presented. The absence of a definitive standard for determining flux distribution and heating boundaries for heaters of this type required the development of a test plan incorporating several thermal mapping techniques. These techniques include the development of heat flux using a multiple calorimeter array for both vacuum and ambient test conditions, and a photographic method for detecting heating boundaries. The test plan and thermal mapping techniques are discussed.

Allen, J. E.↗