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 1,099 records · Page 61

Stress Relaxation Related to Spontaneous Thin Film Buckling: Correlation between Finite Element Calculations and Micro Diffraction Analysis

Compressive residual stresses generated during thin film deposition may lead to undesirable film damage, such as delamination, buckling, and flaking, ultimately leading to the failure of the device employing the film. Understanding the residual stress generation and role in these damage mechanisms is necessary to preserve thin film integrity and optimize its functional properties. Thin shell theory has been used for decades to predict buckling but the results have not yet been correlated with experimental data since the techniques used to measure stress in metallic films were not able to do so at the required micron scale until recently. Micro scanning X-ray diffraction now enables the direct mapping of the local stress of metallic films. In this paper, finite element method based on thin shell theory and synchrotron X-ray micro diffraction have been used to determine stress maps of thin film buckling patterns. Calculations of the stress distribution in the metallic films have been performed taking into account the buckling geometry determined from optical measurements. Stress distributions over gold blisters and tungsten wrinkles obtained with the two techniques are in fair agreement and allow for the accurate determination of the stress relaxation profile from the bottom to the top of the buckling, validating the thin shell theory model.

36 MATERIALS SCIENCE↗

Determination of the mechanical properties of blisters by nanoindentation

In light-water reactors, a zirconium oxide layer will be formed at the outer surface of the fuel cladding as a consequence of the chemical reaction which takes place between zirconium and the cooling water. This zirconium oxide layer is a good thermal insulator, but at some point, it may spall-out. As a consequence, a cold spot will be formed at the outer surface of the cladding. In this situation, the hydrogen will quickly migrate to this spot to form a highly-enriched hydride zone, known as blister. In this work, an experimental procedure has been devised to produce blisters with a similar shape to the ones reported in the literature for irradiated cladding with severe hydriding. Nano-indentation tests were performed on blisters prepared in laboratory. A numerical simulation by the finite element method was developed to fit the load versus penetration depth data from the aforementioned nano-indentation tests. By means of an inverse process which combines experimental data and numerical simulations, the constitutive equation of the blister was determined. The process has been used to calculate the plastic stress-strain curve of the blister under study at room temperature. The results show an almost perfect numerical fitting of the experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Gold grating assisted SPR based D-shaped single mode fiber for detection of liquid refractive index

In this paper, a surface plasmon resonance based D-shaped single mode fiber sensor with a gold grating over the polished fiber surface, is analyzed by finite element method for refractive index sensing. Benefitted from the localized surface plasmon resonance effect of gold grating, the proposed sensor achieves a reasonably good wavelength sensitivity of 7590 nm/RIU for the refractive index range varying from 1.34 to 1.38 and the corresponding resolution is 1.31 × 10{sup −5} RIU. Furthermore, the maximum amplitude sensitivity of the sensor is calculated to be 676.08 RIU{sup −1}. The good sensing performance, low cost and simple structure makes this sensor a potential candidate for application in the areas of food processing, environmental monitoring and biochemical analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields

An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.

Fahrenthold, Eric P.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗