SEARCH · Engineering Papers

Results for “Finite element methods”

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.

SourceResearch connectionFrom yearThrough year

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

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Application of the hypercube parallel processor to a large-scale moment method code

The applicability of a parallel computing architecture to the solution of a large-scale moment-method code is investigated. Specifically, the NEC (Numerical Electromagnetics Code) method-of-moments scattering program is implemented on a hypercube parallel processor. The accuracy and the increase in the speed of execution on this parallel architecture are demonstrated. The results show a very large reduction in execution time for large problems. The great potential of this parallel processor is shown for interactive solution of large NEC problems as well as other moment-method techniques such as the finite-element method.

Manshadi, Farzin↗

NASA NTRS2000

Magnetic Field Suppression of Flow in Semiconductor Melt

One of the most promising approaches for the reduction of convection during the crystal growth of conductive melts (semiconductor crystals) is the application of magnetic fields. Current technology allows the experimentation with very intense static fields (up to 80 KGauss) for which nearly convection free results are expected from simple scaling analysis in stabilized systems (vertical Bridgman method with axial magnetic field). However, controversial experimental results were obtained. The computational methods are, therefore, a fundamental tool in the understanding of the phenomena accounting during the solidification of semiconductor materials. Moreover, effects like the bending of the isomagnetic lines, different aspect ratios and misalignments between the direction of the gravity and magnetic field vectors can not be analyzed with analytical methods. The earliest numerical results showed controversial conclusions and are not able to explain the experimental results. Although the generated flows are extremely low, the computational task is a complicated because of the thin boundary layers. That is one of the reasons for the discrepancy in the results that numerical studies reported. Modeling of these magnetically damped crystal growth experiments requires advanced numerical methods. We used, for comparison, three different approaches to obtain the solution of the problem of thermal convection flows: (1) Spectral method in spectral superelement implementation, (2) Finite element method with regularization for boundary layers, (3) Multiquadric method, a novel method with global radial basis functions, that is proven to have exponential convergence. The results obtained by these three methods are presented for a wide region of Rayleigh and Hartman numbers. Comparison and discussion of accuracy, efficiency, reliability and agreement with experimental results will be presented as well.

Fedoseyev, A. I.↗

NASA NTRS2010

Probabilistic Structural Analysis Program

NASA/NESSUS 6.2c is a general-purpose, probabilistic analysis program that computes probability of failure and probabilistic sensitivity measures of engineered systems. Because NASA/NESSUS uses highly computationally efficient and accurate analysis techniques, probabilistic solutions can be obtained even for extremely large and complex models. Once the probabilistic response is quantified, the results can be used to support risk-informed decisions regarding reliability for safety-critical and one-of-a-kind systems, as well as for maintaining a level of quality while reducing manufacturing costs for larger-quantity products. NASA/NESSUS has been successfully applied to a diverse range of problems in aerospace, gas turbine engines, biomechanics, pipelines, defense, weaponry, and infrastructure. This program combines state-of-the-art probabilistic algorithms with general-purpose structural analysis and lifting methods to compute the probabilistic response and reliability of engineered structures. Uncertainties in load, material properties, geometry, boundary conditions, and initial conditions can be simulated. The structural analysis methods include non-linear finite-element methods, heat-transfer analysis, polymer/ceramic matrix composite analysis, monolithic (conventional metallic) materials life-prediction methodologies, boundary element methods, and user-written subroutines. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. NASA/NESSUS 6.2c is structured in a modular format with 15 elements.

Pai, Shantaram S.↗

DOE OSTI2016

Finite element analysis of true and pseudo surface acoustic waves in one-dimensional phononic crystals

In this paper, we report a theoretical investigation of surface acoustic waves propagating in one-dimensional phononic crystal. Using finite element method eigenfrequency and frequency response studies, we develop two model geometries suitable to distinguish true and pseudo (or leaky) surface acoustic waves and determine their propagation through finite size phononic crystals, respectively. The novelty of the first model comes from the application of a surface-like criterion and, additionally, functional damping domain. Exemplary calculated band diagrams show sorted branches of true and pseudo surface acoustic waves and their quantified surface confinement. The second model gives a complementary study of transmission, reflection, and surface-to-bulk losses of Rayleigh surface waves in the case of a phononic crystal with a finite number of periods. Here, we demonstrate that a non-zero transmission within non-radiative band gaps can be carried via leaky modes originating from the coupling of local resonances with propagating waves in the substrate. Finally, we show that the transmission, reflection, and surface-to-bulk losses can be effectively optimised by tuning the geometrical properties of a stripe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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.↗

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 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 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↗