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 991 records · Page 55

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Air-Coupled Nondestructive Evaluation Dissected

This work develops a two-dimensional theoretical model to simulate the behavior of a fully non-contact air-coupled nondestructive evaluation system for a thin isotropic plate. The model is divided into transmission, guided wave propagation and reception phase. The validation of the complete model was carried out by modeling the same system by means of finite element method using a Multiphysics software. In addition, the dependency of the generated Lamb waves on different transmitter’s parameters and incidence angle is thoroughly investigated. The results of the acoustic pressure excited by the transducer, the out-of-plane velocity amplitudes for the generated first antisymmetric Lamb wave mode, and the radiated pressure from the plate caused by the leaky Lamb wavefield were all compared between the two models and a reasonable degree of similarity was found.

42 ENGINEERING↗

DOE OSTI2016

Heat and mass transfer in porous cavity: Assisting flow

In this paper, investigation of heat and mass transfer in a porous cavity is carried out. The governing partial differential equations are non-dimensionalised and solved using finite element method. The left vertical surface of the cavity is maintained at constant temperature and concentration which are higher than the ambient temperature and concentration applied at right vertical surface. The top and bottom walls of the cavity are adiabatic. Heat transfer is assumed to take place by natural convection and radiation. The investigation is carried out for assisting flow when buoyancy and gravity force act in same direction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Adhesive and Cohesive Strength in FeB/Fe{sub 2}B Systems

In this work, FeB/Fe{sub 2}B systems were evaluated by the scratch test. The powder-pack boriding process was performed on the surface of AISI M2 steel. The mechanical parameters, such as yield stress and Young’s modulus of the boride layer, were obtained by the instrumented indentation technique. Residual stresses produced on the boride layer were estimated by using the x-ray diffraction (XRD) technique. The scratch test was performed in order to evaluate the cohesive/adhesive strength of the FeB/Fe{sub 2}B coating. In addition, a numerical evaluation of the scratch test on boride layers was performed by the finite element method. Maximum principal stresses were related to the failure mechanisms observed by the experimental scratch test. Shear stresses at the interfaces of the FeB/Fe{sub 2}B/substrate system were also evaluated. Finally, the results obtained provide essential information about the effect of the layer thickness, the residual stresses, and the resilience modulus on the cohesive/adhesive strength in FeB/Fe{sub 2}B systems.

36 MATERIALS SCIENCE↗

DOE OSTI2016

One-step fabrication of submicrostructures by low one-photon absorption direct laser writing technique with local thermal effect

In this work, local thermal effect induced by a continuous-wave laser has been investigated and exploited to optimize the low one-photon absorption (LOPA) direct laser writing (DLW) technique for fabrication of polymer-based microstructures. It was demonstrated that the temperature of excited SU8 photoresist at the focusing area increases to above 100 °C due to high excitation intensity and becomes stable at that temperature thanks to the use of a continuous-wave laser at 532 nm-wavelength. This optically induced thermal effect immediately completes the crosslinking process at the photopolymerized region, allowing obtain desired structures without using the conventional post-exposure bake (PEB) step, which is usually realized after the exposure. Theoretical calculation of the temperature distribution induced by local optical excitation using finite element method confirmed the experimental results. LOPA-based DLW technique combined with optically induced thermal effect (local PEB) shows great advantages over the traditional PEB, such as simple, short fabrication time, high resolution. In particular, it allowed the overcoming of the accumulation effect inherently existed in optical lithography by one-photon absorption process, resulting in small and uniform structures with very short lattice constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Sierra/SolidMechanics 4.54. Capabilities in Development

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 4.54 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as peridynamics and the reproducing kernel particle method (RKPM), numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and /-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Sierra/SolidMechanics 4.58. Capabilities In Development

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 4.58 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as peridynamics and the reproducing kernel particle method (RKPM), numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and /-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

A simple structure of all circular-air-holes photonic crystal fiber for achieving high birefringence and low confinement loss

We propose a simple structure of photonic crystal fibers (PCFs) with high birefringence and low confinement loss based on one rectangular centric ring of smaller circular air holes (CAHs) in the fiber core, and three rings of larger CAHs in the fiber cladding. This simple geometry (using all CAHs with two different air hole sizes) is capable of achieving a flexible control of the birefringence, B = 5.501 × 10{sup −3}, and ultra-low confinement loss, 7.30 × 10{sup −5 }dB/km, at an excitation wavelength of λ = 1550 nm. The birefringence value is ∼5.0 times greater than that obtained for conventional CAH PCF. This simple structure has the added advantage from the view point of easy fabrication, robustness, and cost. A full-vector finite element method combined with anisotropic perfectly matched layers was used to analyze the various fiber structures. We have analyzed four cases of CAH PCFs, focusing on the core asymmetry design as opposed to the conventional approach of CAHs or elliptical air holes on the cladding and core. The robustness against manufacturing inaccuracies of the proposed structure has also been further investigated in this work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1987

Shape design sensitivity analysis and optimal design of structural systems

The material derivative concept of continuum mechanics and an adjoint variable method of design sensitivity analysis are used to relate variations in structural shape to measures of structural performance. A domain method of shape design sensitivity analysis is used to best utilize the basic character of the finite element method that gives accurate information not on the boundary but in the domain. Implementation of shape design sensitivty analysis using finite element computer codes is discussed. Recent numerical results are used to demonstrate the accuracy obtainable using the method. Result of design sensitivity analysis is used to carry out design optimization of a built-up structure.

Choi, Kyung K.↗

NASA NTRS1982

Finite element analysis of aerodynamic heating in three dimensional viscous high speed compressible flow: An assessment

The current capability of the finite element method for solving problems of viscous flow is reviewed. Much work has been directed to the simulation of incompressible flows and the relevant features are described. The methods available for, and the problems associated with, the finite element solution of high speed viscous compressible flows are analyzed. A plan for developing finite element research in this area with experimental support is presented.

Morgan, K.↗

DOE OSTI2015

Localized surface phonon polariton resonances in polar gallium nitride

We demonstrate the excitation of localized surface phonon polaritons in an array of sub-diffraction pucks fabricated in an epitaxial layer of gallium nitride (GaN) on a silicon carbide (SiC) substrate. The array is characterized via polarization- and angle-dependent reflection spectroscopy in the mid-infrared, and coupling to several localized modes is observed in the GaN Reststrahlen band (13.4–18.0 μm). The same structure is simulated using finite element methods and the charge density of the modes are studied; transverse dipole modes are identified for the transverse electric and magnetic polarizations and a quadrupole mode is identified for the transverse magnetic polarization. The measured mid-infrared spectrum agrees well with numerically simulated spectra. This work could enable optoelectronic structures and devices that support surface modes at mid- and far-infrared wavelengths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Band structure of cavity-type hypersonic phononic crystals fabricated by femtosecond laser-induced two-photon polymerization

The phononic band diagram of a periodic square structure fabricated by femtosecond laser pulse-induced two photon polymerization is recorded by Brillouin light scattering (BLS) at hypersonic (GHz) frequencies and computed by finite element method. The theoretical calculations along the two main symmetry directions quantitatively capture the band diagrams of the air- and liquid-filled structure and moreover represent the BLS intensities. The theory helps identify the observed modes, reveals the origin of the observed bandgaps at the Brillouin zone boundaries, and unravels direction dependent effective medium behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

A comprehensive simulation model of the performance of photochromic films in absorbance-modulation-optical-lithography

Optical lithography is the most prevalent method of fabricating micro-and nano-scale structures in the semiconductor industry due to the fact that patterning using photons is fast, accurate and provides high throughput. However, the resolution of this technique is inherently limited by the physical phenomenon of diffraction. Absorbance-Modulation-Optical Lithography (AMOL), a recently developed technique has been successfully demonstrated to be able to circumvent this diffraction limit. AMOL employs a dual-wavelength exposure system in conjunction with spectrally selective reversible photo-transitions in thin films of photochromic molecules to achieve patterning of features with sizes beyond the far-field diffraction limit. We have developed a finite-element-method based full-electromagnetic-wave solution model that simulates the photo-chemical processes that occur within the thin film of the photochromic molecules under illumination by the exposure and confining wavelengths in AMOL. This model allows us to understand how the material characteristics influence the confinement to sub-diffraction dimensions, of the transmitted point spread function (PSF) of the exposure wavelength inside the recording medium. The model reported here provides the most comprehensive analysis of the AMOL process to-date, and the results show that the most important factors that govern the process, are the polarization of the two beams, the ratio of the intensities of the two wavelengths, the relative absorption coefficients and the concentration of the photochromic species, the thickness of the photochromic layer and the quantum yields of the photoreactions at the two wavelengths. The aim of this work is to elucidate the requirements of AMOL in successfully circumventing the far-field diffraction limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1977

Displacement incrementation in non-linear structural analysis by the self-correcting method

The prediction of nonlinear structural behavior by the finite element method wherein buckling does not occur has received considerable attention and, with it, reasonable success has been achieved. However, the post-buckling problem has been less actively pursued probably because of the inherent numerical difficulties encountered. This note reviews very briefly the numerical methods currently being used for pre- and post-buckling analysis and presents a self-correcting approach based on load and displacement incrementation which is shown to be efficient, reliable, and easy to program. Numerical solutions are presented which demonstrate the effectiveness of the method.

Haisler, W. E.↗

NASA NTRS2001

Discontinuous Galerkin Methods for NonLinear Differential Systems

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗

DOE OSTI2016

Evaluation of CFD Methods for Simulation of Two-Phase Boiling Flow Phenomena in a Helical Coil Steam Generator

The U.S. Department of Energy, Office of Nuclear Energy charges participants in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program with the development of advanced modeling and simulation capabilities that can be used to address design, performance and safety challenges in the development and deployment of advanced reactor technology. The NEAMS has established a high impact problem (HIP) team to demonstrate the applicability of these tools to identification and mitigation of sources of steam generator flow induced vibration (SGFIV). The SGFIV HIP team is working to evaluate vibration sources in an advanced helical coil steam generator using computational fluid dynamics (CFD) simulations of the turbulent primary coolant flow over the outside of the tubes and CFD simulations of the turbulent multiphase boiling secondary coolant flow inside the tubes integrated with high resolution finite element method assessments of the tubes and their associated structural supports. This report summarizes the demonstration of a methodology for the multiphase boiling flow analysis inside the helical coil steam generator tube. A helical coil steam generator configuration has been defined based on the experiments completed by Polytecnico di Milano in the SIET helical coil steam generator tube facility. Simulations of the defined problem have been completed using the Eulerian-Eulerian multi-fluid modeling capabilities of the commercial CFD code STAR-CCM+. Simulations suggest that the two phases will quickly stratify in the slightly inclined pipe of the helical coil steam generator. These results have been successfully benchmarked against both empirical correlations for pressure drop and simulations using an alternate CFD methodology, the dispersed phase mixture modeling capabilities of the open source CFD code Nek5000.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2018

The effect of cell-size dispersity on the mechanical properties of closed-cell aluminum foam

Highlights: • The effect of cell-size dispersity on the Young’s modulus of foam was investigated through experiments and simulations • Crushing experiments were conducted on 27 foam specimens in order measure the relative Young’s modulus • Three types of numerical models, including monodispersed, bidispersed and polydispersed models were studied • The numerical convergence on the model size for different models was also studied • The difference among numerical results obtained from the three models was insignificant - Abstract: The mechanical properties of open-celled and closed-celled metallic foam are strongly influenced by the cellular microstructure and material properties of the base material from which the foam is made. The present study investigated the microstructure and mechanical properties of ALPORAS closed-cell Al foam. In particular, the effect of cell-size dispersity on the Young's modulus of foam was investigated through crushing experiments and numerical simulations. Crushing experiments were first conducted on two mutually perpendicular directions. Numerical models of varying cell-size dispersity and complexity were subsequently developed, including the monodispersed model that is based on the Kelvin structure, the bidispersed model, and the polydispersed model that is based on the random Laguerre tessellation algorithm. The numerical models were discretized using shell elements, and the simulations were conducted using the ABAQUS finite element method software. Relative Young's moduli with various relative densities were predicted, and the numerical convergence on the model size was studied accordingly. Numerical results obtained from all three numerical models developed were found to overestimate the experimental results by a factor of approximately 3. The difference between the numerical and experimental results can be attributed to the difference in detailed microstructural characteristics. Furthermore, the difference among numerical results obtained from the three models was insignificant.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Interfacial Strength Evaluation of Oxide Films on Carbon Steel by Using the Laser Shock Adhesion Test

We investigate the adhesion strength (interfacial strength) of the oxide film formed on a carbon steel substrate by using a laser shock adhesion test (LaSAT). A rolled carbon steel (SS400) plate is first heated at high temperature (700 °C) under the atmosphere (heat treatment) to cause the formation of oxide films about 20-50 μm thick on the steel surface. The film consisted of multiple layers, including wustite (FeO), magnetite (Fe{sub 3}O{sub 4}), and hematite (Fe{sub 2}O{sub 3}). Next, the x-ray diffraction method was used to measure the residual stress in the film. We find that compression stress develops in the film because of a mismatch of thermal strain in between the film and substrate after the heat treatment. With LaSAT, strong elastic wave generates by a pulsed laser ablation of the substrate backside, and results in interfacial fractures in the oxide films. To monitor the normal (out-of-plane) displacement of the top surface film, a laser ultrasonic interferometer is simultaneously utilized, such that we detect film delamination and determine the critical laser energy for interfacial fracture to obtain the critical stress (interfacial strength). Furthermore, additional pulsed laser irradiation is conducted against a pre-delaminated specimen in order to examine the growth of delamination area and evaluate the interfacial fracture toughness. Finally, elastic wave propagation in the sample is computed using the finite element method (FEM) in order to evaluate the interfacial stress field, accurately. For the FEM, multiple layers in the oxide film are modeled, and the residual stress due to thermal strain mismatch is introduced. The FEM computations reveal the stress distribution around the interface and evaluate the critical stress and critical stress intensity factor for delamination. We thus quantitatively evaluate the interfacial strength and interfacial fracture toughness of the oxide film formed on carbon steel. The mechanism of interfacial fracture is also discussed on the basis of microstructures in a film with multiple layers. This result may be useful for understanding the adhesion quality of oxide film during the hot-rolling process.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Brillouin light scattering study of transverse mode coupling in confined yttrium iron garnet/barium strontium titanate multiferroic

Using the space-resolved Brillouin light scattering spectroscopy we study the transformation of dynamic magnetization patterns in a bilayer multiferroic structure. We show that in the comparison with a single yttrium iron garnet (YIG) film magnetization distribution is transformed in the bilayer structure due to the coupling of waves propagating both in an YIG film (magnetic layer) and in a barium strontium titanate slab (ferroelectric layer). We present a simple electrodynamic model using the numerical finite element method to show the transformation of eigenmode spectrum of confined multiferroic. In particular, we demonstrate that the control over the dynamic magnetization and the transformation of spatial profiles of transverse modes in magnetic film of the bilayer structure can be performed by the tuning of the wavevectors of transverse modes. The studied confined multiferroic stripe can be utilized for fabrication of integrated dual tunable functional devices for magnonic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗