Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matched”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Servo-aero-gravo-elastic (SAGE) scaling and its application to a 13-MW downwind turbine

Reduced scale wind turbines can be extremely cost-effective to test new rotor concepts since prototype costs are heavily dependent on the rotor diameter. Ideally, the scaled model would have the same non-dimensional deflections, dynamics, and control behavior as the full-scale model. This would provide a high-fidelity demonstration of the full-scale performance, which is ideal if the full-scale turbine has significant aeroelastic interactions. To this end, servo-aero-gravo-elastic (SAGE) scaling is developed and applied to a 13-MW turbine that is scaled to a 20% scale model. The scaling preserves the tip-speed ratio, the rotor speed normalized by the flapping frequency, and the tip deflections normalized by the blade length. In addition, the controller employs the same control structure (gain-scheduled pitch control and variable speed torque control) and is scaled dynamically (e.g., matching non-dimensional time constant of the pitch angle, etc.). Furthermore, the thrust, gravity, and centrifugal moments are scaled such that the load angles are preserved as a function of a non-dimensional wind speed. However, the environmental scaling must consider differences in Reynolds number (since this parameter cannot be held constant) and subsequent changes in the axial induction factor. While the presented results showcase these differences during operational conditions, the non-dimensional tip deflections remain comparable through all wind speed ranges, indicating the viability of the SAGE scaling method in matching full-scale aeroelastic responses.

Kaminski, Meghan (ORCID:0000000322334784)↗
DOE OSTI2022

Shock-driven dispersal of a corrugated finite-thickness particle layer

A research area emerging in the multiphase flow community is the study of shock-driven multiphase instability (SDMI), a gas–particle analog of the traditional fluid-fluid Richtmyer–Meshkov instability (RMI). In this work, we study the interaction of planar air shocks with corrugated glass particle curtains through the use of numerical simulations with an Eulerian–Lagrangian approach. One objective of this study is to compare the simulated particle curtains to a comparable set of shock tube experiments performed to analyze traditional RMI of a gas curtain. The simulations are set to match the experimental shock Mach numbers and perturbation wavelengths (3.6 and 7.2 mm) while also matching the Atwood number of the experiments to the multiphase Atwood number of the simulations. Varying particle diameters are tested in the simulations to explore the impact of particle diameter on the evolution of the particle curtain. This simulation setup allows for a one-to-one comparison between RMI and SDMI under comparable conditions while also allowing for a separate study into the validity of the use of the multiphase Atwood number to compare the single-phase and multiphase instabilities. In particular, we show that the comparison depends on the diameter of the particles (thus, dependent on the Stokes number of the flow). A second objective of this study is to analyze the effect of the initial particle volume fraction on the evolution of the curtain and the behavior of the instability. This is done through analyzing the effect of the multiphase terms of the vorticity evolution equation on the vorticity deposition in SDMI. Also discussed is the effect of the particle diameter on the multiphase generation terms as well as in the baroclinic vorticity generation term in SDMI as the shock passes over the curtain.

42 ENGINEERING↗
DOE OSTI2022

Yield degradation mechanisms for two-shock capsules evaluated through simulations

An investigation of twenty two-shock campaign indirectly driven capsules on the National Ignition Facility was conducted using the xRAGE computer code. The two-shock platform was developed to look at the sensitivity of fuel–ablator mix with shock timing, asymmetry, surface roughness, and convergence on roughly ignition size scale capsules. This platform used CH/CD (plastic/deuterated plastic) shell capsules that were about 685-μm outer radius and filled with D2 or hydrogen-tritium (HT) gas. The experimental radius and velocity vs time, neutron yield, burn averaged ion temperature (Tion), burn width, and self-emission image size were compared to one-dimensional (1D) and two-dimensional (2D) simulations. Our 2D simulations suggest that the mixing of glass from the fill tube was the dominant source of impurity in the gas region of the capsule during burn, along with fuel–ablator mix. The mass of glass mixed in is about 5–10 ng. Our 2D simulations capture most of the yield trends from different degradation mechanisms, and they match the observed burn width and Tion measurements. Our 2D models match all the available data to within 2.5 times the normalized experimental error for 19 of 20 capsules.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Fast uncertainty estimates in deep learning interatomic potentials

Deep learning has emerged as a promising paradigm to give access to highly accurate predictions of molecular and material properties. A common short-coming shared by current approaches, however, is that neural networks only give point estimates of their predictions and do not come with predictive uncertainties associated with these estimates. Existing uncertainty quantification efforts have primarily leveraged the standard deviation of predictions across an ensemble of independently trained neural networks. This incurs a large computational overhead in both training and prediction, resulting in order-of-magnitude more expensive predictions. Here, we propose a method to estimate the predictive uncertainty based on a single neural network without the need for an ensemble. This allows us to obtain uncertainty estimates with virtually no additional computational overhead over standard training and inference. We demonstrate that the quality of the uncertainty estimates matches those obtained from deep ensembles. We further examine the uncertainty estimates of our methods and deep ensembles across the configuration space of our test system and compare the uncertainties to the potential energy surface. Finally, we study the efficacy of the method in an active learning setting and find the results to match an ensemble-based strategy at order-of-magnitude reduced computational cost.

Chemistry↗
DOE OSTI2023

Theoretical investigation of the triggering of neoclassical tearing modes by transient resonant magnetic perturbations in NSTX

The extended perturbed equilibrium code asymptotic matching code is used to simulate the triggering of n = 1 neoclassical tearing modes (NTMs) by a pulsed, rotating, n = 1, resonant magnetic perturbation (RMP) in two example NSTX discharges. Although the two discharges are significantly different from one another, the results of the two sets of simulations are quite similar. The critical n = 1 RMP pulse amplitude required to trigger an n = 1 NTM is minimized when the RMP pulse rotation frequency matches the linear natural frequency of an n = 1 tearing mode, resonant within the plasma, that is metastable to an NTM. However, if there is a frequency mismatch, then the seed magnetic island chain driven at the relevant resonant surface is forced to rotate with respect to the RMP, because the RMP pulse amplitude is nowhere near sufficient to lock the island chain to the RMP. This rotation causes the critical RMP pulse amplitude required to trigger an NTM to oscillate as the RMP pulse duration is varied. Furthermore, the critical amplitude is minimized when the RMP pulse duration is such that seed island chain executes a half-integer number of rotations with respect to the pulse. All of the minima have the same value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Identification and characterization of deep nitrogen acceptors in β-Ga 2 O 3 using defect spectroscopies

The ability to achieve highly resistive beta-phase gallium oxide (β-Ga 2 O 3 ) layers and substrates is critical for β-Ga 2 O 3 high voltage and RF devices. To date, the most common approach involves doping with iron (Fe), which generates a moderately deep acceptor-like defect state located at E C -0.8 eV in the β-Ga 2 O 3 bandgap. Recently, there has been growing interest in alternative acceptors, such as magnesium (Mg) and nitrogen (N), due to their predicted deeper energy levels, which could avoid inadvertent charge modulation during device operation. In this work, a systematic study that makes direct correlations between the introduction of N using ion implantation and the observation of a newly observed deep level at E C -2.9 eV detected by deep-level optical spectroscopy (DLOS) is presented. The concentration of this state displayed a monotonic dependence with N concentration over a range of implant conditions, as confirmed by secondary ion mass spectrometry (SIMS). With a near 1:1 match in absolute N and E C -2.9 eV trap concentrations from SIMS and DLOS, respectively, which also matched the measured removal of free electrons from capacitance-voltage studies, this indicates that N contributes a very efficiently incorporated compensating defect. Density functional theory calculations confirm the assignment of this state to be an N (0/-1) acceptor with a configuration of N occupying the oxygen site III [N O(III) ]. The near ideal efficiency for this state to compensate free electrons and its location toward the midgap region of the β-Ga 2 O 3 bandgap demonstrates the potential of N doping as a promising approach for producing semi-insulating β-Ga 2 O 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2023

Validity condition for the local sheath impedance boundary condition and a nonlocal generalization

ICRF sheaths can cause unwanted interactions of high-power RF waves with material surfaces in magnetic fusion devices. In previous work, a local RF sheath impedance boundary condition (BC) was derived for use in ICRF codes together with a microscale (i.e., Debye or sheath width scale) model for obtaining the sheath impedance used in that BC. This local RF sheath BC matches the normal component of current and electrostatic potential across the sheath-plasma interface. Collapsing the matching conditions at the sheath-plasma interface to a BC depends on the assumption of scale separation, which can be violated when conditions along the local radius of curvature of the surface vary sufficiently rapidly. The validity condition is explored in this contribution, with special attention to the case where the magnetic field approaches being tangent to the surface. When the local sheath BC no longer applies, a non-local sheath BC is developed under the assumption of a more relaxed scale separation assumption. It is shown that the non-local sheath BC reduces to the previous local sheath BC under appropriate conditions. Furthermore, a surface-integrated sheath admittance parameter describes the 2D physics in the new BC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simultaneous reproduction of experimental profiles, fluxes, transport coefficients, and turbulence characteristics via nonlinear gyrokinetic profile predictions in a DIII-D ITER similar shape plasma

Experimental conditions obtained on the DIII-D tokamak in the ITER Similar Shape (ISS) have been compared extensively with nonlinear gyrokinetic simulation using the CGYRO code [Candy et al., J. Comput. Phys. 324, 73–93 (2016)] with comparisons spanning ion and electron heat fluxes, electron and impurity particle transport, and turbulent fluctuation levels and characteristics. Bayesian optimization techniques [Rodriguez-Fernandez et al., Nucl. Fusion 62(7), 076036 (2022)], combined with nonlinear gyrokinetics, have been used to obtain simultaneously Qi, Qe, and Γe flux-matched profiles that are found to be in good agreement with experimental profile measurements. Synthetic diagnostics were used to compare measured beam emission spectroscopy and correlation electron cyclotron emission turbulent fluctuations with nonlinear simulation. Although some disagreements exist, nonlinear simulations are found to be in generally good agreement with measured fluctuation levels, spectral shapes, and measured radial trends in low-k δne/ne and δTe/Te. Low (Li and C) and mid-Z (Ca) impurity transport was also compared with these flux-matched simulations. Fully stripped, low-Z impurities are well reproduced by the gyrokinetic modeling while clear disagreement exists in comparisons with mid-Z impurities. Nonlinear gyrokinetic investigation into the Z dependence of impurity transport in the ISS conditions is also performed, demonstrating clear trends of impurity diffusion with impurity Z (both D∝Z and D∝1/Z) that vary with the radial location studied. These trends are shown to result from the local dominance of ion temperature gradient or ∇n driven trapped electron mode turbulence and may contribute to the disagreement between simulation and experiment in mid-Z impurity transport. The results of this work represent one of the most complete validation studies of the gyrokinetic model performed to date and provide an example of new capabilities for predicting performance in future fusion devices.

Physics↗
DOE OSTI2024

Reconciling calculations and measurements of inverse bremsstrahlung absorption

It was recently shown that the use of Coulomb logarithms appropriate for bremsstrahlung radiation (rather than transport processes) along with corrections for the Langdon effect and ion screening reproduced measurements of collisional absorption in well-characterized underdense plasmas [D. Turnbull et al., Phys. Rev. Lett. 130, 145103 (2023)]. However, it was recognized at the time that the use of the standard absorption-reduction factor from Langdon's seminal paper was inconsistent with the use of Coulomb logarithms that are thermally averaged over a Maxwellian electron-velocity distribution function. A more accurate approach would be to average over the expected super-Gaussian distribution function while accounting for the Gaunt factor's velocity dependence, which somewhat mitigates the Langdon effect; however, at that time, this theory matched the data less well. This conflict is now eliminated with the additional insight that the ionization state of our mid-Z ion species (when present) was lower than had been assumed, as evidenced by the Thomson-scattering data and time-dependent Cretin simulations. We are now able to show that an improved treatment of the Langdon effect provides the best match to data. Otherwise, the prior conclusions remain unchanged. We also show an example of the substantial expected impact to the absorption rate in calculations of indirect-drive hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Interpolation of compound semiconductor alloy parameters from those of their constituents

Several methods have been proposed for interpolation of the value of physical parameters of quaternary alloys from those of their constituent ternary and binary sub-alloys. These expressions agree when non-linear bowing terms are not required; they differ in how the bowing terms of the bounding ternaries should be utilized. Common interpolation expressions for quaternaries can be generalized into two groups: (1) those that use a linear interpolation of the nearest ternary parameter values and (2) those that interpolate over binary values with a bowing term derived from the bounding ternaries. The second group of methods is equivalent to a polynomial expansion over the alloy’s interpolation space. For compound semiconductor alloys, the geometry of the composition space is the direct sum of the group-III and group-V mixture sub-spaces. The mixture sub-spaces are best described using barycentric coordinates on a regular simplex. A general polynomial expansion of the value of an alloy parameter using barycentric coordinates for the group-III and group-V simplex spaces is described along with an algorithm to generate interpolation expressions for alloys with arbitrary numbers of elements, including quinary and senary alloys. It is shown that a polynomial expansion produces values in closer agreement with the direct gap of quaternaries lattice-matched to common substrates than do approaches using an interpolation of the ternary values, despite a prominent recommendation to the contrary. Finally, a quaternary correction term is described that improves the predicted direct bandgap energies of GaInAsSb for compositions near those lattice matched to InP, InAs, and GaSb.

Olesberg, Jonathon T. [Sandia National Laboratorie↗
DOE OSTI2024

Three-wave coupling observed between a shear Alfvén wave and a kink-unstable magnetic flux rope

Results from a laboratory experiment are presented in which, for the first time, a shear Alfvén wave is launched using an antenna in a current-carrying plasma column that is tailored to be either stable or unstable to the kink oscillation. As the plasma is driven kink unstable, the frequency power spectrum of the Alfvén wave evolves from a single peak to a peak with multiple sidebands separated by integer multiples of the kink frequency. The main sidebands (one on either side of the launched wave peak in the power spectrum) are analyzed using azimuthal wavenumber matching, perpendicular and parallel wavenumber decomposition, and bispectral time series analysis. The dispersion relation and three-wave matching conditions are satisfied, given each sideband is a propagating Alfvén wave that results from the interaction of the pump Alfvén wave and the co-propagating component of a half-wavelength, standing kink mode. The interaction is shown to generate smaller perpendicular wavelength Alfvén waves that drive energy transport to scales that will approach the dissipation scale of k⊥ρs=1, with k⊥ being the perpendicular wavenumber and ρs being the ion gyroradius at the electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Neutron-producing gas puff Z-pinch experiments on a fast, low-impedance, 0.5 MA linear transformer driver

A study on the neutron production from single and double gas puff Z-pinches on the CESZAR linear transformer driver with ~0.45 MA current and 170 ns rise time is presented. Total neutron yield measurements made with a LaBr activation detector are compared for three configurations, using a double nozzle setup. When a single, hollow, deuterium gas shell was used, reliable implosions could only be attained at higher load mass than the optimal value to match implosion time with the driver rise time, with neutron yields of ~10 6 per pulse. The use of a double gas puff configuration with a deuterium center jet allowed a reduction in the shell density and operation closer to machine-matched conditions, recording up to (4.1 ± 0.3) × 10 7 neutrons/pulse when either Kr or D 2 was used in the shell. For a comparable mass and implosion time, using a higher atomic-number gas in the outer shell results in more unstable plasma surface and smaller plasma radius at the location of instability bubbles, which, however, do not seem to consistently correlate with a higher neutron yield. Comparing implosion dynamics with models and neutron yields with literature scaling suggests that the machine current is not well coupled to the plasma during the final stages of compression. Optimizing current and energy coupling to the pinched plasma is critical to improving performance, particularly in low-impedance drivers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Results from a synthetic model of the ITER XRCS-Core diagnostic based on high-fidelity x-ray ray tracing

A high-fidelity synthetic diagnostic has been developed for the ITER core x-ray crystal spectrometer diagnostic based on x-ray ray tracing. This synthetic diagnostic has been used to model expected performance of the diagnostic, to aid in diagnostic design, and to develop engineering tolerances. The synthetic model is based on x-ray ray tracing using the recently developed xicsrt ray tracing code and includes a fully three-dimensional representation of the diagnostic based on the computer aided design. The modeled components are: plasma geometry and emission profiles, highly oriented pyrolytic graphite pre-reflectors, spherically bent crystals, and pixelated x-ray detectors. Plasma emission profiles have been calculated for Xe 44+ , Xe 47+ , and Xe 51+ , based on an ITER operational scenario available through the Integrated Modelling & Analysis Suite database, and modeled within the ray tracing code as a volumetric x-ray source; the shape of the plasma source is determined by equilibrium geometry and an appropriate wavelength distribution to match the expected ion temperature profile. All individual components of the x-ray optical system have been modeled with high-fidelity producing a synthetic detector image that is expected to closely match what will be seen in the final as-built system. Particular care is taken to maintain preservation of photon statistics throughout the ray tracing allowing for quantitative estimates of diagnostic performance.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Electronic properties of corundum-like Ir 2 O 3 and Ir 2 O 3 -Ga 2 O 3 alloys

In the hexagonal, corundum-like structure, α-Ga 2 O 3 has a bandgap of ~5.1 eV, which, combined with its relatively small electron effective mass, high Baliga's figure of merit, and high breakdown field, makes it a promising candidate for power electronics. Ga 2 O 3 is easy to dope n-type, but impossible to dope p-type, impeding the realization of some electronic device designs. Developing a lattice-matched p-type material that forms a high-quality heterojunction with n-type Ga 2 O 3 would open new opportunities in electronics and perhaps optoelectronic devices. In this work, we studied Ir 2 O 3 as a candidate for that purpose. Using hybrid density functional theory calculations we predict the electronic band structure of α-Ir 2 O 3 and compare that to α-Ga 2 O 3 , and study the stability and electronic properties of α-(Ir x Ga 1–x ) 2 O 3 alloys. We discuss the band offset between the two materials and compare it with recently available experimental data. We find that the Ir d bands that compose the top of the valence band in α-Ir 2 O 3 are much higher in energy than O p bands in α-Ga 2 O 3 , possibly enabling effective p-type doping. Finally, our results provide an insight into using the Ir 2 O 3 or Ir 2 O 3 -Ga 2 O 3 alloys as p-type material lattice-matched to α-Ga 2 O 3 for the realization of p–n heterojunctions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A comparison of probabilistic generative frameworks for molecular simulations

Generative artificial intelligence is now a widely used tool in molecular science. Despite the popularity of probabilistic generative models, numerical experiments benchmarking their performance on molecular data are lacking. Here, in this work, we introduce and explain several classes of generative models, broadly sorted into two categories: flow-based models and diffusion models. We select three representative models: neural spline flows, conditional flow matching, and denoising diffusion probabilistic models, and examine their accuracy, computational cost, and generation speed across datasets with tunable dimensionality, complexity, and modal asymmetry. Our findings are varied, with no one framework being the best for all purposes. In a nutshell, (i) neural spline flows do best at capturing mode asymmetry present in low-dimensional data, (ii) conditional flow matching outperforms other models for high-dimensional data with low complexity, and (iii) denoising diffusion probabilistic models appear the best for low-dimensional data with high complexity. Our datasets include a Gaussian mixture model and the dihedral torsion angle distribution of the Aib9 peptide, generated via a molecular dynamics simulation. We hope our taxonomy of probabilistic generative frameworks and numerical results may guide model selection for a wide range of molecular tasks.

Artificial intelligence↗
DOE OSTI2025

Evaluating nonlocal heat transport in directly driven chromium spheres using x-ray spectroscopy

We report on experiments investigating heat transport in laser-generated plasmas using directly driven chromium spheres. The spheres are fielded at the OMEGA laser facility and are driven with laser intensities of 5×10 14 Wcm −2 . Plasma conditions in the corona and scattered light are measured experimentally and compared against predictions from two-dimensional (2D) radiation-hydrodynamic simulations using different heat transport models. Spectroscopic analysis of x-ray self-emission is used as an additional diagnostic. X-ray emission is integrated over a large region of the plasma, probing regions that are not observed by localized optical Thomson scattering. In particular, x-ray emission peaks near the plasma critical density, so emission from optically thin lines provides information on plasma conditions where nonlocal transport is most likely to be significant. Three common heat transport models are considered: local transport with flux limiters f = 0.15 and f = 0.03, and the nonlocal Schurtz–Nicolai–Busquet (SNB) model. Consistent with previous work, both the high-flux (f = 0.15) and SNB models show good agreement with experimentally measured plasma conditions in the corona despite overpredicting laser absorption, whereas the low-flux (f = 0.03) model fails to match any experimental data. Conditions inferred from x-ray self-emission line ratios support this conclusion during the period of laser peak power, although synthetic spectra for all models fail to match the experiment during the transient portions of the pulse. For these reasons, the low-flux model is again rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Calibration of a broadband x-ray crystal spectrometer using a continuum x-ray source and photon-counting detectors

We report measurements of sensitivity of a broadband (≃20–30 keV) x-ray crystal spectrometer using a bremsstrahlung continuum x-ray source and several detectors, including an energy-discriminating photon-counting point detector (Si drift diode, Amptek Inc.), imaging hybrid photon-counting detectors (Eiger2-Si and Eiger2-CdTe, DECTRIS Ltd.), and image plates (SR type). Sensitivity is defined as a ratio of the crystal-reflected energy-position-dispersed spectrum measured by a given detector to the incident non-dispersed spectrum measured using the energy-discriminating point detector. The sensitivity derived from data measured exclusively by the point detector is considered detector-independent and serves as a reference. The sensitivities derived from data collected with the hybrid photon-counting detectors were matched to the reference using a scale factor of 1.05. The image plate-derived sensitivities required a scale factor of 1.16 to match the reference. Furthermore, the resulting mismatch in the shapes of all measured scaled sensitivities in the range of the spectrometer was ≲±10%, while the mismatch between the shapes of the sensitivities corresponding to the Si-based detectors was ≲±2.5%.

Crystal optics↗