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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1980

An approach to simulation of the finite element machine

The paper describes the simulation approach for the finite element machine (FEM), a special-purpose computer for structural analysis calculations under development at the NASA-Langley Research Center. The FEM consists of an array (1000 or more) of general-purpose microcomputers performing structural analysis in an asychronous parallel manner. A simulator of sufficient flexibility to model the behavior of the FEM on many levels has been designed and coded, and verification has begun using the experimental FEM hardware. The structure of the simulator is described, and preliminary simulation results are presented.

Sawyer, P. L.↗

DOE OSTI2018

Microstructural study and simulation of intrinsic two-way shape memory behavior of functionally graded Ni-rich/NiTiCu thin film

Highlights: • Ni-rich/Ti-rich bi-layer thin film was deposited on Si by using DC magnetron sputtering. • The bi-layer exhibits the gradient structure and mechanical properties through the thickness. • The bi-layer shows intrinsic two-way shape memory effect due to functionally graded structure. • The intrinsic two-way shape memory effect was confirmed by finite element simulation. - Abstract: In this research, the gradient structure and mechanical properties of an austenitic/martensitic (Ni{sub 50.8}Ti/Ni{sub 45}TiCu{sub 5} (at.%)) bi-layer thin film was investigated. The bi-layer was deposited on a Si substrate using DC magnetron sputtering. After crystallization of the film at 773 K for 60 min, the microstructure and mechanical properties of the bi-layer thin film were characterized using X-ray diffraction (XRD), transmission electron microscopy (STEM and HRTEM) and nanoindentation, respectively. The diffraction pattern illustrated that the crystallized bi-layer was combined of martensitic and austenitic layers while secondary ion mass spectroscopy (SIMS) and high resolution transmission electron microscopy analysis demonstrated the existence of a compositional gradient through the thickness and residual strain in the interface of the bi-layer, respectively. Furthermore, the compositional gradient in the bi-layer led to gradual variations in the structure, hardness and Young's modulus through the thickness of the bi-layer. The intrinsic two-way shape memory effect due to its functionally graded structure and the induced stress was confirmed by experimental test and finite element simulation.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus

A particle-in-cell algorithm is derived with a canonical Poisson structure in the formalism of finite element exterior calculus. The resulting method belongs to the class of gauge-compatible splitting algorithms, which exactly preserve gauge symmetries and their associated conservation laws via the momentum map. We numerically demonstrate this time invariance of the momentum map and its usefulness in establishing precise initial conditions with a desired initial electric field and fixed background charge. The restriction of this canonical, finite element Poisson structure to the 1X2P $1\frac {1}{2}$ -dimensional phase space is also considered and simulated numerically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Instrumented Photovoltaic Modules for Environmental Characterization and Simulation Model Validation

Photovoltaic modules are subjected to various mechanical stressors in their deployment environments, ranging from installation handling to wind and snow loads. Damage incurred during these mechanical events has the potential to initiate subsequent degradation mechanisms, reducing useful module lifespan. Thus, characterizing the mechanical state of photovoltaic modules is pertinent to the development of reliable packaging designs. In this work, photovoltaic modules with strain gauges directly incorporated into the module laminate were fabricated and subjected to mechanical loading to characterize internal strains within the module when under load. These experimental measurements were then compared against results obtained by high-fidelity finite-element simulations. The simulation results showed reasonable agreement in the strain values over time; however, there were large discrepancies in the magnitudes of these strains. Both the instrumentation technique and the finite-element simulations have areas where they can improve. These areas of improvement have been documented. Despite the observed discrepancies between the experimental and simulated results, the module instrumentation proved to be a useful gauge in monitoring and characterizing the mechanical state. With some process improvements, this method could potentially be applied to other environments that a photovoltaic module will encounter in its lifetime that are known to cause damage and degrade performance.

14 SOLAR ENERGY↗

DOE OSTI2014

Pressure-assisted fabrication of organic light emitting diodes with MoO{sub 3} hole-injection layer materials

In this study, pressures of ∼5 to ∼8 MPa were applied to organic light emitting diodes containing either evaporated molybdenum trioxide (MoO{sub 3}) or spin-coated poly(3,4-ethylene dioxythiophene) doped with poly(styrene sulphonate) (PEDOT:PSS) hole-injection layers (HILs). The threshold voltages for both devices were reduced by about half, after the application of pressure. Furthermore, in an effort to understand the effects of pressure treatment, finite element simulations were used to study the evolution of surface contact between the HIL and emissive layer (EML) under pressure. The blister area due to interfacial impurities was also calculated. This was shown to reduce by about half, when the applied pressures were between ∼5 and 8 MPa. The finite element simulations used Young's modulus measurements of MoO{sub 3} that were measured using the nanoindentation technique. They also incorporated measurements of the adhesion energy between the HIL and EML (measured by force microscopy during atomic force microscopy). Within a fracture mechanics framework, the implications of the results are then discussed for the pressure-assisted fabrication of robust organic electronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Inverse design of hypoeutectoid pearlite steel microstructures using a deep learning and genetic algorithm optimization framework

Goal-oriented microstructure design in metallic materials is a challenging task due to complex structure-property relationships. Traditional experimental and computational approaches are time-intensive and economically inefficient, limiting their applicability for large-scale design space exploration. Here, in this work, we propose an end-to-end framework that integrates deep learning models with genetic optimization to design microstructures with targeted mechanical properties. Deep learning models enable accurate forward design, while their integration with genetic optimization enables efficient inverse design within a few hours, compared to days or weeks using conventional finite element simulations. The framework combines experimental characterization and finite element modeling to analyze the influence of microstructural features on the mechanical behavior of hypoeutectoid steels. Data from both experiments and simulations are used to train the deep learning models. To demonstrate its effectiveness, we apply the framework to 0.63% carbon steel with proeutectoid ferrite and pearlite phases, commonly used in industrial applications. In this study, 2D microstructures were used for modeling, selected primarily for computational efficiency and to establish proof of concept. The framework successfully optimizes microstructures for targeted yield strength, ultimate strength, and stress concentration factors while significantly reducing computational time. Beyond hypoeutectoid steels, this scalable framework can be extended to other material systems and integrated with additive manufacturing, offering an efficient approach for accelerating microstructure design for specific engineering applications.

ConvLSTM↗

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↗

NASA NTRS1991

FIESTA ROC: A new finite element analysis program for solar cell simulation

The Finite Element Semiconductor Three-dimensional Analyzer by Ralph O. Clark (FIESTA ROC) is a computational tool for investigating in detail the performance of arbitrary solar cell structures. As its name indicates, it uses the finite element technique to solve the fundamental semiconductor equations in the cell. It may be used for predicting the performance (thereby dictating the design parameters) of a proposed cell or for investigating the limiting factors in an established design.

Clark, Ralph O.↗

DOE OSTI2019

Dynamic mechanical behavior of lithium-ion pouch cells subjected to high-velocity impact

With broad-scale deployment of lithium-ion batteries, mechanical-abuse induced failure events such as short circuits and thermal runaway have made safety a focus of attention in battery design. Most previous studies have focused on quasi-static and low-velocity impact analyses, which could not represent the mechanical behavior of lithium-ion batteries in realistic situations. Here we report a combined experimental and computational study on the dynamic response of lithium-ion pouch cells subjected to high-velocity (200 m/s~1000 m/s) impact. Dynamic finite element simulations were performed to study the effects of internal interfacial behavior and external loading and boundary conditions on the dynamic mechanical behavior of pouch cells. Specifically, mechanical resistance across interfaces between different cell components have a pronounced impact on energy distribution and absorption, thereby determining the evolution of subsequent failure events. These results have practical implications for designing cells with liquid versus gel-type electrolytes, for example. Furthermore, finite element simulations on different external loading and boundary conditions indicate that a direct interaction between the impactor and the cell boundaries lead to a significant change in residual velocity, while the response of the pouch cell is independent of the prescribed far-field boundary conditions. In other words, there is not a significant difference in the response of pouch-format cells that have stacked electrodes versus wound electrodes at such high strain rates. These results are significantly different from quasi-static simulations wherein the evolution of failure events follows a very different pathway. In addition, the impact velocity profiles under different initial impact velocities offer the ballistic limit for the pouch cell and inform subsequent propagation analyses at the battery module level. Collectively, the findings presented here not only offer new perspectives on the role of the interfacial resistances and external boundary conditions on the mechanical safety of lithium-ion pouch cells, but also have implications for further investigations on cell design, modeling approaches, and test method development to study safety of lithium-ion batteries under extreme loading conditions.

25 ENERGY STORAGE↗

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↗

NASA NTRSDate not supplied

Evaluation of Human Spaceflight-Related Tissue Weight Relief Using Whole Body Finite Element Model Simulations

Tissue Weight Relief (TWR) is a physiological condition observed in human spaceflight. It not only impacts the injury biomechanics of soft tissue but also the physiological responses of the cardiovascular system both due to fluid redistribution and the effect of tissue-related transmural pressure on the large venous blood vessels. Understanding the effects of tissue weight relief is especially important because of the role it may play in understanding the cause of Space Associated Neuro-Ocular Syndrome (SANS). SANS can be characterized by a number of ocular changes which reduce visual acuity and SANS related symptoms occur in up to 51% of astronauts. A prevailing theory for the causation of SANS is that of headward (cephalad) fluid shift and a prolonged increase of Intracranial Pressure (ICP) similar to intracranial hypertension, which is not fully supported by the experimental data or astronaut symptom reporting. However, it is still believed that SANS is caused by a pressure change in the eye and the surrounding tissues. It has been proposed that TWR plays a substantial role in affecting internal pressures and fluid shifts in microgravity. In this effort, two whole-body Finite Element (FE) models – Elemance and THUMS – are used to ascertain the microgravity-associated TWR of the musculature surrounding the lower body veins. Elemance and THUMS are physics-based computational models that have been validated and verified for several automotive and domestic applications, and as such, can simulate the relief of soft tissue weight due to changes in the gravitation vector. Specifically, the current effort modeled the transition of the gravitational vector from 1G to 0G, applied across the whole-body model in a supine position. For each 1G to 0G transition simulation, the lower body vein’s transmural pressure-time profile was extracted and averaged around the anterior portion of the thigh muscle. The ascertained transmural pressure changes from 1G to 0G transition are given in Figure 1 for the Elemance and the THUMS FE models. The transmural pressure changes of 10 mmHg and 21 mmHg are in the same order of magnitude as Lu’s value of 44 mmHg. It is to be noted that Lu implemented a 0D to 1D lumped parameter model and the Elemance and THUMS are 3D higher order computational models. This proof-of-concept approach demonstrates that TWR pressure can be adequately estimated with in silico techniques however, further in silico investigations need to be conducted to address the unique contributions to the transmural pressure from each of the computational models.

Finite element modeling↗

DOE OSTI2015

The dynamics of Al/Pt reactive multilayer ignition via pulsed-laser irradiation

Reactive multilayers consisting of alternating layers of Al and Pt were irradiated by single laser pulses ranging from 100 μs to 100 ms in duration, resulting in the initiation of rapid, self-propagating reactions. The threshold intensities for ignition vary with the focused laser beam diameter, bilayer thickness, and pulse length and are affected by solid state reactions and conduction of heat away from the irradiated regions. High-speed photography was used to observe ignition dynamics during irradiation and elucidate the effects of heat transfer into a multilayer foil. For an increasing laser pulse length, the ignition process transitioned from a more uniform to a less uniform temperature profile within the laser-heated zone. A more uniform temperature profile is attributed to rapid heating rates and heat localization for shorter laser pulses, and a less uniform temperature profile is due to slower heating of reactants and conduction during irradiation by longer laser pulses. Finite element simulations of laser heating using measured threshold intensities indicate that micron-scale ignition of Al/Pt occurs at low temperatures, below the melting point of both reactants.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Direct evidence of strain transfer for InAs island growth on compliant Si substrates

Semiconductor heteroepitaxy on top of thin compliant layers has been explored as a path to make inorganic electronics mechanically flexible as well as to integrate materials that cannot be grown directly on rigid substrates. Here, we show direct evidences of strain transfer for InAs islands on freestanding Si thin films (7 nm). Synchrotron X-ray diffraction measurements using a beam size of 300 × 700 nm{sup 2} can directly probe the strain status of the compliant substrate underneath deposited islands. Using a recently developed diffraction mapping technique, three-dimensional reciprocal space maps were reconstructed around the Si (004) peak for specific illuminated positions of the sample. The strain retrieved was analyzed using continuous elasticity theory via Finite-element simulations. The comparison of experiment and simulations yields the amount of strain from the InAs islands, which is transferred to the compliant Si thin film.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2015

Deformation mechanisms in a precipitation-strengthened ferritic super alloy revealed by in situ neutron dffraction studies at elevated temperatures

The ferritic superalloy Fe–10Ni–6.5Al–10Cr–3.4Mo strengthened by ordered (Ni,Fe)AlB2-type precipitates is a candidate material for ultra-supercritical steam turbine applications above 923 K. Despite earlier success in improving its room-temperature ductility, the creep resistance of this material at high temperatures needs to be further improved, which requires a fundamental understanding of the high-temperature deformation mechanisms at the scales of individual phases and grains. In situ neutron diffraction has been utilized to investigate the lattice strain evolution and the microscopic load-sharing mechanisms during tensile deformation of this ferritic superalloy at elevated temperatures. Finite-element simulations based on the crystal plasticity theory are employed and compared with the experimental results, both qualitatively and quantitatively. Based on these interphase and intergranular load-partitioning studies, it is found that the deformation mechanisms change from dislocation slip to those related to dislocation climb, diffusional flow and possibly grain boundary sliding, below and above 873 K, respectively. Insights into microstructural design for enhancing creep resistance are also discussed.

nuclear (including radiation effects), defects, me↗