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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Modelling Sawing of Metal Tubes Through FEM Simulation

The paper presents the development of a numerical model of the sawing process of AISI 304 thin tubes, which is cut through a circular blade with alternating roughing and finishing teeth. The numerical simulation environment is the three-dimensional FEM software Deform v.10.1. The teeth actual trajectories were determined by a blade kinematics analysis developed in Matlab. Due to the manufacturing rolling steps and subsequent welding stage, the tube material is characterized by a gradient of properties along its thickness. Consequently, a simplified cutting test was set up and carried out in order to identify the values of relevant material parameters to be used in the numerical model. The dedicated test was the Orthogonal Tube Cutting test (OTC), which was performed on an instrumented lathe. The proposed numerical model was validated by comparing numerical results and experimental data obtained from sawing tests carried out on an industrial machine. The following outputs were compared: the cutting force, the chip thickness, and the chip contact area.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Fourier Analyses of High-Order Continuous and Discontinuous Galerkin Methods

In this paper, we present a Fourier analysis of wave propagation problems subject to a class of continuous and discontinuous discretizations using high-degree Lagrange polynomials. This allows us to obtain explicit analytical formulas for the dispersion relation and group velocity and, for the first time to our knowledge, characterize analytically the emergence of gaps in the dispersion relation at specific wavenumbers, when they exist, and compute their specific locations. Wave packets with energy at these wavenumbers will fail to propagate correctly, leading to significant numerical dispersion. We also show that the Fourier analysis generates mathematical artifacts, and we explain how to remove them through a branch selection procedure conducted by analysis of eigenvectors and associated reconstructed solutions. The higher frequency eigenmodes, named erratic in this study, are also investigated analytically and numerically.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Dynamics, Radiation, and Overall Energy Budget of Earthquake Rupture With Coseismic Off-Fault Damage

Earthquake ruptures dynamically activate coseismic off–fault damage around fault cores. Systematic field observation efforts have shown the distribution of off–fault damage around main faults, while numerical modeling using elastic–plastic off–fault material models have demonstrated the evolution of coseismic off–fault damage during earthquake ruptures. Laboratory–scale microearthquake experiments have pointed out the enhanced high–frequency radiation due to the coseismic off–fault damage. However, the detailed off–fault fracturing mechanisms, subsequent radiation, and its contribution to the overall energy budget remain to be fully understood because of limitations of current observational techniques and model formulations. Here, we constructed a new physics–based dynamic earthquake rupture modeling framework, based on the combined finite–discrete element method, to investigate the fundamental mechanisms of coseismic off–fault damage, and its effect on the rupture dynamics, the radiation and the overall energy budget. We conducted a 2–D systematic case study with depth and showed the mechanisms of dynamic activation of the coseismic off–fault damage. We found the decrease in rupture velocity and the enhanced high–frequency radiation in near field due to the coseismic off–fault damage. We then evaluated the overall energy budget, which shows a significant contribution of the coseismic off–fault damage to the overall energy budget even at depth, where the damage zone width becomes narrower. Here, the present numerical framework for the dynamic earthquake rupture modeling thus provides new insights into earthquake rupture dynamics with the coseismic off–fault damage.

58 GEOSCIENCES↗
DOE OSTI2011

A Modified Johnson-Cook Model for Sheet Metal Forming at Elevated Temperatures and Its Application for Cooled Stress-Strain Curve and Spring-Back Prediction

In this study, a modified Johnson-Cook (J-C) model and an innovated method to determine (J-C) material parameters are proposed to predict more correctly stress-strain curve for tensile tests in elevated temperatures. A MATLAB tool is used to determine material parameters by fitting a curve to follow Ludwick's hardening law at various elevated temperatures. Those hardening law parameters are then utilized to determine modified (J-C) model material parameters. The modified (J-C) model shows the better prediction compared to the conventional one. As the first verification, an FEM tensile test simulation based on the isotropic hardening model for boron sheet steel at elevated temperatures was carried out via a user-material subroutine, using an explicit finite element code, and compared with the measurements. The temperature decrease of all elements due to the air cooling process was then calculated when considering the modified (J-C) model and coded to VUMAT subroutine for tensile test simulation of cooling process. The modified (J-C) model showed the good agreement between the simulation results and the corresponding experiments. The second investigation was applied for V-bending spring-back prediction of magnesium alloy sheets at elevated temperatures. Here, the combination of proposed J-C model with modified hardening law considering the unusual plastic behaviour for magnesium alloy sheet was adopted for FEM simulation of V-bending spring-back prediction and shown the good comparability with corresponding experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The estimation of Poisson’s ratio by time-averaging and Cornu’s method for isotropic beams

Prior research suggests that direct static (e.g., uniaxial testing) and dynamic (e.g., ultrasonic wave propagation analysis) measurements differ in the estimation of Poisson’s ratio because anisotropies and heterogeneities in the sample material affect the two types of tests differently. Even assuming isotropic and homogeneous material properties, prior research further suggests that discrepancies between static/dynamic test results will exist because the error of the diagnostic techniques for the measurand are inherently different. Finally, thermodynamic effects are not present in static tests but can significantly affect dynamic test results. Given the potential for all these variables to produce discrepancies, it would be helpful to have the measurement of Poisson’s ratio obtainable from the same theory and experimental measurements by either static or dynamic testing methods. Our finite element calculations show that by combining time-averaged scanning digital holography with Cornu’s method, it is theoretically possible to estimate the effective Poisson’s ratio from the anticlastic contours at the antinode of the first out-of-plane bending mode shape. This is true regardless of frequency and therefore applicable for both static and dynamic measurements. Our results show that the estimate of Poisson’s ratio by Cornu’s method using data from simulations of mode shapes approaches the true value of Poisson’s ratio. Additionally, our research suggests that beam geometry and boundary conditions are fundamental factors limiting the convergence of the estimate of Poisson’s ratio to the true value of Poisson’s ratio regardless of performing a static or dynamic test.

42 ENGINEERING↗
DOE OSTI2017

BeamDyn: a high-fidelity wind turbine blade solver in the FAST modular framework

Here, this paper presents a numerical implementation of the geometrically exact beam theory based on the Legendre-spectral-finite-element (LSFE) method. The displacement-based geometrically exact beam theory is presented, and the special treatment of three-dimensional rotation parameters is reviewed. An LSFE is a high-order finite element with nodes located at the Gauss-Legendre-Lobatto points. These elements can be an order of magnitude more computationally efficient than low-order finite elements for a given accuracy level. The new module, BeamDyn, is implemented in the FAST modularization framework for dynamic simulation of highly flexible composite-material wind turbine blades within the FAST aeroelastic engineering model. The framework allows for fully interactive simulations of turbine blades in operating conditions. Numerical examples are provided to validate BeamDyn and examine the LSFE performance as well as the coupling algorithm in the FAST modularization framework. BeamDyn can also be used as a stand-alone high-fidelity beam tool.

17 WIND ENERGY↗
DOE OSTI2011

Nanoscale silicon-on-insulator deformation induced by stressed liner structures

Rotation and strain fields were mapped across silicon-on-insulator (SOI) regions induced by overlying stressed Si{sub 3}N{sub 4} features using x-ray nanobeam diffraction. The distribution in SOI tilt exhibited an antisymmetric distribution with a maximum magnitude of 7.9 milliradians, representing one of the first direct measurements of the lattice tilt conducted in situ within buried layers using a spot size of less than 100 nm. The measured rotation distribution corresponds to simulated values generated by boundary element method modeling, indicating that the strain transfer into the underlying SOI primarily induces elastic deformation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Formability analysis of aluminum alloy sheets at elevated temperatures with numerical simulation based on the M-K method

The effect of temperature on formability of aluminum alloy sheet is studied by developing the Forming Limit Diagrams, FLD, for aluminum alloy 3000-series using the Marciniak and Kuczynski technique by numerical simulation. The numerical model is conducted in LS-DYNA and incorporates the Barlat's YLD2000 anisotropic yield function and the temperature dependant Bergstrom hardening law. Three different temperatures; room temperature, 250 deg. C and 300 deg. C, are studied. For each temperature case, various loading conditions are applied to the M-K defect model. The effect of the material anisotropy is considered by varying the defect angle. A simplified failure criterion is used to predict the onset of necking. Minor and major strains are obtained from the simulations and plotted for each temperature level. It is demonstrated that temperature improves the forming limit of aluminum 3000-series alloy sheet.

36 MATERIALS SCIENCE↗
DOE OSTI2014

SU-F-19A-11: Retrospective Evaluation of Thermal Coverage by Thermobrachytherapy Seed Arrangements of Clinical LDR Prostate Implants

Purpose: The superiority of a properly-administered combination of radiation therapy and hyperthermia over radiation alone in treatment of human cancers has been demonstrated in multiple studies examining radiobiology, local control, and survival. Unfortunately, hyperthermia is not yet a common modality in oncology practice, due in part to the technical difficulty of heating a deep-seated target volume to sufficient temperature. To address this problem, our group has invented a thermobrachytherapy (TB) seed based on a commonly-used low dose-rate permanent brachytherapy seed for implant in solid tumors. Instead of the tungsten radiographic marker of the standard seed, the TB seed contains one of a self-regulating ferromagnetic alloy. Placement of a patient implanted with such seeds in an oscillating magnetic field generates heat via induction of eddy currents. We present the results of studies of the capability of clinically-realistic TB seed arrangements to adequately heat defined target volumes. Methods: Seed distributions for several past LDR prostate permanent implant brachytherapy patients were reproduced in the finite element analysis software package COMSOL Multiphysics 4.4, with the difference that TB seeds were modelled, rather than the radiation-only seeds actually used for their treatments. The implant geometries were mainly of the modified peripheral loading type; a range of prostatic volumes and blood perfusion rates likely to be seen in a clinical setting were examined. Results: According to the simulations, when distributed to optimize radiation dose, TB seeds also produce sufficient heat to provide thermal coverage of the target given proper selection of the magnetic field strength. However, the thermal distributions may be improved by additional use of hyperthermia-only seeds. Conclusion: A dual-modality seed intended as an alternative to and using the same implantation apparatus and technique as the standard LDR permanent implant seed has been successfully evaluated for its ability to provide sufficient hyperthermia in clinically-realistic implants. This work was partially supported by the National Institutes of Health (NIH) STTR Grant No. R41 CA153631-01A1.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Two-Level Least Squares Methods in Krylov Subspaces

Two-level least squares acceleration approaches are applied to the Chebyshev acceleration method and the restarted conjugate residual method in solving systems of linear algebraic equations with sparse unsymmetric coefficient matrices arising from finite volume or finite element approximations of boundary-value problems on irregular grids. Application of the proposed idea to other iterative restarted processes also is considered. The efficiency of the algorithms suggested is investigated numerically on a set of model Dirichlet problems for the convection-diffusion equation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Machine learning simulation of finite element analysis in augmented reality

Media, method and system for approximating a finite element analysis texture map for an object. To accomplish this, the object is converted to a computer generated model and finite element analysis is performed for a plurality of different simulated inputs to generate a plurality of simulated mappings. Each simulated mapping is converted into a simulated texture map. A machine learning model is trained on the simulated inputs and simulated texture maps to generate a texture map which approximates a finite element analysis. The machine learning model receives a user input and generates the texture map therefrom. The texture map is then wrapped to the object and displayed.

Scherer, Derek Carl↗
DOE OSTI2017

Multislice does it all—calculating the performance of nanofocusing X-ray optics

Here, we describe an approach to calculating the optical performance of a wide range of nanofocusing X-ray optics using multislice scalar wave propagation with a complex X-ray refractive index. This approach produces results indistinguishable from methods such as coupled wave theory, and it allows one to reproduce other X-ray optical phenomena such as grazing incidence reflectivity where the direction of energy flow is changed significantly. Just as finite element analysis methods allow engineers to compute the thermal and mechanical responses of arbitrary structures too complex to model by analytical approaches, multislice propagation can be used to understand the properties of the real-world optics of finite extent and with local imperfections, allowing one to better understand the limits to nanoscale X-ray imaging.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Microstructural effects in the simulation of creep of concrete

The influence of the microstructure on the visco-elastic properties of concrete is investigated using finite element simulations at the meso-scale. We first derive a constitutive law for the creep of the cement paste which accounts for both recoverable and permanent deformations, as well as the influence of temperature and internal relative humidity. The model is calibrated on a set of experiments at the cement paste scale, and then validated after upscaling to the concrete scale. The model is then applied to study the influence of the microstructure on the macroscopic creep of concrete. We show that materials with finer particles exhibit less creep, and that the anisotropy of creep can be explained with the shape and orientation of the aggregates. Furthermore, the acceleration of the stress relaxation in the presence of damage is explained by the micro-mechanical interaction between the aggregates, the cement paste, and the micro-cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Coupling extended magnetohydrodynamic fluid codes with radiofrequency ray tracing codes for fusion modeling

Neoclassical tearing modes are macroscopic ($L$ ~ 1 m) instabilities in magnetic fusion experiments; if unchecked, these modes degrade plasma performance and may catastrophically destroy plasma confinement by inducing a disruption. Fortunately, the use of properly tuned and directed radiofrequency waves ($λ$ ~ 1 mm) can eliminate these modes. Numerical modeling of this difficult multiscale problem requires the integration of separate mathematical models for each length and time scale [T. G. Jenkins, S. E. Kruger, Phys. Plasmas 19 (2012) 122508]; the extended MHD model captures macroscopic plasma evolution while the RF model tracks the flow and deposition of injected RF power through the evolving plasma profiles. The scale separation enables use of the eikonal (ray-tracing) approximation to model the RF wave propagation. In this work we demonstrate a technique, based on methods of computational geometry, for mapping the ensuing RF data (associated with discrete ray trajectories) onto the finite-element/pseudospectral grid that is used to model the extended MHD physics. In the new representation, the RF data can then be used to construct source terms in the equations of the extended MHD model, enabling quantitative modeling of RF-induced tearing mode stabilization. Though our specific implementation uses the NIMROD extended MHD [C. R. Sovinec et al., J. Comp. Phys. 195 (2004) 355] and GENRAY RF [A. P. Smirnov, R. W. Harvey, K. Kupfer, Bull. Amer. Phys. Soc. 39 (1994) 1626] codes, the approach presented can be applied more generally to any code coupling requiring the mapping of ray tracing data onto Eulerian grids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Relaxation Phenomena of Binary Gas-Mixture with Sizable Difference in their Temperatures and Velocities on the basis of Boltzmann Equations

Relaxation phenomena in the binary gas-mixture with different temperature and different velocities are discussed on the basis of two Boltzmann equations. The Hermite expansion method, extended by H.Grad to multidimensional space, is applied to express distribution functions and the Galerkin method is used to solve two Boltzmann equations. Thus, a system of differential equations for the expansion coefficients is obtained. The time development of the system is calculated numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Characterization Of Flow Stress Of Different AA6082 Alloys By Means Of Hot Torsion Test

FEM simulations are become the most powerful tools in order to optimize the different aspects of the extrusion process and an accurate flow stress definition of the alloy is a prerequisite for a reliable effectiveness of the simulation. In the paper the determination of flow stress by means of hot torsion test is initially presented and discussed: the several approximations that are usually introduced in flow stress computation are described and computed for an AA6082 alloy in order to evidence the final effect on curves shapes. The procedure for regressing the parameters of the sinhyperbolic flow stress definition is described in detailed and applied to the described results. Then four different alloys, extracted by different casting batches but all namely belonging to the 6082 class, were hot torsion tested in comparable levels of temperature and strain rate up to specimen failure. The results are analyzed and discussed in order to understand if a mean flow stress behavior can be identified for the whole material class at the different tested conditions or if specific testing conditions (chemical composition of the alloy, specimen shape, etc) influence the materials properties to a higher degree.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Numerical simulation of linear fiction welding (LFW) processes

Solid state welding processes are becoming increasingly important due to a large number of advantages related to joining ''unweldable'' materials and in particular light weight alloys. Linear friction welding (LFW) has been used successfully to bond non-axisymmetric components of a range of materials including titanium alloys, steels, aluminum alloys, nickel, copper, and also dissimilar material combinations. The technique is useful in the research of quality of the joints and in reducing costs of components and parts of the aeronautic and automotive industries.LFW involves parts to be welded through the relative reciprocating motion of two components under an axial force. In such process the heat source is given by the frictional forces work decaying into heat determining a local softening of the material and proper bonding conditions due to both the temperature increase and the local pressure of the two edges to be welded. This paper is a comparative test between the numerical model in two dimensions, i.e. in plane strain conditions, and in three dimensions of a LFW process of AISI1045 steel specimens. It must be observed that the 3D model assures a faithful simulation of the actual threedimensional material flow, even if the two-dimensional simulation computational times are very short, a few hours instead of several ones as the 3D model. The obtained results were compared with experimental values found out in the scientific literature.

36 MATERIALS SCIENCE↗