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

Evolution of linear perturbations in Lemaítre-Tolman-Bondi void models

We study the evolution of linear perturbations in a Lemaítre-Tolman-Bondi (LTB) void model with realistic cosmological initial conditions. Linear perturbation theory in LTB models is substantially more complicated than in standard Friedmann universes as the inhomogeneous background causes gauge-invariant perturbations couple at first order. As shown by Clarkson et al. (2009) [1], the evolution is constrained by a system of linear partial differential equations which need to be integrated numerically. We present a new numerical scheme using finite element methods to solve this equation system and generate scalar initial conditions based on Gaussian random fields with an underlying power spectrum for the Bardeen potential. After spherical harmonic decomposition, the initial fluctuations are propagated in time and estimates of angular power spectra of each gauge invariant variable are computed as functions of redshift. This allows to analyse the coupling strength in a statistical way. We find significant couplings up to 25% for large and deep voids of Gpc scale as required to fit the distance redshift relations of SNe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Surface-plasmon enhanced photodetection at communication band based on hot electrons

Surface plasmons can squeeze light into a deep-subwavelength space and generate abundant hot electrons in the nearby metallic regions, enabling a new paradigm of photoconversion by the way of hot electron collection. Unlike the visible spectral range concerned in previous literatures, we focus on the communication band and design the infrared hot-electron photodetectors with plasmonic metal-insulator-metal configuration by using full-wave finite-element method. Titanium dioxide-silver Schottky interface is employed to boost the low-energy infrared photodetection. The photodetection sensitivity is strongly improved by enhancing the plasmonic excitation from a rationally engineered metallic grating, which enables a strong unidirectional photocurrent. With a five-step electrical simulation, the optimized device exhibits an unbiased responsivity of ∼0.1 mA/W and an ultra-narrow response band (FWHM = 4.66 meV), which promises to be a candidate as the compact photodetector operating in communication band.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Preliminary analysis of FCM fuel thermal-mechanical performance under normal and transient conditions

Fully ceramic microencapsulated (FCM) fuel is a candidate fuel form to replace commercial UO{sub 2} fuel rod in LWR to enhance fuel reliability under accident conditions. The FCM pellet consists of TRISO particles embedded inside the dense SiC matrix. This fuel design is more complex than conventional LWR fuel, which needs new fuel performance models or analysis tools to assist fuel design and concept evaluation. In this paper an attempt was made to develop a finite element method to analyze FCM fuel pellet thermal-mechanical performance. Three idealized cases of normal operation, power ramp, and LOCA were analyzed. The results showed that the temperature distribution and stress distribution of FCM pellet exhibited a strong non-uniformity. The irradiation strain and creep of buffer layer and PyC (pyrolytic carbon) layers had a strong effect on the mechanical performance of FCM pellet. The hoop stresses of SiC layers and SiC matrix increased significantly during the LOCA and varied with radial positions. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Modified D-shaped SPR PCF polarization filter at telecommunication wavelengths

A compact surface plasmon photonic crystal fiber (PCF) polarization filter is proposed and analyzed at the two communication windows of 1.3 and 1.52 µm. The y-polarized mode can be filtered out at wavelength of 1.3 µm while the x-polarized mode can be excluded at λ = 1.52 µm using the same design. The effects of the geometrical parameters of the PCF on the polarization filter performance are studied by full vectorial finite element method. The numerical simulations reveal that the suggested polarization filter with a compact device length of 23 µm can achieve good crosstalk (CT) of − 46.4 and − 41.2 dB at the communication wavelengths of 1.3 and 1.52 µm, respectively. Further, bandwidths of 45 nm and 110 nm can be obtained around λ = 1.3 µm and 1.52 µm, respectively where the CT is better than − 20 dB. Therefore, the reported design has advantages in terms of low crosstalk, large bandwidth, short device length and multi-functionality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Multiwall carbon nanotube microcavity arrays

Periodic highly dense multi-wall carbon nanotube (MWCNT) arrays can act as photonic materials exhibiting band gaps in the visible regime and beyond terahertz range. MWCNT arrays in square arrangement for nanoscale lattice constants can be configured as a microcavity with predictable resonance frequencies. Here, computational analyses of compact square microcavities (≈0.8 × 0.8 μm{sup 2}) in MWCNT arrays were demonstrated to obtain enhanced quality factors (≈170–180) and narrow-band resonance peaks. Cavity resonances were rationally designed and optimized (nanotube geometry and cavity size) with finite element method. Series (1 × 2 and 1 × 3) and parallel (2 × 1 and 3 × 1) combinations of microcavities were modeled and resonance modes were analyzed. Higher order MWCNT microcavities showed enhanced resonance modes, which were red shifted with increasing Q-factors. Parallel microcavity geometries were also optimized to obtain narrow-band tunable filtering in low-loss communication windows (810, 1336, and 1558 nm). Compact series and parallel MWCNT microcavity arrays may have applications in optical filters and miniaturized optical communication devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A finite-element visualization of quantum reactive scattering. II. Nonadiabaticity on coupled potential energy surfaces

This is the second in a series of papers detailing a MATLAB based implementation of the finite element method applied to collinear triatomic reactions. Here, we extend our previous work to reactions on coupled potential energy surfaces. The divergence of the probability current density field associated with the two electronically adiabatic states allows us to visualize in a novel way where and how nonadiabaticity occurs. A two-dimensional investigation gives additional insight into nonadiabaticity beyond standard one-dimensional models. We study the F({sup 2}P) + HCl and F({sup 2}P) + H{sub 2} reactions as model applications. Our publicly available code (http://www2.chem.umd.edu/groups/alexander/FEM) is general and easy to use.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Diffusion Limit–Preserving Lumped DFEMs on AMR Meshes

Here, we present sweep-compatible, novel upwinding recipes for the bilinear discontinuous (BLD) finite element method (FEM) that allows lumped BLD to be used on adaptive mesh refinement (AMR) meshes for thick transport applications without adding additional degrees of freedom at hanging nodes that exist on refinement boundaries. We analyze the properties of the upwinding and lumping that are needed for BLD to get the thick diffusion limit on such meshes, present results demonstrating locking with the wrong recipe, and present results showing error convergence and robustness properties for two diffusive problems on a variety of AMR meshes.

42 ENGINEERING↗
DOE OSTI2012

ELECTRON TRANSPORT IN THE FAST SOLAR WIND

The electron velocity distribution function is studied in the extended solar corona above coronal holes (i.e., the inner part of the fast solar wind) from the highly collisional corona close to the Sun to the weakly collisional regions farther out. The electron kinetic equation is solved with a finite-element method in velocity space using a linearized Fokker-Planck collision operator. The ion density and temperature profiles are assumed to be known and the electric field and electron temperature are determined self-consistently. The results show quantitatively how much lower the electron heat flux and the thermal force are than predicted by high-collisionality theory. The sensitivity of the particle and heat fluxes to the assumed ion temperature profile and the applied boundary condition at the boundary far from the Sun is also studied.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Structural FEM analysis of the strut-to-fuselage joint of a two-seat composite aircraft

An analysis of a strut-to-fuselage joint is realized in order to evaluate the zones with a high probability of failure by means of a safety factor. The whole section is analyzed using the Finite Element Method (FEM) so as to estimate static resistance behavior, therefore it is necessary a numerical mock-up of the section, the mechanical properties of the Carbon-Epoxy (C-Ep) material, and to evaluate the applied loads. Results of the analysis show that the zones with higher probability of failure are found around the wing strut and the fuselage joint, with a safety factor lower than expected in comparison with the average safety factor used on aircrafts built mostly with metals.

42 ENGINEERING↗
DOE OSTI2022

Time domain probabilistic seismic risk analysis using ground motion prediction equations of Fourier amplitude spectra

Modeling of Fourier amplitude spectra (FAS) of seismic motions has gained much attention in engineering seismology. In the past few years, several ground motion prediction equations (GMPEs) and inter-frequency correlation structure of FAS have been established. Due to many preferable characteristics of FAS, probabilistic seismic hazard/risk analysis is rapidly changing from ergodic, spectrum acceleration Sa(T 0 )-based approach to non-ergodic, site-specific, FAS-based approach. This paper presents time domain intrusive framework for probabilistic seismic risk analysis using GMPE of FAS. Herein, methodology for time domain stochastic ground motion modeling based on GMPEs of FAS is presented in some detail. The simulated uncertain motions are modeled as a random process and represented by polynomial chaos Karhunen-Loève expansion. The random process excitations are further propagated into the uncertain structural system using Galerkin stochastic finite element method (SFEM). Probabilistic evolution of structural response is solved, and such solution is used to develop seismic risk for any damage state. The presented framework is illustrated through seismic risk analysis of a four-story building subjected to possible earthquakes from two strike slip faults. The influences of the epistemic uncertainties in source stress drop Δσ and site attenuation κ0 on seismic risk are investigated. The need for non-ergodic seismic risk analysis with source-specific and site specific characterizations is emphasized.

58 GEOSCIENCES↗
DOE OSTI2014

Simulation and analysis on ultrasonic testing for the cement grouting defects of the corrugated pipe

The defects exist in the cement grouting process of prestressed corrugated pipe may directly impair the bridge safety. In this paper, sound fields propagation in concrete structures with corrugated pipes and the influence of various different defects are simulated and analyzed using finite element method. The simulation results demonstrate a much complex propagation characteristic due to multiple reflection, refraction and scattering, where the scattering signals caused by metal are very strong, while the signals scattered by an air bubble are weaker. The influence of defect both in time and frequency domain are found through deconvolution treatment. In the time domain, the deconvolution signals correspond to larger defect display a larger head wave amplitude and shorter arrive time than those of smaller defects; in the frequency domain, larger defect also shows a stronger amplitude, lower center frequency and lower cutoff frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Conduction Heat Transfer through Solid in Porous Materials: A Comparative Study by Finite-Element Simulations and Effective Medium Approximations

The development of insulation materials with low effective thermal conductivity is essential for energy savings in various applications, including buildings, food services, pipe insulation, and refrigeration. Such materials can be developed by using micro- or nanoporous structures, as well as low-thermal conductivity gases and/or reduced pressure. A variety of effective medium approximation (EMA) models have been developed to study thermal transport through solid in porous structures. However, in many cases, the impacts of porosity on solid conduction are not well predicted by EMA models because of the assumptions made for simplification. Furthermore, the results can vary by a factor of up to 1.5, depending on the morphology of the pores. Hence, proper guidance is needed to choose the appropriate EMA model for a given morphology. This work presents a finite element method study using COMSOL Multiphysics software for various geometries, including hollow cubes, overlapping hollow spheres (normally stacked and tightly packed), and Voronoi structures, and compares the results to those obtained from various EMA models. Simulation results for the hollow cubes and Voronoi structures matched well with the Maxwell-Eucken and Russell model. The results for hollow spheres closely matched the Bauer model and the Glicksman model with fitting parameters. This work provides guidance on analyzing and designing insulation materials for energy savings in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Design considerations of highly efficient D-shaped plasmonic biosensor

Design considerations of photonic crystal fiber (PCF) surface plasmon biosensor with high sensitivity to monitor glucose concentration is reported and studied. Based on a well-known large mode area (LMA) single mode PCF, two different configurations have been studied to investigate the impact of the etching process on the biosensor sensitivity. Furthermore, the possibility of infiltrating one hole by a plasmonic material to increases the biosensor sensitivity is studied. A full vectorial finite element method (FVFEM) is used to carry out the analysis of the reported biosensor. In addition, the geometrical parameters of the suggested biosensors have been optimized to obtain the highest sensitivity. The suggested biosensor based on a D-shaped PCF with plasmonic rod achieves sensitivity as high as 13,600 nm/RIU with corresponding resolution of 7.35×10{sup −6} RIU.The analysis also reveals that the proposed biosensor has a linear performance which is needed practically. Therefore, the reported biosensor has advantages in terms of fabrication feasibility and high linear sensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Onset of transition from laminar to chaos in MHD mixed convection of a lid-driven trapezoidal cavity filled with Cu-water nanofluid

The present study investigates the thermal mixing scenarios of steady magneto-hydrodynamic (MHD) mixed convection in a two-dimensional lid-driven trapezoidal cavity filled with Cu-water nanofluid. The top wall of the cavity slides with a uniform velocity from left to right direction, while the other walls are fixed. The bottom wall is kept with a constant higher temperature than the top one. The governing mass, momentum and energy equations are expressed in non-dimensional forms and Galerkin finite element method has been employed to solve these equations. Special attention is paid on investigating the onset of transition from laminar to chaos at pure mixed convection case. Hence, the computations are carried out for a wide range of Reynolds numbers (Re = 0.1 − 400) and Grashof numbers (Gr = 10{sup −2} − 1.6 × 10{sup 5}) at unity Richardson number and fixed Hartmann number (Ha = 10). The variation of average Nusselt number of the bottom heated wall indicates the influence of governing parameters (Re and Gr) on heat transfer characteristics. The results are presented and explained through the visualisation of isotherms, streamlines and heatlines.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Dielectric relaxation analysis of Pb(Zr{sub 0.54},Ti{sub 0.46})O{sub 3} thin films: Electric field dependence

350 nm-thick Perovskite PbZr{sub 0.54}Ti{sub 0.46}O{sub 3} (PZT) thin films were deposited on Al{sub 2}O{sub 3} substrates by sputtering with and without an additional 10-nm-thick TiO{sub x} buffer layer. X-ray diffraction patterns showed that in presence of TiO{sub x} buffer layer, PZT film was highly oriented along the (111) direction film, whereas the unbuffered, counterpart was polycrystalline. A full wave electromagnetic analysis using a vector finite element method was performed to determine the tunability and the complex permittivity up to 67 GHz. A comparison between the electromagnetic analysis and Cole-Cole relaxation model was proposed. Through an original study of the relaxation time as a function of the electric field, values, such as 2 ps and 0.6 ps, were estimated for E{sub DC} = 0 kV/cm and 235 kV/cm, respectively, and in both cases (111)-PZT and polycrystalline-PZT. The distribution of relaxation times is found to be larger for (111)-PZT film, which is probably related to the film microstructure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Highly non-linear solid core photonic crystal fiber with one nano hole

The numerical study of newly designed solid core photonic crystal fiber (SCPCF) having three hexagonal air hole rings in cladding region and one small nano hole at the center are presented. By using full vectorial finite element method (FV-FEM), we analyses the optical properties like effective area, nonlinearity and confinement loss of the proposed PCF. Results show that the change in core diameter controls the effective area, nonlinearity and confinement loss. A low effective area (3.34 µm{sup 2}), high nonlinearity (36.34 W{sup −1}km{sup −1}) and low confinement loss (0.00106 dB/km) are achieved at the communication wavelength 1.55 µm for the SCPCF having core air hole diameter 0.10 µm, cladding air holes diameter 1.00 µm and pitch 2.50 µm. This type of PCF is very useful in non-linear applications such as supercontinuum generation, four wave mixing, second harmonic generation etc.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Ultrafast spectroscopy of super high frequency mechanical modes of doubly clamped beams

We use ultrafast pump-probe spectroscopy to study the mechanical vibrations in the time domain of doubly clamped silicon nitride beams. Beams with two different clamping conditions are investigated. Finite element method calculations are performed to analyse the mode spectra of both structures. By calculating the strain integral on the surface of the resonators, we are able to reproduce the effect of the detection mechanism and identify all the measured modes. We show that our spectroscopy technique combined with our modelling tools allow the investigation of several different modes in the super high frequency range (3-30 GHz) and above, bringing more information about the vibration modes of nanomechanical resonators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Temperature compensated transfer function for probability of detection (POD) in high temperature ultrasonic NDE using low temperature signals - 249

Ultrasonic piezoelectric based transducers offer the potential for in-coolant NDT measurements and imaging required at high temperature for the safe operation of advanced small modular reactors (SMRs). Current work estimates a significant 97% reduction in the probability of detection (POD) for an artificial flaw size of 2{lambda} due to degradation of the transducer performance at elevated temperature. The source of variability for POD in the current work is the temperature dependence of material matrix of piezoelectric material PZT-5A used in the sensor. The physics based model developed using finite element method quantifies reduction in the scattering amplitude for the side drilled hole (SDH) from 15 C to 195 C. Using model assisted probability of detection (MAPOD) approach, sensor temperature effect on the POD is demonstrated for SDHs in the low carbon steel specimen. Scattering amplitudes for the room temperature simulated data is validated with the experimental data. A temperature correction factor and transfer factors are reported to transform the simulated temperature effect from physics based model to the experimental data for estimation of high temperature POD using low temperature ultrasonic measurements. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗