Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “toroidal grids”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Verification of the ideal magnetohydrodynamic response at rational surfaces in the VMEC code

The VMEC nonlinear ideal MHD equilibrium code [S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553 (1983)] is compared against analytic linear ideal MHD theory in a screw-pinch-like configuration. The focus of such analysis is to verify the ideal MHD response at magnetic surfaces which possess magnetic transform (ι) which is resonant with spectral values of the perturbed boundary harmonics. A large aspect ratio circular cross section zero-beta equilibrium is considered. This equilibrium possess a rational surface with safety factor q = 2 at a normalized flux value of 0.5. A small resonant boundary perturbation is introduced, exciting a response at the resonant rational surface. The code is found to capture the plasma response as predicted by a newly developed analytic theory that ensures the existence of nested flux surfaces by allowing for a jump in rotational transform (ι=1/q). The VMEC code satisfactorily reproduces these theoretical results without the necessity of an explicit transform discontinuity (Δι) at the rational surface. It is found that the response across the rational surfaces depends upon both radial grid resolution and local shear (dι/dΦ, where ι is the rotational transform and Φ the enclosed toroidal flux). Calculations of an implicit Δι suggest that it does not arise due to numerical artifacts (attributed to radial finite differences in VMEC) or existence conditions for flux surfaces as predicted by linear theory (minimum values of Δι). Scans of the rotational transform profile indicate that for experimentally relevant levels of transform shear the response becomes increasing localised. Furthermore, careful examination of a large experimental tokamak equilibrium, with applied resonant fields, indicates that this shielding response is present, suggesting the phenomena is not limited to this verification exercise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Pressure flux-split viscous solutions for subsonic diffusers

A method for predicting flow in diffusers with inlet swirl has been developed. Solutions of the pressure flux-split Reduced Navier Stokes (RNS) equations are obtained for flows in axisymmetric diffusers. Viscous-Inviscid interactions and flow fields with toroidal recirculation regions are efficiently captured. The computational model is verified by comparision with experimental data and other computations. The extreme sensitivity to grid, turbulent closure model and inlet profiles is discussed.

Pordal, H. S.↗
DOE OSTI2014

Scientific Discovery through Advanced Computing

My work during the past year has involved application and development of the M3D code. M3D (Multi-level 3D) carries out simulation studies of plasmas using multiple levels of physics, geometry, and grid models. M3D combines a two dimensional unstructured mesh with finite element discretization in poloidal planes, with a pseudo spectral representation in the toroidal direction. Two critical MHD issues have emerged with the ITER effort: the suppression of edge localized modes (ELMs) by resonant magnetic perturbations (RMPs) and disruptions. Large ELMs can cause an unacceptable heat and particle flux to the wall and divertor. The application of RMPs to ITER is now being planned, although there are still some questions about how RMPs work. Disruptions can produce a much larger heat and particle than ELMs, and in addition can cause unacceptable mechanical stress on the wall and blanket. I have been carrying out extended MHD simulations which address these issues.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Continuum kinetic modeling of the tokamak plasma edge

The first 4D (axisymmetric) high-order continuum gyrokinetic transport simulations that span the magnetic separatrix of a tokamak are presented. The modeling is performed with the COGENT code, which is distinguished by fourth-order finite-volume discretization combined with mapped multiblock grid technology to handle the strong anisotropy of plasma transport and the complex X-point divertor geometry with high accuracy. The calculations take into account the effects of fully nonlinear Fokker-Plank collisions, electrostatic potential variations, and anomalous radial transport. Topics discussed include: (a) ion orbit loss and the associated toroidal rotation and (b) edge plasma relaxation in the presence of anomalous radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development and characterization of a laser-gated, high resolution x-ray radiography platform for high energy density experiments using toroidally bent crystals

Bent crystal x-ray imagers are a well-established diagnostic tool to study high energy density (HED) objects by acquiring two-dimensional x-ray radiographs. Often, studying these HED objects requires very high spatial resolution, which is limited by astigmatism when using spherically bent crystals. By using toroidal-shaped crystals instead, astigmatism can be reduced and the overall spatial resolution of the instrument improved. Here, the development and characterization of a laser-gated x-ray radiography platform at the National Ignition Facility using a toroidal crystal is presented. This includes measurements of the spatial and temporal profile, which is determined by the x-ray source. In conclusion, to properly validate the reduced astigmatism, a new multi-plane grid approach was implemented, which, for the first time, allows a depth of field measurement, alignment verification and x-ray source size estimate in a single shot, demonstrating the current platform’s capability to provide sub-10 μm resolution over a 0.9 mm depth of field with a temporal resolution of 110 ps.

Crystal optics↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Functional Plasma-Sprayed Ceramic Coatings for Use in Fusion Energy Systems

There is growing demand for new, cleaner ways to supply the rapidly growing global energy demand, as concerns grow over climate change and declining supplies of fossil fuels. Nuclear fusion, the phenomenon in which atomic nuclei collide together and release massive quantities of energy, is one of the most promising options for generating large amounts of carbon-free energy in the future.1 To initiate a common nuclear fusion reaction of hydrogen isotopes – deuterium and tritium, a hot plasma state is needed at incredibly high temperatures approaching 100 million degrees centigrade. One practical way to achieve these conditions is to use magnetic confinement of the plasma in a toroidal nuclear reactor known as a tokamak (Figure 1). A commercial fusion power station will then harvest the energy by converting the heat provided by the fusion chain reaction into steam to power gas turbines, and generate electricity on to the grid, while the excess neutrons also produced in the process is captured by a blanket made of dense material surrounding the tokomak.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Characterization of fast magnetosonic waves driven by compact toroid plasma injection along a magnetic field

Magnetosonic waves are low-frequency, linearly polarized magnetohydrodynamic (MHD) waves commonly found in space, responsible for many well-known features, such as heating of the solar corona. In this work, we report observations of interesting wave signatures driven by injecting compact toroid (CT) plasmas into a static Helmholtz magnetic field at the Big Red Ball Facility at Wisconsin Plasma Physics Laboratory. By comparing the experimental results with the MHD theory, we identify that these waves are the fast magnetosonic modes propagating perpendicular to the background magnetic field. Additionally, we further investigate how the background field, preapplied poloidal magnetic flux in the CT injector, and the coarse grid placed in the chamber affect the characteristics of the waves. Since this experiment is part of an ongoing effort of creating a target plasma with tangled magnetic fields as a novel fusion fuel for magneto-inertial fusion (MIF), our current results could shed light on future possible paths of forming such a target for MIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

NASA Tech Briefs, October 2003

Topics covered include: Cryogenic Temperature-Gradient Foam/Substrate Tensile Tester; Flight Test of an Intelligent Flight-Control System; Slat Heater Boxes for Thermal Vacuum Testing; System for Testing Thermal Insulation of Pipes; Electrical-Impedance-Based Ice-Thickness Gauges; Simulation System for Training in Laparoscopic Surgery; Flasher Powered by Photovoltaic Cells and Ultracapacitors; Improved Autoassociative Neural Networks; Toroidal-Core Microinductors Biased by Permanent Magnets; Using Correlated Photons to Suppress Background Noise; Atmospheric-Fade-Tolerant Tracking and Pointing in Wireless Optical Communication; Curved Focal-Plane Arrays Using Back-Illuminated High-Purity Photodetectors; Software for Displaying Data from Planetary Rovers; Software for Refining or Coarsening Computational Grids; Software for Diagnosis of Multiple Coordinated Spacecraft; Software Helps Retrieve Information Relevant to the User; Software for Simulating a Complex Robot; Software for Planning Scientific Activities on Mars; Software for Training in Pre-College Mathematics; Switching and Rectification in Carbon-Nanotube Junctions; Scandia-and-Yttria-Stabilized Zirconia for Thermal Barriers; Environmentally Safer, Less Toxic Fire-Extinguishing Agents; Multiaxial Temperature- and Time-Dependent Failure Model; Cloverleaf Vibratory Microgyroscope with Integrated Post; Single-Vector Calibration of Wind-Tunnel Force Balances; Microgyroscope with Vibrating Post as Rotation Transducer; Continuous Tuning and Calibration of Vibratory Gyroscopes; Compact, Pneumatically Actuated Filter Shuttle; Improved Bearingless Switched-Reluctance Motor; Fluorescent Quantum Dots for Biological Labeling; Growing Three-Dimensional Corneal Tissue in a Bioreactor; Scanning Tunneling Optical Resonance Microscopy; The Micro-Arcsecond Metrology Testbed; Detecting Moving Targets by Use of Soliton Resonances; and Finite-Element Methods for Real-Time Simulation of Surgery.

Source record↗
DOE OSTI2021

A Broad Grid of 2D Kilonova Emission Models

Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of 2-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1 M ⊙ and 0.05 to 0.3c, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (Y e ) ejecta with a robust r-process composition and either a spherical or lobed high-Y e ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code SuperNu using a full suite of lanthanide and 4th row element opacities. We examine the trends of these models with parameter variation, show how it can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

ATEQ: Adaptive toroidal equilibrium code

A radially adaptive numerical scheme is developed to solve the Grad–Shafranov equation for axisymmetric magnetohydrodynamic equilibrium. Here, a decomposition with independent solutions is employed in the radial direction, and Fourier decomposition is used in the poloidal direction. The independent solutions are then obtained using an adaptive shooting scheme together with the multi-region matching technique in the radial direction. Accordingly, the adaptive toroidal equilibrium (ATEQ) code is constructed for axisymmetric equilibrium studies. The adaptive numerical scheme in the radial direction improves considerably the accuracy of the equilibrium solution. The decomposition with independent solutions effectively reduces the matrix size in solving the magnetohydrodynamic equilibrium problem. The reduction of the matrix size is about an order of magnitude as compared with the conventional radially grid-based numerical schemes. Also, in this ATEQ numerical scheme, no matter how accuracy in the radial direction is imposed, the size of matrices basically does not change. The small matrix size scheme gives ATEQ more flexibility to address the requirement of the number of Fourier components in the poloidal direction in tough equilibrium problems. These two unique features, the adaptive shooting and small matrix size, make ATEQ useful to improve tokamak equilibrium solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing with 3D Adaptive Mesh Refinement and Local Adaptive Time Stepping

General relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multifaceted optimizations that, collectively, accelerate computation by 2–5 orders of magnitude for a wide range of applications. First, it introduces a spherical grid with 3D adaptive mesh refinement that operates in each of the three dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-resolution computational performance. Second, we demonstrate that local adaptive time stepping on a logarithmic spherical-polar grid accelerates computation by a factor of ≲10 compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of ~10 9 zone cycles per second per node on 5400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We illustrate H-AMR's computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of 13,440 × 4608 × 8092 cells and a total of ≲22 billion cells and ~0.65 × 10 8 time steps, it is among the largest astrophysical simulations ever performed. We find that frame dragging by the black hole tears up the disk into two independently precessing subdisks. The innermost subdisk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Gyroaveraging operations using adaptive matrix operators

A new adaptive scheme to be used in particle-in-cell codes for carrying out gyroaveraging operations with matrices is presented. This new scheme uses an intermediate velocity grid whose resolution is adapted to the local thermal Larmor radius. The charge density is computed by projecting marker weights in a field-line following manner while preserving the adiabatic magnetic moment μ. These choices permit to improve the accuracy of the gyroaveraging operations performed with matrices even when strong spatial variation of temperature and magnetic field is present. Accuracy of the scheme in different geometries from simple 2D slab geometry to realistic 3D toroidal equilibrium has been studied. As a result, a successful implementation in the gyrokinetic code XGC is presented in the delta-f limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

EPSI 2013-2014. Annual Progress Report

Initially a separate orthogonal expansion for the low-k components and FEM for high-k components was considered, or a multi-scale FEM method, but a more attractive hybrid FDM/FEM model for parallel computations has been devised with minimal impact to formerly developed code. This system involves using Lagrange polynomials along magnetic field lines to establish a FDM discretization of the toroidal contributions to the perpendicular Laplacian, which is the operator of importance in gyrokinetic solution of the vector potential. The method was found analytically to be stable, due to some dissipation from varying field-line lengths. The accuracy of this method has been evaluated using the method of manufactured solutions (with Bei Wang of Princeton University), whereby we find that for a Poisson problem the maximum error roughly diminishes linearly as a function of the grid resolution (uniformly refined), whereas the L2 norm as square root of the resolution improvement. Most of the error appears to be contributed near boundaries and where the FEM grid has anisotropy, so we are trying to improve this analysis using improved tessellations obtained from Mark Shepard, and publish the results as soon as the problems with the grid have been eliminated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Comparison of methods for numerical calculation of continuum damping

Continuum resonance damping is an important factor in determining the stability of certain global modes in fusion plasmas. A number of analytic and numerical approaches have been developed to compute this damping, particularly, in the case of the toroidicity-induced shear Alfvén eigenmode. This paper compares results obtained using an analytical perturbative approach with those found using resistive and complex contour numerical approaches. It is found that the perturbative method does not provide accurate agreement with reliable numerical methods for the range of parameters examined. This discrepancy exists even in the limit where damping approaches zero. When the perturbative technique is implemented using a standard finite element method, the damping estimate fails to converge with radial grid resolution. The finite elements used cannot accurately represent the eigenmode in the region of the continuum resonance, regardless of the number of radial grid points used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of an implicit electromagnetic capability for a hybrid gyrokinetic ion-fluid electron model

Here we report on the development and implementation of a hybrid kinetic ion–fluid electron model for electromagnetic COGENT simulations of edge plasmas. COGENT is a finite-volume gyrokinetic code that employs a locally field-aligned coordinate system combined with a mapped multi-block grid technology to handle strongly anisotropic edge plasma turbulence. The simulation model involves the long-wavelength limit of the ion gyrokinetic equation coupled to the vorticity and Ohm's law equations for the electromagnetic field perturbations. In order to handle the fast Alfvén wave time scales, an implicit-explicit time integration approach with a physics-based preconditioner is used. The model is successfully applied to the simulations of ion-scale resistive-drift ballooning turbulence in a toroidal annulus geometry. Substantial speed-up over a fully explicit time integration approach is observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗