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.

SourceResearch connectionFrom yearThrough year

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

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

NASA NTRS1978

Two-dimensional computer code for design assessment of SRM aft closure environment

A general interpolants method for constructing numerical analogs of the partial differential equations of continuum mechanics and which combines the best features of both finite element and finite difference methods was used in a two dimensional axisymmetric analysis of the space shuttle solid motor aft closure region flowfield. The numerical technique used is discussed as well as the results of the analysis. A steady state solution for the submerged nozzle nose region indicates the development of an area gas flow recirculation near the propellant burning surface boundary and in the region of the underside of the nozzle nose. This recirculation region remains spatially fixed for successive iterations (greater than 800) of the solution.

Freeman, J. A.↗

NASA NTRS1985

On the Interconnection of Incompatible Solid Finite Element Meshes Using Multipoint Constraints

Incompatible meshes, i.e., meshes that physically must have a common boundary, but do not necessarily have coincident grid points, can arise in the course of a finite element analysis. For example, two substructures may have been developed at different times for different purposes and it becomes necessary to interconnect the two models. A technique that uses only multipoint constraints, i.e., MPC cards (or MPCS cards in substructuring), is presented. Since the method uses only MPC's, the procedure may apply at any stage in an analysis; no prior planning or special data is necessary.

Fox, G. L.↗

NASA NTRS1987

Probabilistic structural analysis to evaluate the structural durability of SSME critical components

NASA Lewis Research Center is currently developing probabilistic structural analysis methods for select Space Shuttle Main Engine (SSME) structural components. Briefly, the deterministic, three-dimensional, inelastic analysis methodology developed under the Hot Section Technology (HOST) and R and T Base Programs is being augmented to accommodate the complex probabilistic loading spectra, the thermoviscoplastic material behavior, and the material degradation associated with the environment of space propulsion system structural components representative of the SSME such as turbine blades, transfer ducts, and liquid-oxygen posts. The development of probabilistic structural analysis methodology consists of the following program elements: (1) composite load spectra; (2) probabilistic structural analysis methods; (3) probabilistic finite element theory - new variational principles; and (4) probabilistic structural analysis application. In addition, the program includes deterministic analysis elements: (1) development of structural tailoring computer codes (SSME/STAEBL); (2) development of dynamic creep buckling/ratcheting theory; (3) evaluation of the dynamic characteristics of single-crystal SSME blades; (4) development of SSME blade damper technology; and (5) development of integrated boundary elements for hotfluid structure interaction.

Chamis, Christos C.↗

NASA NTRS1986

A Petrov-Galerkin finite element method for the compressible Euler and Navier-Stokes equations

An overview of a new finite element method for the compressible Euler and Navier-Stokes equations is presented. The discretization is based on entropy variables. The method is developed within the framework of a Petrov-Galerkin formulation. Two perturbations are added to the weighting function; one is a generalization of the SUPG operator and the other is designed to enhance shock capturing capability. The treatment of boundary conditions and the consistent calculation of boundary fluxes are addressed. Results of numerical tests are presented which confirm the robustness and wide applicability of the method.

Franca, L. P.↗

NASA NTRS1987

Model-size reduction for the analysis of symmetric structures with asymmetric boundary conditions

A simple computational procedure is presented for reducing the size of the analysis model for a symmetric structure with asymmetric boundary conditions to that of the corresponding structure with symmetric boundary conditions. The procedure is based on approximating the asymmetric response of the structure by a linear combination of symmetric and antisymmetric global approximation vectors (or modes). The key elements of the procedure are (1) restructuring the governing finite-element equations to delineate the contributions to the symmetric and antisymmetric components of the asymmetric response, (2) successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite-element method is first used to generate a few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The effectiveness of the computational procedure is demonstrated by means of numerical examples of linear static problems of shells, and its potential for solving nonlinear problems is discussed.

Noor, Ahmed K.↗

NASA NTRS1984

Large-Scale Advanced Prop-Fan (LAP) blade design

This report covers the design analysis of a very thin, highly swept, propeller blade to be used in the Large-Scale Advanced Prop-Fan (LAP) test program. The report includes: design requirements and goals, a description of the blade configuration which meets requirements, a description of the analytical methods utilized/developed to demonstrate compliance with the requirements, and the results of these analyses. The methods described include: finite element modeling, predicted aerodynamic loads and their application to the blade, steady state and vibratory response analyses, blade resonant frequencies and mode shapes, bird impact analysis, and predictions of stalled and unstalled flutter phenomena. Summarized results include deflections, retention loads, stress/strength comparisons, foreign object damage resistance, resonant frequencies and critical speed margins, resonant vibratory mode shapes, calculated boundaries of stalled and unstalled flutter, and aerodynamic and acoustic performance calculations.

Violette, John A.↗

NASA NTRS1988

Efficient vibration mode analysis of aircraft with multiple external store configurations

A coupling method for efficient vibration mode analysis of aircraft with multiple external store configurations is presented. A set of low-frequency vibration modes, including rigid-body modes, represent the aircraft. Each external store is represented by its vibration modes with clamped boundary conditions, and by its rigid-body inertial properties. The aircraft modes are obtained from a finite-element model loaded by dummy rigid external stores with fictitious masses. The coupling procedure unloads the dummy stores and loads the actual stores instead. The analytical development is presented, the effects of the fictitious mass magnitudes are discussed, and a numerical example is given for a combat aircraft with external wing stores. Comparison with vibration modes obtained by a direct (full-size) eigensolution shows very accurate coupling results. Once the aircraft and stores data bases are constructed, the computer time for analyzing any external store configuration is two to three orders of magnitude less than that of a direct solution.

Karpel, M.↗

NASA NTRS1997

The Reduction of Ducted Fan Engine Noise Via A Boundary Integral Equation Method

The development of a Boundary Integral Equation Method (BIEM) for the prediction of ducted fan engine noise is discussed. The method is motivated by the need for an efficient and versatile computational tool to assist in parametric noise reduction studies. In this research, the work in reference 1 was extended to include passive noise control treatment on the duct interior. The BEM considers the scattering of incident sound generated by spinning point thrust dipoles in a uniform flow field by a thin cylindrical duct. The acoustic field is written as a superposition of spinning modes. Modal coefficients of acoustic pressure are calculated term by term. The BEM theoretical framework is based on Helmholtz potential theory. A boundary value problem is converted to a boundary integral equation formulation with unknown single and double layer densities on the duct wall. After solving for the unknown densities, the acoustic field is easily calculated. The main feature of the BIEM is the ability to compute any portion of the sound field without the need to compute the entire field. Other noise prediction methods such as CFD and Finite Element methods lack this property. Additional BIEM attributes include versatility, ease of use, rapid noise predictions, coupling of propagation and radiation both forward and aft, implementable on midrange personal computers, and valid over a wide range of frequencies.

Tweed, J.↗

NASA NTRSDate not supplied

Transonic Cruise Slotted Wing Design for Commercial Transport Aircraft using CDISC

A knowledge-based aerodynamic design method, CDISC, has been extended to enable the computational design of transonic cruise slotted wings for commercial transport aircraft. The cruise slotted wing is a multielement wing concept with a forward main element and an aft flap element, separated to form an intermediate slot. This slot favorably redirects airflow from the main element lower surface toward the low-momentum, upper-surface boundary layer of the flap. Relative to the supercritical wing, the cruise slotted wing enables greater aft loading that helps to reduce shock strength and transonic pressure drag. The cruise slotted wing may be considered a passive drag reduction technology with potential fuel-burn savings and increased vehicle range for next-generation aircraft. The current paper seeks to quantify the drag-saving benefits of cruise slotted wing technology in application to single-aisle commercial transport aircraft. A series of multielement design constraints were developed within CDISC for the design of a partial-span, cruise slotted wing for a Mach-0.8 variant of the Common Research Model. To mitigate the skin-friction drag penalty associated with cruise slotted wings, the design features a flap with natural laminar flow over approximately 75% of the flap surface area. Cruise drag estimates from the NASA USM3D-ME flow solver between a supercritical wing and the partial-span, cruise slotted wing design were within one drag count. Near-cruise, off-design analyses showed limited laminar-flow sensitivity to angle of attack and a more gradual drag rise compared to a supercritical wing. Based on this observed benefit, future work is motivated to design a cruise slotted wing with laminar flow on both wing elements to achieve a significant reduction in cruise drag while delaying drag divergence.

CDISC↗

NASA NTRSDate not supplied

Transonic Cruise Slotted Wing Design for Commercial Transport Aircraft using CDISC

A knowledge-based aerodynamic design method, CDISC, has been extended to enable the computational design of transonic cruise slotted wings for commercial transport aircraft. The cruise slotted wing is a multielement wing concept with a forward main element and an aft flap element, separated to form an intermediate slot. This slot favorably redirects airflow from the main element lower surface toward the low-momentum, upper-surface boundary layer of the flap. Relative to the supercritical wing, the cruise slotted wing enables greater aft loading that helps to reduce shock strength and transonic pressure drag. The cruise slotted wing may be considered a passive drag-reduction technology with potential fuel burn savings and increased vehicle range for next-generation aircraft. The current paper seeks to quantify the drag-saving benefits of cruise slotted wing technology in application to single-aisle commercial transport aircraft. A series of multielement design constraints were developed within CDISC for the design of a partial-span, cruise slotted wing for a Mach-0.8 variant of the Common Research Model. To mitigate the skin-friction drag penalty associated with cruise slotted wings, the design features a flap with natural laminar flow over approximately 75% of the flap surface area. Cruise drag estimates from the NASA USM3D flow solver between a supercritical wing and the partial-span, cruise slotted wing design were within one drag count. Near-cruise, off-design analyses showed limited laminar flow sensitivity to angle of attack and a more gradual drag rise compared to a supercritical wing. Based on this demonstrated benefit, future work is motivated to design a cruise slotted wing with natural laminar flow on both wing elements to achieve a significant reduction in cruise drag while providing delayed drag divergence.

CDISC↗

NASA NTRS1984

Multi-grid technique for the solution of incompressible Navier-Stokes equations

A computationally efficient method is proposed for obtaining fine-mesh solutions of the Navier-Stokes equations for high-Re flow, including separation. The method involves implementing a multi-grid solution procedure with suitably chosen elements, such that the resulting overall solution technique is efficient as well as robust. The robustness and efficiency of the solution technique are demonstrated by applying it to three model problems: flow in a driven cavity, downstream asymptotic flow in curved ducts of square and polar sections, and Newmann boundary-value problem in clustered curvilinear orthogonal coordinates.

Ghia, K. N.↗

NASA NTRS1992

Contact thermal shock test of ceramics

A novel quantitative thermal shock test of ceramics is described. The technique employs contact between a metal-cooling rod and hot disk-shaped specimen. In contrast with traditional techniques, the well-defined thermal boundary condition allows for accurate analyses of heat transfer, stress, and fracture. Uniform equibiaxial tensile stresses are induced in the center of the test specimen. Transient specimen temperature and acoustic emission are monitored continuously during the thermal stress cycle. The technique is demonstrated with soda-lime glass specimens. Experimental results are compared with theoretical predictions based on a finite-element method thermal stress analysis combined with a statistical model of fracture. Material strength parameters are determined using concentric ring flexure tests. Good agreement is found between experimental results and theoretical predictions of failure probability as a function of time and initial specimen temperature.

Rogers, W. P.↗

NASA NTRS1990

Approximate minimum-time trajectories for two-link flexible manipulators

The method of recursive quadratic programming was used to generate approximate minimum-time tip trajectories for two-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 bending only in the horizontal plane. 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 to compare a semi-rigid configuration to a flexible one. The level of bending flexibility necessary to excite significant modal behavior is demonstrated. Applied torques for minimum-time maneuvers are shown to be very similar between configurations and retain much of the qualitative character of rigid-body slewing motion.

Eisler, G. R.↗

NASA NTRS2014

Flutter Analysis for Turbomachinery Using Volterra Series

The objective of this paper is to describe an accurate and efficient reduced order modeling method for aeroelastic (AE) analysis and for determining the flutter boundary. Without losing accuracy, we develop a reduced order model based on the Volterra series to achieve significant savings in computational cost. The aerodynamic force is provided by a high-fidelity solution from the Reynolds-averaged Navier-Stokes (RANS) equations; the structural mode shapes are determined from the finite element analysis. The fluid-structure coupling is then modeled by the state-space formulation with the structural displacement as input and the aerodynamic force as output, which in turn acts as an external force to the aeroelastic displacement equation for providing the structural deformation. NASA's rotor 67 blade is used to study its aeroelastic characteristics under the designated operating condition. First, the CFD results are validated against measured data available for the steady state condition. Then, the accuracy of the developed reduced order model is compared with the full-order solutions. Finally the aeroelastic solutions of the blade are computed and a flutter boundary is identified, suggesting that the rotor, with the material property chosen for the study, is structurally stable at the operating condition, free of encountering flutter.

computational fluid dynamics↗

NASA NTRSDate not supplied

Supersonic Traveling Crossflow Wave Characteristics in Ground and Flight Tests

This paper continues analysis of data acquired in a series of supersonic crossflow transition experiments on a 65° swept- wing model tested in both wind tunnel and flight. Flow visualization data are summarized to provide context for the swept-wing instability mechanisms of interest. Unsteady pressure measurements obtained near the wing model leading edge (LE) are studied using 2nd- and 3rd-order statistical methods to extract traveling crossflow wave characteristics from these data. Comparisons between tunnel and flight data with boundary-layer (BL) Linear Stability Theory (LST) predictions provide insights into BL transition phenomena similarities and differences observed in tests using the same swept-wing model configurations. Comparisons of unsteady pressure data results for wing LE configurations with and without distributed-roughness-element (DRE) patterns installed, suggest traveling crossflow disturbance growth is reduced when this type of flow control is applied. A brief stud y of different background surface roughness levels on the wing LE suggests the potential for another flow control approach to promote laminar flow on swept wings using streamwise-biased surface finishes.

Lewis R Owens↗

DOE OSTI2024

A Simplified Method for Predicting Shaker Voltage in IMMATs

Impedance Matched Multi-Axis Tests (IMMATs) can replicate in-service vibration induced stress more accurately than single axis shaker table tests as they can better match a part’s operational boundary conditions and excite it in multiple degrees of freedom simultaneously. Here, the shakers used in IMMATs are less powerful than shaker tables, so shaker force limits can be exceeded during tests if they are not placed adequately for the desired environment. The ability to predict shaker voltage and force before performing a test is, therefore, helpful in selecting shaker locations so that their limits are not exceeded. In this study, electrodynamic shakers were modeled as discrete electromechanical systems, and the shaker parameters were chosen to match experimentally obtained acceleration/voltage frequency response functions (FRFs). These models were coupled to a finite element model of the device under test (DUT) via dynamic substructuring, and the substructured model was demonstrated to accurately predict shaker voltage as well as the error in reproducing the environment at multiple accelerometer locations. A simple method called the FRF Multiplication method, in which the FRF of the substructured system is approximated as the product of two separate FRFs of the shaker and DUT respectively, was proposed and applied to the same system, yielding similar voltage and error predictions to those obtained using substructuring. Simple case studies were presented to explore the applicability of the proposed method, and it was demonstrated to have similar accuracy to the substructuring method in a range of cases. Additionally, we showed that while it was not possible to derive a unique model of the shakers from acceleration/voltage FRFs alone, the models that could be obtained were sufficient to predict test error almost perfectly and shaker voltage with less than 40 percent error.

42 ENGINEERING↗

DOE OSTI2025

Tunable magnons of an antiferromagnetic Mott insulator via interfacial metal-insulator transitions

Antiferromagnetic insulators present a promising alternative to ferromagnets due to their ultrafast spin dynamics essential for low-energy terahertz spintronic device applications. Magnons, i.e., quantized spin waves capable of transmitting information through excitations, serve as a key functional element in this paradigm. However, identifying external mechanisms to effectively tune magnon properties has remained a major challenge. Here we demonstrate that interfacial metal-insulator transitions offer an effective method for controlling the magnons of Sr 2 IrO 4 , a strongly spin-orbit coupled antiferromagnetic Mott insulator. Resonant inelastic x-ray scattering experiments reveal a significant softening of zone-boundary magnon energies in Sr 2 IrO 4 films epitaxially interfaced with metallic 4 d transition-metal oxides. Therefore, the magnon dispersion of Sr 2 IrO 4 can be tuned by metal-insulator transitions of the 4 d transition-metal oxides. We tentatively attribute this non-trivial behavior to a long-range phenomenon mediated by magnon-acoustic phonon interactions. Our experimental findings introduce a strategy for controlling magnons and underscore the need for further theoretical studies to better understand the underlying microscopic interactions between magnons and phonons.

magnetic properties and materials↗

NASA NTRS1977

Structural analysis of ultra-high speed aircraft structural components

The buckling characteristics of a hypersonic beaded skin panel were investigated under pure compression with boundary conditions similar to those found in a wing mounted condition. The primary phases of analysis reported include: (1) experimental testing of the panel to failure; (2) finite element structural analysis of the beaded panel with the computer program NASTRAN; and (3) summary of the semiclassical buckling equations for the beaded panel under purely compressive loads. A comparison of each of the analysis methods is also included.

Lenzen, K. H.↗