Search NASA⌕ Search

SEARCH · Search NASA

Results for “triangularity”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Face- and Cell-Averaged Nodal-Gradient Approach to Cell-Centered Finite-Volume Method on Mixed Grids

In this paper, the averaged nodal-gradient approach previously developed for triangular grids is extended to mixed triangular-quadrilateral grids. It is shown that the face- averaged approach leads to deteriorated iterative convergence on quadrilateral grids. To develop a convergent solver, we consider cell-averaging instead of face-averaging for quadri- lateral cells. We show that the cell-averaged approach leads to a convergent solver and can be efficiently combined with the face-averaged approach on mixed grids. The method is demonstrated for various inviscid and viscous problems from low to high Mach numbers on two-dimensional mixed grids.

Nishikawa, Hiroaki↗
NASA NTRSDate not supplied

A Simplified FANG Cell-Centered Finite-Volume Method and Comparison with Other Methods for Trouble-Prone Grids

We propose a simplication of the face-averaged nodal-gradient (FANG) method for a cell-centered finite-volume Euler/Navier-Stokes solver on arbitrary grids, and compare it with other gradient methods for trouble-prone grids in two dimensions. The implementation of the FANG method is simplified by adding the face-neighbor cells of the cells around a node to a least-squares gradient stencil. The resulting method is stable for both triangular and quadrilateral grids. Although it increases the residual stencil for triangular grids, it allows the solver to work seamlessly for mixed grids and greatly simplifies the implementation, especially in three dimensions. For comparison, only explicit weighted/unweighted least-squares cell-centered and nodal gradient methods are considered. These gradients are used in both inviscid and viscous schemes, and we investigate their impact on the iterative convergence of an implicit defect-correction solver on difficult grids such as highly-curved-and-thin grids and highly distorted anisotropic grids. Finally, we will also consider a face-stencil-based limiter and compare it with a conventional cell-stencil-based limiter.

Hiroaki Nishikawa↗
DOE OSTI2022

Neumann Series in MGS-GMRES and Inner-Outer Iterations: Preprint

A low-synchronization MGS-GMRES Krylov solver employing a truncated Neumann series for the inverse compact WY MGS correction matrix T is presented. A corollary to the backward stability result of Paige et al. [1] establishes that T = I - Lk is sufficient for convergence of GMRES when kLkp F = O("p)_p F (B), where the strictly lower triangular matrix L is defined by the inner products of Krylov vectors V T 1:k-2 vk-1. The preconditioner is the classical Ruge-Stuben AMG algorithm with compatible relaxation and inner-outer Gauss-Seidel smoother. This smoother may also be expressed as a truncated Neumann series. Drop tolerances are applied to the lower triangular matrices arising in the smoother in order to reduce the number of non-zeros and accelerate the time to solution. The number of small matrix elements are found to increase from fine to coarse levels and thus the effciency gains are greater for large problems with many levels in the V -cycle. The solver is applied to the pressure continuity equation for the incompressible Navier-Stokes equations. Unlike the inner-outer iteration, the solver convergence rate with the standard Gauss-Seidel smoother deteriorates with dropping. The solver compute time is reduced by up to 50% without a change in the convergence rate.

Gauss-Seidel smoother↗
DOE OSTI2022

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗
DOE OSTI2023

Experiments to Measure the Effect of Tantalum on Critical Systems

Sandia National Laboratories (SNL) and Oak Ridge National Laboratory (ORNL) collaborated to develop a capability to test the epithermal/intermediate cross sections of materials at the SNL critical experiment facility using the Seven Percent Critical Experiment (7uPCX) fuel. As a result, a new set of critical experiments has been designed to target the epithermal cross sections of tantalum (Ta) and is scheduled to be performed at SNL in 2023. These critical experiments will be evaluated for inclusion in the 2024 edition of the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook. The focus of these critical experiments is to measure the effects of Ta on the critical array size. The critical array size will be determined by an approach-to-critical experiment with the number of fuel rods in the array as a free parameter. The core configurations are designed to optimize the reactivity worth of Ta and the overall percentage of Ta absorption rates in the epithermal/intermediate energy range (0.625 eV – 100 keV). The baseline core configuration includes 7uPCX fuel rods set at a triangular pitch of 1.016 cm and a central dry test region that utilizes a cadmium liner for filtering out thermal neutrons. The central test region has locations for 85 Ta rods set at a triangular pitch of 0.813 cm. The Ta reactivity worth for the case with 85 Ta rods is approximately 2.55%, with the percentage of Ta absorption rates within the intermediate energy range at nearly 90%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Hole‐Carrier‐Dominant Transport in 2D Single‐Crystal Copper

Abstract In 2D noble metals like copper, the carrier scattering at grain boundaries has obscured the intrinsic nature of electronic transport. However, it is demonstrated that the intrinsic nature of transport by hole carriers in 2D copper can be revealed by growing thin films without grain boundaries. As even a slight deviation from the twin boundary is perceived as grain boundaries by electrons, it is only through the thorough elimination of grain boundaries that the hidden hole‐like attribute of 2D single‐crystal copper can be unmasked. Two types of Fermi surfaces, a large hexagonal Fermi surface centered at the zone center and the triangular Fermi surface around the zone corner, tightly matching to the calculated Fermi surface topology, confirmed by angle‐resolved photoemission spectroscopy (ARPES) measurements and vivid nonlinear Hall effects of the 2D single‐crystal copper account for the presence of hole carriers experimentally. This breakthrough suggests the potential to manipulate the majority carrier polarity in metals by means of grain boundary engineering in a 2D geometry.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗
DOE OSTI2018

Contribution to Spin Tubes Study

In this work, we aim to study both analytically and numerically the behavior of triangular spin tubes made by three chains which are connected to one another by exchange, in the presence of single-ion anisotropy. By diagonalizing the spin Hamiltonian, the partition function was exactly calculated within the transfer matrix approach and several relevant physical quantities were derived and analyzed. Depending on the sign and strengths of the exchange and anisotropy constants, the system reveals a rich variety of features such as frustration and steplikes in the magnetization curves arising from both geometry and conflicting interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Increasing Biomass Production on Limited Land Area Through an Optimal Planting Arrangement

Planting density is a primary consideration in silviculture; however, planting arrangement is often ignored. Most, if not all, forest plantations are arranged in rectangular or square lattices (i.e., grids). Using a simple mathematical model, we investigate the potential influence of planting arrangement on planting density, biomass yield, and rotation period by assuming that efficiently arranging trees is similar to packing congruent circles on a plane. The hexagonal lattice achieves the densest circle packing on a plane; therefore, a hexagonal or triangular lattice arrangement of stems provides the highest planting density for a given spacing. Using packing density to quantify arrangement efficiency, tree crowns in a hexagonal lattice fill approximately 90.7% of available area at initial canopy contact, while tree crowns in a square lattice fill approximately 78.5% of available area at initial canopy contact. The hexagonal lattice permits about a 15% higher density than the square lattice, which allows canopy closure to occur earlier without any change in individual tree growth. Short rotation woody crop (SRWC) systems are excellent candidates under the model's assumptions of level stand with even-age monoculture. If belowground resources are non-limiting, a hexagonal lattice arrangement shortens rotation period and thus optimizes the biomass yield per land area over time. Lastly, higher productivity over time is central to sustainable and efficient use of limited area for bioenergy and biomass products.

09 BIOMASS FUELS↗
DOE OSTI2018

A spin glass state in Ba 3 TiRu 2 O 9

The magnetic properties of Ba 3 TiRu 2 O 9 , whose crystal structure is based on stacked triangular planar lattices of M 2 O 9 dimers and single MO 6 octahedra, are reported. The system is magnetically disturbed by a substantial amount of Ti/Ru chemical disorder. The Weiss temperature and effective magnetic moment were found to be –29.5 K and 1.82 μ B /f.u. respectively, and a bifurcation in the zero field cooled and field cooled magnetic susceptibility is observed below 4.7 K, suggesting that this is a compositionally-disordered spin-glass system. Finally, the material is a semiconductor with an activation energy for charge transport of approximately 0.14 eV.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Static and dynamic magnetic properties of the spiin-$\frac{5}{2}$triangle lattice antiferromagnet Na 3 Fe(PO 4 ) 2 studied by 31 P NMR

We report 31 P nuclear magnetic resonance (NMR) measurements have been carried out to investigate the magnetic properties and spin dynamics of Fe 3+ (S = 5/2) spins in the two-dimensional triangular lattice (TL) compound Na 3 Fe(PO 4 ) 2 . The temperature (T) dependence of nuclear spin-lattice relaxation rates (1/T 1 ) shows a clear peak around Néel temperature, T N = 10.9 K, corresponding to an antiferromagnetic (AFM) transition. From the temperature dependence of NMR shift (K) above T N , an exchange coupling between Fe 3+ spins was estimated to be J/k B ≃ 1.9 K using the spin-5/2 Heisenberg isotropic-TL model. The temperature dependence of 1/T 1 T divided by the magnetic susceptibility (χ), 1/T 1 T χ , above T N proves the AFM nature of spin fluctuations below ~50 K in the paramagnetic state. In the magnetically ordered state below T N , the characteristic rectangular shape of the NMR spectra is observed, indicative of a commensurate AFM state in its ground state. The strong temperature dependence of 1/T 1 in the AFM state is well explained by the two-magnon (Raman) process of the spin waves in a 3D antiferromagnet with a spin-anisotropy energy gap of 5.7 K. The temperature dependence of sublattice magnetization is also well reproduced by the spin waves. Those results indicate that the magnetically ordered state of Na 3 Fe(PO 4 ) 2 is a conventional 3D AFM state, and no obvious spin frustration effects were detected in its ground state

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

A data-driven approach to real-time vertical position estimation for NSTX-U vertical stability control

In this paper, a database of 77 996 plasma equilibrium reconstructions from 727 discharges during the initial operation of the NSTX-U spherical tokamak is analyzed to develop a statistically robust model of the plasma vertical position for real-time control. A variety of regression models are developed and tested, ranging in complexity from linear models to deep neural networks, and including input signals ranging from the four pairs of flux loops used historically on NSTX-U up to the full set of 389 real-time signals available to the plasma control system. A linear model based on 140 real-time magnetics signals is found to offer excellent accuracy, with a coefficient of determination R 2 = 0.906. The robustness of this model to limited training data, new operating scenarios, and signal errors is tested, and a procedure is demonstrated to tune the model parameters to optimize its robustness. A time-dependent plasma equilibrium solver, TokaMaker, is used to simulate vertical stability control in NSTX-U, demonstrating that it should be possible to iteratively tune the parameters of a linear vertical position model to stabilize both positive and negative triangularity plasmas in future experiments.

magnetic diagnostics↗
DOE OSTI2022

Radiative pulsed L-mode operation in ARC-class reactors

Abstract A new ARC-class, highly-radiative, pulsed, L-mode, burning plasma scenario is developed and evaluated as a candidate for future tokamak reactors. Pulsed inductive operation alleviates the stringent current drive requirements of steady-state reactors, and operation in L-mode affords ELM-free access to ∼ 90 % core radiation fractions, significantly reducing the divertor power handling requirements. In this configuration the fusion power density can be maximized despite L-mode confinement by utilizing high-field to increase plasma densities and current. This allows us to obtain high gain in robust scenarios in compact devices with P fus > 1000 MW despite low confinement. We demonstrate the feasibility of such scenarios here; first by showing that they avoid violating 0D tokamak limits, and then by performing self-consistent integrated simulations of flattop operation including neoclassical and turbulent transport, magnetic equilibrium, and radiofrequency current drive models. Finally we examine the potential effect of introducing negative triangularity with a 0D model. Our results show high-field radiative pulsed L-mode scenarios are a promising alternative to the typical steady state advanced tokamak scenarios which have dominated tokamak reactor development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MHD stability analysis against pressure and current-driven modes in the SMall Aspect Ratio Tokamak

Linear magnetohydrodynamic (MHD) simulations for the SMall Aspect Ratio Tokamak (SMART) have been carried out for the first time, for both positive (PT) and negative triangularity (NT) shaped plasmas using the MARS-F code. The MHD stability of projected SMART plasmas against internal kinks, infernal modes and edge peeling-ballooning modes have been analyzed for a wide range of realistic equilibria. A stabilization of internal kinks and infernal modes is observed when increasing the safety factor profile and reducing plasma beta. PT shaped plasmas are more stable against both internal kinks and infernal modes than their counterpart NT shaped plasmas. Toroidal flows have little impact on the MHD stability of the internal kinks, but they have a strong stabilizing effect on infernal modes, which can be further mitigated in NT shaped plasmas. The MHD stability of peeling-ballooning modes is reduced in NT shaped plasmas, as observed in conventional tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Novel intrinsic helical cores and MHD dynamo flux pumping evidence in DIII-D

Evidence of a MHD dynamo brought about by intrinsic helical core structures in negative triangularity (NT) and ITER baseline scenario (IBS) plasmas demonstrates paths to operating tokamaks without seeds for deleterious global MHD modes. A stationary sawtooth- and edge localized mode-free state maintains high performance for long durations in a strong NT shape in the DIII-D tokamak. Disparity between minimum safety factor values predicted by magnetic flux diffusion theory and derived from experiment suggests anomalous flux diffusion by an MHD dynamo, driven by flows associated with quasi-interchange modes. Stability calculations with the GATO code find NT and IBS experimental equilibria unstable to quasi-interchange modes, in agreement with observed saturated activity. IBS discharges with similar magnetic winding structure occupy robust helical core states following a transition from sawtooth instabilities. These states agree with prior modeling of helical core bifurcation thresholds with the VMEC equilibrium code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Magnetism and quantum melting in moiré-material Wigner crystals

Recent experiments have established that semiconductor-based moir´e materials can host incompressible states at a series of fractional moir´e-miniband fillings. These states have been identified as generalized Wigner crystals in which electrons localize on a subset of the available triangular-lattice moir´e superlattice sites. In this article, we use momentum-space exact diagonalization to investigate the many-body ground state evolution at rational fillings from the weak-hopping classical lattice gas limit, in which only spin degrees-of-freedom are active at low energies, to the strong-hopping metallic regime where the Wigner crystals melt. We specifically address the nature of the magnetic ground states of the generalized Wigner crystals at fillings ν = 1/3 and ν = 2/3.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effects of pressure on the electronic and magnetic properties of bulk NiI 2

Transition metal dihalides have recently garnered interest in the context of two-dimensional van der Waals magnets as their underlying geometrically frustrated triangular lattice leads to interesting competing exchange interactions. In particular, NiI 2 is a magnetic semiconductor that has been long known for its exotic helimagnetism in the bulk. Recent experiments have shown that the helimagnetic state survives down to the monolayer limit with a layer-dependent magnetic transition temperature that suggests a relevant role of the interlayer coupling. Here, we explore the effects of hydrostatic pressure as a means to enhance this interlayer exchange and ultimately tune the electronic and magnetic response of NiI 2 . We study first the evolution of the structural parameters as a function of external pressure using first-principles calculations combined with x-ray diffraction measurements. We then examine the evolution of the electronic structure and magnetic exchange interactions via first-principles calculations and Monte Carlo simulations. Further, we find that the leading interlayer coupling is an antiferromagnetic second-nearest-neighbor interaction that increases monotonically with pressure. The ratio between isotropic third- and first-nearest-neighbor intralayer exchanges, which controls the magnetic frustration and determines the magnetic propagation vector q of the helimagnetic ground state, is also enhanced by pressure. As a consequence, our Monte Carlo simulations show a monotonic increase in the magnetic transition temperature, indicating that pressure is an effective means to tune the magnetic response of NiI 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Wavelet-based local mesh refinement for rainfall–runoff simulations

A wavelet-based local mesh refinement (wLMR) strategy is designed to generate multiresolution and unstructured triangular meshes from real digital elevation model (DEM) data for efficient hydrological simulations at the catchment scale. The wLMR strategy is studied considering slope- and curvature-based refinement criteria to analyze DEM inputs: the slope-based criterion uses bed elevation data as input to the wLMR strategy, whereas the curvature-based criterion feeds the bed slope data into it. The performance of the wLMR meshes generated by these two criteria is compared for hydrological simulations; first, using three analytical tests with the systematic variation in topography types and then by reproducing laboratory- and real-scale case studies. The bed elevation on the wLMR meshes and their simulation results are compared relative to those achieved on the finest uniform mesh. Repeated analytical tests show that the slope- and curvature-based criteria are equally effective with the wLMR strategy, and that it is easier to decide which criterion to take in relation to the (regular) shape of the topography. For the realistic case studies: (i) slope analysis provides a better metric to assess the correlation of a wLMR mesh to the fine uniform mesh and (ii) both criteria predict outlet hydrographs with a close predictive accuracy to that on the uniform mesh, but the curvature-based criterion is found to slightly better capture the channeling patterns of real DEM data.

54 ENVIRONMENTAL SCIENCES↗