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 829 records · Page 46

Discretized energy minimization in a wave guide with point sources

An anti-noise problem on a finite time interval is solved by minimization of a quadratic functional on the Hilbert space of square integrable controls. To this end, the one-dimensional wave equation with point sources and pointwise reflecting boundary conditions is decomposed into a system for the two propagating components of waves. Wellposedness of this system is proved for a class of data that includes piecewise linear initial conditions and piecewise constant forcing functions. It is shown that for such data the optimal piecewise constant control is the solution of a sparse linear system. Methods for its computational treatment are presented as well as examples of their applicability. The convergence of discrete approximations to the general optimization problem is demonstrated by finite element methods.

Propst, G.↗

Implementation of Combined Feather and Surface-Normal Ice Growth Models in LEWICE/X

Experimental observations have shown that discrete rime ice growths called feathers, which grow in approximately the direction of water droplet impingement, play an important role in the growth of ice on accreting surfaces for some thermodynamic conditions. An improved physical model of ice accretion has been implemented in the LEWICE 2D panel-based ice accretion code maintained by the NASA Lewis Research Center. The LEWICE/X model of ice accretion explicitly simulates regions of feather growth within the framework of the LEWICE model. Water droplets impinging on an accreting surface are withheld from the normal LEWICE mass/energy balance and handled in a separate routine; ice growth resulting from these droplets is performed with enhanced convective heat transfer approximately along droplet impingement directions. An independent underlying ice shape is grown along surface normals using the unmodified LEWICE method. The resulting dual-surface ice shape models roughness-induced feather growth observed in icing wind tunnel tests. Experiments indicate that the exact direction of feather growth is dependent on external conditions. Data is presented to support a linear variation of growth direction with temperature and cloud water content. Test runs of LEWICE/X indicate that the sizes of surface regions containing feathers are influenced by initial roughness element height. This suggests that a previous argument that feather region size is determined by boundary layer transition may be incorrect. Simulation results for two typical test cases give improved shape agreement over unmodified LEWICE.

Matthew T Velazquez↗

Prediction of Sound Waves Propagating Through a Nozzle Without/With a Shock Wave Using the Space-Time CE/SE Method

The benchmark problems in Category 1 (Internal Propagation) of the third Computational Aeroacoustics (CAA) Work-shop sponsored by NASA Glenn Research Center are solved using the space-time conservation element and solution element (CE/SE) method. The first problem addresses the propagation of sound waves through a nearly choked transonic nozzle. The second one concerns shock-sound interaction in a supersonic nozzle. A quasi one-dimension CE/SE Euler solver for a nonuniform mesh is developed and employed to solve both problems. Numerical solutions are compared with the analytical solution for both problems. It is demonstrated that the CE/SE method is capable of solving aeroacoustic problems with/without shock waves in a simple way. Furthermore, the simple nonreflecting boundary condition used in the CE/SE method which is not based on the characteristic theory works very well.

Wang, Xiao-Yen↗

Advancing xEMU Lunar Dust Mitigation Devices

Fine, electrically charged, glass like dust particles caused significant damage to the Apollo EMU9 during lunar EVAs, identified as one of the greatest challenges to future exploration. Passive Lunar Dust Mitigation Devices (LDMD) were developed, within a SBIR Phase II, to prohibit this dust from interrupting venting space suit component operation. A Computational Fluid-Dynamics and Discrete Element Method (CFD-DEM) Simulation Tool was developed at the University of Colorado, Boulder to predict venting gas flow ability to self-clean adhered dust particles from LDMD surfaces. Lunar dust properties (i.e., adhesion and cohesion strengths) required to complete Simulation Tool analysis are relatively unknown due to considerable differences between the Earth and the Moon (i.e., gravity, humidity) and due to an absence of dust particles in their native state. Analytically determining gas flow velocity, density and direction within the fluid Boundary Layer, microns from LDMD surfaces presented a second challenge. Dusty Plasma Chamber testing is being performed at Auburn University, Auburn to observe electrostatically charged dust behavior as it adheres to LDMD prototypes and specific geometry and is then blown away by metered gas flow. Test articles were developed to offer insight into the impact of different flow geometries, surface roughness and dust removal within the gas Boundary-Layer. Observed dust behavior is currently being developed to support the CFD-DEM analysis. Many Simulation Tool analytical cases have been processed to support the intention of completing sensitivity studies to assess how different dust adherence values and Boundary Layer fluid properties impact LDMD self-cleaning effectivity.

Thomas J Stapleton↗

Worst case estimation of homology design by convex analysis

The methodology of homology design is investigated for optimum design of advanced structures. for which the achievement of delicate tasks by the aid of active control system is demanded. The proposed formulation of homology design, based on the finite element sensitivity analysis, necessarily requires the specification of external loadings. The formulation to evaluate the worst case for homology design caused by uncertain fluctuation of loadings is presented by means of the convex model of uncertainty, in which uncertainty variables are assigned to discretized nodal forces and are confined within a conceivable convex hull given as a hyperellipse. The worst case of the distortion from objective homologous deformation is estimated by the Lagrange multiplier method searching the point to maximize the error index on the boundary of the convex hull. The validity of the proposed method is demonstrated in a numerical example using the eleven-bar truss structure.

Yoshikawa, N.↗

Safer Systems: A NextGen Aviation Safety Strategic Goal

The Joint Planning and Development Office (JPDO), is charged by Congress with developing the concepts and plans for the Next Generation Air Transportation System (NextGen). The National Aviation Safety Strategic Plan (NASSP), developed by the Safety Working Group of the JPDO, focuses on establishing the goals, objectives, and strategies needed to realize the safety objectives of the NextGen Integrated Plan. The three goal areas of the NASSP are Safer Practices, Safer Systems, and Safer Worldwide. Safer Practices emphasizes an integrated, systematic approach to safety risk management through implementation of formalized Safety Management Systems (SMS) that incorporate safety data analysis processes, and the enhancement of methods for ensuring safety is an inherent characteristic of NextGen. Safer Systems emphasizes implementation of safety-enhancing technologies, which will improve safety for human-centered interfaces and enhance the safety of airborne and ground-based systems. Safer Worldwide encourages coordinating the adoption of the safer practices and safer systems technologies, policies and procedures worldwide, such that the maximum level of safety is achieved across air transportation system boundaries. This paper introduces the NASSP and its development, and focuses on the Safer Systems elements of the NASSP, which incorporates three objectives for NextGen systems: 1) provide risk reducing system interfaces, 2) provide safety enhancements for airborne systems, and 3) provide safety enhancements for ground-based systems. The goal of this paper is to expose avionics and air traffic management system developers to NASSP objectives and Safer Systems strategies.

Darr, Stephen T.↗

Mosaic wedge-and-strip arrays for large format microchannel plate detectors

The authors present a novel method for joining wedge-and-strip patterns on single anodes in a mosaic array. With only a modest increase in complexity over three-conductor anodes currently in use, the ultimate detector position resolution can be significantly improved, and large-format microchannel plate detectors with pore-size-limited resolution are made possible. The problem of the transition from one anode to the next has been solved with a novel linear encoding scheme, which exhibits essentially distortionless behavior at boundaries parallel to the conducting elements and only slight distortion at the orthogonal boundaries. The ultimate resolution for two anode designs, one designed for large-format imaging and the other for high-resolution spectroscopy, is also predicted.

Martin, Christopher↗

Wallops Ship Surveillance System

Approved as a Wallops control center backup system, the Wallops Ship Surveillance Software is a day-of-launch risk analysis tool for spaceport activities. The system calculates impact probabilities and displays ship locations relative to boundary lines. It enables rapid analysis of possible flight paths to preclude the need to cancel launches and allow execution of launches in a timely manner. Its design is based on low-cost, large-customer- base elements including personal computers, the Windows operating system, C/C++ object-oriented software, and network interfaces. In conformance with the NASA software safety standard, the system is designed to ensure that it does not falsely report a safe-for-launch condition. To improve the current ship surveillance method, the system is designed to prevent delay of launch under a safe-for-launch condition. A single workstation is designated the controller of the official ship information and the official risk analysis. Copies of this information are shared with other networked workstations. The program design is divided into five subsystems areas: 1. Communication Link -- threads that control the networking of workstations; 2. Contact List -- a thread that controls a list of protected item (ocean vessel) information; 3. Hazard List -- threads that control a list of hazardous item (debris) information and associated risk calculation information; 4. Display -- threads that control operator inputs and screen display outputs; and 5. Archive -- a thread that controls archive file read and write access. Currently, most of the hazard list thread and parts of other threads are being reused as part of a new ship surveillance system, under the SureTrak project.

Smith, Donna C.↗

A rigorous solution of the Navier-Stokes equations for unsteady viscous flow at high Reynolds numbers around oscillating airfoils

A method based on the Navier-Stokes equations was developed for analyzing the unsteady incompressible viscous flow around oscillating airfoils at high Reynolds numbers. The Navier-Stokes equations have been integrated in their classical Helmholtz vorticity transport equation form, and the instantaneous velocity field at each time step was determined by the solution of Poisson's equation. A refined finite element was utilized to allow for a conformable solution of the stream function and its first space derivatives at the element interfaces. A corresponding set of accurate boundary conditions was applied; thus obtaining a rigorous solution for the velocity field. The details of the computational procedure and examples of computed results describing the unsteady flow characteristics around the airfoil are presented.

Bratanow, T.↗

Test and Analysis of Buckling-Critical Stiffened Metallic Launch Vehicle Cylinders

A summary of a NASA design analysis and test program on the buckling of large-scale, integrally-stiffened metallic cylinders is presented. The test article designs were developed based on proposed NASA launch vehicle cylinder designs and span a significant portion of the design space. Various loading conditions were applied to the cylinders and include different combinations of axial compression, bending, and internal pressure loads that simulate typical launch vehicle loading scenarios. The data gathered from this test program is being used to develop and validate new analysis-based knockdown factors and design guidelines for these stiffened metallic cylinders. In this paper, the test article designs and fabrication methods are described along with the test facilities and instrumentation. Selected test and finite element analysis results are presented and compared and are used to illustrate the typical response characteristics of the stiffened metallic cylinders considered. Overall, good qualitative agreement is found, however, several discrepancies in the results were identified across several of the tests. The discrepancies were investigated thoroughly and can be attributed to variations in the as-built skin and stiffener geometry, variations in the measured geometric imperfection, and modeling assumptions associated with the boundary conditions, and loading imperfections. Based on these findings, several refinements were made to the finite element models which significantly improved the correlation between the test and analysis results. These modeling refinements are described and the updated analysis results are presented.

Hilburger, Mark W.↗

F-ANG+: A 3-D Augmented-Stencil Face-Averaged Nodal-Gradient Cell-Centered Finite-Volume Method for Hypersonic Flows

We describe the extension of a 2-D simplified face-averaged nodal-gradient (F-ANG) method to 3-D and demonstrate that the 3-D simplified F-ANG method is accomplished by augmenting the nodecentered gradient least squares stencil. This augmented stencil F-ANG method is shown to result in advection and diffusion schemes that are stable for hexahedral, prismatic, pyramidal and tetrahedral cells without having to resort to cell-averaged nodal gradients. In addition, we describe the modifications to the augmented stencil required to support the use of wall function boundary conditions. Finally we describe a consistent, face-stencil based multi-dimensional limiter procedure (MLP), and show it to be fully consistent and compatible with the linearity-preserving unstructured- MUSCL (LP-U-MUSCL) scheme for all values of kappa. These methods and schema are implemented in the cell-centered finite-volume code VULCAN-CFD, which is then used to investigate whether the robustness improvements demonstrated in 2-D carry over to 3-D by computing hypersonic flows using mixed-element grids as well as highly adapted tetrahedral grids.

Weighted Least-Squares↗

The FEM-R-Matrix Approach: Use of Mixed Finite Element and Gaussian Basis Sets for Electron Molecule Collisions

For the calculation of electron molecule collision cross sections R-matrix methods automatically take advantage of the division of configuration space into an inner region (I) bounded by radius tau b, where the scattered electron is within the molecular charge cloud and the system is described by an correlated Configuration Interaction (CI) treatment in close analogy to bound state calculations, and an outer region (II) where the scattered electron moves in the long-range multipole potential of the target and efficient analytic methods can be used for solving the asymptotic Schroedinger equation plus boundary conditions.

Thuemmel, Helmar T.↗

Unstructured viscous grid generation by advancing-front method

A new method of generating unstructured triangular/tetrahedral grids with high-aspect-ratio cells is proposed. The method is based on new grid-marching strategy referred to as 'advancing-layers' for construction of highly stretched cells in the boundary layer and the conventional advancing-front technique for generation of regular, equilateral cells in the inviscid-flow region. Unlike the existing semi-structured viscous grid generation techniques, the new procedure relies on a totally unstructured advancing-front grid strategy resulting in a substantially enhanced grid flexibility and efficiency. The method is conceptually simple but powerful, capable of producing high quality viscous grids for complex configurations with ease. A number of two-dimensional, triangular grids are presented to demonstrate the methodology. The basic elements of the method, however, have been primarily designed with three-dimensional problems in mind, making it extendible for tetrahedral, viscous grid generation.

Pirzadeh, Shahyar↗

Unstructured viscous grid generation by advancing-layers method

A new method of generating unstructured triangular/tetrahedral grids with high-aspect-ratio cells is proposed. The method is based on new grid-marching strategy referred to as 'advancing-layers' for construction of highly stretched cells in the boundary layer and the conventional advancing-front technique for generation of regular, equilateral cells in the inviscid-flow region. Unlike the existing semi-structured viscous grid generation techniques, the new procedure relies on a totally unstructured advancing-front grid strategy resulting in a substantially enhanced grid flexibility and efficiency. The method is conceptually simple but powerful, capable of producing high quality viscous grids for complex configurations with ease. A number of two-dimensional, triangular grids are presented to demonstrate the methodology. The basic elements of the method, however, have been primarily designed with three-dimensional problems in mind, making it extendible for tetrahedral, viscous grid generation.

Pirzadeh, Shahyar↗

Thermal Buckling Analysis of Rectangular Panels Subjected to Humped Temperature Profile Heating

This research investigates thermal buckling characteristics of rectangular panels subjected to different types of humped temperature profile heating. Minimum potential energy and finite-element methods are used to calculate the panel buckling temperatures. The two methods give fairly close thermal buckling solutions. 'Buckling temperature magnification factor of the first kind, eta' is established for the fixed panel edges to scale up the buckling solution of uniform temperature loading case to give the buckling solution of the humped temperature profile loading cases. Also, 'buckling temperature magnification factor of the second kind, xi' is established for the free panel edges to scale up the buckling solution of humped temperature profile loading cases with unheated boundary heat sinks to give the buckling solutions when the boundary heat sinks are heated up.

Ko, William I.↗

3-D-Observation of Matrix of MIL 090657 Meteorite by Absorption-Phase Tomography

MIL 090657 meteorite (CR2.7) is one of the least altered primitive carbonaceous c hondrites [1]. This meteorite has amorphous silicates like GEMS (glass with embedded metal and sulfide), which are characteristically contained in cometary dust, in matrix [2,3] as with the Paris meteorite [4]. Three lithologies have been recognized; lithology-1 (L 1) dominated by submicron anhydrous silicates, lithology-2 (L2) by GEMS-like amorphous silicates and lithology-3 (L3) by phyllosilicates [2]. Organic materials are abundant in L 1 and L2 [2,3]. L 1 and L2 were further divided into sub-lithology respectively based on their textures and compositions [5]. These studies were performed by 2D SEM and TEM observations of sample surfaces and thin sections that are unable to reveal what constitute each lithology and how these lithologies are distributed and related to each other. This information will provide important insights into alteration and aggregation processes on asteroids and in the early solar nebula. In this study, MIL 090657 matrix was examined in 3D using two types of X-ray tomography; DET (dual-energy tomography) [6] and SIXM (scanning-imaging X-ray microscopy) [7]. Mineral phases can be discriminated based on absorption contrasts at two different X-ray energies in DET. In SIXM, materials composed of light elements such as water or organic materials can be identified based on phase and absorption contrasts. By combining these methods, we can discriminate not only organic materials from voids but also hydrous alteration products, such as hydrated silicates and carbonates, from anhydrous minerals [8]. In this study, we first observed cross sections of MIL 090657 matrix fragments C1 00 mm) in detail using FE-SEM/ EDS. Based on the results, three house-shaped samples (3 0 -50 mm) were extracted from L 1, L2 and their boundary (H1, H3 and H5, respectively) using FIB. 3D imaging of these samples were conducted at BL47XU of SPring-8, a synchrotron radiation facility, with -30-40 nm/ voxel and - 70-80 nm/ voxel at 7keV and 8keV in DET and -100 nm/ voxel at 8keV in SIXM. We found new lithologies that we named L4, L5 and L6 in H1 and H3 in addition to L 1 and L2. L4. L5 and L6 are mainly composed of probably phyllosilicates with different Fe contents. Sulfide and framboidal magnetite were recognized in L4. L5 includes magnetite and carbonate and L6 includes anhydrous silicates having cracks inside. L 1, L2, L4 and L5 are porous while few voids were observed in L6. L4 adjoins to L 1 with boundary, which is not very distinct. L2, L5 and L6 adjoin to each other, and the boundaries of L6 with L2 and L5 are clear. In H5, coarse mineral grains C5-1 0 mm) such as Fe-metal and enstatite are present in L 1 and L2. L 1-L2 boundary is not sharp in 3 D. In conclusion, we found a variety of lithologies by 3D observation for the first time, suggesting that t he MIL 0 90657 meteorite experienced complex alteration and aggregation histories. As L2 is dominated by amorphous silicates, which are extremely susceptible to aqueous alteration, t his is p resumed to be t he most primitive lithology. The contact between L2 and phyllosilicate-bearing lithologies (L5 and L6) with clear boundaries indicates that they were aggregated after aqueous alteration of L5 and L6. The indistinct boundary between L 1 and L2 is suggesting that these two lithologies might originally be the same aggregate composed of amorphous silicates and coarse mineral grains. L 1 might have experienced weak aqueous alteration followed by mild thermal alteration [2], while L2 did not undergo aqueous alteration.

Miyama, Sugimoto↗

Automatic Incorporation of Flaws into Cad and FEA Models Based On Measured NDE Data

Automatic incorporation of non-destructive evaluation (NDE) data into CAD and finite element (FE) models for structural analysis is essential for accurate residual-life and strength analyses of composite structures. Currently, two significant but independent challenges exist for incorporating NDE data in such analyses: interpretation of the NDE signals to locate and characterize the structural defects and the incorporation of that information into structural analysis tools. The first challenge will be addressed through linear inversion of calculated synthetic responses generated by a simulation engine to evaluate flaw size, location, and other characteristics. The second challenge will be addressed by building on prior and parallel defect modeling efforts to integrate suitable flaw models into the CAD/finite-element representation of a specimen. Ultrasonic NDE information is usually represented using C-Scans of time-of-flight or intensity information. Other modalities such as thermography and computed tomography (CT) could also be used to characterize the structure. This information will first be projected into material coordinates and then aligned with the specimen’s finite element mesh. The NDE data can then be incorporated by inserting flaw geometry into the CAD and FE analysis model and automatic modification of local material properties and boundary conditions. This process will enable the use of structural analysis tools even for practical modalities such as ultrasound and thermography which may not directly give accurate-enough geometric information to fully define the flaw. The full flaw characteristics will be inferred using prior expert and modeling knowledge to categorize flaw signatures. Outputs for this project will include demonstrated methods for creating structural analysis geometric models that incorporate flaw characterization information; these models can then be used directly in popular commercial FEA software. In addition, the primary deliverable for the project will include a software that will be able to automatically incorporate different kinds of flaws into an existing FE models. The models will be generated in standard file formats to enable direct import into the analysis software.

Adarsh Krishnamurthy↗

Analysis of three-dimensional unsteady flow around oscillating wings

A method based on the Navier-Stokes equations was developed for determining analytically the three-dimensional unsteady flow patterns around oscillating wings. The Helmholz vorticity transport equations were discretized in three-dimensional finite element form from a variational formulation and integrated numerically. At each time step of the numerical integration the velocity field was calculated from the representation of the three-dimensional wing by a system of optimized distribution of vortices in space. During the numerical integration of the vorticity transport equations the time-dependent boundary conditions on the wing were specified as external constraint conditions. Examples of obtained results describing the three-dimensional unsteady flow around a wing were presented.

Bratanow, T.↗