Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite differences”

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

Robust long-range magnetic correlation across antiphase domain boundaries in Sr 2 CrReO 6

Here, we report a resonant elastic x-ray scattering study of a thin-film sample of Sr 2 CrReO 6 , which has one of the highest ferrimagnetic transition temperatures among ordered double perovskites.We found resonantly enhanced magnetic Bragg peaks at Q = (odd, odd, odd) at both the rhenium L 2 and L 3 edges, which coincide with the structural Bragg peaks of Sr 2 CrReO 6 . By analyzing the widths of these Bragg peaks, we extracted very different structural and magnetic correlation lengths. The former is about 15 nm, while the latter is constrained by the instrumental resolution to be at least 90 nm. We argue that a finite structural correlation length is consistent with the existence of antiphase nanodomains in our sample. On the other hand, a much larger magnetic correlation length indicates that the magnetic correlation extends far beyond the boundaries of individual domains and is consistent with strong antiferromagnetic coupling between different antiphase domains. Last, from the azimuthal dependence of the magnetic intensity, we show that the magnetic moments lie perpendicular to the c axis, which explains the earlier bulk magnetization data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A transfer learning approach for acoustic emission zonal localization on steel plate-like structure using numerical simulation and unsupervised domain adaptation

The detection and localization of damage in metallic structures using acoustic emission (AE) monitoring and artificial intelligence technology such as deep learning has been widely studied. However, a current challenge of this approach is the difficulty of obtaining sufficient labeled historical AE signals for the training process of deep learning models. This problem can be approached through the implementation of transfer learning. The innovation of this paper lies in the development of a transfer learning approach for AE source localization on a stainless-steel structure when no historical labeled AE signals are available for training. A finite element model is developed to generate numerical AE signals for the training. Unsupervised domain adaptation (UDA) technology is utilized to reduce the distribution difference between the numerical and the realistic AE signals and to derive the localization results of the unlabeled realistic AE signals. Finally, the results suggest that the proposed approach is capable of localizing AE signals with high accuracy in the absence of labeled training data.

42 ENGINEERING↗
DOE OSTI2022

Quantum membrane phases in synthetic lattices of cold molecules or Rydberg atoms

Here, we calculate properties of dipolar interacting ultracold molecules or Rydberg atoms in a semisynthetic three-dimensional configuration—one synthetic dimension plus a two-dimensional real-space optical lattice or periodic microtrap array—using the stochastic Green's function quantum Monte Carlo method. Through a calculation of thermodynamic quantities and appropriate correlation functions, along with their finite-size scalings, we show that there is a second-order transition to a low-temperature phase in which two-dimensional “sheets” form in the synthetic dimension of internal rotational or electronic states of the molecules or Rydberg atoms, respectively. Simulations for different values of the interaction V, which acts between atoms or molecules that are adjacent both in real and synthetic space, allow us to compute a phase diagram. We find a finite-temperature transition at sufficiently large V as well as a quantum phase transition—a critical value V c below which the transition temperature vanishes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Novel Transformer with Variable Leakage and Magnetizing Inductances

This paper introduces the concept of a novel variable inductance transformer (VIT) whose leakage and magnetizing inductances can be varied dynamically and independently to aid the advancement of research in isolated power electronic converters. In this transformer, the leakage inductance can be varied by varying the length of overlap between the primary and secondary windings, while the magnetizing inductance can be varied by varying the air gap between the core legs. Individual controls for the two inductances ensure that one can be varied independent of the other, as needed to satisfy different power conversion objectives. In this paper, the analytical models of the variable leakage and magnetizing inductances are presented. Analytical results are compared with those obtained from Finite Element Method (FEM) and an experimental prototype of a VIT. Analytical solutions of the variable leakage inductance present an error of less than 2.35 % when compared to FEM solutions, while the analytical solutions of the variable magnetizing inductance present a maximum error of 1.82 % when compared to experimental measurements.

image method, leakage inductance, magnetizing indu↗
DOE OSTI2023

Calculation fo the Reflection and Transmission Coefficients for Gradual Tapered Waveguide Structures

It is required to find the reflection and transmission coefficients for a tapered waveguide. It is assumed that the taper is very long and gradual such that no higher order modes are excited from the discontinuity of the waveguide. Three different methods are presented for finding the reflection coefficient and two methods are presented for finding the transmission coefficient. A traditional differential equation method which finds only the reflection coefficient is presented. An additional set of differential equations are presented which finds both the reflection and transmission coefficients. A FEM (finite element method), which solves for both the reflection and transmission coefficients, is presented with its limitations. An example of a tapered coaxial line and waveguide is presented with the different methods compared and contrasted.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

A smooth contact algorithm for the combined finite discrete element method

From its inception, the combined finite discrete element method has used a distributed potential contact force algorithm to resolve interaction between finite elements. The contact interaction algorithm relies on evaluation of the contact force potential field. The problem with existing algorithms is that the potential field introduces artificial numerical non-smoothness in the contact force. This work introduces a smooth potential field based on the finite element topology, and a generalized contact interaction law is constructed on top of the smooth potential field. Further, a number of validation cases for the proposed algorithm, considering different shapes of discrete elements, are presented, and detailed aspects of the proposed contact interaction law are tested with numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Non-intrusive reduced order modeling of natural convection in porous media using convolutional autoencoders: Comparison with linear subspace techniques

Natural convection in porous media is a highly nonlinear multiphysical problem relevant to many engineering applications (e.g., the process of CO 2 sequestration). Here, we extend and present a non-intrusive reduced order model of natural convection in porous media employing deep convolutional autoencoders for the compression and reconstruction and either radial basis function (RBF) interpolation or artificial neural networks (ANNs) for mapping parameters of partial differential equations (PDEs) on the corresponding nonlinear manifolds. To benchmark our approach, we also describe linear compression and reconstruction processes relying on proper orthogonal decomposition (POD) and ANNs. Further, we present comprehensive comparisons among different models through three benchmark problems. The reduced order models, linear and nonlinear approaches, are much faster than the finite element model, obtaining a maximum speed-up of 7 × 10 6 because our framework is not bound by the Courant–Friedrichs–Lewy condition; hence, it could deliver quantities of interest at any given time contrary to the finite element model. Our model’s accuracy still lies within a relative error of 7% in the worst-case scenario. We illustrate that, in specific settings, the nonlinear approach outperforms its linear counterpart and vice versa. We hypothesize that a visual comparison between principal component analysis (PCA) and t-Distributed Stochastic Neighbor Embedding (t-SNE) could indicate which method will perform better prior to employing any specific compression strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Reversible Switch in Charge Storage Enabled by Selective Ion Transport in Solid Electrolyte Interphase

Solid-electrolyte interphases (SEIs) in advanced rechargeable batteries ensure reversible electrode reactions at extreme potentials beyond the thermodynamic stability limits of electrolytes by insulating electrons while allowing working ions to transport. Such selective ion transport occurs naturally in biological cell membranes as a ubiquitous prerequisite of many life processes and a foundation of biodiversity. In addition, cell membranes can selectively open and close the ion channels in response to external stimuli (e.g., electrical, chemical, mechanical, thermal), giving rise to “gating” mechanisms that help manage intracellular reactions. We wondered whether the chemistry and structure of SEIs can mimic cell membranes, such that ion gating can be replicated. That is, can SEIs realize a reversible switching between two electrochemical behaviors, i.e., the ion intercalation chemistry of batteries and the ion adsorption of capacitors? Herein, we report such SEIs that result in thermally activated selective ion transport. The function of open/close gate switches is governed by the chemical and structural dynamics of SEIs under different thermal conditions, with precise behaviors as conducting and insulating interphases that enable battery and capacitive processes within a finite temperature window. Such an ion gating function is synergistically contributed by Arrhenius-activated ion transport and SEI dissolution/regrowth. Following the understanding of this new mechanism, we then develop an electrochemical method to heal the SEI layer in situ. As a result, the knowledge acquired in this work reveals the possibility of hitherto unknown biomimetic properties of SEIs, which will guide us to leverage such complexities to design better SEIs for future battery chemistries.

25 ENERGY STORAGE↗
DOE OSTI2020

Toroidal modeling of runaway electron loss due to 3-D fields in DIII-D and COMPASS

The 3-D field induced relativistic runaway electron (RE) loss has been simulated for DIII-D, COMPASS and ITER plasmas, utilizing the MARS-F code incorporated with the recently developed and updated RE orbit module (REORBIT). Modeling shows effectively 100% loss of a post-disruption, high-current runaway beam in DIII-D, due to the 1 kG level of magnetic field perturbation produced by a fast growing n = 1 resistive kink instability. The RE loss is shown to be independent of the particle energy or the initial location of particles in the configuration space. Applied resonant magnetic perturbation (RMP) fields from in-vessel coils are not effective for RE beam mitigation in DIII-D, but do produce finite (>10%) RE loss in COMPASS, consistent with experimental observations in above two devices. The major reasons for this difference in RE control by RMP between these two devices are (i) the coil proximity to the RE beam and (ii) the effective coil current scaling versus the machine size and the toroidal magnetic field. About 10% RE loss fraction is also predicted for an ITER 15 MA scenario with pre-disruption plasma, highlighting the role of the plasma response. Up to 30% loss is computed, however, by artificially scaling the equilibrium pressure to zero. This is due to the more resistive plasma response and stronger resulting field line stochasticity. Distributions of the lost REs to the limiting surface show poloidally peaked profile near the high-field-side in both DIII-D and COMPASS, covering about 100o poloidal angle. Higher perturbation field level and/or higher particle energy also result in REs being lost to the low-field-side of the limiting surface of these two devices, increasing the effective wetted area. Finally, for the modeled ITER plasmas, the applied RMP field, with optimal poloidal spectrum that maximizes the edge localized mode control, leads to REs being lost to the lower divertor region of the limiting surface within a narrow poloidal band.

3-D field↗
DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Ultrafast field-driven valley polarization of transition metal dichalcogenide quantum dots

Abstract We study theoretically the electron dynamics of transition metal dichalcogenide (TMDC) quantum dots (QDs) in the field of an ultrashort and ultrafast circularly polarized optical pulse. The QDs have the shape of a disk and their electron systems are described within an effective model with infinite mass boundary conditions. Similar to TMDC monolayers, a circularly polarized pulse generates ultrafast valley polarization of such QDs. The dependence of the valley polarization on the size of the dot is sensitive to the dot material and, for different materials, show both monotonic increase with the dot radius and nonmonotonic behavior with a local maximum at a finite dot radius.

Physics↗
DOE OSTI2024

Topological Green’s Function Zeros in an Exactly Solved Model and Beyond

The interplay of topological electronic band structures and strong interparticle interactions provides a promising path towards the constructive design of robust, long-range entangled many-body systems. As a prototype for such systems, we here study an exactly integrable, local model for a fractionalized topological insulator. Using a controlled perturbation theory about this limit, we demonstrate the existence of topological bands of zeros in the exact fermionic Green’s function and show that in this model they do affect the topological invariant of the system, but not the quantized transport response. Close to (but prior to) the Higgs transition signaling the breakdown of fractionalization, the topological bands of zeros acquire a finite “lifetime.” We also discuss the appearance of edge states and edge zeros at real space domain walls separating different phases of the system. This model provides a fertile ground for controlled studies of the phenomenology of Green’s function zeros and the underlying exactly solvable lattice gauge theory illustrates the synergetic cross pollination between solid-state theory, high-energy physics, and quantum information science.

Bollmann, Steffen↗
DOE OSTI2026

Optimal Litz Wire Design for High-Power Wireless Charging System

In this study, an optimized Litz wire and polyphase coil winding design is proposed for high-power wireless charging systems operating at 85 kHz frequency. Litz wire in a wireless charging system contributes to more than half of the losses of the system, making it a critical design component. Here, the Litz wire is optimized through Finite Element Analysis (FEA) considering electrical and mechanical design aspects and the design was validated through experimental tests on different Litz wires with various strand-diameters, twisting ratio, and number of strands. The proposed optimal Litz wire maximizes the system efficiency considering the tradeoff between coil resistance and mutual inductance. The optimal Litz wire shows that Litz wire constructed with the 63µm strands significantly outperform the traditionally used Litz wire with 100µm strand diameter. Additionally, the coil with optimal Litz wire design shows 60% reduction in the copper volume for the same loss, and at the same time increase the coil-to-coil efficiency by more than 2% by increasing the mutual inductance and reducing coil current. The optimized Litz wire was implemented in a 22-kW wireless charging system achieving 94.7% efficiency.

Aydin, Emrullah [ORNL]↗
DOE OSTI2020

Graphite-Embedded High-Performance Insulated Metal Substrate for Wide-Bandgap Power Modules

Emerging wide-bandgap (WBG) semiconductor devices such as silicon carbide (SiC) metal–oxide semiconductor field-effect transistors (MOSFETs) and gallium nitride high-electron-mobility transistors can handle high power in reduced semiconductor areas better than conventional Si-based devices owing to superior material properties. With increased power loss density in a WBG-based converter and reduced die size in power modules, thermal management of power devices must be optimized for high performance. Here, we present a graphite-embedded insulated metal substrate (thermally-annealed-pyrolytic-graphite-embedded insulated metal substrate—IMSwTPG) designed for WBG power modules. Theoretical thermal performance analysis of graphite-embedded metal cores is presented, with design details for IMSwTPG with embedded graphite to replace a direct-bonded copper (DBC) substrate. The proposed IMSwTPG is compared with an aluminum nitride-based DBC substrate using finite-element thermal analysis for steady-state and transient thermal performance. The solutions’ thermal performances are compared under different coolant temperature and thermal loading conditions, and the proposed substrate's electrical performance is validated with static and dynamic characterization. Using graphite-embedded substrates, junction-to-case thermal resistance of SiC MOSFETs can be reduced up to 17%, and device current density can be increased by 10%, regardless of the thermal management strategy used to cool the substrate. Reduced transient thermal impedance of up to 40% of dies owing to increased heat capacity is validated in transient thermal simulations and experiments. The half-bridge power module's electrical performance is evaluated for on -state resistance, switching performance, and switching loss at three junction temperature conditions. The proposed substrate solution has minimal impact on conduction and switching performance of SiC MOSFETs.

42 ENGINEERING↗
DOE OSTI2021

A Hybrid Finite-Volume, Discontinuous Galerkin Discretization for the Radiative Transport Equation

In this report we propose a hybrid spatial discretization for the radiative transport equation that combines a second-order discontinuous Galerkin (DG) method and a second-order finite-volume (FV) method. The strategy relies on a simple operator splitting that has been used previously to combine different angular discretizations. Unlike standard FV methods with upwind fluxes, the hybrid approach is able to accurately simulate problems in scattering dominated regimes. However, it requires less memory and yields a faster computational time than a uniform DG discretization. In addition, the underlying splitting allows naturally for hybridization in both space and angle. Numerical results are given to demonstrate the efficiency of the hybrid approach in the context of discrete ordinate angular discretizations and Cartesian spatial grids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Selective Heating Through Y-Junction Waveguide Designed by Acoustic Shape Optimization

Unlike phononic crystals or systems designed by topology optimization, waveguides designed by shape optimization do not have voids or internal defects, making the fabrication process more suitable for additive manufacturing. By designing a Y-junction waveguide through shape optimization, an ultrasonic wave can be controlled so that it propagates to a predetermined location just by adjusting its frequency. These demultiplexed ultrasonic waves can be used to transport signals or stimulate nearby materials. As an example, the ultrasonic wave is converted to heat at different locations, which can be applied to mechanisms that can take advantage of heating. First, shape optimization is performed on a cylindrical structure to selectively propagate ultrasonic waves of a particular frequency while attenuating others, which is analyzed through a finite element model. The numerical study results are compared with experimental measurements from samples fabricated through additive manufacturing methods. After verifying the concept, the Y-junction waveguide is fabricated to demultiplex the wave and selectively heat different locations. Here the results show that the method of combining shape optimization with additive manufacturing is exceptionally simple and capable of demultiplexing ultrasonic waves, which can replace complex electrical components with single-material waveguides.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Brown-York charges at null boundaries

The Brown-York stress tensor provides a means for defining quasilocal gravitational charges in subregions bounded by a timelike hypersurface. We consider the generalization of this stress tensor to null hypersurfaces. Such a stress tensor can be derived from the on-shell subregion action of general relativity associated with a Dirichlet variational principle, which fixes an induced Carroll structure on the null boundary. The formula for the mixed-index tensor T i j takes a remarkably simple form that is manifestly independent of the choice of auxiliary null vector at the null surface, and we compare this expression to previous proposals for null Brown-York stress tensors. The stress tensor we obtain satisfies a covariant conservation equation with respect to any connection induced from a rigging vector at the hypersurface, as a result of the null constraint equations. For transformations that act covariantly on the boundary structures, the Brown-York charges coincide with canonical charges constructed from a version of the Wald-Zoupas procedure. For anomalous transformations, the charges differ by an intrinsic functional of the boundary geometry, which we explicity verify for a set of symmetries associated with finite null hyper-surfaces. Applications of the null Brown-York stress tensor to symmetries of asymptotically flat spacetimes and celestial holography are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗