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

Probing the properties and mechanisms of failure waves in soda-lime glass

Soda-lime glass (SLG) and other silica glasses exhibit the failure wave phenomenon under shock compression. The mechanism responsible for this peculiar behavior of glasses is still unresolved. In this study, a series of plate impact experiments was performed at three different impact stresses of 6.4, 8.3, and 10.8 GPa to better understand the mechanisms underlying the failure wave phenomenon. Specifically, spall experiments were conducted to probe the speed and existence of failure waves at different stresses in SLG. A layered glass target was used to probe the possibility of a “renucleation” of the wave at the SLG–SLG interface. When it existed, the failure wave was inferred to propagate at a speed of 1.3 km/s. However, it was observed that the failure wave phenomenon ceases to exist for impact stresses higher than 10 GPa. In experiments with a 6.4 GPa impact stress, the peak free surface velocity was significantly less than what is predicted by stress-Hugoniot calculations. This velocity deficit and other important features of the measured free surface velocity profiles were simulated using finite element analysis by incorporating an abrupt densification of SLG at a critical stress in the equation of state. Furthermore, this densification feature is similar to what would be expected of a phase transition. Although unable to unambiguously reveal the mechanism causing the failure wave phenomenon, the results of the present work clearly indicate that the failure wave causes a secondary compression and densification in SLG.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Low spontaneous Brillouin scattering in anti-resonant hollow-core fibers in GHz frequency range

Brillouin light scattering (BLS) is a powerful experimental tool that can be used to gain insights into the fundamental and applied properties of matter, like dispersions of quasiparticles in a solid, as well as their spatiotemporal dynamics. Many applications of light scattering favor the use of optical fibers in place of free-space optics. In this study, we compare the performance of anti-resonant hollow core fibers to that of conventional solid core fused silica fibers for BLS experiments in the GHz frequency range. Conventional fibers are barely suitable for low-noise measurements because of the spontaneous scattering of photons on various phononic modes present in the core and cladding. In the case of the hollow-core fiber, we identify a range of discrete phononic modes and associate them with the various acoustic modes of the structure surrounding the hollow core using finite-element numerical simulations. The measured relative intensity of the spontaneous BLS signal from these modes is orders of magnitude smaller than that of a solid-core fiber, making anti-resonant hollow-core fibers one of the best solutions for the single-mode light guidance for BLS and potentially other low-noise photonic experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗
DOE OSTI2018

Strain-Induced Spin-Resonance Shifts in Silicon Devices

In spin-based quantum-information-processing devices, the presence of control and detection circuitry can change the local environment of a spin by introducing strain and electric fields, altering its resonant frequencies. These resonance shifts can be large compared to intrinsic spin linewidths, and it is therefore important to study, understand, and model such effects in order to better predict device performance. We investigate a sample of bismuth donor spins implanted in a silicon chip, on top of which a superconducting aluminum microresonator is fabricated. The on-chip resonator provides two functions: it produces local strain in the silicon due to the larger thermal contraction of the aluminum, and it enables sensitive electron spin-resonance spectroscopy of donors close to the surface that experience this strain. Through finite-element strain simulations, we are able to reconstruct key features of our experiments, including the electron spin-resonance spectra. Our results are consistent with a recently observed mechanism for producing shifts of the hyperfine interaction for donors in silicon, which is linear with the hydrostatic component of an applied strain.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Periodic chiral magnetic domains in single-crystal nickel nanowires

We report here on experimental and computational investigations of the domain structure of single-crystal Ni nanowires (NWs). The ${\sim}200\times{}200\times{}8000\mathrm{n}{\mathrm{m}}^{3}$ Ni NWs were grown by a thermal chemical vapor deposition technique that results in single-crystal structures. Magnetoresistance measurements of individual NWs suggest the average magnetization points largely off the NW long axis at zero field. X-ray photoemission electron microscopy images obtained at room temperature show a well-defined periodic magnetization pattern along the surface of the nanowires with a period of ${{\lambda}}_{\mathrm{avg}}=239\pm{}37\mathrm{nm}$. Finite element micromagnetic simulations reveal that when the material parameters of the modeled system match those of nickel crystal at $T=10\mathrm{K}$, an oscillatory magnetization configuration with a period closely matching experimental observation (${\lambda}=222\mathrm{nm}$) is obtainable at remanence. This magnetization configuration involves a periodic array of alternating chirality vortex domains distributed along the length of the NW. Vortex formation is attributable to the relatively high cubic anisotropy of the single crystal Ni NW system at $T=10\mathrm{K}$ and its reduced structural dimensions. The periodic alternating chirality vortex state is a topologically protected metastable state, analogous to an array of 360° domain walls in a thin strip. Simulations show that other remanent states are also possible, depending on the field history. At room temperature ($T=273\mathrm{K}$), simulations show vortices are no longer stable due to the expected reduced cubic anisotropy of the system, suggesting a disparity between the fabricated and modeled nanowires. Negative uniaxial anisotropy and magnetoelastic effects in the presence of compressive biaxial strain are shown to promote and restore formation of vortices at room temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Hydroforming of elliptical cavities

Activities of the past several years in developing the technique of forming seamless (weldless) cavity cells by hydroforming are summarized. An overview of the technique developed at DESY for the fabrication of single cells and multicells of the TESLA cavity shape is given and the major rf results are presented. The forming is performed by expanding a seamless tube with internal water pressure while simultaneously swaging it axially. Prior to the expansion the tube is necked at the iris area and at the ends. Tube radii and axial displacements are computer controlled during the forming process in accordance with results of finite element method simulations for necking and expansion using the experimentally obtained strain-stress relationship of tube material. In cooperation with industry different methods of niobium seamless tube production have been explored. The most appropriate and successful method is a combination of spinning or deep drawing with flow forming. Several single-cell niobium cavities of the 1.3 GHz TESLA shape were produced by hydroforming. They reached accelerating gradients E acc up to 35 MV/m after buffered chemical polishing (BCP) and up to 42 MV/m after electropolishing (EP). More recent work concentrated on fabrication and testing of multicell and nine-cell cavities. Several seamless two- and three-cell units were explored. Accelerating gradients E acc of 30–35 MV/m were measured after BCP and E acc up to 40 MV/m were reached after EP. Nine-cell niobium cavities combining three three-cell units were completed at the company E. Zanon. These cavities reached accelerating gradients of E acc = 30–35 MV/m. One cavity is successfully integrated in an XFEL cryomodule and is used in the operation of the FLASH linear accelerator at DESY. Additionally the fabrication of bimetallic single-cell and multicell NbCu cavities by hydroforming was successfully developed. Several NbCu clad single-cell and double-cell cavities of the TESLA shape have been fabricated. The clad seamless tubes were produced using hot bonding or explosive bonding and subsequent flow forming. The thicknesses of Nb and Cu layers in the tube wall are about 1 and 3 mm respectively. The rf performance of the best NbCu clad cavities is similar to that of bulk Nb cavities. The highest accelerating gradient achieved was 40 MV/m. The advantages and disadvantages of hydroformed cavities are discussed in this paper.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shield Design for 50 kW Three-Phase Wireless Charging System

In this paper, magnetic field emissions (MFE) from a 50 kW bipolar coil-based three-phase wireless charging system (WCS) is investigated. Then, a shield design is proposed to suppress the MFE below the safety limit identified by the International Commission of Non-Ionized Radiation Protection (ICNIRP) guidelines. Bipolar coil-based three-phase (3φ) WCSs provide extremely high power-density; therefore, the 3φ WCSs are highly promising for high-power fast-wireless-charging of the electric vehicles. For the high-power EV application, it is a critical challenge to design a shield to suppress the MFE below the safety limit. The traditional aluminum shield is found ineffective for the bipolar 3φ-WCS. The proposed shield is designed with high-permeability magnetic material, i.e. ferrite, nanocrystalline material, etc. The proposed shield is simulated in finite element analysis (FEA) and the FEA results show that the proposed magnetic shield effectively suppresses the magnetic field emission below the ICNIRP limit. The generalized design of the proposed shield can be used for even higher power bipolar coil based WCSs.

Mohammad, Mostak↗
DOE OSTI2023

Additively Manufactured Copper Windings with Hilbert Structure

In this paper, finite element method simulations and frequency analysis of additively manufactured windings with a Hilbert–cross section pattern were performed over a wide frequency range. The current density distributions for the proposed winding conductor and for a square-wire winding conductor at selected frequencies were compared. An approximately 25% winding loss reduction at higher frequencies and an approximately 5% increase at low frequencies were observed for the analyzed Hilbert structure.

Wojda, Rafal↗
DOE OSTI2023

Design and Optimization of High-frequency Transformer for Isolated Single-Stage Three-Phase AC/DC Converter using Bidirectional Switches

This paper presents the design and optimization of the high-frequency transformer for an isolated single-stage threephase AC/DC converter enabled by the 4H-SiC BiDirectional Field Effect Transistor (BiDFET). Non-dominated Sorting Genetic Algorithm II (NSGA-II) is used to perform a multiobjective optimization to minimize loss, weight, and volume. A 20 kW, 50 kHz hardware prototype based on the selected optimal point is characterized by a power amplifier. The finite element analysis simulation results are also presented to validate the core magnetic flux density and winding electric field strength. The experimental prototype was tested at 800 Vdc/480 Vrms at 10 kW rated power.

BiDFET, High-frequency transformer, Multiobjective↗
DOE OSTI2020

Comparison of Magnetic Field Emission from Unipolar and Bipolar Coil-Based Wireless Charging Systems

In this study, the leakage magnetic field (LMF) from the unipolar (circular, rectangular, etc.) and bipolar (double-D) coil-based wireless charging systems (WCSs) are investigated and compared. The unipolar and bipolar pads, which are the two most widely adopted charging-pads, have significantly different magnetic field emission characteristics. In this paper, the spatial-distribution and vector-pattern of the leakage magnetic field are studied and compared for similar 11 kW unipolar and bipolar coil-based WCSs. The field emissions are simulated using finite element analysis (FEA) and tested experimentally with an 11 kW bipolar coil-based WCS prototype. The results show that the leakage magnetic field of the unipolar and bipolar pads have widely different vector-orientation and spatial distributions of the leakage magnetic field. The results show that the unipolar coil based WCS has a dominant vertical and the bipolar pads coil-based WCS has dominant horizontally aligned leakage magnetic field. Considering the shielding, the aluminum plate is found highly effective for the unipolar pad and almost ineffective for bipolar pads.

Mohammad, Mostak↗
DOE OSTI2020

Computationally Efficient AC Resistance Model for Stator Winding With Rectangular Conductors

The resistance of the electric machine’s windings is an important parameter when analyzing the machine’s performance, as conduction losses in the winding are often the largest loss component. However, this resistance changes with frequency due to the skin and proximity effects, making it challenging to accurately estimate these losses. In this article, a computationally efficient, general method of calculating the ac resistance of an electric machine winding with rectangular conductors is presented. The method uses a 1-D magnetoquasi-static model of the stator slot region to develop boundary transfer relations between individual turns. These relations are then used to create an impedance matrix from which the ac resistance can be calculated. This technique is more flexible than previous solutions since it can accommodate different winding configurations such as fractional-pitch windings. It is shown that a fractional-pitch winding has a lower ac resistance than a full-pitch winding. Finally, finite-element analysis simulations and experimental results are used to validate the proposed method.

42 ENGINEERING↗
DOE OSTI2022

Laminated Permanent Magnets Enable Compact Magnetic Components in Current-Source Converters

Magnetic components have become one of the primary barriers to high power density converters, especially to current-source converters (CSCs). In CSCs, the inductors/transformers have a predominantly dc flux to store energy, which can be offset by standard permanent magnets (PMs). However, eddy currents in standard PMs induce significant losses and thermal stress for medium-frequency applications. This article proposes using laminated PMs to offset the dc flux while reducing eddy currents, leading to significant reductions in the size, cost, and losses of inductors/transformers. Furthermore, the laminated PMs’ optimal location, orientation, and distribution are investigated to generalize this approach for maximum benefits. Here, three-dimensional finite-element analysis simulation and hardware experiments are presented to validate the effectiveness of the proposed approach in a medium-frequency transformer (MFT) for CSCs. Compared with standard PMs, the proposed use of laminated PMs reduces aggregate core-plus-PM losses by 85% in experiments, and thus, relaxes the MFT thermal design. Finally, the proposed approach is experimentally validated in a 40 kVA flyback-type MFT for a soft-switching solid-state transformer. Compared to the traditional design without any PMs, the proposed design increases the saturation current by 46% while inducing only 5% more losses, leading to significant savings in magnetics cost and size.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2018

Nanoplastic removal function and the mechanical nature of colloidal silica slurry polishing

The nanomechanical deformations on a broad range of optical material surfaces (single crystals of Al 2 O 3 [sapphire], SiC, Y 3 Al 5 O 12 [YAG], CaF 2 , and LiB 3 O 5 [LBO]; a SiO 2 –Al 2 O 3 –P 2 O 5 –Li 2 O glass-ceramics [Zerodur]; and glasses of SiO 2 :TiO 2 [ULE], SiO 2 [fused silica], and P 2 O 5 –Al 2 O 3 –K 2 O–BaO [Phosphate]) near the elastic-plastic load boundary have been measured by nanoindentation and nanoscratching to mimic the nanoplastic removal caused by a single slurry particle during polishing. Nanoindenation in air was performed to determine the workpiece hardness at various loads using a commercial nanoindenter with a Berkovich tip. Similarly, an atomic force microscope (AFM) with a stiff diamond coated tip (150 nm radius) was used to produce nanoplastic scratches in air and aqueous environments over a range of applied loads (~20-170 μN). The resulting nanoplastic deformation of the nanoscratches were used to calculate the removal function (i.e., depth per pass) which ranged from 0.18 to 3.6 nm per pass for these materials. A linear correlation between the nanoplastic removal function and the polishing rate (using a fixed polishing process with colloidal silica slurry on a polyurethane pad) of these materials was observed implying that: (a) the polishing mechanism using colloidal silica slurry can be dominated by mechanical rather than chemical interactions; and (b) the nanoplastic removal function, as opposed to interface particle interactions, is the controlling factor for the polishing material removal rate. Furthermore, this correlation is consistent with the Ensemble Hertzian Multi-Gap (EHMG) microscopic material removal rate model described previously. The nanoplastic removal depth was also found to correlate to the measured nanoindentation hardness (H1) of the optical material, scaling as H 1 -3.5 . Two-dimensional (2D) finite element analysis simulations of nanoindentation showed a similar nonlinear dependence of plastic deformation with the workpiece material hardness. The findings of this study are used to determine an effective Preston coefficient for the material removal rate expression and enhance the predictive nature of the nanoplastic polishing rate for various materials utilizing their material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mechanism-Based Energy Regularization in Computational Modeling of Quasibrittle Fracture

Abstract Quasibrittle materials are featured by a strain-softening constitutive behavior under many loading scenarios, which could eventually lead to localization instability. It has long been known that strain localization would result in spurious mesh sensitivity in finite element (FE) simulations. Previous studies have shown that, for the case of fully localized damage, the mesh sensitivity can be mitigated through energy regularization of the material constitutive law. However, depending on the loading configuration and structural geometry, quasibrittle structures could exhibit a complex damage process, which involves both localized and diffused damage patterns at different stages of loading. This study presents a generalized energy regularization method that considers the spatial and temporal evolution of damage pattern. The method introduces a localization parameter, which describes the local damage pattern. The localization parameter governs the energy regularization of the constitutive model, which captures the transition from diffused to localized damage during the failure process. The method is cast into an isotropic damage model, and is further extended to rate-dependent behavior. The energy regularization scheme is directly incorporated into the kinetics of damage growth. The model is applied to simulate static and dynamic failures of ceramic specimens. It is shown that the present model is able to effectively mitigate the spurious mesh sensitivity in FE simulations of both types of failure. The present analysis demonstrates the essential role of mechanism-based energy regularization of constitutive relation in FE simulations of quasibrittle fracture.

Mechanics↗
DOE OSTI2021

Design of a Single Layer Metamaterial for Pressure Offloading of Transtibial Amputees

While using a prosthesis, transtibial amputees may experience pain and discomfort brought on by large pressure gradients at the interface between the residual limb and prosthetic socket. Current prosthetic interface solutions attempt to alleviate these pressure gradients by using soft homogenous liners to reduce and distribute pressures. This research investigates an additively manufactured metamaterial inlay with a tailored mechanical response to reduce peak pressure gradients around the limb. The inlay uses a hyperelastic behaving metamaterial (US10244818) comprised of triangular pattern unit cells which can be 3D printed with walls of various thicknesses controlled by draft angles. The hyperelastic material properties are modeled using a Yeoh 3rd order model. The 3rd order coefficients can be adjusted and optimized, which corresponds to a change in the unit cell wall thickness to create an inlay that can meet the unique offloading needs of an amputee. Finite element analysis simulations evaluated the pressure gradient reduction from: 1) a common homogenous silicone liner, 2) a prosthetist's inlay prescription that utilizes three variations of the metamaterial, and 3) a metamaterial solution with optimized Yeoh 3rd order coefficients. When compared to a traditional homogenous silicone liner for two unique limb loading scenarios, the prosthetist prescribed inlay and the optimized material inlay can achieve equal or greater pressure gradient reduction capabilities. These results show the potential feasibility of implementing this metamaterial as a method of personalized medicine for transtibial amputees by creating a customizable interface solution to meet the unique performance needs of an individual patient.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Impact of static magnetic fields on the radial line slot antenna plasma source

The radial line slot antenna plasma source is used in semiconductor device fabrication. As is the case for all plasma sources, ever more strict uniformity control requirements are driven by the precision demands of new device technologies. Large volume diffusion plasmas, of which the radial line slot antenna source is one type, must overcome transport effects or diffusion modes that tend to 'center peak' the plasma density near the wafer being processed. One way to resolve problematic transport effects is the insertion of magnetic fields into the plasma region. In this paper, the impact of the magnetic field on plasma properties is parameterized as a function of slot configuration. The magnetic field orientation and the magnitude of magnetic field are varied in a computational study in which the source is modeled as a two-dimensional axisymmetric quasineutral plasma. This work employs a finite element model simulation. The magnitude of magnetic fields considered is 50 Gauss maximum with a microwave power of 3000 W at a pressure of 20 mTorr. 20 mTorr is chosen as this is a condition where diffusion effects are challenging to counteract. The study showed that there are specific conditions for slot configuration and magnetic field that improve the plasma controllability and some that do not. Plasma property modulation is most effective when the plasma source region is placed at large radius with axial magnetic fields. There are synergistic effects between the slot location and magnetic field that are important and placing large magnetic fields at the chamber edge alone does not result in improved uniformity. Electron cyclotron resonance (ECR) heating and the impact of pulsing the magnetic fields are presented. ECR heating is not important for the conditions relevant to this paper and pulsing is shown to have benefit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

MO-E-17A-12: Direct Realization of the CT Dose to Phantom: Energy to Heat Conversion in Polyethylene Using Calorimetry

Purpose: To develop a primary reference standard for absorbed dose to phantom for medical CT dosimetry. CT dosimetry relies on the implementation of the CTDI standard based on air kerma. We are taking a step toward an absorbed dose to water standard by first investigating the dose in a solid phantom. By directly measuring the heat converted from the incident radiation, the absorbed dose in the phantom at a point can be assessed with primary methods without relying on indirect conversions. Methods: The calorimeter contains two small thermistors embedded in a removable PE “core” inserted into the cylindrical HDPE phantom. A core made with polystyrene (PS) was also tested because of its purportedly negligible heat defect. Measurements were made using the two cores and with a calibrated ionization chamber in a CT beam. The air chamber values were converted to the dose to medium using appropriate stopping-power ratios from the literature, and then compared to the thermal data. Results: The PS core data yielded a dose of 1.3 times (4-run average, 3% std. dev.) higher than the converted chamber value, whereas the PE core data were inexplicably higher. The possible systematic errors include 1) excess heat from the thermistors, 2) in PE the exothermic chemical reactions, 3) uncertainties of the specific heat capacities of the materials, 4) thermal drift, and 5) theoretical conversion of chamber values. Monte Carlo simulations and finite element heat transfer calculations were performed to address some of these issues. The general validity was assessed in a 6 MV photon beam with an entirely different calibration scheme. Conclusion: This study demonstrates the feasibility but also revealed the difficulty in developing a new primary reference standard for absorbed dose to material for CT. Additional experimental and theoretical work is planned to achieve our goal.

07 ISOTOPE AND RADIATION SOURCES↗