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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

From force chains to nonclassical nonlinear dynamics in cemented granular materials

In this letter, we present evidence for a mechanism responsible for the nonclassical nonlinear dynamics observed in many cemented granular materials that are generally classified as mesoscopic nonlinear elastic materials. We demonstrate numerically that force chains are created within the complex grain-pore network of these materials when subjected to dynamic loading. The interface properties between grains along with the sharp and localized increase of the stress occurring at the grain-grain contacts leads to a reversible decrease of the elastic properties at macroscopic scale and peculiar effects on the propagation of elastic waves when grain boundary properties are appropriately considered. These effects are observed for relatively small amplitudes of the elastic waves, i.e., within tens of microstrain, and relatively large wavelengths, i.e., orders of magnitude larger than the material constituents. The mechanics are investigated numerically using the hybrid finite-discrete-element method and match those observed experimentally using nonlinear resonant ultrasound spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

In situ measurement and modelling of the growth and length scale of twins in α -uranium

Twinning is an important deformation mechanism in crystalline materials. Twins nucleate, then grow, forming a discrete pattern of deformation bands with a certain thickness. The mechanism for twin growth is not completely understood. We present a simple phase field model which captures the physics of twin growth and reproduces the details observed during in situ electron backscatter diffraction experiments on α-uranium. The interaction between residual dislocations at twin interfaces and mobile dislocations in untwinned regions increases the stress needed for twin growth. This is described by a nonlocal term in the proposed constitutive equations: the nucleation stress of a twin increases proportionally to the twin phase field in a neighborhood. Competition between slip and twinning favors the nucleation of a new twin rather than the growth of a preexisting twin. The phase field model is coupled with a crystal plasticity finite element solver that includes the plastic deformation due to twinning. Simulations are able to reproduce the number of twins and their thicknesses as a function of the strain observed during in situ electron backscatter diffraction experiments. The stress-strain curve is also reproduced. This comparison allows the values of the nucleation stress and the interaction strength between residual and mobile dislocations to be found. It also gives an estimation of the density of residual dislocations that is consistent with the observed twin thickness. Finally, this model can be applied to understand microstructural effects in materials using twinning as a strengthening mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal plasticity finite element simulation of lattice rotation and x-ray diffraction during laser shock compression of tantalum

Here we present a crystal plasticity model tailored for high-pressure, high-strain-rate conditions that uses a multiscale treatment of dislocation-based slip kinetics. We use this model to analyze the pronounced plasticity-induced lattice rotations observed in shock-compressed polycrystalline tantalum via in situ x-ray diffraction. By making direct comparisons between experimentally measured and simulated texture evolution, we can explain how the details of the underlying slip kinetics control the degree of lattice rotation that ensues. Specifically, we show that only the highly nonlinear kinetics caused by dislocation nucleation can explain the magnitude of the rotation observed under shock compression. We demonstrate a good fit between our crystal plasticity model and x-ray diffraction data and exploit the data to quantify the dislocation nucleation rates that are otherwise poorly constrained by experiment in the dynamic compression regime.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evaluation of Surface Crack Formation in Photovoltaic Backsheets Using Fragmentation and Finite Element Simulations

Backsheet cracking is among the most commonly observed degradation modes of photovoltaic modules in the field. Cracks can reduce the ability of backsheets to fulfill their functions, such as protection of the modules from the environment or electrical insulation. Here, this work presents an evaluation of the degradation and cracking propensity of two backsheets during accelerated IEC TS 62788-7-2 (International Electrotechnical Commission) aging under the A3 condition with a fragmentation test using a coextruded polyamide backsheet (AAA), and a laminated multilayer backsheet with a polyethylene terephthalate core and outer layer and ethylene vinyl acetate inner layer (PPE). Results show the surface embrittlement of the AAA outer layer during exposure. A longer time of exposure causes cracks to form at lower strains during stretching and creates deeper cracks. Accordingly, Young's modulus of the outer layer increases, as measured by cross-sectional nanoindentation. PPE exhibits cracking after exposure as well. While no similar increase of modulus or crack depth can be observed, the outer layer of PPE exhibits more obvious signs of erosion during exposure, including progressive morphological changes and thickness losses. A finite element model is developed to simulate surface crack formation, based on initially zero-thickness decohesion elements. Decohesion criteria define a critical stress, at which these elements grow, and cracks can begin to form. These criteria are obtained via parameter optimization by comparison between simulation and experiment. The model is used to interpret crack formation in both backsheets.

14 SOLAR ENERGY↗
DOE OSTI2021

Computational Modeling of Photovoltaic Mini-Modules Undergoing Accelerated Stress Testing

A finite element model of a four-cell photovoltaic mini-module was developed and compared to experimental results from an accelerated stress test protocol in order to validate that computational models can accurately represent their physical counterparts when subjected to mechanical loading and to assess mini-module representativeness against full scale photovoltaic modules. Deflected shapes across the simulated mini-modules were compared to measured mini-module shapes when subjected to various pressure loads. Displaced mini-module shape results constrained to the experimental protocols of 0.4 mm and 1.1 mm of displacement at the mini-module center were compared to experimental results of full-size modules subjected to module qualification test load levels of 1.0 kPa and 2.4 kPa, to assess if the bending of mini-modules was representative of full-sized modules under the load. Temperature cycling was incorporated into the model to simulate the impacts of stress due to thermal expansion of the backsheet and cells. A preliminary uncertainty analysis was performed to show how variations in material properties and geometric parameters change the simulation results.

computational modeling↗
DOE OSTI2025

Design of a Structure for Assembly and Cooling the Magnet of the Next-Generation 45 GHz ECR Ion Source MARS-D

Here, the current Electron Cyclotron Resonance Ion Sources (ECRISs), constructed with Nb-Ti wires and the conventional racetrack-and-solenoid structure, have achieved operating frequencies up to 28 GHz and utilized about 90% of the critical current of the Nb-Ti wire. A Mixed Axial and Radial field System Demonstrator (MARS-D) is being developed at Lawrence Berkeley National Laboratory (LBNL). This system, which consists of an innovative hexagonal Closed-Loop Coil (CLC) and a set of solenoids, can generate higher magnetic fields (up to 150% ) while requiring only about 50% of the superconducting wire, enabling Nb-Ti wires to be used in the next-generation 45 GHz ECRIS. However, the assembly and cooling of such an efficient and compact magnet are particularly challenging due to the small radial gap between the CLC and solenoids, as well as the tight operating temperature margin. To address these challenges, a structure was developed that combines a three-section radially split solenoid mandrel with a series of shrink-fit reinforcement rings and cooling channels. This paper presents the detailed structure, manufacturing method, assembly procedure, impregnation method, mechanical Finite Element Analysis (FEA) comparison, and thermal FEA comparison.

ECR ion sources↗
DOE OSTI2021

Microstructure based modeling of friction stir welded joint of dissimilar metals using crystal plasticity

We report the friction stir welding (FSW) process shows promising results in joining dissimilar metals which are otherwise almost impossible to join using traditional welding techniques. Being a new technique, the deformation and the failure mechanism of the joints made by the FSW process needs to be investigated. In this work, a joint between AZ31 Mg alloy and DP590 steel is modeled using phenomenological crystal plasticity formulation on the mesoscale in the form of a representative volume element (RVE). The interface of the two materials is modeled using a cohesive zone model. A parametric study has been performed to understand the effect of grain size and interface fracture toughness as well as strength on the mechanical performance of the joint. It was found that the grain size of AZ31 Mg alloy, as well as DP590 steel, has little effect on the overall joint performance. On the other hand, interface fracture toughness and strength have a significant impact on the mechanical properties of the joint

36 MATERIALS SCIENCE↗
DOE OSTI2021

AIR Algebraic Multigrid for a Space-Time Hybridizable Discontinuous Galerkin Discretization of Advection(-Diffusion)

This article investigates the efficiency, robustness, and scalability of approximate ideal restriction (AIR) algebraic multigrid as a preconditioner in the all-at-once solution of a space-time hybridizable discontinuous Galerkin discretization of advection-dominated flows. The motivation for this study is that the time-dependent advection-diffusion equation can be seen as a “steady” advection-diffusion problem in $(d+1)$-dimensions and AIR has been shown to be a robust solver for steady advection-dominated problems. Numerical examples demonstrate the effectiveness of AIR as a preconditioner for advection-diffusion problems on fixed and time-dependent domains, using both slab-by-slab and all-at-once space-time discretizations, and in the context of uniform and space-time adaptive mesh refinement. A closer look at the geometric coarsening structure that arises in AIR also explains why AIR can provide robust, scalable, space-time convergence on advective and hyperbolic problems, while most multilevel parallel-in-time schemes struggle with such problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Contributions to MoDELib SOFTWARE

The purpose of the current request is to enable LANL employees to contribute computer source code to the existing public repository of the MoDELib software package. This software implements discrete dislocation dynamics (DDD) and finite element (FEM) methods and is currently a vital component of an ongoing DR project at LANL, in collaboration with its original author and maintainer Giacomo Po. Contributions from LANL employees would aim to enhance the reliability, accuracy, and performance of MoDELib simulations using LANL's high performance computing platforms through bug fixes, algorithmic refinements, and parallelization.

Julian, Nicholas↗
DOE OSTI2026

HygroThermFEM v1.0

HygroThermFEM is a Finite Element Method-based numerical calculation engine for solving 2-D heat and moisture transfer problems. This numerical engine is used in the THERM software tool, and its primary purpose is for the analysis of building envelopes (e.g., windows, walls, roofs, foundations, etc.). However, the engine can also be used for any heat and moisture transfer problems that require solving fundamental 2-D energy and mass transfer equations. Fluid flow solutions (Navier-Stokes momentum equations) are not included, but the correlations for various convection heat transfer situations are provided, including the translation of complex cavity geometries into those for which correlations are applicable. The calculation engine is written in C++ and includes an API for connecting to third-party tools.

Vidanovic, Dragan [Lawrence Berkeley National Labo↗
DOE OSTI2026

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗
DOE OSTI2024

EchemFEM: A Firedrake-based Python package for electrochemical transport

The transition from fossil fuels to renewable energy has brought about a rapid increase in the availability of clean electricity. However, electricity generated from sources such as wind and solar are limited to intermittent operation due to daily and seasonal variation. One solution is to utilize electrochemical devices in energy storage and electrochemical manufacturing applications, where they can harness surplus energy and decarbonize chemical industries traditionally reliant on petrochemical feedstocks. Managing the growing prevalence of renewable energy underscores the importance of developing and scaling up these technologies, which can in turn facilitate the achievement of carbon emission reduction commitments of companies and developed economies. Likewise, the electrification of transport creates an increasing need for energy-dense electrochemical energy storage devices such as batteries and supercapacitors. Naturally, simulation tools are required to assist in the design of efficient and industrial-scale electrochemical devices.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2019

Combustion Research I.1 Sprays, Flow, Heat Transfer and Turbulent Mixing

LANL is working with the Partnership for Advance Combustion Engines (PACE) a DOE sponsored consortium of national laboratories for the research and development for creating smarter engine design tools for sustainable transportation. Much of that effort is focused on heat, mass and momentum transfer within and engine. LANL’s focus in 2019, in particular under PACE has been spray modeling and engine mixture formulation, development of conjugate heat transfer analysis (heat transfer to the environment surrounding the combustion chamber). Our efforts relay on the Fast, Easy, Accurate and Robust Continuum Engineering (FEARCE) software developed at LANL as the next generation of engine modeling software for turbulent reactive and multiphase flow, a 2019 R&D 100 winner. This work is relevant to the DOE Vehicle Technologies Office efforts at addressing national energy security as well. Less dependence on petroleum products leads to greater energy security. By Environmental Protection Agency standards, some vehicles are now reaching 42–50 mpg mark. These are conventional gasoline engines. With continued investment and research into new technical innovations, the potential exists to save more than 4 million barrels of oil per day or approximately $200 to $400 million per day. This would be a significant decrease in emissions and use of petroleum and a very large stimulus to the U.S. economy. Better understanding of fuel injection and fuel–air mixing, thermodynamic combustion losses, and combustion/emission formation processes enhances our ability to minimize fuel use and unwanted emissions. Helping to accomplish this understanding, the FEARCE modeling program is working to provide state-of-theart capability for accurately simulating combustion processes: to have a predictive methodology in a software helping industry and researchers not only meet national goals on fuel usage and emissions, but global goals. In addition, a predictive, robust, and accurate capability for simulating the engine combustion process helps to minimize time and labor for development of new engine technology.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗