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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Illusion optics: Optically transforming the nature and the location of electromagnetic emissions

Complex electromagnetic structures can be designed by using the powerful concept of transformation electromagnetics. In this study, we define a spatial coordinate transformation that shows the possibility of designing a device capable of producing an illusion on an antenna radiation pattern. Indeed, by compressing the space containing a radiating element, we show that it is able to change the radiation pattern and to make the radiation location appear outside the latter space. Both continuous and discretized models with calculated electromagnetic parameter values are presented. A reduction of the electromagnetic material parameters is also proposed for a possible physical fabrication of the device with achievable values of permittivity and permeability that can be obtained from existing well-known metamaterials. Following that, the design of the proposed antenna using a layered metamaterial is presented. Full wave numerical simulations using Finite Element Method are performed to demonstrate the performances of such a device.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1980

Finite element models and feedback control of flexible aerospace structures

Large flexible aerospace structures, such as the solar power satellite, are distributed parameter systems with very complex continuum descriptions. This paper investigates the use of finite element methods to produce reduced-order models and finite dimensional feedback controllers for these structures. The main results give conditions under which stable control of the finite element model will produce stable control of the actual structure.

Balas, M. J.↗
DOE OSTI2014

Optical and mechanical mode tuning in an optomechanical crystal with light-induced thermal effects

We report on the modification of the optical and mechanical properties of a silicon 1D optomechanical crystal cavity due to thermo-optic effects in a high phonon/photon population regime. The cavity heats up due to light absorption in a way that shifts the optical modes towards longer wavelengths and the mechanical modes to lower frequencies. By combining the experimental optical results with finite-difference time-domain simulations, we establish a direct relation between the observed wavelength drift and the actual effective temperature increase of the cavity. By assuming that the Young's modulus decreases accordingly to the temperature increase, we find a good agreement between the mechanical mode drift predicted using a finite element method and the experimental one.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2004

Nonlinear Analysis of Bonded Composite Single-LAP Joints

This study presents a semi-analytical solution method to analyze the geometrically nonlinear response of bonded composite single-lap joints with tapered adherend edges under uniaxial tension. The solution method provides the transverse shear and normal stresses in the adhesive and in-plane stress resultants and bending moments in the adherends. The method utilizes the principle of virtual work in conjunction with von Karman s nonlinear plate theory to model the adherends and the shear lag model to represent the kinematics of the thin adhesive layer between the adherends. Furthermore, the method accounts for the bilinear elastic material behavior of the adhesive while maintaining a linear stress-strain relationship in the adherends. In order to account for the stiffness changes due to thickness variation of the adherends along the tapered edges, their in-plane and bending stiffness matrices are varied as a function of thickness along the tapered region. The combination of these complexities results in a system of nonlinear governing equilibrium equations. This approach represents a computationally efficient alternative to finite element method. Comparisons are made with corresponding results obtained from finite-element analysis. The results confirm the validity of the solution method. The numerical results present the effects of taper angle, adherend overlap length, and the bilinear adhesive material on the stress fields in the adherends, as well as the adhesive, of a single-lap joint

Oterkus, E.↗
DOE OSTI2016

Radiation and viscous dissipation effect on square porous annulus

The present study is carried out to investigate the effect of radiation and viscous dissipation in a square porous annulus subjected to outside hot T{sub h} and inside cold T{sub c} temperature. The square annulus has a hollow section of dimension D×D at the interior of annulus. The flow is assumed to obey Darcy law. The governing equations are non-dimensionalised and solved with the help of finite element method. Results are discussed with respect to viscous dissipation parameter, radiation parameter and size of the hollow section of annulus.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

FEM design and simulation of a short, 10 MV, S-band Linac with Monte Carlo dose simulations

Purpose: Current commercial 10 MV Linac waveguides are 1.5 m. The authors’ current 6 MV linear accelerator–magnetic resonance imager (Linac–MR) system fits in typical radiotherapy vaults. To allow 10 MV treatments with the Linac–MR and still fit within typical vaults, the authors design a 10 MV Linac with an accelerator waveguide of the same length (27.5 cm) as current 6 MV Linacs. Methods: The first design stage is to design a cavity such that a specific experimental measurement for breakdown is applicable to the cavity. This is accomplished through the use of finite element method (FEM) simulations to match published shunt impedance, Q factor, and ratio of peak to mean-axial electric field strength from an electric breakdown study. A full waveguide is then designed and tuned in FEM simulations based on this cavity design. Electron trajectories are computed through the resulting radio frequency fields, and the waveguide geometry is modified by shifting the first coupling cavity in order to optimize the electron beam properties until the energy spread and mean energy closely match values published for an emulated 10 MV Linac. Finally, Monte Carlo dose simulations are used to compare the resulting photon beam depth dose profile and penumbra with that produced by the emulated 10 MV Linac. Results: The shunt impedance, Q factor, and ratio of peak to mean-axial electric field strength are all matched to within 0.1%. A first coupling cavity shift of 1.45 mm produces an energy spectrum width of 0.347 MeV, very close to the published value for the emulated 10 MV of 0.315 MeV, and a mean energy of 10.53 MeV, nearly identical to the published 10.5 MeV for the emulated 10 MV Linac. The depth dose profile produced by their new Linac is within 1% of that produced by the emulated 10 MV spectrum for all depths greater than 1.5 cm. The penumbra produced is 11% narrower, as measured from 80% to 20% of the central axis dose. Conclusions: The authors have successfully designed and simulated an S-band waveguide of length of 27.5 cm capable of producing a 10 MV photon beam. This waveguide operates well within the breakdown threshold determined for the cavity geometry used. The designed Linac produces depth dose profiles similar to those of the emulated 10 MV Linac (waveguide-length of 1.5 m) but yields a narrower penumbra.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Risk-Averse Control of Fractional Diffusion with Uncertain Exponent

In this paper, we introduce and analyze a new class of optimal control problems constrained by elliptic equations with uncertain fractional exponents. We utilize risk measures to formulate the resulting optimization problem. We develop a functional analytic framework, study the existence of solution and rigorously derive the first-order optimality conditions. Additionally, we employ a sample-based approximation for the uncertain exponent and the finite element method to discretize in space. Further, we prove the rate of convergence for the optimal risk neutral controls when using quadrature approximation for the uncertain exponent and conclude with illustrative examples.

finite element method↗
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 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↗