Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element methods”

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

An immersed boundary method for fluid–structure interaction with compressible multiphase flows

This paper presents a two-dimensional immersed boundary method for fluid–structure interaction with compressible multiphase flows involving large structure deformations. This method involves three important parts: flow solver, structure solver and fluid–structure interaction coupling. In the flow solver, the compressible multiphase Navier–Stokes equations for ideal gases are solved by a finite difference method based on a staggered Cartesian mesh, where a fifth-order accuracy Weighted Essentially Non-Oscillation (WENO) scheme is used to handle spatial discretization of the convective term, a fourth-order central difference scheme is employed to discretize the viscous term, the third-order TVD Runge–Kutta scheme is used to discretize the temporal term, and the level-set method is adopted to capture the multi-material interface. In this work, the structure considered is a geometrically non-linear beam which is solved by using a finite element method based on the absolute nodal coordinate formulation (ANCF). The fluid dynamics and the structure motion are coupled in a partitioned iterative manner with a feedback penalty immersed boundary method where the flow dynamics is defined on a fixed Lagrangian grid and the structure dynamics is described on a global coordinate. We perform several validation cases (including fluid over a cylinder, structure dynamics, flow induced vibration of a flexible plate, deformation of a flexible panel induced by shock waves in a shock tube, an inclined flexible plate in a hypersonic flow, and shock-induced collapse of a cylindrical helium cavity in the air), and compare the results with experimental and other numerical data. The present results agree well with the published data and the current experiment. Finally, we further demonstrate the versatility of the present method by applying it to a flexible plate interacting with multiphase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Calculation of reinforced-concrete frame strength under a simultaneous static cross section load and a column lateral impact

When designing buildings with reinforced concrete that are planned to resist dynamic loads it is necessary to calculate this structural behavior under operational static and emergency impact and blast loads. Calculations of the structures under shock-wave loads can be performed by solving dynamic equations that do not consider static loads. Due to this fact the calculation of reinforced concrete frame under a simultaneous static and dynamic load in full 3d settings becomes a very non trivial and resource consuming problem. This problem can be split into two tasks. The first one is a shock-wave problem that can be solved using software package RANET-3, which allows solving the problem using finite elements method adapted for dynamic task. This method calculates strain-stress state of the material and its dynamic destruction, which is considered as growth and consolidation of micro defects under loading. On the second step the results of the first step are taken as input parameters for quasi static calculation of simultaneous static and dynamic load using finite elements method in AMP Civil Engineering-11.

36 MATERIALS SCIENCE↗

Prediction of high temperature metal matrix composite ply properties

The application of the finite element method (superelement technique) in conjunction with basic concepts from mechanics of materials theory is demonstrated to predict the thermomechanical behavior of high temperature metal matrix composites (HTMMC). The simulated behavior is used as a basis to establish characteristic properties of a unidirectional composite idealized an as equivalent homogeneous material. The ply properties predicted include: thermal properties (thermal conductivities and thermal expansion coefficients) and mechanical properties (moduli and Poisson's ratio). These properties are compared with those predicted by a simplified, analytical composite micromechanics model. The predictive capabilities of the finite element method and the simplified model are illustrated through the simulation of the thermomechanical behavior of a P100-graphite/copper unidirectional composite at room temperature and near matrix melting temperature. The advantage of the finite element analysis approach is its ability to more precisely represent the composite local geometry and hence capture the subtle effects that are dependent on this. The closed form micromechanics model does a good job at representing the average behavior of the constituents to predict composite behavior.

Caruso, J. J.↗

An interactive graphics system to facilitate finite element structural analysis

The characteristics of an interactive graphics systems to facilitate the finite element method of structural analysis are described. The finite element model analysis consists of three phases: (1) preprocessing (model generation), (2) problem solution, and (3) postprocessing (interpretation of results). The advantages of interactive graphics to finite element structural analysis are defined.

Burk, R. C.↗

Finite element analyses of plane thermal entry-length flows

It is pointed out that the design of actively cooled engine and airframe structures for hypersonic flight requires strong interaction between fluid, thermal, and structural analyses. Research programs are being conducted to develop an integrated thermal-structural analysis capability based upon the finite element method. In connection with these programs, a finite element engineering approach has been developed based upon a number of assumptions customarily used in practical heat transfer analysis. An evaluation is conducted of the finite element engineering approach for entry-length flows. The evaluation makes use of an analysis of plane thermal entry-length flows by a finite element approach utilizing mean fluid temperatures with a convection coefficient. Another analysis is based on a finite element continuum approach, taking into account the rigorous momentum, mass, and energy equations of viscous, incompressible flow. The engineering approach was found to give excellent agreement with the rigorous continuum approach for Peclet number values not less than 100.

Thornton, E. A.↗

A spectral element method for fluid dynamics - Laminar flow in a channel expansion

A spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible Navier-Stokes equations. In the spectral element discretization, the computational domain is broken into a series of elements, and the velocity in each element is represented as a high-order Lagrangian interpolant through Chebyshev collocation points. The hyperbolic piece of the governing equations is then treated with an explicit collocation scheme, while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle. The implementation of the technique is demonstrated on a one-dimensional inflow-outflow advection-diffusion equation, and the method is then applied to laminar two-dimensional (separated) flow in a channel expansion. Comparisons are made with experiment and previous numerical work.

Patera, A. T.↗

Hierarchical Strategy for Rapid Analysis Environment

A new philosophy is developed wherein the hierarchical definition of data is made use of in creating a better environment to conduct analyses of practical problems. This system can be adapted to conduct virtually any type of analysis, since this philosophy is not bound to any specific kind of analysis. It provides a framework to manage different models and its results and more importantly, the interaction between the different models. Thus, it is ideal for many types of finite element analyses like globalAoca1 analysis and those that involve multiple scales and fields. The system developed during the course of this work is just a demonstrator of the basic concepts. A complete implementation of this strategy could potentially make a major impact on the way analyses are conducted. It could considerably reduce the time frame required to conduct the analysis of real-life problems by efficient management of the data involved and reducing the human effort involved. It also helps in better decision making because of more ways to interpret the results. The strategy has been currently implemented for structural analysis, but with more work it could be extended to other fields of science when the finite element method is used to solve the differential equations numerically. This report details the work that has been done during the course of this project and its achievements and results. The following section discusses the meaning of the word hierarchical and the different references to the term in the literature. It talks about the development of the finite element method, its different versions and how hierarchy has been used to improve the methodology. The next section describes the hierarchical philosophy in detail and explains the different concepts and terms associated with it. It goes on to describe the implementation and the features of the demonstrator. A couple of problems are analyzed using the demonstrator program to show the working of the system. The two problems considered are two dimensional plane stress analysis problems. The results are compared with those obtained using conventional analysis. The different challenges faced during the development of this system are discussed. Finally, we conclude with suggestions for future work to add more features and extend it to a wider range of problems.

Whitcomb, John↗

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗

Demonstration Of Ultra HI-FI (UHF) Methods

Computational aero-acoustics (CAA) requires efficient, high-resolution simulation tools. Most current techniques utilize finite-difference approaches because high order accuracy is considered too difficult or expensive to achieve with finite volume or finite element methods. However, a novel finite volume approach (Ultra HI-FI or UHF) which utilizes Hermite fluxes is presented which can achieve both arbitrary accuracy and fidelity in space and time. The technique can be applied to unstructured grids with some loss of fidelity or with multi-block structured grids for maximum efficiency and resolution. In either paradigm, it is possible to resolve ultra-short waves (less than 2 PPW). This is demonstrated here by solving the 4th CAA workshop Category 1 Problem 1.

Dyson, Rodger W.↗

Planning, creating and documenting a NASTRAN finite element model of a modern helicopter

Mathematical models based on the finite element method of structural analysis as embodied in the NASTRAN computer code are widely used by the helicopter industry to calculate static internal loads and vibration of airframe structure. The internal loads are routinely used for sizing structural members. The vibration predictions are not yet relied on during design. NASA's Langley Research Center sponsored a program to conduct an application of the finite element method with emphasis on predicting structural vibration. The Army/Boeing CH-47D helicopter was used as the modeling subject. The objective was to engender the needed trust in vibration predictions using these models and establish a body of modeling guides which would enable confident future prediction of airframe vibration as part of the regular design process.

Gabal, R.↗

Planning, creating and documenting a NASTRAN finite element model of a modern-helicopter

Mathematical models based on the finite element method of structural analysis as embodied in the NASTRAN computer code are widely used by the helicopter industry to calculate static internal loads and vibration of airframe structure. The internal loads are routinely used for sizing structural members. The vibration predictions are not yet relied on during design. NASA's Langley Research Center sponsored a program to conduct an application of the finite element method with emphasis on predicting structural vibration. The Army/Boeing CH-47D helicopter was used as the modeling subject. The objective was used to engender the needed trust in vibration predictions using these models and establish a body of modeling guides which would enable confident future prediction of airframe vibration as part of the regular design process.

Gabel, R.↗

A Hybrid Numerical Analysis Method for Structural Health Monitoring

A new hybrid surface-integral-finite-element numerical scheme has been developed to model a three-dimensional crack propagating through a thin, multi-layered coating. The finite element method was used to model the physical state of the coating (far field), and the surface integral method was used to model the fatigue crack growth. The two formulations are coupled through the need to satisfy boundary conditions on the crack surface and the external boundary. The coupling is sufficiently weak that the surface integral mesh of the crack surface and the finite element mesh of the uncracked volume can be set up independently. Thus when modeling crack growth, the finite element mesh can remain fixed for the duration of the simulation as the crack mesh is advanced. This method was implemented to evaluate the feasibility of fabricating a structural health monitoring system for real-time detection of surface cracks propagating in engine components. In this work, the authors formulate the hybrid surface-integral-finite-element method and discuss the mechanical issues of implementing a structural health monitoring system in an aircraft engine environment.

Forth, Scott C.↗

Overcoming element quality dependence of finite elements with adaptive extended stencil FEM (AES‐FEM)

Summary The finite element methods (FEMs) are important techniques in engineering for solving partial differential equations, but they depend heavily on element shape quality for stability and good performance. In this paper, we introduce the Adaptive Extended Stencil Finite Element Method (AES‐FEM) as a means for overcoming this dependence on element shape quality. Our method replaces the traditional basis functions with a set of generalized Lagrange polynomial basis functions , which we construct using local weighted least‐squares approximations. The method preserves the theoretical framework of FEM and allows imposing essential boundary conditions and integrating the stiffness matrix in the same way as the classical FEM. In addition, AES‐FEM can use higher‐degree polynomial basis functions than the classical FEM, while virtually preserving the sparsity pattern of the stiffness matrix. We describe the formulation and implementation of AES‐FEM and analyze its consistency and stability. We present numerical experiments in both 2D and 3D for the Poisson equation and a time‐independent convection–diffusion equation. The numerical results demonstrate that AES‐FEM is more accurate than linear FEM, is also more efficient than linear FEM in terms of error versus runtime, and enables much better stability and faster convergence of iterative solvers than linear FEM over poor‐quality meshes. Copyright © 2016 John Wiley & Sons, Ltd.

Conley, Rebecca↗

Control rod treatment for FEM based radiation transport methods

This paper presents a novel control rod treatment for finite element method (FEM) based reactor physics simulations. The new method falls into the class of homogenization methods which treat material interfaces within mesh elements by defining a representative, homogenized nuclear cross section set. In contrast to existing homogenization based cusping treatments that attempt to approximate flux-volume weighting of the cross sections, we propose to exactly integrate terms appearing in the variational form of the transport or diffusion equation over mixed material elements. This goal is achieved by defining smoothly varying cross sections in elements that straddle material boundaries. It is shown that the smoothly varying cross sections can be obtained by projecting the piecewise set of cross sections onto the set of Legendre polynomials up to a finite order. Implementation of smoothly varying cross sections is supported within the Rattlesnake radiation multiphysics code that is currently under development at the Idaho National Laboratory. Homogenization methods traditionally suffer from cusping effects that originate from overestimation of absorption cross sections in strongly neutron-absorbing materials as found in control rods. We demonstrate that the presented method is (1) highly effective in removing cusping effects and (2) provides accurate results based on the C5G7-TD benchmark, a time-dependent version of the C5G7 benchmark.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Finite element applications to explore the effects of partial bonding on metal matrix composite properties

The mechanics of materials approach (definition of E, G, Nu, and Alpha) and the finite element method are used to explore the effects of partial bonding and fiber fracture on the behavior of high temperature metal matrix composites. Composite ply properties are calculated for various degrees of disbonding to evaluate the sensitivity of these properties to the presence of fiber/matrix disbonding and fiber fracture. The mechanics of materials approach allows for the determination of the basic ply material properties needed for design/analysis of composites. The finite element method provides the necessary structural response (forces and displacements) for the mechanics of materials equations. Results show that disbonding of fractured fibers affect only E sub (111) and alpha sub (111) significantly.

Caruso, J. J.↗

Finite element applications to explore the effects of partial bonding on metal matrix composite properties

The mechanics of materials approach (definition of E, G, nu, and alpha) and the finite element method are used to explore the effects of partial bonding and fiber fracture on the behavior of high temperature metal matrix composites. Composite ply properties are calculated for various degrees of disbonding to evaluate the sensitivity of these properties to the presence of fiber/matrix disbonding and fiber fracture. The mechanics of materials approach allows for the determination of the basic ply material properties needed for design/analysis of composites. The finite element method provides the necessary structural response (forces and displacements) for the mechanics of materials equations. Results show that disbonding of fractured fibers affect only E-l(11) and alpha-l(11) significantly.

Caruso, J. J.↗

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗

Reducing Transient Vibrations due to Rotating Shafts

An approach is proposed for calculating the transient torque imbalance of rotating shafts in the diesel generators of autonomous power units. The maximum shear strain of the rotor shaft in the generator is assessed by the finite-element method.

42 ENGINEERING↗