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 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS1993

A unit-cell model of textile composite beams for predicting stiffness properties

Flexural stiffness properties of a textile composite beam are obtained from a finite-element model of the unit cell. Three linearly independent deformations, namely, pure extension, pure bending and pure shear, are applied to the unit cell. The top and bottom surfaces of the beam are assumed to be traction free. Periodic boundary conditions on the lateral boundaries of the unit cell are enforced by multi-point constraint elements. From the forces acting on the unit cell, the flexural stiffness coefficients of the composite beam are obtained. The difficulties in determining the transverse shear stiffness are discussed, and a modified approach is presented. The methods are first verified by applying them to isotropic and bimaterial beams for which the results are known, and then illustrated for a simple plain-weave textile composite.

Sankar, Bhavani V.↗

NASA NTRS1993

Elastic interactions of a fatigue crack with a micro-defect by the mixed boundary integral equation method

In this paper, the mixed boundary integral equation method is developed to study the elastic interactions of a fatigue crack and a micro-defect such as a void, a rigid inclusion or a transformation inclusion. The method of pseudo-tractions is employed to study the effect of a transformation inclusion. An enriched element which incorporates the mixed-mode stress intensity factors is applied to characterize the singularity at a moving crack tip. In order to evaluate the accuracy of the numerical procedure, the analysis of a crack emanating from a circular hole in a finite plate is performed and the results are compared with the available numerical solution. The effects of various micro-defects on the crack path and fatigue life are investigated. The results agree with the experimental observations.

Lua, Yuan J.↗

NASA NTRS1999

Parallel Aeroelastic Analysis Using ENSAERO and NASTRAN

A high fidelity parallel static structural analysis capability is created and interfaced to the multidisciplinary analysis package ENSAERO-MPI of Ames Research Center. This new module replaced ENSAERO's lower fidelity simple finite element and modal modules. Full aircraft structures may be more accurately modeled using the new finite element capability. Parallel computation is performed by breaking the full structure into multiple substructures. This approach is conceptually similar to ENSAERO's multi-zonal fluid analysis capability. The new substructure code is used to solve the structural finite element equations for each substructure in parallel. NASTRAN/COSMIC is utilized as a front end for this code. Its full library of elements can be used to create an accurate and realistic aircraft mode. It is used to create the stiffness matrices for each sub-structure. The new parallel code then uses an iterative preconditioned conjugate gradient method to solve the global structural equations for the sub-structure boundary nodes. Results are presented for a wing-body configuration.

Eldred, Lloyd B.↗

NASA NTRS1976

Convection cells in the early lunar magma ocean - Trace-element evidence

A method was proposed for recognizing suites of samples which are related to the differentiation of distinct and regionally separated magmas in the early lunar surface ocean. The sample suites were identified by three Cl-P2O5 compositional groups. New experimental data continue to add samples to these groups. Each magma was associated with a convection cell the boundaries of which are proposed. Independent geochemical evidence supporting a limited number of early magmas is cited. Two new possible Cl-P2O5 groups are identified. Samples not indigenous to a given site and their possible source regions are identified. It is proposed that many Apollo 17 basalts are not locally derived.

Jovanovic, S.↗

NASA NTRS2022

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.

Boundary Layer Stability↗

NASA NTRS1988

Bio-, Magneto- and event-stratigraphy across the K-T boundary

Determining the time and the time structure of rare events in geology can be accomplished by applying three different and independent stratigraphic methods: Biostratigraphy, magneto-stratigraphy and event-stratigraphy. The optimal time resolution of the two former methods is about 1000 years, while by means of event-stratigraphy a resolution of approximately one year can be achieved. For biostratigraphy across the Cretaceous-Tertiary (K-T) boundary micro- and nannofossils have been found best suited. The qualitative and quantitative analyses of minerals and trace elements across the K-T boundary show anomalies on a millimeter scale and permit conclusions regarding the time structure of the K-T event itself. The results of the analyses find a most consistent explanation by the assumption of an extraterrestrial impact. The main portion of the material rain from the atmosphere evidently was deposited within a short time. The long-time components consist of the finest portion of the material rain from the atmosphere and the transported and redeposited fall-out.

Preisinger, A.↗

NASA NTRS2002

Discrete Roughness Effects on Shuttle Orbiter at Mach 6

Discrete roughness boundary layer transition results on a Shuttle Orbiter model in the NASA Langley Research Center 20-Inch Mach 6 Air Tunnel have been reanalyzed with new boundary layer calculations to provide consistency for comparison to other published results. The experimental results were previously obtained utilizing the phosphor thermography system to monitor the status of the boundary layer via global heat transfer images of the Orbiter windward surface. The size and location of discrete roughness elements were systematically varied along the centerline of the 0.0075-scale model at an angle of attack of 40 deg and the boundary layer response recorded. Various correlative approaches were attempted, with the roughness transition correlations based on edge properties providing the most reliable results. When a consistent computational method is used to compute edge conditions, transition datasets for different configurations at several angles of attack have been shown to collapse to a well-behaved correlation.

Berry, Scott A.↗

NASA NTRSDate not supplied

Instabilities in the Wake of Roughness on a Flat Plate in a Quiet Supersonic Tunnel

Roughness-induced transition is an unavoidable reality in practical high-speed vehicles. Typical prediction of transition due to roughness include algebraic correlations and, more recently, semi-empirical methods. In the NASA Langley Research Center Supersonic Low Disturbance Tunnel, a Mach 3.5 quiet tunnel, several transition experiments have been performed in the past decade to better understand the mechanisms by which small roughness causes transition in a supersonic boundary layer. The study started first with isolated roughness elements of different planforms and shapes and progressed to increasingly more complicated geometries before arriving at a pseudorandom roughness, defined by an analytic function. Concurrent computational efforts progressed with these studies as well, starting with the use of linear stability theory and progressing to the use of harmonic linearized Navier-Stokes to predict the growth of boundary layer stabilities.

Amanda Chou↗

NASA NTRS1983

Moving finite elements in 2-D

The mathematical background regarding the moving finite element (MFE) method of Miller and Miller (1981) is discussed, taking into account a general system of partial differential equations (PDE) and the amenability of the MFE method in two dimensions to code modularization and to semiautomatic user-construction of numerous PDE systems for both Dirichlet and zero-Neumann boundary conditions. A description of test problem results is presented, giving attention to aspects of single square wave propagation, and a solution of the heat equation.

Gelinas, R. J.↗

NASA NTRS2011

Robust and Accurate Shock Capturing Method for High-Order Discontinuous Galerkin Methods

A simple yet robust and accurate approach for capturing shock waves using a high-order discontinuous Galerkin (DG) method is presented. The method uses the physical viscous terms of the Navier-Stokes equations as suggested by others; however, the proposed formulation of the numerical viscosity is continuous and compact by construction, and does not require the solution of an auxiliary diffusion equation. This work also presents two analyses that guided the formulation of the numerical viscosity and certain aspects of the DG implementation. A local eigenvalue analysis of the DG discretization applied to a shock containing element is used to evaluate the robustness of several Riemann flux functions, and to evaluate algorithm choices that exist within the underlying DG discretization. A second analysis examines exact solutions to the DG discretization in a shock containing element, and identifies a "model" instability that will inevitably arise when solving the Euler equations using the DG method. This analysis identifies the minimum viscosity required for stability. The shock capturing method is demonstrated for high-speed flow over an inviscid cylinder and for an unsteady disturbance in a hypersonic boundary layer. Numerical tests are presented that evaluate several aspects of the shock detection terms. The sensitivity of the results to model parameters is examined with grid and order refinement studies.

Atkins, Harold L.↗

NASA NTRS2015

Tetrahedral-Mesh Simulation of Turbulent Flows with the Space-Time Conservative Schemes

Direct numerical simulations of turbulent flows are predominantly carried out using structured, hexahedral meshes despite decades of development in unstructured mesh methods. Tetrahedral meshes offer ease of mesh generation around complex geometries and the potential of an orientation free grid that would provide un-biased small-scale dissipation and more accurate intermediate scale solutions. However, due to the lack of consistent multi-dimensional numerical formulations in conventional schemes for triangular and tetrahedral meshes at the cell interfaces, numerical issues exist when flow discontinuities or stagnation regions are present. The space-time conservative conservation element solution element (CESE) method - due to its Riemann-solver-free shock capturing capabilities, non-dissipative baseline schemes, and flux conservation in time as well as space - has the potential to more accurately simulate turbulent flows using unstructured tetrahedral meshes. To pave the way towards accurate simulation of shock/turbulent boundary-layer interaction, a series of wave and shock interaction benchmark problems that increase in complexity, are computed in this paper with triangular/tetrahedral meshes. Preliminary computations for the normal shock/turbulence interactions are carried out with a relatively coarse mesh, by direct numerical simulations standards, in order to assess other effects such as boundary conditions and the necessity of a buffer domain. The results indicate that qualitative agreement with previous studies can be obtained for flows where, strong shocks co-exist along with unsteady waves that display a broad range of scales, with a relatively compact computational domain and less stringent requirements for grid clustering near the shock. With the space-time conservation properties, stable solutions without any spurious wave reflections can be obtained without a need for buffer domains near the outflow/farfield boundaries. Computational results for the isotropic turbulent flow decay, at a relatively high turbulent Mach number, show a nicely behaved spectral decay rate for medium to high wave numbers. The high-order CESE schemes offer very robust solutions even with the presence of strong shocks or widespread shocklets. The explicit formulation in conjunction with a close to unity theoretical upper Courant number bound has the potential to offer an efficient numerical framework for general compressible turbulent flow simulations with unstructured meshes.

Chang, Chau-Lyan↗

NASA NTRS1979

Multi-element airfoil viscous-inviscid interactions

Subsonic viscous-inviscid interactions for multi-element airfoils are predicted by iterating between inviscid and viscous solutions until the performance coefficients converge. Inviscid flow is modelled by using distributed source-vortex singularities on configuration surface panels. Viscous effects are calculated by an existing laminar separation bubble model and a NASA-Lockheed boundary layer-wake method. Numerical formulations and example calculations are presented.

Gross, L. W.↗

NASA NTRS1983

Focal axis resolver for offset reflector antennas

Method and apparatus for determining the focal axis of an asymmetrical antenna such as an offset paraboloid reflector whose physical rim is not coincident with the boundary of the electrical aperture but whose focal point is known is provided. A transmitting feed horn array consisting of at least two feed horn elements is positioned asymmetrically on either side of an estimated focal axis which is generally inclined with respect to the boresight axis of the antenna. The feed horn array is aligned with the estimated focal axis so that the phase centers (CP sub 1, CP sub 2) of the two feed horn elements are located on a common line running through the focal point (F) orthogonally with respect to the estimated focal axis.

Schmidt, R. F.↗

NASA NTRS1985

Three-dimensional hybrid-stress finite element analysis of composite laminates with cracks and cutouts

A three-dimensional hybrid-stress finite element analysis of composite laminates containing cutouts and cracks is presented. Fully three-dimensional, hexahedral isoparametric elements of the hybrid-stress model are formulated on the basis of the Hellinger-Reissner variational principle. Traction-free edges, cutouts, and crack surfaces are modeled by imposition of exact traction boundary conditions along element surfaces. Special boundary and surface elements are constructed by introducing proper constraints on assumed stress functions. The Lagrangian multiplier technique is used to enforce ply-interface continuity conditions in hybrid bimaterial composite elements for modeling the interface region in a composite laminate. Two examples are given to illustrate the capability of the present method of approach: (1) the well-known delamination problem in an angle-ply laminate, and (2) the important problem of a composite laminate containing a circular hole. Results are presented in detail for each case. Implications of interlaminar and intralaminar crack initiation, growth and fracture in composites containing cracks and cutouts are discussed.

Wang, S. S.↗

NASA NTRS1987

Design of Software for Design of Finite Element for Structural Analysis

The concepts of software engineering which allow a user of the finite element method to describe a model, to collect and to check the model data in a data base as well as to form the matrices required for a finite element calculation are examined. Next the components of the model description are conceived including the mesh tree, the topology, the configuration, the kinematic boundary conditions, the data for each element, and the loads. The possibilities for description and review of the data are considered. The concept of the segments for the modularization of the programs follows the components of the model description. The significance of the mesh tree as a globular guiding structure will be understood in view of the principle of the unity of the model, the mesh tree, and the data base. The user-friendly aspects of the software system will be summarized: the principle of language communication, the data generators, error processing, and data security.

Helfrich, Reinhard↗

NASA NTRS1992

Coupled finite-difference/finite-element approach for wing-body aeroelasticity

Computational methods using finite-difference approaches for fluids and finite-element approaches for structures have individually advanced to solve almost full-aircraft configurations. However, coupled approaches to solve fluid/structural interaction problems are still in their early stages of development, particularly for complex geometries using complete equations such as the Euler/Navier-Stokes equations. Earlier work demonstrated the success of coupling finite-difference and finite-element methods for simple wing configurations using the Euler/Navier-Stokes equations. In this paper, the same approach is extended for general wing-body configurations. The structural properties are represented by beam-type finite elements. The flow is modeled using the Euler/Navier-Stokes equations. A general procedure to fully couple structural finite-element boundary conditions with fluid finite-difference boundary conditions is developed for wing-body configurations. Computations are made using moving grids that adapt to wing-body structural deformations. Results are illustrated for a typical wing-body configuration.

Guruswamy, Guru P.↗

NASA NTRS2007

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗

NASA NTRS2018

Unstructured Grid Adaptation and Solver Technology for Turbulent Flows

Unstructured grid adaptation is a tool to control Computational Fluid Dynamics (CFD) discretization error. However, adaptive grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Issues that prevent the use of adaptive grid methods are identified by applying unstructured grid adaptation methods to a series of benchmark cases. Once identified, these challenges to existing adaptive workflows can be addressed. Unstructured grid adaptation is evaluated for test cases described on the Turbulence Modeling Resource (TMR) web site, which documents uniform grid refinement of multiple schemes. The cases are turbulent flow over a Hemisphere Cylinder and an ONERA M6Wing. Adaptive grid force and moment trajectories are shown for three integrated grid adaptation processes with Mach interpolation control and output error based metrics. The integrated grid adaptation process with a finite element (FE) discretization produced results consistent with uniform grid refinement of fixed grids. The integrated grid adaptation processes with finite volume schemes were slower to converge to the reference solution than the FE method. Metric conformity is documented on grid/metric snapshots for five grid adaptation mechanics implementations. These tools produce anisotropic boundary conforming grids requested by the adaptation process.

Park, Michael A.↗