Search NASA⌕ Search

SEARCH · Search NASA

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effective sensor location for detection of change in structural dynamic response

This paper investigates the sensitivity of structural system response to the sensor location by investigating consequences of small changes in the location to the structural system response. The paper discusses how maximum observability (based on mode shape and the participation of that mode in the input provided) drives optimal location. The structural responses were investigated in terms of the g-rms response for various low-frequency inputs (pure sinusoids and real-life inputs such as an earthquake and trains). Results were then analyzed in the context of Modal Contributions Factors (MCF) and changes to the Force-to-response Transfer Functions (TRFs). A modal-matching process is first presented using a Matlab TM -based Finite Element Method (FEM) model of a cantilever beam and instrumentation to determine the location of a small mass based on three different criteria. Subsequently, the structural response is investigated using experiments and the FEM model. The accelerometer of small mass (at 1/3 height) was moved up or down to obtain changes in the structural response (TRF) to various realistic low-frequency inputs. Modal Contribution Factor (MCF) and derivative (slope) of the associated mode-shapes were correlated to the observed changes in TRFs. Results show how optimal sensor locations for detecting change in structural response can be based on the MCFs and the associated mode-shapes.

Woodall, J.↗
DOE OSTI2021

A simplified framework for prediction of sensor network coverage in real-time structural health monitoring of plate-like structures

In this article, the coverage area prediction of piezoelectric sensor network for detecting a specific type of under-surface crack in plate-like structures is addressed. In particular, this article proposes a simplified framework to estimate the coverage of any given sensor network arrangement when a critical defect is known. Based on numerical results from finite element methods (FEM), a simplified framework to estimate coverage area of any given network arrangement is developed. Using such a simplified framework, one can avoid time-consuming procedure of evaluating numerous FEM models in estimating sensor network coverage. Back-scatter fields of partial cracks are estimated using a proposed function, whose parameters are estimated from the results of a limited number of FEM simulations. The proposed function efficiently predicts the back-scattered field of any combination of transmitters and receivers for a given crack geometry. Superposition is used to estimate the coverage area of an arbitrary piezoelectric (e.g., PZT) sensor network. It is shown that the coverage area of a sensor network depends on both sensor network geometry and defect properties (e.g., crack inclination) and it is not necessarily a linear function of the number of sensors. Furthermore, it is shown that the network arrangement has an important effect on the geometry of the coverage area. Experimental results of a network of 14 PZTs in two clusters confirm the accuracy of the method.

Engineering↗
DOE OSTI2019

Computational modeling of continuum scale constitutive equations to improve biomass feedstock material handling and conveying systems

Existing design practices for common feedstock handling equipment such as hoppers and screw feeders have not been successful in many situations. The lack of an appropriate toolset to characterize and model biomass flow behavior is deemed the root cause of plugging, stalls, and flow-related issues in biomass handling systems design and operations. New standardized characterization protocols in conjunction with appropriate constitutive models will lead to successful biomass handling system design and operations. This study investigates constitutive models for bulk solid flow based on the first principles to overcome issues in biomass handling. The focus is on two major issues in biomass handling, i.e., initiating and sustaining of flow in the hopper that is widely used in biomass conveying operation. Accordingly, this article focuses on analytical models of biomass hopper flow. Mohr-Coulomb, Drucker-Prager, and modified Cam-clay models are investigated as constitutive bulk solid flow models. Parameters of those models are determined using a set of triaxial tests. Triaxial tests are conducted using the cubical triaxial tester (CTT) in the Food and Biomaterials Lab of the Pennsylvania State University. Milled 1mm corn stover and 2mm wood chip biomass materials having different flow characteristics are examined, and their respective flow behavior out of a hopper are simulated using the Finite Element method implementing identified constitutive models and determined fundamental mechanical properties. It is found that, with the same hopper geometry, ground corn stover 1mm tends to exhibit funnel flow whereas ground Douglas fir 2 mm tends to exhibit mass flow. These simulation results are validated with hopper flow experiments. This study demonstrates how biomass flow problems can be predicted and mitigated by numerically modeled fundamental constitutive equations.

09 BIOMASS FUELS↗
DOE OSTI2019

Thermal behavior characterization of a kilowatt-power-level cryogenically cooled Yb:YAG active mirror laser amplifier

We present an analysis of the thermal behavior of a high-energy kilowatt-average-power diode-pumped cryogenically cooled Yb:YAG active mirror laser amplifier based on measurements and simulations. Maps of the temperature distribution of the laser material at pump powers up to 1 kW were obtained for the first time by spatially and spectrally resolving the fluorescence induced by a scanning probe beam. The wavefront distortion resulting from the front surface deformation and the overall deformation of the gain medium assembly were measured using a Mach–Zehnder interferometer. The measured deformations agree well with the results of thermomechanical modeling using finite element method simulations, and with the results of focal length shift measurements. The relative contributions to the optical path difference (OPD) of the mechanical deformations, refractive index changes, and electronic contribution are discussed. We propose that the Cr 4+ :YAG cladding plays a significant role in both the temperature distribution and the overall OPD changes. Moreover, the pump-induced mechanical deformations of the assembly dominate the OPD changes in this kilowatt-average-pump-power cryogenically cooled Yb:YAG active mirror laser.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces

The scattering-type scanning near-field optical microscope (s-SNOM) has emerged as a powerful tool for resolving nanoscale inhomogeneities in laterally heterogeneous samples. However, most analytical models used to predict the scattering near-field signals are assuming homogenous landscapes (bulk materials), resulting in inconsistencies when applied to samples with more complex configurations. In this work, we combine the point-dipole model (PDM) to the finite-element method (FEM) to account for the lateral and vertical heterogeneities while keeping the computation time manageable. Full images, spectra, or hyperspectral line profiles can be simulated by calculating the self-consistent dipole radiation demodulated at higher harmonics of the tip oscillation, mimicking real experimental procedures. Using this formalism, we clarify several important yet puzzling experimental observations in near-field images on samples with rich typography and complex material compositions, heterostructures of two-dimensional material flakes, and plasmonic antennas. The developed method serves as a basis for future investigations of nano-systems with nontrivial topography.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Parametric study of optical transmission through plasmonic hole arrays modulated by the phase transition of vanadium dioxide

We have performed comprehensive electromagnetic simulations and preliminary experiments to explore the effects of geometrical and material parameters on the extraordinary optical transmission (EOT) through periodic arrays of subwavelength holes in a bilayer stack consisting of a gold or silver film atop a vanadium dioxide film (Au/Ag + VO 2 ), where the latter undergoes a semiconductor-to-metal phase transition. Using the finite-difference time-domain (FDTD) and finite-element methods (FEM), we vary iteratively the array periodicity, VO 2 film thickness and hole diameters, as well as the refractive index inside the VO 2 -layer holes and the VO 2 optical constants. For each variation, we compare the metallic-to-semiconducting ratios of the zero-order transmission ( T 00 ) peaks and find sharp maxima in these ratios within narrow parameter ranges. The maxima arise from Fabry-Perot and Fano-type resonances that minimize T 00 in the semiconducting phase of the perforated bilayers. At a fixed array period, the primary factors controlling the VO 2 -enabled EOT modulation are the VO 2 thickness, diameter of the VO 2 -layer holes, and absorption in the two VO 2 phases. Besides uncovering the origins of the higher metallic-phase T 00 , this study provides a protocol for optimizing the performance of the bilayer hole arrays for potential uses as dynamically tunable nano-optical devices.

Donev, Eugenii U. (ORCID:000000030555219X)↗
DOE OSTI2021

Self-assembled multifunctional nanostructures for surface passivation and photon management in silicon photovoltaics

This work reports the fabrication and characterization of multifunctional, nanostructured passivation layers formed using a self-assembly process that provide both surface passivation and improved light trapping in crystalline silicon photovoltaic (PV) cells. Scalable block copolymer self-assembly and vapor phase infiltration processes are used to form arrays of aluminum oxide nanostructures (Al 2 O 3 ) on crystalline silicon without substrate etching. The Al 2 O 3 nanostructures are characterized using scanning electron microscopy (SEM), transmission electron microscopy (TEM), and spectroscopic ellipsometry. Injection-level dependent photoconductance measurements are used to determine the effective carrier lifetime of the samples to confirm the nanostructures successfully passivate the Si surface. Finite element method simulations and reflectance measurement show that the nanostructures increase the internal rear reflectance of the PV cell by suppressing the parasitic optical losses in the metal contact. An optimized morphology of the structures is identified for their potential use in PV cells as multifunctional materials providing surface passivation, photon management, and carrier transport pathways.

14 SOLAR ENERGY↗
DOE OSTI2015

Indentation-derived elastic modulus of multilayer thin films: Effect of unloading induced plasticity

Nanoindentation is useful for evaluating the mechanical properties, such as elastic modulus, of multilayer thin film materials. A fundamental assumption in the derivation of the elastic modulus from nanoindentation is that the unloading process is purely elastic. In this work, the validity of elastic assumption as it applies to multilayer thin films is studied using the finite element method. The elastic modulus and hardness from the model system are compared to experimental results to show validity of the model. Plastic strain is shown to increase in the multilayer system during the unloading process. Additionally, the indentation-derived modulus of a monolayer material shows no dependence on unloading plasticity while the modulus of the multilayer system is dependent on unloading-induced plasticity. Lastly, the cyclic behavior of the multilayer thin film is studied in relation to the influence of unloading-induced plasticity. Furthermore, it is found that several cycles are required to minimize unloading-induced plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Adamantine 1.0: A Thermomechanical Simulator for Additive Manufacturing

Adamantine is a thermomechanical simulation code that is written in C++ and built on top of deal.II (Arndt et al., 2023), p4est (Burstedde et al., 2011), ArborX (Lebrun-Grandié et al., 2020), Trilinos (The Trilinos Project Team, 2020), and Kokkos (Trott et al., 2022). Adamantine was developed with additive manufacturing in mind and it is particularly well adapted to simulate fused filament fabrication, directed energy deposition, and powder bed fusion. Adamantine employs the finite element method with adaptive mesh refinement to solve a nonlinear anisotropic heat equation, enabling support for various additive manufacturing processes. It can also perform elastoplastic and thermoelastoplastic simulations. It can handle materials in three distinct phases (solid, liquid, and powder) to accurately reflect the physical state during different stages of the manufacturing process. To enhance simulation accuracy, adamantine incorporates data assimilation techniques (Asch et al., 2016). This allows it to integrate experimental data from sensors like thermocouples and infrared (IR) cameras. This combined approach helps account for errors arising from input parameters, material properties, models, and numerical calculations, leading to more realistic simulations that reflect what occurs in a particular print.

36 MATERIALS SCIENCE↗
DOE OSTI2011

EMPHASIS/Nevada UTDEM user guide. Version 2.0.

The Unstructured Time-Domain ElectroMagnetics (UTDEM) portion of the EMPHASIS suite solves Maxwell's equations using finite-element techniques on unstructured meshes. This document provides user-specific information to facilitate the use of the code for applications of interest. UTDEM is a general-purpose code for solving Maxwell's equations on arbitrary, unstructured tetrahedral meshes. The geometries and the meshes thereof are limited only by the patience of the user in meshing and by the available computing resources for the solution. UTDEM solves Maxwell's equations using finite-element method (FEM) techniques on tetrahedral elements using vector, edge-conforming basis functions. EMPHASIS/Nevada Unstructured Time-Domain ElectroMagnetic Particle-In-Cell (UTDEM PIC) is a superset of the capabilities found in UTDEM. It adds the capability to simulate systems in which the effects of free charge are important and need to be treated in a self-consistent manner. This is done by integrating the equations of motion for macroparticles (a macroparticle is an object that represents a large number of real physical particles, all with the same position and momentum) being accelerated by the electromagnetic forces upon the particle (Lorentz force). The motion of these particles results in a current, which is a source for the fields in Maxwell's equations.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2011

EMPHASIS/Nevada CABANA user Guide. Version 2.0.

The CABle ANAlysis (CABANA) portion of the EMPHASIS{trademark} suite is designed specifically for the simulation of cable system-generated electromagnetic pulse (SGEMP). The code can be used to evaluate the response of a specific cable design to threat or to compare and minimize the relative response of difference designs. This document provides user-specific information to facilitate the application of the code to cables of interest. It solves the electrical portion of a cable SGEMP simulation. It takes specific results from the deterministic radiation-transport code CEPTRE as sources and computes the resulting electrical response to an arbitrary cable load. The cable geometry itself is also arbitrary and is limited only by the patience of the user in meshing and by the available computing resources for the solution. The CABANA simulation involves solution of the quasi-static Maxwell equations using finite-element method (FEM) techniques.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2012

h and hp -adaptive finite elements with local ALE for modeling turbulent reactive flow in engines with injection

LANL and our collaborators are facilitating the effort for predictive internal combustion engine modeling. We are developing an hp -adaptive finite element methods (FEM) software that provides a high degree of accuracy and solution robustness. This FEM projection method, a predictor-corrector scheme (PCS) has excellent capability over all flow regimes. The system is applicable to Newtonian and non-Newtonian fluids, turbulent reactive flows with sprays and using a new local ALE method for immersed moving parts. These algorithms are capable of representing the physics within an engine. T-3 and its collaborators provide software that many use directly for engine simulation or that they may alter with various models. KIVA software is in use worldwide by hundreds of Universities and dozens of corporations involved in research, engine design and manufacturing. We are working to deliver a more predictive capability to help in the understanding of combustion engines, thereby providing for greater efficiency and lower unwanted emissions.

02 PETROLEUM↗
DOE OSTI2012

Logistical simulation of spent nuclear fuel disposal in a salt repository with low temperature limits

Two logistical modeling cases evaluated the relative cost, schedule, and thermal results for strategies that would limit peak salt temperature in a spent-fuel repository to 90°C and 150°C. Runs with the code CALVIN 4.0 show that these low thermal limits are logistically feasible for a salt repository using relatively small waste packages (4-PWR/9-BWR or 12-PWR/24-BWR sizes). Attainment of low-thermal goals was represented using instantaneous waste package thermal power limits at emplacement (3,200 W and 6,800 W) which came from a correlation between package power and peak salt temperature that was developed by a previous study using the finite element method.

Kalinina, Elena Arkadievna↗
DOE OSTI2013

Modeling electron transport in the presence of electric and magnetic fields.

This report describes the theoretical background on modeling electron transport in the presence of electric and magnetic fields by incorporating the effects of the Lorentz force on electron motion into the Boltzmann transport equation. Electromagnetic fields alter the electron energy and trajectory continuously, and these effects can be characterized mathematically by differential operators in terms of electron energy and direction. Numerical solution techniques, based on the discrete-ordinates and finite-element methods, are developed and implemented in an existing radiation transport code, SCEPTRE.

Fan, Wesley C.↗
DOE OSTI2013

SciDAC Center for Gyrokinetic Particle Simulation of Turbulent Transport in Burning Plasmas

During the first year of the SciDAC gyrokinetic particle simulation (GPS) project, the GPS team (Zhihong Lin, Liu Chen, Yasutaro Nishimura, and Igor Holod) at the University of California, Irvine (UCI) studied the tokamak electron transport driven by electron temperature gradient (ETG) turbulence, and by trapped electron mode (TEM) turbulence and ion temperature gradient (ITG) turbulence with kinetic electron effects, extended our studies of ITG turbulence spreading to core-edge coupling. We have developed and optimized an elliptic solver using finite element method (FEM), which enables the implementation of advanced kinetic electron models (split-weight scheme and hybrid model) in the SciDAC GPS production code GTC. The GTC code has been ported and optimized on both scalar and vector parallel computer architectures, and is being transformed into objected-oriented style to facilitate collaborative code development. During this period, the UCI team members presented 11 invited talks at major national and international conferences, published 22 papers in peer-reviewed journals and 10 papers in conference proceedings. The UCI hosted the annual SciDAC Workshop on Plasma Turbulence sponsored by the GPS Center, 2005-2007. The workshop was attended by about fifties US and foreign researchers and financially sponsored several gradual students from MIT, Princeton University, Germany, Switzerland, and Finland. A new SciDAC postdoc, Igor Holod, has arrived at UCI to initiate global particle simulation of magnetohydrodynamics turbulence driven by energetic particle modes. The PI, Z. Lin, has been promoted to the Associate Professor with tenure at UCI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Technical Spotlight: NEAMS Structural Mechanics with Diablo

The Diablo code being developed at Lawrence Livermore National Laboratory (LLNL) uses implicit, Lagrangian finite element methods for the simulation of solid mechanics and multi-physics events over moderate to long time frames. A recent effort examined technologies needed to extend Diablo for use in seismic response incorporating soil-structure interaction (SSI). Our SSI study was intended to illustrate how advanced time-domain simulation approaches compare with and complement the standard frequency-domain approach for SSI.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Failure Predictions for VHTR Core Components using a Probabilistic Contiuum Damage Mechanics Model

The proposed work addresses the key research need for the development of constitutive models and overall failure models for graphite and high temperature structural materials, with the long-term goal being to maximize the design life of the Next Generation Nuclear Plant (NGNP). To this end, the capability of a Continuum Damage Mechanics (CDM) model, which has been used successfully for modeling fracture of virgin graphite, will be extended as a predictive and design tool for the core components of the very high- temperature reactor (VHTR). Specifically, irradiation and environmental effects pertinent to the VHTR will be incorporated into the model to allow fracture of graphite and ceramic components under in-reactor conditions to be modeled explicitly using the finite element method. The model uses a combined stress-based and fracture mechanics-based failure criterion, so it can simulate both the initiation and propagation of cracks. Modern imaging techniques, such as x-ray computed tomography and digital image correlation, will be used during material testing to help define the baseline material damage parameters. Monte Carlo analysis will be performed to address inherent variations in material properties, the aim being to reduce the arbitrariness and uncertainties associated with the current statistical approach. The results can potentially contribute to the current development of American Society of Mechanical Engineers (ASME) codes for the design and construction of VHTR core components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

EMPHASIS(TM)/Nevada Unstructured FEM Implementation Version 2.1.1.

EMPHASIS TM /NEVADA is the SIERRA/NEVADA toolkit implementation of portions of the EMP HASIS TM code suite. The purpose of the toolkit i m- plementation is to facilitate coupling to other physics drivers such as radi a- tion transport as well as to better manage code design, implementation, co m- plexity, and important verification and validation processes. This document describes the theory and implementation of the unstructured finite - element method solver , associated algorithms, and selected verification and valid a- tion . Acknowledgement The author would like to recognize all of the ALEGRA team members for their gracious and willing support through this initial Nevada toolkit - implementation process. Although much of the knowledge needed was gleaned from document a- tion and code context, they were always willing to consult personally on some of the less obvious issues and enhancements necessary.

Turner, C. David↗