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 217 records · Page 12

Enhanced Inductance and Winding Loss Model for Coupled Inductors

An enhanced model for coupled inductors with significant fringing effect is developed in the form of equivalent circuit. The equivalent circuit is derived from a physical model that captures the flux paths through a leakage inductor and two mutual inductors on the primary and secondary side. A mutual winding resistor is added in parallel with the mutual inductor to model the winding loss affected by the fringing flux that penetrates the winding. The equivalent circuit is verified by finite-element simulation (FES) and successfully predicts the open-circuit winding loss and winding loss variation with phase-shift impact.

Xue, Lincoln↗

Development of dual phase soft magnetic laminate for advanced electric machines

The development of dual phase soft magnetic laminate for advanced electric machines at General Electric is summarized, including alloy development, properties, manufacturing trials, prototype experience, potential applications, and future material development. Dual phase soft magnetic laminates provide the flexibility to tune magnetic/non-magnetic state at the desired locations, using solution nitriding treatment with protective mask/coating. This enabling material technology benefits electric machine performance by reducing flux leakage in the laminates. It also decouples the electromagnetic and mechanical design, which opens up the design space for advanced electric machines. Magnetic and mechanical properties of dual phase laminates were compared to commercial FeSi and FeCo soft magnetic laminates. Driven by the market needs for high power density and high-speed electric machines, a synchronous reluctance (SynRel) machine without PMs using dual phase soft magnetic laminates was assessed in machine electromagnetic performance as well as rotor stress analysis by finite element simulation. Successful prototype build and testing further retired manufacturing risks and advanced the technology maturity level. Future material development with higher mechanical strength, higher saturation magnetization, higher permeability, and lower loss are highly desired for ground and air transportation electrification.

Rallabandi, Vandana↗

Design and Engineering of LUPIN: A Test-Bed Radiofrequency Ion Source for Enhanced Neutral Beam Injection on DIII-D

The Large, Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma (ICP) chamber for demonstrating plasma performance of an RF ICP positive ion source upgrade for the DIII-D Neutral Beam Injection (NBI) system. LUPIN will be used to investigate ion source physics, including neutral gas dynamics, plasma density uniformity, interactions with Faraday shields, and power coupling to novel RF antenna designs. LUPIN has an RF generator capable of delivering 20 kW of power at 2 MHz, which is coupled into a cylindrical quartz vessel measuring 20 cm in length and 10 cm in radius. This configuration matches the power density requirements for a full-scale ion source. Target hydrogen and deuterium plasma densities exceeding 10 18 m -3 would relate to extracted ion current densities of 2100 A/m 2 for 10 s. Vacuum conductance and gas flow calculations predict a maximum achievable neutral gas flow rate of 1675 Pa·L/s at 5 Pa of He, which mimics the gas flow of the DIII-D NBI system. Designs have been developed for an internal Faraday shield to mitigate heat flux and ion sputtering on the dielectric vessel. Thermomechanical finite element simulations demonstrated the Faraday shield design to be capable of withstanding anticipated heat loads from worst-case operation scenarios. Results of upcoming experimental investigations on LUPIN will guide the design of a full-scale prototype for DIII-D integration.

Faraday shield↗

Coupling Finite Element and Finite Volume Simulation Within MOOSE

This work focuses on adding capability to couple finite element (FE) variables into finite volume (FV) physics within INL’s Multiphysics Object-Oriented Simulation Environment (MOOSE). This coupling can allow for improvement of multiphysics simulations where one set of physics is best suited for an FE discretization, while another set is best suited for an FV discretization. Electrohydrodynamics, which govern liquid metal reactor concepts and plasma dynamics, is a fitting example where the electromagnetic field equations are solved with FE and the fluid flow is solved with FV. The new FE to FV coupling method can be summarized as taking the element or face average of the FE variable value or gradient and applying that quantity directly in the FV equation objects.

97 MATHEMATICS AND COMPUTING↗

COUPLING SMOOTHED PARTICLE HYDRODYNAMICS WITH FINITE ELEMENT METHOD TO SIMULATE RESIDUAL STRESSES FROM FRICTION STIR PROCESSING

Friction stir processing (FSP) is a solid-state material processing technique that locally modifies the microstructure but also induces undesirable residual stresses. A robust numerical model for the FSP can help in mitigating these residual stresses. Heat source models within a finite element method (FEM) framework suffer from inaccuracies. In contrast, smoothed particle hydrodynamics (SPH) model that explicitly captures the material flow near the tool and the associated heat generation are accurate. However, the computational expense of SPH simulations can be prohibitive. In this work, we propose a coupled SPH-FEM framework. SPH is used to model the heat generation accurately near the tool and which is then inserted into to the FEM model as a heat source. To verify this proposed coupling approach, a test case is set up with typical FSP conditions and it is modeled in both SPH and SPH-FEM. The temperatures profiles were compared after the simulations have reached steady-state temperatures. The similarity of the temperature profiles from SPH-FEA and SPH validated the proposed coupling approach. This proposed approach achieves the accuracy of the SPH method while potentially retaining the low computational expense of FEM.

smoothed particle hydrodynamics, Finite Element Me↗

Implementation of Creep-Damage Model for Concrete Fracture in MOOSE

Numerical implementation of an isotropic creep-damage model for concrete in multiphysics object-oriented simulation environment (MOOSE) finite element framework is presented in this paper. The constitutive model considers the combined effect of instantaneous and delayed strains on damage propagation. The implementation considers creep using generalized Maxwell or Kelvin-Voigt models. Using strain splitting assumptions, the total mechanical strains are split into elastic and creep components. Damage is considered to evolve as a function of the elastic and creep strains. This work considers damage as a function of fracture energy using the characteristic length of each finite element. This approach preserves the energy release rate of each element and avoids vanishing energy dissipation as the mesh is refined. Further, a creep-damage parameter is used to quantify the effect of creep strain on damage. The model is tested against published results on notched three-point bending specimens involving non-linear creep and predicts that about a third of the creep strain contributes towards damage for the experiments simulated. Results show that the proposed framework has predictive capabilities, and the model can be extended for more complex systems.

36 MATERIALS SCIENCE↗

Droplet formation simulation using mixed finite elements

Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The importance of sub-meter-scale snow roughness on conductive heat flux of Arctic sea ice

Abstract The conductive heat flux through the snow and ice is a critical component of the mass and energy budgets in the Arctic sea ice system. We use high horizontal resolution (3–15 cm) measurements of snow topography to explore the impacts of sub-meter-scale snow surface roughness on heat flux as simulated by the Finite Element method. Simulating horizontal heat flux in a variable snow cover modestly increases the total simulated heat flux. With horizontal heat flux, as opposed to simple 1D-vertical heat flux modeling, the simulated heat flux is 10% greater than that for uniform snow with the same mean snow thickness for a 31.5 × 21 m region of sea ice (the largest region we studied). Vertical-only (1D) heat flux simulates just a 6% increase for the same region. However, this is highly dependent on observation resolution. Had we measured the snow cover at 1 m horizontal spacing or greater, simulating horizontal heat flux would not have changed the net heat flux from that simulated with vertical-only heat flux. These findings suggest that measuring and modeling snow roughness at sub-meter horizontal scales may be necessary to accurately represent horizontal heat flux on level Arctic sea ice.

Geology↗

Investigation of Cutter–Woodchip Contact Pressure in a New Biomass Comminution System

A Crumbler ® rotary shear system was recently developed to process fuel-grade woodchips into feedstocks having more uniform particle sizes for improved flowability and higher thermochemical output yield than a traditional hammer mill. It represents a significant innovation in the state of the art of biomass comminution equipment. However, the milling unit has experienced significant abrasive wear of the feeding teeth and cutters when processing hard and dusty feedstocks, such as logging residue and corn stover. Here, we present initial simulation results from a finite-element stress analysis of the rotary shear, which is used to investigate the dynamics of interaction between the cutters and woodchips. This is the first step for optimizing the cutter design to improve the system lifetime and processing efficiency, thus reducing downtime and improving overall productivity.

42 ENGINEERING↗

ExaAM: Metal additive manufacturing simulation at the fidelity of the microstructure

Additive manufacturing (AM), or 3D printing, of metals is transforming the fabrication of components, in part by dramatically expanding the design space, allowing optimization of shape and topology. However, although the physical processes involved in AM are similar to those of welding, a field with decades of experimental, modeling, simulation, and characterization experience, qualification of AM parts remains a challenge. The availability of exascale computational systems, particularly when combined with data-driven approaches such as machine learning, enables topology and shape optimization as well as accelerated qualification by providing process-aware, locally accurate microstructure and mechanical property models. We describe the physics components comprising the Exascale Additive Manufacturing simulation environment and report progress using highly resolved melt pool simulations to inform part-scale finite element thermomechanics simulations, drive microstructure evolution, and determine constitutive mechanical property relationships based on those microstructures using polycrystal plasticity. We report on implementation of these components for exascale computing architectures, as well as the multi-stage simulation workflow that provides a unique high-fidelity model of process–structure–property relationships for AM parts. In addition, we discuss verification and validation through collaboration with efforts such as AM-Bench, a set of benchmark test problems under development by a team led by the National Institute of Standards and Technology.

3D printing↗

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗

Uncertainty quantification in elastic constants of SiC f /SiC m tubular composites using global sensitivity analysis

Silicon carbide fiber and silicon carbide matrix (SiC f /SiC m ) tubes produced through the chemical vapor infiltration process have become a candidate cladding material in nuclear applications. The performance of this composite is influenced by many variables such as braiding angle, porosity, material properties, etc., which vary over a range of values due to the inherent fluctuations in the manufacturing process. In this study, the variability in elastic constants of SiC f /SiC m composite has been quantified through multiscale finite element (FE) simulations, variable screening, and high-fidelity surrogate modeling. The key variables dominantly affecting the elastic constants of SiC f /SiC m tubes were identified using global sensitivity analysis. A surrogate to the high-fidelity FE-based model was used in Monte Carlo simulations to generate a hundred thousand samples from which the uncertainty in elastic constants was assessed. It turned out that the coefficient of variation was less than 10%.

Materials Science↗

Tailored Graphene Micropatterns by Wafer-Scale Direct Transfer for Flexible Chemical Sensor Platform

2D materials, such as graphene, exhibit great potential as functional materials for numerous novel applications due to their excellent properties. The grafting of conventional micropatterning techniques on new types of electronic devices is required to fully utilize the unique nature of graphene. However, the conventional lithography and polymer-supported transfer methods often induce the contamination and damage of the graphene surface due to polymer residues and harsh wet-transfer conditions. In this work, a novel strategy to obtain micropatterned graphene on polymer substrates using a direct curing process is demonstrated. Employing this method, entirely flexible, transparent, well-defined self-activated graphene sensor arrays, capable of gas discrimination without external heating, are fabricated on 4 in. wafer-scale substrates. Finite element method simulations show the potential of this patterning technique to maximize the performance of the sensor devices when the active channels of the 2D material are suspended and nanoscaled. This study contributes considerably to the development of flexible functional electronic devices based on 2D materials.

2D materials, Graphene, Microscale patterning, Che↗

High Throughput Computational Framework of Materials Properties for Extreme Environments

This project aims to establish a framework capable of efficiently predicting the properties of structural materials for service in harsh environments over a wide range of temperatures and over long periods of time. The approach is to develop and integrate high throughput first-principles calculations in combination with machine learning (ML) methods, perform high throughput CALPHAD (calculations of phase diagrams) modeling, and carry out finite element method (FEM) simulations. Relevant to high temperature service in fossil power system, nickel-based superalloys such as Inconel 740 and Haynes 282 as well as the associated (Ni-Cr-Co)-Al-C-Fe-Mn-Mo-Nb-Si-Ti system, were investigated. The present framework was built on the concept of phase-based property data, in which properties of individual phases are modeled as a function of internal and external independent variables. This project established an open-source infrastructure with the following capabilities: (1) High throughput implementation of first-principles calculations at finite temperatures and variable compositions using both accurate phonon calculations and the efficient Debye model for thermodynamic properties, elastic constants, diffusion coefficients, vacancy formation, stacking and twin faults, and dislocation mobility; i.e., using the developed code DFTTK; (2) Machine learning capabilities to predict the above properties so that the number of first-principles calculations can be significantly reduced; e.g., using the developed code SIPFENN; (3) High throughput CALPHAD modeling of the above properties as a function of temperature and composition using our unique capability based on ESPEI and PyCalphad; (4) New capabilities to predict the stress-strain behavior of individual phases; and (5) New models for tensile strength prediction in common FEM software with the crystal plasticity finite element simulations (CPFEM).

, Ni-based superalloys↗