SEARCH · Engineering Papers

Results for “reach code”

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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2024

SN 2021adxl: A luminous nearby interacting supernova in an extremely low-metallicity environment

SN 2021adxl is a slowly evolving, luminous, Type IIn supernova with asymmetric emission line profiles, similar to the well-studied SN 2010jl. We present extensive optical, near-ultraviolet, and near-infrared photometry and spectroscopy covering ∼1.5 years post discovery. SN 2021adxl occurred in an unusual environment, atop a vigorously star-forming region that is offset from its host galaxy core. The appearance of Ly α and O II , as well as the compact core, would classify the host of SN 2021adxl as a “Blueberry” galaxy, analogous to higher redshift, low-metallicity, star-forming dwarf “Green Pea” galaxies. Using several abundance indicators, we find a metallicity of the explosion environment of only ∼0.1 Z ⊙ , the lowest reported metallicity for a Type IIn SN environment. SN 2021adxl reaches a peak magnitude of M r ≈ −20.2 mag and since discovery, SN 2021adxl has faded by only ∼4 magnitudes in the r band with a cumulative radiated energy of ∼1.5 × 10 50 erg over 18 months. SN 2021adxl shows strong signs of interaction with a complex circumstellar medium, seen by the detection of X-rays, revealed by the detection of coronal emission lines, and through multi-component hydrogen and helium profiles. In order to further understand this interaction, we model the H α profile using a Monte Carlo electron scattering code. The blueshifted high-velocity component is consistent with emission from a radially thin spherical shell resulting in the broad emission components due to electron scattering. Using the velocity evolution of this emitting shell, we find that the SN ejecta collide with circumstellar material of at least ∼5 M ⊙ assuming a steady-state mass-loss rate of ∼4 − 6 × 10 −3 M ⊙ yr −1 for the first ∼200 days of evolution. SN 2021adxl was last observed to be slowly declining at ∼0.01 mag d −1 , and if this trend continues, SN 2021adxl will remain observable after its current solar conjunction. Continuing the observations of SN 2021adxl may reveal signatures of dust formation or an infrared excess, similar to that seen for SN 2010jl.

Brennan, S. J. (ORCID:0000000313256235)↗

DOE OSTI2015

An Approach to Enhance pnetCDF Performance in Environmental Modeling Applications

I/O has been considered as a bottleneck in parallel applications. The software package, pnetCDF which works with parallel file systems, was developed to address this issue and provide parallel I/O capability. This study examines the performance of a novel approach which performs data aggregation along either row or column dimension of the spatial domain, and then applies the pnetCDF parallel I/O paradigm. The test was done with three different domain sizes which represents small, moderately large and large data domains, using a small scale Community Multi-scale Air Quality model (CMAQ) mocked up code. The examination includes comparing I/O performance with traditional serial I/O technique, straight application of pnetCDF, and the data aggregation along row and column dimension before applying pnetCDF. After the comparison, “optimal” I/O configurations for this new novel approach were quantified. Data aggregation along the row dimension (pnetCDFcr) works better than along the column dimension (pnetCDFcc) although it may perform slightly worse than straight the pnetCDF method with a small number of processors. When the number of processors becomes larger, pnetCDFcr out performs pnetCDF significantly. If the number of processors keeps increasing, pnetCDF reaches a point that the performance is even worse than the serial I/O technique. This new approach has also been tested on a real application where it performs two times better than the straight pnetCDF paradigm.

I/O↗

DOE OSTI2015

The Acoustic Emission signal acquired by the microphones placed in the CABRI test device along the fourteen last R.I.A. experiments: an example of reproducible research in nuclear science

The CABRI facility is an experimental nuclear reactor of the French Atomic Energy Commission. It is located at the Cadarache Research Centre in southern France and it is designed to act as a support to the French nuclear infrastructure. The purpose of the new testing programme termed, 'CABRI International Programme' (CIP) is to study the behaviour of PWR-type fuel rods at high burnup, equipped with an 'advanced' cladding, under Reactivity Initiated Accident (RIA) conditions (such as the scenario of a control rod ejection). Within the framework of this programme, piloted and funded by the French Institute of Nuclear Radioprotection and Safety (IRSN), ten tests are to be conducted with a frequency of two tests per year. The LPRE laboratory of the CEA which is in charge of the Preparation, realisation and breakdown of the test results studies the possibility to set up a new test analysis based on the processing of signals coming from sensors placed within the test equipment. During the experimental phase, the behaviour of the fuel element generates acoustic waves which can be detected by two microphones placed upstream and downstream of the test device. Studies showed the interest to realize temporal and spectral analyses on these signals by showing the existence of signatures which can be correlated with physical phenomena as the rod failure or the test shutdown (i.e. scram). The work presented in this article results from the will to consolidate these studies. Since the main phenomenon to be tracked is the fuel rod failure, the aim would be to highlight specific events which would have been precursors of the rod failure in order to use in the future these signals for further interpretation. Such an antecedent works resumption leads to a better understanding of the experimental needs and constitutes a good initial state to prepare the new very fast digital data acquisition systems. Although all the raw data are accessible in the form of text files, analyses and graphics representations were not straightforward to reproduce from the ancient studies since that, on one hand, people who were in charge of the original work left the laboratory and on the other hand because it is not easy when the time passes, even with our own work, to be able to remember the steps of data manipulations and the exact setup: - During the ancient experiments the use of analog data acquisition systems required to digitize tapes to be able to realize computer treatments. That had had for consequence to lose the initial dating. This one must be correctly edited to do temporal comparisons. - Analyses require functions for calculations whose parameters has to be well-known to reach the same results. We thus wished to manage our workflow in the idea that it can be easily reproducible on all the experiments. The object of the work presented in this article was to put in practice this strong bind between the data, treatments and generation of the document in order not to hesitate to do the iteration principle in action. We do not have to be afraid by the data driven analyses. According to the philosophy of the literate programming, the text of the technical document is woven with the computer code that produces all the printed output as tables, graphs for the study eliminating hence the unrealistic cut and paste. This difficulty is not specific to the nuclear domain. For many years, researchers have been worked out solutions to this mundane issue. And, presently, new technologies and high-level programming languages offer us actual answers. We will firstly present the tools applied in our laboratory to implement this workflow, then we will describe the global perception carried out to continue the study of the Acoustic Emission signals recorded by the two microphones during the fourteen last CABRI R.I.A. test.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2014

Electron-optical systems for planar gyrotrons

The methodology of designing an electron-optical system (EOS) that forms sheet helical electron beams (HEBs) for high-power gyrotrons is developed. As an example, we consider the EOS for a 140-GHz gyrotron operated at the first harmonic of the cyclotron frequency with an accelerating voltage of 50 kV, a beam current of 30 A, and a magnetic field compression of 36. A planar geometry of the magnetron-injection gun (MIG) is suggested. The adiabatic theory of MIGs modified for the planar geometry of EOS is used for preliminary estimations of MIG parameters. Numerical simulation of the HEB properties based on the CST STUDIO SUITE 3D code is performed to find the optimal configuration of a planar MIG. The accuracy of the calculated data is discussed. The main factors that affect the HEB quality are considered. It is shown that a sheet HEB with a pitch-factor of 1.3 and velocity spread not exceeding 25%–30% can be formed; this is quite acceptable for high-efficiency operation of modern gyrotrons. Calculation of the beam-wave interaction with the obtained HEB parameters proved that a high output power with a sufficiently good efficiency of about 20% can be reached. Simulations show the feasibility of the experimental implementation of a novel planar EOS and its use in short-wave planar gyrotrons. The developed technique can be used for the study and optimization of planar gyrotrons of different frequency bands and power levels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Modeling jet and outflow feedback during star cluster formation

Powerful jets and outflows are launched from the protostellar disks around newborn stars. These outflows carry enough mass and momentum to transform the structure of their parent molecular cloud and to potentially control star formation itself. Despite their importance, we have not been able to fully quantify the impact of jets and outflows during the formation of a star cluster. The main problem lies in limited computing power. We would have to resolve the magnetic jet-launching mechanism close to the protostar and at the same time follow the evolution of a parsec-size cloud for a million years. Current computer power and codes fall orders of magnitude short of achieving this. In order to overcome this problem, we implement a subgrid-scale (SGS) model for launching jets and outflows, which demonstrably converges and reproduces the mass, linear and angular momentum transfer, and the speed of real jets, with ∼1000 times lower resolution than would be required without the SGS model. We apply the new SGS model to turbulent, magnetized star cluster formation and show that jets and outflows (1) eject about one-fourth of their parent molecular clump in high-speed jets, quickly reaching distances of more than a parsec, (2) reduce the star formation rate by about a factor of two, and (3) lead to the formation of ∼1.5 times as many stars compared to the no-outflow case. Most importantly, we find that jets and outflows reduce the average star mass by a factor of ∼ three and may thus be essential for understanding the characteristic mass of the stellar initial mass function.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Pre-Equilibrium De-Excitations for Neutrino-Nucleus Interactions

The Deep Underground Neutrino Experiment (DUNE) is sensitive to MeV-scale energy depositions from low-energy astrophysical neutrinos, including those from core-collapse supernovae. Interpreting these detector signals requires accurate modeling of the nuclear de-excitation that follows from the neutrino-nucleus interaction. The MARLEY (Model of Argon Reaction Low-Energy Yields) event generator specializes in the low-energy regime. MARLEY currently assumes the residual nucleus equilibrates immediately after the primary interaction. This omits the intermediate pre-equilibrium stage in which energy redistributes among nucleons until statistical equilibrium is reached. While pre-equilibrium effects are well established for nucleon-induced reactions, they have not previously been studied for neutrino-nucleus interactions. This work addresses that gap by implementing a two-component exciton model, which is the first dedicated treatment of pre-equilibrium de-excitation for neutrino-nucleus interactions, restructured around MARLEY's existing class hierarchy to prepare for direct integration, including particle-hole state densities, internal transition rates, and pre-equilibrium particle emission. The calculations show encouraging agreement with the TALYS-2.2 nuclear reaction code for neutron-nucleus interactions. We further propose a concrete integration path into the full MARLEY event generator, including derived class structure and an extended event record for pre-equilibrium vertices in support of future reweighting. Remaining work focuses on refining the emission width calculation, adding $\gamma$-ray emission, and completing this integration to quantify the impact of pre-equilibrium effects on the expected low-energy neutrino signals in DUNE and similar experiments.

Visser, Erin [Michigan State U., East Lansing (mai↗

DOE OSTI2024

Integrated modeling of boron powder injection for real-time plasma-facing component conditioning

An integrated modeling framework for investigating the application of solid boron (B) powder injection for real-time surface conditioning of plasma-facing components (PFCs) in tokamak environments is presented. Utilizing the DIII-D impurity powder dropper (IPD) setup, this study simulates B powder injection scenarios ranging from milligrams to tens of milligrams per second, corresponding to boron flux rates of 10 20 –10 21 B/s in standard L-mode conditions. The comprehensive modeling approach combines EMC3-EIRENE for simulating the deuterium plasma background and the Dust Injection Simulator (DIS) for the ablation and transport of the boron powder particles. EMC3 trace impurity fluid modeling results show substantial boron transport to the inboard lower divertor, predominantly influenced by the main ion plasma flow. The dependency on powder particle size (5-250 µm) was found to be insignificant for the scenario considered. The effects of erosion and redeposition were considered to reconcile the discrepancies with experimental observations, which saw substantial deposition on the outer divertor plasma-facing components. For this purpose, the WallDYN3D code was updated to include boron sources within the plasma domain and integrated into the modeling framework. The mixed-material migration modeling shows evolving boron deposition patterns, suggesting the formation of mixed B-C layers or predominantly B coverage depending on the powder mass flow rate. While the modeling outcomes at lower B injection rates tend to align with DIII-D experimental observations, the prediction of near-pure boron layers at higher rates has yet to be experimentally verified in the carbon environment of the DIII-D tokamak. The extensive reach of boron layers found in the modeling suggests the need for modeling that encompasses the entire wall geometry for more accurate experimental correlations. This integrated approach sets a precedent for analyzing and applying real-time in-situ boron coating techniques in advanced tokamak scenarios, potentially extendable to ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

The structural, electronic and phonon behavior of CsPbI{sub 3}: A first principles study

Metal halide perovskites are optoelectronic materials that have attracted enormous attention as solar cells with power conversion efficiencies reaching 20%. The benefit of using hybrid compounds resides in their ability to combine the advantage of these two classes of compounds: the high mobility of inorganic materials and the ease of processing of organic materials. In spite of the growing attention of this new material, very little is known about the electronic and phonon properties of the inorganic part of this compounds. A theoretical study of structural, electronic and phonon properties of metal-halide cubic perovskite, CsPbI{sub 3} is presented, using first-principles calculations with planewave pseudopotential method as personified in PWSCF code. In this approach local density approximation (LDA) is used for exchange-correlation potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗

DOE OSTI2017

Calculations vs. Measurements for Remnant Dose Rates from SNS Spent Structures

The Spallation Neutron Source (SNS) in Oak Ridge, Tennessee, is an accelerator driven neutron scattering facility for materials research. Presently SNS is capable to operate at 1.4 MW proton beam power incident on a mercury target with a proton beam energy of 1 GeV and 60 Hz repetition rate. SNS target system components are periodically replaced because they reach their end-of-life due to radiation induced material damage. Target vessel, which houses mercury target, is exchanged about two-three times per year and the proton beam window (PBW) is exchanged every two – three years.Each spent structure that leaves the SNS site requires supporting documentation with radionuclide inventory and dose rate prediction for the time of the transportation. Neutronics analyses are performed, assuming realistic irradiation history and decay case to ensure that the container/package, housing the structure, is compliant with the waste management regulations. Analyses are complex due to geometry, multi-code usage and following data treatment.To validate analyses, measurements of dose rates from the spent target vessel # 13 and PBW module #5 were performed. Neutronics analyses were performed to calculate residual dose rates from both structures for the time of measurements.

Popova, Irina I.↗

DOE OSTI2014

Transport simulations of linear plasma generators with the B2.5-Eirene and EMC3-Eirene codes

Linear plasma generators are cost effective facilities to simulate divertor plasma conditions of present and future fusion reactors. For this research, the codes B2.5-Eirene and EMC3-Eirene were extensively used for design studies of the planned Material Plasma Exposure eXperiment (MPEX). Effects on the target plasma of the gas fueling and pumping locations, heating power, device length, magnetic configuration and transport model were studied with B2.5-Eirene. Effects of tilted or vertical targets were calculated with EMC3-Eirene and showed that spreading the incident flux over a larger area leads to lower density, higher temperature and off-axis profile peaking in front of the target. In conclusion, the simulations indicate that with sufficient heating power MPEX can reach target plasma conditions that are similar to those expected in the ITER divertor. B2.5-Eirene simulations of the MAGPIE experiment have been carried out in order to establish an additional benchmark with experimental data from a linear device with helicon wave heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

A structural model of the long-term degradation of the concrete biological shield

The concrete biological shield (CBS) of light water reactors is exposed to high neutron radiation dose in the long term, which may lead to the degradation of the concrete’s mechanical properties. Given the important shielding role of the CBS, it is necessary to investigate the irradiation effects at the structural scale and provide estimates of the damage extent from the wall’s inner surface to study potential license renewals. For this purpose, we developed a mechanical model accounting for radiation-induced expansion, creep, and damage in concrete using the Grizzly finite element code, informed by ex-core neutron flux calculations using the VERA tool. The model was applied to a 3D CBS structure represented by the CBS wall, a steel liner, reinforcement bars, and a concrete base mat and evaluated damage at 40, 60, and 80 years of operation. The VERA model predicted a maximum fluence of approximately 2 x 10 19 ncm -2 at 80 years of operation. The results showed that damage is highest at the inner surface of the CBS wall and gradually decreases with depth. It extends beyond the rebar after 60 years and reaches a depth of approximately 12 cm at 80 years.

42 ENGINEERING↗

DOE OSTI2014

The coyote universe extended: Precision emulation of the matter power spectrum

Modern sky surveys are returning precision measurements of cosmological statistics such as weak lensing shear correlations, the distribution of galaxies, and cluster abundance. To fully exploit these observations, theorists must provide predictions that are at least as accurate as the measurements, as well as robust estimates of systematic errors that are inherent to the modeling process. In the nonlinear regime of structure formation, this challenge can only be overcome by developing a large-scale, multi-physics simulation capability covering a range of cosmological models and astrophysical processes. As a first step to achieving this goal, we have recently developed a prediction scheme for the matter power spectrum (a so-called emulator), accurate at the 1% level out to k ∼ 1 Mpc{sup –1} and z = 1 for wCDM cosmologies based on a set of high-accuracy N-body simulations. It is highly desirable to increase the range in both redshift and wavenumber and to extend the reach in cosmological parameter space. To make progress in this direction, while minimizing computational cost, we present a strategy that maximally reuses the original simulations. We demonstrate improvement over the original spatial dynamic range by an order of magnitude, reaching k ∼ 10 h Mpc{sup –1}, a four-fold increase in redshift coverage, to z = 4, and now include the Hubble parameter as a new independent variable. To further the range in k and z, a new set of nested simulations run at modest cost is added to the original set. The extension in h is performed by including perturbation theory results within a multi-scale procedure for building the emulator. This economical methodology still gives excellent error control, ∼5% near the edges of the domain of applicability of the emulator. A public domain code for the new emulator is released as part of the work presented in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

The properties and fate of the galactic center G2 cloud

The object G2 was recently discovered descending into the gravitational potential of the supermassive black hole (BH) Sgr A*. We test the photoionized cloud scenario, determine the cloud properties, and estimate the emission during the pericenter passage. The incident radiation is computed starting from the individual stars at the locations of G2. The radiative transfer calculations are conducted with CLOUDY code and 2011 broadband and line luminosities are fitted. The spherically symmetric, tidally distorted, and magnetically arrested cloud shapes are tested with both the interstellar medium dust and 10 nm graphite dust. The best-fitting magnetically arrested model has the initial density n {sub init} = 1.8 × 10{sup 5} cm{sup –3}, initial radius R {sub init} = 2.2 × 10{sup 15} cm = 17 mas, mass m {sub cloud} = 4 M {sub Earth}, and dust relative abundance A = 0.072. It provides a good fit to 2011 data, is consistent with the luminosities in 2004 and 2008, and reaches an agreement with the observed size. We revise down the predicted radio and X-ray bow shock luminosities to be below the quiescent level of Sgr A*, which readily leads to non-detection in agreement to observations. The magnetic energy dissipation in the cloud at the pericenter coupled with more powerful irradiation may lead to an infrared source with an apparent magnitude m{sub L{sup ′}}≈13.0. No shock into the cloud and no X-rays are expected from cloud squeezing by the ambient gas pressure. Larger than previously estimated cloud mass m {sub cloud} = (4-20) M {sub Earth} may produce a higher accretion rate and a brighter state of Sgr A* as the debris descend onto the BH.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Use of Hygrothermal Models for Understanding Water Transport in Corrosion Under Insulation Applications

Thermal insulation is used in operating facilities to conserve heat and protect against freezing amongst others. A consequence of insulating the equipment is the necessity to manage the introduced threat of corrosion under insulation (CUI). For CUI to occur, water and oxygen must enter and migrate through the insulation to reach the external surface of the equipment. The water transport characteristics are dependent on several factors such as type of insulation, type of jacketing, equipment operating temperature, external weather, water entry/leakage rate and cyclic service. In hot piping, there are competing water transport characteristics, as in water vapor moves outwards away from the equipment as water enters insulation. Knowing the water transport and the parameters that influence the time of wetness at the metal surface helps in understanding conditions favoring CUI.The use of transient hygrothermal models for moisture control is well established in the building insulation design codes and standards. The building designs naturally shed the liquid water to minimize entry and facilitate vapor management so that moisture doesn’t accumulate within building. Several building industry hygrothermal models have been developed and are available for commercial use. One such commercial model has been used to understand water transport in a CUI application. The case study involves evaluation of piping and equipment installed with a closed cell polyurethane insulation. The hygrothermal model provided insights on the parameters influencing the time of wetness and the ease of water escaping the equipment-insulation-jacketing system.

Venkateswaran, Sai P.↗

DOE OSTI2013

Modeling of plasma processes in the slowly diverging magnetic fields at the exit of an applied-field magnetoplasmadynamic thruster

The performance of plasma thrusters with applied electric and magnetic fields can be enhanced by increasing the magnetic field strength, which is applied in the thrust chamber and the exit region propulsive plume. The ejected plasma which passes through a slowly diverging magnetic field will expand but can be restricted within the magnetic nozzle fields. To examine in detail the processes that occur, a new method with Particle-in-cell calculations is applied here. A two-dimensional axisymmetric particle dynamic code is used to model an AF-MPDT (Applied-field Magnetoplasmadynamic Thruster) for which extensive experimental data are available; it used Ar propellant and had applied magnetic coils of 101.5 mm radius and 153 mm length. From the results of the simulation study, it is found that total thrust increases linearly with magnetic field strength in the range of 0–0.1 T, but it decreases with increasing applied magnetic field up to 0.6 T. Thrust efficiency is found to increase to a maximum of 8.4% when B = 0.1 T; further, the peak value of nozzle efficiency reaches 91% at a moderate magnetic field (0.3 T). In detail, it is found that distributions of plasma density (10{sup 14}–10{sup 15} m{sup −3}) that form in the magnetic nozzle demonstrate a significant pattern of concentration up to fields of B = 0.3 T where ions begin to be magnetized. However, azimuthal velocities of ions behave differently with different degrees of magnetization, i.e., weakly magnetized ions follow rotating electrons in a right-handed direction, while fully magnetized ions revolve in left-handed direction due to electromagnetic forces. Notably, a feedback effect on total magnetic field due to plasma motion identified in other studies is not found to be present in the working conditions of the AF-MPDT examined here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

The magnetic furnace: intense core dynamos in B stars

The dynamo action achieved in the convective cores of main-sequence massive stars is explored here through three-dimensional (3D) global simulations of convective core dynamos operating within a young 10 M{sub ⊙} B-type star, using the anelastic spherical harmonic code. These simulations capture the inner 65% of this star by radius, encompassing the convective nuclear-burning core (about 23% by radius) and a portion of the overlying radiative envelope. Eight rotation rates are considered, ranging from 0.05% to 16% of the surface breakup velocity, thereby capturing both convection that barely senses the effects of rotation and other situations in which the Coriolis forces are prominent. The vigorous dynamo action realized within all of these turbulent convective cores builds magnetic fields with peak strengths exceeding a megagauss, with the overall magnetic energy (ME) in the faster rotators reaching super-equipartition levels compared to the convective kinetic energy (KE). The core convection typically involves turbulent columnar velocity structures roughly aligned with the rotation axis, with magnetic fields threading through these rolls and possessing complex linkages throughout the core. The very strong fields are able to coexist with the flows without quenching them through Lorentz forces. The velocity and magnetic fields achieve such a state by being nearly co-aligned, and with peak magnetic islands being somewhat displaced from the fastest flows as the intricate evolution proceeds. As the rotation rate is increased, the primary force balance shifts from nonlinear advection balancing Lorentz forces to a magnetostrophic balance between Coriolis and Lorentz forces.

79 ASTRONOMY AND ASTROPHYSICS↗