Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element simulation”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Characterization of Fe-6Si Soft Magnetic Alloy Produced by Laser-Directed Energy Deposition Additive Manufacturing

Abstract Commercial electrical steels, Fe-Si alloys with < 4 wt.% Si, are inexpensive and efficient materials for electrical power conversion. Further efficiency improvements require increasing the silicon concentration to 6 wt.%, at which point the material becomes brittle and difficult to form by conventional rolling and sheet fabrication methods. Additive manufacturing stands to overcome challenges with commercial manufacturing techniques by leveraging near-net-shape fabrication. The wide array of process conditions provides additive manufacturing with increased flexibility, enabling control over the microstructure and mechanical properties. This work explores the microstructures and magnetic properties of ring-shaped Fe-Si alloys produced using concentric and cross-hatch tool paths on a laser-directed energy deposition additive manufacturing system. Concentric-built samples exhibit elongated grain structures while cross-hatch-built samples comprise lower aspect ratio grain structures. Thermal finite element analysis simulations model the stress conditions produced by the different scan path geometries. Microhardness measurements probe the mechanical properties as a function of anneal temperature, providing a qualitative understanding of the intergranular defect density. Soft magnetic properties measured under quasistatic and AC conditions show frequency- and microstructure-dependent coercivity and permeability. Finally, analysis of the core loss quantifies how the build strategies and thermal treatments influence efficiency in electrical power conversion applications. Understanding the influences of scan path geometry and thermal treatment provides a pathway towards application of additively manufactured soft magnetic materials.

Adamczyk, Jesse M.↗
DOE OSTI2023

Contact pressure analysis to allow improved design for the clearing plate in a biomass comminution system

For biomass processing systems, feedstocks with small but equal sizes can improve flowability and thermochemical outputs. The Crumbler rotary shear system was designed to achieve such ideal feedstocks to reduce costs and energy throughout the biomass supply chain. Improving the wear resistance of major components in this system is important to decrease the cost. Accelerating the feedstock flow through the rotary shear machine plays an equivalent role of improving its efficiency. An analytical analysis and a finite-element analysis of the stationary clearing plate are presented in this work. Here, these tools are used to investigate the effect of contact surface curvature on the impact pressure. The optimized new clearing plate design is expected to improve the component’s lifetime and wood chips flowability. The field tests confirm that the trend of the new clearing plate design guided by FEA. The optimization method, model verification, and validation experience applied in this work can also be applied to other static components with similar simple contact wear problems.

09 BIOMASS FUELS↗
DOE OSTI2018

Investigation on Vibration Response of Aluminum Foam Beams Using Speckle Interferometry

This paper firstly presents a brief review on the representations of speckle fringe pattern in vibration measurement using time-averaged electronic speckle pattern interference (ESPI) method. Based on laser phase noise, a new representation of speckle fringe intensity generated by real-time subtraction ESPI method for vibrating measurement is then proposed. In the experimental performance, the real-time subtraction ESPI method is employed to inspect the vibrating response of a closed-cell aluminum foam cantilever beam. In contract with the results of the finite elemental method simulation, the vibration mode shapes obtained by the ESPI method are well agreed with numerical prediction. The optical vibration analysis is also carried out to determine the effective Young’s modulus of aluminum foam, and the results verify the validation of the ESPI method for investigation on mechanical properties of metal foam materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Statistical analysis of spatial distribution of pores in metal additive manufacturing

We report mechanical properties of additively manufactured parts are sensitive to the presence of pores form during the manufacturing process. The impact of pores on the mechanical performance has been investigated extensively with respect to different parameters, such as pore volume fraction, shape, and size. However, statistical investigations focusing on the relationships between the spatial distribution of pores and process parameters; and consequently, the performance of the manufactured parts are scattered and limited. Also, these sparse investigations usually suffer from an ambiguous definition of terminologies. For instance, the required criteria to consider a point pattern as complete spatial randomness (CSR) are generally not clarified. Moreover, no numerical formalism is yet developed to show how much the observed spatial results are statistically significant. To address these shortcomings, the statistical definition of CSR in a point pattern and the procedure to quantitatively determine the deviation of a pattern from CSR were discussed. The explained statistical approach was used to investigate the effect of scanning speed parameter on the spatial distribution of gas pores in laser powder bed fusion manufactured stainless steel parts. Furthermore, and to highlight the impact of the spatial distribution of pores on mechanical properties, fatigue performances of parts with clustered and randomly distributed pores were simulated by finite element analysis. It is shown that by reducing the scanning speed, the spatial distribution of gas pores deviated more from CSR, and correspondingly fatigue performance deteriorated.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Reduction of baseplate distortion during directed energy deposition using compliant features

Distortion in additive manufacturing (AM) remains a barrier to its adoption in precision industries. Baseplate warpage is one such issue which compromises the feasibility of post-process precision machining. The restriction to thermal contraction of the deposited part, imposed by the baseplate, generates bending moments, that in turn causes warpage. A novel distortion mitigation strategy using baseplates with integrated compliant features, which enables thermal contraction of the build is presented. Six unique design concepts are evaluated, including a solid reference, using two deposition geometries. A laser, hot-wire, directed energy deposition (DED) process is used to deposit a symmetric cylindrical part (C-part) and a T-shaped asymmetric part (T-part). Flatness deviation of baseplates is measured using structured light 3D scanning. Measurements reveal a reduction in net flatness deviation of 58.8% for the C-part and 40.9% for the T-part, compared to the solid reference. While most designs yielded reductions exceeding 30% and 20% for the C- and T-parts, respectively, one configuration resulted in increased deviation. Finite element (FE) simulations are performed to elucidate the underlying mechanisms affecting distortion of compliant baseplates during DED. Despite variations between predictions and measurements, agreement in the general trend is observed. It also revealed that initial flatness errors in baseplates significantly affect its deviation during deposition. Predictions indicate that compliant features significantly affect the thermal distribution as well as the evolution of flatness deviation in the baseplate during deposition. Notably, one design exhibited a reduction in distortion during cooling, following its initial increase during deposition. FE predictions show a maximum reduction of 60.9% and 38.8% in net flatness deviation for the C- and T-parts, respectively. The performance of compliant baseplates is found to be governed by both its the thermal and mechanical characteristics, which are crucial factors to be considered during design.

Mathews, Ritin [ORNL] (ORCID:0000000301440828)↗
DOE OSTI2025

Hydrogen charging and desorption from microstructural viewpoint: A method for deconvoluting hydrogen desorption spectra and unveiling the hydrogen-microstructure interaction

Understanding the interaction of hydrogen with microstructural features in metallic materials is crucial for designing hydrogen-resistant alloys. Although thermal desorption spectroscopy (TDS) is widely used for investigating the hydrogen binding behavior of various microstructural features, its application to face-centered cubic (fcc) metals and alloys that exhibit low hydrogen diffusivity is limited due to the lumped TDS desorption signals. This paper shows that, by coupling a Sofronis–McMeeking type hydrogen transport model with a microstructure-informed finite-element model, TDS data can be deconvoluted to reveal the underlying adsorption–diffusion–desorption processes, hydrogen diffusivity, and trap-binding energies. In conclusion, the austenitic steel SS316L in solution-annealed condition is used as a demonstration material, and we focused on investigating the interaction of deuterium (hydrogen isotope) with grain boundaries, which is difficult to investigate from experiments alone but critical for design of alloys for hydrogen infrastructure.

Finite element simulation, Polycrystalline Microst↗
DOE OSTI2026

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗
DOE OSTI2025

Co-training of multiple neural networks for simultaneous optimization and training of physics-informed neural networks for composite curing

This paper introduces a Physics-Informed Neural Network (PINN) technique that co-trains neural networks (NNs) that represent each function in a system of equations to simultaneously solve equations representing an out-of-autoclave (OOA) cure process while conducting optimization in adherence to process requirements. Specifically, this co-training approach benefits from using NNs to represent OOA inputs (air temperature profile) and outputs (part and tool temperature profiles and degree of cure). Production requirements can then be levied on the inputs, such as maximum air temperature and minimum cure cycle, and simultaneously on the outputs, such as degree of cure, maximum part temperature, and part temperature rate limits. The technique is validated with finite element (FE) simulations and physical experiments for curing a Toray T830H-6 K/3900-2D composite panel. Furthermore, this novel approach efficiently models and optimizes the OOA cure process.

Composite curing↗
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

A Physically-Based Model for Thermo-Oxidative and Hydrolytic Aging of Elastomers

A computationally efficient model is proposed to capture the loss of mechanical performance due to chemical aging that are formed as the competition of chain scission and cross-link formation/dissolution, such as thermo-oxidative aging or hydrolytic aging. The model should be considered an extension of our recent models [1, 2, 3] which further simplifies the matrix behavior based on the assumption of independence of environmental and mechanical damage. The model uses this assumption to reduce the necessary material parameters needed to model constitutive and inelastic behavior of elastomers during aging. To this end, the model can provide accurate predictions of the material performance with the significantly fewer number of fitting parameters. The model is relevant for all decay mechanisms formed by the occurrence of two simultaneous micro-mechanisms; (i) formation/reduction of the cross-links, and (ii) chain scission, both of which are present in thermo-oxidation and hydrolytic aging. Assuming the alteration of the chain density along the aging trajectory is identical to the peroxide cross-link density for thermo-oxidation, and the change of the average molecular weight for hydrolysis, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time and aging temperature. Next, the modified network alteration model is formulated for implementation into Finite Element (FE) simulations. The model is built on the presumption of homogeneous and consistent oxygen/water absorption and thus is mainly relevant for relatively thin samples exposed to environmental loads for a long time. The proposed model includes only six physically inspired material parameters. Therefore, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is the right choice for advanced implementations in FE programs.

42 ENGINEERING↗
DOE OSTI2019

Predicting forming limit diagrams for magnesium alloys using crystal plasticity finite elements

The crystal plasticity finite element (CPFE) method was used to predict forming limit strains of hexagonal close-packed (HCP) polycrystalline magnesium alloy sheets, namely AZ31 and ZE10, under an isothermal temperature condition. The strain rate-dependent uniaxial tensile test data along various loading directions and an initial texture were used to determine the constitutive parameters of the crystal plasticity model. A hybrid representative volume element approach, which combines the CPFE and the Marciniak–Kuczynski model, was developed to obtain the forming limit diagram of the magnesium alloys. The predicted forming limits were in excellent agreement with the Nakazima test results. Moreover, the microscopic responses such as slip/twin activation and deformation texture changes during various loading paths from uniaxial tension to the balanced biaxial tension were comprehensively analyzed and discussed from the CPFE results to understand the underlying micro-mechanism for the macro-mechanical responses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling inter- and intra-granular dislocation transport using crystal plasticity

Here, this work presents the development of a crystal plasticity material model that incorporates both dislocation transport within grains and dislocation transfer across grain boundaries. This model has been implemented in the open-source finite element code MOOSE. In addition, a novel geometry-based criterion is developed to determine the direction of dislocation transfer across grain boundaries. The transfer criterion incorporates the geometric features of the grain boundary, such as the grain boundary plane normal, and its misorientation, which is accounted for through the orientation of the incoming and outgoing slip systems. The model is tested with several cases, including a copper single crystal, bi-crystal, and polycrystal. The development of the transfer criterion, implementation of the model, and its application to these test cases are discussed in detail.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Scalable preconditioning for the stabilized contact mechanics problem

We present a family of preconditioning strategies for the contact problem in fractured and faulted porous media. We combine low-order continuous finite elements to simulate the bulk deformation with piecewise constant Lagrange multipliers to impose the frictional contact constraints. This formulation is not uniformly inf-sup stable and requires stabilization. We improve previous work by Franceschini et al. (2020) by introducing a novel jump stabilization technique that requires only local geometrical and mechanical properties. We then design scalable preconditioning strategies that take advantage of the block structure of the Jacobian matrix using a physics-based partitioning of the unknowns by field type, namely displacement and Lagrange multipliers. The key to the success of the proposed preconditioners is a pseudo-Schur complement obtained by eliminating the Lagrange multiplier degrees of freedom, which can then be efficiently solved using an optimal multigrid method. Numerical results, including complex real-world problems, are presented to illustrate theoretical properties, scalability and robustness of the preconditioner. A comparison with other approaches available in the literature is also provided.

58 GEOSCIENCES↗
DOE OSTI2023

Corrosion of Al-Fe self-pierce riveting joints with multiphysics-based modeling and experiments

Self-piercing riveting (SPR) is an extensively used joining technique to assemble dissimilar materials. However, this joining of the dissimilar materials can generate galvanic/crevice corrosion, which can drastically impact the mechanical properties and the service life of the joint. In this study, a multiphysics-based hybrid modeling approach is developed for the galvanic corrosion of the Al-Fe SPR joints, which can consider the corrosion initiation and corrosion evaluation jointly. Experimental studies are first performed to extract information regarding the corrosion initiation sites, corrosion evolution and overall corrosion induced material loss in SPR joints. The initiation information is then passed onto the multiphysics FE model thus making it a hybrid model. This model can help understand the influences of metal microstructure on the corrosion propagation while the morphology changes can also be analyzed. Using the developed hybrid modeling approach, thorough parametric studies can be performed to explore the coupled impacts of multiple corrosion factors on the corrosion behavior of the joints. The developed hybrid model is validated on the prediction of galvanic corrosion for the SPR joints by comparing with experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Fabrication and characterization of powder metallurgy tantalum components prepared by high compaction pressure technique

The present study has investigated the consolidation behaviors of tantalum powders during compaction and sintering, and the characteristics of sintered components. For die compaction, the densification behaviors of the powders are simulated by finite element analyses based on the yield function proposed by Shima and Oyane. Accordingly, the green density distribution for coarser particles is predicted to be more uniform because they exhibits higher initial relative tap density owing to lower interparticle friction. It is also found that cold isostatic pressing is capable of producing higher dense compacts compared to the die pressing. However, unlike the compaction behavior, the sintered density of smaller particles is found to be higher than those of coarser ones owing to their higher specific surface area. The maximum sintered density was found to be 0.96 of theoretical density where smaller particles were pressed isostatically at 400 MPa followed by sintering at 2000 °C. Moreover, the effects of processing conditions on grain size and texture were also investigated. The average grain size of the sintered specimen is 30.29 μm and its texture is less than 2 times random intensity. Consequently, it is concluded that the higher pressure compaction technique is beneficial to produce high dense and texture-free tantalum components compared to hot pressing and spark plasma sintering. - Highlights: • Higher Ta density is obtained from higher pressure and sintering temperature. • High compaction method enables P/M Ta to achieve the density of 16.00 g·cm{sup −3}. • A P/M Ta component with fine microstructure and random orientation is developed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

An isothermal forming process with multi-stage variable speed for magnesium component assisted by sensitivity analysis

We have developed a multi-stage variable speed (MSV) isothermal forming process based on strain rate sensitivity (SRS) analysis, and successfully produced a large size of rib-web component of magnesium alloy from a cylindrical billet. The obtained rib-web component has been processed into parts for mechanical properties examination and the examination results met the standard of being safely used in aircrafts. The process parameters of the forming technology and microstructure change were explicitly linked by the SRS coefficient. The finite element (FE) simulation and SRS distribution were implemented to identify the relationship among process parameters, shaping characteristics, and microstructure development. The MSV speed process was employed for reducing the forming load, refining grain sizes and increasing the strength of the component. With the selected process parameters the maximum forming load was reduced to be less than 20 MN, which is 55.3% lower than that used in the same process with the initial constant speed unchanged. Fine dynamic recrystallization (DRX) were observed from all parts of the forged cover plate including the convex, ribs, and web region. It was also found that higher strain resulted in more uniform distribution of fine ß-Mg17Al12 precipitates and larger DRX volume fraction. Importantly the developed isothermal forming process with severe plastic deformation can be applied for obtaining grain refinement so improving the mechanical properties of Mg products.

Magnesium alloy, Isothermal forming, Multi-stage v↗
DOE OSTI2022

The effects of microstructure on deformation twinning in Mg WE43

The interplay between microstructure and deformation twinning in a WE43-T6 Mg alloy under uniaxial compression was investigated using a combination of scanning electron microscopy with digital image correlation (SEM-DIC) and crystal plasticity finite element (CPFE) simulation. To improve understanding of the statistical characteristics of deformation twin formation, microstructural effects were characterized in over 1000 grains through metrics including the nominal Schmid Factor, grain size, geometric compatibility factor (m'), residual Burgers vector, and the strain accommodated by neighboring grains. There was a strong correlation between the nominal Schmid Factor and both twin activation and variant selection, but this was not fully deterministic. Deformation twinning also exhibited a strong dependence on existing slip and twinning in the neighboring grain and on the m' value. Within the range of grain sizes present in this condition, grain size was determined to have minimal effect on deformation twinning. Finally, statistical analysis of CPFE simulations was used to further investigate microstructural effects on twinning, and qualitatively captured the effect of the nominal Schmid Factor.

36 MATERIALS SCIENCE↗