Search NASA⌕ Search

SEARCH · Search NASA

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Development on electromagnetic impedance function modeling and its estimation

Today the Electromagnetic methods such as magnetotellurics (MT) and controlled sources audio MT (CSAMT) is used in a broad variety of applications. Its usefulness in poor seismic areas and its negligible environmental impact are integral parts of effective exploration at minimum cost. As exploration was forced into more difficult areas, the importance of MT and CSAMT, in conjunction with other techniques, has tended to grow continuously. However, there are obviously important and difficult problems remaining to be solved concerning our ability to collect process and interpret MT as well as CSAMT in complex 3D structural environments. This talk aim at reviewing and discussing the recent development on MT as well as CSAMT impedance functions modeling, and also some improvements on estimation procedures for the corresponding impedance functions. In MT impedance modeling, research efforts focus on developing numerical method for computing the impedance functions of three dimensionally (3-D) earth resistivity models. On that reason, 3-D finite elements numerical modeling for the impedances is developed based on edge element method. Whereas, in the CSAMT case, the efforts were focused to accomplish the non-plane wave problem in the corresponding impedance functions. Concerning estimation of MT and CSAMT impedance functions, researches were focused on improving quality of the estimates. On that objective, non-linear regression approach based on the robust M-estimators and the Hilbert transform operating on the causal transfer functions, were used to dealing with outliers (abnormal data) which are frequently superimposed on a normal ambient MT as well as CSAMT noise fields. As validated, the proposed MT impedance modeling method gives acceptable results for standard three dimensional resistivity models. Whilst, the full solution based modeling that accommodate the non-plane wave effect for CSAMT impedances is applied for all measurement zones, including near-, transition-as well as the far-field zones, and consequently the plane wave correction is no longer needed for the impedances. In the resulting robust impedance estimates, outlier contamination is removed and the self consistency between the real and imaginary parts of the impedance estimates is guaranteed. Using synthetic and real MT data, it is shown that the proposed robust estimation methods always yield impedance estimates which are better than the conventional least square (LS) estimation, even under condition of severe noise contamination. A recent development on the constrained robust CSAMT impedance estimation is also discussed. By using synthetic CSAMT data it is demonstrated that the proposed methods can produce usable CSAMT transfer functions for all measurement zones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Variational multiscale element-free Galerkin method combined with the moving Kriging interpolation for solving some partial differential equations with discontinuous solutions

Solving partial differential equations with discontinuous solutions is an important challenging problem in numerical analysis. To this end, there are some methods such as finite volume method, discontinuous Galerkin approach and particle technique that are able to solve these problems. In the current paper, the moving Kriging element-free Galerkin method has been combined with the variational multiscale algorithm to obtain acceptable and high-resolution solutions. For testing this technique, we select some PDEs with discontinuous solution such as Burgers’, Sod’s shock tube, advection–reaction–diffusion, Kuramoto–Sivashinsky, Boussinesq and shallow water equations. First, we obtain a time-discrete scheme by approximating time derivative via finite difference technique. Then we introduce the moving Kriging interpolation and also obtain their shape functions. We use the element-free Galerkin method for approximating the spatial derivatives. This method uses a weak form of the considered equation that is similar to the finite element method with the difference that in the classical element-free Galerkin method test and trial functions are moving least squares approximation (MLS) shape functions. Since the shape functions of moving least squares (MLS) approximation do not have Kronecker delta property, we cannot implement the essential boundary condition, directly. Thus, we employ the shape functions of moving Kriging interpolation and radial point interpolation technique which have the mentioned property. Also, in the element-free Galekin method, we do not use any triangular, quadrangular or other type of meshes. The element-free Galerkin method is a global method while finite element method is a local one. This technique employs a background mesh for integration which makes it different from the truly mesh procedures. Several test problems are solved and numerical simulations are reported which confirm the efficiency of the proposed schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

The evolution history of superalloy powders during hot consolidation and plastic deformation

Highlights: • Element segregations in a powder metallurgy superalloy at different states are identified. • The surface energy drives precipitation of oxides, carbides, and γ' particles on powder surface during hot consolidation. • Prior particle boundary structure is distorted and more uniformly refined grains are obtained after hot extrusion. - Abstract: Powder metallurgy superalloys are widely used to produce critical rotating parts in turbine engines. However, the issues like what is the driving force of precipitation on prior particle boundary (PPB) during consolidation and how the PPB structure evolves during plastic deformation still remain to be resolved. In this study, with the characterization of the microstructure evolution from powder, consolidation, to plastic deformation, the segregation of alloying elements on the particle boundaries at elevated temperature was analyzed by equilibrium thermodynamic calculations, and the distortion of PPB structure under large plastic deformation was quantitively simulated by finite element (FE) methods. The results indicated that surface energy serves as the main driving force of PPB precipitation and hot extrusion could distort the PPB structure and yield more uniformly refined grains.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Predicting Deformation Limits of Dual-Phase Steels Under Complex Loading Paths

The deformation limits of various DP980 steels are examined in this study with deformation instability theory. Under uniaxial tension, overall stress-strain curves of the material are estimated based on simple rule of mixture (ROM) with both iso-strain and iso-stress assumptions. Under complex loading paths, actual microstructure-based finite element (FE) method is used to explicitly resolve the deformation incompatibilities between the soft ferrite and hard martensite phases. The results show that, for uniaxial tension, the deformation instability theory with iso-strain-based ROM can be used to provide the lower bound estimate of the uniform elongation (UE) for the various DP980 considered. Under complex loading paths, the deformation instability theory with microstructure-based FE method can be used in examining the effects of various microstructural features on the deformation limits of DP980 steels.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Galerkin formulation of path integrals in lattice field theory

We present a mathematical framework for Galerkin formulations of path integrals in lattice field theory. The framework is based on using the degrees of freedom (DOFs) associated to a Galerkin discretization as the fundamental lattice variables. We formulate standard concepts in lattice field theory, such as the partition function and correlation functions, in terms of the DOFs. For example, using continuous finite element spaces, we show that the two-point spatial correlation function can be defined between any two points on the domain (as opposed to at just lattice sites) and furthermore, this two-point function satisfies a weak propagator (or Green’s function) identity, in analogy to the continuum case, as well as a convergence estimate obtained from the standard finite element techniques. Furthermore, this framework leads naturally to higher-order formulations of lattice field theories by considering higher-order finite element spaces for the Galerkin discretization. We consider analytical and numerical examples of scalar field theory to investigate how increasing the order of piecewise polynomial finite element spaces affect the approximation of lattice observables. Finally, we sketch an outline of this Galerkin framework in the context of gauge field theories.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Auxiliary Space Preconditioning of Finite Element Equations Using a Nonconforming Interior Penalty Reformulation and Static Condensation

We modify the well-known interior penalty finite element discretization method so that it allows for element-by-element assembly. This is possible due to the introduction of additional unknowns associated with the interfaces between neighboring elements. The resulting bilinear form, and a Schur complement (reduced) version of it, are utilized in a number of auxiliary space preconditioners for the original conforming finite element discretization problem. Furthermore, these preconditioners are analyzed on the fine scale, and their performance is illustrated on model second order scalar elliptic problems discretized with high order elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

PROTEUS-MOC User Manual

The PROTEUS-MOC code is a three-dimensional (3D) neutron transport code based on the finite element mesh and the method of characteristics (MOC) which combines a 2D MOC with the discontinuous Galerkin finite element method for the axial direction. Thus, for PROTEUS-MOC, a 3D geometry is represented with by the implicit extrusion of a single 2D planar geometry with different material assignments allowable on each plane. The code requires four input files for a steady-state calculation without the thermal feedback: a driver input, a mesh input, a cross section input, and a material assignment input. With thermal feedback, a T/H input file is needed. For kinetics, a kinetics driver input file is required. The PROTEUS-MOC software produces a text output file as well as a data file which includes integral and average quantities such on fluxes, powers, temperatures, densities, etc. The detailed output file can be converted to the VTK format using the post processing code so that the outputs can be visualized using a visualization software such as VisIt.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Selective growth of Ni4Ti3 precipitate variants induced by complicated cyclic stress during laser additive manufacturing of NiTi-based composites

Highlights: • The NiTi-based composites were in-situ synthesized by SLM. • Ni{sub 4}Ti{sub 3} could tailor the phase transformation temperature of TiC/NiTi composites. • A three-dimensional stress field modeling of TiC/NiTi composites was established. • The relationship between growth of Ni{sub 4}Ti{sub 3} variants and cyclic stress was disclosed. - Abstract: NiTi-based composites were in-situ synthesized by selective laser melting (SLM) additive manufacturing technology with formation of a great many Ni{sub 4}Ti{sub 3} precipitates. The differential scanning calorimetry (DSC) were carried out to exhibit the similar phase transformation temperature (A{sub s} = 51.82 °C, A{sub f} = 94.02 °C, M{sub s} = 55.87 °C and M{sub f} = 18.93 °C) by comparing with that of SLM-fabricated equi-atomic NiTi alloy, which could be attributed to the compensation of Ni{sub 4}Ti{sub 3} precipitates for the effect of an addition of TiC particles. Three different Ni{sub 4}Ti{sub 3} precipitate variants ((111){sub R}//(1 1-bar 1-bar ){sub B2}, (111){sub R}//( 1-bar 1-bar 1){sub B2}, and (111){sub R}//( 1-bar 1 1-bar ){sub B2}) with various precipitate diameters were observed within SLM fabricated part, showing apparently selective growth behavior. By high-resolution transmission electron microscopy, severe lattice defects formed surrounding the smaller variant, while relatively intact lattice structure was found around the larger one, which indicated that high strain energy was produced in the interface between the matrix and the smaller variant. Furthermore, finite element simulation method was applied to disclose the evolution features of complicated cyclic stress during SLM processing of TiC/NiTi composites. By predicting, the main constraint mode within as-fabricated track was tensile stress effect. Based on the redistribution behavior of thermal stress and the coupling effect between the thermal stress and strain field induced by Ni{sub 4}Ti{sub 3} precipitate, the selective growth mechanism of Ni{sub 4}Ti{sub 3} variant was revealed.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A Modified Isotropic–Kinematic Hardening Model to Predict the Defects in Tube Hydroforming Process

Numerical simulations of tube hydroforming process of hollow crankshafts were conducted by using finite element analysis method. Moreover, the modified model involving the integration of isotropic–kinematic hardening model with ductile criteria model was used to more accurately optimize the process parameters such as internal pressure, feed distance and friction coefficient. Subsequently, hydroforming experiments were performed based on the simulation results. The comparison between experimental and simulation results indicated that the prediction of tube deformation, crack and wrinkle was quite accurate for the tube hydroforming process. Finally, hollow crankshafts with high thickness uniformity were obtained and the thickness distribution between numerical and experimental results was well consistent.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Deformation heterogeneity and intragrain lattice misorientation in high strength contrast, dual-phase bridgmanite/periclase

A bridgmanite/periclase aggregate is investigated due to its prevalence in the Earth’s lower mantle and its importance for understanding geodynamic processes. In dual-phase polycrystalline aggregates, both strength contrast between phases and single crystal elastic and plastic anisotropy are known to influence the development of intragrain misorientation, and the evolution of crystallographic texture. In this study, full-field crystal plasticity simulations are performed with a finite element solution method, and applied to an aggregate of a mechanically strong orthorhombic phase with perovskite structure (bridgmanite, MgSiO 3 ), and a relatively weaker cubic phase (periclase, MgO). The relative strengths of the phases are parameterized to elucidate the effect that strength contrast has on texture evolution and single-phase simulations are performed for comparison. Altogether, results indicate that the relative strength between the two phases influences the development of plasticity, and the overall texture evolution. Results are discussed in light of trends related to the evolution of plasticity and misorientation in the aggregate, and their dependence on both the strength contrast between phases as well as the spatial distribution of grains and phases.

36 MATERIALS SCIENCE↗
DOE OSTI2016

SU-F-J-88: Comparison of Two Deformable Image Registration Algorithms for CT-To-CT Contour Propagation

Purpose: To compare the contour propagation accuracy of two deformable image registration (DIR) algorithms in the Raystation treatment planning system – the “Hybrid” algorithm based on image intensities and anatomical information; and the “Biomechanical” algorithm based on linear anatomical elasticity and finite element modeling. Methods: Both DIR algorithms were used for CT-to-CT deformation for 20 lung radiation therapy patients that underwent treatment plan revisions. Deformation accuracy was evaluated using landmark tracking to measure the target registration error (TRE) and inverse consistency error (ICE). The deformed contours were also evaluated against physician drawn contours using Dice similarity coefficients (DSC). Contour propagation was qualitatively assessed using a visual quality score assigned by physicians, and a refinement quality score (0 0.9 for lungs, > 0.85 for heart, > 0.8 for liver) and similar qualitative assessments (VQS < 0.35, RQS > 0.75 for lungs). When anatomical structures were used to control the deformation, the DSC improved more significantly for the biomechanical DIR compared to the hybrid DIR, while the VQS and RQS improved only for the controlling structures. However, while the inclusion of controlling structures improved the TRE for the hybrid DIR, it increased the TRE for the biomechanical DIR. Conclusion: The hybrid DIR was found to perform slightly better than the biomechanical DIR based on lower TRE while the DSC, VQS, and RQS studies yielded comparable results for both. The use of controlling structures showed considerable improvement in the hybrid DIR results and is recommended for clinical use in contour propagation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Instrument design optimization with computational methods

Using Finite Element Analysis to approximate the solution of differential equations, two different instruments in experimental Hall C at the Thomas Jefferson National Accelerator Facility are analyzed. The time dependence of density uctuations from the liquid hydrogen (LH2) target used in the Q wea k experiment (2011-2012) are studied with Computational Fluid Dynamics (CFD) and the simulation results compared to data from the experiment. The 2.5 kW liquid hydrogen target was the highest power LH2 target in the world and the first to be designed with CFD at Jefferson Lab. The first complete magnetic field simulation of the Super High Momentum Spectrometer (SHMS) is presented with a focus on primary electron beam deflection downstream of the target. The SHMS consists of a superconducting horizontal bending magnet (HB) and three superconducting quadrupole magnets. The HB allows particles scattered at an angle of 5:5 deg to the beam line to be steered into the quadrupole magnets which make up the optics of the spectrometer. Without mitigation, remnant fields from the SHMS may steer the unscattered beam outside of the acceptable envelope on the beam dump and limit beam operations at small scattering angles. A solution is proposed using optimal placement of a minimal amount of shielding iron around the beam line.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

A Simulation Study of a Pixelated CZT Detector Combined Energy-Independent Collimator Using Finite Element and Monte Carlo Methods for SPECT Applications

Compared to other semiconductor materials, Cadmium Zinc Telluride (CZT) is currently considered a leading material for use in clinical SPECT applications [1]. Over the last few decades, simulation and experimental studies related to the charge sharing effect of pixelated CZT detectors have been performed in order to estimate the charge sharing effect between two neighboring pixels and to determine the optimum pixel size of CZT detectors. The goal of this work is to investigate the charge-transport effect in our pixelated CZT detector design using various radionuclides covering a broad range of gamma emission energies from Tl-201 to I-131 that are used for clinical SPECT applications and to estimate physical effects in the detectors combined with our energy-optimized high-sensitivity pixel-matching parallel-hole collimator [2].

36 MATERIALS SCIENCE↗