Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element Method”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Modelling DC responses of 3D complex fracture networks

Here, the determination of the geometrical properties of fractures plays a critical role in many engineering problems to assess the current hydrological and mechanical states of geological media and to predict their future states. However, numerical modeling of geoelectrical responses in realistic fractured media has been challenging due to the explosive computational cost imposed by the explicit discretizations of fractures at multiple length scales, which often brings about a tradeoff between computational efficiency and geologic realism. Here, we use the hierarchical finite element method to model electrostatic response of realistically complex 3D conductive fracture networks with minimal computational cost.

58 GEOSCIENCES↗
DOE OSTI2019

Electromagnetic wave propagation based upon spectral-element methodology in dispersive and attenuating media

SUMMARY We build on mathematical equivalences between Maxwell’s wave equations for an electromagnetic medium and elastic seismic wave equations. This allows us to readily model Maxwell’s wave propagation in the spectral-element codes SPECFEM2D and SPECFEM3D, written for acoustic, viscoelastic and poroelastic seismic wave propagation, providing the ability to handle complex geometries, inherent to finite-element methods and retaining the strength of exponential convergence and accuracy due to the use of high-degree polynomials to interpolate field functions on the elements, characteristic to spectral-element methods (SEMs). Attenuation and dispersion processes related to the frequency dependence of dielectric permittivity and conductivity are also included using a Zener model, similar to shear attenuation in viscoelastic media or viscous diffusion in poroelastic media, and a Kelvin–Voigt model, respectively. Ability to account for anisotropic media is also discussed. Here, we limit ourselves to certain dielectric permittivity tensor geometries, in order to conserve a diagonal mass matrix after discretization of the system of equations. Doing so, simulation of Maxwell’s wave equations in the radar frequency range based on SEM can be solved using explicit time integration schemes well suited for parallel computation. We validate our formulation with analytical solutions. In 2-D, our implementation allows for the modelling of both a transverse magnetic (TM) mode, suitable for surface based reflection ground penetration radar type of applications, and a transverse electric (TE) mode more suitable for crosshole and vertical radar profiling setups. Two 2-D examples are designed to demonstrated the use of the TM and TE modes. A 3-D example is also presented, which allows for the full TEM solution, different antenna orientations, and out-of-plane variations in material properties.

58 GEOSCIENCES↗
DOE OSTI2026

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗
DOE OSTI2025

Parameter extraction for a SPICE model of an hTron superconducting thermal switch

Efficiently simulating large circuits is crucial to the development of superconducting nanowire-based electronics. However, current simulation tools for this technology are not adapted to the scaling of circuit size and complexity. We focus on the multilayered heater-nanocryotron (hTron), a promising superconducting nanowire-based switch used in applications such as superconducting nanowire single-photon detector readout. Previously, the hTron was modeled using traditional finite-element methods, which fall short in simulating systems at a larger scale. An empirical-based method would be better adapted to this task, enhancing both simulation speed and agreement with experimental data. In this work, we perform switching current and activation delay measurements on 17 hTron devices. We then develop a method for extracting physical fitting parameters used to characterize the devices. We build a SPICE behavioral model that reproduces the static and transient device behavior using these parameters, and validate it by comparing its performance to a model developed in prior work, showing an improvement in simulation time by several orders of magnitude. Furthermore, our model provides circuit designers with a tool to help understand the hTron’s behavior during all design stages, thus enabling broader use of the hTron across various new areas of application.

Caloritronics↗
DOE OSTI2024

Conceptual design of a high reactive-power ferroelectric fast reactive tuner

We present a novel design of a ferroelectric fast reactive tuner (FE-FRT) capable of modulating mega-VAR reactive power on a submicrosecond timescale. The high reactive power capability of our design extends the range of applications of reactive tuners to numerous applications. We present a detailed analytical model of the performance of a megawatt-class reactive power device and benchmark it against finite-element method eigenmode and frequency domain electromagnetic simulations. We introduce new features, including an annulus design for the ferroelectric capacitors and capacitive window coupling to the cavity. We consider thermal design issues and nonlinear effects in the ferroelectric. The model covers several configurations, allowing control of the frequency of superconducting and normal-conducting cavities in a variety of applications and frequencies. We calculate that the FE-FRT designed should be capable of handling around 0.45 MVAR of reactive power with around 3 kW of resistive losses, providing a frequency tuning range of 8 kHz in an example of 400 MHz cavity geometry. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Long lifetime of thermally excited magnons in bulk yttrium iron garnet

Spin currents are generated within the bulk of magnetic materials due to heat flow, an effect called intrinsic spin Seebeck. This bulk bosonic spin current consists of a diffusing thermal magnon cloud, parametrized by the magnon chemical potential (μ m ), with a diffusion length of several microns in yttrium iron garnet (YIG). Transient optothermal measurements of the spin-Seebeck effect (SSE) as a function of temperature reveal the time evolution of μ m due to intrinsic SSE in YIG. The interface SSE develops at times <2 ns while the intrinsic SSE signal continues to evolve at times >500 μs, dominating the temperature dependence of SSE in bulk YIG. Time-dependent SSE data are fit to a multitemperature model of coupled spin/heat transport using the finite-element method (FEM), where the magnon spin lifetime (τ) and magnon-phonon thermalization time (τ mp ) are used as fit parameters. From 300 to 4 K, τ mp ) varies from 1 to 10 ns, whereas tau varies from 2 to 60 μ s with the spin lifetime peaking at 90 K. At low temperature, a reduction in τ is observed consistent with impurity relaxation reported in ferromagnetic resonance measurements. Finally, these results demonstrate that the thermal magnon cloud in YIG contains extremely low-frequency magnons (similar to 10 GHz), providing spectral insight to the microscopic scattering processes involved in magnon spin/heat diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Eigenmodes of a quartz tuning fork and their application to photoinduced force microscopy

Here, we examine the mechanical eigenmodes of a quartz tuning fork (QTF) for the purpose of facilitating its use as a probe for multifrequency atomic force microscopy (AFM). We perform simulations based on the three-dimensional finite element method and compare the observed motions of the beams with experimentally measured resonance frequencies of two QTF systems. The comparison enabled us to assign the first seven asymmetric eigenmodes of the QTF. We also find that a modified version of single beam theory can be used to guide the assignment of mechanical eigenmodes of QTFs. The usefulness of the QTF for multifrequency AFM measurements is demonstrated through photoinduced force microscopy measurements. By using the QTF in different configurations, we show that the vectorial components of the photoinduced force can be independently assessed and that lateral forces can be probed in true noncontact mode.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Operator product expansion for radial lattice quantization of 3D ϕ 4 theory

At its critical point, the three-dimensional lattice Ising model is described by a conformal field theory (CFT), the 3D Ising CFT. Instead of carrying out simulations on Euclidean lattices, we use the quantum finite elements method to implement radially quantized critical ϕ 4 theory on simplicial lattices approaching R × S 2 . Computing the four-point function of identical scalars, we demonstrate the power of radial quantization by the accurate determination of the scaling dimensions Δ ε and Δ T as well as ratios of the operator product expansion coefficients f σ σ ε and f σ σ T of the first spin-0 and spin-2 primary operators ε and T of the 3D Ising CFT. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2025

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Advancing material modeling in hydrocodes using a concurrent finite-element and molecular dynamics multiscale framework

We present a multiscale simulation framework that couples the finite-element method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99% efficiency. These results highlight the framework's potential for large-scale multiscale modeling across a broad range of materials and conditions.

Computer science↗
DOE OSTI2023

Genetic Algorithm-Based Commutation Angle Control for Torque Ripple Mitigation in Switched Reluctance Motor Drives

This work addresses the application of the Genetic Algorithm (GA) technique to optimize the commutation angles of a 2 kW 8/6 switched reluctance machine (SRM). The primary goal is to reduce the well-known drawback of SRMs: the torque ripple. Firstly, the machine was modeled in Matlab /Simulink ® using lookup tables obtained via finite element method (FEM) simulations. Subsequently, the model was used to perform the GA routine aiming to find the optimal phase commutation angles that minimize the torque ripple factor. Notably, the torque performance of the SRM was significantly affected by the commutation angles during the search for the optimal solution. Afterwards, the GA results for four different operation points were verified experimentally through a developed drive platform with digital signal processor-based (DSP) control and an asymmetric bridge converter. As showed by the experiments, the proposed approach was suitable to reduce the torque ripple by more than 50% for one of the evaluated operating points. Furthermore, it was confirmed that the torque ripple mitigation led to acoustic noise improvement.

42 ENGINEERING↗
DOE OSTI2025

Flexible PCB Windings Size Optimization for Winding AC Resistance Minimization Under the Sinusoidal Voltage Regime

In this paper, optimization of the winding AC resistance of the flexible printed circuit board (FPCB) wire windings is performed. A one-dimensional model of the FPCB winding is introduced and an equation for the FPCB winding AC resistance and its low and medium frequency approximation are derived. The approximate FPCB winding AC resistance equation is used to derive optimum thickness of the FPCB winding conductor thickness at which minimum of winding AC resistance (global optimum) is achieved. The Finite Element Method analysis and experimental verification of derived equations (winding resistance, impedance, power loss, and temperature measurements) is performed in order validate derived model and winding resistance equation.

42 ENGINEERING↗
DOE OSTI2019

Design of a loop-gap resonator with bimodal uniform fields using finite element analysis

The loop-gap resonator (LGR) was originally developed to provide a uniform microwave magnetic field on a sample for electron spin resonance (ESR) experiments. The LGR is composed of one or more loops and gaps acting as inductances and capacitances respectively. Typical LGR designs produce a uniform field on a sample at a single resonant frequency, but for certain experiments it is necessary to study the response of a material to uniform fields at multiple frequencies applied simultaneously. In this study we develop an empirical design procedure using finite element method calculations to design an asymmetric loop-gap resonator with uniform fields at two frequencies in the same sample volume and analyze the field uniformity, frequency tunability and filling factors, providing comparison to a manufactured device.

Computational electromagnetics↗
DOE OSTI2014

Design and development of permanent magnet based focusing lens for J-Band Klystron

Applying permanent magnet technology to beam focusing in klystrons can reduce their power consumption and increase their reliability of operation. Electromagnetic design of the beam focusing elements, for high frequency travelling wave tubes, is very critical. The magnitude and profile of the magnetic field need to match the optics requirement from beam dynamics studies. The rise of the field from cathode gun region to the uniform field region (RF section) is important as the desired transition from zero to peak axial field must occur over a short axial distance. Confined flow regime is an optimum choice to minimize beam scalloping but demands an axial magnetic field greater than 2 - 3 times the Brillouin flow field. This necessitates optimization in the magnet design achieve high magnetic field within given spatial constraints. Electromagnetic design and simulations were done using 3D Finite element method (FEM) analysis software. A permanent magnet based focusing lens for a miniature J-Band klystron has been designed and developed at Control Instrumentation Division, BARC. This paper presents the design, simulation studies, beam transmission and RF tests results for J Band klystron with permanent magnet focusing lens. (author)

43 PARTICLE ACCELERATORS↗
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

Mechanical Analysis of Carbon Fiber Retaining Sleeve for A High-speed Outer Rotor SPM Electric Motor Design

Outer rotor motors can be designed with a larger airgap diameter than inner rotor motors for the same overall diameter and therefore provide higher torque which makes them ideal for high torque space constrained applications. This advantage fits well with the increasing demand of high-speed electric machines in many applications due to their inherent high-power density. However, increasing the motor speed results in significant centrifugal loads for the outer rotor motor design. Layers of carbon fiber sleeve winding on the rotor outer surface is necessary to guarantee its structural integrity without introducing too much weight on the motor itself. In this paper, we studied the mechanical stress of a carbon fiber retaining sleeve in a surface permanent magnet outer rotor motor that spins up to 20,000 RPM. Analytical results from finite element method (FEM) confirmed that the maximum stress experienced by the carbon fiber sleeve was under the material’s yield strength. The proposed sleeve will allow this outer rotor motor to operate at high speeds.

Lin, Lianshan↗