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 1,351 records · Page 75

All sources NASA NTRS DOE OSTI

DOE OSTI2012

Computation of multi-region relaxed magnetohydrodynamic equilibria

We describe the construction of stepped-pressure equilibria as extrema of a multi-region, relaxed magnetohydrodynamic (MHD) energy functional that combines elements of ideal MHD and Taylor relaxation, and which we call MRXMHD. The model is compatible with Hamiltonian chaos theory and allows the three-dimensional MHD equilibrium problem to be formulated in a well-posed manner suitable for computation. The energy-functional is discretized using a mixed finite-element, Fourier representation for the magnetic vector potential and the equilibrium geometry; and numerical solutions are constructed using the stepped-pressure equilibrium code, SPEC. Convergence studies with respect to radial and Fourier resolution are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Computation of Multi-region Relaxed Magnetohydrodynamic Equilibria

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

A point-centered arbitrary Lagrangian Eulerian hydrodynamic approach for tetrahedral meshes

We present a three dimensional (3D) arbitrary Lagrangian Eulerian (ALE) hydrodynamic scheme suitable for modeling complex compressible flows on tetrahedral meshes. The new approach stores the conserved variables (mass, momentum, and total energy) at the nodes of the mesh and solves the conservation equations on a control volume surrounding the point. This type of an approach is termed a point-centered hydrodynamic (PCH) method. The conservation equations are discretized using an edge-based finite element (FE) approach with linear basis functions. All fluxes in the new approach are calculated at the center of each tetrahedron. A multidirectional Riemann-like problem is solved at the center of the tetrahedron. The advective fluxes are calculated by solving a 1D Riemann problem on each face of the nodal control volume. A 2-stage Runge–Kutta method is used to evolve the solution forward in time, where the advective fluxes are part of the temporal integration. The mesh velocity is smoothed by solving a Laplacian equation. The details of the new ALE hydrodynamic scheme are discussed. Results from a range of numerical test problems are presented.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1996

Computational Analysis of Static and Dynamic Behaviour of Magnetic Suspensions and Magnetic Bearings

Static modelling of magnetic bearings is often carried out using magnetic circuit theory. This theory cannot easily include nonlinear effects such as magnetic saturation or the fringing of flux in air-gaps. Modern computational tools are able to accurately model complex magnetic bearing geometries, provided some care is exercised. In magnetic suspension applications, the magnetic fields are highly three-dimensional and require computational tools for the solution of most problems of interest. The dynamics of a magnetic bearing or magnetic suspension system can be strongly affected by eddy currents. Eddy currents are present whenever a time-varying magnetic flux penetrates a conducting medium. The direction of flow of the eddy current is such as to reduce the rate-of-change of flux. Analytic solutions for eddy currents are available for some simplified geometries, but complex geometries must be solved by computation. It is only in recent years that such computations have been considered truly practical. At NASA Langley Research Center, state-of-the-art finite-element computer codes, 'OPERA', 'TOSCA' and 'ELEKTRA' have recently been installed and applied to the magnetostatic and eddy current problems. This paper reviews results of theoretical analyses which suggest general forms of mathematical models for eddy currents, together with computational results. A simplified circuit-based eddy current model proposed appears to predict the observed trends in the case of large eddy current circuits in conducting non-magnetic material. A much more difficult case is seen to be that of eddy currents in magnetic material, or in non-magnetic material at higher frequencies, due to the lower skin depths. Even here, the dissipative behavior has been shown to yield at least somewhat to linear modelling. Magnetostatic and eddy current computations have been carried out relating to the Annular Suspension and Pointing System, a prototype for a space payload pointing and vibration isolation system, where the magnetic actuator geometry resembles a conventional magnetic bearing. Magnetostatic computations provide estimates of flux density within airgaps and the iron core material, fringing at the pole faces and the net force generated. Eddy current computations provide coil inductance, power dissipation and the phase lag in the magnetic field, all as functions of excitation frequency. Here, the dynamics of the magnetic bearings, notably the rise time of forces with changing currents, are found to be very strongly affected by eddy currents, even at quite low frequencies. Results are also compared to experimental measurements of the performance of a large-gap magnetic suspension system, the Large Angle Magnetic Suspension Test Fixture (LAMSTF). Eddy current effects are again shown to significantly affect the dynamics of the system. Some consideration is given to the ease and accuracy of computation, specifically relating to OPERA/TOSCA/ELEKTRA.

Britcher, Colin P.↗

DOE OSTI2019

Three-dimensional explicit finite element formulation for shear localization with global tracking of embedded weak discontinuities

Shear localization is a damage mechanism characterized by the concentration of plastic deformation within thin bands of materials, often observed as a precursor to failure in metals. In this study, we present a three-dimensional explicit finite element formulation with embedded weak discontinuities for the treatment of shear localization under dynamic loading conditions. In this formulation, the onset of localization is detected via a material stability analysis suitable for rate-sensitive materials, and the continuity of propagating shear bands is ensured using a global tracking strategy involving a heat-conduction type boundary value problem, solved for a scalar level set function over the global domain of the problem. In addition, we propose a modified quadrature rule that is used to compute the contribution of an individual element to the global finite element arrays, taking into account the position of the embedded shear band within that element.

42 ENGINEERING↗

DOE OSTI2018

Optimized Electromagnetic Fields Levitate Bulk Metallic Materials

The bulk metallic electromagnetic levitation technique has a promising application prospect in condensed physics and materials science. In order to resolve the problem of weak levitation force and the small size of samples, an optimized solution based on finite element analysis (FEA) is developed to seek the best electromagnetic field distribution in the levitation zone. The influence of material properties on the distribution of magnetic flux density, surface pressure, and the Lorentz force were studied by the FEA method with various sample sizes of titanium and copper balls. The levitation capability of cylindrical coil system which consists of a single-layer coil with an inner diameter of 30 mm and the levitated object with a diameter of 20 mm was analyzed. To keep a balance between the turn number and increment efficiency of force, the best number of 7 turns is confirmed. Then, two types of optimized coils were predicted by numerical simulations, and the levitation force is improved drastically. The double-layer coil displays an outstanding levitation capability owing to the strong electromagnetic field generated by multilayer turns. In addition, 5 to 50 g aluminum, titanium, copper, and nickel are levitated by the above three kinds of coils and the levitation capacities are validated from the recorded data of the minimum operating current that can levitate the samples. The experimental results show good agreements with the simulation results, and the minimum current for optimized coils is reduced efficiently.

36 MATERIALS SCIENCE↗

DOE OSTI2015

An adhesive contact mechanics formulation based on atomistically induced surface traction

In this work, we have developed a novel multiscale computational contact formulation based on the generalized Derjuguin approximation for continua that are characterized by atomistically enriched constitutive relations in order to study macroscopic interaction between arbitrarily shaped deformable continua. The proposed adhesive contact formulation makes use of the microscopic interaction forces between individual particles in the interacting bodies. In particular, the double-layer volume integral describing the contact interaction (energy, force vector, matrix) is converted into a double-layer surface integral through a mathematically consistent approach that employs the divergence theorem and a special partitioning technique. The proposed contact model is formulated in the nonlinear continuum mechanics framework and implemented using the standard finite element method. With no large penalty constant, the stiffness matrix of the system will in general be well-conditioned, which is of great significance for quasi-static analysis. Three numerical examples are presented to illustrate the capability of the proposed method. Results indicate that with the same mesh configuration, the finite element computation based on the surface integral approach is faster and more accurate than the volume integral based approach. In addition, the proposed approach is energy preserving even in a very long dynamic simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

High-order gradients with the shifted boundary method: An embedded enriched mixed formulation for elliptic PDEs

Here, we propose an extension of the embedded boundary method known as “shifted boundary method” to elliptic diffusion equations in mixed form (e.g., Darcy flow, heat diffusion problems with rough coefficients, etc.). Our aim is to obtain an improved formulation that, for linear finite elements, is at least second-order accurate for both flux and primary variable, when either Dirichlet or Neumann boundary conditions are applied. Following previous work of Nishikawa and Mazaheri in the context of residual distribution methods, we consider the mixed form of the diffusion equation (i.e., with Darcy-type operators), and introduce an enrichment of the primary variable. This enrichment is obtained exploiting the relation between the primary variable and the flux variable, which is explicitly available at nodes in the mixed formulation. The proposed enrichment mimics a formally quadratic pressure approximation, although only nodal unknowns are stored, similar to a linear finite element approximation. We consider both continuous and discontinuous finite element approximations and present two approaches: a non-symmetric enrichment, which, as in the original references, only improves the consistency of the overall method; and a symmetric enrichment, which enables a full error analysis in the classical finite element context. Combined with the shifted boundary method, these two approaches are extended to high-order embedded computations, and enable the approximation of both primary and flux (gradient) variables with second-order accuracy, independently on the type of boundary conditions applied. We also show that the primary variable is third-order accurate, when pure Dirichlet boundary conditions are embedded.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

On optimal control of hybrid dynamical systems using complementarity constraints

Optimal control for switch-based dynamical systems is a challenging problem in the process control literature. In this study, we model these systems as hybrid dynamical systems with finite number of unknown switching points and reformulate them using non-smooth and non-convex complementarity constraints as a mathematical program with complementarity constraints (MPCC). We utilize a moving finite element based strategy to discretize the differential equation system to accurately locate the unknown switching points at the finite element boundary and achieve high-order accuracy at intermediate non-collocation points. We propose a globalization approach to solve the discretized MPCC problem using a mixed NLP/MILP-based strategy to converge to a non-spurious first-order optimal solution. The method is tested on three dynamic optimization examples, including a gas–liquid tank model and an optimal control problem with a sliding mode solution.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

An analytical and numerical study of the nonlinear reflection at a stress-free surface

Implementation of the ultrasonic second harmonic generation has typically been restricted to simple setups such as through-transmission or Rayleigh surface waves. Recent research has evaluated the second harmonic waves generation in P- and SV- waves reflected from a stress-free surface to enable the single-sided interrogation of a specimen. This research considers the second harmonic generation in an aluminum specimen, which is analytically evaluated using an approach based on a perturbation method. Here, the model is chosen to mimic an experimental setup where the longitudinal wave is generated at oblique angle using a wedge transducer. Due to the mode conversion at the interface of the wedge and the specimen, it is necessary to evaluate longitudinal and shear waves, determining all second harmonic waves generated in the bulk and at the stress-free boundary. The theoretically developed model is then implemented in a commercial finite element code, COMSOL, using increasing fundamental wave amplitudes for different values of third order elastic constants. The results of this computational model verify the analytical approach and the proposed measurement setup, taking into account assumptions and approximations of the solution procedure. Furthermore, the computational model is used to draw important conclusions relevant to the experimental setup, including the need to avoid interaction with diffracted waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Planetary normal mode computation: Parallel algorithms, performance and reproducibility

This report is an extension of work entitled “Computing planetary interior normal modes with a highly parallel polynomial filtering eigensolver.” by Shi et al., [1] originally presented at the SC18 conference. A highly parallel polynomial filtered eigensolver was developed and exploited to calculate the planetary normal modes. The proposed method is ideally suited for computing interior eigenpairs for large-scale eigenvalue problems as it greatly enhances memory and computational efficiency. In this article, the second-order finite element method is used to further improve the accuracy as only the first-order finite element method was deployed in the previous work. The parallel algorithm, its parallel performance up to 20k processors, and the great computational accuracy are illustrated. The reproducibility of the previous work was successfully performed on the Student Cluster Competition at the SC19 conference by several participant teams using a completely different Mars-model dataset on different clusters. Both weak and strong scaling performances of the reproducibility by the participant teams were impressive and encouraging. The analysis and reflection of their results are demonstrated and future direction is discussed.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1973

A thermal, thermoelastic, and wear simulation of a high-energy sliding contact problem

This paper describes an investigation of the sliding contact problem encountered in high-energy disk brakes. The analysis includes a simulation modeling, using the finite element method, of the thermoelastic instabilities that cause transient changes in contact to occur on the friction surface. In order to include the effect of wear of the concentrated contacts on the friction surface, a wear criterion is proposed that results in prediction of wear rates for disk brakes that are quite close to experimentally determined wear rates. The thermal analysis shows that the transient temperature distribution in a disk brake can be determined more accurately by use of this thermomechanical analysis than by a more conventional analysis that assumes constant contact conditions. It is also shown that lower, more desirable, temperatures in disk brakes can be attained by increasing the volume, the thermal conductivity, and especially, the heat capacity of the brake components.

Kennedy, F. E., Jr.↗

NASA NTRS1976

Hybrid finite element methods

The purpose of this paper is to show how Lagrange multipliers can be used with finite elements to achieve a number of desirable properties in the underlying approximation. For elliptic boundary value problems, variational principles can be developed in which all boundary conditions are natural. In fluid flow problems, one can endow the approximations with physically essential conservation laws.

Fix, G. M.↗

NASA NTRS1976

Thermal stresses in a spherical pressure vessel having temperature-dependent, transversely isotropic, elastic properties

Rayleigh-Ritz and modified Rayleigh-Ritz procedures are used to construct approximate solutions for the response of a thick-walled sphere to uniform pressure loads and an arbitrary radial temperature distribution. The thermoelastic properties of the sphere are assumed to be transversely isotropic and nonhomogeneous; variations in the elastic stiffness and thermal expansion coefficients are taken to be an arbitrary function of the radial coordinate and temperature. Numerical examples are presented which illustrate the effect of the temperature-dependence upon the thermal stress field. A comparison of the approximate solutions with a finite element analysis indicates that Ritz methods offer a simple, efficient, and relatively accurate approach to the problem.

Tauchert, T. R.↗

NASA NTRS1978

An advanced structural analysis/synthesis capability - ACCESS 2

An advanced automated design procedure for minimum weight design of structures (ACCESS 2) is reported. Design variable linking, constraint deletion, and explicit constraint approximation are used to effectively combine finite element and nonlinear mathematical programming techniques. The approximation concepts approach to structural synthesis is extended to problems involving fiber composite structure, thermal effects and natural frequency constraints in addition to the usual static stress and displacement limitations. Sample results illustrating these new features are given.

Schmit, L. A.↗

NASA NTRS1981

Development of a unified numerical procedure for free vibration analysis of structures

This paper presents the details of a unified numerical algorithm and the associated computer program developed for the efficient determination of natural frequencies and modes of free vibration of structures. Both spinning and nonspinning structures with or without viscous and/or structural damping may be solved by the current procedure; in addition, the program is capable of solving static problems with multiple load cases as well as the quadratic matrix eigenproblem associated with a finite dynamic element structural discretization.

Gupta, K. K.↗

NASA NTRS1981

The stability of pseudospectral-Chebyshev methods

The pseudospectral-Chebyshev methods are shown to be convergent in variable coefficient problems and, in some cases, hyperbolic problems. The analysis demonstrates that the rate of convergence is greater for finite difference methods or the finite element method. For a single first-order hyperbolic equation, the method is seen as remaining stable even when the coefficient changes sign, although in this case it is specified that care must be taken to have adequate spatial resolution. It is noted that this fact, combined with the fact that collocation methods are easy to apply in the nonlinear case, shows that the pseudospectral method is in general preferable to the Galerkin or Tau methods.

Gottlieb, D.↗

NASA NTRS1981

An adaptive, multi-level method for elliptic boundary value problems

Subroutine PLTMG is a FORTRAN program for solving self-adjoint elliptic boundary value problems in general regions of R-squared. It is based on a piecewise linear triangle finite element method, an adaptive grid refinement procedure, and a multi-level iterative method to solve the resulting sets of linear equations. This paper describes the method and presents some numerical results and comparisons.

Bank, R. E.↗