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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

In situ characterization of laser-generated melt pools using synchronized ultrasound and high-speed X-ray imaging

Metal additive manufacturing is a fabrication method that forms a part by fusing layers of powder to one another. An energy source, such as a laser, is commonly used to heat the metal powder sufficiently to cause a molten pool to form, which is known as the melt pool. The melt pool can exist in the conduction or the keyhole mode where the material begins to rapidly evaporate. The interaction between the laser and the material is physically complex and difficult to predict or measure. Here, high-speed X-ray imaging was combined with immersion ultrasound to obtain synchronized measurements of stationary laser-generated melt pools. Furthermore, two-dimensional and three-dimensional finite-element simulations were conducted to help explain the ultrasonic response in the experiments. In particular, the time-of-flight and amplitude in pulse-echo configuration were observed to have a linear relationship to the depth of the melt pool. These results are promising for the use of ultrasound to characterize the melt pool behavior and for finite-element simulations to aid in interpretation.

42 ENGINEERING↗
DOE OSTI2024

Facile Tensile Testing Platform for In Situ Transmission Electron Microscopy of Nanomaterials

In situ tensile testing using transmission electron microscopy (TEM) is a powerful technique to probe structure-property relationships of materials at the atomic scale. In this work, a facile tensile testing platform for in situ characterization of materials inside a transmission electron microscope is demonstrated. The platform consists of: 1) a commercially available, flexible, electron-transparent substrate (e.g., TEM grid) integrated with a conventional tensile testing holder, and 2) a finite element simulation providing quantification of specimen-applied strain. The flexible substrate (carbon support film of the TEM grid) mitigates strain concentrations usually found in free-standing films and enables in situ straining experiments to be performed on materials that cannot undergo localized thinning or focused ion beam lift-out. The finite element simulation enables direct correlation of holder displacement with sample strain, providing upper and lower bounds of expected strain across the substrate. The tensile testing platform is validated for three disparate material systems: sputtered gold-palladium, few-layer transferred tungsten disulfide, and electrodeposited lithium, by measuring lattice strain from experimentally recorded electron diffraction data. The results show good agreement between experiment and simulation, providing confidence in the ability to transfer strain from holder to sample and relate TEM crystal structural observations with material mechanical properties.

2D materials↗
DOE OSTI2021

Physics-Informed Machine Learning Models for Predicting the Progress of Reactive-Mixing

This paper presents a physics-informed machine learning (ML) framework to construct reduced-order models (ROMs) for reactive-transport quantities of interest (QoIs) based on high-fidelity numerical simu-lations. QoIs include species decay, product yield, and degree of mixing. The ROMs for QoIs are applied to quantify and understand how the chemical species evolve over time. First, high-resolution datasets for constructing ROMs are generated by solving anisotropic reaction-di?usion equations using a non-negative finite element formulation for di?erent input parameters. The reactive-mixing model input parameters are: time-scale associated with flipping of velocity, spatial-scale controlling small/large vortex structures of velocity, perturbation parameter of the vortex-based velocity, anisotropic dispersion strength/contrast, and molecular diffusion. Second, random forests, F-test, and mutual information criterion are used to evaluate the importance of model inputs/features with respect to QoIs. We observed that anisotropic dispersion strength/contrast is the most important feature and time-scale associated with flipping of velocity is the least important feature. Third, Support Vector Machines (SVM) and Support Vector Regression (SVR) are used to construct ROMs based on the model inputs. The constructed SVR-ROMs are then used to predict scaling of QoIs. We also present estimates and inequalities on the QoIs, which inform that the species decay, mix, and produce in an exponential fashion. These inequalities also inform that a radial basis function is the most suitable kernel for the SVM/SVR models for QoIs. It is observed that R2-score for SVR-ROMs on unseen data is greater than 0.9, implying that the SVR-ROMs are able to predict the reaction-diffusion system state reasonably well. Finally, in terms of the computational cost, the proposed SVM-ROMs are O(107) times faster than running a high-fidelity finite element simulation for evaluating QoIs. This makes the proposed ML-based ROMs attractive for reactive-transport sensing and real-time monitoring applications as they are significantly faster yet reasonably accurate.

Mudunuru, Maruti K.↗
DOE OSTI2020

Global System Reduction Order Modeling for Localized Feature Inclusion

The development of reduced-order models remains an active research area, despite advances in computational resources. Here we develop a novel order-reduction approach that is designed to incorporate isolated regions that contain, for example, nonlinearitites or accumulating damage. The approach is designed to use global modes of the overall system response, which are then naturally coupled to the response in the isolated region of interest. Two examples are provided to demonstrate both the accuracy and the computational efficiency of the proposed approach. The performance of this approach is compared to the exact response corresponding to a finite element simulation for the chosen problems. In addition, the accuracy and computational efficiency are shown relative to a standard Galerkin reduction based on the linear normal modes. It is found that the proposed reduction offer computational efficiency comparable to a Galerkin reduction, but more accurately represents the response of the system when both are compared to the finite element simulation.

42 ENGINEERING↗
DOE OSTI2020

Iodine Effects on Medical Isotope Producer Delay & Guard Bed Temperature and Performance in Support of NA-23 Program

In support of NA-23’s objective to accelerate the establishment of a reliable commercial United States (U.S.) domestic supply of 99Mo produced without the use of highly enriched uranium (HEU), the Pacific Northwest National Laboratory (PNNL) was tasked with modeling the performance of activated carbon delay and guard beds that may be typical for non-HEU 99Mo production processes. The specific tasks accomplished by PNNL included calculations of long-term effects of iodine on delay bed performance and development of a finite element simulation to better understand delay bed performance under various scenarios specifically with radioxenon and radioiodine decay heating of the delay beds. The main issue of interest in this study is understanding the effect of iodine on the abatement performance of the beds if iodine is not sufficiently removed from gases prior to reaching the delay beds. Finite element simulations of delay bed performance were developed using COMSOL Multiphysics® Software. These simulations are based on partial differential equations that describe xenon and iodine adsorption, heat transfer and decay of radioxenon and radioiodine isotopes. Heat transfer in the activated carbon beds is simulated with a heat source tied to the decay of radioactive xenon and iodine species. These simulations were used to investigate delay bed temperature and hold up time using radioxenon and radioiodine activities and flowrates expected from a non-HEU production process. A bounding case scenario was simulated, where all the iodine was released into the delay beds due to insufficient guarding of the delay beds. In addition, preliminary investigations of the guard beds were conducted.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

Effect of Diphenyl Content on Viscoelasticity of Poly(dimethyl- co -diphenyl)siloxane Melt and Network

Polydimethylsiloxane (PDMS) is one of the most widely used polymeric materials for sealants, adhesives, lubricants, and thermal as well as electrical insulation. At low temperatures, however, PDMS is subject to crystallization that can cause deterioration in mechanical function. A common way to suppress such crystallization is through the incorporation of phenylsiloxane into the backbone of polysiloxane. Nevertheless, the introduction of phenyl components, even in small quantities, could potentially change the properties of the siloxane in a significant way. In this work, a series of mechanical tests and finite element simulations were performed to study the macroscale viscoelasticity of two poly(dimethyl-co-diphenyl)siloxane formulations in order to understand the effects of a few percent diphenyl contents on the viscoelasticity of the polysiloxane material. Here, we utilized the small-angle X-ray scattering to investigate the microscopic structures of the copolymers and broadband dielectric spectroscopy and rheology to probe the chain dynamics at the microscale. The results of these characterizations were used to inform the finite element simulations. We found that the degree of cross-linking does not significantly alter the microstructure but can profoundly affect the viscoelastic response of the copolymer networks. The corresponding hysteretic behavior is interpreted in terms of reptation-like motion and relaxation of the effective free chains in the cured polymer network. The relaxation of the copolymer chains is slowed significantly by even a small increase in the molar ratio of the diphenyl component.

36 MATERIALS SCIENCE↗
DOE OSTI2019

AC Resistance Reduction Using Orthogonal Air Gaps in High-Frequency Inductors: Preprint

This paper presents a relatively simple technique to reduce winding losses due to fringing fields in high-frequency inductors. In high-frequency power electronics, ac inductor winding losses are affected by skin and proximity effects, including uneven current distribution due to fringing magnetic fields around airgaps. It is well known how fringing effects can be mitigated using distributed airgaps, at the expense of non-standard core or winding geometry. In planar magnetics, the orthogonal-airgap approach proposed in this paper combines airgaps in core segments parallel with the windings with conventionally placed airgaps in segments perpendicular to the windings. The approach is developed using a 1D analytical framework and validated by 2D finite-element simulations. Analytical guidelines are presented to optimize the airgaps to achieve minimum ac resistance. As a case study, a planar inductor is designed for an 8~kW SiC-based buck converter operating at 250~kHz. It is shown how the orthogonal airgaps result in more than 45\% reduction in ac resistance and substantially reduced inductor losses compared to the design using standard airgaps. The results are verified by loss measurements on an experimental converter prototype.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Ultrasonic waves from radial mode excitation of a piezoelectric disc on the surface of an elastic solid

In this work, we discuss the ultrasonic wave propagation characteristics in an isotropic elastic solid due to the radial mode excitation of a piezoelectric disc actuator bonded to its surface. Finite element simulations using coupled electromechanical modeling is employed to investigate the wave propagation behavior. We find that the radial mode vibrations on the surface of the elastic solid generates all the three types of ultrasonic waves: longitudinal, shear, and surface waves. The waves in the solid are comprised of a central lobe and multiple side lobes based on the frequency of the radial mode excitation. The central lobe is predominantly composed of longitudinal waves and the side lobes are composed of shear waves. While the longitudinal waves have a strong central lobe and weak side lobes that are not fully developed, shear waves consists of fully developed side lobes. In addition, we observe that the longitudinal waves have fewer side lobes inside the solid compared to the shear waves. A semi-analytical approach is presented to explain the above observations. The interface traction boundary condition between the piezoelectric disc and the elastic solid is approximated as a truncated 'Bessel' excitation over the area of the piezoelectric disc in contact with the solid. A Fourier–Bessel series expansion technique is used to investigate the wave propagation behavior from the above traction boundary condition. The results obtained explain the observations from the finite element simulations with regard to the number and distribution of side lobes pertaining to the longitudinal and shear waves generated from the radial modes of the piezoelectric disc. The proposed semi-analytical approach is general and can be applied to any arbitrary axisymmetric excitations on the surface of an elastic solid.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High strain-rate compression behavior of polymeric rod and plate Kelvin lattice structures

The compressive high strain-rate behavior of polymeric Kelvin lattice structures with rod-based or plate-based unit cells was investigated through experimental techniques and finite element simulations. Polymeric lattice structures with 5x5x5 unit cell geometries were manufactured on the millimeter scale using vat polymerization additive manufacturing and tested at low (0.001/s) and high (1000/s) strain-rates. High strain-rate experiments were performed and validated for a viscoelastic split-Hopkinson (Kolsky) pressure bar system (SHPB) coupled with high-speed imaging and digital image correlation (DIC). Experimental results at both low and high strain-rates show the formation of a localized deformation band which was more prevalent in low relative density specimens and low strain-rate experiments. Strain-rate effects of lattice specimens strongly correlate with effects of the base polymer material; both bulk polymer and lattice specimen demonstrated strain-rate hardening, strain-rate stiffening, and decreased fracture strain under dynamic loading. Results show mechanical failure properties and energy absorption depended strongly on the relative density of the lattice specimen and exhibited distinct scaling between relative density and geometry type (rod, plate) and loading rate. High relative density plate-lattices demonstrated inferior mechanical properties to rod-lattices; however, there exists a critical relative density for a given mechanical property (17%- 28%) below which plate-lattices outperform rod-lattices of similar mass. As a result, high strain-rate explicit finite element simulations were performed and showed good agreement with the mechanical failure trends and deformation modes observed in the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Improving Computational Efficiency of Mechanical Finite Element Method Simulations for PV Modules: Preprint

In this work we have elucidated the trade off between structural mechanics FEM model accuracy and computation time by employing lower fidelity viscoelastic models. Results indicate that computation time can easily be cut in half while only expecting a potential maximum error of 10 % by considering lower fidelity models. A novel approach to produce the Prony Series fit for viscoelastic characterization has also been presented. Employing this approach, we were able to achieve a further 10 % reduction in computation time without further sacrificing simulation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

Identification Uncertainty in Inverse Material Model Parameter Determination: A Sensitivity‐Based Decision Process for Load Path Selection

This research proposes a sensitivity-based framework for selecting the optimal prescribed loading path for a biaxial cruciform specimen. Optimality here is determined by the direction and magnitude of the prescribed displacement that minimizes the influence of random noise on the material model parameter identification. Using simulated experimental data based on finite element simulation, in this work, we identify the material model parameters of a Ludwik hardening model and plane stress implementation of the Hill-48 yield criterion using finite element model updating (FEMU). Our analysis reveals that the identification (or estimator) uncertainty of model parameters depends on the displacement boundary conditions (i.e., loading sequence) and the ground-truth value of the individual parameters. Optimal experimental design (OED) criteria based on the Fisher information matrix were investigated to mitigate indecision in the choice of optimal load path when the identification uncertainty of different material model parameters optimized at different load paths. The determinant of the Fisher information matrix was chosen here as the more useful metric due to its ability to capture uncertainty of the most influential material model parameters. The proposed framework demonstrates potential for real-time automated load step selection using scalar criteria derived prior to mechanical loading. The framework can be generalized to other geometries, boundary conditions and material models, allowing this procedure to be utilized for different experimental configurations and materials.

Fayad, Samuel S. [University of Illinois at Urbana↗
DOE OSTI2025

Prelude to Simulation Diagnostics in MOOSE: Diagnosing Meshes

As finite element simulations increase in complexity, so does the need for high quality meshes. Without proper diagnostics, issues pertaining to poor mesh quality or improper mesh construction may go undetected, causing simulations to fail unexpectedly. This paper details the new mesh diagnostic system available in MOOSE, which provides users with the tools to inspect mesh quality for both externally and internally generated meshes. This system can detect issues arising from element volume size, non-conformal elements, intersecting edges, missing sidesets, as well as other common mesh issues. As meshes increase in size and complexity, such issues become more difficult to diagnose through visual inspection alone.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗
DOE OSTI2022

On Residual Stresses and Reference Temperatures in Thermomechanical Simulations of Photovoltaic Modules Using the Finite Element Method

Thermomechanical simulation of photovoltaic (PV) modules using the finite element method (FEM) is a useful tool to evaluate module design features in terms of structural integrity, reliability, and durability. One of the main challenges in the numerical modeling of a PV module is the incorporation of residual stresses induced by the manufacturing process. Modeling assumptions and abstractions are necessary to limit the model complexity and reduce the computational time. However, oversimplifications and incorrect assumptions can lead to erroneous numerical results. Unfortunately, much simulation work still neglects process-induced stresses. This can lead to incorrect predictions of the stress-strain history and erroneous conclusions during the design process. In this work, we review current modeling practices for incorporating process-induced stresses, and contrast numerical models that consider residual stresses with those that neglect them. We find that the simulation objective and available material properties dictate which process steps need to be modeled, and explore in depth the modeling of residual stresses induced by the lamination process. We demonstrate that a simplified cooldown procedure at the beginning of the simulation can increase the model accuracy and discuss appropriate choices for starting and reference temperatures in the finite element model.

14 SOLAR ENERGY↗
DOE OSTI2024

Fabrication and Performance Evaluation of Double-Sided Copper Nanowire Arrays as Thermal and Electrical Interfacial Layers

Reducing contact interface thermal and electrical resistances is in great demand across various industries, particularly in the semiconductor industry. This study introduces an approach using double-sided copper nanowire (Cu NW) arrays on copper sheets as both thermal and electrical interfacial layers, designed to effectively accommodate the topographical inconsistencies between contact surfaces. Further, experimental outcomes reveal a significant reduction in thermal contact resistance (TCR), with a value of 2.5 mm 2 K W –1 , thereby exceeding the efficiency of reported nanostructural thermal interface materials (TIMs). Additionally, when utilized as an electrical interfacial layer, these double-sided Cu NWs arrays dramatically reduced electrical contact resistance (ECR), outperforming traditional conductive grease in applications necessitating separable bonding, though showing comparable performance to costly silver-based conductive epoxies required for permanent, inseparable bonds. The promising results of the double-sided Cu NWs arrays in reducing both TCR and ECR, confirmed by finite element simulation, highlight their substantial potential in advancing TIMs and electrical interconnection applications across various sectors.

42 ENGINEERING↗
DOE OSTI2022

Performance of a MQXF Nb 3 Sn Quadrupole Magnet Under Different Stress Level

In a dipole or in a quadrupole accelerator magnet, the displacement of the coil turns induced by the electromagnetic forces can cause quenches limiting the magnet performance. For this reason, an azimuthal preload is applied to avoid azimuthal movements of the coil up to the required operational current. However, several tests showed that accelerator magnets can operate with a partial preload, i.e. that coil unloading during the ramp does not prevent reaching higher currents. This issue is particularly relevant for Nb 3 Sn magnets, where the loads applied to the Nb3Sn filaments can reach the degradation limits of critical current. In order to investigate the impact of coil preload on the quench performance, the MQXFS6 short model quadrupole for the High Luminosity Upgrade was tested under an azimuthal preload at 80% of the short sample current, reaching 93% of short sample current at 1.9 K. The preload was then released to 60%, still showing ability to operate in the range of 80-85% of short sample current as required by HL-LHC project. With this lower preload, the ability of going above 90% of short sample was lost, and a significant training appeared above 85%. When the preload was restored to the original 80% value, the magnet reached with few quenches 95% of short sample (13.4 T peak field). Magnetic measurements confirm the larger movement of the coil in the case with lower preload, and agree with finite element simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Thermodynamic theory of crystal plasticity: Formulation and application to polycrystal fcc copper

We present a thermodynamic description of crystal plasticity. Our formulation is based on the Langer-Bouchbinder-Lookman thermodynamic dislocation theory (TDT), which asserts the fundamental importance of an effective temperature that describes the state of configurational disorder and therefore the dislocation density of the crystalline material. We extend the TDT description from isotropic plasticity to crystal plasticity with many slip systems. Finite-element simulations show favourable comparison with experiments on polycrystal fcc copper under uniaxial compression, tension, and simple shear. Here, the thermodynamic theory of crystal plasticity thus provides a thermodynamically consistent and physically rigorous description of dislocation motion in crystals. We also discuss new insights about the interaction of dislocations belonging to different slip systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A high-frequency, low-power resonant radio-frequency neutron spin flipper for high-resolution spectroscopy

Here, we present a resonant-mode, transverse-field, radio-frequency (rf) neutron spin flipper design that uses high-temperature superconducting films to ensure sharp transitions between uniform magnetic field regions. Resonant mode allows for low-power, high-frequency operation but requires strict homogeneity of the magnetic fields inside the device. This design was found to efficiently flip neutrons at 96.6 ± 0.6% at an effective frequency of 4 MHz in bootstrap configuration with a beam size of 2.4 × 2.5 cm 2 and a wavelength of 0.4 nm. The high frequency and efficiency enable this device to perform high-resolution neutron spectroscopy with comparable performance with currently implemented rf flipper designs. The limitation of the maximum frequency was found due to the field homogeneity of the device. We numerically analyze the maximum possible efficiency of this design using a Bloch solver simulation with magnetic fields generated from finite-element simulations. We also discuss future improvements of the efficiency and frequency to the design based on the experimental and simulation results.

47 OTHER INSTRUMENTATION↗