SEARCH · Engineering Papers

Results for “finite element problems”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

High-Order Shock Fitting with Finite Element Methods

A moving-grid, shock-tracking, finite element method has been implemented that can achieve high-order accuracy for flow simulations with shocks. In this approach, element edges in the computational mesh are fitted to the shock front and moved with the shock throughout the simulation. The Euler equations are solved on the moving mesh in an arbitrary Lagrangian-Eulerian framework. Three different methods for specifying the shock motion have been tested,and the order of accuracy of the resulting solutions was verified. It was found that one of the methods is applicable only for problems with a constant upstream flow state, whereas the other two methods converge to the expected solution for general shock problems. When using a finite element method with polynomial degree p, the two general shock motion methods converged with the expected order of accuracy of (p+1). Finally, accurate solutions are obtained for supersonic nozzle and blunt body bow shock problems. Based on the results, it is concluded that the shock tracking method can preserve high-order solution accuracy in the vicinity of a shock wave.

Luke D'Aquila↗

DOE OSTI2018

Optimal Control of a Degenerate PDE for Surface Shape

Controlling the shapes of surfaces provides a novel way to direct self-assembly of colloidal particles on those surfaces and may be useful for material design. This motivates the investigation of an optimal control problem for surface shape in this paper. Specifically, we consider an objective (tracking) functional for surface shape with the prescribed mean curvature equation in graph form as a state constraint. The control variable is the prescribed curvature. We prove existence of an optimal control, and using improved regularity estimates, we show sufficient differentiability to make sense of the first order optimality conditions. This allows us to rigorously compute the gradient of the objective functional for both the continuous and discrete (finite element) formulations of the problem. Numerical results are shown to illustrate the minimizers and optimal controls on different domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Multi-fidelity Uncertainty Quantification for Homogenization Problems in Structure-Property Relationships from Crystal Plasticity Finite Elements

Crystal plasticity finite element method (CPFEM) has been an integrated computational materials engineering (ICME) workhorse to study materials behaviors and structure-property relationships for the last few decades. These relations are mappings from the microstructure space to the materials properties space. Due to the stochastic and random nature of microstructures, there is always some uncertainty associated with materials properties, for example, in homogenized stress-strain curves. For critical applications with strong reliability needs, it is often desirable to quantify the microstructure-induced uncertainty in the context of structure-property relationships. However, this uncertainty quantification (UQ) problem often incurs a large computational cost because many statistically equivalent representative volume elements (SERVEs) are needed. In this article, we apply a multi-level Monte Carlo (MLMC) method to CPFEM to study the uncertainty in stress-strain curves, given an ensemble of SERVEs at multiple mesh resolutions. By using the information at coarse meshes, we show that it is possible to approximate the response at fine meshes with a much reduced computational cost. We focus on problems where the model output is multi-dimensional, which requires us to track multiple quantities of interest (QoIs) at the same time. In conclusion, our numerical results show that MLMC can accelerate UQ tasks around 2.23x, compared to the classical Monte Carlo (MC) method, which is widely known as ensemble average in the CPFEM literature.

36 MATERIALS SCIENCE↗

NASA NTRS1992

On development of a finite dynamic element and solution of associated eigenproblem by a block Lanczos procedure

The paper first presents the details of the development of a new six-noded plane triangular finite dynamic element. A block Lanczos algorithm is developed next for the accurate and efficient solution of the quadratic matrix eigenvalue problem associated with the finite dynamic element formulation. The resulting computer program fully exploits matrix sparsity inherent in such a discretization and proves to be most efficient for the extraction of the usually required first few roots and vectors, including repeated ones. Most importantly, the present eigenproblem solution is shown to be comparable to that of the corresponding finite element analysis, thereby rendering the associated dynamic element method rather attractive owing to superior convergence characteristics of such elements, presented herein.

Gupta, K. K.↗

NASA NTRS1993

Development of an hp-version finite element method for computational optimal control

The purpose of this research effort was to begin the study of the application of hp-version finite elements to the numerical solution of optimal control problems. Under NAG-939, the hybrid MACSYMA/FORTRAN code GENCODE was developed which utilized h-version finite elements to successfully approximate solutions to a wide class of optimal control problems. In that code the means for improvement of the solution was the refinement of the time-discretization mesh. With the extension to hp-version finite elements, the degrees of freedom include both nodal values and extra interior values associated with the unknown states, co-states, and controls, the number of which depends on the order of the shape functions in each element. One possible drawback is the increased computational effort within each element required in implementing hp-version finite elements. We are trying to determine whether this computational effort is sufficiently offset by the reduction in the number of time elements used and improved Newton-Raphson convergence so as to be useful in solving optimal control problems in real time. Because certain of the element interior unknowns can be eliminated at the element level by solving a small set of nonlinear algebraic equations in which the nodal values are taken as given, the scheme may turn out to be especially powerful in a parallel computing environment. A different processor could be assigned to each element. The number of processors, strictly speaking, is not required to be any larger than the number of sub-regions which are free of discontinuities of any kind.

Hodges, Dewey H.↗

DOE OSTI2015

Aleph code electrostatic solver verification

Aleph is an electrostatic particle-in-cell code which uses the finite element method to solve for the electric potential and field based on external potentials and discrete charged particles. The field solver in Aleph was verified for two problems and matched the analytic theory for finite elements. The first problem showed the mesh-refinement convergence for a nonlinear field with no particles within the domain. This matched the theoretical convergence rates of second order for the potential field and first order for the electric field. Then the solution for a single particle in an infinite domain was compared to the analytic solution. This also matched the theory of first order convergence in both the potential and electric fields for both problems over a refinement factor of 16. These solutions give confidence that the field solver and charge weighting schemes are implemented correctly. This page intentionally left blank.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Development and applications of two computational procedures for determining the vibration modes of structural systems

Two computational procedures for analyzing complex structural systems for their natural modes and frequencies of vibration are presented. Both procedures are based on a substructures methodology and both employ the finite-element stiffness method to model the constituent substructures. The first procedure is a direct method based on solving the eigenvalue problem associated with a finite-element representation of the complete structure. The second procedure is a component-mode synthesis scheme in which the vibration modes of the complete structure are synthesized from modes of substructures into which the structure is divided. The analytical basis of the methods contains a combination of features which enhance the generality of the procedures. The computational procedures exhibit a unique utilitarian character with respect to the versatility, computational convenience, and ease of computer implementation. The computational procedures were implemented in two special-purpose computer programs. The results of the application of these programs to several structural configurations are shown and comparisons are made with experiment.

Kvaternik, R. G.↗

DOE OSTI2016

Formulation and numerical analysis of a fully-coupled dynamically deforming electromagnetic wire

Here, an electromagnetic beam model is developed for the simulation of actuated electronic textiles. The beam is solved using a nonlinear director-based kinematic description with additional temperature and electric potential fields along its length. The three fields are fully coupled by mutual dependences on the deformation, Lorenz force, back electromotive force, temperature dependent constitutive responses, and the Seebeck effect. Instead of solving Maxwell's equations in full detail, a quasistatic approximation is used to solve the electric potential in the presence of a moving material medium. The current-carrying beam approximation is used to further simplify the solution space for the potential. While this formulation alleviates the spatial and temporal discretization restrictions, the coupled problem is an index-1 semi-explicit Differential Algebraic Equation requiring special treatment. The time dependent problem is solved using different Runge-Kutta methods. Diagonally implicit Runge-Kutta methods and explicit Runge-Kutta methods using implicit solution of the electric potential problem are explored. The finite element model is implemented using the open source package FEniCS, which is able to automatically generate the linearizations of the multiphysics equations required for the implicit solutions. A model problem is constructed with which to test and analyze the physical formulation and numerical solution techniques. The time stepping methods are verified using the convergence orders of the higher-order Runge-Kutta methods. Runtime comparisons show that the explicit methods are generally more computationally efficient than the implicit schemes used for this problem. For the implicit schemes, a staggered solution is significantly faster than a monolithic solution at most time step sizes. However, at very large time steps, such as those that would be used for dynamic relaxation, the monolithic solution can be more efficient than the staggered solution.

36 MATERIALS SCIENCE↗

NASA NTRS1989

A combined finite element and boundary integral formulation for solution via CGFFT of 2-dimensional scattering problems

A new technique is presented for computing the scattering by 2-D structures of arbitrary composition. The proposed solution approach combines the usual finite element method with the boundary integral equation to formulate a discrete system. This is subsequently solved via the conjugate gradient (CG) algorithm. A particular characteristic of the method is the use of rectangular boundaries to enclose the scatterer. Several of the resulting boundary integrals are therefore convolutions and may be evaluated via the fast Fourier transform (FFT) in the implementation of the CG algorithm. The solution approach offers the principle advantage of having O(N) memory demand and employs a 1-D FFT versus a 2-D FFT as required with a traditional implementation of the CGFFT algorithm. The speed of the proposed solution method is compared with that of the traditional CGFFT algorithm, and results for rectangular bodies are given and shown to be in excellent agreement with the moment method.

Collins, Jeffery D.↗

NASA NTRS1990

Periodic trim solutions with hp-version finite elements in time

Finite elements in time as an alternative strategy for rotorcraft trim problems are studied. The research treats linear flap and linearized flap-lag response both for quasi-trim and trim cases. The connection between Fourier series analysis and hp-finite elements for periodic a problem is also examined. It is proved that Fourier series is a special case of space-time finite elements in which one element is used with a strong displacement formulation. Comparisons are made with respect to accuracy among Fourier analysis, displacement methods, and mixed methods over a variety parameters. The hp trade-off is studied for the periodic trim problem to provide an optimum step size and order of polynomial for a given error criteria. It is found that finite elements in time can outperform Fourier analysis for periodic problems, and for some given error criteria. The mixed method provides better results than does the displacement method.

Peters, David A.↗

NASA NTRS1991

Periodic trim solutions with hp-version finite elements in time

Hou, Lin-Jun↗

NASA NTRS1982

Reliability of the finite element method for calculating free edge stresses in composite laminates

The interlaminar normal stress distributions along the interface between the +45 deg and -45 deg plies of a graphite/epoxy laminate, obtained by various investigators, were found to disagree in both magnitude and sign. The reliability of the displacement-formulated finite element method in analyzing the edge-stress problem of a composite laminate is investigated. The history of the edge-stress problem is reviewed, and two well-known elasticity problems, one involving a stress discontinuity and one a singularity, are analyzed. The finite element analysis in these problems yields accurate stress distributions everywhere except in two elements closest to the stress discontinuity or singularity. Stress distributions for a + or -45 deg ply laminate near the singularity were similar to those of the two elasticity problems, demonstrating the methods, accuracy for calculating interlaminar stresses in composite laminates. The disagreement between the numerical methods was attributed to the unsymmetric stress tensor at singularity.

Whitcomb, J. D.↗

NASA NTRS1993

Development of an hp-version finite element method for computational optimal control

The purpose of this research effort is to develop a means to use, and to ultimately implement, hp-version finite elements in the numerical solution of optimal control problems. The hybrid MACSYMA/FORTRAN code GENCODE was developed which utilized h-version finite elements to successfully approximate solutions to a wide class of optimal control problems. In that code the means for improvement of the solution was the refinement of the time-discretization mesh. With the extension to hp-version finite elements, the degrees of freedom include both nodal values and extra interior values associated with the unknown states, co-states, and controls, the number of which depends on the order of the shape functions in each element.

Hodges, Dewey H.↗

NASA NTRS1990

Finite element two-dimensional panel flutter at high supersonic speeds and elevated temperature

A consistent finite element formulation for flutter of two-dimensional panels including effects of temperature variation along panel length is presented. The von Karman nonlinear strain-displacement relation is used to account for the large-deflections, and the quasi-steady first-order piston theory is employed for aerodynamic loading in the finite-element formulation. The panel flutter under combined thermal/aerodynamic loading can be separated into solving two problems in sequence. They are: (1) postbuckling deflection, stresses and buckling stability boundary, and (2) critical dynamic pressure and linear flutter boundary. The solution procedures in solving these problems are presented in detail. Finite element results of two-dimensional panels with uniform temperature change are compared with classic analytical solutions. Influence of nonuniform temperature distributions on critical buckling temperature, flutter deflection and stresses, and flutter boundary is also given.

Xue, David Y.↗

DOE OSTI2024

SCEPTRE 2.7 User's Guide

Sandia’s Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) is a computer code that solves the linear Boltzmann transport equation, particularly targeting coupled photon-electron problems. It uses unstructured finite element meshes in space, multigroup in energy, and discrete ordinates (Sn) or other methods in angle. SCEPTRE uses an xml-based input file to specify the problem. This report documents the options and syntax of that input file.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

A new numerically stable sequential algorithm for coupled finite-strain elastoplastic geomechanics and flow

In this work, we propose a new numerically stable sequential implicit method for coupled flow and finite-strain multiplicative elastoplastic geomechanics. We find from stability analysis that the sequential method that solves the flow problem first by fixing the first Piola-Kirchhoff total stress, solving the geomechanics (solid deformation) problem at the next step, is unconditionally stable (contractive and B-stable). In this sense, this method named the fixed first Piola-Kirchhoff stress method is an extension of the fixed stress method in coupled flow and infinitesimal geomechanics. We also study the fixed second Piola-Kirchhoff stress method, comparing it with the fixed first Piola-Kirchhoff stress method, because the constitutive relations are formulated by the second Piola-Kirchhoff total stress, although fixing the second Piola-Kirchhoff total stress field does not provide theoretical unconditional stability. In space discretization, we use the finite element method for the geomechanics problem with the total Lagrangian approach, while employing the finite volume method for the flow problem. Geometrical nonlinearity from the total Lagrangian approach results in full-tensor permeability even if the initial permeability is isotropic. To consider the full-tensor permeability accurately, we employ the multipoint flux approximation in solving the flow problem. In time discretization, the backward Euler method is used. Then, we perform numerical experiments of the two fixed-stress sequential methods with various scenarios, and find superiority of their numerical stability. Specifically, the two fixed-stress sequential methods provide almost identical results for elasticity. For J2 plasticity, the two sequential methods yield slightly different numerical results although the trends of the results are similar. For the Drucker-Prager plasticity, we find that the results from the two methods are almost identical. For all the test cases, the two sequential methods are numerically stable.

58 GEOSCIENCES↗

NASA NTRS1982

A wave envelope finite element scheme for acoustical radiation

The aeroacoustic problem associated with the radiation of fan noise from the inlet of a turbofan aircraft engine, the dimensions of which are generally many times larger than the acoustical wavelengths of the major energy-carrying frequencies, is considered. In the present approach, a conventional finite element solution in the inner region is compatibly matched to a 'wave envelope' finite element solution in a large but finite outer region. The inclusion of a wavelike variation with the correct asymptotic decay in the shape functions for the outer region preserves the correct behavior of the solution at large distances. The method is initially presented for a simple one-dimensional model based on the solution of Webster's horn equation. Results are presented for the specific case of a uniform cylindrical section joined to a conical expansion, and also for a simple axisymmetric test case of the calculation of acoustical pressure generated by a vibrating circular piston located at the center of an infinite rigid baffle.

Astley, R. J.↗

DOE OSTI2023

SCEPTRE 2.6 User's Guide

Sandia’s Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) is a computer code that solves the linear Boltzmann transport equation, particularly targeting coupled photon-electron problems. It uses unstructured finite element meshes in space, multigroup in energy, and discrete ordinates (S n ) or other methods in angle. SCEPTRE uses an xml-based input file to specify the problem. This report documents the options and syntax of that input file.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗