Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Cell migration and organization in three-dimensional in vitro culture driven by stiffness gradient: Cell Migration Driven by Stiffness Gradient

Durotaxis, a phenomenon that cells move according to changes in stiffness of the extra cellular matrix, has emerged as a crucial parameter controlling cell migration behavior. The current study provides a simple method to generate three-dimensional continuous stiffness variations without changing other physical characteristics of the extra cellular environment. Using Finite Element simulations, the stiffness and the stiffness gradient variations are evaluated quantitatively, leading to an analysis of the dependence of cell migration behavior on the substrate stiffness parameters. We tested various cell lines on several 3-D environments. In this work, the durotaxis results show that the cell migration velocity does not have any consistency with the stiffness of the substrate, rather it is more related to the stiffness gradient of the substrate. This finding suggests a new mechanism underlying the durotaxis phenomenon, highlighting the importance of the substrate stiffness gradient, rather than the stiffness itself.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A performance spectrum for parallel computational frameworks that solve PDEs

Summary Important computational physics problems are often large‐scale in nature, and it is highly desirable to have robust and high performing computational frameworks that can quickly address these problems. However, it is no trivial task to determine whether a computational framework is performing efficiently or is scalable. The aim of this paper is to present various strategies for better understanding the performance of any parallel computational frameworks for solving PDEs. Important performance issues that negatively impact time‐to‐solution are discussed, and we propose a performance spectrum analysis that can enhance one's understanding of critical aforementioned performance issues. As proof of concept, we examine commonly used finite element simulation packages and software and apply the performance spectrum to quickly analyze the performance and scalability across various hardware platforms, software implementations, and numerical discretizations. It is shown that the proposed performance spectrum is a versatile performance model that is not only extendable to more complex PDEs such as hydrostatic ice sheet flow equations but also useful for understanding hardware performance in a massively parallel computing environment. Potential applications and future extensions of this work are also discussed.

Chang, J.↗
DOE OSTI2023

Modeling of the effectiveness of novel edge seal designs for fast, low‐Cap‐Ex manufacturing

Abstract In the manufacturing of photovoltaic (PV) modules, the lamination process can take up to 20 min to complete. In this work, new lamination processes are being developed, and have been prototyped, which hope to be able to cut this time down to as little as 30 s. This could provide significant savings in the cost of lamination equipment, floor space, and energy. PV modules are expected to have a lifespan exceeding 20 to 30 years. For moisture‐sensitive PV technologies, the edge seal between the two layers of glass can be the weakest point of its reliability. There is an inherent challenge when evaluating edge seal materials due to their low permeation rates. As part of Colorado State University's Photovoltaic Research and Development 2, work at the National Renewable Energy Laboratory has developed models to evaluate edge seal configurations in glass‐glass PV modules. Here, this new manufacturing process is evaluated for long‐term moisture durability. Different edge seal design options within glass–glass PV modules are explored. Most of these designs are targeting a superstrate on glass configuration, e.g CdTe, but some designs could be used on conventional crystalline Si cells. Using COMSOL finite element simulation software, we investigated the edge seal and interlayer design configurations containing silicone perimeter edge adhesive, desiccated polyisobutylene‐based edge seal, air, and polyolefin while integrating climate conditions equivalent to a hot and humid environment such as Miami, Florida. We found optimized configurations that will allow the module to prevent moisture ingress over 50 years minimizing the amount of time and material used while utilizing polymers that are easily dispensed.

14 SOLAR ENERGY↗
DOE OSTI2016

Anomalous Oxidative Diffusion in Titanium Pyrotechnic Powders

It has long been observed that oxidation processes in metals tend to follow a parabolic rate law associated with the growth of a surface oxide layer. Here we observe that for certain titanium powders, the expected parabolic law (∝t 1/2 ) is recovered, yet for others, the exponent differs significantly. One explanation for this non-parabolic, anomalous diffusion arises from fractal geometry. Theoretical considerations indicate that the time response of diffusion-limited processes in an object closely follow a power-law in time (t n ) with n=(E–D)/2, where E is the object's Euclidean dimension and D is its boundary's Hausdorff dimension. Non-integer, (fractal) values of D will result in n≠1/2. Finite element simulations of several canonical fractal objects were performed to verify the application of this theory; the results matched the theory well. Two different types of titanium powder were tested in isothermal thermogravimetric tests under dilute oxygen. Time-dependent mass uptake data were fit with power-law forms and the associated exponents were used to determine an equivalent fractal dimension. One Ti powder type has an implied surface dimension of ca. 2.3 to 2.5, suggesting fractal geometry may be operative. Finally, the other has a dimension near 2.0, indicating it behaves like traditional material.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Grain Orientation Dependence of the Residual Lattice Strain in a Cold Rolled Interstitial-Free Steel

The experimentally measured grain-orientation-dependent residual lattice strains, evolved in an interstitia-free steel after 70% cold rolling reduction, are studied by means of crystal elastic visco-plastic finite element simulations, which provides a very satisfactory prediction of deformation texture. The calculated residual lattice strain pole figure matches well with the experimentally measured counterpart within the highest density regions of major texture components observed. Both experimental evidence and results of modeling clearly indicate that the residual lattice strain is orientation dependent, based on comprehensive information on the evolution of residual lattice strain in various crystallographic orientations during plastic deformation. It appears that in a cold rolled material, there is a general correlation between the stresses developed just prior to unloading and the residual lattice strains in particular directions. Here, it is also shown that the cumulative plastic shear does not reveal a clear correlation with the components of residual lattice strain while presented in the normal correlation plot, however, this relationship can be better understood by means of the orientation distribution function of residual lattice strain, which can be derived from the neutron or X-ray diffraction experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Input specific neural networks

Neural networks have emerged as powerful tools for mapping between inputs and outputs. However, their black-box nature limits the ability to encode or impose specific structural relationships between inputs and outputs. Many scientific and engineering problems, such as constitutive modeling in solid mechanics, require networks that can enforce convexity, monotonicity, or other structural constraints to ensure physical consistency. Here, we introduce the Input Specific Neural Network (ISNN), a new architecture that enables multiple, distinct constraints to be imposed on different input subsets for scalar-valued outputs. This framework unifies convex, monotone–convex, monotone, and arbitrary mappings within a single network for the first time. Two ISNN architectures with analytical first- and second-order derivatives are developed. We demonstrate the performance on synthetic toy problems, inverse problems in isotropic hyperelasticity, and finite element simulations. ISNNs achieve improved extrapolation behavior, require fewer invariant inputs than standard input convex networks for polyconvex potentials, and enable significant computational savings via manual differentiation. We also show how ISNNs can be used to learn structural relationships between inputs and outputs via a binary gating mechanism. Particularly, ISNNs are employed to model a homogenized anisotropic free energy potential in a decoupled multiscale setting. The network learns whether or not the potential should be modeled as polyconvex and retains only the relevant layers while using the minimum number of inputs. ISNNs provide a flexible foundation for embedding structural priors into neural networks, enhancing both interpretability and stability. They are broadly applicable across computational mechanics and other scientific domains requiring constrained functional relationships.

Jadoon, Asghar A. [Univ. of Texas, Austin, TX (Uni↗
DOE OSTI2015

On-demand in situ generation of oxygen in a nanofluidic embedded planar microband electrochemical reactor

In situ generation of reagents and their subsequent use downstream presents new opportunities to amplify the utility of nanofluidic devices by exploiting the confined geometry to address mass transport limitations on reaction kinetics and efficiency. Oxygen, an inherently valuable reactant, can be produced from electrolysis of water, a process that can be conveniently integrated within a nanofluidic system. Here, we construct and characterize a nanofluidic device consisting of a planar microband electrode embedded within a nanochannel for in situ electrochemical generation and optical monitoring of O 2 . Fluorescein, a dye with a pH-sensitive emission intensity, was used to monitor the spatiotemporal characteristics of the oxidation of H 2 O, using the co-produced H + . Application of anodic potentials at the nanochannel-embedded electrode results in a decrease in fluorescence intensity, which reflects the decreasing solution pH. A combination of fluorescence intensity and chronoamperometric response was used to quantitatively determine proton generation, and the H + /O 2 stoichiometry was then used to determine the concentration of the O 2 in the channel. Comparison of the experimental results to finite element simulations validates the use of fluorescein emission intensity to spectroscopically determine the local oxygen concentration in the nanochannel. By varying the applied potential, spatially averaged O 2 concentrations ranging from 0.13 to 0.41 mM were generated. In conclusion, the results demonstrate a convenient route to in situ modulation of the dissolved O 2 level in a nanofluidic device and the use of an optical probe to monitor its spatial and temporal distribution under flow conditions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Determining anisotropic slip system rate sensitivities of Ti-6Al-4V using high-energy X-ray diffraction microscopy

Variation of strain rate sensitivity among different families of slip systems in the hexagonal close-packed (α) phase of titanium (Ti) alloys has the potential to alter microscale load redistribution during both creep and dwell fatigue loading. However, existing literature contains conflicting reports regarding the degree of anisotropy present in the strain rate sensitives between α slip system families across Ti alloys. Here, we quantify the strain rate sensitivity of α slip system families in Ti-6Al-4V using high-energy X-ray diffraction microscopy (HEDM). We present a novel procedure in which we utilize the HEDM-measured grain-scale stress states during stress relaxation to determine the strain rate sensitivities of the basal $\langle a \rangle$ , prismatic $\langle a \rangle$, and first-order pyramidal $\langle c + a \rangle$ slip system families. In addition, we measure the strain rate sensitivity exponent of the macroscopic response and the α phase for comparison. We find rate sensitivities for the different slip system families to range from 0.02 to 0.04, which—while varied—are relatively isotropic in comparison with some values presented in the literature. We also demonstrate the effects of the measured anisotropic rate sensitivities using crystal plasticity finite element simulations.

Peterson, Kenneth M. [Pennsylvania State Univ., Un↗
DOE OSTI2014

Three-Dimensional Characterization and Modeling of Microstructural Weak Links for Spall Damage in FCC Metals

Local microstructural weak links for spall damage were investigated using three-dimensional (3-D) characterization in multicrystalline copper samples (grain size ≈ 450 µm) shocked with laser-driven plates at low pressures (2 to 4 GPa). The thickness of samples and flyer plates, approximately 1000 and 500 µm respectively, led to short pressure pulses that allowed isolating microstructure effects on local damage characteristics. Electron Backscattering Diffraction and optical microscopy were used to relate the presence, size, and shape of porosity to local microstructure. The experiments were complemented with 3-D finite element simulations of individual grain boundaries (GBs) that resulted in large damage volumes using crystal plasticity coupled with a void nucleation and growth model. Results from analysis of these damage sites show that the presence of a GB-affected zone, where strain concentration occurs next to a GB, correlates strongly with damage localization at these sites, most likely due to the inability of maintaining strain compatibility across these interfaces, with additional effects due to the inclination of the GB with respect to the shock. Results indicate that strain compatibility plays an important role on intergranular spall damage in metallic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Insights from the Lattice-Strain Evolution on Deformation Mechanisms in Metallic-Glass-Matrix Composites

In this paper, in situ high-energy synchrotron X-ray diffraction experiments and micromechanics-based finite element simulations have been conducted to examine the lattice-strain evolution in metallic-glass-matrix composites (MGMCs) with dendritic crystalline phases dispersed in the metallic-glass matrix. Significant plastic deformation can be observed prior to failure from the macroscopic stress–strain curves in these MGMCs. The entire lattice-strain evolution curves can be divided into elastic–elastic (denoting deformation behavior of matrix and inclusion, respectively), elastic–plastic, and plastic–plastic stages. Characteristics of these three stages are governed by the constitutive laws of the two phases (modeled by free-volume theory and crystal plasticity) and geometric information (crystalline phase morphology and distribution). The load-partitioning mechanisms have been revealed among various crystalline orientations and between the two phases, as determined by slip strain fields in crystalline phase and by strain localizations in matrix. Finally, implications on ductility enhancement of MGMCs are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Triaxial Constraint and Tensile Strength Enhancement in Brazed Joints

A brazed joint consists of a low-melting point and thin interlayer sandwiched between the high-melting-point base materials, in which the interlayer strength is typically lower than that of the base material. When this butt-joined composite is loaded uniaxially in the direction perpendicular to the plane of the brazing layer, the tensile strength is found to be much higher than that of the braze. This seems to violate the iso-stress condition in such a butt-joint serial configuration. Furthermore, the stress triaxiality has been usually ascribed, but without a quantitative rationalization, as being responsible for this tensile strength enhancement. Here a complete finite element simulation has been conducted to study the dependence of triaxiality and strength enhancement on geometric and material parameters. Two asymptotic limit solutions (based on Bridgman and Xia–Shih solutions, respectively) have been identified to understand the simulation results. The critical role of void evolution has been revealed when making a quantitative comparison to available experiments. In addition, ductility of the brazed joint, which has not been fully addressed in literature, is investigated by the Gurson–Tvergaard–Needleman model.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Intergranular Strain Evolution During Biaxial Loading: A Multiscale FE-FFT Approach

Predicting the macroscopic and microscopic mechanical response of metals and alloys subjected to complex loading conditions necessarily requires a synergistic combination of multiscale material models and characterization techniques. This article focuses on the use of a multiscale approach to study the difference between intergranular lattice strain evolution for various grain families measured during in situ neutron diffraction on dog bone and cruciform 316L samples. At the macroscale, finite element simulations capture the complex coupling between applied forces and gauge stresses in cruciform geometries. The predicted gauge stresses are used as macroscopic boundary conditions to drive a mesoscale full-field elasto-viscoplastic fast Fourier transform crystal plasticity model. The results highlight the role of grain neighborhood on the intergranular strain evolution under uniaxial and equibiaxial loading.

36 MATERIALS SCIENCE↗
DOE OSTI2016

On capturing the grain-scale elastic and plastic anisotropy of alpha-Ti with spherical nanoindentation and electron back-scattered diffraction

Here, spherical nanoindentation combined with electron back-scattered diffraction has been employed to characterize the grain-scale elastic and plastic anisotropy of single crystal alpha-Ti of two different compositions (in two different titanium alloys). Data analyses protocols needed to reliably extract the desired properties of interest are extended and demonstrated in this paper. Specifically, the grain-scale mechanical response is extracted in the form of indentation stress-strain curves for commercially pure (CP-Ti) alpha-Ti and alloyed (Ti-64) titanium from measurements on polycrystalline samples. The results are compared with responses of single crystals and nanoindentation tests (hardness and modulus) from the literature, and the measured indentation moduli are validated using crystal-elastic finite element simulations. The results obtained in this study show that (i) it is possible to characterize reliably the elastic and plastic anisotropy of alpha-Ti (hcp) of varying alloying contents with spherical nanoindentation stress-strain curves, (ii) the indentation modulus of alpha-Ti-64 is 5–10% less than CP-Ti, and (iii) the indentation yield strength of alpha-Ti-64 is 50–80% higher than CP-Ti.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Modeling slip system strength evolution in Ti-7Al informed by in-situ grain stress measurements

In this work, far-field high-energy X-ray diffraction microscopy is used to asses the evolution of slip system strengths in hexagonal close-packed (HCP) Ti-7Al during tensile deformation in-situ. The following HCP slip system families are considered: basal $\langle$a$\rangle$, prismatic $\langle$a$\rangle$, pyramidal$\langle$a$\rangle$, and first-order pyramidal $\langle$c + a$\rangle$. A 1 mm length of the specimen's gauge section, marked with fiducials and comprised of an aggregate of over 500 grains, is tracked during continuous deformation. The response of each slip system family is quantified using ‘slip system strength curves’ that are calculated from the average stress tensors of each grain over the applied deformation history. These curves, which plot the average resolved shear stress for each slip system family versus macroscopic strain, represent a mesoscopic characterization of the aggregate response. A short time-scale transient softening is observed in the basal $\langle$a$\rangle$, prismatic $\langle$a$\rangle$, and pyramidal $\langle$a$\rangle$ slip systems, while a long time-scale transient hardening is observed in the pyramidal $\langle$c + a$\rangle$ slip systems. These results are used to develop a slip system strength model as part of an elasto-viscoplastic constitutive model for the single crystal behavior. A suite of finite element simulations is performed on a virtual polycrystal to demonstrate the relative effects of the different parameters in the slip system strength model. The model is shown to accurately capture the macroscopic stress-strain response using parameters that are chosen to capture the mesoscopic slip system responses.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Mechanical properties of metal-ceramic nanolaminates: Effect of constraint and temperature

Al/SiC nanolaminates with equal nominal thicknesses of the Al and SiC layers (10, 25, 50 and 100 nm) were manufactured by magnetron sputtering. The mechanical properties were measured at 25 °C and 100 °C by means of nanoindentation and micropillar compression tests and the deformation mechanisms were analyzed by in situ micropillar compression tests in the transmission electron microscope. In addition, finite element simulations of both tests were carried out to ascertain the role played by the strength of the Al layers and by the elastic constraint of the ceramic layers on the plastic flow of Al in the mechanical response. It was found that the mechanical response was mainly controlled by the constraint during nanoindentation or micropillar compression tests of very thin layered (≈10 nm) laminates, while the influence of the strength of Al layers was not as critical. This behavior was reversed, however, for thick layered laminates (100 nm). Here, these mechanisms point to the different effects of layer thickness during nanoindentation and micropillar compression, at both temperatures, and showed the critical role played by constraint on the mechanical response of nanolaminates made of materials with a very large difference in the elasto-plastic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Particle size effects in metallic microparticle impact-bonding

In cold spray processing of metals, the critical velocity for particle bonding is dependent on both intrinsic material properties, such as density and (spall) strength, and extrinsic process parameters, such as the powder particle size. Here, we specifically isolate and investigate particle size effects on the critical velocity for bonding through laser- induced single particle impact experiments and finite element simulations. We also present a predictive framework to correlate particle size and critical velocity. We show that an increase in particle size leads to an increase in the temperature of the jet formed at the interface of the particle and the substrate. This increased temperature locally decreases the spall strength of the material which, in turn, decreases the critical velocity for larger particles.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Moisture, anisotropy, stress state, and strain rate effects on bighorn sheep horn keratin mechanical properties

Our paper investigates the effects of moisture, anisotropy, stress state, and strain rate on the mechanical properties of the bighorn sheep (Ovis Canadensis) horn keratin. The horns consist of fibrous keratin tubules extending along the length of the horn and are contained within an amorphous keratin matrix. We tested samples in the rehydrated (35 wt.% water) and ambient dry (10 wt.% water) conditions along the longitudinal and radial directions under tension and compression. Increased moisture content was found to increase ductility and decrease strength, as well as alter the stress state dependent nature of the material. Furthermore, the horn keratin demonstrates a significant strain rate dependence in both tension and compression, and also showed increased energy absorption in the hydrated condition at high strain rates when compared to quasi-static data, with increases of 114% in tension and 192% in compression. Compressive failure occurred by lamellar buckling in the longitudinal orientation followed by shear delamination. Tensile failure in the longitudinal orientation occurred by lamellar delamination combined with tubule pullout and fracture. Finally, the structure-property relationships quantified here for bighorn sheep horn keratin can be used to help validate finite element simulations of ram’s impacting each other as well as being useful for other analysis regarding horn keratin on other animals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiscale mechanical design of the lightweight, stiff, and damage-tolerant cuttlebone: A computational study

Cuttlebone, the endoskeleton of cuttlefish, offers an intriguing biological structural model for designing low-density cellular ceramics with high stiffness and damage tolerance. Cuttlebone is highly porous (porosity ~93%) and lightweight (density less than 20% of seawater), constructed mainly by brittle aragonite (95 wt%), but capable of sustaining hydrostatic water pressures over 20 atmospheres and exhibits energy absorption capability under compression comparable to many metallic foams (~4.4 kJ/kg). Here, in this work, we computationally investigate how such remarkable mechanical efficiency is enabled by the multiscale structure of cuttlebone. Using the common cuttlefish, Sepia Officinalis, as a model system, we first conducted high-resolution synchrotron micro-computed tomography (µ-CT) and quantified the cuttlebone's multiscale geometry, including the 3D asymmetric shape of individual walls, the wall assembly patterns, and the long-range structural gradient of walls across the entire cuttlebone (ca. 38 chambers). The acquired 3D structural information enables systematic finite-element simulations, which further reveal the multiscale mechanical design of cuttlebone: at the wall level, wall asymmetry provides optimized energy absorption while maintaining high structural stiffness; at the chamber level, variation of walls (number, pattern, and waviness amplitude) contributes to progressive damage; at the entire skeletal level, the gradient of chamber heights tailors the local mechanical anisotropy of the cuttlebone for reduced stress concentration. Our results provide integrated insights into understanding the cuttlebone's multiscale mechanical design and provide useful knowledge for the designs of lightweight cellular ceramics.

36 MATERIALS SCIENCE↗