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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Transonic flow computations in cascades using finite element method

Analysis of transonic flow through a cascade of airfoils is investigated using the finite element method. Development of a computational grid suitable for complex flow structures and different types of boundary conditions is presented. An efficient pseudo-time integration scheme is developed for the solution of equations. Modeling of the shock and the convergence characteristics of the developed scheme are discussed. Numerical results include a 45 deg staggered cascade of NACA 0012 airfoils with inlet flow Mach number of 0.8 and angles of attack 1, 0, and 1 deg.

Akay, H. U.↗

NASA NTRS1985

Nonlinear quasi-fracture behavior of polymers

This paper deals with the nonlinear quasi-fracture behavior of polymers. Using the finite element method, the quasi-fracture model composed of an isolated craze opening with a nonlinear boundary envelope forces represented by oriented fibril domains is analyzed in simply-stressed rectangular plate. The stress distribution, the displacement profile, as well as the yield or plastic or flow zone adjacent to the craze tips are obtained. It is found that the variation of the stress distribution near the craze tip is not as sharp as that for an elastic crack. The initial load required for yielding or flowing is greater for the quasifracture as compared with that for a classical fracture. Under the same load the yield or flow zone is somewhat smaller as compared with that obtained on the basis of linear elastic-plastic fracture mechanics theory.

Sun, B. N.↗

NASA NTRS1987

Calculation of stresses and forces between the skin and stiffener in composite panels

A methodology for computing the stresses at the interface between the skin and stiffener in stiffened composite panels is described. The methodology is based on finite-element analyses and on elasticity solution. The finite-element analyses are standard, while the elasticity solution is based on eigenvalue expansions of the stress function. The eigenvalue expansions are assumed to be valid in the local region where the stiffener flange terminates. The local elasticity solution is coupled to the global finite-element analysis using collocation on the boundary of the local region. Accuracy and convergence of the method are discussed and several examples of its utility are presented.

Hyer, M. W.↗

NASA NTRS2007

Compatibility Condition in Theory of Solid Mechanics (Elasticity, Structures, and Design Optimization)

The strain formulation in elasticity and the compatibility condition in structural mechanics have neither been understood nor have they been utilized. This shortcoming prevented the formulation of a direct method to calculate stress. We have researched and understood the compatibility condition for linear problems in elasticity and in finite element analysis. This has lead to the completion of the method of force with stress (or stress resultant) as the primary unknown. The method in elasticity is referred to as the completed Beltrami-Michell formulation (CBMF), and it is the integrated force method (IFM) in structures. The dual integrated force method (IFMD) with displacement as the primary unknown has been formulated. IFM and IFMD produce identical responses. The variational derivation of the CBMF yielded the new boundary compatibility conditions. The CBMF can be used to solve stress, displacement, and mixed boundary value problems. The IFM in structures produced high-fidelity response even with a modest finite element model. The IFM has influenced structural design considerably. A fully utilized design method for strength and stiffness limitation has been developed. The singularity condition in optimization has been identified. The CBMF and IFM tensorial approaches are robust formulations because of simultaneous emphasis on the equilibrium equation and the compatibility condition.

Patnaik, Surya N.↗

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 NTRS1999

A High Order Discontinuous Galerkin Method for 2D Incompressible Flows

In this paper we introduce a high order discontinuous Galerkin method for two dimensional incompressible flow in vorticity streamfunction formulation. The momentum equation is treated explicitly, utilizing the efficiency of the discontinuous Galerkin method The streamfunction is obtained by a standard Poisson solver using continuous finite elements. There is a natural matching between these two finite element spaces, since the normal component of the velocity field is continuous across element boundaries. This allows for a correct upwinding gluing in the discontinuous Galerkin framework, while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability The method is suitable for inviscid or high Reynolds number flows. Optimal error estimates are proven and verified by numerical experiments.

Liu, Jia-Guo↗

NASA NTRS2004

Boundary Layer Control for Hypersonic Airbreathing Vehicles

Active and passive methods for tripping hypersonic boundary layers have been examined in NASA Langley Research Center wind tunnels using a Hyper-X model. This investigation assessed several concepts for forcing transition, including passive discrete roughness elements and active mass addition (or blowing), in the 20-Inch Mach 6 Air and the 31-Inch Mach 10 Air Tunnels. Heat transfer distributions obtained via phosphor thermography, shock system details, and surface streamline patterns were measured on a 0.333-scale model of the Hyper-X forebody. The comparisons between the active and passive methods for boundary layer control were conducted at test conditions that nearly match the Hyper-X nominal Mach 7 flight test-point of an angle-of-attack of 2-deg and length Reynolds number of 5.6 million. For passive roughness, the primary parametric variation was a range of trip heights within the calculated boundary layer thickness for several trip concepts. The passive roughness study resulted in a swept ramp configuration, scaled to be roughly 0.6 of the calculated boundary layer thickness, being selected for the Mach 7 flight vehicle. For the active blowing study, the manifold pressure was systematically varied (while monitoring the mass flow) for each configuration to determine the jet penetration height, with schlieren, and transition movement, with the phosphor system, for comparison to the passive results. All the blowing concepts tested, which included various rows of sonic orifices (holes), two- and three-dimensional slots, and random porosity, provided transition onset near the trip location with manifold stagnation pressures on the order of 40 times the model surface static pressure, which is adequate to ensure sonic jets. The present results indicate that the jet penetration height for blowing was roughly half the height required with passive roughness elements for an equivalent amount of transition movement.

Berry, Scott A.↗

NASA NTRS2004

A Critical Review of the Transport and Decay of Wake Vortices in Ground Effect

This slide presentation reviews the transport and decay of wake vortices in ground effect and cites a need for a physics-based parametric model. The encounter of a vortex with a solid body is always a complex event involving turbulence enhancement, unsteadiness, and very large gradients of velocity and pressure. Wake counter in ground effect is the most dangerous of them all. The interaction of diverging, area-varying, and decaying aircraft wake vortices with the ground is very complex because both the vortices and the flow field generated by them are altered to accommodate the presence of the ground (where there is very little room to maneuver) and the background turbulent flow. Previous research regarding vortex models, wake vortex decay mechanisms, time evolution within in ground effect of a wake vortex pair, laminar flow in ground effect, and the interaction of the existing boundary layer with a convected vortex are reviewed. Additionally, numerical simulations, 3-dimensional large-eddy simulations, a probabilistic 2-phase wake vortex decay and transport model and a vortex element method are discussed. The devising of physics-based, parametric models for the prediction of (operational) real-time response, mindful of the highly three-dimensional and unsteady structure of vortices, boundary layers, atmospheric thermodynamics, and weather convective phenomena is required. In creating a model, LES and field data will be the most powerful tools.

Sarpkaya, T.↗

NASA NTRS1978

Two-dimensional separated wake modeling and its use to predict maximum section lift coefficient

A technique for computing the lift of separating multielement airfoils in incompressible flow is presented. The procedure employs repeated application of a panel method to solve for the separated wake displacement surface using entirely inviscid boundary conditions. Results are presented that compare computed pressure distributions with those measured in the wind tunnel for airfoils with one, two, and four elements with separation on each element. A method employing this technique is presented which shows promise in predicting airfoil section lift through stall.

Henderson, M. L.↗

NASA NTRS1982

Noise transmission into semicylindrical enclosures through discretely stiffened curved panels

An analytical study of sound transmission into semicylindrical enclosures through discretely stiffened curved elastic panels is presented. The transmitted sound is estimated by solving the acoustic wave equation for the interior acoustic field, a Galerkin-like method being used. This solution is then coupled to the vibration of the stiffened panels. The response characteristics of these panels are determined by using a modal analysis where the modes are obtained by the finite element-strip method. Numerical results include spectra of the interior sound pressure due to white noise, turbulent boundary layer and propeller noise inputs.

Chang, M. T.↗

NASA NTRS1990

A new finite element method for solving compressible Navier-Stokes equations based on an operator splitting method and h-p adaptivity

A new finite element method solving compressible Navier-Stokes equations is proposed. The method is based on a version of Strang's operator splitting and an h-p adaptive finite element approximation in space. This paper contains the formulation of the method with a detailed discussion of boundary conditions, a sample adaptive strategy and numerical examples involving compressible viscous flow over a flat plate with Reynolds number Re = 1000 and Re = 10,000.

Demkowicz, L.↗

NASA NTRS1991

A geometrically nonlinear analysis of interlaminar stresses in unsymmetrically laminated plates subjected to uniform thermal loading

An analytical study of interlaminar stresses in unsymmetrically laminated plates is presented. The study examines the linear elastic, large deflection response of square laminated composite plates subjected to uniform thermal loading. Both cross-ply and angle-ply, symmetric and unsymmetric, laminates are examined to evaluate the effects of mismatch between adjacent layers in elastic constants and coefficients of thermal expansion. A geometrically nonlinear kinematic description is used to predict the large out-of-plane (transverse) deflections. The nonlinear, three-dimensional boundary value problems are formulated from elasticity theory and approximate solutions are determined using the finite element method. A global/local analysis procedure is used to obtain improved free edge stress predictions. For the laminates and loading conditions considered, the results indicate that the out-of-plane deflections of the unsymmetric laminates reduce interlaminar shear stresses. These deflections also reduce interlaminar normal stresses in some laminates and increase these stresses for others. The results indicate that as the out-of-plane deflections become large, the differences in interlaminar stress predictions between linear and nonlinear theory can become quite large.

Norwood, D. S.↗

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.↗

NASA NTRS2003

Electromagnetic Interactions GenERalized (EIGER): Applications at Johnson Space Center

This slide presentation reviews the software framework for the analysis and design of complex electromagnetic systems. The system called Electromagnetic Interactions GenERalized (EIGER) is a hybrid finite element solution for wave equations and deru=ives boundary element solutions for integral equation formulations. Written in Fortran 90 Eiger uses object-oriented design methods to abstract key analysis components. EIGER has been applied to higher order modeling for analysis of antennas to assist in the design of the antennas of a mini-AERCam. EIGER might also be used to model the Ion Cyclotron Resonant heating stage for the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR).

Dobbins, Justin A.↗

NASA NTRS1974

Boundary value problems with incremental plasticity in granular media

Discussion of the critical state concept in terms of an incremental theory of plasticity in granular (soil) media, and formulation of the governing equations which are convenient for a computational scheme using the finite element method. It is shown that the critical state concept with its representation by the classical incremental theory of plasticity can provide a powerful means for solving a wide variety of boundary value problems in soil media.

Chung, T. J.↗

NASA NTRS1993

Multigrid methods and high order finite difference for flow in transition - Effects of isolated and distributed roughness elements

The high order finite difference and multigrid methods have been successfully applied to direct numerical simulation (DNS) for flow transition in 3D channels and 3D boundary layers with 2D and 3D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semicoarsening multigrid method associated with line distributive relaxation scheme, and a new treatment of the outflow boundary condition, which needs only a very short buffer domain to damp all wave reflection, are developed. These approaches make the multigrid DNS code very accurate and efficient. This makes us not only able to do spatial DNS for the 3D channel and flat plate at low computational costs, but also able to do spatial DNS for transition in the 3D boundary layer with 3D single and multiple roughness elements. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments.

Liu, C.↗

NASA NTRS2019

Investigation of Transitional Shock-Wave/Boundary Layer Interactions Using Direct Numerical Simulations

Interaction of a transitional boundary layer with a normal-shock is investigated using unstructured tetrahedral meshes under the numerical framework of the space-time conservation element, solution element (CESE) method. The computations mimic recent experimental efforts at the University of Tennessee Space Institute, where a Mach 2.0 flow interacts with a tall cylinder attached to a flat plate. The location of the cylinder with respect to the flat plate leading edge determines if the incoming boundary layer is laminar, transitional or fully turbulent. Four representative flow conditions exemplifying laminar and transitional boundary layers are analyzed by direct numerical simulations. Similar to what was observed in the experiments for the case of transitional interaction, the computations reveal an intermittent upstream influence (UI) shock that repeatedly travels upstream from the lambda-foot toward the leading edge before vanishing. Through detailed unsteady flow analysis obtained using Fourier analysis and dynamic mode decomposition techniques, the presence of disturbances with similar frequencies as those measured in experiments were identified in the flow along with locations that appear to influence the dynamics of the flow.

Venkatachari, Balaji Shankar↗

NASA NTRS1985

Advances and trends in structures and dynamics; Proceedings of the Symposium, Washington, DC, October 22-25, 1984

Among the topics discussed are developments in structural engineering hardware and software, computation for fracture mechanics, trends in numerical analysis and parallel algorithms, mechanics of materials, advances in finite element methods, composite materials and structures, determinations of random motion and dynamic response, optimization theory, automotive tire modeling methods and contact problems, the damping and control of aircraft structures, and advanced structural applications. Specific topics covered include structural design expert systems, the evaluation of finite element system architectures, systolic arrays for finite element analyses, nonlinear finite element computations, hierarchical boundary elements, adaptive substructuring techniques in elastoplastic finite element analyses, automatic tracking of crack propagation, a theory of rate-dependent plasticity, the torsional stability of nonlinear eccentric structures, a computation method for fluid-structure interaction, the seismic analysis of three-dimensional soil-structure interaction, a stress analysis for a composite sandwich panel, toughness criterion identification for unidirectional composite laminates, the modeling of submerged cable dynamics, and damping synthesis for flexible spacecraft structures.

Noor, A. K.↗