Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,243 records · Page 69

Characterizing Tradeoffs in Memory, Accuracy, and Speed for Chemistry Tabulation Techniques

Chemistry tabulation is a common approach in practical simulations of turbulent combustion at engineering scales. Linear interpolants have traditionally been used for accessing precomputed multidimensional tables but suffer from large memory requirements and discontinuous derivatives. Higher-degree interpolants address some of these restrictions but are similarly limited to relatively low-dimensional tabulation. Artificial neural networks (ANNs) can be used to overcome these limitations but cannot guarantee the same accuracy as interpolants and introduce challenges in reproducibility and reliable training. These challenges are enhanced as the physics complexity to be represented within the tabulation increases. Here, we assess the efficiency, accuracy, and memory requirements of Lagrange polynomials, tensor product B-splines, and ANNs as tabulation strategies. We analyze results in the context of nonadiabatic flamelet modeling where higher dimension counts are necessary. While ANNs do not require structuring of data, providing benefits for complex physics representation, interpolation approaches often rely on some structuring of the table. Interpolation using structured table inputs that are not directly related to the variables transported in a simulation can incur additional query costs. This is demonstrated in the present implementation of heat losses. We show that ANNs, despite being difficult to train and reproduce, can be advantageous for high-dimensional, unstructured datasets relevant to nonadiabatic flamelet models. Furthermore we demonstrate that Lagrange polynomials show significant speedup for similar accuracy compared to B-splines.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MOOSE Reactor Module: An Open-Source Capability for Meshing Nuclear Reactor Geometries

The U.S. Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has developed numerous physics solvers utilizing the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) framework for multiphysics reactor analysis. These solvers require input finite element meshes representing the discretized spatial domain. Typically, reactor analysts turn to licensed tools for the creation of reactor geometry meshes. Recently, open-source functionality has been added to the MOOSE framework to mesh common reactor geometries and improve MOOSE-based nuclear reactor application user workflows. The new functionality is primarily contained in the new Reactor module of MOOSE and includes support for hexagonal pins, assemblies, and cores, extended Cartesian geometry support, options for modeling static and rotating control drums within a hexagonal assembly, core periphery triangulation, and automatic tagging of pin, assembly, plane, and depletion regions for easier post processing of physics results. A set of reactor geometry mesh builder objects further streamlines the construction of hexagonal and Cartesian cores and allows mapping of materials to regions during mesh generation. The meshes produced with the MOOSE Reactor module may be used directly within MOOSE-based applications or exported as Exodus II files for use in other finite element solvers. The tools have been demonstrated and verified using a variety of NEAMS physics solvers on a range of reactor applications, including a sodium-cooled fast reactor core analysis using Griffin, a fast reactor assembly thermal deformation analysis using MOOSE Tensor Mechanics, and a heat pipe–cooled microreactor coupled analysis using Griffin, Bison, and Sockeye. MOOSE’s Reactor module provides significant advantages compared to the use of external meshing tools when analyzing Cartesian and hexagonal reactor lattices using MOOSE-based applications: immediate accessibility (open-source) to the end user, low barrier to entry for new users, speed of mesh generation, volume preservation of meshed fuel pins, and simplification of analysis workflow when used in conjunction with MOOSE-based applications.

99 GENERAL AND MISCELLANEOUS↗

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗

Discovering Active Subspaces for High-Dimensional Computer Models

Dimension reduction techniques have long been an important topic in statistics, and active subspaces (AS) have received much attention this past decade in the computer experiments literature. The most common approach towards estimating the AS is to use Monte Carlo with numerical gradient evaluation. While sensible in some settings, this approach has obvious drawbacks. Recent research has demonstrated that active subspace calculations can be obtained in closed form, conditional on a Gaussian process (GP) surrogate, which can be limiting in high-dimensional settings for computational reasons. In this paper, we produce the relevant calculations for a more general case when the model of interest is a linear combination of tensor products. These general equations can be applied to the GP, recovering previous results as a special case, or applied to the models constructed by other regression techniques including multivariate adaptive regression splines (MARS). Furthermore, using a MARS surrogate has many advantages including improved scaling, better estimation of active subspaces in high dimensions and the ability to handle a large number of prior distributions in closed form. In one real-world example, we obtain the active subspace of a radiation-transport code with 240 inputs and 9,372 model runs in under half an hour.

97 MATHEMATICS AND COMPUTING↗

Accelerating boundary analog of a Kerr black hole

An accelerated boundary correspondence (i.e. a flat spacetime accelerating mirror trajectory) is derived for the Kerr spacetime, with a general formula that ranges from the Schwarzschild limit (zero angular momentum) to the extreme maximal spin case (yielding asymptotic uniform acceleration). Here, the beta Bogoliubov coefficients reveal the particle spectrum is a Planck distribution at late times with temperature cooler than a Schwarzschild black hole, due to the 'spring constant' analog of angular momentum. The quantum stress tensor indicates a constant emission of energy flux at late times consistent with eternal thermal equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗

Semiclassical dynamics of Hawking radiation

Here, we consider gravity in 3+1 spacetime dimensions coupled to N scalar matter fields in a semiclassical limit where $N\to\infty$. The dynamical evolution of a black hole including the back-reaction of the Hawking radiation on the metric is formulated as an initial-value problem. The quantum stress-energy tensor is evaluated using a point-splitting regularization along spacelike geodesics. To account for the quantum entanglement of the matter fields, they are treated as a set of bilocal collective fields defined on spacelike hypersurfaces. The resulting semiclassical field equations include terms up to fourth order in derivatives that can be treated in a perturbative $\hbar$ expansion. The formulation we arrive at should be amenable to numerical simulation of time-dependent semiclassical spacetime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Path integral factorization and the gravitational effective action

Abstract We discuss the factorization and continuity properties of fields in the Euclidean gravitational path integral with higher dimension operators constructed from powers of the Riemann tensor. We construct the boundary terms corresponding to the microcanonical ensemble and show that the saddle point approximation to the path integral with a quasilocal energy constraint generally yields a saddle point with discontinuous temperature. This extends a previous result for the Euclidean Schwarzschild–de Sitter geometry in Einstein gravity and shows that it is robust against at least some types of quantum corrections from heavy fields. As an application, we compute the entropy of SdS in D = 4 using the BTZ method. Our result matches the entropy calculated using Wald’s formula.

Astronomy & Astrophysics↗

Inflation alternative via the gravitational field of a singularity

Abstract We explore the scenario that the observable Universe emerged from the vicinity of a negative mass ring singularity, and all content of the Universe travels at the same group velocity close to the speed of light on a geodesic trajectory along the axis of rotation of the singularity. In appropriate coordinate parametrization and evaluated on the trajectory, we find that the metric tensor in the vicinity of the trajectory exhibits a conformal scale factor a ( η ) with contraction and subsequent expansion properties that solve the horizon problem. We then introduce a static flow of gravitating radiation along the trajectory (perturbatively with respect to the mass scale of the singularity) to model a homogeneous radiation dominated Universe. Solving the Einstein field equations with a physically motivated ansatz of metric perturbation then reveals that the effective conformal scale factor indeed grows asymptotically with the same power law as expected in a conventional radiation dominated Universe.

Astronomy & Astrophysics↗

Inertial and anisotropic pressure effects on cross-field electron transport in low-temperature magnetized plasmas

In this paper, a one-dimensional (1D) particle-in-cell Monte Carlo collision (PIC-MCC) model is developed to investigate the effects of anisotropic pressure and inertial terms due to non-Maxwellian velocity distribution functions on cross-field electron transport. The conservation of momentum is evaluated by taking the moments of the first-principles gas-kinetic equation. A steady-state discharge is obtained without any low-frequency ionization oscillations by considering an anomalous electron scattering profile. Further, the results obtained from the 1D PIC-MCC model are compared with fluid models, including the quasi-neutral drift-diffusion (DD), non-neutral DD, and full fluid moment models. The discharge current obtained from the PIC-MCC model is in good agreement with the fluid models. The cross-field electron transport due to the inertial terms, i.e. the gradient of axial and azimuthal drift, is evaluated. Moreover, PIC-MCC simulation results show non-zero, anisotropic, off-diagonal pressure tensor terms due to asymmetric non-Maxwellian electron velocity distribution function, potentially contributing to cross-field electron transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Changes of deformed shell gaps at N ~ 100 in light rare-earth, neutron-rich nuclei

There have been compelling evidences indicating that deformed single-particle states in neutron- rich regions are different from those in the stable region. It was pointed out that for the light rare- earth (60Nd, 62Sm, and 64Gd), neutron-rich (N = 98-102) nuclei, the Woods-Saxon potential, the Nilsson modified oscillator potential with "universal" parameters, and the folded Yukawa potential all failed to describe the correct ordering of neutron single-particle states. The location and size of deformed shell gaps in this mass region are under current debate. We propose a modification for the "standard" Nilsson parameters of Bengtsson and Ragnarsson introduced in 1985. The proposed N-dependent spin-orbit interaction brings an effect of the tensor force for the deformed potential, which causes dynamic changes in deformed shell gaps with neutron number and deformation. By applying the modified Nilsson parameters to generate deformed bases for the Projected Shell Model, we demonstrate that the calculation can explain consistently the current experimental data, including the ground state configuration in odd-neutron nuclei, upbending of the yrast moment of inertia at higher spins and the energies of 2-quasineutron 6- and 4- isomers in even-even nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetic properties of the Ising-like rare earth pyrosilicate: D-Er 2 Si 2 O 7

Ising-like spin-1/2 magnetic materials are of interest for their ready connection to theory, particularly in the context of quantum critical behavior. In this work we report detailed studies of the magnetic properties of a member of the rare earth pyrosilicate family, D-Er 2 Si 2 O 7 , which is known to display a highly anisotropic Ising-like g-tensor and effective spin-1/2 magnetic moments. We used powder neutron diffraction, powder inelastic neutron spectroscopy (INS), and single crystal AC susceptibility to characterize its magnetic properties. Neutron diffraction enabled us to determine the magnetic structure below the known transition temperature (T N = 1.9 K) in zero field, confirming that the magnetic state is a four-sublattice antiferromagnetic structure with two non-collinear Ising axes, as was previously hypothesized. Our powder INS data revealed a gapped excitation at zero field, consistent with anisotropic (possibly Ising) exchange. An applied field of 1 T produces a mode softening, which is consistent with a field-induced second order phase transition. To assess the relevance of D-Er 2 Si 2 O 7 to the transverse field Ising model, we performed AC susceptibility measurements on a single crystal with the magnetic field oriented in the direction transverse to the Ising axes. This revealed a transition at 2.65 T at 0.1 K, a field significantly higher than the mode-softening field observed by powder INS, showing that the field-induced phase transitions are highly field-direction dependent as expected. Furthermore, these measurements suggest that D-Er 2 Si 2 O 7 may be a candidate for further exploration related to the transverse field Ising model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Low-frequency Raman active modes of twisted bilayer MoS 2

Abstract We study the low-frequency Raman active modes of twisted bilayer MoS 2 for several twist angles using a force-field approach and a parametrized bond polarizability model. We show that twist angles near high-symmetry stacking configurations exhibit stacking frustration that leads to significant buckling of the moiré superlattice. We find that atomic relaxation due to the twist is of prime importance. The periodic displacement of the Mo atoms shows the realization of a soliton network, and in turn, leads to the emergence of a number of frequency modes not seen in the high-symmetry stacking systems. Some of the modes are only seen in the XZ Raman polarization setup while others are seen in the XY setup. The symmetry of the normal modes, and how this affects the Raman tensors is examined in detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cryogenic platform to investigate strong microwave cavity-spin coupling in correlated magnetic materials

We present a comprehensive exploration of loop-gap resonators for electron spin resonance (ESR) studies, enabling investigations into the hybridization of solid-state magnetic materials with microwave polariton modes. The experimental setup, implemented in a Physical Property Measurement System by Quantum Design, allows for measurements of ESR spectra at temperatures as low as 2 Kelvin. The versatility of continuous wave ESR spectroscopy is demonstrated through experiments on CuSO 4 •5H 2 O and MgCr 2 O 4 , showcasing the g-tensor and magnetic susceptibilities of these materials. The study delves into the challenges of fitting spectra under strong hybridization conditions and underscores the significance of proper calibration and stabilization. The detailed guide provided serves as a valuable resource for laboratories interested in exploring hybrid quantum systems through microwave resonators.

electron spin resonance spectroscopy↗

An assessment of full-wave effects on Maxwellian lower-hybrid wave damping

Lower-hybrid current drive (LHCD) actuators are important components of modern day fusion experiments as well as proposed fusion reactors. However, simulations of LHCD often differ substantially from experimental results, and from each other, especially in the inferred power deposition profile shape. Here we investigate some possible causes of this discrepancy; ‘full-wave’ effects such as interference and diffraction, which are omitted from standard raytracing simulations and the breakdown of the raytracing near reflections and caustics. We compare raytracing simulations to state-of-the-art full-wave simulations using matched hot-plasma dielectric tensors in realistic tokamak scenarios for the first time. We show that differences between full-wave simulations and raytracing in previous work were primarily due to numerical and physical inconsistencies in the simulations, and we demonstrate that quantitative agreement between raytracing and converged full-wave simulations can be obtained in reactor relevant-scenarios and qualitative agreement can be obtained in situations with weak damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Closure theory for high-collisionality multi-ion plasmas

A general formalism is developed to construct and solve a system of linearized moment equations for parallel and perpendicular closures in high-collisionality plasmas. It is applicable for multiple ion species with arbitrary masses, temperatures, charges, and densities. The convergence of closure coefficients is evaluated by increasing the number of moments from 2 to 32 for scalar, vector, and rank-2 tensor moments. As an example, the complete set of closure coefficients for a deuterium-carbon plasma over the entire Hall parameter range is presented. Furthermore, the closure coefficients at various temperature ratios show that the one-temperature closure coefficients can differ significantly from the two-temperature coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The magnetic gradient scale length explains why certain plasmas require close external magnetic coils

Abstract The separation between the last closed flux surface of a plasma and the external coils that magnetically confine it is a limiting factor in the construction of fusion-capable plasma devices. This plasma-coil separation must be large enough so that components such as a breeding blanket and neutron shielding can fit between the plasma and the coils. Plasma-coil separation affects reactor size, engineering complexity, and particle loss due to field ripple. For some plasmas it can be difficult to produce the desired flux surface shaping with distant coils, and for other plasmas it is infeasible altogether. Here, we seek to understand the underlying physics that limits plasma-coil separation and explain why some configurations require close external coils. In this paper, we explore the hypothesis that the limiting plasma-coil separation is set by the shortest scale length of the magnetic field as expressed by the ∇ B tensor. We tested this hypothesis on a database of > 40 stellarator and tokamak configurations. Within this database, the coil-to-plasma distance compared to the minor radius varies by over an order of magnitude. The magnetic scale length is well correlated to the coil-to-plasma distance of actual coil designs generated using the REGCOIL method (Landreman 2017 Nucl. Fusion 57 046003). Additionally, this correlation reveals a general trend that larger plasma-coil separation is possible with a small number of field periods.

Physics↗

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗

Electron kinetics in a positive column of AC discharges in a dynamic regime

Abstract We have performed hybrid kinetic-fluid simulations of a positive column in alternating current (AC) argon discharges over a range of driving frequencies f and gas pressure p for the conditions when the spatial nonlocality of the electron energy distribution function (EEDF) is substantial. Our simulations confirmed that the most efficient conditions of plasma maintenance are observed in the dynamic regime when time modulations of mean electron energy (temperature) are substantial. The minimal values of the root mean square electric field and the electron temperature have been observed at f/p values of about 3 kHz Torr −1 in a tube of radius R = 1 cm. The ionization rate and plasma density reached maximal values under these conditions. The numerical solution of a kinetic equation allowed accounting for the kinetic effects associated with spatial and temporal nonlocality of the EEDF. Using the kinetic energy of electrons as an independent variable, we solved an anisotropic tensor diffusion equation in phase space. We clarified the role of different flux components during electron diffusion in phase space over surfaces of constant total energy. We have shown that the kinetic theory uncovers a more exciting and rich physics than the classical ambipolar diffusion (Schottky) model. Non-monotonic radial distributions of excitation rates, metastable densities, and plasma density have been observed in our simulations at pR > 6 Torr cm. The predicted off-axis plasma density peak in the dynamic regime has never been observed in experiments so far. We hope our results stimulate further experimental studies of the AC positive column. The kinetic analysis could help uncover new physics even for such a well-known plasma object as a positive column in noble gases.

Physics↗