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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Design and Evaluation of an Experimental Technique for Mechanical and Fatigue Testing of Sub-Sized Samples

Some recently developed processing techniques allow fabrication of small specimens. The fatigue behavior of these small samples cannot be characterised according to the existing standards due to size limitations. The main objective of the present work is to develop an original technique for both mechanical and fatigue testing of sub-sized specimens produced by a high pressure torsion (HPT) method. The fatigue behavior of the sub-sized Ti-6Al-4 V samples, made from HPT-processed disks, is compared with that of the standard samples prepared according to the ASTM standards. Finite element method (FEM) analysis of the stress distribution under applied static load, for both sub-sized specimens and grips designed for their testing, is performed. Good agreement between the outcomes of tensile testing of the standard ASTM samples and sub-sized samples is observed. Fatigue data, represented by Wöhler curves for both types of samples, showed significant divergence of high cycle fatigue properties caused by incompliance of specimen geometries and surface/defect size ratio. The obtained results can be used to scale up the fatigue behavior from sub-sized specimens to the standard ones, as well as to predict the fatigue behavior of the bulk nanostructured metals and alloys produced by HPT.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Microwave ablation at 915 MHz vs 2.45 GHz: A theoretical and experimental investigation

Purpose: The relationship between microwave ablation system operating frequency and ablation performance is not currently well understood. The objective of this study was to comparatively assess the differences in microwave ablation at 915 MHz and 2.45 GHz. Methods: Analytical expressions for electromagnetic radiation from point sources were used to compare power deposition at the two frequencies of interest. A 3D electromagnetic-thermal bioheat transfer solver was implemented with the finite element method to characterize power deposition and thermal ablation with asymmetrical insulated dipole antennas (single-antenna and dual-antenna synchronous arrays). Simulation results were validated against experiments in ex vivo tissue. Results: Theoretical, computational, and experimental results indicated greater power deposition and larger diameter ablation zones when using a single insulated microwave antenna at 2.45 GHz; experimentally, 32 ± 4.1 mm and 36.3 ± 1.0 mm for 5 and 10 min, respectively, at 2.45 GHz, compared to 24 ± 1.7 mm and 29.5 ± 0.6 mm at 915 MHz, with 30 W forward power at the antenna input port. In experiments, faster heating was observed at locations 5 mm (0.91 vs 0.49 °C/s) and 10 mm (0.28 vs 0.15 °C/s) from the antenna operating at 2.45 GHz. Larger ablation zones were observed with dual-antenna arrays at 2.45 GHz; however, the differences were less pronounced than for single antennas. Conclusions: Single- and dual-antenna arrays systems operating at 2.45 GHz yield larger ablation zone due to greater power deposition in proximity to the antenna, as well as greater role of thermal conduction.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A mathematical model of the passage of an asteroid-comet body through the Earth’s atmosphere

In the paper, a mathematical model and a numerical algorithm are proposed for modeling the complex of phenomena which accompany the passage of a friable asteroid-comet body through the Earth’s atmosphere: the material ablation, the dissociation of molecules, and the radiation. The proposed model is constructed on the basis of the Navier-Stokes equations for viscous heat-conducting gas with an additional equation for the motion and propagation of a friable lumpy-dust material in air. The energy equation is modified for the relation between two its kinds: the usual energy of the translation of molecules (which defines the temperature and pressure) and the combined energy of their rotation, oscillation, electronic excitation, dissociation, and radiation. For the mathematical model of atmosphere, the distribution of density, pressure, and temperature in height is taken as for the standard atmosphere. An asteroid-comet body is taken initially as a round body consisting of a friable lumpy-dust material with corresponding density and significant viscosity which far exceed those for the atmosphere gas. A numerical algorithm is proposed for solving the initial-boundary problem for the extended system of Navier-Stokes equations. The algorithm is the combination of the semi-Lagrangian approximation for Lagrange transport derivatives and the conforming finite element method for other terms. The implementation of these approaches is illustrated by a numerical example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Heat transfer in porous medium embedded with vertical plate: Non-equilibrium approach - Part A

Heat transfer in a porous medium embedded with vertical flat plate is investigated by using thermal non-equilibrium model. Darcy model is employed to simulate the flow inside porous medium. It is assumed that the heat transfer takes place by natural convection and radiation. The vertical plate is maintained at isothermal temperature. The governing partial differential equations are converted into non-dimensional form and solved numerically using finite element method. Results are presented in terms of isotherms and streamlines for various parameters such as heat transfer coefficient parameter, thermal conductivity ratio, and radiation parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

A THREE-DIMENSIONAL NUMERICAL SOLUTION FOR THE SHAPE OF A ROTATIONALLY DISTORTED POLYTROPE OF INDEX UNITY

We present a new three-dimensional numerical method for calculating the non-spherical shape and internal structure of a model of a rapidly rotating gaseous body with a polytropic index of unity. The calculation is based on a finite-element method and accounts for the full effects of rotation. After validating the numerical approach against the asymptotic solution of Chandrasekhar that is valid only for a slowly rotating gaseous body, we apply it to models of Jupiter and a rapidly rotating, highly flattened star ({alpha} Eridani). In the case of Jupiter, the two-dimensional distributions of density and pressure are determined via a hybrid inverse approach by adjusting an a priori unknown coefficient in the equation of state until the model shape matches the observed shape of Jupiter. After obtaining the two-dimensional distribution of density, we then compute the zonal gravity coefficients and the total mass from the non-spherical model that takes full account of rotation-induced shape change. Our non-spherical model with a polytropic index of unity is able to produce the known mass of Jupiter with about 4% accuracy and the zonal gravitational coefficient J {sub 2} of Jupiter with better than 2% accuracy, a reasonable result considering that there is only one parameter in the model. For {alpha} Eridani, we calculate its rotationally distorted shape and internal structure based on the observationally deduced rotation rate and size of the star by using a similar hybrid inverse approach. Our model of the star closely approximates the observed flattening.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Scalability Tests of a Finite Element Code on Hundreds of Thousands Cores and Heterogeneous Architecture

A multi-scale finite element method code, msFEM, is tested on Jaguar and Nebulae, two petaflops computers that were listed as #1 and #2 on the Top500 list of June 2010 at the time of the tests. The flat MPI version of msFEM is scaled from 20K up to 200K CPU cores on Jaguar, delivering 70% parallel efficiency at the 200K cores with a finite element model of eight millions of degrees of freedom. GPU versions, in both double precision and mixed precision coded through MPI+OpenMP+CUDA hybrid programming, 900 GPU nodes on Jaguar and 1500 GPU nodes on Nebulae, achieving remarkable 90 + % parallel efficiency on the systems. The mixed-precision GPU version delivers further 1.5 times of speedup over the fully double precision version with no significant implementational cost. The large-scale tests support that the msFEM runs efficiently on petaflops computers and is highly potential for domain applications at extreme-scale.

finite element method↗
DOE OSTI2013

Hamiltonian discontinuous Galerkin FEM for linear, rotating incompressible Euler equations: Inertial waves

A discontinuous Galerkin finite element method (DGFEM) has been developed and tested for the linear, three-dimensional, rotating incompressible Euler equations. These equations admit complicated wave solutions, which poses numerical challenges. These challenges concern: (i) discretisation of a divergence-free velocity field; (ii) discretisation of geostrophic boundary conditions combined with no-normal flow at solid walls; (iii) discretisation of the conserved, Hamiltonian dynamics of the inertial-waves; and, (iv) large-scale computational demands owing to the three-dimensional nature of inertial-wave dynamics and possibly its narrow zones of chaotic attraction. These issues have been resolved, for example: (i) by employing Dirac’s method of constrained Hamiltonian dynamics to our DGFEM for linear, compressible flows, thus enforcing the incompressibility constraints; (ii) by enforcing no-normal flow at solid walls in a weak form and geostrophic tangential flow along the wall; and, (iii) by applying a symplectic time discretisation. We compared our simulations with exact solutions of three-dimensional incompressible flows, in (non) rotating periodic and partly periodic cuboids (Poincaré waves). Additional verifications concerned semi-analytical eigenmode solutions in rotating cuboids with solid walls. Finally, a simulation in a tilted rotating tank, yielding more complicated wave dynamics, demonstrates the potential of our new method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

WE-G-17A-09: Novel Magnetic Shielding Design for Inline and Perpendicular Integrated 6 MV Linac and 1.0 T MRI Systems

Purpose: The influence of fringe magnetic fields delivered by magnetic resonance imaging (MRI) on the beam generation and transportation in Linac is still a major challenge for the integration of linear accelerator and MRI (Linac-MRI). In this study, we investigated an optimal magnetic shielding design for Linac-MRI and further characterized the beam trajectory in electron gun. Methods: Both inline and perpendicular configurations were analyzed in this study. The configurations, comprising a Linac-MRI with a 100cm SAD and an open 1.0 T superconductive magnet, were simulated by the 3D finite element method (FEM). The steel shielding around the Linac was included in the 3D model, the thickness of which was varied from 1mm to 20mm, and magnetic field maps were acquired with and without additional shielding. The treatment beam trajectory in electron gun was evaluated using OPERA 3d SCALA with and without shielding cases. Results: When Linac was not shielded, the uniformity of diameter sphere volume (DSV) (30cm) was about 5 parts per million (ppm) and the fringe magnetic fields in electron gun were more than 0.3 T. With shielding, the magnetic fields in electron gun were reduced to less than 0.01 T. For the inline configuration, the radial magnetic fields in the Linac were about 0.02T. A cylinder steel shield used (5mm thick) altered the uniformity of DSV to 1000 ppm. For the perpendicular configuration, the Linac transverse magnetic fields were more than 0.3T, which altered the beam trajectory significantly. A 8mm-thick cylinder steel shield surrounding the Linac was used to compensate the output losses of Linac, which shifted the magnetic fields' uniformity of DSV to 400 ppm. Conclusion: For both configurations, the Linac shielding was used to ensure normal operation of the Linac. The effect of magnetic fields on the uniformity of DSV could be modulated by the shimming technique of the MRI magnet. NIH/NIGMS grant U54 GM104944, Lincy Endowed Assistant Professorship.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Screening sensitivity analysis of a PWR fuel assembly FEM structural model

Within the present project, computational methods to describe the static mechanical response of PWR fuel assemblies (FA) in the reactor core are developed. For this purpose a structural model of a PWR fuel assembly has been implemented with the Finite Element Method (FEM) using the ANSYS APDL Mechanical code. The purpose of the model is to investigate the long-term creep deformation of the fuel assemblies in the reactor core, commonly known as fuel assembly bow. FA bow is considered to be a complex process with a large number of influencing mechanisms and many only roughly known boundary conditions. Uncertainty and sensitivity analyses are a common way to investigate the predictive power and reliability of such complex models. Hence, as a first approach to understand the importance of the many model parameters, a qualitative sensitivity analysis of the model was conducted using the screening method proposed by Morris. In the first part of the paper, the FA structural model and its principal parameters are presented. The second part presents first the different simulation experiments used for the screening analyses. Based on the model description, relevant input parameters are then identified and a typical range of values is determined. The sensitivity analysis is performed in two subsequent steps. First, the model is screened for influential stiffness parameters with the objective of stiffening the FA structure. Second, the influence of the uncertain boundary conditions and material models on the FA creep deformation is tested. For the considered FA design it is concluded that the influence of the involved uncertainties is comparable to the influence of the parameters which can increase the stiffness. An integral approach accounting for these uncertainties is hence necessary to predict FA bow reliably in order to take additional measures for a further reduction of FA bow amplitudes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Preparation of room temperature terahertz detector with lithium tantalate crystal and thin film

Research on room temperature terahertz (THz) detector is essential for promoting the application of THz science and technology. Both lithium tantalate crystal (LiTaO{sub 3}) and lithium tantalate thin film were used to fabricate the THz detector in this paper. Polishing process were used to reduce the thickness of LiTaO{sub 3} crystal slice by chemical mechanical polishing techniques and an improved sol-gel process was used to obtain high concentration LiTaO{sub 3} precursor solution to fabricate LiTaO{sub 3} thin film. Three dimension models of two THz detectors were set up and the temperature increasing map of two devices were simulated using finite element method. The lowest noise equivalent power value for terahertz detector using pyroelectric material reaches 6.8 × 10{sup −9} W at 30 Hz operating frequency, which is suitable for THz imaging application.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Newmark local time stepping on high-performance computing architectures

In multi-scale complex media, finite element meshes often require areas of local refinement, creating small elements that can dramatically reduce the global time-step for wave-propagation problems due to the CFL condition. Local time stepping (LTS) algorithms allow an explicit time-stepping scheme to adapt the time-step to the element size, allowing near-optimal time-steps everywhere in the mesh. We develop an efficient multilevel LTS-Newmark scheme and implement it in a widely used continuous finite element seismic wave-propagation package. In particular, we extend the standard LTS formulation with adaptations to continuous finite element methods that can be implemented very efficiently with very strong element-size contrasts (more than 100x). Capable of running on large CPU and GPU clusters, we present both synthetic validation examples and large scale, realistic application examples to demonstrate the performance and applicability of the method and implementation on thousands of CPU cores and hundreds of GPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Multiscale Capability in Rattlesnake using Contiguous Discontinuous Discretization of Self-Adjoint Angular Flux Equation

In this report, we present a new upwinding scheme for the multiscale capability in Rattlesnake, the MOOSE based radiation transport application. Comparing with the initial implementation of multiscale utilizing Lagrange multipliers to impose strong continuity of angular flux on interface of in-between subdomains, this scheme does not require the particular domain partitioning. This upwinding scheme introduces discontinuity of angular flux and resembles the classic upwinding technique developed for solving first order transport equation using discontinuous finite element method (DFEM) on the subdomain interfaces. Because this scheme restores the causality of radiation streaming on the interfaces, significant accuracy improvement can be observed with moderate increase of the degrees of freedom comparing with the continuous method over the entire solution domain. Hybrid SN-PN is implemented and tested with this upwinding scheme. Numerical results show that the angular smoothing required by Lagrange multiplier method is not necessary for the upwinding scheme.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Strength of concrete structures under dynamic loading

The use of elastic supports is one the efficient methods of decreasing the dynamic loading. The paper describes the influence of elastic supports on the stress-strain state of steel concrete structures exposed to one-time dynamic loading resulting in failure. Oblique bending beams on elastic supports and their elastic, elastoplastic, and elastoplastic consolidation behavior are considered in this paper. For numerical calculations the developed computer program is used based on the finite element method. Research findings prove high efficiency of elastic supports under dynamic loading conditions. The most effective behavior of elastic supports is demonstrated at the elastoplastic stage. A good agreement is observed between the theoretical and experimental results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Sensitivity analysis of static resistance of slender beam under bending

The paper deals with statical and sensitivity analyses of resistance of simply supported I-beams under bending. The resistance was solved by geometrically nonlinear finite element method in the programme Ansys. The beams are modelled with initial geometrical imperfections following the first eigenmode of buckling. Imperfections were, together with geometrical characteristics of cross section, and material characteristics of steel, considered as random quantities. The method Latin Hypercube Sampling was applied to evaluate statistical and sensitivity resistance analyses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Stress field and interaction forces between dislocations and precipitate distributions

A computational method is developed for calculation of the stress field and interaction forces between dislocations and precipitates of arbitrary shape and distribution. The internal stress generated by precipitates due to coherency strain is implemented within the discrete dislocation dynamics (DDD) framework. The s-version finite element method (s-FEM), which models a precipitate of arbitrary shape using a local mesh is used to calculate coherency stress fields. The method facilitates meshing precipitate volumes of arbitrary geometry, and multiple local meshes can be superimposed at various positions of the global mesh. Accuracy and convergence conditions of the method are established. For a single precipiate, the method is shown to be 3.5 times faster than a standard FEM calculation for the same accuracy, and the gain in speed increases with the number of precipitates. The CRSS of spherical precipitates is found to be higher than disk-shaped precipitates of the same volume fraction.

42 ENGINEERING↗
DOE OSTI2014

Application of stochastic Galerkin FEM to the complete electrode model of electrical impedance tomography

The aim of electrical impedance tomography is to determine the internal conductivity distribution of some physical body from boundary measurements of current and voltage. The most accurate forward model for impedance tomography is the complete electrode model, which consists of the conductivity equation coupled with boundary conditions that take into account the electrode shapes and the contact resistances at the corresponding interfaces. If the reconstruction task of impedance tomography is recast as a Bayesian inference problem, it is essential to be able to solve the complete electrode model forward problem with the conductivity and the contact resistances treated as a random field and random variables, respectively. In this work, we apply a stochastic Galerkin finite element method to the ensuing elliptic stochastic boundary value problem and compare the results with Monte Carlo simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A New Configuration of Paralleled Modular ANPC Multilevel Converter Controlled by an Improved Modulation Method for 1 MHz, 1 MW EV Charger

In this work, a new configuration of the modular multilevel converter (MLC) based on the parallel connection of three-level active-neutral-point-clamped (3L-ANPC) cells as well as its improved modulation method is proposed for 1 MHz, 1 MW electric vehicle (EV) megacharger. In the proposed paralleled modular ANPC-MLC, only six high-frequency silicon carbide (SiC) power switches operating at 333 kHz are required to generate 1 MHz switching frequency spectrum. Moreover, the operating voltage of all power devices is halved, the magnitude of the first switching frequency harmonic cluster is decreased by the factor of five, and the load current is equally distributed between the 3L-ANPC legs by employing the proposed improved modulation method. Hence, the modularity, efficiency, and power density of the proposed converter are notably increased, whereas the value of passive components and the overall switching loss are remarkably decreased. In addition, an optimized design of the one 3L-ANPC cell of the proposed paralleled modular ANPC-MLC for 1 MHz, 1MW EV megacharger using Ansys SIwave, Icepak, and Q3D finite element method platforms is presented and analyzed in detail. The provided experimental results of the down-scaled setup verify the feasibility and viability of the proposed configuration as well as its improved switching pattern.

42 ENGINEERING↗
DOE OSTI2018

Effect of Tempering and Rolling on Fatigue Crack Growth Behavior of Modified 9Cr-1Mo Steel

In the present work, fatigue crack growth (FCG) behavior of modified 9Cr-1Mo steel has been investigated near fatigue threshold and Paris regime at room temperature (RT) and 650 °C. The modified 9Cr-1Mo steel was subjected to normalizing at 1050 °C followed by tempering (650-750 °C) and warm rolling at 550 °C followed by tempering in order to improve its precipitation state and refine microstructure. The microstructural parameters such as grain size, block size, lath width and precipitate size were examined and correlated with FCG and tensile properties of investigated samples. The results showed that decreasing tempering temperature significantly reduces the fatigue threshold at RT while increases the yield strength at elevated temperature (i.e., 650 °C). Increase in fatigue threshold is found to increase with the increase in block size, whereas high density of fine precipitates contributed the high temperature strength. The FCG rate in Paris regime was less affected by heat treatment and rolling process at RT. However, rolled sample tested at 650 °C shows a decrease in FCG rate as compared to heat-treated samples. In terms of higher FCG resistance and yield strength, warm rolling followed by tempering at low temperature (i.e., 700 °C) can be considered optimum. For the numerical simulation of FCG behavior, extended finite element method (XFEM) along with Paris law has been implemented in ABAQUS using python script. The simulated FCG behavior is found in good agreement with the experimental results.

36 MATERIALS SCIENCE↗