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,135 records · Page 63

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↗

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↗

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↗

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↗

Experimental and FEM Study of Residual Stresses During Ultra-Precision Turning of Aluminum 2024-T3

Aluminum 2024-T3 is commonly used in ultra-precision turning. In this study, the residual stresses of aluminum 2024-T3 were investigated by the finite element method (FEM) and a series of cutting experiments. The effects of depth of cut, cutting speed, rake angle and tool edge radius on the residual stresses were discussed by FEM. Moreover, an orthogonal simulation experiment was also conducted. The order of the strong-to-weak sequence of factors on residual stresses in cutting direction is the rake angle, tool edge radius, depth of cut and cutting speed. The cutting experiments were conducted and the surface residual stresses in the cutting direction were measured by X-ray diffraction. The experimental results showed that increasing cutting speed and decreasing tool edge radius can reduce the surface residual stresses.

36 MATERIALS SCIENCE↗

A simple design of highly birefringent and nonlinear photonic crystal fiber with ultra-flattened dispersion

We present a compressed hexagonal photonic crystal fiber (CH-PCF) based on a simple structure with ultra-flattened dispersion, high birefringence and nonlinearity. The propagation characteristics are studied accurately by utilizing full vector finite-element method (FEM). Numerical results indicate that the birefringence of 1.59 × 10{sup −2} was obtained with high nonlinear coefficient of 42.58 W{sup −1} km{sup −1}(X) and 52.80 W{sup −1} km{sup −1}(Y) at an effective wavelength of 1.55 μm. Furthermore, CH-PCF presents an ultra-flattened dispersion with variation of ± 3.4 ps/nm/km(X) and ± 1.6 ps/nm/km(Y) over 1.3 μm to 1.8 μm (0.5 μm bandwidth), respectively. In addition, this circular-hole structure can be fabricated conveniently at most easy compared with complex designs, which is an advantage for practical applications such as dispersion compensation and nonlinear optics application.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetic field calculation for a 10 MeV positron emission tomography cyclotron

The magnetic field calculation and correction for a 10 MeV positron emission tomography cyclotron is presented. 3D TOSCA analysis results are compared with the measured data, and the calculation error is used to calibrate the B-H curve to obtain a very precise finite element method estimator, which is used to predict the correction of the magnet pole for achieving the isochronous field. The isochronous field error is approximated with the effects of a set of standard patches. On the assumption that the effect of each small patch is proportional to its surface, the correction of the magnet pole is found by solving a system of equations using the least square scheme. The magnet shimming is performed and the measured magnetic field is found in good agreement with the prediction, with an error less than 2 G.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗