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.

SourceResearch connectionFrom yearThrough year

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

At least 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

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 OSTI2014

Non destructive technique for cracks detection by an eddy current in differential mode for steel frames

In this paper, a nondestructive technique is used as a tool to control cracks and microcracks in materials. A simulation by a numerical approach such as the finite element method is employed to detect cracks and eventually; to study their propagation using a crucial parameter such as the stress intensity factor. This approach has been used in the aircraft industry to control cracks. Besides, it makes it possible to highlight the defects of parts while preserving the integrity of the controlled products. On the other side, it is proven that the reliability of the control of defects gives convincing results for the improvement of the quality and the safety of the material. Eddy current testing (ECT) is a standard technique in industry for the detection of surface breaking flaws in magnetic materials such as steels. In this context, simulation tools can be used to improve the understanding of experimental signals, optimize the design of sensors or evaluate the performance of ECT procedures. CEA-LIST has developed for many years semi-analytical models embedded into the simulation platform CIVA dedicated to non-destructive testing. The developments presented herein address the case of flaws located inside a planar and magnetic medium. Simulation results are obtained through the application of the Volume Integral Method (VIM). When considering the ECT of a single flaw, a system of two differential equations is derived from Maxwell equations. The numerical resolution of the system is carried out using the classical Galerkin variant of the Method of Moments. Besides, a probe response is calculated by application of the Lorentz reciprocity theorem. Finally, the approach itself as well as comparisons between simulation results and measured data are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

DOE OSTI2018

Modeling Resonances of Dielectric Spheres in the Terahertz Range

We have studied the resonance properties of dielectric spheres with subwavelength dimensions in the optical and terahertz frequency range. Distribution of the electric field vector in a sphere irradiated by a plane polarized wave has been numerically simulated using the finite element method. The obtained results clearly demonstrate that it is in principle possible to use resonances of various orders in dielectric spheres for the amplification of current in a conducting antenna with dimensions several times smaller than the radiation wavelength. In the case of resonances within short spectral intervals, it is possible to provide a current density gain up to 25–30 dB as compared to that in the absence of a dielectric sphere surrounding the antenna.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2016

Influence of radiation on double conjugate diffusion in a porous cavity

The current work highlights the effect of radiation on the conjugate heat and mass transfer in a square porous cavity having a solid wall. The solid wall is placed at the center of cavity. The left surface of cavity is maintained at higher temperature T{sub w} and concentration C{sub w} whereas the right surface is maintained at T{sub c} and C{sub c} such that T{sub w}>T{sub c} and Cw>Cc. The top and bottom surfaces are adiabatic. The governing equations are solved with the help of finite element method by making use of triangular elements. The results are discussed with respect to two different heights of solid wall inside the porous medium along with the radiation parameter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

Diffusion Acceleration Schemes for Self-Adjoint Angular Flux Formulation with a Void Treatment

A Galerkin weak form for the monoenergetic neutron transport equation with a continuous finite element method and discrete ordinate method is developed based on self-adjoint angular flux formulation. This weak form is modified for treating void regions. A consistent diffusion scheme is developed with projection. Correction terms of the diffusion scheme are derived to reproduce the transport scalar flux. A source iteration that decouples the solution of all directions with both linear and nonlinear diffusion accelerations is developed and demonstrated. One-dimensional Fourier analysis is conducted to demonstrate the stability of the linear and nonlinear diffusion accelerations. Numerical results of these schemes are presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2019

Plasmon Enhancement of the Electric Field in Mid-Infrared Ge/Si Quantum-Dot Photodetectors with Different Thicknesses of the Active Region

The spatial distribution of the electric field in Ge/Si photodetector heterostructures coated with a gold film containing a regular two-dimensional array of subwavelength apertures is calculated by the finite-element method. The array period and aperture diameter are 1.2 and 0.7 μm, respectively. The efficiency of field enhancement is determined for different thicknesses of the active region occupied by quantum dots. It is demonstrated that the field-enhancement factor for an electromagnetic wave incident on the structure from the side of the substrate is ~3.5 times larger than that for a wave incident from the opposite side. In the first case, the field-enhancement factor varies nonmonotonically with the thickness of the active region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Influence of AC Loss on Stress and Strain of Superconducting Coils

The second-generation (2G) high-temperature superconducting (HTS) coated conductors (CC) are increasingly used in power systems recently, especially in large-capacity superconducting magnetic energy storage (SMES). HTSCC in superconducting energy storage coil is subjected to thermal stress which is caused by thermal contraction due to AC loss. The thermal stress will degrade the critical currents and may cause the displacement of the superconducting coil and insulating epoxy resin. In this work, the AC loss of the superconducting coil is firstly obtained by finite element method (FEM) software COMSOL with the method of iterative conductivity. Furthermore, the effect of AC loss on the superconducting coil is thus analyzed by solving the coupling problems which include magnetic field, mechanical behaviors, and heat transfer in solids. The result shows that the AC loss of the superconducting coil is concentrated in the middle of the coil and the maximum AC loss density is 5.96 × 10{sup 5} W/m{sup 3}; the maximum stress laid at the inner radius of the coil and the maximum stress is 0.01 MPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Influence of cementation on the yield surface of rocks numerically determined from digital microstructures

Digital Rock Physics has reached a level of maturity on the characterisation of primary properties that depend on the microstructure – such as porosity, permeability or elastic moduli – by numerically solving field equations on μCT scan images of rock. After small deformations or at depth though, most rocks eventually reach their limit of elasticity and the complementary plastic properties are needed to describe the full mechanical behaviour. Currently, determination of a rock’s yield surface from its microstructure is often restricted to semi-analytical criteria derived by limit analysis or numerical simulations performed on idealised geometries. Such simplification lacks representativeness, particularly for processes that affect directly the pore-grain interface such as the cementation phenomenon, happening during diagenesis. Eventually, only direct numerical simulation of elasto-plasticity performed on digitalised microstructures can be used to assess the strength of different cemented materials and its evolution with the alteration of the microstructure. In this study, we provide a comprehensive parametric study on the impact of cementation on rock strength for real microstructures of cemented granular materials. Compared to most previous studies, the whole yield surface is determined numerically (using Finite Element Method) in order to assess the influence of cementation for different stress-paths. The previously known tendency of rock to strengthen with increasing cementation volume is verified. New results on the influence of cement property namely Young’s modulus, friction and cohesion on the rock’s yield surface are explored. The envelopes obtained are compared to the ones obtained by experimental data and existing models. Finally, the framework presented in this study showcases the wider possibility of determining any rock’s or porous material’s yield surface from its microstructure.

36 MATERIALS SCIENCE↗

DOE OSTI2013

Features of the stress-strain state of Si/SiO{sub 2}/Ge heterostructures with germanium nanoislands of a limited density

Within the elastic continuum model, with the use of the finite-element method, the stress-strain state of silicon-germanium heterostructures with semispherical germanium islands grown on an oxidized silicon surface is calculated. It is shown that as the density of islands is increased to limiting values, in the SiGe structure with open quantum dots the value and spatial distribution of the elastic-strain fields significantly change. The results of theoretical calculation allow the heterostructure portions with the maximum variation in the stress-strain state to be determined. The position of such a portions can be controlled by changing the density of islands.

36 MATERIALS SCIENCE↗

DOE OSTI2018

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↗

DOE OSTI2016

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↗

DOE OSTI2019

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↗

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↗