Engineering Papers⌕ Search

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.

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

At least 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Sierra/SolidMechanics 5.6 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 5.6 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 the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-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 OSTI2022

Sierra/SolidMechanics 5.8 In-Development Manual

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 5.8 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 the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

42 ENGINEERING↗
DOE OSTI2022

Sierra/SolidMechanics 5.10 In-Development Manual

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 5.10 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 the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

42 ENGINEERING↗
DOE OSTI2023

Sierra/SolidMechanics 5.16 Capabilities in Development Manual

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 5.16 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 the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-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 OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1982

Creep crack-growth: A new path-independent integral (T sub c), and computational studies

The development of valid creep fracture criteria is considered. Two path-independent integral parameters which show some degree of promise are the C* and (Delta T)sub c integrals. The mathematical aspects of these parameters are reviewed by deriving generalized vector forms of the parameters using conservation laws which are valid for arbitrary, three dimensional, cracked bodies with crack surface tractions (or applied displacements), body forces, inertial effects, and large deformations. Two principal conclusions are that (Delta T)sub c has an energy rate interpretation whereas C* does not. The development and application of fracture criteria often involves the solution of boundary/initial value problems associated with deformation and stresses. The finite element method is used for this purpose. An efficient, small displacement, infinitesimal strain, displacement based finite element model is specialized to two dimensional plane stress and plane strain and to power law creep constitutive relations. A mesh shifting/remeshing procedure is used for simulating crack growth. The model is implemented with the quartz-point node technique and also with specially developed, conforming, crack-tip singularity elements which provide for the r to the n-(1+n) power strain singularity associated with the HRR crack-tip field. Comparisons are made with a variety of analytical solutions and alternate numerical solutions for a number of problems.

Stonesifer, R. B.↗
NASA NTRS1986

A domain method for shape design sensitivity analysis of built-up structures

A numerical method of obtaining accurate shape design sensitivity information for built-up structures is developed and demonstrated through examples. The method utilizes the basic character of the finite element method that gives accurate information not on the boundary, but in the domain. The method is shown to drastically simplify the derivation of shape design sensitivity formulas for complex built-up structures.

Choi, Kyung K.↗
NASA NTRS1988

Semi-discrete approximations to nonlinear systems of conservation laws; consistency and L(infinity)-stability imply convergence

A convergence theory for semi-discrete approximations to nonlinear systems of conservation laws is developed. It is shown, by a series of scalar counter-examples, that consistency with the conservation law alone does not guarantee convergence. Instead, a notion of consistency which takes into account both the conservation law and its augmenting entropy condition is introduced. In this context it is concluded that consistency and L(infinity)-stability guarantee for a relevant class of admissible entropy functions, that their entropy production rate belongs to a compact subset of H(loc)sup -1 (x,t). One can now use compensated compactness arguments in order to turn this conclusion into a convergence proof. The current state of the art for these arguments includes the scalar and a wide class of 2 x 2 systems of conservation laws. The general framework of the vanishing viscosity method is studied as an effective way to meet the consistency and L(infinity)-stability requirements. How this method is utilized to enforce consistency and stability for scalar conservation laws is shown. In this context we prove, under the appropriate assumptions, the convergence of finite difference approximations (e.g., the high resolution TVD and UNO methods), finite element approximations (e.g., the Streamline-Diffusion methods) and spectral and pseudospectral approximations (e.g., the Spectral Viscosity methods).

Tadmor, Eitan↗
NASA NTRS1994

Stress intensity factors for surface and corner cracks emanating from a wedge-loaded hole

To assist analysis of riveted lap joints, stress intensity factors are determined for surface and corner cracks emanating from a wedge-loaded hole by using a 3-D weight function method in conjunction with a 3-D finite element method. A stress intensity factor equation for surface cracks is also developed to provide a closed-form solution. The equation covers commonly-encountered geometrical ranges and retains high accuracy over the entire range.

Zhao, W.↗
NASA NTRS2016

Extracting Damping Ratio from Dynamic Data and Numerical Solutions

There are many ways to extract damping parameters from data or models. This Technical Memorandum provides a quick reference for some of the more common approaches used in dynamics analysis. Described are six methods of extracting damping from data: the half-power method, logarithmic decrement (decay rate) method, an autocorrelation/power spectral density fitting method, a frequency response fitting method, a random decrement fitting method, and a newly developed half-quadratic gain method. Additionally, state-space models and finite element method modeling tools, such as COMSOL Multiphysics (COMSOL), provide a theoretical damping via complex frequency. Each method has its advantages which are briefly noted. There are also likely many other advanced techniques in extracting damping within the operational modal analysis discipline, where an input excitation is unknown; however, these approaches discussed here are objective, direct, and can be implemented in a consistent manner.

Casiano, M. J.↗
DOE OSTI2017

Multiresolution molecular mechanics: Surface effects in nanoscale materials

Surface effects have been observed to contribute significantly to the mechanical response of nanoscale structures. The newly proposed energy-based coarse-grained atomistic method Multiresolution Molecular Mechanics (MMM) (Yang, To (2015), ) is applied to capture surface effect for nanosized structures by designing a surface summation rule SR{sup S} within the framework of MMM. Combined with previously proposed bulk summation rule SR{sup B}, the MMM summation rule SR{sup MMM} is completed. SR{sup S} and SR{sup B} are consistently formed within SR{sup MMM} for general finite element shape functions. Analogous to quadrature rules in finite element method (FEM), the key idea to the good performance of SR{sup MMM} lies in that the order or distribution of energy for coarse-grained atomistic model is mathematically derived such that the number, position and weight of quadrature-type (sampling) atoms can be determined. Mathematically, the derived energy distribution of surface area is different from that of bulk region. Physically, the difference is due to the fact that surface atoms lack neighboring bonding. As such, SR{sup S} and SR{sup B} are employed for surface and bulk domains, respectively. Two- and three-dimensional numerical examples using the respective 4-node bilinear quadrilateral, 8-node quadratic quadrilateral and 8-node hexahedral meshes are employed to verify and validate the proposed approach. It is shown that MMM with SR{sup MMM} accurately captures corner, edge and surface effects with less 0.3% degrees of freedom of the original atomistic system, compared against full atomistic simulation. The effectiveness of SR{sup MMM} with respect to high order element is also demonstrated by employing the 8-node quadratic quadrilateral to solve a beam bending problem considering surface effect. In addition, the introduced sampling error with SR{sup MMM} that is analogous to numerical integration error with quadrature rule in FEM is very small. - Highlights: • Surface effect captured by Multiresolution Molecular Mechanics (MMM) is presented. • A novel surface summation rule within the framework of MMM is proposed. • Surface, corner and edges effects are accuterly captured in two and three dimension. • MMM with less 0.3% degrees of freedom of atomistics reproduces atomistic results.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS2012

Using Polarimetric SAR Data to Infer Soil Moisture from Surfaces with Varying Subsurface Moisture Profiles

A time-series approach is used to estimate the moisture content-based on polarimetric SAR data. It is found that under the assumption of constant soil moisture, empirically observed relationships between radar backscatter and moisture are only half as sensitive to moisture as compared to actual radar data. A numerical finite element method is used to calculate the radar backscatter for rough soils with arbitrarily varying soil moisture as a function of depth. Several instance of drying and wetting moisture profiles are considered and the radar backscatter is calculated in each case. Radar backscatter is found to crucially depend on the soil moisture variation in the top half wavelength of soil.

subsurface sensing↗
DOE OSTI2015

Direct measurement of electron beam quality conversion factors using water calorimetry

Purpose: In this work, the authors describe an electron sealed water calorimeter (ESWcal) designed to directly measure absorbed dose to water in clinical electron beams and its use to derive electron beam quality conversion factors for two ionization chamber types. Methods: A functioning calorimeter prototype was constructed in-house and used to obtain reproducible measurements in clinical accelerator-based 6, 9, 12, 16, and 20 MeV electron beams. Corrections for the radiation field perturbation due to the presence of the glass calorimeter vessel were calculated using Monte Carlo (MC) simulations. The conductive heat transfer due to dose gradients and nonwater materials was also accounted for using a commercial finite element method software package. Results: The relative combined standard uncertainty on the ESWcal dose was estimated to be 0.50% for the 9–20 MeV beams and 1.00% for the 6 MeV beam, demonstrating that the development of a water calorimeter-based standard for electron beams over such a wide range of clinically relevant energies is feasible. The largest contributor to the uncertainty was the positioning (Type A, 0.10%–0.40%) and its influence on the perturbation correction (Type B, 0.10%–0.60%). As a preliminary validation, measurements performed with the ESWcal in a 6 MV photon beam were directly compared to results derived from the National Research Council of Canada (NRC) photon beam standard water calorimeter. These two independent devices were shown to agree well within the 0.43% combined relative uncertainty of the ESWcal for this beam type and quality. Absorbed dose electron beam quality conversion factors were measured using the ESWcal for the Exradin A12 and PTW Roos ionization chambers. The photon-electron conversion factor, k{sub ecal}, for the A12 was also experimentally determined. Nonstatistically significant differences of up to 0.7% were found when compared to the calculation-based factors listed in the AAPM’s TG-51 protocol. General agreement between the relative electron energy dependence of the PTW Roos data measured in this work and a recent MC-based study are also shown. Conclusions: This is the first time that water calorimetry has been successfully used to measure electron beam quality conversion factors for energies as low as 6 MeV (R{sub 50} = 2.25 cm)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Enhanced sensitivity of hemoglobin sensor using dual-core photonic crystal fiber

A Hemoglobin (Hb) biosensor based on dual-core photonic crystal fiber is proposed and analyzed. In this paper, the effective refractive index dependency on the Hb concentration within a blood sample is utilized in obtaining the transmission spectrum variation. The results are calculated using full vectorial finite element method. Through this study, the effect of the structure geometrical parameters on the sensor performance is optimized to maximize the sensor sensitivity. The numerical results show a sensitivity of 8.013 nm/g/dL for the X-polarized mode at 3.7 cm fiber length and 7.68 nm/g/dL for the Y-polarized mode at 3.2 cm fiber length of the proposed sensor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Characterization of latent image of electron beam resist via critical-dimension resonant soft X-ray scattering

Metrology plays a crucial role in semiconductor manufacturing by providing accurate and precise measurement and characterization of critical parameters. With the development of high-resolution extreme ultraviolet lithography (EUVL) processes, critical dimensions are shrinking to sub-10 nm. Resist materials encounter the challenge of providing heightened sensitivity and a handle on exacerbating stochastic variations. A comprehensive understanding of the chemical profile of the latent image is pivotal for mitigating stochastic effects and optimizing pattern quality. However, the subtle differences in chemistry between the exposed and unexposed regions of the resists make it extremely challenging to characterize the latent images with sub-nanometer precision. Here, we develop the metrology with critical-dimension resonant soft X-ray scattering (CD-RSoXS) to probe the chemical profiles of latent images stored in resist after exposure. The combination of absorption spectroscopy and enhanced scattering contrast makes it possible to characterize the subtle structural and chemical variations in the latent image. Moreover, the results of the measurements are compared with the simulations with a finite element method–based Maxwell solver to extract a detailed profile of the latent and developed images. We demonstrate that the CD-RSoXS technique can provide valuable insights into the high spatial resolution and local chemical sensitivity simultaneously, which is crucial to understanding the resolution limits and stochastic effects in EUVL processes.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Adhesive-Bonded Composite Joint Analysis with Delaminated Surface Ply Using Strain-Energy Release Rate

This paper presents an analytical model to determine the strain energy release rate due to an interlaminar crack of the surface ply in adhesively bonded composite joints subjected to axial tension. Single-lap shear-joint standard test specimen geometry with thick bondline is followed for model development. The field equations are formulated by using the first-order shear-deformation theory in laminated plates together with kinematics relations and force equilibrium conditions. The stress distributions for the adherends and adhesive are determined after the appropriate boundary and loading conditions are applied and the equations for the field displacements are solved. The system of second-order differential equations is solved to using the symbolic computation tool Maple 9.52 to provide displacements fields. The equivalent forces at the tip of the prescribed interlaminar crack are obtained based on interlaminar stress distributions. The strain energy release rate of the crack is then determined by using the crack closure method. Finite element analyses using the J integral as well as the crack closure method are performed to verify the developed analytical model. It has been shown that the results using the analytical method correlate well with the results from the finite element analyses. An attempt is made to predict the failure loads of the joints based on limited test data from the literature. The effectiveness of the inclusion of bondline thickness is justified when compared with the results obtained from the previous model in which a thin bondline and uniform adhesive stresses through the bondline thickness are assumed.

Chadegani, Alireza↗
DOE OSTI2019

Experimental and Numerical Investigations of Cyclic Plastic Deformation of Al-Mg Alloy

In this work, low-cycle fatigue (LCF) tests are performed on as-received and solution-treated Al-Mg alloy (AA 5754). Solution heat treatment of the as-received alloy is performed at 530 °C for 2 h and then quenched in water at room temperature. The low-cycle fatigue response of the as-received and solution-treated alloy is evaluated at five different strain amplitudes from 0.2 to 0.6%. The strain rate of 0.5 × 10{sup −3} s{sup −1} is kept constant in all the LCF tests. The LCF behavior of AA 5754 is modeled by finite element method (FEM) using Chaboche kinematic hardening cyclic plasticity model. Hardening behavior or fatigue life of as-received and solution-treated alloy is successfully simulated using this model. The results obtained by FEM simulations are compared with those obtained by LCF experiments. Microstructural studies are also performed for as-received and solution-treated alloy. The fracture surfaces obtained from LCF tests are analyzed by scanning electron microscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effect of viscous dissipation and radiation in an annular cone

The viscous dissipation is an effect due to which heat is generated inside the medium. The presence of radiation further complicates the heat transfer behavior inside porous medium. The present paper discusses the combined effect of viscous dissipation and radiation inside a porous medium confined in an annular cone with inner radius r{sub i}. The viscous dissipation and radiation terms are included in the energy equation thereby solving the coupled momentum and energy equations with the help of finite element method. The results are presented in terms of isothermal and streamline indicating the thermal and fluid flow behavior of porous medium. It is found that the combination of viscous dissipation and radiation parameter and the cone angle has significant effect on the heat transfer and fluid flow behavior inside the porous medium. The fluid velocity is found to increase with the increase in Raleigh number.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗