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 1,027 records · Page 57

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

Identification of High-Frequency Transverse Acoustic Modes in Multi-Nozzle Can Combustors

Abstract High frequency thermoacoustic instabilities are problematic for lean-premixed gas turbines. Identifying which acoustic mode is being excited is important, in that it provides insight into potential mitigation measures and mechanical stress/life. However, the frequency spacing between modes becomes quite close for high frequency instabilities in a can combustor. This makes it difficult to distinguish between the modes (e.g., the first transverse mode vs. a higher order axial/mixed mode) based upon frequency calculations alone, which inevitably have uncertainties in boundary conditions, temperature profiles, and combustion response. This paper presents a methodology to simultaneously identify the acoustic mode shapes in the axial and azimuthal directions from acoustic pressure measurements. Multiple high temperature pressure transducers, located at distinct axial and azimuthal positions, are flush mounted in the combustor wall. The measured pressure oscillations from each sensor are then used to reconstruct the pressure distributions by using a least squares method in conjunction with a solution of a three dimensional wave equation. In order to validate the methodology, finite element method (FEM) calculations with estimated post-flame temperature is used to provide the candidate frequencies and corresponding mode shapes. The results demonstrate the reconstructed mode shapes and standing/spinning character of transverse waves, as well as the associated frequencies, both of which are consistent with the FEM predictions. Nodal line location was also extracted from the experimental data during the instabilities in the pressure data.

Kim, J.↗

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

Magneto-Structural Combined Dimensional and Topology Optimization of Interior Permanent Magnet Synchronous Machine Rotors

A magneto-structural combined dimensional and topology optimization approach for interior permanent magnet synchronous machine (IPMSM) rotor design is proposed using the Solid Isotropic Material Penalization (SIMP) density-based topology optimization method. This method optimizes the location and dimensions of the permanent magnets (PMs) embedded in the rotor, while simultaneously optimizing the lamination flux barriers using topology optimization.To allow for the location and dimensions of the PMs to change while retaining the same number of mesh elements and connectivity, the mesh must be deformed. Mesh deformation is accomplished using a Laplace smoothing technique. To increase the range of PM displacement and size changes, the non-PM areas of the rotor design domain are divided into sub-regions along whose boundaries the mesh nodes are redistributed also using a Laplace smoothing algorithm. Three examples of IPMSM rotor optimizations are presented to demonstrate the proposed approach.

33 ADVANCED PROPULSION SYSTEMS↗

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 the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems

In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗

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

Probabilistic Structural Analysis Methods (PSAM) for Select Space Propulsion System Components

Probabilistic Structural Analysis Methods (PSAM) are described for the probabilistic structural analysis of engine components for current and future space propulsion systems. Components for these systems are subjected to stochastic thermomechanical launch loads. Uncertainties or randomness also occurs in material properties, structural geometry, and boundary conditions. Material property stochasticity, such as in modulus of elasticity or yield strength, exists in every structure and is a consequence of variations in material composition and manufacturing processes. Procedures are outlined for computing the probabilistic structural response or reliability of the structural components. The response variables include static or dynamic deflections, strains, and stresses at one or several locations, natural frequencies, fatigue or creep life, etc. Sample cases illustrates how the PSAM methods and codes simulate input uncertainties and compute probabilistic response or reliability using a finite element model with probabilistic methods.

Source record↗