Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element analysis”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Development of a thermal creep model for aluminum alloy 6061 cladding in U-10Mo monolithic fuel plates

Plate-type fuel elements consisting of a high-density, low-enriched uranium (LEU) U–10Mo-based fuel foil encapsulated in an aluminum alloy (AA) cladding are fabricated using the hot isostatic pressing (HIP) technique. During the HIP process, the fuel plate system is heated to 560 °C, then cooled to room temperature. This heat cycle significantly affects the mechanical properties of the aluminum cladding, and experimental investigations have shown that, post-HIP bonding, the mechanical properties of the aluminum cladding transition from those of AA 6061-T6 to something closer to the O temper. More specifically, the ultimate strength of the cladding decreases while its ductility increases, making it challenging to capture the changes in mechanical behavior and material properties. Understanding the residual stresses generated during the HIP process is critical for assessing the fuel plate’s integrity under various temperature, pressure, and irradiation. To simulate the HIP bonding process, the elastic, plastic, and thermal properties of the cladding are assumed to be similar to those of AA 6061-O temper. However, the primary challenge lies in the lack of available data for the creep model of the AA 6061 cladding during this transient process of HIP. The present study focuses on developing a computational model that predicts the creep behavior of the aluminum cladding in the fuel plates during the HIP process, as cladding creep significantly influences the residual stresses generated in U-10Mo fuel plates during HIP fabrication. Furthermore, as HIP bonding occurs at high temperatures that are nearing the melting point of aluminum, the present work considered a temperature-dependent Arrhenius-type creep model. In particular, a hyperbolic sine creep model is employed to estimate the creep properties of the as-fabricated aluminum cladding. In conclusion, the residual stresses predicted in the U-10Mo fuel when using the newly calibrated creep model closely align with the experimental measurements, validating the model’s accuracy.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

High temperature fracture toughness of single-layer CrAlN and CrAlSiN hard coatings

Hardness and toughness are among of the most important physical properties of hard coatings and it is important to analyze them at temperatures present during coating applications. CrAlN and CrAlSiN are hard coating materials frequently used in high temperature applications. Here, while both possess high oxidation resistance and thermal stability, they differ in their microstructure. In this work, the fracture toughness of magnetron-sputtered single-layer CrAlN and CrAlSiN hard coatings was determined up to 700 °C by in-situ scanning electron microscopy micro-cantilever deflection tests. The results show a decreasing trend of the fracture toughness with increasing temperature. While both coatings' fracture toughness trends are similar up to 400 °C, they diverge at higher temperatures. In general, the evolution of the fracture toughness with increasing temperature shows similarities to the observed changes in hardness in the same temperature range.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Regiospecific Cellulose Orientation and Anisotropic Mechanical Property in Plant Cell Walls

Cellulose microfibrils (CMFs) are a major load-bearing component in plant cell walls. Thus, their structures have been studied extensively with spectroscopic and microscopic characterization methods; but the findings from these two approaches were inconsistent, which hampers mechanistic understanding of cell wall mechanics. Here, we report the regiospecific assembly of CMFs in the periclinal wall of plant epidermal cells. Using sum frequency generation spectroscopic imaging, we found that CMFs are highly aligned in the cell edge region where two cells form a junction, whereas they are mostly isotropic on average throughout the wall thickness in the flat face region of the epidermal cell. This subcellular-level heterogeneity in CMF alignment provided insights into tissue-level anisotropy in tensile modulus of cell wall materials. Furthermore, this finding also has resolved a previous contradiction between the spectroscopic and microscopic imaging studies, which paves a foundation for better understanding of cell wall artichecture, especially structure-geometry relationships.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Effect of H + Exchange and Surface Impurities on Bulk and Interfacial Electrochemistry of Garnet Solid Electrolytes

Contact loss and current constriction pose significant challenges at the Li metal interface of solid-state batteries. For garnet-structured Li 7 La 3 Zr 2 O 12 (LLZO), these effects are amplified by Li + /H + exchange and surface contamination reactions, which lead to conductivity losses and poor Li wetting. In this study, we utilize a variety of surface treatment processes across 37 cells to selectively induce proton exchange and contamination reactions in LLZO. The resulting bulk and surface chemistry is systematically characterized and correlated to changes in electrochemical properties. Additionally, we combine impedance analysis and finite element method modeling to deconvolute sources of impedance contributions at the Li metal interface. Specifically, we show that constriction impedance at the Li metal interface arises not solely from voids, but also from ionically-resistive surface contaminants. Further, these findings emphasize the connection between ionic conductivity and constriction, demonstrating that micron-scale ionically-resistive components increase constriction even with identical contact geometries. Finally, we leverage our comprehensive dataset to highlight unstable overpotential growth as a failure mechanism, additionally showing that the phase of a cell’s impedance is a sensitive indicator for the onset of interfacial instability. Overall, this study clarifies the impacts of proton exchange and surface contamination on electrochemical properties at the Li|solid electrolyte interface and elucidates insights that are generalizable to other solid-state battery systems.

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

Intrinsic thermal interfacial resistance measurement in bonded metal–polymer foils

Heat conduction through bonded metal–polymer interfaces often limits the overall heat transfer in electronic packaging, batteries, and heat recovery systems. To design the thermal circuit in such systems, it is essential to measure the thermal interfacial resistance (TIR) across ~1 µm to 100 µm junctions. Previously reported TIR of metal–polymer junctions utilize ASTM E1530-based two-block systems that measure the TIR by applying pressure across the interface through external heating and cooling blocks. Here, we report a novel modification of the ASTM-E1530 technique that employs integrated heaters and sensors to provide an intrinsic TIR measurement of an adhesively bonded metal–polymer junction. We design the measurement technique using finite element simulations to either passively suppress or actively compensate the lateral heat diffusion through the polymer, which can minimize the systematic error to ≲5%. Through proof-of-concept experiments, we report the TIR of metal–polymer interfaces made from DuPont’s Pyralux double-side copper-clad laminates, commonly used in flexible printed circuit boards. Our TIR measurement errors are <10%. We highlight additional sources of errors due to non-idealities in the experiment and discuss possible ways to overcome them. Our measurement technique is also applicable to interfaces that are electrically insulating such as adhesively joined metal–metal junctions and sputter-coated or welded metal–polymer junctions. Altogether, the technique is capable of measuring TIR ≳10 –5 m 2 KW –1 in bonded metal–polymer foils and can be tailored for in situ measurements in flexible electronics, circuit packaging, and other hybrid metal–polymer systems.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Dynamically morphing microchannels in liquid crystal elastomer coatings containing disclinations

Liquid crystal elastomers (LCEs) hold a major promise as a versatile material platform for smart soft coatings, since their orientational order can be predesigned to program a desired dynamic profile. In this work, we introduce temperature-responsive dynamic coatings based on LCEs with arrays of singular defects-disclinations that run parallel to the surface. The disclinations form in response to antagonistic patterns of the molecular orientation at the top and bottom surfaces, imposed by the plasmonic mask photoalignment. Upon heating, an initially flat LCE coating develops linear microchannels located above each disclination. The stimulus that causes a non-flat profile of LCE coatings upon heating is the activation force induced by the gradients of molecular orientation around disclinations. To describe the formation of microchannels and their thermal response, we adopt a Frank-Oseen model of disclinations in patterned director field and propose a linear elasticity theory to connect the complex spatially-varying molecular orientation to the displacements of the LCE. The thermo-responsive surface profiles predicted by the theory and by the finite element modeling are in good agreement with the experimental data; in particular, higher gradients of molecular orientation produce a stronger modulation of the coating profile. The elastic theory and the finite element simulations allow us to estimate the material parameter that characterizes the elastomer coating’s response to the thermal activation. The disclinations containing LCEs show potential as soft dynamic coatings with a predesigned responsive surface profile.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Directional shock diode behavior through the interaction of geometric voids in engineered polymer assemblies

With the advent of additive manufacturing (AM) techniques, a new class of shockwave mitigation and structural supports has been realized through the hierarchical assembly of polymer materials. To date, there have been a limited number of studies investigating the role of structure on shockwave localization and whether AM offers a means to tailor shockwave behavior. Of particular interest is whether the mesoscopic structure can be tailored to achieve shockwave properties in one direction of impact vs the other. In this work, we illustrate directional response in engineered polymer foams. In situ time-resolved x-ray phase contrast imaging at the Advanced Photon Source was used to characterize these diode-like structures. This work offers a breakthrough in materials technology for the development of protective structures that require augmentation of shock in one direction while diminishing transmission in the opposite direction.

42 ENGINEERING↗
DOE OSTI2021

Time integrator agnostic charge conserving finite element PIC

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a longstanding problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leapfrog time stepping method. While this is a significant advance, leapfrog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss’ laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the Debye length has been resolved, and real mass ratios are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Thermal conductivity modeling of monodispersed microspheres using discrete element method

Particle beds are widely used in various systems and processes, such as particle heat exchangers, granular flow reactors, and additive manufacturing. Accurate modeling of the thermal conductivity of particle beds and understanding of their heat transfer mechanisms are important. However, previous models were based on simple cubic packing of particles, which could not accurately represent the actual heat transfer processes under certain conditions. Here, we examine the effect of the packing structure on the thermal conductivity of particle beds. We use monodispersed silica microspheres with average particle sizes ranging from 23 to 330 μm as a model material. We employ a transient hot-wire technique to measure the thermal conductivity of the particle beds with packing density of 43%–57% within a temperature range of room temperature to 500 °C and under N2 gaseous pressures of 20–760 Torr. We then use a discrete element method (DEM) to obtain the realistic packing structure of the particles, which is then fed into a finite-element model (FEM) to calculate the thermal conductivity, with the consideration of solid conduction, gas conduction, and radiation heat transfer. Our results show that the thermal conductivity model based on the more realistic random packing structure derived from the DEM shows better agreement with the experimental data compared to that based on the simple cubic-packing structure. The combined DEM and FEM methodology can serve as a useful tool to predict the effective thermal conductivity of particle beds and to quantify different heat transfer mechanisms under various conditions.

14 SOLAR ENERGY↗
DOE OSTI2021

A semi-analytical model for laser induced densification in fused silica at elevated temperatures

A semi-analytical model for deterministic surface smoothing of fused silica below the evaporation temperature threshold is presented. Using this model, a method for calculating deterministic track depths for laser polishing fused silica is derived and validated with experiments. In this work, model predictions were in good agreement with experiment across varying laser preheating temperature and power. With the semi-analytical method presented here, deterministic laser smoothing can be achieved in the densification regime to close the process loop of laser polishing.

36 MATERIALS SCIENCE↗
DOE OSTI2021

THz time-domain characterization of amplifying quantum-cascade metasurface

An amplifying quantum-cascade (QC) metasurface, the key component of the QC vertical-external-cavity surface-emitting-laser (VECSEL), is studied as a function of injected current density using reflection-mode terahertz time domain spectroscopy. Nearly perfect absorption is measured at zero bias, which is associated with the transition from the weak to strong coupling condition between the metasurface resonance and an intersubband transition within the QC material. An increase in reflectance is observed as the device is biased, both due to reduction in intersubband loss and the presence of intersubband gain. Significant phase modulation associated with the metasurface resonance is observed via electrical control, which may be useful for electrical tuning of QC-VECSEL. These results provide insight into the interaction between the intersubband QC-gain material and the metasurface and modify the design rules for QC-VECSELs for both biased and unbiased regions.

42 ENGINEERING↗
DOE OSTI2022

An E and B gyrokinetic simulation model for kinetic Alfvén waves in tokamak plasmas

The gyrokinetic particle simulation is a powerful tool for studies of transport, nonlinear phenomenon, and energetic particle physics in tokamak plasmas. While most gyrokinetic simulations make use of the scalar and vector potentials, a new model (GK-E&B) has been developed by using the E and B field in a general form and has been implemented in simulating kinetic Alfvén waves in uniform plasma. In our work, the Chen et al. GK-E&B model has been expressed, in general, tokamak geometry using the local orthogonal coordinates and general tokamak coordinates. Its reduction for uniform plasma is verified, and the numerical results show good agreement with the original work. The theoretical dispersion relation and numerical results in the local model in screw pinch geometry are also in excellent agreement. Numerical results show excellent performance in a realistic parameter regime of burning plasmas with high values of β/(M e k$^{2}_{⊥}$ρ$^{2}_{i}$), which is a challenge for traditional methods due to the “cancellation” problem. As one application, the GK-E&B model is implemented with kinetic electrons in the local single flux surface limit. With the matched International Tokamak Physics Activity-Toroidicity-induced Alfvén Eigenmodes parameters adopted, numerical results show the capability of the GK-E&B in treating the parallel electron Landau damping for realistic tokamak plasma parameters. As another application, the global GK-E&B model has been implemented with the dominant electron contribution in the cold electron limit. Its capability in simulating the finite E || due to the finite electron mass is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Controlling superconductivity of CeIrIn 5 microstructures by substrate selection

Superconductor/metal interfaces are usually fabricated in heterostructures that join these dissimilar materials. A conceptually different approach has recently exploited the strain sensitivity of heavy-fermion superconductors, selectively transforming regions of the crystal into the metallic state by strain gradients. The strain is generated by differential thermal contraction between the sample and the substrate. Here, we present an improved finite-element model that reliably predicts the superconducting transition temperature in CeIrIn 5 even in complex structures. Different substrates are employed to tailor the strain field into the desired shapes. Using this approach, both highly complex and strained as well as strain-free microstructures are fabricated to validate the model. This enables a high degree of control over the microscopic strain fields and forms the basis for more advanced structuring of superconductors as in Josephson junctions yet also finds natural use cases in any material class in which a modulation of the physical properties on a chip is desirable.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2022

Sensing depths in frequency domain thermoreflectance

In this work, a method is developed to calculate the length into a sample to which a Frequency Domain Thermoreflectance (FDTR) measurement is sensitive. Sensing depth and sensing radius are defined as limiting cases for the spherically spreading FDTR measurement. A finite element model for FDTR measurements is developed in COMSOL multiphysics and used to calculate sensing depth and sensing radius for silicon and silicon dioxide samples for a variety of frequencies and laser spot sizes. The model is compared to experimental FDTR measurements. Design recommendations for sample thickness are made for experiments where semi-infinite sample depth is desirable. For measurements using a metal transducer layer, the recommended sample thickness is three thermal penetration depths, as calculated from the lowest measurement frequency.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulating the non-monotonic strain response of nanoporous multiferroic composites under electric field control

In this work, we simulate and analyze the mechanical response of a class of multiferroic materials consisting of a templated porous nanostructure made out of cobalt ferrite (CFO) partially filled by atomic layer deposition (ALD) with a ferroelectric phase of lead zirconate titanate (PZT). The strain in the device is measured when an electric field is applied for varying ALD thicknesses, displaying a non-monotonic dependence with a maximum strain achieved for a coating thickness of 3 nm. To understand this behavior, we apply finite element modeling to the smallest repeatable unit of the nanoporous template and simulate the mechanical response as a function of PZT coating thickness. We find that this non-monotonic response is caused by the interplay between two driving forces opposing one another. First, increased porosity works toward increasing the strain due to a reduced system stiffness. Second, decreased porosity involves a larger mass fraction of PZT, which drives the electro-mechanical response of the structure, thus leading to a larger strain. The balance between these two driving forces is controlled by the shear coupling at the CFO/PZT interface and the effective PZT cross section along the direction of the applied electric field. Here, our numerical results show that considering a nonlinear piezoelectric response for PZT leads to an improved agreement with the experimental data, consistent with ex situ poling of the nanostructure prior to magnetic measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Charge carrier motion and effect of fixed oxide charge in a microstructured silicon radiation detector

Signal formation in a microstructured semiconductor neutron detector is more complex than in planar diode geometry. Three-dimensional microstructures are laterally smaller than the ionization cloud length, and the electric fields may be weak enough to exhibit plasma time effects. This work is the first detailed treatment of charge carrier motion in these complex semiconductor devices to replicate the time profile and signal magnitude. Simulations were performed using COMSOL Multiphysics to investigate various parameters that affect the propagation of the charge cloud. It was observed that the size of the simulated three-dimensional structure had an impact on the induced current pulse, indicating the importance of simulation geometry optimization to accurately simulate charge cloud expansion. COMSOL Multiphysics was used to replicate accurate charge creation profiles using energy deposition information imported from radiation transport codes. A detailed simulation methodology is presented to benchmark preamplifier event pulses along with complexities in modeling the charge carrier motion along the etched microstructured trenches with Si–SiO 2 boundary conditions, including fixed oxide charge and interface trapping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗