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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
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 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 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↗
DOE OSTI2021

Tuning Electrogenerated Chemiluminescence Intensity Enhancement Using Hexagonal Lattice Arrays of Gold Nanodisks

Electrogenerated chemiluminescence (ECL) microscopy shows promise as a technique for mapping chemical reactions on single nanoparticles. The technique’s spatial resolution is limited by the quantum yield of the emission and the diffusive nature of the ECL process. To improve signal intensity, ECL dyes have been coupled with plasmonic nanoparticles, which act as nanoantennas. Here, we characterize the optical properties of hexagonal arrays of gold nanodisks and how they impact the enhancement of ECL from the coreaction of tris(2,2'-bipyridyl)dichlororuthenium(II) hexahydrate and tripropylamine. Here, we find that varying the lattice spacing results in a 23-fold enhancement of ECL intensity because of increased dye-array near-field coupling as modeled using finite element method simulations.

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

Curvilinear Magnonic Crystal Based on 3D Hierarchical Nanotemplates

Curvilinear magnetic nanostructures enable control of magnetization dynamics through geometry-induced anisotropy and chiral interactions, as well as magnetic field modulation. In this work, we report a curvilinear magnonic crystal based on large-area square arrays of truncated nanospikes fabricated by conformal coating of 3D hierarchical templates with permalloy thin films. Brillouin light scattering spectroscopy reveals an anisotropic band structure with multiple dispersive and folded Bloch-type dispersive spin-wave modes as well as nondispersive modes exhibiting direction-dependent frequency shifts and intensity asymmetries along lattice principal axes. Finite element micromagnetic simulations indicate that curvature-induced variations of the demagnetizing field govern the magnonic response, enabling the identification of modes propagating in nanochannels and others localized on nanospike apexes or along the ridges connecting adjacent nanospikes. The combination of geometric curvature and optical probing asymmetry produces directional dependence of magnonic bands, establishing 3D hierarchical templates as a versatile platform for curvature-engineered magnonics.

3D hierarchical templates↗
DOE OSTI2021

Trap-Assisted Dopant Compensation Prevents Shunting in Poly-Si Passivating Interdigitated Back Contact Silicon Solar Cells

Using a trap-assisted compensation model, we explain why polycrystalline Si (poly-Si) passivating contacts are able to achieve low leakage current between the doped fingers of interdigitated back contact (IBC) monocrystalline Si solar cells despite mixing of boron and phosphorus dopants in the isolation region. The fill factor of IBC solar cells is strongly affected by the electrical isolation region between n- and p-type fingers, as this region is critical in minimizing shunting losses. During fabrication of monocrystalline Si solar cells with poly-Si passivating contacts, the intrinsic poly-Si isolation region inevitably gets contaminated with both p- and n-type dopants. Using dopant profiles measured with time-of-flight secondary ion mass spectrometry and scanning probe measurements of the isolation region, we demonstrate that despite the dopant spreading during cell processing, a well-compensated region between the doped fingers exists that prevents shunting. The trap-assisted dopant compensation mechanism significantly widens the compensated region to tens of microns, where the residual dopant densities are below the trap density. This enables a high-resistivity region, resulting in low shunt current. Using one-dimensional (1-D) finite element diode simulations, we identify the design parameters and experimental conditions under which a sufficiently resistive region can form. Furthermore, our measurements of 2-D local resistivity and work function maps across the isolation region using scanning spreading resistance microscopy and Kelvin probe force microscopy demonstrate the existence of a highly resistive, wide compensated region and confirm our proposed compensation mechanism. For our structures, this region is ~25 µm in width within a ~150 µm wide finger isolation region with nearly 3 orders of magnitude higher resistivity than the regions dominated by a single type of dopant.

14 SOLAR ENERGY↗
DOE OSTI2024

Distinct Mechanism of Anti-Corrosion and Swelling-Adhesion Modeling of Low-Dimensional Nylon-Fluoropolymer Composite Coatings

Compared to other polymers, composite coatings with fluoropolymers as the matrix have attracted considerable interest due to their mechanical performance, chemical stability, and low surface energy. Herein, we present an anticorrosion mechanism of fluoropolymer composite coatings achieved by dispersing varying amounts of polyamide 12 (PA-12) particles over a poly(vinylidene fluoride-co-hexafluoropropylene) (PVDF-HFP) matrix. The choice of PVDF-HFP as a coating matrix was driven by its superior processability and mechanical properties over Teflon, while PA-12 afforded swelling capacity, preventing electrolyte permeation, and inhibiting matrix stress crack propagation. The corrosion resistance of the coatings was assessed mainly by potentiodynamic polarization and impedance measurements while being immersed in a NaCl solution. Our electrochemical measurement findings showed that 0.75% w/w PA-12 in the PVDF-HFP matrix significantly enhanced corrosion protection even after 21 days of immersion. Microscopy, dielectric spectroscopy, surface analysis, and mechanical testing (American Society for Testing and Materials standards) corroborated the results. PVDF-HFP coatings containing a minimum percolation threshold of 0.75 wt % PA-12 with 1 mil thickness on mild steel exhibited a remarkable increase in resistance and film adhesion, and a notable corrosion rate decrease. Finite element analysis simulation with a bilinear traction-separation law confirmed the swelling mechanism and adhesion behavior. Finally, this study highlights the potential of nylon particle dispersion in PVDF-HFP matrices as a practical approach to developing advanced anticorrosion coatings with promising practical applications in various industries.

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

An active learning framework for the rapid assessment of galvanic corrosion

The current present in a galvanic couple can define its resistance or susceptibility to corrosion. However, as the current is dependent upon environmental, material, and geometrical parameters it is experimentally costly to measure. To reduce these costs, Finite Element (FE) simulations can be used to assess the cathodic current but also require experimental inputs to define boundary conditions. Due to these challenges, it is crucial to accelerate predictions and accurately predict the current output for different environments and geometries representative of in-service conditions. Machine learned surrogate models provides a means to accelerate corrosion predictions. However, a one-time cost is incurred in procuring the simulation and experimental dataset necessary to calibrate the surrogate model. Therefore, an active learning protocol is developed through calibration of a low-cost surrogate model for the cathodic current of an exemplar galvanic couple (AA7075-SS304) as a function of environmental and geometric parameters. The surrogate model is calibrated on a dataset of FE simulations, and calculates an acquisition function that identifies specific additional inputs with the maximum potential to improve the current predictions. This is accomplished through a staggered workflow that not only improves and refines prediction, but identifies the points at which the most information is gained, thus enabling expansion to a larger parameter space. The protocols developed and demonstrated in this work provide a powerful tool for screening various forms of corrosion under in-service conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Single-particle detection of enhanced polarizability in Au-decorated semiconducting nanorods via scanning dielectric microscopy

Hybrid nanostructures that combine semiconducting and metallic components offer great potential for photothermal therapy, optoelectronics, and sensing, by integrating tunable optical properties with enhanced light absorption and charge transport. Boosting the integrated performance of these hybrid systems demands techniques capable of probing local variations of the physical properties inaccessible to bulk analysis. Here, we report the single-particle dielectric characterization of hybrid, semiconducting bismuth sulfide (Bi 2 S 3 ) nanorods (NR) decorated with metallic Au nanoparticles (NP), employing scanning dielectric microscopy, which uses electrostatic force microscopy in combination with finite-element numerical simulations. We reveal a pronounced enhancement in the local dielectric response of Bi2S3 upon Au decoration, attributed to interfacial polarization and electron transfer from Au to the Bi 2 S 3 matrix, thus suggesting a enhanced metallic-like polarizability at the single-particle level. Numerical simulations show that the response is dominated by the vertical component of the permittivity and that the decorating metallic Au NP produce only moderate shielding of the semiconductor Bi 2 S 3 NR core, indicating that the large increase in the dielectric response originates primarily from intrinsic modifications within the NR. Overall, these findings provide direct insight into structure–property relationships at the single-particle level, supporting the rational design of advanced hybrid nanostructures with tailored electronic functionalities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Probing the properties and mechanisms of failure waves in soda-lime glass

Soda-lime glass (SLG) and other silica glasses exhibit the failure wave phenomenon under shock compression. The mechanism responsible for this peculiar behavior of glasses is still unresolved. In this study, a series of plate impact experiments was performed at three different impact stresses of 6.4, 8.3, and 10.8 GPa to better understand the mechanisms underlying the failure wave phenomenon. Specifically, spall experiments were conducted to probe the speed and existence of failure waves at different stresses in SLG. A layered glass target was used to probe the possibility of a “renucleation” of the wave at the SLG–SLG interface. When it existed, the failure wave was inferred to propagate at a speed of 1.3 km/s. However, it was observed that the failure wave phenomenon ceases to exist for impact stresses higher than 10 GPa. In experiments with a 6.4 GPa impact stress, the peak free surface velocity was significantly less than what is predicted by stress-Hugoniot calculations. This velocity deficit and other important features of the measured free surface velocity profiles were simulated using finite element analysis by incorporating an abrupt densification of SLG at a critical stress in the equation of state. Furthermore, this densification feature is similar to what would be expected of a phase transition. Although unable to unambiguously reveal the mechanism causing the failure wave phenomenon, the results of the present work clearly indicate that the failure wave causes a secondary compression and densification in SLG.

36 MATERIALS SCIENCE↗