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

A generalized 3D elastic model for nanoscale, self-assembled oxide-metal thin films with pillar-in-matrix configurations

In recent years, functional oxide-metal based vertically aligned nanocomposite (VAN) thin films have gained interest due to their intriguing physical properties and multifunctionalities stemming from the complex interactions between the two phases in the film and the substrate. In this work, we develop a model for studying the energetics of these thin film systems, including the effects of both lattice mismatch and capillary forces due to interface curvature. Each phase is incorporated into the model using a phase indicator function, and we introduce the capillary forces as body forces using a vector density representation of the interface. The model is implemented using the finite element method to study the deformation of the thin film which is composed of Au nanopillars embedded in a La 0.7 Sr 0.3 MnO 3 (LSMO) matrix on an SrTiO 3 (STO) substrate. Further, the results suggest that the total energy is lowest for random configurations of pillars compared to ordered square and hexagonal lattice configurations, consistent with the random distribution of pillars found in experiments. Furthermore, we find that the interfacial energy dominates the total energy of each configuration, suggesting that interfacial energy in the system is an important design parameter for nanocomposite growth, along with the lattice mismatch.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of the plastic deformation microstructure in metals on the kinetics of recrystallization: A phase-field study

The sensitivity of recrystallization kinetics in metals to the heterogeneity of microstructure and deformation history is a widely accepted experimental fact. However, most of the available recrystallization models employ either a mean field approach or use grain-averaged parameters, and thus neglecting the mesoscopic heterogeneity induced by prior deformation. Here in the present study, we investigate the impact of deformation-induced dislocation (subgrain) structure on the kinetics of recrystallization in metals using the phase-field approach. The primary focus here is upon the role of dislocation cell boundaries. The free energy formulation of the phase-field model accounts for the heterogeneity of the microstructure by assigning localized energy to the resulting dislocation microstructure realizations generated from experimental data. These microstructure realizations are created using the universal scaling laws for the spacing and the misorientation angles of both the geometrically necessary and incidental dislocation boundaries. The resulting free energy is used into an Allen-Cahn based model of recrystallization kinetics, which are solved using the finite element method. The solutions thus obtained shed light on the critical role of the spatial heterogeneity of deformation in the non-smooth growth of recrystallization nuclei and on the final grain structure. The results showed that, in agreement with experiment, the morphology of recrystallization front exhibits protrusions and retrusions. By resolving the subgrain structure, the presented algorithm paves the way for developing predictive kinetic models that fully account for the deformed state of recrystallizing metals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A model for oblique impacts on material surfaces

Many practical situations of material damage, wear, and erosion involve collisions between small particles and surfaces at inclined angles. While there are many well-validated models of normal incidence impact situations, elastic-plastic models for oblique incidence impact events are lacking. Here the finite element method is used to predict the normal and tangential coefficient of restitution in oblique impacts for hard, elastic spheres impacting an elastic-perfectly plastic material surface. The proposed model covers various impact angles ranging from 0° to 45°, within a limiting impact velocity below which the effects of heating are negligible. The normal coefficient of restitution follows power-laws with respect to normalized values of the impact velocity. Interestingly, the tangential coefficient of restitution follows a linear relationship with impact velocity. Together, these results provide a semi-empirical set of equations predicting oblique impact rebounds (both velocity and trajectory) for a wide range of conditions and material properties, with which experimental results can be rapidly interpreted. As a result, Laser-Induced Particle Impact Test (LIPIT) data are also presented for aluminum particles impacting aluminum substrates, at impact angles of 25° and 40°; the results compare favorably with the model and validate the general use of such models for the analysis of experimental data.

Coefficient of restitution↗
DOE OSTI2026

Transient Multiphysics Simulations with Pin Power Reconstruction in the Griffin Reactor Physics Code

This work introduces the pin power reconstruction capability available in the Griffin reactor physics code. This capability is implemented in an unstructured mesh framework, and the methods introduced are applied to the 2D SIMBA reactor core, which has assemblies and pins arranged in a hexagonal lattice. Since this reactor has a non-Cartesian geometry and also operates in the thermal spectrum, a general approach to pin power reconstruction is adopted, where SPH-based equivalence is leveraged to preserve assembly-wise reaction rates, while computing full-core form functions to preserve pin-wise fission production rates within the fuel pins of the reactor core. In a 2D microreactor benchmark problem, this pin power reconstruction approach was shown to reproduce pin powers compared to the Serpent2 Monte Carlo code for fixed temperature conditions and control drum rotation angles, yielding a core-wide RMS error level of 0.6\% and a maximum absolute pin error of 2.3\%. In addition, a tabulated library of multigroup cross sections, SPH factors, and form functions was generated to demonstrate the applicability of pin power reconstruction to a thermal feedback problem. Finally, a control drum transient was successfully simulated, showcasing the application of pin power reconstruction in a transient multiphysics feedback problem.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Simultaneous shape and topology optimization of inflatable soft robots

Simultaneous shape and topology optimization is used to design pressure-activated inflatable soft robots. The pressure loaded boundary is meshed conformingly and shape optimized, while the morphology of the robot is topology optimized. The design objective is to exert maximum force on an object, i.e. to produce soft “grippers”. The robot’s motion is modeled using nearly incompressible finite deformation hyperelasticity. To ensure stability of the robot, the buckling load factors obtained via linearized buckling analyses are constrained. The finite element method is used to evaluate the optimization cost and constraint functions and the adjoint method is employed to compute their sensitivities. The numerical examples produce pressure-driven soft robots with varying complexity. We also compare our simultaneous optimization results to those obtained via sequential topology and then shape optimization.

42 ENGINEERING↗
DOE OSTI2025

Influence of sizing concentration on strength, stiffness, and porosity in textile grade carbon fiber (TCF)-Epoxy composites: Revealing inverse trends

The effect of fiber sizing (i.e., surface treatment) concentration (0 %, 1.36 %, 1.52 %, 1.94 %, and 2.13 %) on the mechanical properties (tensile, flexural, interlaminar shear strength (ILSS), and low velocity impact) of textile grade carbon fiber (TCF)-epoxy composite is examined. An inverse relationship between the strength and stiffness of the composite is observed with increased sizing concentration. The root mean square (RMS) roughness of the fiber surface increased from 17.8 nm (unsized) to 22.7 nm with 2.13 % sizing concentration. It was found that the tensile strength increased by 131 % from 221.4 ± 18.5 MPa (unsized) to 510.8 ± 28.05 MPa (for 1.36 % sizing) and further by 155 %–563.7 ± 14.95 MPa at 2.13 % sizing. On the contrary, the stiffness is initially increased by 126 % from 33.52 ± 7.80 GPa (unsized) to 75.9 ± 3.21 GPa (for 1.36 % sizing) but reduced with further increase in the sizing concentration. A single fiber pull-out test is simulated using the finite element method to validate the reverse trend in strength and stiffness. The varying sizing concentration is simulated by introducing an interface of varying thickness between fibers and matrix. Simulation results confirm that a thicker interface, corresponding to a higher sizing concentration, decreases interfacial shear stress, enhancing material strength while decreasing stiffness. The reverse trend in strength and stiffness with the sizing concentration aligns with experimental observations. In conclusion, the present study emphasizes the importance of sizing concentration for mechanical properties and provide a design criterion for customized high-strength and high-stiffness applications.

Porosity↗
DOE OSTI2023

Simulation toolkit for digital material characterization of large image-based microstructures

In this paper, an efficient image-based simulation toolkit for material characterization is presented, which is scalable to work from personal computers to workstations. The effective thermal conductivity, elasticity, and permeability are evaluated employing a computational homogenization framework based on the Finite Element Method (FEM). Two complementary open-source packages are presented: one developed in Python, which can convert digital images into voxel meshes (pyTomoviewer); the other developed in Julia, that can run numerical simulations to compute effective material properties (chpack). Also, a CUDA C version of chpack is provided (chfem_gpu). They were designed to deal with large multi-phase models, so strategies were devised to minimize their memory footprint, while avoiding a high toll on execution time. The voxel-based approach significantly simplifies the FEM meshes and allows efficient matrix-free implementations. In that sense, to handle large linear systems of equations, the element-by-element (EBE) technique is adopted, in conjunction with a low-memory implementation of the Preconditioned Conjugate Gradient (PCG) method. Finally, the code was thoroughly tested on an artificial geometry made of a square array of cylinders, for which analytical solutions exist, as well as on a real micro-tomographic reconstruction of FiberForm TM , a carbon preform commonly used in thermal protection systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Cost-Performance Trade-off in thermoelectric air conditioning system with graded and constant material properties

Thermoelectric (TE) air cooling is a solid-state technology that has the potential to replace conventional vapor compression-based air conditioning. In this paper, we present a detailed system-level modeling for thermoelectric air conditioning system with position-dependent (graded) and constant material properties. Strategies for design optimization of the system are provided in terms of cost-performance trade-off. Realistic convection heat transfer at both sides of the system are taken into account in our modeling. Here, the effects of convection heat transfer coefficients, air flowrate, and thermoelectric material properties are investigated with varying key parameters such as TE leg thickness, module fill factor, and input current. Both constant material properties and graded properties are considered for the TE materials, and they are compared in terms of the degree of cooling, coefficient of performance (COP), and power consumption. For graded materials, we employ one-dimensional finite element methods to solve the coupled electrical and thermal current equations with arbitrary profile of material properties varying with position along the TE legs. We find that graded materials can enhance the degree of cooling, but only at the expense of COP, compared to the case of constant property materials. With constant material properties of ZT = 1 and relatively low electric current, the power consumption of TE coolers can be lower than those of conventional air conditioners at an equivalent cooling capacity. Considering additional advantages such as demand-flexible operation, low noise, and high scalability, thermoelectric cooling could be a competitive technology for future air conditioning applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Environmental influence on cracking and debonding of electrically conductive adhesives

Electrically conductive adhesives (ECAs) are starting to replace metallic solders in recent designs of photovoltaic (PV) modules. This transition represents a significant material change, and a proper understanding of the durability and reliability of the new interconnect needs to be established. This paper presents our continued work on developing a degradation model for ECA interconnects in PV modules. Here, we characterize the fracture mechanics properties of an epoxy-based ECA, for both critical and subcritical, mode I and mode II loading conditions. Emphasis is on the influence of different environmental conditions such as temperature and humidity. We use the Finite Element Method to account for residual stresses, induced by temperature changes and moisture absorption, and correct the apparent fracture toughness. We found that high moisture levels not only can weaken the fracture resistance of the ECA interconnect, but can also promote subcritical debonding at significantly lower driving forces than in dry environments.

42 ENGINEERING↗
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 OSTI2023

Modeling the 4D discharge of lithium-ion batteries with a multiscale time-dependent deep learning framework

The lithium-ion battery (LIB) field is moving towards the direction of investigating spatially resolved physical phenomena in the 3D porous microstructure of electrodes. These pore-scale simulations give new insights into the local dynamics of lithiation/de-lithiation and charge transport, Nevertheless, the computational time of these simulations limits the integration of these models in optimization workflows of cycling conditions or electrode manufacturing processes. Machine learning models present a way of assessing in real-time the performance of materials. While several successful techniques for replicating simulations with machine learning have been proposed, this case study presents a more demanding problem, due to the necessity of understanding the behavior of heterogeneous 3D local data, as it evolves in time: this poses both a scientific and a technical challenge. To this end, we propose an autoregressive multiscale convolutional neural network model to predict relevant quantities at the pore-scale in the solid phase: the lithium concentration (in the active material) and potential (in the active material and carbon binder). Here, these are ultimately used to reconstruct the battery discharge curve. 3D images of the electrode microstructures are the input to the network, trained with a dataset of finite element method simulations to predict the discharge behavior of the cathode side in lithium ion batteries. We propose this machine learning model as a proof-of-concept of the applicability of multiscale networks for time-dependent physics problems. The trained model exhibits very high accuracy (with errors lower than 2 %) in forecasting the discharge behavior of new unseen cathodes.

25 ENERGY STORAGE↗
DOE OSTI2022

Microstructure-Sensitive modeling of surface roughness and notch effects on extreme value fatigue response

Fatigue life knockdown due to surface roughness induces local stress concentrations. Physical fatigue experiments are crucial to qualify materials for use but are limited in their capability to examine the mechanisms underlying the fatigue crack formation driving force due to these stress concentrations. Microstructure-sensitive computational models can aid in this regard and are employed in this work to examine the effects of surface roughness in Al 7075-T6. Digital microstructure models are subjected to crystal plasticity finite element method (CPFEM) simulations to examine the combined effects of microstructure and grain scale asperities on the driving force for fatigue crack formation. Following elastic–plastic shakedown, mesoscale volume-averaged fatigue indicator parameters (FIPs) are computed within fatigue damage process zones of grains. An increase in the intensity of realistic surface roughness profiles corresponds accordingly to larger FIPs but it is difficult to precisely interpret the mechanisms and effects of this intensification. We thus parametrically investigate surface asperities that couple with microstructure and lack of constraint on slip at the free surface to describe microstructure-sensitive surface roughness knockdown effects on fatigue resistance. The depth of a single valley (i.e., notch) due to surface roughness is found to be more detrimental to fatigue resistance than the radius (i.e., notch width, not to be confused with notch acuity). Finally, we also find that different crystallographic textures, grain morphologies, and number of sampled grains have a limited effect on the depth of influence of these notches.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗
DOE OSTI2020

Composite Structured Surfaces for Durable Dropwise Condensation

Dropwise condensation of steam on hydrophobic substrates has a 10X higher heat transfer coefficient compared to filmwise condensation. To promote dropwise condensation, low surface energy hydrophobic coatings (polymers) are typically utilized. The low intrinsic thermal conductivity (k < 1 W/(m · K)) of polymers, coupled with high heat transfer coefficient of dropwise condensation (100 kW/(m 2 · K)), necessitates that the coating be thin (< 1µm) in order to avoid reducing the overall heat exchanger conductance. However, thin polymeric films easily degrade. The two opposing requirements result in the need for optimization between the durability (thick coating) and the heat transfer (thin coating). To enable high thermal conductivity in thicker coatings, we develop metal-polymer structured surfaces. By using porous structures as inter-connected heat-conducting backbones that are filled with hydrophobic materials, we enable tuning of the coating effective thermal conductivity and surface energy. Three metal structures were studied; micro/nanowires, inverse opals, and sintered spheres. Heat transfer performance was calculated using three-dimensional finite element method simulations with two distinct boundary conditions; convection at the walls and isothermal walls. Interestingly, the overall conductance shows up to 40% difference depending on the boundary condition used in calculating the composite coating effective thermal conductivity. We use our model to predict the heat transfer performance as a function of metal fraction by volume and by surface area for condensation. By coupling our thermal simulations with a previously verified analytical model for predicting wetting behavior on heterogeneous surfaces, we propose a regime map to predict dropwise-to-filmwise transition. Furthermore, our work not only forms a starting point for the development of durable dropwise condensing surfaces, it identifies important considerations needed for computing effective thermal conductivity of composites.

42 ENGINEERING↗
DOE OSTI2020

Modeling The Effects of Loading Scenario and Thermal Expansion Coefficient on Potential Failure of Cryo-compressed Hydrogen Vessels

A multiscale thermomechanical model for a simplified Type-3 cryogenic compressed hydrogen (H2) storage vessel is developed in this paper. The model accounts for the temperature-dependent elastic-plastic behavior of the vessel carbon/epoxy composite overwrap and aluminum alloy liner. The homogenized thermo-elastic-plastic behavior for the individual laminas of the vessel layup is obtained by an incremental Eshelby-Mori-Tanka approach associated with a micromechanical failure criterion to predict lamina failure while a standard elastic-plastic constitutive model is used to describe the behavior of a typical aluminum alloy assumed for the liner. The vessel response to external loadings is achieved by a finite element method. Four loading scenarios representing four thermomechanical cycles applied to the vessel are analyzed to evaluate constituent and lamina stresses as well as the associate failure criterion during the cycle according to these scenarios. The model can provide helpful guidance to mitigate thermal stresses by an adequate selection of loading scenario, optimizing the layup and by tailoring thermomechanical properties of the resin matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗