Search NASA⌕ Search

SEARCH · Search NASA

Results for “field line tracing”

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 OSTI2021

Model-based real-time surface heat flux and temperature estimation for the DIII-D tokamak

A control-oriented model for monitoring of wall power flux densities on the DIII-D tokamak has been successfully implemented and validated experimentally. Future reactors will have to withstand severe steady state high heat flux loads on plasma-facing components (PFCs). Due to the difficulty of directly-measuring local heat fluxes on these components, monitoring and protection of PFCs during the plasma discharge can benefit from simplified physics-based real-time (RT) functional models to estimate and guide heat load control. As a first step into the development, a control-oriented model for monitoring of wall power flux densities and temperatures on DIII-D tokamak has been successfully implemented. The paper discusses the experimental demonstration and comparison of the 2-D model-based wall heat flux algorithm on the DIII-Dinner wall limiter (IWL) against infra-red camera heat flux measurements for limited plasma configurations. The paper also reports on the benchmarking of the field line tracing environment, SMITER, developed at ITER organization on DIII-D tokamak against experimental IR diagnostic data and the derivation of the component shaping weighting factors for the 2-D model-based approach. Here, the extension of the model-based approach for surface temperature estimation on the DIII-D IWL is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

X-point and divertor filament dynamics from gas puff imaging on TCV

Abstract A new gas puff imaging diagnostic has been installed on the TCV tokamak, providing two-dimensional insights into scrape-off-layer (SOL) turbulence dynamics above, at and below the magnetic X-point. A detailed study in L-mode, attached, lower single-null discharges shows that statistical properties have little poloidal variations, while vast differences are present in the 2D behaviour of intermittent filaments. Strongly elongated filaments, just above the X-point and in the divertor far-SOL, show a good consistency in shape and dynamics with field-line tracing from filaments at the outboard midplane, highlighting their connection. In the near-SOL of the outer divertor leg, short-lived, high frequency and more circular (diameter ∼15 sound Larmor radii) filaments are observed. These divertor-localised filaments appear born radially at the position of maximum density and display a radially outward motion with velocity ≈400 m s −1 that is comparable to radial velocities of upstream-connected filaments. Conversely, in these discharges ( B → × ∇ B pointing away from the divertor), these divertor filaments’ poloidal velocities differ strongly from those of upstream-connected filaments. The importance of divertor-localised filaments upon radial transport and profile broadening is explored using filament statistics and in situ kinetic profile measurements along the divertor leg. This suggests that these filaments contribute significantly to electron density profile broadening in the divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development and validation of non-axisymmetric heat flux simulations with 3D fields using the HEAT code

A new comprehensive module to simulate heat fluxes from three-dimensional (3D) magnetic fields has been implemented in the HEAT code. Especially compact tokamaks like SPARC require tools to predict and manage large heat fluxes. Existing release versions of HEAT can only simulate axisymmetric heat flux on 3D plasma facing components. The new module uses an M3D-C1 perturbed equilibrium and the MAFOT code to trace field lines of the perturbed 3D magnetic field. Heat flux is then assigned to the resulting footprints via a 3D layer model. The model distinguishes between the scrape-off layer, the magnetic lobes and the private flux region, and employs only 0D parameters like the layer width, diffusive spread and the last closed flux surface position in the perturbed edge to generate a heat flux profile. The magnitude is normalized to the total input power. Resulting heat flux simulations are compared and validated against infrared measurements in the DIII-D tokamak with applied 3D fields; good agreement is found for several cases. The new module can now be applied to the SPARC tokamak; a preliminary result for applied rotating 3D fields is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Observation of stationary filaments with resonant magnetic perturbations in NSTX

The application of resonant magnetic perturbations (RMPs) induces striated heat fluxes on the divertor plate in tokamaks. These striated heat fluxes were previously considered strike point splitting by the intersection between the divertor and predicted magnetic lobe structures. Recently, the existence of stationary 3D filaments has been observed for the first time in National Spherical Torus Experiment (NSTX) plasmas with toroidal mode number n = 3 RMPs. Magnetic field line tracing of these filaments from the midplane to the divertor is consistent with the quantity and the radial and toroidal location of striated heat fluxes observed on the divertor. This study suggests that divertor striated heat flux patterns are linked to these stationary filaments rather than to strike point splitting or magnetic lobe structures for NSTX. This finding will help guide future efforts to understand the heat flux striations in NSTX with RMPs in order to control and reduce the divertor peak heat flux.

RMPs↗
DOE OSTI2026

Effect of controlled magnetic island bifurcation on electron diffusion

Magnetic islands strongly influence cross-field electron transport in magnetized plasmas. In particular, bifurcations of the island topology modify the number and location of O-points, X-points, and separatrix boundaries, thereby altering diffusion pathways. In recent DIII-D experiments, external magnetic perturbations were used to rotate and periodically bifurcate the island on the q = 2 surface, causing a switchback between a q = 2/1-dominated structure and a narrower q = 4/2-dominated structure. To investigate how this topological change affects electron transport, we employ the field line tracing code TRIP3D with an implemented collisional operator. Thermal, tracer electrons launched from O-points, X-points, and outside separatrix boundaries reveal distinct diffusion regimes, including classical, subdiffusive, and superdiffusive behavior, depending on both the dominant island mode and launch location. These results suggest that island bifurcation can alter electron diffusion across rational surfaces, with direct implications for particle confinement. While the present work emphasizes diffusion as a general framework, the findings provide insight into the conditions under which electron trapping into an island or stochastization of the island's separatrix can enable additional mechanisms, such as the generation of energetic electrons.

DIII-D↗
DOE OSTI2022

Reassessing energy deposition for the ITER 5 MA vertical displacement event with an improved DINA model

The beryllium (Be) main chamber wall interaction during a 5 MA/1.8 T upward, unmitigated VDE scenario, first analysed in [J. Coburn et al., Phys. Scr. T171 (2020) 014076] for ITER, has been re-evaluated using the latest energy deposition analysis software. Updates to the DINA disruption model are summarized, including an improved numerical convergence for the 0D power balance, limitations on the safety factor within the plasma core, and the choice to maintain a constant plasma +halo poloidal cross-section. Such updates result in a broad halo region and higher radiated power fractions compared to previous models. The new scenario lasts for ~75 ms and deposits ~29 MJ of energy, with the radial distribution of parallel heat flux q‖(r)resembling an exponential falloff with an effective λ_q=75-198 mm. A maximum halo width w_h of 0.52 m at the outboard midplane is observed. SMITER field line tracing and energy deposition simulations calculate a q_{⊥,max} of ~83 MW/m^2 on the upper first wall panels (FWP). Heat transfer calculations with the MEMOS-U code show that the FWP surface temperature reaches ~1000 K, well below the Be melt threshold. Variations of this 5 MA scenario with Be im-purity densities from 0 to 3∙10^19 m^-3 also remain below the melt threshold despite differences in energy deposition and duration. These results are in contrast to the early study which predicted melt damage to the first wall [J. Coburn et al., Phys. Scr. T171 (2020) 014076], and emphasize the importance of accurate models for the halo width w_h and the heat flux distribution q‖(r)within that halo width. The 2020 halo model in DINA has been compared with halo current experiments on COMPASS, JET, and Alcator C-Mod, and the preliminary results build confidence in the broad halo width predictions. Results for the 5 MA VDE are compared with those for a 15 MA equivalent, generated using the new DINA model. At the higher current, significant melting of the upper FWP is to be expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

Magnetospheric mapping with quantitative geomagnetic field models

The Mead-Fairfield geomagnetic field models were used to trace field lines between the outer magnetosphere and the earth's surface. The results are presented in terms of ground latitude and local time contours projected to the equatorial plane and into the geomagnetic tail. With these contours various observations can be mapped along field lines between high and low altitudes. Low altitudes observations of the polar cap boundary, the polar cusp, the energetic electron trapping boundary and the sunward convection region are projected to the equatorial plane and compared with the results of the model and with each other. The results provide quantitative support to the earlier suggestions that the trapping boundary is associated with the last closed field line in the sunward hemisphere, the polar cusp is associated with the region of the last closed field line, and the polar cap projects to the geomagnetic tail and has a low latitude boundary corresponding to the last closed field line.

Fairfield, D. H.↗
NASA NTRS1978

Linear force-free fields in the lower corona

Boundary conditions are imposed on the normal component of B on two parallel planes which delineate the force-free volume of linear magnetic fields in low coronal levels of the solar atmosphere. Green's function is employed to ensure that equations for the magnetic field energy remain bounded and smooth. Field line tracings of a simple bipolar source distribution are developed and compared with empirical data obtained from visual observations in H(alpha), EUV, and X-rays.

Barbosa, D. D.↗
NASA NTRS1980

Spacelab energetic ion mass spectrometer

Basic design criteria are given for an ion mass spectrometer for use in studying magnetospheric ion populations. The proposed instrument is composed of an electrostatic analyzer followed by a magnetic spectrometer and simultaneously measures the energy per unit and mass per unit charge of the ion species. An electromagnet is used for momentum analysis to extend the operational energy range over a much wider domain than is possible with the permanent magnets used in previous flights. The energetic ion source regions, ion energization mechanisms, field line tracing, coordinated investigations, and orbit considerations are discussed and operations of the momentum analyzer and of the electrostatic energy analyzer are examined.

Whalen, B. A.↗
NASA NTRS1992

Probing the magnetosphere using chemical releases from the Combined Release and Radiation Effects Satellite

An overview is presented of the chemical release experiments from NASA's Combined Release and Radiation Effects Satellite (CRRES) program. Preliminary results are given for the CRRES investigations of (1) stimulated electron and ion precipitation, (2) ion transport in the magnetotail, (3) critical ionization velocity, (4) field line tracing and parallel acceleration, (5) diamagnetic cavity formation and collapse, and (6) plasma instabilities. The chemical vapor properties from a thermite release mechanism are also briefly described.

Bernhardt, P. A.↗
NASA NTRS1995

Interactive data exploration and particle tracking for general circulation models

The SCENE environment for interactive visualization of complex data sets is discussed. This environment is used to create tools for graphical exploration of atmospheric flow models. These tools may be extended by the user in a seamless manner, so that no programming is required. A module for accurately tracing field lines and particle trajectories in SCENE is presented. This is used to examine the flowfield qualitatively with streamlines and pathlines and to identify critical points in the velocity field. The paper also describes a visualization tool for general circulation models on which the primary features of the environment are demonstrated.

Rosenbaum, R. I.↗
NASA NTRSDate not supplied

Magnetic Mapping in the Inner Magnetosphere using Kamodo

Many models require specialized access and interpolation schemes to effectively extract and interpolate their outputs. In particular, the Block-Adaptive Tree Solarwind Roe Upwind Scheme (BATSRUS) component of the Space Weather Modeling Framework (SWMF) requires Kamodo to take advantage of its block-based adaptive grid structure, and the Lyon-Fedder Mobarry magnetosphere model (or its successor GAMERA) needs a scheme that appreciates the distorted spherical arrangement of grid vertices on a non-orthogonal grid. With the flythrough layer developed by Ringuette et al. (SH42E-2337), the underlying model readers have been adapted to use multiple time steps in a single Python session to perform 4- dimensional interpolations in time and space. Kamodo now utilizes lazy interpolation that loads data only when needed. We present the successful integration of SWMF/BATSRUS magnetosphere access and interpolation into the new 4D Kamodo framework utilizing an external library of C code. Through function composition, Kamodo facilitates the calculation of derived quantities and the transformation of positions and vectors into different coordinate systems. This work is a significant step towards performing field line tracing in Kamodo with SWMF magnetosphere outputs.

Lutz Rastaetter↗
NASA NTRS2023

Prediction and Verification of Parker Solar Probe Solar Wind Sources at 13.3 Rs

Drawing connections between heliospheric spacecraft and solar wind sources is a vital step in understanding the evolution of the solar corona into the solar wind and contextualizing in situ timeseries. Furthermore, making advanced predictions of this linkage for ongoing heliospheric missions, such as Parker Solar Probe (Parker), is necessary for achieving useful coordinated remote observations and maximizing scientific return. The general procedure for estimating such connectivity is straightforward (i.e., magnetic field line tracing in a coronal model) but validating the resulting estimates is difficult due to the lack of an independent ground truth and limited model constraints. In its most recent orbits, Parker has reached perihelia of 13.3R⊙ and moreover travels extremely fast prograde relative to the solar surface, covering over 120° longitude in 3 days. Here we present footpoint predictions and subsequent validation efforts for Parker Encounter 10, the first of the 13.3R⊙ orbits, which occurred in November 2021. We show that the longitudinal dependence of in situ plasma data from these novel orbits provides a powerful method of footpoint validation. With reference to other encounters, we also illustrate that the conditions under which source mapping is most accurate for near-ecliptic spacecraft (such as Parker) occur when solar activity is low, but also require that the heliospheric current sheet is strongly warped by mid-latitude or equatorial coronal holes. Lastly, we comment on the large-scale coronal structure implied by the Encounter 10 mapping, highlighting an empirical equatorial cut of the Alfvèn surface consisting of localized protrusions above unipolar magnetic separatrices.

S. T. Badman↗
NASA NTRSDate not supplied

On the Solar Open Magnetic Flux Problem

The solar magnetic fields emerging from the photosphere into the chromosphere and corona are comprised of a combination of “closed” (field lines with both ends rooted at the Sun) and “open” (field lines with only one end at the Sun) fields. Since the early 2000’s, the magnitude of total unsigned open magnetic flux estimated by coronal models has been in significant disagreement with in situ spacecraft observations, especially during solar maximum. Estimates of total open unsigned magnetic flux using coronal hole observations (e.g., using extreme ultraviolet (EUV) or Helium (He) are in general agreement with the coronal model results and thus show similar disagreements with in situ observations. This paper provides a brief overview of the problem, summarizes the proposed explanations for the discrepancies, and presents results that strongly support the explanation that the discrepancy is due to dynamics at the open-closed boundary. These results are derived from determination of the total unsigned open magnetic flux, utilizing the Wang-Sheeley-Arge (WSA) model at a particular spatial resolution and different field-line tracing methods. One of these methods produces excellent agreement with in situ observations. Our results imply that active regions residing near the boundaries of mid-latitude coronal holes are the primary source of the missing open flux. Furthermore, they resolve many of the seemingly contradictory facts that have made resolving the open-flux problem so difficult.

Charles Arge↗
NASA NTRS1972

ALLMAG, GDALMG, LINTRA: Computer programs for geomagnetic field and field-line calculations

A set of computer programs has been developed for the calculation of the geomagnetic field and the tracing of field lines in space. The basic subroutine, geocentric ALLMAG, contains coefficients for seven recently-published field models as built-in data statements. At execution time the user can vary the model and/or the time period by simply changing input parameters. Subroutine GDALMG is adapted for input and output in geodetic coordinates. ALLMAG and GDALMG are equivalent to Cain's FIELD and FIELDG, with the added flexibility of the choice of seven models. LINTRA traces field lines from any point in space to a specified altitude intersect in the same or opposite hemisphere, using any of the models contained in ALLMAG. Input is in either geocentric or geodetic coordinates, and output is returned in both. McIlwain's INVAR package, which calculates B and L, has been adapted to use ALLMAG. All programs are described in detail, and sample calculations are given.

Stassinopoulos, E. G.↗
NASA NTRS1988

Mapping of auroral kilometric radiation sources to the aurora

Auroral kilometric radiation (AKR) and optical auroral emissions are observed simultaneously using plasma wave instrumentation and auroral imaging photometers carried on the DE 1 spacecraft. The DE 1 plasma wave instrument measures the relative phase of signals from orthogonal electric dipole antennas, and from these measurements, apparent source directions can be determined with a high degree of precision. Wave data are analyzed for several strong AKR events, and source directions are determined for several emission frequencies. By assuming that the AKR originates at cyclotron resonant altitudes, a candidate source field line is identified. When the selected source field line is traced down to auroral altitudes on the concurrent DE 1 auroral image, a striking correspondece between the AKR source field line and localized auroral features is produced. The magnetic mapping study provides strong evidence that AKR sources occur on field lines associated with discrete auroral arcs, and it provides confirmation that AKR is generated near the electron cyclotron frequency.

Huff, R. L.↗
NASA NTRS1990

Magnetic field configuration of the theta aurora

A magnetic configuration of the open magnetosphere is described which is conducive to the formation of the theta aurora when the IMF has a significant northward component. A magnetic field topology and polar cap configuration, derived from a quantitative model of the open magnetosphere that incorporates Crooker's antiparallel merging hypothesis, are presented. Under this hypothesis, when the IMF has a northward component, the dayside merging line bifurcates, leaving a large fraction of the subsolar magnetopause untouched by the merging process. The polar cap, defined by tracing magnetic field lines that connect from the solar wind to the earth, is similarly bifurcated, leaving a sun-aligned stagnation region that is not magnetically connected to the solar wind and may plausibly be associated with the sun-aligned 'bar' of the theta aurora. The model provides testable predictions with regard to the position of this 'convection gap' in both Northern and Southern Hemispheres as functions of IMF direction.

Toffoletto, F. R.↗
NASA NTRS1975

Magnetospheric mapping with a quantitative geomagnetic field model

Mapping the magnetosphere on a dipole geomagnetic field model by projecting field and particle observations onto the model is described. High-latitude field lines are traced between the earth's surface and their intersection with either the equatorial plane or a cross section of the geomagnetic tail, and data from low-altitude orbiting satellites are projected along field lines to the outer magnetosphere. This procedure is analyzed, and the resultant mappings are illustrated. Extension of field lines into the geomagnetic tail and low-altitude determination of the polar cap and cusp are presented. It is noted that while there is good agreement among the various data, more particle measurements are necessary to clear up statistical uncertainties and to facilitate comparison of statistical models.

Fairfield, D. H.↗