Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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

OpenTurbine [SWR-23-07]

OpenTurbine is an open-source code designed for the simulation of wind turbine structural dynamics. It is a flexible multi-body dynamics solver wherein the wind turbine is modeled as a collection of beam finite elements, rigid bodies, and constraints. OpenTurbine is equipped with an application programming interface primarily designed for coupling to computational fluid dynamics codes, but it is also suited for coupling to aerodynamics solvers based on blade element momentum theory.

Sprague, MichaelA.↗
DOE OSTI2020

An Initial Assessment of the Design Margins of Different ASME Section III, Division 5 Design Rules

This report develops a method for assessing the design margin of the ASME Section III, Division 5 Class A design rules, focusing on a definition of margin as the ratio between the actual, expected component life and the ASME design life. Full inelastic finite element simulations with a complete damage model capturing the available creep, creep-fatigue, fatigue, and tension test failure data produce the expected component lives, which can then be compared to corresponding ASME design calculations. The focus of this particular work is on Alloy 617 and the ASME creep-fatigue design rules, but the approach is general and could be applied to other design limits and to other materials. A margin assessment of a representative Alloy 617 component shows the ASME creep-fatigue design rules applied to this component have a margin of 10 on design life. The report describes several challenges in completing a more comprehensive margin assessment of the ASME rules, in particular challenges in developing probabilistic assessment methods for high temperature structural components. The report details these challenge and discusses the immediate possibilities of the deterministic margin assessment method developed here.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Experimental Study of Mechanical Response of Polyurethane Rubber at Varying Strain Rates for W80-4 System

This report describes an experimental study to determine the mechanical behavior of the polyurethane rubber material that was used in the W80-4 MCC Shock/Breach Phase 1 test series. Compression experiments were conducted on cylindrical specimens over a wide range of loading rates to characterize the material over the range of strain rates that were experienced in the shock/breach testing. Additionally, specimen diameter was varied to determine the effect of confinement on the material response and was found to be significant. This data is used to populate a material model to enable accurate analyses and finite element simulations of the shock/breach test series.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice

The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells’ location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells’ diffusion and immune cells’ influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes’ diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Removal of Non-Specifically Bound Proteins Using Rayleigh Waves Generated on ST-Quartz Substrates

Label-free biosensors are plagued by the issue of non-specific protein binding which negatively affects sensing parameters such as sensitivity, selectivity, and limit-of-detection. In the current work, we explore the possibility of using the Rayleigh waves in ST-Quartz devices to efficiently remove non-specifically bound proteins via acoustic streaming. A coupled-field finite element (FE) fluid structure interaction (FSI) model of a surface acoustic wave (SAW) device based on ST-Quartz substrate in contact with a liquid loading was first used to predict trends in forces related to SAW-induced acoustic streaming. Based on model predictions, it is found that the computed SAW body force is sufficient to overcome adhesive forces between particles and a surface while lift and drag forces prevent reattachment for a range of SAW frequencies. We further performed experiments to validate the model predictions and observe that the excitation of Rayleigh SAWs removed non-specifically bound (NSB) antigens and antibodies from sensing and non-sensing regions, while rinsing and blocking agents were ineffective. An amplified RF signal applied to the device input disrupted the specific interactions between antigens and their capture antibody as well. ST-quartz allows propagation of Rayleigh and leaky SH-SAW waves in orthogonal directions. Thus, the results reported here could allow integration of three important biosensor functions on a single chip, i.e., removal of non-specific binding, mixing, and sensing in the liquid phase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Power Electronics Materials and Bonded Interfaces - Reliability and Lifetime

The emergence of wide-bandgap (WBG) devices in power electronics packages has brought significant attention to the reliability and thermal performance of the bonded interfaces. These interfaces, such as die-attach and substrate-attach, need to withstand the severe stress induced on them by the high-temperature (> 200°C) operation of the WBG devices. Sintered silver and transient liquid phase alloys are key bonded interfaces that can enable the high-temperature operation of WBG devices, however it is essential to evaluate their reliability and study the possible failure mechanisms before they can be used in power electronics packages. Sintered silver is already being adopted in commercial packages as a die-attach, but large-area substrate-attach offers additional reliability challenges. In this project, we investigate the reliability of these interface materials under accelerated thermal cycling and correlate these results with finite element simulations to develop lifetime prediction models. The failure mechanisms of these bonded interfaces are also studied.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Adaptive time scaling for accelerating explicit finite element analysis

A method for accelerating an explicit finite element analysis (FEA) simulation of a modeled system or process includes performing an initial iteration of the FEA simulation according to a baseline time interval via an FEA computing network, and calculating a criteria ratio of a predetermined set of scaling criteria for the modeled system or process. The method includes determining a time-scaling factor using the criteria ratio via the FEA computing network as a function of the criteria ratio, and then applying the time-scaling factor to the baseline time interval to generate a scaled time interval. The scaled time interval accelerates simulation time of the FEA simulation. The method includes performing a subsequent iteration of the explicit FEA simulation at the scaled time interval using the FEA computing network. The process continues for subsequent iterations, with the time-scaling factor adapting with each iteration.

Chen, Jian↗
DOE OSTI2021

Conduction Heat Transfer through Solid in Porous Materials: A Comparative Study by Finite-Element Simulations and Effective Medium Approximations

The development of insulation materials with low effective thermal conductivity is essential for energy savings in various applications, including buildings, food services, pipe insulation, and refrigeration. Such materials can be developed by using micro- or nanoporous structures, as well as low-thermal conductivity gases and/or reduced pressure. A variety of effective medium approximation (EMA) models have been developed to study thermal transport through solid in porous structures. However, in many cases, the impacts of porosity on solid conduction are not well predicted by EMA models because of the assumptions made for simplification. Furthermore, the results can vary by a factor of up to 1.5, depending on the morphology of the pores. Hence, proper guidance is needed to choose the appropriate EMA model for a given morphology. This work presents a finite element method study using COMSOL Multiphysics software for various geometries, including hollow cubes, overlapping hollow spheres (normally stacked and tightly packed), and Voronoi structures, and compares the results to those obtained from various EMA models. Simulation results for the hollow cubes and Voronoi structures matched well with the Maxwell-Eucken and Russell model. The results for hollow spheres closely matched the Bauer model and the Glicksman model with fitting parameters. This work provides guidance on analyzing and designing insulation materials for energy savings in the future.

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 OSTI2020

Volumetric analysis and mesh generation of real and artificial microstructural geometries

Producing a viable finite element mesh of realistic microstructural structural geometry is a critical step in analyzing the thermo-mechanical behavior of complex multi-material composites. Advancements in imaging technology such as micro computed tomography have allowed modelers to access high resolution mesoscale geometries for direct numerical simulation. However, converting from voxel based 3D images to usable finite element meshes has been challenging. A robust method including algorithms and software scripts for generating finite element meshes from 3D imaged microstructures is presented in this paper. It includes a routine for inserting cohesive elements around material interfaces to enable modeling of interface properties including delamination and damage. The algorithms and procedures presented in this method leverage currently available software packages for processing surface based geometry into volume based meshes. In addition to converting real geometry from physical imaging systems, algorithms for producing numerically generated and statistically equivalent microstructural geometry are also included. These artificial microstructures can be a valuable resource for modelers when physical specimens do not exist or are limited in quantity: Method establishes a workflow from voxel data to viable finite element mesh including interface information, and Includes a method for synthetic geometry generation based on metrics from real microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling and analyzing parasitic parameters in high frequency converters

This research focuses on electromagnetic interference (EMI) / electromagnetic compatibility (EMC) design and analysis in power electronics systems. To limit the EMI under the standards, different methods and strategies are investigated. Parasitic parameters of high frequency (HF) transformer are analyzed using a novel analytical method, finite element method (FEM), and experimental measurements for different structures and windings arrangements. Also, the magnetic field, electric field, electric displacement, and electric potential distribution are simulated and analyzed. Moreover, a high voltage system is considered and analyzed to improve the EMC. The EMI propagation paths are analyzed. The EMI noise level of the system is obtained and compared to the IEC61800-3 standard. To improve the EMC, the parasitic parameters of the transformer, as the main path of EMI circulation, are analyzed and optimized to block the propagation. Furthermore, the geometry structure of the HF transformer is optimized to lower the parasitics in the system. Three pareto-optimal techniques are investigated for the optimization. The models and results are verified by 3D-FEM and experimental results for several given scenarios. Furthermore, the EMC modeling and conducted EMI analysis are developed for a system including an AC-DC-DC power supply (rectifier and dual active bridge (DAB) converter). Moreover, the common mode (CM) EMI noise propagation through the system is discussed and the noise sources and effect of components on the noise are analyzed. Additionally, the CM impedance of different parts of the system and the noise levels are discussed. Finally, EMI attenuation techniques were applied to the system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗
DOE OSTI2022

Improved Axisymmetric and High Temperature Material Structural Modeling in MOOSE and NEML

This report describes improvements made to the solid mechanics formulation in the MOOSE open source finite element simulation environment and the open source Nuclear Engineering Material model Library (NEML) for mechanical constitutive models. The focus of these improvements is to improve the usability and performance of simulations involving one or both pieces of software. Specifically, this work completes a new system for solid mechanics simulations in the MOOSE ecosystem providing exact linearizations and optimal (quadratic) convergence, for a variety of coordinate systems and material types, including large deformation simulations. This work then provides users a framework to build highly efficient mechanical simulations of structures or materials or to couple in additional MOOSE physics modules to build complex, scalable multiphysics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗
DOE OSTI2023

Modeling of the Advanced Test Reactor Using OpenMC, Cubit and Griffin

In the pursuit of the ability to perform multiphysics simulations of the Advanced Test Reactor, it is crucial to have a fast and highly accurate deterministic model. To achieve this, a contemporary two-step method is employed. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor multiphysics application based on the Multiphysics Object-Oriented Simulation Environment. To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material IDs are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Preliminary comparisons indicate good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 50 pcm in the two-dimensional geometry configuration. However, in three-dimensional calculations, an unacceptably large error is found in the Griffin solution. Future work is planned to resolve this discrepancy.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Stress intensity factor models using mechanics-guided decomposition and symbolic regression

The finite element method can be used to compute accurate stress intensity factors (SIFs) for cracks with complex geometries and boundary conditions. In contrast, handbook solutions act as surrogate SIF models that provide significantly faster evaluation times. However, the development of conventional surrogate SIF models relies on manual development based on low-order parameterizations. This limits surrogate model accuracy and generalizability. Here, in this paper, we develop a framework for the automated development of mechanics-guided handbook SIF solutions by using interpretable machine learning via genetic programming for symbolic regression (GPSR). Formalizing the mechanics-based approach of Raju and Newman, SIF training data is decomposed into multiple subsets. This decomposition enables parallel GPSR model development of subfunctions, each of which accounts for specific geometrical corrections with respect to a known analytical model. Using this mechanics-based approach with GPSR allows for equations to be learned with improved accuracy and reduced complexity relative to the Raju Newman equations while maintaining the inherent interpretability of mathematical expressions. In this paper, we present equations that match the complexity of the Raju Newman equations while having reduced error, as well as equations with similar errors and reduced complexity.

42 ENGINEERING↗
DOE OSTI2022

A concurrent irradiation-mechanics multiscale coupling model

The vast majority of our current knowledge regarding the basic mechanisms controlling irradiation effect on mechanical properties is based almost entirely on the results of post-irradiation experiments or theoretical models. However, the concurrent effects of irradiation, mechanical stress, and thermal damage on the failure phenomena of materials and components remain largely unexplored due to its internal multiscale-multiphysics coupling nature. Here, we present here a concurrent irradiation-mechanics multiscale coupling model. The concurrent evolutions of nanoscale irradiation defect clusters, microscale dislocation configurations, and mechanical responses are well captured based on coupling cluster dynamics, discrete dislocation dynamics, and the finite element methods using an effective time marching scheme. Model predictions of defect densities and size are in general agreement with experimental observations. Irradiation hardening is shown to take place also in samples undergoing concurrent irradiation-mechanical loading, similar to samples tested post-irradiation. However, the occurrence of plastic flow localization and dislocation channel formation is not accompanied with apparent yield drop (softening) under concurrent irradiation-mechanical loading conditions, which is different from the post-irradiation case.

42 ENGINEERING↗
DOE OSTI2024

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗