Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Evolution of A-Type Macrosegregation in Large Size Steel Ingot After Multistep Forging and Heat Treatment

A-type macrosegregation refers to the channel chemical heterogeneities that can be formed during solidification in large size steel ingots. In this research, a combination of experiment and simulation was used to study the influence of open die forging parameters on the evolution of A-type macrosegregation patterns during a multistep forging of a 40 metric ton (MT) cast, high-strength steel ingot. Macrosegregation patterns were determined experimentally by macroetch along the longitudinal axis of the forged and heat-treated ingot. Mass spectroscopy, on more than 900 samples, was used to determine the chemical composition map of the entire longitudinal sectioned surface. FORGE NxT 1.1 finite element modeling code was used to predict the effect of forging sequences on the morphology evolution of A-type macrosegregation patterns. For this purpose, grain flow variables were defined and implemented in a large scale finite element modeling code to describe oriented grains and A-type segregation patterns. Examination of the A-type macrosegregation showed four to five parallel continuous channels located nearly symmetrical to the axis of the forged ingot. In some regions, the A-type patterns became curved or obtained a wavy form in contrast to their straight shape in the as-cast state. Mass spectrometry analysis of the main alloying elements (C, Mn, Ni, Cr, Mo, Cu, P, and S) revealed that carbon, manganese, and chromium were the most segregated alloying elements in A-type macrosegregation patterns. The observed differences were analyzed using thermodynamic calculations, which indicated that changes in the chemical composition of the liquid metal can affect the primary solidification mode and the segregation intensity of the alloying elements. Finite element modeling simulation results showed very good agreement with the experimental observations, thereby allowing for the quantification of the influence of temperature and deformation on the evolution of the shape of the macrosegregation channels during the open die forging process.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gyrokinetic electromagnetic particle simulations in triangular meshes with C1 finite elements

The triangular mesh-based gyrokinetic scheme enables comprehensive axis-to-edge studies across the entire plasma volume. Our approach employs triangular finite elements with first-derivative continuity (C1), building on previous work to facilitate gyrokinetic simulations. Additionally, we have adopted the mixed variable/pullback scheme for gyrokinetic electromagnetic particle simulations. The filter-free treatment in the poloidal cross-section with triangular meshes introduces unique features and challenges compared to previous treatments using structured meshes. Our implementation has been validated through benchmarks using ITPA-toroidicity-induced Alfvén eigenmode parameters, showing its capability in moderate to small electron skin depth regimes. Additional examinations using experimental parameters confirm its applicability to realistic plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Material Variability Effects on Damage Development Within Composite Adhesively Bonded Joints

Use of composite adhesively bonded joints (ABJ) is of critical importance to the adoption of composite materials in the automotive industry, as ABJ enable lower stress concentrations as compared to conventional mechanically fastened joints. ABJ are better suited for joining composite materials as compared to fastened joints because fastened joints require drilling of holes which locally affect composite material structure. Composite materials and adhesives are subject to unavoidable stochastic local material variations which make different failure scenarios possible. An experimentally tested ABJ configuration is simulated using finite element analysis (FEA). Experimentally, under tension, the joints failed by three major failure modes, with peak loads ranging from 13.0-16.1 kips. Progressive failure analysis tools are used to simulate damage development within each material within the joint. The simulation agreed well with the average experimental peak load. Stochastically occurring adhesive porosity and matrix-fiber micro-disbonding were numerically simulated. The simulations revealed a similar trend as observed experimentally: joints which failed at higher peak loads had lower levels of damage within the face-sheets of the composite panels which were adhesively bonded; these joints which failed at higher peak loads also had greater damage in the doubler of the experimentally tested double lap joint configuration.

Richard Larson↗
DOE OSTI2015

Multi-element least square HDMR methods and their applications for stochastic multiscale model reduction

Stochastic multiscale modeling has become a necessary approach to quantify uncertainty and characterize multiscale phenomena for many practical problems such as flows in stochastic porous media. The numerical treatment of the stochastic multiscale models can be very challengeable as the existence of complex uncertainty and multiple physical scales in the models. To efficiently take care of the difficulty, we construct a computational reduced model. To this end, we propose a multi-element least square high-dimensional model representation (HDMR) method, through which the random domain is adaptively decomposed into a few subdomains, and a local least square HDMR is constructed in each subdomain. These local HDMRs are represented by a finite number of orthogonal basis functions defined in low-dimensional random spaces. The coefficients in the local HDMRs are determined using least square methods. We paste all the local HDMR approximations together to form a global HDMR approximation. To further reduce computational cost, we present a multi-element reduced least-square HDMR, which improves both efficiency and approximation accuracy in certain conditions. To effectively treat heterogeneity properties and multiscale features in the models, we integrate multiscale finite element methods with multi-element least-square HDMR for stochastic multiscale model reduction. This approach significantly reduces the original model's complexity in both the resolution of the physical space and the high-dimensional stochastic space. We analyze the proposed approach, and provide a set of numerical experiments to demonstrate the performance of the presented model reduction techniques. - Highlights: • Multi-element least square HDMR is proposed to treat stochastic models. • Random domain is adaptively decomposed into some subdomains to obtain adaptive multi-element HDMR. • Least-square reduced HDMR is proposed to enhance computation efficiency and approximation accuracy in certain conditions. • Integrating MsFEM and multi-element least square HDMR can significantly reduce computation complexity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2019

A Multi-Physics Study on High-Specific Power Li-O2 Batteries for Electric Aircraft

Commercialization of lithium-air batteries faces many challenges, such as electrolyte decomposition, short cycle life, low energy and power density, etc. However, commercialization of Li-O2 batteries for aeronautics is much more challenging due to additional safety constraints on cyclability and performance (high specific power and specific energy). For this presentation, we will discuss inter-related aspects of physics-based modeling of a pack: cell and battery model calibration. In addition, we will evaluate and present optimal battery designs for high discharge current density, high discharge time, and low battery mass using simulation-based optimization.The Finite Element Model (FEM) used to simulate a Li-O2 cell is based on the work of Bevara [1]. The different aspects of the model are based on: porous electrode theory and concentrated electrolyte theory; quantum tunneling model for the resistance of conformal layer of discharge product (Li2O2) [1]; Butler-Volmer kinetics for electrochemical reaction; Fick's diffusion for oxygen transport; and an oxygen dissolution model is applied at the air/electrolyte interface [2]. The electrolyte properties such as ion conductivity, ion diffusion, oxygen diffusion, and mass density of the electrolyte were taken from Molecular Dynamics (MD) simulations [3]; while the other model parameters, which includes mass of cell components, were calibrated to match experiments at high discharge current densities. The cell mass includes the anode, cathode, separator, electrolyte, and other components (such as current collector). This calibrated model is used to perform parametric studies on cathode thickness, porosity, tortuosity, carbon particle size, electrolyte transport and material properties, partial pressure of oxygen, discharge time, and discharge current density to study optimal designs for high specific power and energy. References:1. Bevara, V. & Andrei, P. (2014), J. Electrochem. Soc. 161 (14), A2068-A2079.2.Mehta, M. & Andrei, P. (2015), J. Power Sources. 286, 299-308.3.Liyana-Arachchi, T.; Haskins, J.; Burke, C.; Diederichsen, K.; McCloskey, B.; & Lawson, J. (2018), J. Phys. Chem. B. 122 (36), 8548 - 8559.4.Choi, W.; Kikumoto, H.; Choudhary, R. & Ooka, R. (2018), Applied Energy, 209, 306-321.

Mehta, Mohit↗
NASA NTRS2002

Comparisons of Kinematics and Dynamics Simulation Software Tools

Kinematic and dynamic analyses for moving bodies are essential to system engineers and designers in the process of design and validations. 3D visualization and motion simulation plus finite element analysis (FEA) give engineers a better way to present ideas and results. Marshall Space Flight Center (MSFC) system engineering researchers are currently using IGRIP from DELMIA Inc. as a kinematic simulation tool for discrete bodies motion simulations. Although IGRIP is an excellent tool for kinematic simulation with some dynamic analysis capabilities in robotic control, explorations of other alternatives with more powerful dynamic analysis and FEA capabilities are necessary. Kinematics analysis will only examine the displacement, velocity, and acceleration of the mechanism without considering effects from masses of components. With dynamic analysis and FEA, effects such as the forces or torques at the joint due to mass and inertia of components can be identified. With keen market competition, ALGOR Mechanical Event Simulation (MES), MSC visualNastran 4D, Unigraphics Motion+, and Pro/MECHANICA were chosen for explorations. In this study, comparisons between software tools were presented in terms of following categories: graphical user interface (GUI), import capability, tutorial availability, ease of use, kinematic simulation capability, dynamic simulation capability, FEA capability, graphical output, technical support, and cost. Propulsion Test Article (PTA) with Fastrac engine model exported from IGRIP and an office chair mechanism were used as examples for simulations.

Yeu-Sheng Paul Shiue↗
NASA NTRS2022

Laser Patterning Assisted Devitrification and Domain Engineering of Amorphous and Nanocrystalline Alloys

We have investigated laser-patterning and devitrification as a novel method to realize emergent properties in modified regions at the ribbon surface consisting of periodic localized laser spots thermally treated under rapid heating and cooling conditions. Fe-based amorphous ribbons were annealed for systematically varied laser power. Extremely high heating and cooling rates on the order of 10(exp 8)-10(exp 10)K/s have been estimated by finite element analysis simulations. Observations suggest surface melting followed by rapid solidification fast enough to quench and form an amorphous structure upon cooling. Diffusion of heat occurs from the laser irradiated surface region so that a finite volume of material surrounding this region rises in temperature above the crystallization temperature relevant for conventional isothermal annealing experiments. The underlying mechanism of laser annealing and consequences of heat transients involved are discussed in terms of impacts upon micro/nanostructure, residual stresses, and magnetic domain structure surrounding the laser irradiated region. The study illustrates the potential to exploit spatially optimized phase transformations in a scalable manufacturing process of amorphous and nanocrystalline alloys to locally access otherwise inaccessible extreme heating and cooling rates.

Amorphous ribbons↗
NASA NTRS1993

A class of hybrid finite element methods for electromagnetics: A review

Integral equation methods have generally been the workhorse for antenna and scattering computations. In the case of antennas, they continue to be the prominent computational approach, but for scattering applications the requirement for large-scale computations has turned researchers' attention to near neighbor methods such as the finite element method, which has low O(N) storage requirements and is readily adaptable in modeling complex geometrical features and material inhomogeneities. In this paper, we review three hybrid finite element methods for simulating composite scatterers, conformal microstrip antennas, and finite periodic arrays. Specifically, we discuss the finite element method and its application to electromagnetic problems when combined with the boundary integral, absorbing boundary conditions, and artificial absorbers for terminating the mesh. Particular attention is given to large-scale simulations, methods, and solvers for achieving low memory requirements and code performance on parallel computing architectures.

Volakis, J. L.↗
DOE OSTI2016

Theoretical study on the top- and enclosed-contacted single-layer MoS{sub 2} piezotronic transistors

Recently, the piezotronic effect has been observed in two-dimensional single-layer MoS{sub 2} materials, which have potential applications in force and pressure triggered or controlled electronic devices, sensors, and human-machine interfaces. However, classical theory faces the difficulty in explaining the mechanism of the piezotronic effect for the top- and enclosed-contacted MoS{sub 2} transistors, since the piezoelectric charges are assumed to exist only at the edge of the MoS{sub 2} flake that is far from the electronic transport pathway. In the present study, we identify the piezoelectric charges at the MoS{sub 2}/metal-MoS{sub 2} interface by employing both the density functional theory and finite element method simulations. This interface is on the transport pathway of both top- and enclosed-contacted MoS{sub 2} transistors, thus it is capable of controlling their transport properties. This study deepens the understanding of piezotronic effect and provides guidance for the design of two-dimensional piezotronic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Mechanical properties of niobium radio-frequency cavities

Radio-frequency cavities made of bulk niobium are one of the components used in modern particle accelerators. The mechanical stability is an important aspect of cavity design, which typically relies on finite-element analysis simulations using material properties from tensile tests on sample. This contribution presents the results of strain and resonant frequency measurements as a function of a uniform pressure up to 722 kPa, applied to single-cell niobium cavities with different crystallographic structure, purity and treatments. In addition, burst tests of high-purity multi-cell cavities with different crystallographic structure have been conducted up to the tensile strength of the material. Finite-element analysis of the single-cell cavity geometry is in good agreement with the observed behavior in the elastic regime assuming a Young's modulus value of 88.5 GPa and a Poisson's ratio of 0.4, regardless of crystallographic structure, purity or treatment. However, the measured yield strength and tensile strength depend on crystallographic structure, material purity and treatment. In particular, the results from this study show that the mechanical properties of niobium cavities with large crystals are comparable to those of cavities made of fine-grain niobium.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulated effects of sample size and grain neighborhood on the modeling of extreme value fatigue response

Assessing the size of representative volume elements (RVEs) for fatigue-related applications is challenging. A RVE relevant to random microstructure requires a volume of material that is sufficiently large to capture the grain/phase heterogeneity that captures all statistical moments of the distribution of the driving force for fatigue crack formation at “hot spot” grains. Consequently, the large size of a microstructure RVE required to study fatigue phenomena is largely computationally intractable and difficult to explore. A more realistic objective in this work is to systematically study, as a function of the size of a statistical sample of microstructure, trends towards convergence of the simulated distribution of driving force for fatigue crack formation. Our present work accordingly leverages the recently developed open-source PRISMS-Fatigue framework to examine the trends in convergence of extreme value distributions (EVD) of Fatigue Indicator Parameters (FIPs) in progressively larger polycrystalline microstructure realizations of FCC Al alloy 7075-T6 using crystal plasticity finite element method simulations. The results are compared to the traditional method in which ensembles of statistical volume elements (SVEs) are simulated to build up statistics intended to approximate those associated with a larger volume of material. The convergence of EVDs with increase of size of a SVE of microstructure is closely related to the extent of grain nearest neighbor (NN) interactions. Accordingly, the sensitivity of the local micromechanical response at hot spot grains is quantitatively investigated by systematically varying the orientations of NN grains. Results indicate that SVEs with cubic crystallographic texture tend towards convergence of the EVD of FIPs with tens of thousands of grains while the random and rolled textures require larger volumes. Simple relationships based on microstructure parameters (e.g., Schmid Factor, grain size, NN misorientation) do not completely correlate to fatigue hot spot grains. Finally, the sensitivity of the extreme value fatigue response at hot spot grains extends to the 3rd NN when a single neighborhood grain orientation is altered.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermal Analysis of a 50 kW Three-Phase Wireless Charging System

In this paper, the thermal analysis of a 50 kW three-phase wireless charging system (WCS) is presented. Addressing the thermal challenge is essential for designing a compact charging pad for 50 kW and higher-power WCS pads. Low thermal conductivity of the Litz wire and ferrite, and uneven distribution of the coil, and core loss cause the high temperature in the charging pads. In this paper, the loss distribution of a 50 kW three-phase WCS is investigated, and the temperature distribution is simulated using finite element analysis (FEA) considering the magnetic and non-magnetic materials of a charging pad. The thermal characteristics of an extremely high-power density 50 kW WCS prototype are tested experimentally for 10 minutes of operation. The simulation and experimental results show that the coil temperature increases to 65°C, the core temperature varies between 50°C to 150°C, and the packaging temperature increases to 65°C after 10 minutes of operation at rated 50 kW output power.

Mohammad, Mostak↗
DOE OSTI2016

GaAs/Ge crystals grown on Si substrates patterned down to the micron scale

Monolithic integration of III-V compounds into high density Si integrated circuits is a key technological challenge for the next generation of optoelectronic devices. In this work, we report on the metal organic vapor phase epitaxy growth of strain-free GaAs crystals on Si substrates patterned down to the micron scale. The differences in thermal expansion coefficient and lattice parameter are adapted by a 2-μm-thick intermediate Ge layer grown by low-energy plasma enhanced chemical vapor deposition. The GaAs crystals evolve during growth towards a pyramidal shape, with lateral facets composed of (111) planes and an apex formed by (137) and (001) surfaces. The influence of the anisotropic GaAs growth kinetics on the final morphology is highlighted by means of scanning and transmission electron microscopy measurements. The effect of the Si pattern geometry, substrate orientation, and crystal aspect ratio on the GaAs structural properties was investigated by means of high resolution X-ray diffraction. The thermal strain relaxation process of GaAs crystals with different aspect ratio is discussed within the framework of linear elasticity theory by Finite Element Method simulations based on realistic geometries extracted from cross-sectional scanning electron microscopy images.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Integrating Experiments, Simulations, and Artificial Intelligence to Accelerate the Discovery of High-Performance Green Composites

The imperative for incorporating greener materials into the aerospace industry necessitates addressing significant challenges associated with the microstructural variability exhibited by recycled and sustainable feedstocks. In this study, we propose an integrated methodology that combines experimental investigations, finite element analysis, and artificial intelligence to develop sustainable composites with consistent properties. Our approach utilizes a pipeline comprising an automated mechanical tester, a finite element method simulator, and a convolutional neural network predictor to identify and optimize fabrication parameters for achieving desired mechanical characteristics in composites. By employing a nested-loop pipeline, our methodology improves sample efficiency, accuracy, and effectively bridges the gap between simulations and real-world performance. This unique methodology offers a promising avenue for facilitating the adoption of aerospace-appropriate green composites.

Athanasiou, Christos↗
DOE OSTI2022

Simulation for the fault current limiting operation of REBCO CORC superconducting cables with different core materials

With high power density and low loss, CORC superconducting cables composed of RE- Ba 2 Cu 3 O 7–δ (REBCO) coated conductors are of great interest for power transmission applications. They can also effectively be designed to serve as fault current limiters, thanks to their sharp superconducting-to-normal transition. Coupled electromagnetic-thermal finite element (FE) simulations implemented in COMSOL Multiphysics were developed and validated against the published results for the fault current limiting performance of two CORC cables of REBCO coasted conductors wound on either copper or stainless steel cores. For improved accuracy, temperature dependencies of electrical and thermal properties of all component materials were considered in the simulations. The simulations were performed in the cross-sections of the cables for every individual layer of each REBCO tape to deliver comprehensive understanding of evolution of the current distribution and temperature rise in those layers. Furthermore, the studies can suggest approaches to optimize cable design for fault current limiting (FCL) applications by assessing the role of individual components. In the paper, possible computational errors, challenges and improvements are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

A hybrid approach for fusing 4D-MRI temporal information with 3D-CT for the study of lung and lung tumor motion

Purpose: Accurate visualization of lung motion is important in many clinical applications, such as radiotherapy of lung cancer. Advancement in imaging modalities [e.g., computed tomography (CT) and MRI] has allowed dynamic imaging of lung and lung tumor motion. However, each imaging modality has its advantages and disadvantages. The study presented in this paper aims at generating synthetic 4D-CT dataset for lung cancer patients by combining both continuous three-dimensional (3D) motion captured by 4D-MRI and the high spatial resolution captured by CT using the authors’ proposed approach. Methods: A novel hybrid approach based on deformable image registration (DIR) and finite element method simulation was developed to fuse a static 3D-CT volume (acquired under breath-hold) and the 3D motion information extracted from 4D-MRI dataset, creating a synthetic 4D-CT dataset. Results: The study focuses on imaging of lung and lung tumor. Comparing the synthetic 4D-CT dataset with the acquired 4D-CT dataset of six lung cancer patients based on 420 landmarks, accurate results (average error <2 mm) were achieved using the authors’ proposed approach. Their hybrid approach achieved a 40% error reduction (based on landmarks assessment) over using only DIR techniques. Conclusions: The synthetic 4D-CT dataset generated has high spatial resolution, has excellent lung details, and is able to show movement of lung and lung tumor over multiple breathing cycles.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Newmark local time stepping on high-performance computing architectures

In multi-scale complex media, finite element meshes often require areas of local refinement, creating small elements that can dramatically reduce the global time-step for wave-propagation problems due to the CFL condition. Local time stepping (LTS) algorithms allow an explicit time-stepping scheme to adapt the time-step to the element size, allowing near-optimal time-steps everywhere in the mesh. We develop an efficient multilevel LTS-Newmark scheme and implement it in a widely used continuous finite element seismic wave-propagation package. In particular, we extend the standard LTS formulation with adaptations to continuous finite element methods that can be implemented very efficiently with very strong element-size contrasts (more than 100x). Capable of running on large CPU and GPU clusters, we present both synthetic validation examples and large scale, realistic application examples to demonstrate the performance and applicability of the method and implementation on thousands of CPU cores and hundreds of GPUs.

97 MATHEMATICS AND COMPUTING↗