Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Definition”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Depletion Benchmark Analysis on a Lead Fast Reactor Using PyARC/OpenMC

PyARC is a user-friendly fast reactor analysis tool that automates multiphysics workflows using the “extended suite” of Argonne Reactor Computation (ARC) codes by providing a single common input for model definition, code execution, and output post-processing. A lead fast reactor (LFR) benchmark model is used to perform depletion calculations using the newly integrated OpenMC depletion capability in PyARC, building on previous analysis using the ARC codes through PyARC and Serpent. Results for core lifetime k-effective, shutdown decay heat, and end-of-life heavy-metal inventory are compared to verify the PyARC/OpenMC integration against the PyARC/ARC workflow and Serpent for depletion analysis of LFR designs. The results show satisfactory agreement among all three methods, with remaining discrepancies largely attributable to differences in nuclear data libraries and decay-chain modeling detail rather than to fundamental modeling limitations.

Kiesling, Kalin R.↗
DOE OSTI2024

A Second Moment Method for k -Eigenvalue Acceleration with Continuous Diffusion and Discontinuous Transport Discretizations

The second moment method is a linear acceleration technique that couples the transport equation to a diffusion equation with transport-dependent additive closures. The resulting low-order diffusion equation can be discretized independent of the transport discretization, unlike diffusion synthetic acceleration, and is symmetric positive definite, unlike quasidiffusion. While this method has been shown to be comparable to quasidiffusion in iterative performance for fixed source and time-dependent problems, it is largely unexplored as an eigenvalue problem acceleration scheme due to the belief that the resulting inhomogeneous source makes the problem ill posed. Recently, a preliminary feasibility study was performed on the second moment method for eigenvalue problems. The results suggested comparable performance to quasidiffusion and more robust performance than diffusion synthetic acceleration. This work extends the initial study to more realistic reactor problems using state-of-the-art discretization techniques. Finally, the results in this paper show that the second moment method is more computationally efficient than its alternatives on complex reactor problems with unstructured meshes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Optical Spectroscopic Investigation of Hexavalent Actinide Ions in n-Dodecane Solutions of Tri-butyl Phosphate

The extraction of hexavalent actinides An(VI) by tri-butyl phosphate (TBP) was investigated by electronic absorption and vibrational spectroscopies. Through a series of spectral subtractions, vibrational spectra associated with TBP, TBP–HNO3 adducts, and An(VI)–TBP complexes could be isolated. Investigation of U(VI) extracts indicated spectral features consistent with the formation of the expected [UO2(NO3)2(TBP)2] complex, but spectral features of other species were clearly evident. Likewise, multiple species were evident in the electronic absorption and vibrational spectra of TBP phases generated by extraction of Pu(VI). Although definitive characterization of the additional species formed could not be achieved in this work, it is hypothesized that they contain 3:1 TBP-to-An(VI) stoichiometry.

Lumetta, Gregg J.↗
DOE OSTI2022

HyKKT: a hybrid direct-iterative method for solving KKT linear systems

Here, we propose a solution strategy for the large indefinite linear systems arising in interior methods for nonlinear optimization. The method is suitable for implementation on hardware accelerators such as graphical processing units (GPUs). The current gold standard for sparse indefinite systems is the LBLT factorization where L is a lower triangular matrix and B is 1×1 or 2×2 block diagonal. However, this requires pivoting, which substantially increases communication cost and degrades performance on GPUs. Our approach solves a large indefinite system by solving multiple smaller positive definite systems, using an iterative solver on the Schur complement and an inner direct solve (via Cholesky factorization) within each iteration. Cholesky is stable without pivoting, thereby reducing communication and allowing reuse of the symbolic factorization. We demonstrate the practicality of our approach on large optimal power flow problems and show that it can efficiently utilize GPUs and outperform LBL T factorization of the full system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Examination of probability distribution of mixture fraction in LES/FDF modelling of a turbulent partially premixed jet flame

An accurate prediction of the probability density function (PDF) of the mixture fraction is crucial to the prediction of combustion since mixing plays an important role in turbulent non-premixed and partially premixed flames. This work provides an assessment of the large-eddy simulation (LES)/filtered density function (FDF) method for the prediction of the PDF of the mixture fraction. The advantage of the LES/FDF method is that it provides the full predictions of the statistical distribution of scalars including the mixture fraction. The predictive accuracy of the method for the PDF is yet to be fully validated. Assessing the prediction of the PDF of the mixture fraction, a conserved scalar, is an important starting point. The Sydney/Sandia inhomogeneous inlet jet flame is used as a test case. A quick comparison shows that the LES/FDF predicted PDF shapes of the mixture fraction deviate significantly from the commonly presumed Beta-PDF as well as from the experimental data in the flame. Here, to examine the source of the discrepancy, we clarify the different PDF definitions used in the comparison among the predictions, measurements, and the presumed shape PDFs. The discrepancy observed from the comparison is largely reconciled by clarifying the difference between the PDFs that are examined. The PDF of the resolved mixture fraction is shown to be close to the Beta-PDF in both the measurements and predictions, while the PDF directly deduced from the LES/FDF particles deviates significantly from the Beta-PDF. A multimodal PDF analysis and a pseudo convergence analysis are conducted to provide plausible evidence to support the predicted multimodal PDF shapes. The sub-filter scale FDF is shown to be close to the Beta-PDF too through the construction of a synthesized PDF, which supports the common presumed Beta-PDF assumption used in the presumed PDF methods when combined with LES.

42 ENGINEERING↗
DOE OSTI2024

The development of Gibbs's dyadic and implications for the gradient of a vector field

In this paper, we review the history of the dyadic as developed by Gibbs. This mathematical construct appeared in the second part of Gibbs's pamphlet on vector analysis (published in 1884), and it represented the first known development of a Cartesian theory of tensors. Gibbs made extensive use of the dyadic to express his theory of linear vector functions, that is, functions that acted on vectors and mapped them to new vectors. The dyadic proved to be a capable vehicle in Gibbs's hands, and his theory for dyadics (which we would now call second-order Cartesian tensors) was relatively advanced. The theory detailed notions such as the decomposition of vectors and conditions under which a tensor would have an inverse. While Gibbs's theory for linear operators expressed by dyadics was robust, it did not seem to garner the attention that the more conventional vector analysis (published in the first half of his pamphlet in 1881) did. Perhaps in part because of the general unfamiliarity with the dyadic, two distinct and conflicting definitions of the gradient of a vector field have arisen in the literature. The details of these differences in notation, possible reasons for the difference, and a potential resolution are proposed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A unified theory of zonal flow shears and density corrugations in drift wave turbulence

A unified theory of zonal flow shears and density corrugations in drift wave turbulence is presented. Polarization and density advection beat excitation are studied in combination with modulational response. Noise is driven by two-time flux correlation. While the effective zonal flow eddy viscosity can go negative, the zonal diffusivity is positive definite. There is no inverse cascade of density corrugation. The connection between avalanches and corrugations is discussed. Here, the zonal cross-correlation is identified and calculated. Conditions for alignment of zonal shears and corrugation gradients are determined, and the implications for staircase structure are discussed. We show that the synergy of beat noise and modulational effects is stronger than either alone. Strong zonal flows can be excited well below the modulational instability threshold. In the context of L–H transition, zonal noise quenches turbulence overshoot by eliminating the threshold for zonal flow excitation. The power threshold for L–H transition is lowered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Observation and characterisation of trapped electron modes in Wendelstein 7-X

In the past, quasi coherent (QC) modes were reported for nearly all tokamaks. The general definition describes modes as QC when the magnitude squared coherence is in the range of 0.3 to 0.6. QC modes are observed in the plasma core as well as in the plasma edge and can have quite different physical origins. The one in the core are observed in plasmas with low collisionality, where the electron temperature exceeds the ion temperature in the plasma core. This is the case for electron cyclotron heating in general. The origin of these modes are electrons trapped within a magnetic mirror, as reported in the past from various fusion devices. The so-called trapped-electron modes (TEMs) belong to drift wave instabilities and can be destabilized by electron-temperature gradients in the plasma core. From the diagnostic point of view, QC modes appear as fluctuations in electron density and temperature. Therefore, the microwave reflectometer is very well suited to monitor these modes. This paper describes experiments, conducted at the Wendelstein 7-X stellarator (W7-X), which aim at detecting QC modes at low wave numbers. A poloidal correlation reflectometer installed at W7-X, is able to measure low wave numbers ($k_\perp\unicode{x2A7D}{3.5}\,\textrm{cm}^{-1}$). For medium line-averaged densities ($\int n_\textrm{e}\unicode{x2A7D}{6 \times10^{19}}\,\textrm{m}^{-2}$) the plasma core is accessible for this diagnostic. For different magnetic configurations and plasma parameters, broad QC structures are observed in the coherence spectra. From the analysis of the rotation and the poloidal structure, these QC modes show the properties of electron-temperature-gradient driven TEMs. A linear relation between the mode velocity and the rotation frequency is found. The relation is uniform and confirms the nature of QC-mode observation as TEM in tokamaks, too.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

TEM turbulence in simulation and experiment in quasisymmetric and degraded-symmetry configurations of the HSX stellarator

The Helically Symmetric Experiment (HSX) has demonstrated reduced neoclassical transport and flow damping with quasisymmetry, but the difference of turbulent transport between quasisymmetric and degraded-symmetry configurations has not been definitively assessed. Profile analysis presented here shows that anomalous transport at the mid-radius is larger with degraded quasisymmetry when temperature and density profiles are matched. An initial comparison (${T_\mathrm{i}/T_\mathrm{e}} = 1$) between the quasisymmetric and degraded symmetry configurations finds that trapped-electron mode (TEM) turbulence can support larger heat fluxes in the degraded symmetry configuration despite smaller linear growth rates. However, further simulations at the experimental temperature ratio ($T_\mathrm{i}/T_\mathrm{e} = 0.2$) do not reproduce this effect. Experimental measurements of heat transport and density fluctuation amplitude are compared to nonlinear gyrokinetic simulations at the experimental temperature ratio. The heat flux and density fluctuation amplitude in simulation reproduce a substantial dependence on the density gradient, and the simulated heat flux matches measurements within experimental uncertainties, indicating that TEM turbulence is responsible for the anomalous transport at the steep density gradient in HSX.

HSX↗
DOE OSTI2024

Protocols and methodologies for acquiring and analyzing critical-current versus longitudinal-strain data in Bi 2 Sr 2 CaCu 2 O 8+x wires

Abstract In the literature on Bi 2 Sr 2 CaCu 2 O 8+ x (Bi-2212) superconducting wires, it is evident that measurement protocols for transport critical-current I c versus longitudinal strain ϵ and definitions of the so-called ‘strain limit’ are generally dissimilar. Yet, values obtained for the ‘strain limit’ are frequently assimilated to being those of the irreversible strain limit ϵ irr , regardless of the I c degradation-criterion used to define it. In effect, ϵ irr should correspond specifically to the I c ( ϵ ) irreversibility onset , where crack formation in Bi-2212 filaments presumably starts. Because I c ( ϵ ) degradation remains progressive over a fairly wide strain range beyond ϵ irr , the different I c degradation-criteria in use do not yield to the same result and, thus, are not equivalent from metrology perspective. Indeed, in studying densified samples of a modern Bi-2212 round wire, we found ϵ irr ≈ 0.4% and ϵ 5% ≈ 0.6% ( ϵ 5% being the strain where I c degrades by 5%). In this paper, we outline and suggest I c ( ϵ )-measurement protocols and data-analysis methodologies in the hope to converge the various approaches taken for studying Bi-2212 strain properties and, thus, remove related result discrepancies. A unified approach would enable more objective data comparisons among laboratories and among different Bi-2212 conductors. It would pave the way for more rigorous studies of effects potentially associated with wire design, powder, heat treatments, and other such parameters on the conductor’s strain properties.

protocols↗
DOE OSTI2025

Generalized fractional approach to solving partial differential equations with arbitrary dispersion relations

The dynamics of physical systems are typically modeled by partial differential equations (PDEs). Because of the complexity in physical systems, PDE models rely on approximations or limiting cases. To model the full complexity of physical systems, it is necessary to use multiscale approaches where appropriate models are used at each level. Alternatively, complexity can be directly tackled by fractional integrodifferential equations that must be derived for each system. Here, we argue that complexity can be restored in PDEs by describing them from a fractional calculus perspective. Rather than deriving a fractional integrodifferential equation, we reinterpret the dispersion relation of the system by use of the Riesz definition, which contains the required information relating the energy and momentum space of the system and thus fully describes their dynamics. The approach is demonstrated by two examples: the Landau–Lifshitz equation in a 1D ferromagnetic chain and a modified KdV equation supporting surface gravity waves or Euler dispersion. The presented approach is applicable to fluids, soft matter, and solid-state matter and can be readily generalized to higher dimensions and more complex systems. While numerical calculations are needed to determine the fractional operator, the approach is analytical and can be utilized to determine analytical solutions and investigate nonlinear problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Weak lensing trispectrum and Kurt-spectra

We introduce two kurt-spectra to probe fourth-order statistics of weak lensing convergence maps. Using state-of-the-art numerical simulations, we study the shapes of these kurt-spectra as a function of source redshifts and smoothing angular scales. We employ a pseudo-Cℓ approach to estimate the spectra from realistic convergence maps in the presence of an observational mask and noise for stage-IV large-scale structure surveys. We compare these results against theoretical predictions calculated using the FFTLog formalism, and find that a simple nonlinear clustering model — the hierarchical ansatz — can reproduce the numerical trends for the kurt-spectra in the nonlinear regime. In addition, we provide estimators for beyond fourth-order spectra where no definitive analytical results are available, and present corresponding results from numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Robustness of cosmic birefringence measurement against Galactic foreground emission and instrumental systematics

The polarization of the cosmic microwave background (CMB) can be used to search for parity-violating processes like that predicted by a Chern-Simons coupling to a light pseudoscalar field. Such an interaction rotates E modes into E modes in the observed CMB signal through an effect known as cosmic birefringence. Even though isotropic birefringence can be confused with the rotation produced by a miscalibration of the detectors' polarization angles, the degeneracy between both effects is broken when Galactic foreground emission is used as a calibrator. Here, in this work, we use realistic simulations of the High-Frequency Instrument of the Planck mission to test the impact that Galactic foreground emission and instrumental systematics have on the recent birefringence measurements obtained through this technique. Our results demonstrate the robustness of the methodology against the miscalibration of polarization angles and other systematic effects, like intensity-to-polarization leakage, beam leakage, or cross-polarization effects. However, our estimator is sensitive to the EB correlation of polarized foreground emission. Here we propose to correct the bias induced by dust EB by modeling the foreground signal with templates produced in Bayesian component-separation analyses that fit parametric models to CMB data. Acknowledging the limitations of currently available dust templates like that of the Commander sky model, high-precision CMB data and a characterization of dust beyond the modified blackbody paradigm are needed to obtain a definitive measurement of cosmic birefringence in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Soft pattern of gravitational Rutherford scattering from heavy target mass expansion

We investigate the soft behavior of the tree-level Rutherford scattering processes mediated via $t$-channel one-graviton exchange. We consider two types of Rutherford scattering processes, {\it e.g.}, a low-energy massless structureless projectile (up to spin-$1$) hits a static massive composite particle carrying various spins (up to spin-$2$), and a slowly-moving light projectile hits a heavy static composite target. The unpolarized cross sections in the first type are found to exhibit universal forms at the first two orders in $1/M$ expansion, yet differ at the next-to-next-to-leading order, though some terms at this order still remain universal or depend on the target spin in a definite manner. The unpolarized cross sections in the second type are universal at the lowest order in projectile velocity expansion and through all orders in $1/M$, independent of the spins of both projectile and target. The universality partially breaks down at relative order-$v^2/M^2$, albeit some terms at this order still depend on the target spin in a specific manner.

Rutherford scattering↗
DOE OSTI2022

Comparison between MAST-U conventional and Super-X configurations through SOLPS-ITER modelling

Abstract MAST-U has recently started operating with a Super-X divertor, designed to increase total flux expansion and neutral trapping, both predicted through simple analytic models and SOLPS calculations to reduce the plasma and impurity density detachment thresholds. In this study, utilising the SOLPS-ITER code, we are quantifying the possible gain allowed by the MAST-U Super-X and neutral baffling geometry, in terms of access to detachment. We show that a significant reduction of the upstream density detachment threshold (up to a factor 1.6) could be achieved in MAST-U, for the Super-X, as opposed to conventional divertor geometry, mainly through an increased total flux expansion, neutral trapping being found very similar between the different configurations. We also show that variations of the strike-point angle are complex to interpret in such a tightly baffled geometry, and that a case in which the target normal points more towards the separatrix does not necessarily imply a lower detachment threshold. As in previous calculations for TCV, we quantify the role of neutral effects through developing and applying a quantitative definition of neutral trapping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗
DOE OSTI2026

Time-dependent scenario modeling for the ST-E1 fusion power plant

ST-E1 is a low aspect ratio fusion power plant being designed by Tokamak Energy targeting 1.5 GW of fusion power. Characterization of the ST-E1 flat-top scenario is described elsewhere McNamara et al (2026 Nucl. Fusion 66 086008); here we focus on addressing the question of how to ramp-up the ST-E1 plasma from an initial state following breakdown and flux-surface formation to the target flat-top state. Being low-aspect ratio, the available solenoid flux of ST-E1 is limited. Therefore, particular consideration is placed on developing ramp-up scenarios that predominantly use inductive flux provided by external vertical field coils. Through time-dependent modeling with METIS, we show that this is possible when the ramp-up is performed at relatively high plasma density: although auxiliary current drive efficiency is reduced, this is significantly outweighed by (1) higher electron-ion collisional equilibration, (2) higher fusion power once ions become sufficiently hot, (3) higher poloidal beta for increased vertical-field flux, and (4) potentially favorable exhaust compatibilities. Ultimately, we show the target ST-E1 flat-top performance can be reached after a ramp-up period lasting 150 s using less than 40 Vs of solenoid flux (with vertical field providing ∼ 90 Vs of flux). The sensitivity to model assumptions are presented, with the general observation that deleterious effects can be mitigated through minor alterations of the auxiliary power temporal waveform and/or total auxiliary power level. The impact of a solenoid and the auxiliary power mix (electron cyclotron heating only versus electron and ion cyclotron heating) on the ST-E1 ramp-up success are also discussed in appendices. On this latter topic, we show that the effect of direct-ion heating during ramp-up is obscured by the uncertainty in the pedestal dynamics, identifying a clear line of future work required to make a definite decision on the ST-E1 auxiliary power mix.

current ramp-up↗