Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-element”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Amplified spontaneous emission from europium-based molecular complexes coupled to photonic crystal cavities

Rare-earth ion-based materials bear many remarkable optical properties that render them highly appealing for lighting and quantum-related applications. However, their small oscillator strength and weak emission often pose limitations. Here, in this work, we synthesize and couple Eu(III)-based molecular complexes to nanobeam photonic crystals supporting air modes. A reasonable spatial overlap between the molecular complexes and cavity modes leads to an average spontaneous emission coupling efficiency of 0.19. Our pump power-dependent photoluminescence measurements evidence amplified spontaneous emission from the molecular complexes with an amplification threshold as low as 4.4 W/cm 2 , likely benefiting from the efficient coupling. These findings suggest that integrating rare-earth ion-based molecular complexes with photonic structures could be a viable approach for regulating their emission characteristics for particular applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Application of the van der Pauw method for electrical conductivity measurements at high temperatures using an insulating compressing ring

High quality data on the high temperature electrical properties of ceramics, particularly oxides, is of great value for material selection, design, and modeling for a broad range of emerging applications. Utilizing the mismatch in the coefficient of thermal expansion between two materials, a purely mechanical method for establishing electrical contact in the van der Pauw geometry to measure the bulk resistivity of ceramic disks at high temperatures is presented. Measurements of a reference material, 20 mol. % Gd-doped cerium oxide, are presented up to 1000 °C. The viability of electrical measurements up to a maximum temperature of 1600 °C is also considered. Measurements are performed using multiple techniques and compared to literature values finding excellent agreement. Finally, the approach described in this work enables the van der Pauw method to be applied to many ceramic materials over a wide range of temperatures and environments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Measurement of helicon waves with phase contrast imaging on DIII-D – A theoretical feasibility study

A DIII-D high-beta H-mode discharge, with I p = 850 kA, B t = 2.1 T n e = 4 10 19 m –3 , has been designed to validate full wave modeling of helicon waves by optimizing the expected response of the Phase Contrast Imaging diagnostic. Helicon waves have been predicted to have high current drive efficiency off-axis without facing the accessibility issues of lower-hybrid waves. To test these predictions experimentally, DIII-D has recently commissioned a high-power helicon antenna. To confidently predict the behavior of helicon waves in future devices, measurements of their fundamental properties and validation against models will be essential. Phase contrast imaging (PCI) is an absolutely calibrated internal reference interferometer able to measure density fluctuations with radial wavenumbers k R between 1.5 cm –1 and 20 cm –1 . For helicon waves 2 cm –1 < k R < 10 cm –1 is expected, allowing PCI to measure their envelope and wavenumber spectrum. This makes PCI a powerful tool for the validation of state-of-the-art models, like the AORSA full wave code. AORSA is used to compute the density perturbations measured by the PCI with 2D calculations corresponding to 11 different toroidal mode numbers combined to resolve the trajectory of the helicon wave in 3D. This is necessary because the waves travel 120 degrees toroidally from the antenna to the PCI. Here, a cold plasma finite element model (CPFEM) [4] is used to predict propagation through the scrape-off layer. The result of the CPFEM model is connected to AORSA by creating an artificial Gaussian antenna on the last closed flux surface. PCI shows best results for waves with small vertical wavenumbers k z . Modeling the helicon waves for several past DIII-D experiments shows that k z is minimized if the intersection of the helicon and the PCI laser beams occurs in the midplane. For such an optimized scenario the predicted signal level is two orders of magnitude larger than the background density fluctuations arising from broadband turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Measurement bias in self-heating x-ray free electron laser experiments from diffraction studies of phase transformation in titanium

X-ray self-heating is a common by-product of X-ray Free Electron Laser (XFEL) techniques that can affect targets, optics, and other irradiated materials. Diagnosis of heating and induced changes in samples may be performed using the x-ray beam itself as a probe. However, the relationship between conditions created by and inferred from x-ray irradiation is unclear and may be highly dependent on the material system under consideration. Here, we report on a simple case study of a titanium foil irradiated, heated, and probed by a MHz XFEL pulse train at 18.1 keV delivered by the European XFEL using measured x-ray diffraction to determine temperature and finite element analysis to interpret the experimental data. We find a complex relationship between apparent temperatures and sample temperature distributions that must be accounted for to adequately interpret the data, including beam averaging effects, multivalued temperatures due to sample phase transitions, and jumps and gaps in the observable temperature near phase transformations. The results have implications for studies employing x-ray probing of systems with large temperature gradients, particularly where these gradients are produced by the beam itself. Finally, this study shows the potential complexity of studying nonlinear sample behavior, such as phase transformations, where biasing effects of temperature gradients can become paramount, precluding clear observation of true transformation conditions.

Crystallography↗
DOE OSTI2025

X-point effects on the ideal MHD modes in tokamaks in the description of dual-poloidal-region safety factor

The flux coordinates with dual-region safety factor (q) in the poloidal direction are developed in this work. The X-point effects on the ideal MHD modes in tokamaks are then analyzed using this coordinate system. Since the X-point effects mainly affect the edge region, the modes localized at the tokamak edge are particularly examined. Two types of modes are studied. The first is related to the conventional peeling or peeling-ballooning modes. The mode existence aligned with the local magnetic field in the poloidally core region as observed experimentally is confirmed. The X points are shown to contribute to a stabilizing effect for the conventionally treated modes with the surface-averaged q and with the tokamak edge portion truncated. The other is the axisymmetric modes localized in the vicinity of X points, which can affect the cross-field-line transport near the X points. The existence of axisymmetric modes points to the possibility of applying a toroidally axisymmetric resonant magnetic perturbation (RMP) in the X-point area for mitigating the edge localized modes, which can be an alternative to the current RMP design. The dual q description also has important implications for the existing non-axisymmetric RMP concept. It helps to understand why the RMP suppression of edge localized modes is difficult to achieve in the double-null tokamak configurations and points to the possibility of further improving the current RMP concept by considering the alignment to the local q.

Fourier analysis↗
DOE OSTI2024

A practical approach to calculating magnetic Johnson noise for precision measurements

Magnetic Johnson noise is an important consideration for many applications involving precision magnetometry, and its significance will only increase in the future with improvements in measurement sensitivity. The fluctuation–dissipation theorem can be utilized to derive analytic expressions for magnetic Johnson noise in certain situations, but when used in conjunction with finite element analysis tools, the combined approach is particularly powerful as it provides a practical means to calculate the magnetic Johnson noise arising from conductors of arbitrary geometry and permeability. In this paper, we demonstrate this method to be one of the most comprehensive approaches presently available to calculate thermal magnetic noise. In particular, its applicability is shown to not be limited to cases where the noise is evaluated at a point in space but also can be expanded to include cases where the magnetic field detector has a more general shape, such as a finite-size loop, a gradiometer, or a detector that consists of a polarized atomic species trapped in a volume. Furthermore, some physics insights gained through studies made using this method are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
DOE OSTI2025

A three-dimensional laser ray-tracing methodology for radiation-hydrodynamics simulations

We report on a methodology for performing laser ray-tracing in three spatial dimensions for radiation-hydrodynamics simulation codes. Our method, which is an extension of that developed in Haines et al., Comput. Fluids 201, 104478 (2020), utilizes an automatically generated separate mesh for the laser ray-tracing from the radiation-hydrodynamics mesh. This enables the laser mesh to be tailored to minimize ray noise with significantly fewer rays than would be required when the ray-tracing is performed on the radiation-hydrodynamics mesh, primarily by allowing the use of high-aspect-ratio cells that are not suitable for hydrodynamics solvers. For a planar target, we show that our method provides a ≈ 100× reduction in computational expense to achieve a fixed level of ray noise relative to ray-tracing directly on the radiation-hydrodynamics mesh. The relatively low ray requirement also enables efficient computation of cross-beam energy transfer. Each cell in the logically cubic laser mesh is a non-convex dodecahedron with triangular sides, and numerical integration of the ray trajectories and inverse bremsstrahlung is performed by mapping each cell to the unit cube. We will describe our methodology in detail as well as its implementation in the xRAGE radiation-hydrodynamics code, discuss performance, and present the results from applying the methodology to test problems with analytic solutions for laser ray-tracing through a quadratic density gradient with an analytic solution as well as for a laser-driven heat front. In 3D radiation-hydrodynamics simulations of laser-driven experiments performed on the National Ignition Facility, laser ray-tracing with our methodology uses less than 1% of total computational time while introducing acceptably low levels of ray noise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Assessment of errors in analytic modeling of permanent magnet electron spectrometers for laser-plasma accelerators

Here, we assess the error in three treatments of a critical electron beam diagnostic for laser-plasma acceleration (LPA) experiments: a permanent magnet electron spectrometer. Since LPA electron spectrometers are often difficult to calibrate due to the scarcity of well-characterized, tunable electron beam sources in the appropriate energy range and the mechanical complexity of electron spectrometers, the standard of calibration is Hall probe measurements. We first compare the electron spectrometer performance between SIMION calculations and the Hall probe measurements. We find up to a 60% (<8%) error in determining the absolute energy for electrons below (above) 80 MeV when using the SIMION-modeled magnetic field vs the measured magnetic field due to SIMION overpredicting fringe field strength. The difference in spectrometer energy resolution is within ∼10% between the modeled and measured fields. We then assess a simple block model, which is commonly used. We find that the simple block model presented here sufficiently captures the predicted as-built spectrometer performance for any application provided that the user can tolerate uncertainties on absolute energy determination up to 5% and errors on energy resolution up to 1%.

Electromagnetic optics↗
DOE OSTI2026

Quantum speed limit for the out-of-time-ordered correlator from an open-system perspective

Scrambling, the delocalization of initially localized quantum information, is commonly characterized by the out-of-time-ordered correlator (OTOC). Employing the OTOC–Renyi-2 entropy theorem, we derive a quantum speed limit for the OTOC, which sets a lower bound for the rate with which information can be scrambled. This bound becomes particularly tractable by describing the scrambling of information in a closed quantum system as an effective decoherence process of an open system interacting with an environment. We prove that decay of the OTOC can be bounded by the strength of the system-environment coupling and two-point environmental correlation functions. We validate our analytic bound numerically using the nonintegrable transverse field Ising model. Furthermore, our results provide a universal and model-agnostic quantitative framework for understanding the dynamical limits of information spreading across quantum many-body physics, condensed matter systems, and engineered quantum platforms.

Fermions↗
DOE OSTI2026

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Thermal and electric field driven rf breakdown precursor formation on metal surfaces

The phenomenon of electric breakdown poses serious challenges to the design of devices that operate in high electric field environments. Experimental evidence points toward breakdown events that are accompanied by elevated temperatures and dark current spikes, which is attributed to high-asperity nanostructure formation that enhances the local electric field and triggers a runaway process. However, the exact mechanistic origin of such nanostructures under typical macroscopic operational conditions of electric field and magnetic-field-mediated heating remains poorly understood. In this work, we simulate the evolution of a copper surface under the combined action of the electric fields and elevated temperatures. Using a mesoscale curvature-driven surface evolution model, we show how a copper surface can undergo a type of dynamical instability that naturally leads to the formation of sharp asperities in realistic experimental conditions. Exploring the combined effect of fields and temperature rise, we identify the critical regimes that allow for the formation of breakdown precursors. The results show that thermoelastic stresses, while not essential, can significantly lower the critical electric field required for runaway surface instability, which is consistent with experimental observations that thermal effects can increase breakdown rates. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Comprehensive model for evaluating voltage losses and performance improvements in thin-film photovoltaic devices

Progress of state-of-the-art and next-generation thin-film photovoltaic devices is often stymied by open-circuit voltage (𝑉 oc ) that is significantly lower than theoretical and practical limits. Yet, effectively diagnosing the primary sources of voltage loss remains challenging. Herein, a sequence of device-level characterization techniques and simulations are employed to identify and rank loss mechanisms. For the research-based Cd⁡(Se,Te) device under study, most of the loss was at the front semiconductor heterointerface due to a clifflike conduction-band offset that lowered the recombination activation energy. Additional losses due to band tails were quantified by photoluminescence analysis. The latter provided the absorption coefficient and activation energy reduction associated with band tails as inputs to device models. Simulations showed that alleviating front-interface issues would improve 𝑉 oc , but it would then be limited by bulk recombination. Further improvement of the bulk would then lead to back-contact limitations. Reducing band tails is beneficial in any circumstance. In conclusion, this analysis provides guidance for reaching toward the radiative 𝑉 oc limit.

14 SOLAR ENERGY↗
DOE OSTI2024

Light-induced transitions of polar state and domain morphology of photoferroelectric nanoparticles

Using the Landau-Ginzburg-Devonshire approach, we study light-induced phase transitions, evolution of polar state and domain morphology in photo-ferroelectric nanoparticles (NPs). Light exposure increases the free carrier density near the NP surface and may in turn induce phase transitions from the nonpolar paraelectric to the polar ferroelectric phase. Using the uniaxial photo-ferroelectric Sn 2 P 2 S 6 as an example, here we show that visible light exposure induces the appearance and vanishing of striped, labyrinthine or curled domains and changes in the polarization switching hysteresis loop shape from paraelectric curves to double, pinched and single loops, as well as the shifting in the position of the tricritical point. Furthermore, we demonstrate that an ensemble of non-interacting photo-ferroelectric NPs may exhibit superparaelectric-like features at the tricritical point, such as strongly frequency-dependent giant piezoelectric and dielectric responses, which can potentially be exploited for piezoelectric applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Bipartite mutual information in classical many-body dynamics

Information theoretic measures have helped to sharpen our understanding of many-body quantum states. As perhaps the most well-known example, the entanglement entropy (or more generally, the bipartite mutual information) has become a powerful tool for characterizing the dynamical growth of quantum correlations. By contrast, although computable, the bipartite mutual information (MI) is almost never explored in classical many particle systems; this owes in part to the fact that computing the MI requires keeping track of the evolution of the full probability distribution, a feat which is rarely done (or thought to be needed) in classical many-body simulations. Here, we utilize the MI to analyze the spreading of information in 1D elementary cellular automata (CA). Broadly speaking, we find that the behavior of the MI in these dynamical systems exhibits a few different types of scaling that roughly correspond to known CA universality classes. Of particular note is that we observe a set of automata for which the MI converges parametrically slowly to its thermodynamic value. We develop a microscopic understanding of this behavior by analyzing a two-species model of annihilating particles moving in opposite directions. Furthermore, our work suggests the possibility that information theoretic tools such as the MI might enable a more fine-grained characterization of classical many-body states and dynamics.

Cellular automata↗
DOE OSTI2023

Crystal plasticity finite element simulation of lattice rotation and x-ray diffraction during laser shock compression of tantalum

Here we present a crystal plasticity model tailored for high-pressure, high-strain-rate conditions that uses a multiscale treatment of dislocation-based slip kinetics. We use this model to analyze the pronounced plasticity-induced lattice rotations observed in shock-compressed polycrystalline tantalum via in situ x-ray diffraction. By making direct comparisons between experimentally measured and simulated texture evolution, we can explain how the details of the underlying slip kinetics control the degree of lattice rotation that ensues. Specifically, we show that only the highly nonlinear kinetics caused by dislocation nucleation can explain the magnitude of the rotation observed under shock compression. We demonstrate a good fit between our crystal plasticity model and x-ray diffraction data and exploit the data to quantify the dislocation nucleation rates that are otherwise poorly constrained by experiment in the dynamic compression regime.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Measurement of Photovoltaic Module Deformation Dynamics During Hail Impact Using Digital Image Correlation

Stereo high-speed video of photovoltaic modules undergoing laboratory hail tests was processed using digital image correlation to determine module surface deformation during and immediately following impact. The purpose of this work was to demonstrate a methodology for characterizing module impact response differences as a function of construction and incident hail parameters. Video capture and digital image analysis were able to capture out-of-plane module deformation to a resolution of ±0.1 mm at 11 kHz on an in-plane grid of 10 × 10 mm over the area of a 1 × 2 m commercial photovoltaic module. With lighting and optical adjustments, the technique was adaptable to arbitrary module designs, including size, backsheet color, and cell interconnection. Furthermore, impacts were observed to produce an initially localized dimple in the glass surface, with peak deflection proportional to the square root of incident energy. Subsequent deformation propagation and dissipation were also captured, along with behavior for instances when the module glass fractured. Natural frequencies of the module were identifiable by analyzing module oscillations postimpact. Limitations of the measurement technique were that the impacting ice ball obscured the data field immediately surrounding the point of contact, and both ice and glass fracture events occurred within 100 μs, which was not resolvable at the chosen frame rate. Increasing the frame rate and visualizing the back surface of the impact could be applied to avoid these issues. Applications for these data include validating computational models for hail impacts, identifying the natural frequencies of a module, and identifying damage initiation mechanisms.

14 SOLAR ENERGY↗