Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous finite element”

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 487 records · Page 27

Computational mechanics needs study

In order to assess the needs in computational mechanics over the next decade, we formulated a questionnaire and contacted computational mechanics researchers and users in industry, government, and academia. As expected, we found a wide variety of computational mechanics usage and research. This report outlines the activity discussed with those contacts, as well as that in our own organizations. It should be noted that most of the contacts were made before the recent decline of the defense industry. Therefore, areas which are strongly defense-oriented may decrease in relative importance. In order to facilitate updating of this study, names of a few key researchers in each area are included as starting points for future literature surveys. These lists of names are not intended to represent those persons doing the best research in that area, nor are they intended to be comprehensive. They are, as previously stated, offered as starting points for future literature searches. Overall, there is currently a broad activity in computational mechanics in this country, with the breadth and depth increasing as more sophisticated software and faster computers become more available. The needs and desires of the workers in this field are as diverse as their background and organizational products. There seems to be some degree of software development in any organization (although the level of activity is highly variable from one organization to another) which has any research component in its mission. It seems, however, that there is considerable use of commercial software in almost all organizations. In most industrial research organizations, it appears that very little actual software development is contracted out, but that most is done in-house, using a mixture of funding sources. Government agencies vary widely in the ratio of in-house to out-house ratio. There is a considerable amount of experimental verification in most, but not all, organizations. Generally, the amount of experimental verification is more than we expected. Of all the survey contacts, one or two believe that the resources they are allocated are sufficient, but most do not. Some believe they have only half the resources they need. Some see their resource deficits as short-term, while others see it as a trend which will continue or perhaps worsen. The pessimism is stronger in the defense and aerospace industry. When considering only the nonlinear development efforts, there appears to be an even mix of geometric and material nonlinearity. There is not much particular emphasis in linear analysis unless it is in extension of current analysis capabilities to larger problems. The primary exception is concern about modeling of composites, where proven methodologies have trailed element and computer hardware development. Most of the people we spoke to use finite element techniques, but there is some finite difference and boundary element work ongoing. There is also some interest in multiple methods. Coupling of finite elements and boundary elements appears to be of high interest, since the two analysis types are complementary.

Griffin, O. Hayden, Jr.↗

The Prediction of Ducted Fan Engine Noise Via a Boundary Integral Equation Method

A computationally efficient Boundary Integral Equation Method (BIEM) for the prediction of ducted fan engine noise is presented. The key features of the BIEM are its versatility and the ability to compute rapidly any portion of the sound field without the need to compute the entire field. Governing equations for the BIEM are based on the assumptions that all acoustic processes are linear, generate spinning modes, and occur in a uniform flow field. An exterior boundary value problem (BVP) is defined that describes the scattering of incident sound by an engine duct with arbitrary profile. Boundary conditions on the duct walls are derived that allow for passive noise control treatment. The BVP is recast as a system of hypersingular boundary integral equations for the unknown duct surface quantities. BIEM solution methodology is demonstrated for the scattering of incident sound by a thin cylindrical duct with hard walls. Numerical studies are conducted for various engine parameters and continuous portions of the total pressure field are computed. Radiation and duct propagation results obtained are in agreement with the classical results of spinning mode theory for infinite ducts.

Boundary integral equation↗

Update on the NASA GRC Stirling Technology Development Project

The Department of Energy, NASA Glenn Research Center (GRC), and Stirling Technology Company (STC) are developing a free-piston Stirling convertor for a Stirling radioisotope power system (SRPS) to provide spacecraft on-board electric power for NASA deep space missions. The SRPS has recently been identified for potential use on the Europa Orbiter and Solar Probe Space Science missions. Stirling is also now being considered for unmanned Mars rovers. NASA GRC is conducting an in-house project to assist in developing the Stirling convertor for readiness for space qualification and mission implementation. As part of this continuing effort, the Stirling convertor will be further characterized under launch environment random vibration testing, methods to reduce convertor electromagnetic interference (EMI) will be developed, and an independent performance verification will be completed. Convertor life assessment and permanent magnet aging characterization tasks are also underway. Substitute organic materials for the linear alternator and piston bearing coatings for use in a high radiation environment have been identified and have now been incorporated in Stirling convertors built by STC for GRC. Electromagnetic and thermal finite element analyses for the alternator are also being conducted. This paper discusses the recent results and status for this NASA GRC in-house project.

Thieme, Lanny G.↗

Three-Dimensional Vibration Isolator for Suppressing High-Frequency Responses for Sage III Contamination Monitoring Package (CMP)

The Stratospheric Aerosol and Gas Experiment (SAGE) III - International Space Station (ISS) instrument will be used to study ozone, providing global, long-term measurements of key components of the Earth's atmosphere for the continued health of Earth and its inhabitants. SAGE III is launched into orbit in an inverted configuration on SpaceX;s Falcon 9 launch vehicle. As one of its four supporting elements, a Contamination Monitoring Package (CMP) mounted to the top panel of the Interface Adapter Module (IAM) box experiences high-frequency response due to structural coupling between the two structures during the SpaceX launch. These vibrations, which were initially observed in the IAM Engineering Development Unit (EDU) test and later verified through finite element analysis (FEA) for the SpaceX launch loads, may damage the internal electronic cards and the Thermoelectric Quartz Crystal Microbalance (TQCM) sensors mounted on the CMP. Three-dimensional (3D) vibration isolators were required to be inserted between the CMP and IAM interface in order to attenuate the high frequency vibrations without resulting in any major changes to the existing system. Wire rope isolators were proposed as the isolation system between the CMP and IAM due to the low impact to design. Most 3D isolation systems are designed for compression and roll, therefore little dynamic data was available for using wire rope isolators in an inverted or tension configuration. From the isolator FEA and test results, it is shown that by using the 3D wire rope isolators, the CMP high-frequency responses have been suppressed by several orders of magnitude over a wide excitation frequency range. Consequently, the TQCM sensor responses are well below their qualification environments. It is indicated that these high-frequency responses due to the typical instrument structural coupling can be significantly suppressed by a vibration passive control using the 3D vibration isolator. Thermal and contamination issues were also examined during the isolator selection period for meeting the SAGE III-ISS instrument requirements.

Li, Y.↗

Mixed finite element models for free vibrations of thin-walled beams

Simple, mixed finite element models are developed for the free vibration analysis of curved thin-walled beams with arbitrary open cross section. The analytical formulation is based on a Vlasov's type thin-walled beam theory with the effects of flexural-torsional coupling, transverse shear deformation and rotary inertia included. The fundamental unknowns consist of seven internal forces and seven generalized displacements of the beam. The element characteristic arrays are obtained by using a perturbed Lagrangian-mixed variational principle. Only C(sup o) continuity is required for the generalized displacements. The internal forces and the Lagrange multiplier are allowed to be discontinuous at interelement boundaries. Numerical results are presented to demonstrate the high accuracy and effectiveness of the elements developed. The standard of comparison is taken to be the solutions obtained by using 2-D plate/shell models for the beams.

Noor, Ahmed K.↗

Mixed finite element models for free vibrations of thin-walled beams

Simple, mixed finite element models are developed for the free vibration analysis of curved thin-walled beams with arbitrary open cross section. The analytical formulation is based on a Vlasov's type thin-walled beam theory with the effects of flexural-torsional coupling, transverse shear deformation and rotary inertia included. The fundamental unknowns consist of seven internal forces and seven generalized displacements of the beam. The element characteristic arrays are obtained by using a perturbed Lagrangian-mixed variational principle. Only C(sup o) continuity is required for the generalized displacements. The internal forces and the Lagrange multiplier are allowed to be discontinuous at interelement boundaries. Numerical results are presented to demonstrate the high accuracy and effectiveness of the elements developed. The standard of comparison is taken to be the solutions obtained by using 2-D plate/shell models for the beams.

Noor, Ahmed K.↗

Meshless Local Petrov-Galerkin Method for Bending Problems

Recent literature shows extensive research work on meshless or element-free methods as alternatives to the versatile Finite Element Method. One such meshless method is the Meshless Local Petrov-Galerkin (MLPG) method. In this report, the method is developed for bending of beams - C1 problems. A generalized moving least squares (GMLS) interpolation is used to construct the trial functions, and spline and power weight functions are used as the test functions. The method is applied to problems for which exact solutions are available to evaluate its effectiveness. The accuracy of the method is demonstrated for problems with load discontinuities and continuous beam problems. A Petrov-Galerkin implementation of the method is shown to greatly reduce computational time and effort and is thus preferable over the previously developed Galerkin approach. The MLPG method for beam problems yields very accurate deflections and slopes and continuous moment and shear forces without the need for elaborate post-processing techniques.

Phillips, Dawn R.↗

Role of Intersections in Fracture Connectivity

Networks of intersecting fractures often provide the flow paths through subsurface reservoirs. Assessing network connectivity is challenging because fracture intersections compose a vanishingly small fraction of the network void volume. In this paper, motivated by 3D X-ray imaging of the simplest element of fracture network, that is, two orthogonal fractures, we perform a percolation and finite-size scaling analysis to study the connectivity provided by fracture intersections. The conditions when an intersection enhances connectivity across a sample depend on spatial correlations in the fracture aperture distributions, on the stress state, and on the direction of flow. Here we consider three flow directions: (a) across intersections, (b) parallel to intersections and (c) around corners. For (a), intersections provide minimal enhancement of connectivity because they contribute little additional void area. For (b), intersections increase the probability of a connected path near threshold by enabling 3D connected pathways that are not possible in parallel fractures. Flow around corners, (c), is fundamentally the result of the intersection connecting two fractures in series and spatial correlations are broken around corners, suppressing the connectivity relative to (a). When the connected fractures are stressed equally, a joint percolation threshold emerges that continues to have scale invariance. However, when the fractures are stressed unequally, the system has mixed percolation without clearly defined percolation thresholds. In all cases, percolation probabilities are found to be scale dependent which has important consequences for the connectivity of larger fracture networks composed of the fundamental element studied here.

02 PETROLEUM↗

The Fractional Step Method Applied to Simulations of Natural Convective Flows

This paper describes research done to apply the Fractional Step Method to finite-element simulations of natural convective flows in pure liquids, permeable media, and in a directionally solidified metal alloy casting. The Fractional Step Method has been applied commonly to high Reynold's number flow simulations, but is less common for low Reynold's number flows, such as natural convection in liquids and in permeable media. The Fractional Step Method offers increased speed and reduced memory requirements by allowing non-coupled solution of the pressure and the velocity components. The Fractional Step Method has particular benefits for predicting flows in a directionally solidified alloy, since other methods presently employed are not very efficient. Previously, the most suitable method for predicting flows in a directionally solidified binary alloy was the penalty method. The penalty method requires direct matrix solvers, due to the penalty term. The Fractional Step Method allows iterative solution of the finite element stiffness matrices, thereby allowing more efficient solution of the matrices. The Fractional Step Method also lends itself to parallel processing, since the velocity component stiffness matrices can be built and solved independently of each other. The finite-element simulations of a directionally solidified casting are used to predict macrosegregation in directionally solidified castings. In particular, the finite-element simulations predict the existence of 'channels' within the processing mushy zone and subsequently 'freckles' within the fully processed solid, which are known to result from macrosegregation, or what is often referred to as thermo-solutal convection. These freckles cause material property non-uniformities in directionally solidified castings; therefore many of these castings are scrapped. The phenomenon of natural convection in an alloy under-going directional solidification, or thermo-solutal convection, will be explained. The development of the momentum and continuity equations for natural convection in a fluid, a permeable medium, and in a binary alloy undergoing directional solidification will be presented. Finally, results for natural convection in a pure liquid, natural convection in a medium with a constant permeability, and for directional solidification will be presented.

Westra, Douglas G.↗

Composite Nozzle/Thrust Chambers Analyzed for Low-Cost Boosters

The Low Cost Booster Technology Program is an initiative to minimize the cost of future liquid engines by using advanced materials and innovative designs, and by reducing engine complexity. NASA Marshall Space Flight Center s 60K FASTRAC Engine is one example where these design philosophies have been put into practice. This engine burns a liquid kerosene/oxygen mixture. It uses a one-piece, polymer composite thrust chamber/nozzle that is constructed of a tape-wrapped silica phenolic liner, a metallic injector interface ring, and a filament-wound epoxy overwrap. A cooperative effort between NASA Lewis Research Center s Structures Division and Marshall is underway to perform a finite element analysis of the FASTRAC chamber/nozzle under all the loading and environmental conditions that it will experience during its lifetime. The chamber/nozzle is a complex composite structure. Of its three different materials, the two composite components have distinctly different fiber architectures and, consequently, require separate material model descriptions. Since the liner is tape wrapped, it is orthotropic in the nozzle global coordinates; and since the overwrap is filament wound, it is treated as a monoclinic material. Furthermore, the wind angle on the overwrap varies continuously along the length of the chamber/nozzle.

Sullivan, Roy M.↗

Elastic-plastic stress concentrations around crack-like notches in continuous fiber reinforced metal matrix composites

Continuous fiber silicon-carbide/aluminum composite laminates with slits were tested statically to failure. Five different layups were examined: (0) sub 8, (0 sub 2/ + or - 45) sub s, (0/90) sub 2s), (0/ + or - 45/90 sub s), and (+ or - 45) sub 2s. Either a 9.5 or a 19 mm slit was machined in the center of each specimen. The strain distribution ahead of the slit tip was found experimentally with a series of strain gages bonded ahead of the slit tip. A three-dimensional finite element program (PAFAC) was used to predict the strain distribution ahead of the slit tip for several layups. For all layups, except the (0) sub 8, the yielding of the metal matrix caused the fiber stress concentration factor to increase with increasing load. This is contrary to the behavior seen in homogeneous materials where yielding causes the stress concentration to drop. For the (0) sub 8 laminate, yielding of the matrix caused a decrease in the fiber stress concentration. The finite element analysis predicted these trends correctly.

Johnson, W. S.↗

Elastic-plastic stress concentrations around crack-like notches in continuous fiber reinforced metal matrix composites

Continuous fiber silicon-carbide/aluminum composite laminates with slits were tested statically to failure. Five different layups were examined: (0) sub 8, (0 sub 2/ + or - 45) sub s, (0/90) sub 2s), (0/ + or - 45/90 sub s), and (+ or - 45) sub 2s. Either a 9.5 or a 19 mm slit was machined in the center of each specimen. The strain distribution ahead of the slit tip was found experimentally with a series of strain gages bonded ahead of the slit tip. A three-dimensional finite element program (PAFAC) was used to predict the strain distribution ahead of the slit tip for several layups. For all layups, except the (0) sub 8, the yielding of the metal matrix caused the fiber stress concentration factor to increase with increasing load. This is contrary to the behavior seen in homogeneous materials where yielding causes the stress concentration to drop. For the (0) sub 8 laminate, yielding of the matrix caused a decrease in the fiber stress concentration. The finite element analysis predicted these trends correctly.

Johnson, W. S.↗

Integrated Thermal Response Tool for Earth Entry Vehicles

A system is presented for multi-dimensional, fully-coupled thermal response modeling of hypersonic entry vehicles. The system consists of a two-dimensional implicit thermal response, pyrolysis and ablation program (TITAN), a commercial finite-element thermal and mechanical analysis code (MARC), and a high fidelity Navier-Stokes equation solver (GIANTS). The simulations performed by this integrated system include hypersonic flow-field, fluid and solid interaction, ablation, shape change, pyrolysis gas generation and flow, and thermal response of heatshield and structure. The thermal response of the ablating and charring heatshield material is simulated using TITAN, and that of the underlying structural is simulated using MARC. The ablating heatshield is treated as an outer boundary condition of the structure, and continuity conditions of temperature and heat flux are imposed at the interface between TITAN and MARC. Aerothermal environments with fluid and solid interaction are predicted by coupling TITAN and GIANTS through surface energy balance equations. With this integrated system, the aerothermal environments for an entry vehicle and the thermal response of both the heatshield and the structure can be obtained simultaneously. Representative computations for a proposed blunt body earth entry vehicle are presented and discussed in detail.

Chen, Y.-K.↗

Debond Analyses for Stitched Composite Structures

The effect of stitching on mode I and mode II strain energy release rates for debond configurations is studied using an analysis based on plate finite elements and the virtual crack closure technique. The stitches were modeled as discrete nonlinear fastener elements with a compliance determined by experiment. The axial and shear behavior of the stitches was considered with both the compliances and failure loads assumed to be independent. The mode I strain energy release rate, G(sub I), was shown to decrease once the debond had grown beyond the first row of stitches and was reduced to zero for long debonds, however, the mode II strain energy release rate, G(sub II), continued to be of significant magnitude over the range of debond lengths considered.

Glaessgen, E. H.↗

Two-Phase Turbulence Statistics from High Fidelity Dispersed Droplet Flow Simulations in a Pressurized Water Reactor (PWR) Sub-Channel with Mixing Vanes

In the dispersed flow film boiling regime (DFFB), which exists under post-LOCA (loss-of-coolant accident) conditions in pressurized water reactors (PWRs), there is a complex interplay between droplet dynamics and turbulence in the surrounding steam. Experiments have accredited particular significance to droplet collision with the spacer-grids and mixing vane structures and their consequent positive feedback to the heat transfer recorded in the immediate downstream vicinity. Enabled by high-performance computing (HPC) systems and a massively parallel finite element-based flow solver—PHASTA (Parallel Hierarchic Adaptive Stabilized Transient Analysis)—this work presents high fidelity interface capturing, two-phase, adiabatic simulations in a PWR sub-channel with spacer grids and mixing vanes. Selected flow conditions for the simulations are informed by the experimental data found in the literature, including the steam Reynolds number and collision Weber number (Wec={40,80}), and are characteristic of the DFFB regime. Data were collected from the simulations at an unprecedented resolution, which provides detailed insights into the continuous phase turbulence statistics, highlighting the effects of the presence of droplets and the comparative effect of different Weber numbers on turbulence in the surrounding steam. Further, axial evolution of droplet dynamics was analyzed through cross-sectionally averaged quantities, including droplet volume, surface area and Sauter mean diameter (SMD). The downstream SMD values agree well with the existing empirical correlations for the selected range of Wec. The high-resolution data repository from the simulations herein is expected to be of significance to guide model development for system-level thermal hydraulic codes.

Saini, Nadish↗

Formulation of full state feedback for infinite order structural systems

The estimation of exact displacement and displacement rate feedback kernels from finite dimensional control solutions based on finite element structural models is discussed. These kernels are then transformed to equivalent curvature and curvature rate feedback kernels. These curvature kernels are augmented with single point displacement and rotation feedback to account for rigid body motions. A growing class of sensors known as area-averaging sensors is used to measure the curvature and curvature rate state functions. The output of area-averaging sensors equals the convolution of all structural curvature states with the spatial sensitivity function of the sensors. Transforming the discrete feedback gains into continuous feedback kernels and employing area-averaging sensors make it possible to implement full state feedback for infinite order structural systems.

Miller, David W.↗

Structural Assessment of Advanced Composite Tow-Steered Shells

The structural performance of two advanced composite tow-steered shells, manufactured using a fiber placement system, is assessed using both experimental and analytical methods. The fiber orientation angles vary continuously around the shell circumference from 10 degrees on the shell crown and keel, to 45 degrees on the shell sides. The two shells differ in that one shell has the full 24-tow course applied during each pass of the fiber placement system, while the second shell uses the fiber placement system s tow drop/add capability to achieve a more uniform shell wall thickness. The shells are tested in axial compression, and estimates of their prebuckling axial stiffnesses and bifurcation buckling loads are predicted using linear finite element analyses. These preliminary predictions compare well with the test results, with an average agreement of approximately 10 percent.

Wu, K. Chauncey↗

Plasticity-mediated deformation instabilities in thin film-compliant substrate systems: direct three-dimensional simulations

Abstract Surface wrinkles driven by mechanical instability commonly form in thin-film structures attached to a compliant substrate. In this study, a recently developed computational approach is employed to simulate the formation and transformation of wrinkles involving plastic yielding of the thin film. The three-dimensional (3D) finite element models contain an embedded imperfection at the film-substrate interface, serving to trigger the bifurcation modes. Successful application of this technique to allow for film plasticity is demonstrated, including the evolution of 3D surface patterns and their correlation with the overall load–displacement response. The simulations reveal that plastic yielding transforms the surface instability patterns into more localized forms. Under uniaxial loading, the sinusoidal elastic wrinkles undergo the wrinkle-to-fold transition. With equi-biaxial loading, the initial square-checkerboard array turns into continuous tall ridges along the 45° directions. In both loading modes, the plasticity-induced instability patterns are only partially relieved upon unloading, leaving permanent features on the surface.

36 MATERIALS SCIENCE↗