Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface fluxes”

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 343 records · Page 19

Low-frequency oscillations in radiative-convective systems

Although eastward propagation has long been considered one of the essential features of the Madden-Julian waves, recent observations have revealed a stationary or quasi-stationary component in the oscillations, particularly in measures of the diabatic heating rate. Wave-CISK theories of the low-frequency oscillations have struggled to explain the observed period and vertical structure of the waves. On the other hand, theoretical and numerical studies have shown that low-frequency waves strongly resembling the observed oscillations can be excited by specified low-frequency oscillations of the convective heating. A problem with the latter set of theories is that the cause of the oscillatory heating has not been satisfactorily explained. It is proposed here that the observed low-frequency wave motions are the response to forcing by an essentially stationary, self-excited oscillating heat source that is produced by nonlinear interactions among radiation, cumulus convection, and the surface fluxes of sensible heat and moisture. Feedback of the large-scale motions on the latent heating is not required. Results from two very different one-dimensional models are presented to support this hypothesis. The physical processes included in the models are essentially the same, that is, radiation, cumulus convection, and the surface fluxes of sensible heat and moisture; the first model is highly simplified, however, while the second includes relatively sophisticated parameterizations of all the relevant physical processes. Results from both models show low-frequency oscillations of the latent heating, temperature, and moisture. Experiments show that the oscillations are favored by a warm sea surface and weak surface wind speeds, consistent with the observed conditions over the Indian Ocean and the tropical western Pacific Ocean.

Hu, QI↗

Force-free fields in thin coronal loops

We solve the force-free equation J x B = 0 for fields which are toroidally symmetric. The technique utilizes an expansion about a cylindrical field and is therefore valid or tori with a large aspect ratio such as long, thin, coronal loops. The calculation is performed in spatial toroidal coordinates, rather than in the flux coordinates used by previous authors; this allows direct calculation of the loci of flux surfaces and of surfaces of constant magnetic pressure. Our solutions differ significantly from toroidal fields in laboratories, which are in general not force-free. They are characterized by field lines whose projections in the poloidal planes are circles with centers displaced by varying distances from the axis of the torus. In general, flux surfaces do not correspond to surfaces of constant magnetic pressure. We have examined solutions corresponding to simple analytic zero-order cylindrical fields. For moderate twists in the zero-order (cylindrical) field, the magnetic pressure is larger on the inner toroidal radius. However, this effect diminishes with twist angle and in fact, for extreme initial twists, the magnetic pressure can be larger on the outer radius. We compare our results with previous work utilizing flux coordinates.

Emslie, A. Gordon↗

Geometric dependencies of the mean E × B shearing rate in negative triangularity tokamaks

Abstract This paper presents a comparative study of the poloidal distribution of the mean E × B shearing rate for positive triangularity (PT) and negative triangularity (NT) tokamaks. The effects of flux surface up–down asymmetry due to asymmetric upper and lower triangularities are also considered. Both direct eddy straining and the effects on Shafranov shift feedback loops are examined. Shafranov shift increases the shearing rate at all poloidal angles for all triangularities, due to flux surface compression. The maximum shearing rate bifurcates at a critical triangularity δ crit ( ≲ 0 ) . Thus, the shearing rate is maximal off the outboard midplane for NT, while it is maximal on the outboard midplane for PT. For up–down asymmetric triangularity, the usual up–down symmetry of the shearing rate is broken. The shearing rate at the outboard midplane is lower for NT than for PT, suggesting that the shearing efficiency in NT is reduced. Implications for turbulence stabilisation and confinement improvement in high- β p NT and internal transport barrier discharges are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stellar chromospheres and coronae in the Ursa Major Cluster stars

IUE spectra of 18 proposed members of the Ursa Major Cluster and Einstein X-ray images of 11 of these stars are discussed. Thirteen stars, six in the nucleus and seven in the extended comoving stream, are probably true members of the UMa Cluster in that their bright ultraviolet and X-ray surface fluxes indicate youth. Four stars, one in the nucleus and three in the stream, exhibit weak ultraviolet and/or X-ray emission suggesting that they are old field stars that have the same space motion as the UMa Cluster. The X-ray surface fluxes of the UMa Cluster stars appear to be brighter than those of the Hyades Cluster stars, consistent with their relative ages. It is argued that chromospheres and transition regions could be present in dwarf stars hotter than B-V of about 0.30, but are unobservable in IUE spectra due to the rapid increase in photospheric flux in hotter stars.

Simon, T.↗

Analysis of the Effects of Vitiates on Surface Heat Flux in Ground Tests of Hypersonic Vehicles

To achieve the high enthalpy conditions associated with hypersonic flight, many ground test facilities burn fuel in the air upstream of the test chamber. Unfortunately, the products of combustion contaminate the test gas and alter gas properties and the heat fluxes associated with aerodynamic heating. The difference in the heating rates between clean air and a vitiated test medium needs to be understood so that the thermal management system for hypersonic vehicles can be properly designed. This is particularly important for advanced hypersonic vehicle concepts powered by air-breathing propulsion systems that couple cooling requirements, fuel flow rates, and combustor performance by flowing fuel through sub-surface cooling passages to cool engine components and preheat the fuel prior to combustion. An analytical investigation was performed comparing clean air to a gas vitiated with methane/oxygen combustion products to determine if variations in gas properties contributed to changes in predicted heat flux. This investigation started with simple relationships, evolved into writing an engineering-level code, and ended with running a series of CFD cases. It was noted that it is not possible to simultaneously match all of the gas properties between clean and vitiated test gases. A study was then conducted selecting various combinations of freestream properties for a vitiated test gas that matched clean air values to determine which combination of parameters affected the computed heat transfer the least. The best combination of properties to match was the free-stream total sensible enthalpy, dynamic pressure, and either the velocity or Mach number. This combination yielded only a 2% difference in heating. Other combinations showed departures of up to 10% in the heat flux estimate.

Cuda, Vincent↗

LDEF environment modeling updates

An updated gas dynamics model for gas interactions around the LDEF is presented that includes improved scattering algorithms. The primary improvement is more accurate predictions of surface fluxes in the wake region. The code used is the Integrated Spacecraft Environments Model (ISEM). Additionally, initial results of a detailed ISEM prediction model of the Solar Array Passive LDEF Experiment (SAMPLE), A0171, is presented. This model includes details of the A0171 geometry and outgassing characteristics of the many surfaces on the experiment. The detailed model includes the multiple scattering that exists between the ambient atmosphere, LDEF outgassing, and atomic oxygen erosion products. Predictions are made for gas densities, surface fluxes and deposition at three different time periods of the LDEF mission.

Gordon, Tim↗

Simulations of plasma current induced by toroidal field ramping down on tokamaks

In tokamak experiments, the BT ramping down can drive off-axis and parallel inductive current. This approach leads to a decrease in q95 and an increase in normalized beta, βN. The off-axis inductive current also broadens the current profile. Typically, the 1.5D transport code is used to simulate the time evolution of the plasma current profile, which is based on the flux-surface-averaged Faraday’s law. The ONETWO code is one of those transport codes. However, this code cannot simulate the situation of the evolving toroidal field BT. This article proposes a new time-dependent model to take into account the BT ramping down situation. A modified formula of flux-surface-averaged Faraday’s law was derived to consider the effect of the BT ramping rate on the current evolution, and it was implemented in the ONETWO code. Then, the modified ONETWO code was used to simulate the BT ramping down experiment on DIII-D. The simulation result of the plasma current evolution with BT ramping down shows a broader current profile with smaller ohmic current induced by the poloidal field, compared to that without BT ramping down.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast physics-based launcher optimization for electron cyclotron current drive

With the increased urgency to design fusion pilot plants, fast optimization of electron cyclotron current drive (ECCD) launchers is paramount. Traditionally, this is done by coarsely sampling the 4D parameter space of possible launch conditions consisting of (1) the launch location (constrained to lie along the reactor vessel), (2) the launch frequency, (3) the toroidal launch angle, and (4) the poloidal launch angle. For each initial condition, a ray-tracing simulation is performed to evaluate the ECCD efficiency. Unfortunately, this approach often requires a large number of simulations (sometimes millions in extreme cases) to build up a dataset that adequately covers the plasma volume, which must then be repeated every time the design point changes. Here we adopt a different approach. Rather than launching rays from the plasma periphery and hoping for the best, we instead directly reconstruct the optimal ray for driving current at a given flux surface using a reduced physics model coupled with a commercial ray-tracing code. Repeating this throughout the plasma volume requires only hundreds of simulations, constituting a significant speedup. The new method is validated on two separate example tokamak profiles, and is shown to reliably drive localized current at the specified flux surface with the same optimal efficiency as obtained from the traditional approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of implementing the Simple Biosphere Model in a general circulation model

The Simple Biosphere Model (SiB) of Sellers et al., was designed to simulate the interactions between the earth's land surface and the atmosphere by treating the vegetation explicitly and realistically, thereby incorporating the biophysical controls on the exchanges of radiation, momentum, sensible and latent heat between the two systems. This paper describes the steps taken to implement SiB in a modified version of the National Meteorological Center's global spectral general circulation model (GCM) and explores the impact of the implementation on the simulated land surface fluxes and near-surface meteorological conditions. The coupled model (SiB-GCM) was used to produce summer and winter simulations. The same GCM was used with a conventional hydrological model (Ctl-GCM) to produce comparable 'control' summer and winter simulations for comparison. It was found that SiB-GCM produced a more realistic partitioning of energy at the land surface than Ctl-GCM. Generally, SiB-GCM, produced more sensible heat flux and less latent heat flux over vegetated land than did Ctl-GCM and this resulted in a much deeper daytime planetary boundary layer and reduced precipitation rates over the continents in SiB-GCM. In the summer simulation, the 200 mb jet stream was slightly weakened in the SiB-GCM relative to the Ctl-GCM results and analyses made from observations.

Sato, N.↗

Visualization techniques for the gyrokinetic tokamak simulation code

Gyrokinetic simulations of plasma microturbulence in tokamaks are challenging to visualize because the compute grid follows the magnetic field lines that spiral around the torus. We have overcome this challenge by developing three new approaches that improve visualization of gyrokinetics. Our techniques work directly with the topology of magnetic flux surfaces where the simulation stores variables in concentric rings on poloidal planes (vertical cross sections of the torus). Our visualization preview step triangulates each consecutive pair of rings to display the data on a poloidal plane. The second visualization technique follows spiral field lines around the torus and constructs polygons to visualize a flux surface. Third, the poloidal triangles are connected between planes to form prisms that compose a 3-D model of the entire torus. The visualization workflow produces detailed geometry that matches the high resolution, irregular compute grid for every time step. The surface and solid models are displayed in scientific visualization programs to effectively explore and communicate the results, including fluctuation of electron density, ion temperature, and electrostatic potential. Highly detailed renderings verify plasma behavior along magnetic field lines over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New linear stability parameter to describe low-β electromagnetic microinstabilities driven by passing electrons in axisymmetric toroidal geometry

In magnetic confinement fusion devices, the ratio of the plasma pressure to the magnetic field energy, β, can become sufficiently large that electromagnetic microinstabilities become unstable, driving turbulence that distorts or reconnects the equilibrium magnetic field. In this paper, a theory is proposed for electromagnetic, electron-driven linear instabilities that have current layers localised to mode-rational surfaces and binormal wavelengths comparable to the ion gyroradius. The model retains axisymmetric toroidal geometry with arbitrary shaping, and consists of orbit-averaged equations for the mode-rational surface layer, with a ballooning space kinetic matching condition for passing electrons. The matching condition connects the current layer to the large scale electromagnetic fluctuations, and is derived in the limit that β is comparable to the square root of the electron-to-ion-mass ratio. Electromagnetic fluctuations only enter through the matching condition, allowing for the identification of an effective β that includes the effects of equilibrium flux surface shaping. The scaling predictions made by the asymptotic theory are tested with comparisons to results from linear simulations of micro-tearing and electrostatic microinstabilities in MAST discharge #6252, showing excellent agreement. In particular, it is demonstrated that the effective β can explain the dependence of the local micro-tearing mode (MTM) growth rate on the ballooning parameter θ 0 – possibly providing a route to optimise local flux surfaces for reduced MTM-driven transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantifying the Role of Ocean Dynamics in SST Variability across GCMs and Observations

Abstract Midlatitude SSTs forced by mesoscale oceanic processes can affect the large-scale atmosphere, pointing to the ocean’s crucial role outside the tropics. Previous studies have shown oceanic mesoscale processes’ effect on global and regional climate variability. This study quantifies the local contribution of ocean dynamics to mixed-layer temperature across the globe by directly estimating the ocean heat flux divergence resolved by state-of-the-art ocean reanalysis, eddy-resolving, and eddy-parameterized versions of two U.S. national climate models and indirectly from air–sea flux satellite-based estimates. Our results show that the eddy-resolving climate simulations resolve mixed-layer temperature variances that are larger and closer to those inferred from observations than both their eddy-parameterized counterparts and ECCO over much of the extratropics. The observations and the eddy-resolving models indicate a more significant role of ocean dynamics in the mixed-layer temperature variability than the surface fluxes over most extratropics compared to their eddy-parameterized versions. A frequency domain analysis shows that the better-resolved ocean mesoscale and thermal gradients enhance the variance over a time scale from 2 months to 30 years. Results show agreement in the ocean’s contribution among satellite-based estimates, ocean reanalysis products, and ocean eddy-resolving simulations. At the same time, differences emerge for ECCO and the eddy-parameterized models, suggesting that surface fluxes account for a larger fraction of the mixed-layer temperature variability in most of the extratropics.

Siqueira, L.↗

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↗

Field lines and magnetic surfaces in a two-component slab/2D model of interplanetary magnetic fluctuations

A two-component model for the spectrum of interplanetary magnetic fluctuations was proposed on the basis of ISEE observations, and has found an intriguing level of application in other solar wind studies. The model fluctuations consist of a fraction of 'slab' fluctuations, varying only in the direction parallel to the locally uniform mean magnetic field B(0) and a complement of 2D (two-dimensional) fluctuations that vary in the directions transverse to B(0). We have developed an spectral method computational algorithm for computing the magnetic flux surfaces (flux tubes) associated with the composite model, based upon a precise analogy with equations for ideal transport of a passive scalar in planar two dimensional geometry. Visualization of various composite models will be presented, including the 80 percent 2D/ 20 percent slab model with delta B/B(0) approximately equals 1 and a minus 5/3 spectral law, that is thought to approximately represent a snapshot of solar wind turbulence. Characteristically, the visualizations show that flux tubes, even when defined as regular on some plane, shred and disperse rapidly as they are viewed along the parallel direction. This diffusive process, which generalizes the standard picture of field line random walk, will be discussed in detail. Evidently, the traditional picture that flux tubes randomize like strands of spaghetti with a uniform tangle along the axial direction is in need of modification.

Matthaeus, W. H.↗

Core plasma fueling by fast inward particle transport after hydrogen pellet injection in Wendelstein 7-X

A large database of more than 1000 individual cryogenic hydrogen pellets injected into Wendelstein 7-X for plasma fueling was analyzed to improve the understanding of the three phases of the process: the ablation, deposition and transport of the pellet material. Kilohertz-sampled electron density and temperature measurements revealed a more complex drift behavior than predicted by numerical code simulation. It could be explained by the poloidal plasma E r x B- drift rotation, which plays a significant role in stellarators, but was not previously considered in pellet injection codes like HPI2. The drift results in a fast poloidal rotation of the pellet material around the plasma core, leading to an almost homogeneous deposition over the involved flux surfaces regardless of magnetic high and low field side injection geometry. Additionally, a novel fast inward directed transport mechanism (‘FIT-effect’) was observed. The effect occurs on timescales of tens of milliseconds and cannot be explained by neoclassical transport or diffusion. It might be linked to the turbulence pinch recently found in Wendelstein 7-X. When the FIT-effect occurs, the pellet particles are rapidly transferred from the deposition flux surfaces to the plasma core, causing the plasma density profile to peak, which is beneficial for confinement in Wendelstein 7-X. The large pellet injection database was statistical analyzed with regard to pellet and plasma parameters, which delivered some starting points towards developing an understanding of the physics behind the FIT-effect. The results indicate, that plasma core fueling via pellet injection is largely independent of the injection geometry in stellarators under certain conditions, reducing the technical complexity of the injection system.

Wendelstein 7-X↗

Nonparametric solutions to the variational principle of ideal magnetohydrodynamics

In an effort to gain a better understanding of MHD equilibria in three dimensions, the lower dimensional cases are studied. The solution of the three-dimensional problem is based on the classical variational principle of ideal magnetohydrodynamics. The crucial assumption for the numerical method is the existence of a nested set of toroidal flux surfaces, which is then used as a coordinate. This paper studies the nonparametric solutions to this variational problem in those cases when the direct solution is known to have islands. A form of the variational principle for the slab geometry is described; the one-dimensional problem is analyzed; and asymptotic expansions and numerical solutions to the two-dimensional problem are discussed. An example is presented which shows that the assumption of nested flux surfaces need not rule out the occurrence of islands.

Betancourt, O. L.↗

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↗

On the enigma of FK Comae

Stellar chromospheric and coronal activity appears ubiquitous among late type stars to the left of the TR-wind boundary line (Linksy and Haisch 1979). The level of activity as measured by the X-ray surface flux is linearly proportional to the stellar angular velocity, with the exception of slowly rotating dwarfs (Walter 1981, 1982; Walter and Bowyer 1981). The peculiar rapidly rotating G giant FK Comae (Merrill 1948) appears to fit into this pattern. Line widths indicate V sin i = 120 + or - 20 km s(-1) (Bopp and Stencel 1981). FK Comae has strong Ca II H and K and H alpha emission, strong transition region UV lines (Bopp and Stencel 1981), and an X-ray surface flux in good agreement with its rapid rotation (Walter 1981). Yet, FK Comae is an enigmatic star. It is a rapid rotator, but it is not clear why it is a rapid rotator. There is no direct evidence for duplicity; indeed, the upper limit of 20 km s(-1) on the K velocity puts tight constraints on any binary configuration, especially if sin i approx 1, as indicated by the large V sin i. Bopp and Stencel (1981) have suggested that FK Comae is an example of a coalesced W UMa system (Webbink 1976), wherein the orbital angular momentum has become rotational angular momentum of the coalesced star.

Walter, F.↗