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 721 records · Page 40

Boundary estimation problems arising in thermal tomography

Problems on the identification of two-dimensional spatial domains arising in the detection and characterization of structural flaws in materials are considered. For a thermal diffusion system with external boundary input, observations of the temperature on the surface are used in a output least squares approach. Parameter estimation techniques based on the method of mappings are discussed and approximation schemes are developed based on a finite element Galerkin approach. Theoretical convergence results for computational techniques are given and the results are applied to experimental data for the identification of flaws in the thermal testing of materials.

Banks, H. T.↗

Three-dimensional finite element progressive failure analysis of composite laminates under axial extension

A three-dimensional (3D) progressive failure algorithm is developed, where the layerwise laminate theory (LWLT) of Reddy is used for kinematic description. The finite element model based on the layerwise theory predicts both inplane and interlaminar stresses with the same accuracy as that of a conventional 3D finite element model and provides a convenient format for modeling the 3D stress fields in composite laminates. A parametric study is conducted to investigate the effect of out-of-plane material properties, 3D stiffness reduction methods, and boundary conditions on the failure loads and strains of a composite laminate under axial extension. The results indicate that different parameters have a different degree of influence on the failure loads and strains. The predictive ability of various phenomenological failure criteria is evaluated in the light of experimental results available in the literature, and the predictions of the LWLT are compared with those of the first-order shear deformation theory. It is concluded that a 3D stress analysis is necessary to predict accurately the failure behavior of composite laminates.

Reddy, Yeruva S.↗

Elasto-Plastic Analysis of Tee Joints Using HOT-SMAC

The Higher Order Theory - Structural/Micro Analysis Code (HOT-SMAC) software package is applied to analyze the linearly elastic and elasto-plastic response of adhesively bonded tee joints. Joints of this type are finding an increasing number of applications with the increased use of composite materials within advanced aerospace vehicles, and improved tools for the design and analysis of these joints are needed. The linearly elastic results of the code are validated vs. finite element analysis results from the literature under different loading and boundary conditions, and new results are generated to investigate the inelastic behavior of the tee joint. The comparison with the finite element results indicates that HOT-SMAC is an efficient and accurate alternative to the finite element method and has a great deal of potential as an analysis tool for a wide range of bonded joints.

Arnold, Steve M.↗

A Wall-Modeled LES Perspective for the High Lift Common Research Model Using LAVA

A new immersed boundary Wall-Modelled Large Eddy Simulation (WMLES) formulation is developed to study high-lift aerodynamics on the NASA High-Lift Common Research Model (HL-CRM). A sequence of Cartesian Octree grids with sizes ranging from 100 Million through 2.02 Billion grid points is utilized to systematically assess grid-sensitivity and convergence for the in-tunnel (QinetiQ) configuration of the model, and remarkable agreement between the immersed boundary and the curvilinear body-aligned WMLES formulations is reported on grids with comparable resolutions. In the free-air configuration of the model, consistent predictions between the Curvilinear Overset and the Cartesian Octree formulations are reported for angles of attack up to C(L,max) at a=19.57. However, some differences in the onset of stall are seen between the two methods for a>20°: while the curvilinear WMLES experiences wing-root separation with increasing angle of attack (Topology A), the Cartesian Octree formulation shows a different flow topology characterized by boundary layer weakness on the main element, emanating from the pylon-wing attachment (Topology B). In order to obtain further insight into the two-distinct topologies, carefully designed numerical experiments to isolate effects of the model standoff and the tunnel wall-boundary layers are conducted using the immersed boundary WMLES formulation. The increased incidence angle-of-attack on the inboard portion of the wing due to the standoff is shown to be sufficient for triggering a switch from Topology B to Topology A in Cartesian WMLES. The role of the floor boundary layer is further examined in detail by identification of additional corner-flow vorticity generated by the viscous juncture flow interactions between the floor boundary layer and the standoff leading to formation of a strong coherent and persistent vortex on the belly-side of the fuselage. The intensity of this vortex is shown to increase with the thickness of the floor boundary layer. A further increase in the incidence angle of attack near the leading-edge strake caused by the presence of this belly-side vortex is quantified for two-distinct floor boundary layers. Both of the floor boundary layers considered result in the onset of large scale wing-root separation at a=21.47in non-confined (free-air) configurations.

TTT↗

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Here, in this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial degree for analytic and cross-code test cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical simulation of hyperbolic heat conduction with convection boundary conditions and pulse heating effects

The paper describes the numerical simulation of hyperbolic heat conduction with convection boundary conditions. The effects of a step heat loading, a sudden pulse heat loading, and an internal heat source are considered in conjunction with convection boundary conditions. Two methods of solution are presened for predicting the transient behavior of the propagating thermal disturbances. In the first method, MacCormack's predictor-corrector method is employed for integrating the hyperbolic system of equations. Next, the transfinite element method, which employs specially tailored elements, is used for accurately representing the transient response of the propagating thermal wave fronts. The agreement between the results of various numerical test cases validate the representative behavior of the thermal wave fronts. Both methods represent hyperbolic heat conduction behavior by effectively modeling the sharp discontinuities of the propagating thermal disturbances.

Glass, David E.↗

Aeroelasticity of wing and wing-body configurations on parallel computers

The objective of this research is to develop computationally efficient methods for solving aeroelasticity problems on parallel computers. Both uncoupled and coupled methods are studied in this research. For the uncoupled approach, the conventional U-g method is used to determine the flutter boundary. The generalized aerodynamic forces required are obtained by the pulse transfer-function analysis method. For the coupled approach, the fluid-structure interaction is obtained by directly coupling finite difference Euler/Navier-Stokes equations for fluids and finite element dynamics equations for structures. This capability will significantly impact many aerospace projects of national importance such as Advanced Subsonic Civil Transport (ASCT), where the structural stability margin becomes very critical at the transonic region. This research effort will have direct impact on the High Performance Computing and Communication (HPCC) Program of NASA in the area of parallel computing.

Byun, Chansup↗

Flutter Analysis of the Transonic Truss-Braced Wing Aircraft Using Transonic Correction

This paper describes a flutter analysis method for the Transonic Truss-Braced Wing aircraft using a vortex lattice method coupled to an unsteady transonic correction method to account for unsteady aerodynamics in transonic flow. A steady-state vortex-lattice model of the Truss-Braced Wing aircraft is developed using vortex lattice code VSPAERO. A transonic and viscous flow correction method is implemented in the VSPAERO model to account for steady-state transonic and viscous flow effects using transonic small disturbance code TSFOIL coupled to an in-house integral boundary layer code. In addition, a wing-strut interference correction method is developed to account for the transonic interference aerodynamics in the strut juncture region using high-fidelity CFD code FUN3D. A structural dynamic finite-element model of the Truss-Braced Wing aircraft is developed using BEAM3D in-house finite-element code and is coupled to the VSPAERO. The BEAM3D model includes a geometric nonlinearity due to the tension in the strut which causes a deflection-dependent nonlinear stiffness. An unsteady transonic correction method is developed to better capture the unsteady aerodynamics in transonic flow. The unsteady transonic correction method makes use of the Theodorsen’s theory to account for the amplitude and phase shift of the unsteady lift coefficient in transonic flow. A preliminary flutter analysis of the Truss-Braced Wing aircraft is conduct to illustrate the unsteady transonic correction approach.

Nhan Nguyen↗

Incompressible Navier-Stokes and parabolized Navier-Stokes solution procedures and computational techniques

Recent developments with finite-difference techniques are emphasized. The quotation marks reflect the fact that any finite discretization procedure can be included in this category. Many so-called finite element collocation and galerkin methods can be reproduced by appropriate forms of the differential equations and discretization formulas. Many of the difficulties encountered in early Navier-Stokes calculations were inherent not only in the choice of the different equations (accuracy), but also in the method of solution or choice of algorithm (convergence and stability, in the manner in which the dependent variables or discretized equations are related (coupling), in the manner that boundary conditions are applied, in the manner that the coordinate mesh is specified (grid generation), and finally, in recognizing that for many high Reynolds number flows not all contributions to the Navier-Stokes equations are necessarily of equal importance (parabolization, preferred direction, pressure interaction, asymptotic and mathematical character). It is these elements that are reviewed. Several Navier-Stokes and parabolized Navier-Stokes formulations are also presented.

Rubin, S. G.↗

Modeling and Analysis of Wrinkled Membranes: An Overview

Thin-film membranes are basic elements of a variety of space inflatable/deployable structures. Wrinkling degrades the performance and reliability of these membrane structures, and hence has been a topic of continued interest. Wrinkling analysis of membranes for general geometry and arbitrary boundary conditions is quite challenging. The objective of this presentation is two-fold. Firstly, the existing models of wrinkled membranes and related numerical solution methods are reviewed. The important issues to be discussed are the capability of a membrane model to characterize taut, wrinkled and slack states of membranes in a consistent and physically reasonable manner; the ability of a wrinkling analysis method to predict the formation and growth of wrinkled regions, and to determine out-of-plane deformation and wrinkled waves; the convergence of a numerical solution method for wrinkling analysis; and the compatibility of a wrinkling analysis with general-purpose finite element codes. According to this review, several opening issues in modeling and analysis of wrinkled membranes that are to be addressed in future research are summarized, The second objective of this presentation is to discuss a newly developed membrane model of two viable parameters (2-VP model) and associated parametric finite element method (PFEM) for wrinkling analysis are introduced. The innovations and advantages of the proposed membrane model and PFEM-based wrinkling analysis are: (1) Via a unified stress-strain relation; the 2-VP model treat the taut, wrinkled, and slack states of membranes consistently; (2) The PFEM-based wrinkling analysis has guaranteed convergence; (3) The 2-VP model along with PFEM is capable of predicting membrane out-of-plane deformations; and (4) The PFEM can be integrated into any existing finite element code. Preliminary numerical examples are also included in this presentation to demonstrate the 2-VP model and PFEM-based wrinkling analysis approach.

Yang, B.↗

The strip method for shape derivatives

Abstract A major challenge in shape optimization is the coupling of finite element method (FEM) codes in a way that facilitates efficient computation of shape derivatives. This is particularly difficult with multiphysics problems involving legacy codes, where the costs of implementing and maintaining shape derivative capabilities are prohibitive. The volume and boundary methods are two approaches to computing shape derivatives. Each has a major drawback: the boundary method is less accurate, while the volume method is more invasive to the FEM code. We introduce the strip method , which computes shape derivatives on a strip adjacent to the boundary. The strip method makes code coupling simple. Like the boundary method, it queries the state and adjoint solutions at quadrature nodes, but requires no knowledge of the FEM code implementations. At the same time, it exhibits the higher accuracy of the volume method. As an added benefit, its computational complexity is comparable to that of the boundary method, that is, it is faster than the volume method. We illustrate the benefits of the strip method with numerical examples.

Hardesty, Sean↗

Critical study of higher order numerical methods for solving the boundary-layer equations

A fourth order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method, which is the natural extension of the second order box scheme to fourth order, was demonstrated with application to the incompressible, laminar and turbulent, boundary layer equations. The efficiency of the present method is compared with two point and three point higher order methods, namely, the Keller box scheme with Richardson extrapolation, the method of deferred corrections, a three point spline method, and a modified finite element method. For equivalent accuracy, numerical results show the present method to be more efficient than higher order methods for both laminar and turbulent flows.

Wornom, S. F.↗

An automatic multigrid method for the solution of sparse linear systems

An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.

Shapira, Yair↗

Uniaxial Tensile Properties of AS4 3D Woven Composites with Four Different Resin Systems: Experimental Results and Analysis: Property Computations

As a part of the NASA Composite Technology for Exploration project, eight different AS4 3D orthogonal woven composite panels were manufactured and were subjected to mechanical testing including uniaxial tension along the weaves' warp direction. Each set, with four different resin systems (KCR-IR6070, EP2400, RTM6, and RS-50), included weave architectures designed using 12K and 6K AS4 carbon fiber yarns. For the tension testing conducted at Room Temperature Ambient (RTA) conditions, the elastic modulus and strength of these eight panels (as-processed and thermally-cycled) were measured and compared while the potential evolution of micro-cracking before and after thermal cycling were monitored via optical microscopy and X-Ray Computed Tomography. The data set also included test results of the as-processed materials at Elevated Temperature Wet (ETW) conditions. In the second part of this study, efforts were made to compute elastic constants for AS4 6K/RTM6 and AS4 12K/RTM6 materials by implementing a finite element approach and the Multiscale Generalized Method of Cells (MSGMC) technique developed at NASA Glenn Research Center. Digimat-FE was used to model the weave architectures, assign properties, calculate yarn properties, create the finite element mesh, and compute the elastic properties by applying periodic boundary conditions to finite element models of each repeating unit cell. The required input data for MSGMC was generated using Matlab® from Digimat exported weave information. Experimental and computational results were compared, and the differences and limitations in correlating to the test data were briefly discussed.

Property Computations↗

Uniaxial Tensile Properties of AS4 3D Woven Composites with Four Different Resin Systems: Experimental Results & Analysis - Property Computations

As a part of the NASA Composite Technology for Exploration project, eight different AS4 3D orthogonal woven composite panels were manufactured and were subjected to mechanical testing including uniaxial tension along the weaves' warp direction. Each set, with four different resin systems (KCR-IR6070, EP2400, RTM6, and RS-50), included weave architectures designed using 12K and 6K AS4 carbon fiber yarns. For the tension testing conducted at Room Temperature Ambient (RTA) conditions, the elastic modulus and strength of these eight panels (as-processed and thermally cycled) were measured and compared while the potential evolution of micro-cracking before and after thermal cycling were monitored via optical microscopy and X-Ray Computed Tomography. The data set also included test results of the as-processed materials at Elevated Temperature Wet (ETW) conditions. In the second part of this study, efforts were made to compute elastic constants for AS4 6K/RTM6 and AS4 12K/RTM6 materials by implementing a finite element approach and the Multiscale Generalized Method of Cells (MSGMC) technique developed at NASA Glenn Research Center. Digimat-FE was used to model the weave architectures, assign properties, calculate yarn properties, create the finite element mesh, and compute the elastic properties by applying periodic boundary conditions to finite element models of each repeating unit cell. The required input data for MSGMC was generated using Matlab® from Digimat exported weave information. Experimental and computational results were compared, and the differences and limitations in correlating to the test data were briefly discussed.

Property Computations↗

Comparative Results from a CFD Challenge Over a 2D Three-Element High-Lift Airfoil

A high-lift workshop was held in May of 1993 at NASA Langley Research Center. A major part of the workshop centered on a blind test of various computational fluid dynamics (CFD) methods in which the flow about a two- dimensional (2D) three-element airfoil was computed without prior knowledge of the experimental data. The results of this 'blind' test revealed: (1) The Reynolds Averaged Navier-Stokes (RANS) methods generally showed less variability among codes than did potential/Euler solvers coupled with boundary-layer solution techniques. However, some of the coupled methods still provided excellent predictions. (2) Drag prediction using coupled methods agreed more closely with experiment than the RANS methods. Lift was more accurately predicted than drag for both methods. (3) The CFD methods did well in predicting lift and drag changes due to changes in Reynolds number, however, they did not perform as well when predicting lift and drag increments due to changing flap gap, (4) Pressures and skin friction compared favorably with experiment for most of the codes. (5) There was a large variability in most of the velocity profile predictions. Computational results predict a stronger siat wake than measured suggesting a missing component in turbulence modeling, perhaps curvature effects.

Klausmeyer, Steven M.↗

Assessment of Preconditioner for a USM3D Hierarchical Adaptive Nonlinear Method (HANIM) (Invited)

Enhancements to the previously reported mixed-element USM3D Hierarchical Adaptive Nonlinear Iteration Method (HANIM) framework have been made to further improve robustness, efficiency, and accuracy of computational fluid dynamic simulations. The key enhancements include a multi-color line-implicit preconditioner, a discretely consistent symmetry boundary condition, and a line-mapping method for the turbulence source term discretization. The USM3D iterative convergence for the turbulent flows is assessed on four configurations. The configurations include a two-dimensional (2D) bump-in-channel, the 2D NACA 0012 airfoil, a three-dimensional (3D) bump-in-channel, and a 3D hemisphere cylinder. The Reynolds Averaged Navier Stokes (RANS) solutions have been obtained using a Spalart-Allmaras turbulence model and families of uniformly refined nested grids. Two types of HANIM solutions using line- and point-implicit preconditioners have been computed. Additional solutions using the point-implicit preconditioner alone (PA) method that broadly represents the baseline solver technology have also been computed. The line-implicit HANIM shows superior iterative convergence in most cases with progressively increasing benefits on finer grids.

Pandya, Mohagna J.↗