Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element Method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

DNS of Flows over Periodic Hills using a Discontinuous-Galerkin Spectral-Element Method

Direct numerical simulation (DNS) of turbulent compressible flows is performed using a higher-order space-time discontinuous-Galerkin finite-element method. The numerical scheme is validated by performing DNS of the evolution of the Taylor-Green vortex and turbulent flow in a channel. The higher-order method is shown to provide increased accuracy relative to low-order methods at a given number of degrees of freedom. The turbulent flow over a periodic array of hills in a channel is simulated at Reynolds number 10,595 using an 8th-order scheme in space and a 4th-order scheme in time. These results are validated against previous large eddy simulation (LES) results. A preliminary analysis provides insight into how these detailed simulations can be used to improve Reynoldsaveraged Navier-Stokes (RANS) modeling

DNS↗
DOE OSTI2018

Multi-objective design of microvascular panels for battery cooling applications

Building on a recently developed optimization method based on an interface-enriched generalized finite element method, multiple objective functions are considered for the optimization of 2D networks of microchannels embedded in battery-cooling panels. The objective functions considered in this study are a differentiable alternative to the maximum temperature (the p-mean of the temperature), the pressure drop and the variance of the temperature. The -constraint method and the normalized normal constraint method are used to generate the pressure-temperature Pareto optimal front of the multi-objective optimization problem. The effects of different operating constraints/conditions such as localization of heat sources, prescribed pump power and imposed flow rate on the optimal designs are investigated. In addition to the topology of the embedded network, the cross sections of the microchannels are also introduced as design parameters to further improve the pressure drop of the designs. The resulting variable-cross-section optimized design is validated with experiment.

energy storage (including batteries and capacitors↗
DOE OSTI2021

Machine learning based inverse modeling of full-field strain distribution for mechanical characterization of a linear elastic and heterogeneous membrane

Heterogeneous membranes or films are thin and soft structures with spatial variations in material property and thickness. Mechanical behavior of heterogeneous membranes is not well understood, mainly due to the difficulty in obtaining accurate and reliable material property data. To understand the mechanical behavior of these materials, accurate and efficient characterization methods for heterogeneous membranes are needed. Here, in this paper, an inverse method based on machine learning is developed to efficiently extract mechanical properties from full-field strain distributions. This approach is demonstrated on a flat heterogeneous membrane with uniform thickness formed by up to four linear elastic synthetic materials in a grid arrangement, and deforming in a moderate strain range (true strain ~10%). The results show that the machine learning method achieves accuracy comparable to the traditional inverse finite element method, and is 6 orders of magnitude faster in the demonstrated case studies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗
NASA NTRS1985

Analysis and test evaluation of the dynamic stability of three advanced turboprop models at zero forward speed

Results of static stability wind tunnel tests of three 62.2 cm (24.5 in) diameter models of the Prop-Fan are presented. Measurements of blade stresses were made with the Prop-Fans mounted on an isolated nacelle in an open 5.5 m (18 ft) wind tunnel test section with no tunnel flow. The tests were conducted in the United Technology Research Center Large Subsonic Wind Tunnel. Stall flutter was determined by regions of high stress, which were compared with predictions of boundaries of zero total viscous damping. The structural analysis used beam methods for the model with straight blades and finite element methods for the models with swept blades. Increasing blade sweep tends to suppress stall flutter. Comparisons with similar test data acquired at NASA/Lewis are good. Correlations between measured and predicted critical speeds for all the models are good. The trend of increased stability with increased blade sweep is well predicted. Calculated flutter boundaries generaly coincide with tested boundaries. Stall flutter is predicted to occur in the third (torsion) mode. The straight blade test shows third mode response, while the swept blades respond in other modes.

Smith, Arthur F.↗
NASA NTRS1989

Comparison of the measured and predicted response of the Earth Radiation Budget Experiment active cavity radiometer during solar observations

The Earth Radiation Budget Experiment consists of an array of radiometric instruments placed in earth orbit by the National Aeronautics and Space Administration to monitor the longwave and visible components of the earth's radiation budget. Presented is a dynamic electrothermal model of the active cavity radiometer used to measure the earth's total radiative exitance. Radiative exchange is modeled using the Monte Carlo method and transient conduction is treated using the finite element method. Also included is the feedback circuit which controls electrical substitution heating of the cavity. The model is shown to accurately predict the dynamic response of the instrument during solar calibration.

Mahan, J. R.↗
NASA NTRS1991

Dynamic analysis of truss-beam system

A simple truss-beam method for determining the dynamic characteristics of the space structures intended to perform various tasks in orbit is presented. Algorithms are provided to determine the flexibility matrix of the truss beam for use in the equation of motion. The natural frequencies obtained through this method are compared with those obtained through the finite element method. An experimental procedure for verifying the theoretical results is described. It is concluded that the truss-beam method is a simple analysis technique that yields reasonably accurate results with a minimum of computational effort, especially for the critical lower bending modes of flexible structures. The method takes less computer time than more conventional methods and can be programed in both BASIC and FORTRAN for use in micro- and mainframe computers.

Abu-Saba, Elias G.↗
NASA NTRS1992

A SINDA thermal model using CAD/CAE technologies

The approach to thermal analysis described by this paper is a technique that incorporates Computer Aided Design (CAD) and Computer Aided Engineering (CAE) to develop a thermal model that has the advantages of Finite Element Methods (FEM) without abandoning the unique advantages of Finite Difference Methods (FDM) in the analysis of thermal systems. The incorporation of existing CAD geometry, the powerful use of a pre and post processor and the ability to do interdisciplinary analysis, will be described.

Rodriguez, Jose A.↗
NASA NTRS1995

P-Finite-Element Program For Analysis Of Plates

BUCKY is p-finite-element computer program for highly accurate analysis of structures. Used to analyze buckling, bending, and in-plane stress-and-strain behaviors of plates. Provides elastic-plastic solutions for isotropic plates in states of plane stress, and axisymmetric solution sequence used to treat three-dimensional problems. Computes response of plate to variety of loading and boundary conditions by use of higher-order displacement function in p-finite-element method. Enables user to obtain results more accurate than obtained by use of traditional h-finite elements. Written in FORTRAN 77.

Smith, James P.↗
NASA NTRS2004

Electron-Hydrogen Elastic Scattering

Scattering by single-electron systems is always of interest because the wave function of the target is known exactly. Various approximations have been employed to take into account distortion produced in the target. Among them are the method of polarized orbitals and the close coupling approximation. Recently, e-H and e-He+ S-wave scattering in the elastic region has been studied using the Feshbach projection operator formalism. In this approach, the usual Hartree-Fock and exchange potentials are augmented by an optical potential and the resulting phase shifts have rigorous lower bounds. Now this method is being applied to the e-H P-wave scattering in the elastic region. The number of terms in the Hylleraas-type wave function for the 1,3 P phase shifts is 84 and the resulting phase shifts (preliminary) are given. The results have been given up to five digits because to that accuracy they are rigorous lower bounds. They are in general agreement with the variational (VAR) results of Armstead, and those obtained from the intermediate energy R-matrix method (RM) of Scholz et al., and the finite element method (FEM) of Botero and Shertzer. The later two methods do not provide any bounds on phase shifts.

Bhatia, A. K.↗
NASA NTRS1977

A finite element simulation of sound attenuation in a finite duct with a peripherally variable liner

Using multimodal analysis, a variational finite element method is presented for analyzing sound attenuation in a three-dimensional finite duct with a peripherally variable liner in the absence of flow. A rectangular element, with cubic shaped functions, is employed. Once a small portion of a peripheral liner is removed, the attenuation rate near the frequency where maximum attenuation occurs drops significantly. The positioning of the liner segments affects the attenuation characteristics of the liner. Effects of the duct termination are important in the low frequency ranges. The main effect of peripheral variation of the liner is a broadening of the attenuation characteristics in the midfrequency range. Because of matrix size limitations of the presently available computer program, the eigenvalue equations should be solved out of core in order to handle realistic sources.

Watson, W. R.↗
DOE OSTI2015

Temperature Profile in Fuel and Tie-Tubes for Nuclear Thermal Propulsion Systems

A finite element method to calculate temperature profiles in heterogeneous geometries of tie-tube moderated LEU nuclear thermal propulsion systems and HEU designs with tie-tubes is developed and implemented in MATLAB. This new method is compared to previous methods to demonstrate shortcomings in those methods. Typical methods to analyze peak fuel centerline temperature in hexagonal geometries rely on spatial homogenization to derive an analytical expression. These methods are not applicable to cores with tie-tube elements because conduction to tie-tubes cannot be accurately modeled with the homogenized models. The fuel centerline temperature directly impacts safety and performance so it must be predicted carefully. The temperature profile in tie-tubes is also important when high temperatures are expected in the fuel because conduction to the tie-tubes may cause melting in tie-tubes, which may set maximum allowable performance. Estimations of maximum tie-tube temperature can be found from equivalent tube methods, however this method tends to be approximate and overly conservative. A finite element model of heat conduction on a unit cell can model spatial dependence and non-linear conductivity for fuel and tie-tube systems allowing for higher design fidelity of Nuclear Thermal Propulsion.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2016

Energy dissipation via acoustic emission in ductile crack initiation

In this article, we present a modeling approach to estimate the energy release due to ductile crack initiation in conjunction to the energy dissipation associated with the formation and propagation of transient stress waves typically referred to as acoustic emission. To achieve this goal, a ductile fracture problem is investigated computationally using the finite element method based on a compact tension geometry under Mode I loading conditions. To quantify the energy dissipation associated with acoustic emission, a crack increment is produced given a pre-determined notch size in a 3D cohesive-based extended finite element model. The computational modeling methodology consists of defining a damage initiation state from static simulations and linking such state to a dynamic formulation used to evaluate wave propagation and related energy redistribution effects. The model relies on a custom traction separation law constructed using full field deformation measurements obtained experimentally using the digital image correlation method. The amount of energy release due to the investigated first crack increment is evaluated through three different approaches both for verification purposes and to produce an estimate of the portion of the energy that radiates away from the crack source in the form of transient waves. The results presented herein propose an upper bound for the energy dissipation associated to acoustic emission, which could assist the interpretation and implementation of relevant nondestructive evaluation methods and the further enrichment of the understanding of effects associated with fracture.

42 ENGINEERING↗
DOE OSTI2015

A Novel Transverse Flux Machine for Vehicle Traction Applications

A novel transverse flux machine topology for electric vehicle traction application using ferrite magnets is presented in this paper. The proposed transverse flux topology utilizes novel magnet arrangements in the rotor that are similar to Halbach-array to boost flux linkage; on the stator side, cores are alternately arranged around a pair of ring windings in each phase to make use of the entire rotor flux that eliminates end windings. Analytical design considerations and finite element methods are used for an optimized design of a scooter in-wheel motor. Simulation results from Finite Element Analysis (FEA) show the motor achieved comparable torque density to conventional rare-earth permanent magnet machines. This machine is a viable candidate for direct drive applications with low cost and high torque density.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1988

Nonconforming mortar element methods: Application to spectral discretizations

Spectral element methods are p-type weighted residual techniques for partial differential equations that combine the generality of finite element methods with the accuracy of spectral methods. Presented here is a new nonconforming discretization which greatly improves the flexibility of the spectral element approach as regards automatic mesh generation and non-propagating local mesh refinement. The method is based on the introduction of an auxiliary mortar trace space, and constitutes a new approach to discretization-driven domain decomposition characterized by a clean decoupling of the local, structure-preserving residual evaluations and the transmission of boundary and continuity conditions. The flexibility of the mortar method is illustrated by several nonconforming adaptive Navier-Stokes calculations in complex geometry.

Maday, Yvon↗
DOE OSTI2015

Novel Transverse Flux Machine for Vehicle Traction Applications: Preprint

A novel transverse flux machine topology for electric vehicle traction applications using ferrite magnets is presented in this paper. The proposed transverse flux topology utilizes novel magnet arrangements in the rotor that are similar to the Halbach array to boost flux linkage; on the stator side, cores are alternately arranged around a pair of ring windings in each phase to make use of the entire rotor flux that eliminates end windings. Analytical design considerations and finite-element methods are used for an optimized design of a scooter in-wheel motor. Simulation results from finite element analysis (FEA) show that the motor achieved comparable torque density to conventional rare-earth permanent magnet (PM) machines. This machine is a viable candidate for direct-drive applications with low cost and high torque density.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

Surface Remeshing with Robust High-Order Reconstruction

Remeshing is an important problem in variety of applications, such as finite element methods and geometry processing. Surface remeshing poses some unique challenges, as it must deliver not only good mesh quality but also good geometric accuracy. For applications such as finite elements with high-order elements (quadratic or cubic elements), the geometry must be preserved to high-order (third-order or higher) accuracy, since low-order accuracy may undermine the convergence of numerical computations. The problem is particularly challenging if the CAD model is not available for the underlying geometry, and is even more so if the surface meshes contain some inverted elements. We describe remeshing strategies that can simultaneously produce high-quality triangular meshes, untangling mildly folded triangles and preserve the geometry to high-order of accuracy. Our approach extends our earlier works on high-order surface reconstruction and mesh optimization by enhancing its robustness with a geometric limiter for under-resolved geometries. We also integrate high-order surface reconstruction with surface mesh adaptation techniques, which alter the number of triangles and nodes. We demonstrate the utilization of our method to meshes for high-order finite elements, biomedical image-based surface meshes, and complex interface meshes in fluid simulations.

Ray, Navamita↗
NASA NTRS1986

Unstructured multigrid methods for the compressible Euler equations

A two-step explicit finite element based solution procedure for the compressible Euler equations is described. Convergence to steady state can be accelerated by using a multigrid technique. The geometric flexibility of the finite element method is retained by using a sequence of unnested grids. The viability of this approach is demonstrated for elliptic problems and initial experiences with the technique applied to the Euler equations are reported.

Loehner, R.↗