Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Magnetic field modeling”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Particle-based modelling of axisymmetric tandem mirror devices

In this work, we describe the use of a 1D-2V quasi-neutral hybrid electrostatic PIC with Monte-Carlo Coulomb collisions and non-uniform magnetic field to model the parallel transport and confinement in an axisymmetric tandem mirror device. End-plugs, based on simple-mirrors, are positioned at each end of the device and fueled with neutral beams (25 and 100 keV) to produce a sloshing ion population and increase the density of the end-plugs relative to the central cell. Results show the formation of a potential difference barrier between the central cell and the end-plugs. This potential confines a large fraction of the low energy thermal ions in the central cell which would otherwise be lost in a simple mirror, demonstrating the advantage of the beam-driven tandem mirror configuration relative to simple mirrors. In addition, we explore the effect of end-plug electron temperature on the confinement time of the device and compare it with theoretical estimates. Finally, we discuss the limitations of the code in its present form and describe the next logical steps to improve its predictive capability such as a fully nonlinear Fokker–Planck collision operator, multiply nested flux surface solutions and modeling the exhaust region up to the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Calculations of synchrotron emission from the terrestrial radiation belts

A theoretical model was developed to allow for the calculation of the synchrotron emission arising from high energy electrons trapped in the Van Allen belts of a planet with a dipole magnetic field. The model is general enough to allow for the calculation of the intensity of radiation received by an observer at any distance from and any latitude about the planet. The model is used to compute the emission from the earth's Van Allen belts that one should expect at various latitudes at a distance of 1.92 earth radii, the position of the Radio Astronomy Explorer satellite that was launched in 1968, for the frequencies 1.3 MHz and 2.2 MHz.

Johnson, R.↗
DOE OSTI2024

Three-dimensional modeling of a plasma in a strong azimuthal magnetic field

Three-dimensional magnetohydrodynamic (MHD) simulations are able to model the generation of disc-shaped plasma, driven by laser ablation from a current carrying rod in a pulsed-power machine producing azimuthal magnetic fields of 2–3 MG. The plasma at such extreme conditions is unique in that the parameter space for the plasma β and Hall parameter χ transition from below unity to greater than unity at different stages of the plasma generation. In simulations, the formation of the plasma disc in the azimuthal direction, is driven by heat flux from the laser spot and depends on the set of transport coefficients used in simulations. The most recent set of transport coefficients leads to the formation of plasma ejecta at the back end of the rod, that qualitatively matches experiments. Specifically, the cross-gradient Nernst effect, which twists magnetic field, is shown to have a large effect on the shape of the back end ejecta. In the direction along the axis of the rod there is propagation of perturbations from the disc as observed in experiments. In simulations the period of temperature perturbations is in good agreement with experimental results. Furthermore, an instability due to coupling of heat flux and the magnetic field advection provides a possible explanation for perturbation growth along the axis of the rod, and the instability growth rate is consistent with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Heating of solar coronal holes by reflected Alfven waves

As a continuation of the work of Moore et al. (1991), who found evidence that coronal holes are heated by Alfven waves that are reflected back down within the coronal holes, this paper shows that to demonstrate this evidence, it is only necessary to consider a subset of the Moore et al. models, namely, those having radial magnetic field. Using these models, it is shown that the Alfven velocity is not constant in the atmosphere of coronal holes, but changes with height (or radius), causing downward reflection of all upward Alfven waves of sufficiently long wavelength (or period).

Moore, R. L.↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

An Analytical Model for Non-Uniform Magnetic Field Effects on Two-Dimensional Laminar Jet Diffusion Flames

In 1846, Michael Faraday found that permanent magnets could cause candle flames to deform into equatorial disks. He believed that the change in flame shape was caused by the presence of charged particles within the flames interacting with the magnetic fields. Later researchers found that the interaction between the flame ions and the magnetic fields were much too small to cause the flame deflection. Through a force analysis, von Engel and Cozens showed that the change in the flame shape could be attributed to the diamagnetic flame gases in the paramagnetic atmosphere. Paramagnetism occurs in materials composed of atoms with permanent magnetic dipole moments. In the presence of magnetic field gradients, the atoms align with the magnetic field and are drawn into the direction of increasing magnetic field. Diamagnetism occurs when atoms have no net magnetic dipole moment. In the presence of magnetic gradient fields, diamagnetic substances are repelled towards areas of decreasing magnetism. Oxygen is an example of a paramagnetic substance. Nitrogen, carbon monoxide and dioxide, and most hydrocarbon fuels are examples of diamagnetic substances. In order to evaluate the usefulness of these magnets in altering flame behavior, a study has been undertaken to develop an analytical model to describe the change in the flame length of a laminar diffusion jet in the presence of a nonuniform magnetic field.

Calvert, M. E.↗
NASA NTRS1995

Coronal roots of solar wind streams: 3-D MHD modeling

Weak (discontinuous) solutions of the 3-D MHD equations look like a promising tool to model the transonic solar wind with structural elements: current sheets, coronal plumes etc. Using the observational information about various coronal emissions one can include these structural elements into the 3-D MHD solar wind model by embedding the discontinuities of given type. Such 3-D MHD structured solar wind is calculated self-consistently: variants are examined via numerical experiments. In particular, the behavior of coronal plumes in the transonic solar wind flow, is modeled. The input information for numerical modeling (for example, the magnetic field map at the very base of the solar corona) can be adjusted so that fast stream arises over the center of the coronal hole, over the coronal hole boundaries and, even, over the region with closed magnetic topology. 3-D MHD equations have the analytical solution which can serve as a model of supersonic trans-alfvenic solar wind in the (5-20) solar radii heliocentric distance interval. The transverse, nonradial total (gas + magnetic field) pressure balance in the flow is the corner-stone of this solution. The solution describes the filamentation (ray-like structure of the solar corona) and streaming (formation of high-speed streams with velocities up to 800 km/sec) as a consequence of the magnetic field spatial inhomogeneous structure and trans-alfvenic character of the flow. The magnetic field works in the model as a 'controller' for the solar wind streaming and filamentation.

Pisanko, Yu. V.↗
NASA NTRS1988

Modelling a C-type flare observed in microwaves and hard X-rays

Using the very large array (VLA) at 6 and 20 cm wavelength and the hard X-ray burst spectrometer on the Solar Maximum Mission, a two-ribbon flare was observed from the onset phase through the maximum and decline on November 14, 1981. Because of the extensive size of the microwave source and the gradual variations in hard X-rays whose spectrum becomes progressively flatter with time, the flare is classified as a C-type flare. Considering the hardening of the X-ray spectrum and its non-impulsive nature, a coronal trap model was invoked for the energetic electrons. The microwave emission is easily accounted for by gyrosynchronous radiation from mildly relativistic electrons. It was found that the source must be optically thick at 20 cm during the maximum phase, but as the source evolved toward an optically thin regime, the intensity decreased while the degree of circular polarization increased. In an initial homogeneous model, we found that the computed microwave spectrum was too narrow to match the patrol spectrum from 606 to 15400 MHz. In the model, the magnetic field consists of a dipolar arcade bridging the H alpha ribbons, and extending to heights of order 40,000 to 50,000 km. The variation of the magnetic field strength from footpoints to apex causes the gyrosynchrotron spectrum to be broader. Preliminary conclusions regarding the electron distributions producing the hard X-rays and the microwaves, and the suitability of this model for C-type flares is presented.

Schmahl, E. J.↗
NASA NTRS1974

The interior structure of Jupiter (consequences of Pioneer 10 data)

Models of the Jovian interiors are based on theoretical equations of state of hydrogen and helium supported by a few experimental points and an observed parameter such as oblateness, gravitational coefficients, heat emission, and magnetic fields. The models fall into three categories: (1) those which assume a uniform and rather low H2/He ratio throughout the planet, (2) those in which this ratio is solar and thus higher and (3) those which take into account the lack of complete miscibility of the two elements in the condensed state. Recent values of the observed parameters obtained by Pioneer 10 permit improvements of the first two models but also pose new questions. In the first category of models the new data indicate that the amount of hydrogen has to be increased, while in the solar models which have a heavy core (made of SiO2, MgO, Fe and Ni), the abundance of hydrogen has to be decreased, both changes pointing in the direction of incomplete miscibility present in the third category of models.

Smoluchowski, R.↗
DOE OSTI2025

Magnetohydrodynamic simulations preliminarily predict the habitability and radio emission of TRAPPIST-1e

TRAPPIST-1e, an Earth-sized exoplanet in the habitable zone of the nearby M dwarf TRAPPIST-1, may experience magnetospheric responses that vary with stellar space weather, which could potentially influence both its habitability and radio emissions. Our objective is to investigate how different Earth-like magnetospheric configurations of TRAPPIST-1e – specifically variations in dipolar magnetic field strength and axial tilt – respond to diverse stellar space weather conditions, including events analogous to coronal mass ejections (CMEs), and to assess their implications for potential habitability and expected radio emissions. We conducted 3D magnetohydrodynamic simulations of the TRAPPIST-1e system using the PLUTO code in spherical coordinates. The planetary magnetic field was modelled as dipolar, with equatorial strengths from Earth-like to several times stronger. The dipole axis spans a representative range of axial tilts. We investigate four stellar wind environments, from sub-Alfvénic flow to CME-like disturbances. Planetary shielding was quantified based on the magnetopause standoff distance, and radio powers were estimated via empirical scaling laws. Our simulations show that both shielding and radio power depend strongly on the magnetic configuration. Stronger fields increase protection, while larger tilts reduce it. Radio power increases with both field strength and tilt across all wind regimes. An Earth-like magnetic field can provide effective shielding even under intense CMEs, whereas high tilts require stronger fields. Predicted radio powers reach ~10 20 erg s −1 during CMEs, making bursts from close-in, magnetised planets more detectable. However, for TRAPPIST-1e, the maximum cyclotron frequency lies below the Earth’s ionospheric cutoff (~10 MHz), making ground-based detection currently infeasible.

Sun: magnetic fields↗
DOE OSTI2020

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Pair Production near Threshold in Pulsar Magnetic Fields

In pulsar polar cap models, curvature radiation gamma rays produce e(+)e(-) pairs in the strong magnetic fields near the surface of the neutron star. While these gamma rays have energies E sub gamma mc(2), they also propagate at very small angles to the fields, such that the threshold condition, E gamma 2mc(2)/sin theta is just barely satisfied when they pair produce. Threshold effects on the pair production attenuation coefficient, which are due to the discreteness of the e(+)e(-) Landau states, must therefore be considered when computing the mean free paths of curvature radiation photons in pulsar magnetic fields. These effects, which are not incorporated in the asymptotic expression for the attenuation coefficient, have some interesting consequences for pulsar models. Since pair production is suppressed near threshold, the photon mean free paths are longer than previously thought. In magnetic fields approximately 6x10 approximately G, the pairs tend to be produced in the ground state Landau level and will not synchrotron radiate. Since synchrotron radiation is an essential ingredient in the electromagnetic cascades which produce low energy pairs above the acceleration region, pulsars with very high magnetic fields may not produce many pairs.

Harding, A. K.↗
DOE OSTI2026

Early Solar Wind and Dynamo Magnetic Field Topology Predictions for (16) Psyche and Other Asteroids

Abstract Asteroid (16) Psyche is a metal‐rich body that might record an ancient coherent magnetization if some relict crust or mantle is preserved. Herein, we use magnetohydrodynamic simulations to predict (16) Psyche's field topology for several distinct pathways: (i) an early solar wind‐induced magnetization imparted after a larger body was impacted, forming the present‐day asteroid, (ii) a core dynamo magnetization imparted in an asteroid that is either presently largely intact or was a rubble pile, and (iii) magnetization in the turbulent solar nebula disk. For pathway (i) we find the field to be predominantly dipolar and spin axis‐aligned. For pathway (ii) we find the field to be either dipolar and spin axis‐misaligned, or highly multipolar. We also find that (iii) a field produced earlier before the solar nebula cleared, would be highly multipolar. In cases (i) and (ii) we also place constraints on the field strength. Simple detection of a magnetic field without constraining its topology and temporal variability would be insufficient to confirm a remanent source, due to the influence of the present‐day solar wind, electromagnetic induction, and (16) Psyche's high obliquity. For sufficiently strong fields however, the field topology and orientation may reveal key observable consequences of the nature and history of (16) Psyche. Our framework is also broadly applicable to the study of magnetic fields from other asteroids.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

The geomagnetic main field and the geodynamo

Information available on the geomagnetic main field and the geodynamo is presented. Attention is given to the process of mapping the magnetic field, the last version of International Geomagnetic Reference Field Model, and maps of the magnetic field at the core-mantle boundary and their interpretation. Particular consideration is given to the existing geodynamo theories, with special relation given to the Braginsky and Meytlis theory of core turbulence in which the turbulence differs fundamentally from classical turbulence of Kolmogorov type.

Bloxham, Jeremy↗
NASA NTRS1980

An equivalent source model of the satellite-altitude magnetic anomaly field over Australia

The low-amplitude, long-wavelength magnetic anomaly field measured between 400 and 700 km elevation over Australia by the POGO satellites is modeled by means of the equivalent source technique. Magnetic dipole moments are computed for a latitude-longitude array of dipole sources on the earth's surface such that the dipoles collectively give rise to a field which makes a least squares best fit to that observed. The distribution of magnetic moments is converted to a model of apparent magnetization contrast in a layer of constant (40 km) thickness, which contains information equivalent to the lateral variation in the vertical integral of magnetization down to the Curie isotherm and can be transformed to a model of variable thickness magnetization. It is noted that the closest equivalent source spacing giving a stable solution is about 2.5 deg, corresponding to about half the mean data elevation, and that the magnetization distribution correlates well with some of the principle tectonic elements of Australia.

Mayhew, M. A.↗
DOE OSTI2021

Magnetospheres of Terrestrial Exoplanets and Exomoons: Implications for Habitability and Detection

Characterizing habitable exoplanets and/or their moons is of paramount importance. Here we show the results of our magnetic field topological modeling, which demonstrate that terrestrial exoplanet–exomoon coupled magnetospheres work together to protect the early atmospheres of both the exoplanet and the exomoon. When exomoon magnetospheres are within the exoplanet's magnetospheric cavity, the exomoon magnetosphere acts like a protective magnetic bubble providing an additional magnetopause confronting the stellar winds when the moon is on the dayside. In addition, magnetic reconnection would create a critical pathway for the atmosphere exchange between the early exoplanet and exomoon. When the exomoon's magnetosphere is outside of the exoplanet's magnetosphere it then becomes the first line of defense against strong stellar winds, reducing the exoplanet's atmospheric loss to space. A brief discussion is given on how this type of exomoon would modify radio emissions from magnetized exoplanets.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

Predicting Hall Thruster Operational Lifetime

A simple analytic model predicted Hall thruster channel erosion based on thruster geometry, operating conditions, and magnetic field configuration. This model relied on a one-dimensional representation of the plasma with a fixed ionization fraction and variable ion energies based on the magnetic field distribution. Sputtering was modeled as the result of elastic scattering of ions by neutrals within the channel. Not all scattered ions and neutrals were assumed to reach the channel walls as a result of additional subsequent scattering events. Incorporating this phenomenon resulted in a greater predicted decrease in erosion rate with time than predicted based only on geometric effects. Results from this model were compared to SPT 100 experimental erosion data.

Manzella, David↗
NASA NTRS1989

Plasma pressure in the environment of Jupiter, inferred from Voyager 1 magnetometer observations

A model combining the internal magnetic field with a self-consistent model of the Jovian magnetodisc was fitted to the Voyager 1 Jovian magnetic field data by means of a generalized inverse technique. The model parameters included the internal field spherical harmonic coefficients as well as with parameters describing the plasma distribution in the magnetosphere. Assuming that the pressure in the middle and outer magnetosphere is related to the unit flux tube volume V through PV exp gamma = const, the model fit yielded a value of 0.88 for gamma. If the hot (30 keV) plasma is transported adiabatically inward under the interchange instability triggered by centrifugal force of the heavy torus ions, losses are not sufficient to account for such a low value of gamma beyond L = 10. Closer to the planet, as the outer edge of the Io plasma torus is approached, PV exp gamma is found to decrease inward, as expected from the particle measurements, which identified an inner boundary of the particle fluxes in that region. With the present stage of the development of magnetodisc models, secular variations of the internal field still remain difficult to estimate.

Caudal, G.↗