S.R. Hudson; D.A. Monticello; A.H. Reiman; D.J. Strickler; S.P. Hirshman; L-P. Ku; E. Lazarus; A. Brooks; M.C. Zarnstorff; A.H. Boozer; G-Y. Fu; and G.H. Neilson
2003-09-15
For the (non-axisymmetric) stellarator class of plasma confinement devices to be feasible candidates for fusion power stations it is essential that, to a good approximation, the magnetic field lines lie on nested flux surfaces; however, the inherent lack of a continuous symmetry implies that magnetic islands responsible for breaking the smooth topology of the flux surfaces are guaranteed to exist. Thus, the suppression of magnetic islands is a critical issue for stellarator design, particularly for small aspect ratio devices. Pfirsch-Schluter currents, diamagnetic currents, and resonant coil fields contribute to the formation of magnetic islands, and the challenge is to design the plasma and coils such that these effects cancel. Magnetic islands in free-boundary high-pressure full-current stellarator magnetohydrodynamic equilibria are suppressed using a procedure based on the Princeton Iterative Equilibrium Solver [Reiman and Greenside, Comp. Phys. Comm. 43 (1986) 157] which iterate s the equilibrium equations to obtain the plasma equilibrium. At each iteration, changes to a Fourier representation of the coil geometry are made to cancel resonant fields produced by the plasma. The changes are constrained to preserve certain measures of engineering acceptability and to preserve the stability of ideal kink modes. As the iterations continue, the coil geometry and the plasma simultaneously converge to an equilibrium in which the island content is negligible, the plasma is stable to ideal kink modes, and the coils satisfy engineering constraints. The method is applied to a candidate plasma and coil design for the National Compact Stellarator Experiment [Reiman, et al., Phys. Plasmas 8 (May 2001) 2083].
Coope, Wilfred Anthony; Ferrando i Margalet, Sergi; Allfrey, Simon J.; Kisslinger, Johann; Wobig, Horst F.G.; Narushima, Yoshiro; Okamura, Shoichi; Suzuki, Chihiro; Watanabe, Kiyomasa Y.; Yamazaki, Kozo; Isaev, Maxim Yu
2004-09-15
The impact of the bootstrap current is investigated on the equilibrium properties of a two-period quasi-axisymmetric stellarator reactor with free boundary and on the corresponding ideal magnetohydrodynamic stability properties. Although the magnetic field strength B spectrum is dominated by a m/n = 1/0 component, the discrete filamentary coils trigger some small-amplitude symmetry-breaking components that can disturb the quasi-symmetry of B. Finite {beta} causes the plasma column to shift outward in the absence of bootstrap current. With a self-consistent bootstrap current in the 1/{nu} regime, the plasma becomes more elongated and more distorted in the horizontally elongated up-down symmetric cross section. At {beta} (approximately equal to) 3.25%, the plasma can be restored to its near-vacuum shape with the application of a vertical field with coil currents 20% of those of the modular coils, but at the expense of a significant mirror component in the B-field spectrum. The bootstrap current causes the rotational transform {iota} profile to increase above the critical resonant value ({iota}{sub c} = 1/2 for {beta} {>=}1.1%) and combines with the Pfirsch-Schlueter current to destabilize a m/n = 2/1 external kink mode for {beta} {>=}1.8%.
S.R. Hudson; D.A. Monticello; A.H. Reiman; A.H. Boozer; D.J. Strickler; S.P. Hirshman; M.C. Zarnstorff
2002-11-19
Magnetic islands in free-boundary stellarator equilibria are suppressed using a procedure that iterates the plasma equilibrium equations and, at each iteration, adjusts the coil geometry to cancel resonant fields produced by the plasma. The coils are constrained to satisfy certain measures of engineering acceptability and the plasma is constrained to ensure kink stability. As the iterations continue, the coil geometry and the plasma simultaneously converge to an equilibrium in which the island content is negligible. The method is applied with success to a candidate plasma and coil design for the National Compact Stellarator eXperiment [Physics of Plasma, 7 (2000) 1911].
Free-boundary magnetohydrodynamic equilibria with flow
Schmitt, R. F.; Park, G. Y.; Guazzotto, L.; Strauss, H.; Chang, C.-S.
2011-02-15
The finite-element M3D code [W. Park et al., Phys. Plasmas 6, 1796 (1999)] has been modified to include a free-boundary equilibrium solver with arbitrary toroidal and poloidal flows. With this modification, the M3D code now has the capability to self-consistently model two essential ingredients necessary for equilibrium calculations in the edge region, namely, free-boundary and arbitrary flow. As a free-boundary code, M3D includes the separatrix and scrape-off layer regions in the equilibrium calculation. Poloidal flows in the subsonic, supersonic, and transonic regimes can be calculated with the new version of the M3D code. Calculation results show that the presence of equilibrium flows, in particular those next to the plasma boundary, can considerably influence the position of the X-point and magnetic separatrix shape/location and hence the position of the strike point on the divertor plates. Moreover, it is shown that poloidal flow is not a rigid-body rotation, with the fastest flows occurring on the inboard side of the plasma. A numerical confirmation of the ''de Laval nozzle'' model of Betti and Freidberg [R. Betti and J. P. Freidberg, Phys. Plasmas 7, 2439 (2000)] for free-boundary equilibrium calculations is obtained, with the formation of the predicted discontinuities between regions of subsonic and supersonic flows (with respect to the poloidal sound speed). Finally, a detailed comparison between isentropic and isothermal equilibria is presented, showing qualitative analogies and quantitative differences.
Free-boundary Full-pressure Island Healing in a Stellarator: Coil-healing
S.R.Hudson; A. Reiman; D. Strickler; A. Brooks; D.A. Monticello; and S.P. Hirshman
2002-05-07
The lack of axisymmetry in stellarators guarantees that in general magnetic islands and chaotic magnetic field lines will exist. As particle transport is strongly tied to the magnetic field lines, magnetic islands and chaotic field lines result in poor plasma confinement. For stellarators to be feasible candidates for fusion power stations it is essential that, to a good approximation, the magnetic field lines lie on nested flux-surfaces, and the suppression of magnetic islands is a critical issue for stellarator coil design, particularly for small aspect ratio devices. A procedure for modifying stellarator coil designs to eliminate magnetic islands in free-boundary full-pressure magnetohydrodynamic equilibria is presented. Islands may be removed from coil-plasma free-boundary equilibria by making small changes to the coil geometry and also by variation of trim coil currents. A plasma and coil design relevant to the National Compact Stellarator Experiment is used to illustrate the technique.
A Priori Estimates for Free Boundary Problem of Incompressible Inviscid Magnetohydrodynamic Flows
NASA Astrophysics Data System (ADS)
Hao, Chengchun; Luo, Tao
2014-06-01
In the present paper, we prove the a priori estimates of Sobolev norms for a free boundary problem of the incompressible inviscid magnetohydrodynamics equations in all physical spatial dimensions n = 2 and 3 by adopting a geometrical point of view used in Christodoulou and Lindblad (Commun Pure Appl Math 53:1536-1602, 2000), and estimating quantities such as the second fundamental form and the velocity of the free surface. We identify the well-posedness condition that the outer normal derivative of the total pressure including the fluid and magnetic pressures is negative on the free boundary, which is similar to the physical condition (Taylor sign condition) for the incompressible Euler equations of fluids.
Global magnetohydrodynamic instabilities in the L-2M stellarator
Mikhailov, M. I.; Shchepetov, S. V.; Nührenberg, C.; Nührenberg, J.
2015-12-15
Analysis of global magnetohydrodynamic (MHD) instabilities in the L-2M stellarator (Prokhorov General Physics Institute, Russian Academy of Sciences) is presented. The properties of free-boundary equilibria states are outlined, the stability conditions for small-scale modes are briefly discussed, and the number of trapped particles is estimated. All the magnetic configurations under study are stable against ballooning modes. It is shown that global ideal internal MHD modes can be found reliably only in Mercier unstable plasmas. In plasma that is stable with respect to the Mercier criterion, global unstable modes that are localized in the vicinity of the free plasma boundary and are not associated with any rational magnetic surface inside the plasma (the so-called peeling modes) can be found. The radial structure of all perturbations under study is almost entirely determined by the poloidal coupling of harmonics. The results of calculations are compared with the available experimental data.
Dislocations in magnetohydrodynamic waves in a stellar atmosphere.
López Ariste, A; Collados, M; Khomenko, E
2013-08-23
We describe the presence of wave front dislocations in magnetohydrodynamic waves in stratified stellar atmospheres. Scalar dislocations such as edges and vortices can appear in Alfvén waves, as well as in general magnetoacoustic waves. We detect those dislocations in observations of magnetohydrodynamic waves in sunspots in the solar chromosphere. Through the measured charge of all the dislocations observed, we can give for the first time estimates of the modal contribution in the waves propagating along magnetic fields in solar sunspots. PMID:24010425
Free boundary skin current MHD (magnetohydrodynamic) equilibria
Reusch, M.F.
1988-02-01
Function theoretic methods in the complex plane are used to develop simple parametric hodograph formulae which generate sharp boundary equilibria of arbitrary shape. The related method of Gorenflo and Merkel is discussed. A numerical technique for the construction of solutions, based on one of the methods is presented. A study is made of the bifurcations of an equilibrium of general form. 28 refs., 9 figs.
Free-boundary ideal MHD stability of W7-X divertor equilibria
NASA Astrophysics Data System (ADS)
Nührenberg, C.
2016-07-01
Plasma configurations describing the stellarator experiment Wendelstein 7-X (W7-X) are computationally established taking into account the geometry of the test-divertor unit and the high-heat-flux divertor which will be installed in the vacuum chamber of the device (Gasparotto et al 2014 Fusion Eng. Des. 89 2121). These plasma equilibria are computationally studied for their global ideal magnetohydrodynamic (MHD) stability properties. Results from the ideal MHD stability code cas3d (Nührenberg 1996 Phys. Plasmas 3 2401), stability limits, spatial structures and growth rates are presented for free-boundary perturbations. The work focusses on the exploration of MHD unstable regions of the W7-X configuration space, thereby providing information for future experiments in W7-X aiming at an assessment of the role of ideal MHD in stellarator confinement.
Free boundary ballooning mode representation
Zheng, L. J.
2012-10-15
A new type of ballooning mode invariance is found in this paper. Application of this invariance is shown to be able to reduce the two-dimensional problem of free boundary high n modes, such as the peeling-ballooning modes, to a one-dimensional problem. Here, n is toroidal mode number. In contrast to the conventional ballooning representation, which requires the translational invariance of the Fourier components of the perturbations, the new invariance reflects that the independent solutions of the high n mode equations are translationally invariant from one radial interval surrounding a single singular surface to the other intervals. The conventional ballooning mode invariance breaks down at the vicinity of plasma edge, since the Fourier components with rational surfaces in vacuum region are completely different from those with rational surfaces in plasma region. But, the new type of invariance remains valid. This overcomes the limitation of the conventional ballooning mode representation for studying free boundary modes.
Free boundaries in problems with hysteresis
Apushkinskaya, D. E.; Uraltseva, N. N.
2015-01-01
Here, we present a survey concerning parabolic free boundary problems involving a discontinuous hysteresis operator. Such problems describe biological and chemical processes ‘with memory’ in which various substances interact according to hysteresis law. Our main objective is to discuss the structure of the free boundaries and the properties of the so-called ‘strong solutions’ belonging to the anisotropic Sobolev class with sufficiently large q. Several open problems in this direction are proposed as well. PMID:26261368
NASA Astrophysics Data System (ADS)
Nicolas, T.; Ichiguchi, K.
2016-02-01
The results of a reduced magnetohydrodynamic model for the interchange mode in a stellarator are reported. The model is based on Strauss equations, with the addition of resistivity, viscosity and perpendicular heat conductivity on the one hand, and ion and electron diamagnetic effects on the other hand. The ideally unstable and stable mode growth rate and frequency are studied. The main result is that the ideally unstable mode can rotate in the electron direction even in the absence of electron diamagnetic effects, due to a linear bifurcation caused by the heat conductivity. The presence of the bifurcation is explained analytically. The diamagnetic effects can be destabilizing. For the ideally stable (resistive) mode, we find that the presence of viscosity and heat conductivity causes rotation in the electron direction in the case where {{T}e}={{T}i} .
RADIATION MAGNETOHYDRODYNAMICS SIMULATION OF PROTO-STELLAR COLLAPSE: TWO-COMPONENT MOLECULAR OUTFLOW
Tomida, Kengo; Tomisaka, Kohji; Ohsuga, Ken; Matsumoto, Tomoaki; Machida, Masahiro N.; Saigo, Kazuya E-mail: tomisaka@th.nao.ac.jp E-mail: masahiro.machida@nao.ac.jp E-mail: matsu@hosei.ac.jp
2010-05-01
We perform a three-dimensional nested-grid radiation magnetohydrodynamics (RMHD) simulation with self-gravity to study the early phase of the low-mass star formation process from a rotating molecular cloud core to a first adiabatic core just before the second collapse begins. Radiation transfer is implemented with the flux-limited diffusion approximation, operator-splitting, and implicit time integrator. In the RMHD simulation, the outer region of the first core attains a higher entropy and its size is larger than that in the magnetohydrodynamics simulations with the barotropic approximation. Bipolar molecular outflow consisting of two components is driven by magnetic Lorentz force via different mechanisms, and shock heating by the outflow is observed. Using the RMHD simulation we can predict and interpret the observed properties of star-forming clouds, first cores, and outflows with millimeter/submillimeter radio interferometers, especially the Atacama Large Millimeter/submillimeter Array.
MAGNETOHYDRODYNAMIC EFFECTS ON PULSED YOUNG STELLAR OBJECT JETS. I. 2.5D SIMULATIONS
Hansen, E. C.; Frank, A.; Hartigan, P.
2015-02-10
In this paper, we explore the dynamics of radiative axisymmetric magnetohydrodynamic (MHD) jets at high resolution using adaptive mesh refinement methods. The goal of the study is to determine both the dynamics and emission properties of such jets. To that end, we have implemented microphysics enabling us to produce synthetic maps of Hα and [S II]. The jets are pulsed either sinusoidally or randomly via a time-dependent ejection velocity which leads to a complicated structure of internal shocks and rarefactions as has been seen in previous simulations. The high resolution of our simulations allows us to explore in great detail the effect of pinch forces (due to the jet's toroidal magnetic field) within the ''working surfaces'' where pulses interact. We map the strong Hα emission marking shock fronts and the strong [S II] emission inside cooling regions behind shocks as observed with high-resolution images of jets. We find that pinch forces in the stronger field cases produce additional emission regions along the axis as compared with purely hydrodynamic runs. These simulations are a first step to understanding the full three-dimensional emission properties of radiative MHD jets.
An overview of unconstrained free boundary problems
Figalli, Alessio; Shahgholian, Henrik
2015-01-01
In this paper, we present a survey concerning unconstrained free boundary problems of type where B1 is the unit ball, Ω is an unknown open set, F1 and F2 are elliptic operators (admitting regular solutions), and is a functions space to be specified in each case. Our main objective is to discuss a unifying approach to the optimal regularity of solutions to the above matching problems, and list several open problems in this direction. PMID:26261367
Cohen, O.; Kashyap, V. L.; Drake, J. J.; Garraffo, C.; Sokolov, I. V.; Gombosi, T. I.
2011-05-20
We carry out the first time-dependent numerical magnetohydrodynamic modeling of an extrasolar planetary system to study the interaction of the stellar magnetic field and wind with the planetary magnetosphere and outflow. We base our model on the parameters of the HD 189733 system, which harbors a close-in giant planet. Our simulation reveals a highly structured stellar corona characterized by sectors with different plasma properties. The star-planet interaction (SPI) varies in magnitude and complexity, depending on the planetary phase, planetary magnetic field strength, and the relative orientation of the stellar and planetary fields. It also reveals a long, comet-like tail which is a result of the wrapping of the planetary magnetospheric tail by its fast orbital motion. A reconnection event occurs at a specific orbital phase, causing mass loss from the planetary magnetosphere that can generate a hot spot on the stellar surface. The simulation also shows that the system has sufficient energy to produce hot spots observed in Ca II lines in giant planet hosting stars. However, the short duration of the reconnection event suggests that such SPI cannot be observed persistently.
The free boundary Euler equations with large surface tension
NASA Astrophysics Data System (ADS)
Disconzi, Marcelo M.; Ebin, David G.
2016-07-01
We study the free boundary Euler equations with surface tension in three spatial dimensions, showing that the equations are well-posed if the coefficient of surface tension is positive. Then we prove that under natural assumptions, the solutions of the free boundary motion converge to solutions of the Euler equations in a domain with fixed boundary when the coefficient of surface tension tends to infinity.
Development of an Auto-Convergent Free-Boundary Axisymmetric Equilibrium Solver
Huang, J.; Menard, J.
2006-01-01
The calculation of the magnetic flux given an assumed value for the current profile in axisymmetric toroidal plasmas is essential in studying the effects of various magnetohydrodynamic (MHD) instabilities upon controlled fusion. To this end, an iterative, modular algorithm coupled with a fast, direct elliptic solver for the Grad-Shafranov equation has been used to reconstruct the desired free-boundary equilibrium solution. This free-boundary Grad-Shafranov (FBGS) equilibrium algorithm is modified with the application of the von Hagenow method for determining the flux on the computational boundary, greatly reducing the time cost from O(N3) to O(N2 ln N) machine operations as compared to current Green’s function methods. The inherent variance in implementing the von Hagenow method gives a mean error bound of 0.1 percent with respect to the normal Green’s method. The improvements will allow the grid resolution to be increased efficiently and automatically to reduce the maximum Grad-Shafranov error to values needed for accurate stability calculations on a more effective time scale.
Free Boundary Problems for a Lotka-Volterra Competition System
NASA Astrophysics Data System (ADS)
Wang, Mingxin; Zhao, Jingfu
2014-09-01
In this paper we investigate two free boundary problems for a Lotka-Volterra type competition model in one space dimension. The main objective is to understand the asymptotic behavior of the two competing species spreading via a free boundary. We prove a spreading-vanishing dichotomy, namely the two species either successfully spread to the entire space as time t goes to infinity and survive in the new environment, or they fail to establish and die out in the long run. The long time behavior of the solutions and criteria for spreading and vanishing are also obtained. This paper is an improvement and extension of J. Guo and C. Wu.
Nonlinear magnetohydrodynamic stability
NASA Technical Reports Server (NTRS)
Bauer, F.; Betancourt, O.; Garabedian, P.
1981-01-01
The computer code developed by Bauer et al. (1978) for the study of the magnetohydrodynamic equilibrium and stability of a plasma in toroidal geometry is extended so that the growth rates of instabilities may be estimated more accurately. The original code, which is based on the variational principle of ideal magnetohydrodynamics, is upgraded by the introduction of a nonlinear formula for the growth rate of an unstable mode which acts as a quantitative measure of instability that is important in estimating numerical errors. The revised code has been applied to the determination of the nonlinear saturation, ballooning modes and beta limits for tokamaks, stellarators and torsatrons.
NASA Technical Reports Server (NTRS)
1995-01-01
This eerie, dark structure, resembling an imaginary sea serpent's head, is a column of cool molecular hydrogen gas (two atoms of hydrogen in each molecule) and dust that is an incubator for new stars. The stars are embedded inside finger-like protrusions extending from the top of the nebula. Each 'fingertip' is somewhat larger than our own solar system. The pillar is slowly eroding away by the ultraviolet light from nearby hot stars, a process called 'photoevaporation.' As it does, small globules of especially dense gas buried within the cloud is uncovered. These globules have been dubbed 'EGGs' -- an acronym for 'Evaporating Gaseous Globules.' The shadows of the EGGs protect gas behind them, resulting in the finger-like structures at the top of the cloud. Forming inside at least some of the EGGs are embryonic stars -- stars that abruptly stop growing when the EGGs are uncovered and they are separated from the larger reservoir of gas from which they were drawing mass. Eventually the stars emerge, as the EGGs themselves succumb to photoevaporation. The stellar EGGS are found, appropriately enough, in the 'Eagle Nebula' (also called M16 -- the 16th object in Charles Messier's 18th century catalog of 'fuzzy' permanent objects in the sky), a nearby star-forming region 7,000 light-years away in the constellation Serpens. The picture was taken on April 1, 1995 with the Hubble Space Telescope Wide Field and Planetary Camera 2. The color image is constructed from three separate images taken in the light of emission from different types of atoms. Red shows emission from singly-ionized sulfur atoms. Green shows emission from hydrogen. Blue shows light emitted by doubly-ionized oxygen atoms.
A semilinear parabolic system with a free boundary
NASA Astrophysics Data System (ADS)
Wang, Mingxin; Zhao, Yonggang
2015-12-01
This paper deals with a semilinear parabolic system with reaction terms {v^p, u^q} and a free boundary {x = s(t)} in one space dimension, where {s(t)} evolves according to the free boundary condition {s'(t) = -μ(u_x + ρ v_x)}. The main aim of this paper was to study the existence, uniqueness, regularity and long-time behavior of positive solution (maximal positive solution). Firstly, we prove that this problem has a unique positive solution when {p, q ≥ 1}, and a (unique) maximal positive solution when {p < 1} or {q < 1}. Then, we study the regularity of {(u,v)} and {s}. At last, we discuss the global existence, finite-time blowup of the unique positive solution (maximal positive solution) and long-time behavior of bounded global solution.
A free boundary approach to shape optimization problems
Bucur, D.; Velichkov, B.
2015-01-01
The analysis of shape optimization problems involving the spectrum of the Laplace operator, such as isoperimetric inequalities, has known in recent years a series of interesting developments essentially as a consequence of the infusion of free boundary techniques. The main focus of this paper is to show how the analysis of a general shape optimization problem of spectral type can be reduced to the analysis of particular free boundary problems. In this survey article, we give an overview of some very recent technical tools, the so-called shape sub- and supersolutions, and show how to use them for the minimization of spectral functionals involving the eigenvalues of the Dirichlet Laplacian, under a volume constraint. PMID:26261362
Free-boundary toroidal Alfvén eigenmodes
NASA Astrophysics Data System (ADS)
Chen, Eugene Y.; Berk, H. L.; Breizman, B.; Zheng, L. J.
2011-05-01
A numerical study is presented for the n = 1 free-boundary toroidal Alfvén eigenmodes (TAE) in tokamaks, which shows that there is considerable sensitivity of n = 1 modes to the position of the conducting wall. An additional branch of the TAE is shown to emerge from the upper continuum as the ratio of conducting wall radius to plasma radius increases. Such phenomena arise in plasma equilibria with both circular and shaped cross sections, where the shaped profile studied here is similar to that found in Alcator C-Mod.
Quiescent H-Mode 3D MHD Free-Boundary Equilibrium
NASA Astrophysics Data System (ADS)
Cooper, W. Anthony; Graves, Jonathan P.; Duval, Basil P.; Porte, Laurie; Sauter, Olivier; Tran, Trach-Minh; Brunetti, Daniele; Pfefferle, David; Raghunathan, Madhusudan; Faustin, Jonathan M.; Patten, Hamish; Kleiner, Andreas; Reimerdes, Holger
2015-11-01
Free boundary magnetohydrodynamic equilibrium states with spontaneous three dimensional deformations of the plasma-vacuum interface are computed with the 3D VMEC solver [Hirshman et al., Comput. Phys. Commun. 43 (1986) 143]. The structures we have obtained have the appearance of saturated ideal external kink/peeling modes. Large edge pressure gradients yield toroidal mode number n = 1 corrugations when the edge bootstrap current is large and n = 4 distortions when this current is small. The deformations of the plasma boundary region induces a nonaxisymmetric Pfirsch-Schlüter current that drives a field-aligned current ribbon which is consistent with experimental observations reported. We claim that the equilibrium states we compute model the Edge Harmonic Oscillation [K.H. Burrell et al., Phys. Plasmas 22 (2005) 021805. W.M. Solomon et al., Phys. Rev. Lett. 113 (2014) 135001] observed on DIII-D and the Outer Mode [E.R. Solano et al., Phys. Rev. Lett. 104 (2014) 135001] found in JET during Quiescent H-mode operation. This work was supported in part by the Swiss National Science Foundation.
Analysis of a free-boundary tumor model with angiogenesis
NASA Astrophysics Data System (ADS)
Friedman, Avner; Lam, King-Yeung
2015-12-01
We consider a free boundary problem for a spherically symmetric tumor with free boundary r < R (t). In order to receive nutrients u the tumor attracts blood vessel at a rate proportional to α (t), so that ∂u/∂r + α (t) (u - u bar) = 0 holds on the boundary, where u bar is the nutrient concentration outside the tumor. A parameter μ in the model is proportional to the 'aggressiveness' of the tumor. When α is a constant, the existence and uniqueness of stationary solution is proved. For the more general situation when α depends on time, we show, under various conditions (that are always satisfied if μ is small), that (i) R (t) remains bounded if α (t) remains bounded; (ii) limt→∞ R (t) = 0 if limt→∞ α (t) = 0; and (iii) lim inf t → ∞ R (t) > 0 if lim inf t → ∞ α (t) > 0. Surprisingly, we exhibit solutions (when μ is not small) where α (t) → 0 exponentially in t while R (t) → ∞ exponentially in t. Finally, we prove the global asymptotic stability of steady state when μ is sufficiently small.
Stalactite growth as a free-boundary problem
NASA Astrophysics Data System (ADS)
Short, Martin B.; Baygents, James C.; Goldstein, Raymond E.
2005-08-01
Stalactites, the most familiar structures found hanging from the ceilings of limestone caves, grow by the precipitation of calcium carbonate from within a thin film of fluid flowing down their surfaces. We have recently shown [M. B. Short, J. C. Baygents, J. W. Beck, D. A. Stone, R. S. Toomey III, and R. E. Goldstein, "Stalactite growth as a free-boundary problem: A geometric law and its Platonic ideal," Phys. Rev. Lett. 94, 018501 (2005)] that the combination of thin-film fluid dynamics, calcium carbonate chemistry, and carbon dioxide diffusion and outgassing leads to a local geometric growth law for the surface evolution which quantitatively explains the shapes of natural stalactites. Here we provide details of this free-boundary calculation, exploiting a strong separation of time scales among that for diffusion within the layer, contact of a fluid parcel with the growing surface, and growth. When the flow rate, the scale of the stalactite, and the chemistry are in the ranges typically found in nature, the local growth rate is proportional to the local thickness of the fluid layer, itself determined by Stokes flow over the surface. Numerical studies of this law establish that a broad class of initial conditions is attracted to an ideal universal shape, whose mathematical form is found analytically. Statistical analysis of stalactite shapes from Kartchner Caverns (Benson, AZ) shows excellent agreement between the average shape of natural stalactites and the ideal shape. Generalizations of these results to nonaxisymmetric speleothems are discussed.
An elliptic hyperbolic free boundary problem modelling cancer therapy
NASA Astrophysics Data System (ADS)
Tao, Youshan; Chen, Miaojun
2006-02-01
In this paper we study a free boundary problem modelling the growth of an avascular tumour with drug application. The tumour consists of two cell populations: live cells and dead cells. The densities of these cells satisfy a system of nonlinear first order hyperbolic equations. The tumour surface is a moving boundary, which satisfies an integro-differential equation. The nutrient concentration and the drug concentration satisfy nonlinear diffusion equations. The nutrient drives the growth of the tumour, whereas the drug is capable of killing cells with Michaelis-Menten kinetics. We prove that this free boundary problem has a unique global solution. Furthermore, we investigate the combined effects of a drug and a nutrient on an avascular tumour growth. We prove that the tumour shrinks to a necrotic core with radius Rs > 0 and that the global solution converges to a trivial steady-state solution under some natural assumptions on the model parameters. We also prove that an untreated tumour shrinks to a dead core or continually grows to an infinite size, which depends on the different parameter conditions.
A non-local free boundary problem arising in a theory of financial bubbles
Berestycki, Henri; Monneau, Regis; Scheinkman, José A.
2014-01-01
We consider an evolution non-local free boundary problem that arises in the modelling of speculative bubbles. The solution of the model is the speculative component in the price of an asset. In the framework of viscosity solutions, we show the existence and uniqueness of the solution. We also show that the solution is convex in space, and establish several monotonicity properties of the solution and of the free boundary with respect to parameters of the problem. To study the free boundary, we use, in particular, the fact that the odd part of the solution solves a more standard obstacle problem. We show that the free boundary is and describe the asymptotics of the free boundary as c, the cost of transacting the asset, goes to zero. PMID:25288815
The diffusive logistic equation with a free boundary and sign-changing coefficient
NASA Astrophysics Data System (ADS)
Wang, Mingxin
2015-02-01
This short paper concerns a diffusive logistic equation with a free boundary and sign-changing coefficient, which is formulated to study the spread of an invasive species, where the free boundary represents the expanding front. A spreading-vanishing dichotomy is derived, namely the species either successfully spreads to the right-half-space as time t → ∞ and survives (persists) in the new environment, or it fails to establish itself and will extinct in the long run. The sharp criteria for spreading and vanishing are also obtained. When spreading happens, we estimate the asymptotic spreading speed of the free boundary.
NASA Astrophysics Data System (ADS)
Garraffo, C.; Cohen, O.; Drake, J. J.; Downs, C.
2013-02-01
We study the influence of the spatial resolution on scales of 5° and smaller of solar surface magnetic field maps on global magnetohydrodynamic solar wind models, and on a model of coronal heating and X-ray emission. We compare the solutions driven by a low-resolution Wilcox Solar Observatory magnetic map, the same map with spatial resolution artificially increased by a refinement algorithm, and a high-resolution Solar and Heliospheric Observatory Michelson Doppler Imager map. We find that both the wind structure and the X-ray morphology are affected by the fine-scale surface magnetic structure. Moreover, the X-ray morphology is dominated by the closed loop structure between mixed polarities on smaller scales and shows significant changes between high- and low-resolution maps. We conclude that three-dimensional modeling of coronal X-ray emission has greater surface magnetic field spatial resolution requirements than wind modeling, and can be unreliable unless the dominant mixed polarity magnetic flux is properly resolved.
Garraffo, C.; Cohen, O.; Drake, J. J.; Downs, C.
2013-02-10
We study the influence of the spatial resolution on scales of 5 Degree-Sign and smaller of solar surface magnetic field maps on global magnetohydrodynamic solar wind models, and on a model of coronal heating and X-ray emission. We compare the solutions driven by a low-resolution Wilcox Solar Observatory magnetic map, the same map with spatial resolution artificially increased by a refinement algorithm, and a high-resolution Solar and Heliospheric Observatory Michelson Doppler Imager map. We find that both the wind structure and the X-ray morphology are affected by the fine-scale surface magnetic structure. Moreover, the X-ray morphology is dominated by the closed loop structure between mixed polarities on smaller scales and shows significant changes between high- and low-resolution maps. We conclude that three-dimensional modeling of coronal X-ray emission has greater surface magnetic field spatial resolution requirements than wind modeling, and can be unreliable unless the dominant mixed polarity magnetic flux is properly resolved.
Spreading and vanishing in the diffusive prey-predator model with a free boundary
NASA Astrophysics Data System (ADS)
Wang, Mingxin
2015-06-01
This paper deals with the diffusive Lotka-Volterra type prey-predator model with a free boundary over a one dimensional habitat. This problem may be used to describe the interaction between indigenous species and invasive species and the spreading of such two species, with the free boundary representing the expanding front. Our main purpose is to study the spreading and vanishing phenomena and long time behaviors of prey and predator.
Boquist, Carl W.; Marchant, David D.
1978-01-01
A ceramic-metal composite suitable for use in a high-temperature environment consists of a refractory ceramic matrix containing 10 to 50 volume percent of a continuous high-temperature metal reinforcement. In a specific application of the composite, as an electrode in a magnetohydrodynamic generator, the one surface of the electrode which contacts the MHD fluid may have a layer of varying thickness of nonreinforced refractory ceramic for electrode temperature control. The side walls of the electrode may be coated with a refractory ceramic insulator. Also described is an electrode-insulator system for a MHD channel.
Magnetohydrodynamic instability
NASA Technical Reports Server (NTRS)
Priest, E. R.; Cargill, P.; Forbes, T. G.; Hood, A. W.; Steinolfson, R. S.
1986-01-01
There have been major advances in the theory of magnetic reconnection and of magnetic instability, with important implications for the observations, as follows: (1) Fast and slow magnetic shock waves are produced by the magnetohydrodynamics of reconnection and are potential particle accelerators. (2) The impulsive bursty regime of reconnection gives a rapid release of magnetic energy in a series of bursts. (3) The radiative tearing mode creates cool filamentary structures in the reconnection process. (4) The stability analyses imply that an arcade can become unstable when either its height or twist of plasma pressure become too great.
NASA Technical Reports Server (NTRS)
Carpenter, Kenneth
2007-01-01
The Stellar Imager (SI) is one of NASA's "Vision Missions" - concepts for future, space-based, strategic missions that could enormously increase our capabilities for observing the Cosmos. SI is designed as a UV/Optical Interferometer which will enable 0.1 milli-arcsecond (mas) spectral imaging of stellar surfaces and, via asteroseismology, stellar interiors and of the Universe in general. The ultra-sharp images of the Stellar Imager will revolutionize our view of many dynamic astrophysical processes by transforming point sources into extended sources, and snapshots into evolving views. SI, with a characteristic angular resolution of 0.1 milli-arcseconds at 2000 Angstroms, represents an advance in image detail of several hundred times over that provided by the Hubble Space Telescope. The Stellar Imager will zoom in on what today-with few exceptions - we only know as point sources, revealing processes never before seen, thus providing a tool as fundamental to astrophysics as the microscope is to the study of life on Earth. SI's science focuses on the role of magnetism in the Universe, particularly on magnetic activity on the surfaces of stars like the Sun. It's prime goal is to enable long-term forecasting of solar activity and the space weather that it drives, in support of the Living With a Star program in the Exploration Era. SI will also revolutionize our understanding of the formation of planetary systems, of the habitability and climatology of distant planets, and of many magneto-hydrodynamically controlled processes in the Universe. Stellar Imager is included as a "Flagship and Landmark Discovery Mission" in the 2005 Sun Solar System Connection (SSSC) Roadmap and as a candidate for a "Pathways to Life Observatory" in the Exploration of the Universe Division (EUD) Roadmap (May, 2005) and as such is a candidate mission for the 2025-2030 timeframe. An artist's drawing of the current "baseline" concept for SI is presented.
Marchant, David D.; Killpatrick, Don H.
1978-01-01
An electrode capable of withstanding high temperatures and suitable for use as a current collector in the channel of a magnetohydrodynamic (MHD) generator consists of a sintered powdered metal base portion, the upper surface of the base being coated with a first layer of nickel aluminide, an intermediate layer of a mixture of nickel aluminide - refractory ceramic on the first layer and a third or outer layer of a refractory ceramic material on the intermediate layer. The sintered powdered metal base resists spalling by the ceramic coatings and permits greater electrode compliance to thermal shock. The density of the powdered metal base can be varied to allow optimization of the thermal conductivity of the electrode and prevent excess heat loss from the channel.
On some free boundary problems of the prey-predator model
NASA Astrophysics Data System (ADS)
Wang, Mingxin
In this paper we investigate some free boundary problems for the Lotka-Volterra type prey-predator model in one space dimension. The main objective is to understand the asymptotic behavior of the two species (prey and predator) spreading via a free boundary. We prove a spreading-vanishing dichotomy, namely the two species either successfully spread to the entire space as time t goes to infinity and survive in the new environment, or they fail to establish and die out in the long run. The long time behavior of solution and criteria for spreading and vanishing are also obtained. Finally, when spreading successfully, we provide an estimate to show that the spreading speed (if exists) cannot be faster than the minimal speed of traveling wavefront solutions for the prey-predator model on the whole real line without a free boundary.
Spreading-vanishing dichotomy in a diffusive logistic model with a free boundary, II
NASA Astrophysics Data System (ADS)
Du, Yihong; Guo, Zongming
We study the diffusive logistic equation with a free boundary in higher space dimensions and heterogeneous environment. Such a model may be used to describe the spreading of a new or invasive species, with the free boundary representing the expanding front. For simplicity, we assume that the environment and the solution are radially symmetric. In the special case of one space dimension and homogeneous environment, this free boundary problem was investigated in Du and Lin (2010) [10]. We prove that the spreading-vanishing dichotomy established in Du and Lin (2010) [10] still holds in the more general and ecologically realistic setting considered here. Moreover, when spreading occurs, we obtain best possible upper and lower bounds for the spreading speed of the expanding front. When the environment is asymptotically homogeneous at infinity, these two bounds coincide. Our results indicate that the asymptotic spreading speed determined by this model does not depend on the spatial dimension.
The minimal habitat size for spreading in a weak competition system with two free boundaries
NASA Astrophysics Data System (ADS)
Wu, Chang-Hong
2015-08-01
In this paper, we focus on the dynamics for a Lotka-Volterra type weak competition system with two free boundaries, where free boundaries which may intersect each other as time evolves are used to describe the spreading of two competing species, respectively. In the weak competition case, the dynamics of this model can be classified into four cases, which forms a spreading-vanishing quartering. The notion of the minimal habitat size for spreading is introduced to determine if species can always spread. Some sufficient conditions for spreading and vanishing are established. Also, when spreading occurs, some rough estimates for spreading speed and the long-time behavior of solutions are established.
The free boundary problem describing information diffusion in online social networks
NASA Astrophysics Data System (ADS)
Lei, Chengxia; Lin, Zhigui; Wang, Haiyan
In this paper we consider a free boundary problem for a reaction-diffusion logistic equation with a time-dependent growth rate. Such a problem arises in the modeling of information diffusion in online social networks, with the free boundary representing the spreading front of news among users. We present several sharp thresholds for information diffusion that either lasts forever or suspends in finite time. In the former case, we give the asymptotic spreading speed which is determined by a corresponding elliptic equation.
NASA Astrophysics Data System (ADS)
Li, Fang; Liang, Xing; Shen, Wenxian
2016-08-01
In this series of papers, we investigate the spreading and vanishing dynamics of time almost periodic diffusive KPP equations with free boundaries. Such equations are used to characterize the spreading of a new species in time almost periodic environments with free boundaries representing the spreading fronts. In the first part of the series, we showed that a spreading-vanishing dichotomy occurs for such free boundary problems (see [16]). In this second part of the series, we investigate the spreading speeds of such free boundary problems in the case that the spreading occurs. We first prove the existence of a unique time almost periodic semi-wave solution associated to such a free boundary problem. Using the semi-wave solution, we then prove that the free boundary problem has a unique spreading speed.
NASA Astrophysics Data System (ADS)
Erdélyi, R.
2007-07-01
The heating of solar atmosphere from chromosphere to corona is one of the key fundamental and yet unresolved questions of modern space and plasma physics. In spite of the multi-fold efforts spanning over half a century including the many superb technological advances and theoretical developments (both analytical and computational) the unveiling of the subtle of coronal heating still remains an exciting job for the 21st century! In the present paper I review the various popular heating mechanisms put forward in the existing extensive literature. The heating processes are, somewhat arbitrarily, classified as hydrodynamic (HD), magnetohydrodynamic (MHD) or kinetic based on the model medium. These mechanisms are further divided based on the time scales of the ultimate dissipation involved (i.e. AC and DC heating, turbulent heating). In particular, attention is paid to discuss shock dissipation, Landau damping, mode coupling, resonant absorption, phase mixing, and, reconnection. Finally, I briefly review the various observational consequences of the many proposed heating mechanisms and confront them with high-resolution ground-based and satellite data currently available.
A priori estimates for the free boundary problem of incompressible neo-Hookean elastodynamics
NASA Astrophysics Data System (ADS)
Hao, Chengchun; Wang, Dehua
2016-07-01
A free boundary problem for the incompressible neo-Hookean elastodynamics is studied in two and three spatial dimensions. The a priori estimates in Sobolev norms of solutions with the physical vacuum condition are established through a geometrical point of view of Christodoulou and Lindblad (2000) [3]. Some estimates on the second fundamental form and velocity of the free surface are also obtained.
Magnetohydrodynamic Turbulence
NASA Astrophysics Data System (ADS)
Montgomery, David C.
2004-01-01
Magnetohydrodynamic (MHD) turbulence theory is modeled on neutral fluid (Navier-Stokes) turbulence theory, but with some important differences. There have been essentially no repeatable laboratory MHD experiments wherein the boundary conditions could be controlled or varied and a full set of diagnostics implemented. The equations of MHD are convincingly derivable only in the limit of small ratio of collision mean-free-paths to macroscopic length scales, an inequality that often goes the other way for magnetofluids of interest. Finally, accurate information on the MHD transport coefficients-and thus, the Reynolds-like numbers that order magnetofluid behavior-is largely lacking; indeed, the algebraic expressions used for such ingredients as the viscous stress tensor are often little more than wishful borrowing from fluid mechanics. The one accurate thing that has been done extensively and well is to solve the (strongly nonlinear) MHD equations numerically, usually in the presence of rectangular periodic boundary conditions, and then hope for the best when drawing inferences from the computations for those astrophysical and geophysical MHD systems for which some indisputably turbulent detailed data are available, such as the solar wind or solar prominences. This has led to what is perhaps the first field of physics for which computer simulations are regarded as more central to validating conclusions than is any kind of measurement. Things have evolved in this way due to a mixture of the inevitable and the bureaucratic, but that is the way it is, and those of us who want to work on the subject have to live with it. It is the only game in town, and theories that have promised more-often on the basis of some alleged ``instability''-have turned out to be illusory.
Qin, H.; Reiman, A.
1996-09-25
An analytic solution is obtained for free-boundary, high-beta equilibria in large aspect ratio tokamaks with a nearly circular plasma boundary. In the absence of surface currents at the plasma-vacuum interface, the free-boundary equilibrium solution introduces constraints arising from the need to couple to an external vacuum field which is physically realizable with a reasonable set of external field coils. This places a strong constraint on the pressure profiles that are consistent with a given boundary shape at high {epsilon}{beta}{sub p}. The equilibrium solution also provides information on the flux surface topology. The plasma is bounded by a separatrix. Increasing the plasma pressure at fixed total current causes the plasma aperture to decrease in a manner that is described.
NASA Astrophysics Data System (ADS)
Kusaka, Yoshiaki
2016-06-01
Stationary free boundary problems related to the formation of the interface between two immiscible fluids are investigated. The corresponding models were introduced originally by Shikhmurzaev to remove singularities arising in contact-line problems. We consider two models which are derived under different assumptions for the structure of the interfacial layer, and prove the existence of rotationally symmetric solutions in Hölder spaces for a sufficiently small rotationally symmetric external force.
Some free boundary problems involving non-local diffusion and aggregation
Carrillo, José Antonio; Vázquez, Juan Luis
2015-01-01
We report on recent progress in the study of evolution processes involving degenerate parabolic equations which may exhibit free boundaries. The equations we have selected follow two recent trends in diffusion theory: considering anomalous diffusion with long-range effects, which leads to fractional operators or other operators involving kernels with large tails; and the combination of diffusion and aggregation effects, leading to delicate long-term equilibria whose description is still incipient. PMID:26261360
A free boundary problem for steady small plaques in the artery and their stability
NASA Astrophysics Data System (ADS)
Friedman, Avner; Hao, Wenrui; Hu, Bei
2015-08-01
Atherosclerosis is a leading cause of death in the United States and worldwide; it originates from a plaque which builds up in the artery. In this paper, we consider a simplified model of plaque growth involving LDL and HDL cholesterols, macrophages and foam cells, which satisfy a coupled system of PDEs with a free boundary, the interface between the plaque and the blood flow. We prove that there exist small radially symmetric stationary plaques and establish a sharp condition that ensures their stability. We also determine necessary and sufficient conditions under which a small initial plaque will shrink and disappear, or persist for all times.
Furth, H.P.; Ludescher, C.
1984-08-01
The present paper briefly reviews the subject of tokamak-stellarator and pinch-stellarator hybrids, and points to two interesting new possibilities: compact-torus-stellarators and mirror-stellarators.
Helical axis stellarator equilibrium model
Koniges, A.E.; Johnson, J.L.
1985-02-01
An asymptotic model is developed to study MHD equilibria in toroidal systems with a helical magnetic axis. Using a characteristic coordinate system based on the vacuum field lines, the equilibrium problem is reduced to a two-dimensional generalized partial differential equation of the Grad-Shafranov type. A stellarator-expansion free-boundary equilibrium code is modified to solve the helical-axis equations. The expansion model is used to predict the equilibrium properties of Asperators NP-3 and NP-4. Numerically determined flux surfaces, magnetic well, transform, and shear are presented. The equilibria show a toroidal Shafranov shift.
Ising antiferromagnet on a finite triangular lattice with free boundary conditions
NASA Astrophysics Data System (ADS)
Kim, Seung-Yeon
2015-11-01
The exact integer values for the density of states of the Ising model on an equilateral triangular lattice with free boundary conditions are evaluated up to L = 24 spins on a side for the first time by using the microcanonical transfer matrix. The total number of states is 2 N s = 2300 ≈ 2.037 × 1090 for L = 24, where N s = L( L+1)/2 is the number of spins. Classifying all 2300 spin states according to their energy values is an enormous work. From the density of states, the exact partition function zeros in the complex temperature plane of the triangular-lattice Ising model are evaluated. Using the density of states and the partition function zeros, we investigate the properties of the triangularlattice Ising antiferromagnet. The scaling behavior of the ground-state entropy and the form of the correlation length at T = 0 are studied for the triangular-lattice Ising antiferromagnet with free boundary conditions. Also, the scaling behavior of the Fisher edge singularity is investigated.
On a Free Boundary Problem for the Curvature Flow with Driving Force
NASA Astrophysics Data System (ADS)
Guo, Jong-Shenq; Matano, Hiroshi; Shimojo, Masahiko; Wu, Chang-Hong
2016-03-01
We study a free boundary problem associated with the curvature dependent motion of planar curves in the upper half plane whose two endpoints slide along the horizontal axis with prescribed fixed contact angles. Our first main result concerns the classification of solutions; every solution falls into one of the three categories, namely, area expanding, area bounded and area shrinking types. We then study in detail the asymptotic behavior of solutions in each category. Among other things we show that solutions are asymptotically self-similar both in the area expanding and the area shrinking cases, while solutions converge to either a stationary solution or a traveling wave in the area bounded case. We also prove results on the concavity properties of solutions. One of the main tools of this paper is the intersection number principle, however in order to deal with solutions with free boundaries, we introduce what we call "the extended intersection number principle", which turns out to be exceedingly useful in handling curves with moving endpoints.
NASA Astrophysics Data System (ADS)
Kent, Stuart; Venkataramani, Shankar C.
2014-07-01
We use conformal maps to study a free boundary problem for a two-fluid electromechanical system, where the interface between the fluids is determined by the combined effects of electrostatic forces, gravity, and surface tension. The free boundary in our system develops sharp corners or singularities in certain parameter regimes, and this is an impediment to using existing "single-scale" numerical conformal mapping methods. The difficulty is due to the phenomenon of crowding, i.e., the tendency of nodes in the preimage plane to concentrate near the sharp regions of the boundary, leaving the smooth regions of the boundary poorly resolved. A natural idea is to exploit the scale separation between the sharp regions and smooth regions to solve for each region separately and then stitch the solutions together. However, this is not straightforward as conformal maps are rigid "global" objects, and it is not obvious how one would patch two conformal maps together to obtain a new conformal map. We develop a "multiscale" (i.e., adaptive) conformal mapping method that allows us to carry out this program of stitching conformal maps on different scales together. We successfully apply our method to the electromechanical model problem.
A Schur complement formulation for solving free-boundary, Stefan problems of phase change
NASA Astrophysics Data System (ADS)
Lun, Lisa; Yeckel, Andrew; Derby, Jeffrey J.
2010-10-01
A Schur complement formulation that utilizes a linear iterative solver is derived to solve a free-boundary, Stefan problem describing steady-state phase change via the Isotherm-Newton approach, which employs Newton's method to simultaneously and efficiently solve for both interface and field equations. This formulation is tested alongside more traditional solution strategies that employ direct or iterative linear solvers on the entire Jacobian matrix for a two-dimensional sample problem that discretizes the field equations using a Galerkin finite-element method and employs a deforming-grid approach to represent the melt-solid interface. All methods demonstrate quadratic convergence for sufficiently accurate Newton solves, but the two approaches utilizing linear iterative solvers show better scaling of computational effort with problem size. Of these two approaches, the Schur formulation proves to be more robust, converging with significantly smaller Krylov subspaces than those required to solve the global system of equations. Further improvement of performance are made through approximations and preconditioning of the Schur complement problem. Hence, the new Schur formulation shows promise as an affordable, robust, and scalable method to solve free-boundary, Stefan problems. Such models are employed to study a wide array of applications, including casting, welding, glass forming, planetary mantle and glacier dynamics, thermal energy storage, food processing, cryosurgery, metallurgical solidification, and crystal growth.
NASA Astrophysics Data System (ADS)
Brun, A. S.; García, R. A.; Houdek, G.; Nandy, D.; Pinsonneault, M.
2015-12-01
We discuss how recent advances in observations, theory and numerical simulations have allowed the stellar community to progress in its understanding of stellar convection, rotation and magnetism and to assess the degree to which the Sun and other stars share similar dynamical properties. Ensemble asteroseismology has become a reality with the advent of large time domain studies, especially from space missions. This new capability has provided improved constraints on stellar rotation and activity, over and above that obtained via traditional techniques such as spectropolarimetry or CaII H&K observations. New data and surveys covering large mass and age ranges have provided a wide parameter space to confront theories of stellar magnetism. These new empirical databases are complemented by theoretical advances and improved multi-D simulations of stellar dynamos. We trace these pathways through which a lucid and more detailed picture of magnetohydrodynamics of solar-like stars is beginning to emerge and discuss future prospects.
Magnetohydrodynamic power generation
NASA Technical Reports Server (NTRS)
Smith, J. L.
1984-01-01
Magnetohydrodynamic (MHD) Power Generation is a concise summary of MHD theory, history, and future trends. Results of the major international MHD research projects are discussed. Data from MHD research is included. Economics of initial and operating costs are considered.
Experiments in Magnetohydrodynamics
ERIC Educational Resources Information Center
Rayner, J. P.
1970-01-01
Describes three student experiments in magnetohydrodynamics (MHD). In these experiments, it was found that the electrical conductivity of the local water supply was sufficient to demonstrate effectively some of the features of MHD flowmeters, generators, and pumps. (LC)
Gyroscopic analog for magnetohydrodynamics
Holm, D.D.
1981-01-01
The gross features of plasma equilibrium and dynamics in the ideal magnetohydrodynamics (MHD) model can be understood in terms of a dynamical system which closely resembles the equations for a deformable gyroscope.
Magnetohydrodynamic fluidic system
Lee, Abraham P.; Bachman, Mark G.
2004-08-24
A magnetohydrodynamic fluidic system includes a reagent source containing a reagent fluid and a sample source containing a sample fluid that includes a constituent. A reactor is operatively connected to the supply reagent source and the sample source. MHD pumps utilize a magnetohydrodynamic drive to move the reagent fluid and the sample fluid in a flow such that the reagent fluid and the sample fluid form an interface causing the constituent to be separated from the sample fluid.
Free boundary problems in shock reflection/diffraction and related transonic flow problems.
Chen, Gui-Qiang; Feldman, Mikhail
2015-09-13
Shock waves are steep wavefronts that are fundamental in nature, especially in high-speed fluid flows. When a shock hits an obstacle, or a flying body meets a shock, shock reflection/diffraction phenomena occur. In this paper, we show how several long-standing shock reflection/diffraction problems can be formulated as free boundary problems, discuss some recent progress in developing mathematical ideas, approaches and techniques for solving these problems, and present some further open problems in this direction. In particular, these shock problems include von Neumann's problem for shock reflection-diffraction by two-dimensional wedges with concave corner, Lighthill's problem for shock diffraction by two-dimensional wedges with convex corner, and Prandtl-Meyer's problem for supersonic flow impinging onto solid wedges, which are also fundamental in the mathematical theory of multidimensional conservation laws. PMID:26261363
Free boundary problems in shock reflection/diffraction and related transonic flow problems
Chen, Gui-Qiang; Feldman, Mikhail
2015-01-01
Shock waves are steep wavefronts that are fundamental in nature, especially in high-speed fluid flows. When a shock hits an obstacle, or a flying body meets a shock, shock reflection/diffraction phenomena occur. In this paper, we show how several long-standing shock reflection/diffraction problems can be formulated as free boundary problems, discuss some recent progress in developing mathematical ideas, approaches and techniques for solving these problems, and present some further open problems in this direction. In particular, these shock problems include von Neumann's problem for shock reflection–diffraction by two-dimensional wedges with concave corner, Lighthill's problem for shock diffraction by two-dimensional wedges with convex corner, and Prandtl-Meyer's problem for supersonic flow impinging onto solid wedges, which are also fundamental in the mathematical theory of multidimensional conservation laws. PMID:26261363
Exact solution of two phase spherical Stefan problem with two free boundaries
NASA Astrophysics Data System (ADS)
Kavokin, Alexey A.; Nauryz, Targyn; Bizhigitova, Nazerke T.
2016-08-01
Solution of the heat equation in a spherical domain with two free boundaries (two-phase Stefan problem) when one of the subdomains degenerates at the initial time is considered. The use of conventional finite-difference methods in these cases is not expedient because of the degenerate domain. The solution is found in the form of combination of Integral Error functions series, [M. Sarsengeldin, and S. Kharin, Filomat, (2016), (in Press)] and then recurrent solvability of nonlinear algebraic equations for determining the coefficients of the series is proved. Such problems are of practical interest for the simulation of laser material processing as well for the modeling of thermal effects of electric arc that ignites during the opening of electric contacts [S. N. Kharin, and M. Sarsengeldin, ï£¡Influence of contact materials on phenomena in a short electrical arcï£¡, in Key Engineering Materials, Trans tech publications, Islamabad, Pakistan, 2012, pp. 321-329].
Development of a free-boundary tokamak equilibrium solver for advanced study of tokamak equilibria
NASA Astrophysics Data System (ADS)
Jeon, Young Mu
2015-09-01
A free-boundary Tokamak equilibrium solver (TES), developed for advanced study of tokamak equilibra, is described with two distinctive features. One is a generalized method to resolve the intrinsic axisymmetric instability, which is encountered in all equilibrium calculations with a freeboundary condition. The other is an extension to deal with a new divertor geometry such as snowflake or X divertors. For validations, the uniqueness of a solution is confirmed by the independence of variations in the computational domain, the mathematical correctness and accuracy of equilibrium profiles are checked by using a direct comparison with an analytic equilibrium known as a generalized Solov'ev equilibrium, and the governing force balance relation is tested by examining the intrinsic axisymmetric instabilities. As an application of an advanced equilibrium study, a snow-flake divertor configuration that requires a second-order zero of the poloidal magnetic flux is discussed in the circumstance of the Korea superconducting tokamak advanced research (KSTAR) coil system.
Casimir force in the O(n→∞) model with free boundary conditions.
Dantchev, Daniel; Bergknoff, Jonathan; Rudnick, Joseph
2014-04-01
We present results for the temperature behavior of the Casimir force for a system with a film geometry with thickness L subject to free boundary conditions and described by the n→∞ limit of the O(n) model. These results extend over all temperatures, including the critical regime near the bulk critical temperature Tc, where the critical fluctuations determine the behavior of the force, and temperatures well below it, where its behavior is dictated by the Goldstone mode contributions. The temperature behavior when the absolute temperature, T, is a finite distance below Tc, up to a logarithmic-in-L proximity of the bulk critical temperature, is obtained both analytically and numerically; the critical behavior follows from numerics. The results resemble-but do not duplicate-the experimental curve behavior for the force obtained for He4 films. PMID:24827202
Comment on "Casimir force in the O (n →∞ ) model with free boundary conditions"
NASA Astrophysics Data System (ADS)
Diehl, H. W.; Grüneberg, Daniel; Hasenbusch, Martin; Hucht, Alfred; Rutkevich, Sergei B.; Schmidt, Felix M.
2015-02-01
In a recent paper by D. Dantchev, J. Bergknoff, and J. Rudnick [Phys. Rev. E 89, 042116 (2014), 10.1103/PhysRevE.89.042116], the problem of the Casimir force in the O (n ) model on a slab with free boundary conditions, investigated earlier by us [Europhys. Lett. 100, 10004 (2012), 10.1209/0295-5075/100/10004], is reconsidered using a mean-spherical model with separate constraints for each layer. The authors (i) question the applicability of the Ginzburg-Landau-Wilson approach to the low-temperature regime, arguing for the superiority of their model compared to the family of ϕ4 models A and B whose numerically exact solutions we determined both for values of the coupling constant 0
Comment on "Casimir force in the O(n→∞) model with free boundary conditions".
Diehl, H W; Grüneberg, Daniel; Hasenbusch, Martin; Hucht, Alfred; Rutkevich, Sergei B; Schmidt, Felix M
2015-02-01
In a recent paper by D. Dantchev, J. Bergknoff, and J. Rudnick [Phys. Rev. E 89, 042116 (2014)], the problem of the Casimir force in the O(n) model on a slab with free boundary conditions, investigated earlier by us [Europhys. Lett. 100, 10004 (2012)], is reconsidered using a mean-spherical model with separate constraints for each layer. The authors (i) question the applicability of the Ginzburg-Landau-Wilson approach to the low-temperature regime, arguing for the superiority of their model compared to the family of ϕ(4) models A and B whose numerically exact solutions we determined both for values of the coupling constant 0
Some free boundary problems in potential flow regime usinga based level set method
Garzon, M.; Bobillo-Ares, N.; Sethian, J.A.
2008-12-09
Recent advances in the field of fluid mechanics with moving fronts are linked to the use of Level Set Methods, a versatile mathematical technique to follow free boundaries which undergo topological changes. A challenging class of problems in this context are those related to the solution of a partial differential equation posed on a moving domain, in which the boundary condition for the PDE solver has to be obtained from a partial differential equation defined on the front. This is the case of potential flow models with moving boundaries. Moreover the fluid front will possibly be carrying some material substance which will diffuse in the front and be advected by the front velocity, as for example the use of surfactants to lower surface tension. We present a Level Set based methodology to embed this partial differential equations defined on the front in a complete Eulerian framework, fully avoiding the tracking of fluid particles and its known limitations. To show the advantages of this approach in the field of Fluid Mechanics we present in this work one particular application: the numerical approximation of a potential flow model to simulate the evolution and breaking of a solitary wave propagating over a slopping bottom and compare the level set based algorithm with previous front tracking models.
Three-dimensional magnetohydrodynamic equilibrium of quiescent H-modes in tokamak systems
NASA Astrophysics Data System (ADS)
Cooper, W. A.; Graves, J. P.; Duval, B. P.; Sauter, O.; Faustin, J. M.; Kleiner, A.; Lanthaler, S.; Patten, H.; Raghunathan, M.; Tran, T.-M.; Chapman, I. T.; Ham, C. J.
2016-06-01
Three dimensional free boundary magnetohydrodynamic equilibria that recover saturated ideal kink/peeling structures are obtained numerically. Simulations that model the JET tokamak at fixed < β > =1.7% with a large edge bootstrap current that flattens the q-profile near the plasma boundary demonstrate that a radial parallel current density ribbon with a dominant m /n = 5/1 Fourier component at {{I}\\text{t}}=2.2 MA develops into a broadband spectrum when the toroidal current I t is increased to 2.5 MA.
Free-Boundary 3D Equilibria and Resistive Wall Instabilities with Extended-MHD
NASA Astrophysics Data System (ADS)
Ferraro, N. M.
2015-11-01
The interaction of the plasma with external currents, either imposed or induced, is a critical element of a wide range of important tokamak phenomena, including resistive wall mode (RWM) stability and feedback control, island penetration and locking, and disruptions. A model of these currents may be included within the domain of extended-MHD codes in a way that preserves the self-consistency, scalability, and implicitness of their numerical methods. Such a model of the resistive wall and non-axisymmetric coils is demonstrated using the M3D-C1 code for a variety of applications, including RWMs, perturbed non-axisymmetric equilibria, and a vertical displacement event (VDE) disruption. The calculated free-boundary equilibria, which include Spitzer resistivity, rotation, and two-fluid effects, are compared to external magnetic and internal thermal measurements for several DIII-D discharges. In calculations of the perturbed equilibria in ELM suppressed discharges, the tearing response at the top of the pedestal is found to correlate with the onset of ELM suppression. Nonlinear VDE calculations, initialized using a vertically unstable DIII-D equilibrium, resolve in both space and time the currents induced in the wall and on the plasma surface, and also the currents flowing between the plasma and the wall. The relative magnitude of these contributions and the total impulse to the wall depend on the resistive wall time, although the maximum axisymmetric force on the wall over the course of the VDE is found to be essentially independent of the wall conductivity. This research was supported by US DOE contracts DE-FG02-95ER54309, DE-FC02-04ER54698 and DE-AC52-07NA27344.
NASA Astrophysics Data System (ADS)
Kawai, Yusuke; Yamada, Yoshio
2016-07-01
This paper deals with a free boundary problem for diffusion equation with a certain class of bistable nonlinearity which allows two positive stable equilibrium states as an ODE model. This problem models the invasion of a biological species and the free boundary represents the spreading front of its habitat. Our main interest is to study large-time behaviors of solutions for the free boundary problem. We will completely classify asymptotic behaviors of solutions and, in particular, observe two different types of spreading phenomena corresponding to two positive stable equilibrium states. Moreover, it will be proved that, if the free boundary expands to infinity, an asymptotic speed of the moving free boundary for large time can be uniquely determined from the related semi-wave problem.
NASA Technical Reports Server (NTRS)
Chiu, H.-Y. (Editor); Muriel, A.
1972-01-01
Aspects of normal stellar evolution are discussed together with evolution near the main sequence, stellar evolution from main sequence to white dwarf or carbon ignition, the structure of massive main-sequence stars, and problems of stellar stability and stellar pulsation. Other subjects considered include variable stars, white dwarfs, close binaries, novae, early supernova luminosity, neutron stars, the photometry of field horizontal-branch stars, and stellar opacity. Transport mechanisms in stars are examined together with thermonuclear reactions and nucleosynthesis, the instability problem in nuclear burning shells, stellar coalescence, and intense magnetic fields in astrophysics. Individual items are announced in this issue.
NASA Astrophysics Data System (ADS)
Wu, Junde
2016-04-01
In this paper we study a free boundary problem modeling tumor growth. The model consists of two elliptic equations describing nutrient diffusion and pressure distribution within tumors, respectively, and a first-order partial differential equation governing the free boundary, on which a Gibbs-Thomson relation is taken into account. We first show that the problem may have none, one or two radial stationary solutions depending on model parameters. Then by bifurcation analysis we show that there exist infinite many branches of non-radial stationary solutions bifurcating from given radial stationary solution. The result implies that cell-to-cell adhesiveness is the key parameter which plays a crucial role on tumor invasion.
MHD kink-driven instabilities in net-current-free stellarators
Rewoldt, G.; Johnson, J.L.
1984-02-01
The Pfirsch-Schlueter current, which is induced in a toroidal device to keep the plasma current diverence-free, is shown to drive a free-boundary instability in a model of a net-current-free ATF-1 stellarator if <..beta..> = 2.6%.
Cornetti, G.M.
1995-12-31
The 3D Navier-Stokes equations are solved via the Characteristic-Galerkin method extended to free boundary problems. A temporal discretization procedure is proposed for the case where a preferential direction to move mesh point exists, as in thin domains. Using a single layer of finite elements, the numerical results cover the so-called shallow water 2D approximation, showing the same wave propagation speed.
Thermoacoustic magnetohydrodynamic electrical generator
Wheatley, John C.; Swift, Gregory W.; Migliori, Albert
1986-01-01
A thermoacoustic magnetohydrodynamic electrical generator includes an intrinsically irreversible thermoacoustic heat engine coupled to a magnetohydrodynamic electrical generator. The heat engine includes an electrically conductive liquid metal as the working fluid and includes two heat exchange and thermoacoustic structure assemblies which drive the liquid in a push-pull arrangement to cause the liquid metal to oscillate at a resonant acoustic frequency on the order of 1,000 Hz. The engine is positioned in the field of a magnet and is oriented such that the liquid metal oscillates in a direction orthogonal to the field of the magnet, whereby an alternating electrical potential is generated in the liquid metal. Low-loss, low-inductance electrical conductors electrically connected to opposite sides of the liquid metal conduct an output signal to a transformer adapted to convert the low-voltage, high-current output signal to a more usable higher voltage, lower current signal.
Thermoacoustic magnetohydrodynamic electrical generator
Wheatley, J.C.; Swift, G.W.; Migliori, A.
1984-11-16
A thermoacoustic magnetohydrodynamic electrical generator includes an intrinsically irreversible thermoacoustic heat engine coupled to a magnetohydrodynamic electrical generator. The heat engine includes an electrically conductive liquid metal as the working fluid and includes two heat exchange and thermoacoustic structure assemblies which drive the liquid in a push-pull arrangement to cause the liquid metal to oscillate at a resonant acoustic frequency on the order of 1000 Hz. The engine is positioned in the field of a magnet and is oriented such that the liquid metal oscillates in a direction orthogonal to the field of the magnet, whereby an alternating electrical potential is generated in the liquid metal. Low-loss, low-inductance electrical conductors electrically connected to opposite sides of the liquid metal conduct an output signal to a transformer adapted to convert the low-voltage, high-current output signal to a more usable higher voltage, lower current signal.
AC magnetohydrodynamic microfluidic switch
Lemoff, A V; Lee, A P
2000-03-02
A microfluidic switch has been demonstrated using an AC Magnetohydrodynamic (MHD) pumping mechanism in which the Lorentz force is used to pump an electrolytic solution. By integrating two AC MHD pumps into different arms of a Y-shaped fluidic circuit, flow can be switched between the two arms. This type of switch can be used to produce complex fluidic routing, which may have multiple applications in {micro}TAS.
Collisional Drift Waves in Stellarator Plasmas
J.L.V. Lewandowski
2003-10-07
A computational study of resistive drift waves in the edge plasma of a stellarator with an helical magnetic axis is presented. Three coupled field equations, describing the collisional drift wave dynamics in the linear approximation, are solved as an initial-value problem along the magnetic field line. The magnetohydrodynamic equilibrium is obtained from a three-dimensional local equilibrium model. The use of a local magnetohydrodynamic equilibrium model allows for a computationally efficient systematic study of the impact of the magnetic field structure on drift wave stability.
Magnetohydrodynamics of fractal media
Tarasov, Vasily E.
2006-05-15
The fractal distribution of charged particles is considered. An example of this distribution is the charged particles that are distributed over the fractal. The fractional integrals are used to describe fractal distribution. These integrals are considered as approximations of integrals on fractals. Typical turbulent media could be of a fractal structure and the corresponding equations should be changed to include the fractal features of the media. The magnetohydrodynamics equations for fractal media are derived from the fractional generalization of integral Maxwell equations and integral hydrodynamics (balance) equations. Possible equilibrium states for these equations are considered.
NASA Astrophysics Data System (ADS)
Li, Zhilin; Xiao, Li; Cai, Qin; Zhao, Hongkai; Luo, Ray
2015-09-01
In this paper, a new Navier-Stokes solver based on a finite difference approximation is proposed to solve incompressible flows on irregular domains with open, traction, and free boundary conditions, which can be applied to simulations of fluid structure interaction, implicit solvent model for biomolecular applications and other free boundary or interface problems. For some problems of this type, the projection method and the augmented immersed interface method (IIM) do not work well or does not work at all. The proposed new Navier-Stokes solver is based on the local pressure boundary method, and a semi-implicit augmented IIM. A fast Poisson solver can be used in our algorithm which gives us the potential for developing fast overall solvers in the future. The time discretization is based on a second order multi-step method. Numerical tests with exact solutions are presented to validate the accuracy of the method. Application to fluid structure interaction between an incompressible fluid and a compressible gas bubble is also presented.
Teaching an Old Dog an Old Trick: FREE-FIX and Free-Boundary Axisymmetric MHD Equilibrium
NASA Astrophysics Data System (ADS)
Guazzotto, Luca
2015-11-01
A common task in plasma physics research is the calculation of an axisymmetric equilibrium for tokamak modeling. The main unknown of the problem is the magnetic poloidal flux ψ. The easiest approach is to assign the shape of the plasma and only solve the equilibrium problem in the plasma / closed-field-lines region (the ``fixed-boundary approach''). Often, one may also need the vacuum fields, i.e. the equilibrium in the open-field-lines region, requiring either coil currents or ψ on some closed curve outside the plasma to be assigned (the ``free-boundary approach''). Going from one approach to the other is a textbook problem, involving the calculation of Green's functions and surface integrals in the plasma. However, no tools are readily available to perform this task. Here we present a code (FREE-FIX) to compute a boundary condition for a free-boundary equilibrium given only the corresponding fixed-boundary equilibrium. An improvement to the standard solution method, allowing for much faster calculations, is presented. Applications are discussed. PPPL fund 245139 and DOE grant G00009102.
NASA Astrophysics Data System (ADS)
Ferraro, N. M.; Jardin, S. C.; Lao, L. L.; Shephard, M. S.; Zhang, F.
2016-05-01
Free-boundary 3D tokamak equilibria and resistive wall instabilities are calculated using a new resistive wall model in the two-fluid M3D-C1 code. In this model, the resistive wall and surrounding vacuum region are included within the computational domain. This implementation contrasts with the method typically used in fluid codes in which the resistive wall is treated as a boundary condition on the computational domain boundary and has the advantage of maintaining purely local coupling of mesh elements. This new capability is used to simulate perturbed, free-boundary non-axisymmetric equilibria; the linear evolution of resistive wall modes; and the linear and nonlinear evolution of axisymmetric vertical displacement events (VDEs). Calculated growth rates for a resistive wall mode with arbitrary wall thickness are shown to agree well with the analytic theory. Equilibrium and VDE calculations are performed in diverted tokamak geometry, at physically realistic values of dissipation, and with resistive walls of finite width. Simulations of a VDE disruption extend into the current-quench phase, in which the plasma becomes limited by the first wall, and strong currents are observed to flow in the wall, in the SOL, and from the plasma to the wall.
Li, Zhilin; Xiao, Li; Cai, Qin; Zhao, Hongkai; Luo, Ray
2016-01-01
In this paper, a new Navier–Stokes solver based on a finite difference approximation is proposed to solve incompressible flows on irregular domains with open, traction, and free boundary conditions, which can be applied to simulations of fluid structure interaction, implicit solvent model for biomolecular applications and other free boundary or interface problems. For some problems of this type, the projection method and the augmented immersed interface method (IIM) do not work well or does not work at all. The proposed new Navier–Stokes solver is based on the local pressure boundary method, and a semi-implicit augmented IIM. A fast Poisson solver can be used in our algorithm which gives us the potential for developing fast overall solvers in the future. The time discretization is based on a second order multi-step method. Numerical tests with exact solutions are presented to validate the accuracy of the method. Application to fluid structure interaction between an incompressible fluid and a compressible gas bubble is also presented. PMID:27087702
Spectrum of anomalous magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Giovannini, Massimo
2016-05-01
The equations of anomalous magnetohydrodynamics describe an Abelian plasma where conduction and chiral currents are simultaneously present and constrained by the second law of thermodynamics. At high frequencies the magnetic currents play the leading role, and the spectrum is dominated by two-fluid effects. The system behaves instead as a single fluid in the low-frequency regime where the vortical currents induce potentially large hypermagnetic fields. After deriving the physical solutions of the generalized Appleton-Hartree equation, the corresponding dispersion relations are scrutinized and compared with the results valid for cold plasmas. Hypermagnetic knots and fluid vortices can be concurrently present at very low frequencies and suggest a qualitatively different dynamics of the hydromagnetic nonlinearities.
Potential vorticity in magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Webb, G. M.; Mace, R. L.
2015-01-01
A version of Noether's second theorem using Lagrange multipliers is used to investigate fluid relabelling symmetries conservation laws in magnetohydrodynamics (MHD). We obtain a new generalized potential vorticity type conservation equation for MHD which takes into account entropy gradients and the J × B force on the plasma due to the current J and magnetic induction B. This new conservation law for MHD is derived by using Noether's second theorem in conjunction with a class of fluid relabelling symmetries in which the symmetry generator for the Lagrange label transformations is non-parallel to the magnetic field induction in Lagrange label space. This is associated with an Abelian Lie pseudo algebra and a foliated phase space in Lagrange label space. It contains as a special case Ertel's theorem in ideal fluid mechanics. An independent derivation shows that the new conservation law is also valid for more general physical situations.
Magnetohydrodynamic inertial reference system
NASA Astrophysics Data System (ADS)
Eckelkamp-Baker, Dan; Sebesta, Henry R.; Burkhard, Kevin
2000-07-01
Optical platforms increasingly require attitude knowledge and optical instrument pointing at sub-microradian accuracy. No low-cost commercial system exists to provide this level of accuracy for guidance, navigation, and control. The need for small, inexpensive inertial sensors, which may be employed in pointing control systems that are required to satisfy angular line-of-sight stabilization jitter error budgets to levels of 1-3 microradian rms and less, has existed for at least two decades. Innovations and evolutions in small, low-noise inertial angular motion sensor technology and advances in the applications of the global positioning system have converged to allow improvement in acquisition, tracking and pointing solutions for a wide variety of payloads. We are developing a small, inexpensive, and high-performance inertial attitude reference system that uses our innovative magnetohydrodynamic angular rate sensor technology.
Stellarator Coil Design and Plasma Sensitivity
Long-Poe Ku and Allen H. Boozer
2010-11-03
The rich information contained in the plasma response to external magnetic perturbations can be used to help design stellarator coils more effectively. We demonstrate the feasibility by first devel- oping a simple, direct method to study perturbations in stellarators that do not break stellarator symmetry and periodicity. The method applies a small perturbation to the plasma boundary and evaluates the resulting perturbed free-boundary equilibrium to build up a sensitivity matrix for the important physics attributes of the underlying configuration. Using this sensitivity information, design methods for better stellarator coils are then developed. The procedure and a proof-of-principle application are given that (1) determine the spatial distributions of external normal magnetic field at the location of the unperturbed plasma boundary to which the plasma properties are most sen- sitive, (2) determine the distributions of external normal magnetic field that can be produced most efficiently by distant coils, (3) choose the ratios of the magnitudes of the the efficiently produced magnetic distributions so the sensitive plasma properties can be controlled. Using these methods, sets of modular coils are found for the National Compact Stellarator Experiment (NCSX) that are either smoother or can be located much farther from the plasma boundary than those of the present design.
Lattice Boltzmann model for resistive relativistic magnetohydrodynamics.
Mohseni, F; Mendoza, M; Succi, S; Herrmann, H J
2015-08-01
In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfvén waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has a negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to σ-1/2, σ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere. PMID:26382548
Lattice Boltzmann model for resistive relativistic magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Mohseni, F.; Mendoza, M.; Succi, S.; Herrmann, H. J.
2015-08-01
In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfvén waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has a negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to σ-1 / 2, σ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere.
New Classes of Quasi-helically Symmetric Stellarators
L.P. Ku and A.H. Boozer
2010-08-09
New classes of quasi-helically symmetric stellarators with aspect ratios ≤ 10 have been found which are stable to the perturbation of magnetohydrodynamic modes at plasma pressures of practical interest. These configurations have large rotational transform and good quality of flux surfaces. Characteristics of some selected examples are discussed in detail. The feasibility of using modular coils for these stellarators has been investigated. It is shown that practical designs for modular coils can be achieved.
Conservation of circulation in magnetohydrodynamics
Bekenstein; Oron
2000-10-01
We demonstrate at both the Newtonian and (general) relativistic levels the existence of a generalization of Kelvin's circulation theorem (for pure fluids) that is applicable to perfect magnetohydrodynamics. The argument is based on the least action principle for magnetohydrodynamic flow. Examples of the new conservation law are furnished. The new theorem should be helpful in identifying new kinds of vortex phenomena distinct from magnetic ropes or fluid vortices. PMID:11089118
NASA Technical Reports Server (NTRS)
Linsky, J. L.
1980-01-01
Developments in the understanding and use of chromospheric diagnostics are discussed with emphasis on the following aspects: (1) trends emerging from semiempirical models of single stars; (2) the validity of claims that theoretical models of chromospheres are becoming realistic; (3) the correlation between the widths of Ca 2 H and K line emission cores and stellar absolute luminosity extending over 15 magnitudes (Wilson-Bappu relation); and (4) the existence of systematic flow patterns in stellar chromospheres.
NASA Astrophysics Data System (ADS)
Peletier, Reynier F.
2013-10-01
This is a summary of my lectures during the 2011 Canary Islands Winter School in Puerto de la Cruz. I give an introduction to the field of stellar populations in galaxies, and highlight some new results. Since the title of the Winter School is Secular Evolution in Galaxies I mostly concentrate on nearby galaxies, which are best suited to study this theme. Of course, the understanding of stellar populations is intimately connected to understanding the formation and evolution of galaxies, one of the great outstanding problems of astronomy. We are currently in a situation where very large observational advances have been made in recent years. Galaxies have been detected up to a redshift of ten. A huge effort has to be made so that stellar population theory can catch up with observations. Since most galaxies are far away, information about them has to come from stellar population synthesis of integrated light. Here I will discuss how stellar evolution theory, together with observations in our Milky Way and Local Group, are used as building blocks to analyse these integrated stellar populations.
Research into the process of impingement of two plane jets of an ideal fluid with free boundaries
NASA Astrophysics Data System (ADS)
Baskakov, V. D.; Karnaukhov, K. A.
2016-07-01
The problem of finding parameters stationary jets, outgoing from a place of impact of two incoming flat ideal jets with free boundaries and possessing the identical speed, but various width, has no decision. Various models are widely used for a conclusion of the missing equation now, but they lead to contradictory results. The new model is offered. Adequacy to the offered model was checked by comparison results with data of the numerical calculations in ANSYS AUTODYN. The dependence approximating results of numerical calculations is developed to increase accuracy in calculations of angular provision of the internal outgoing jet.First of all, the executed researches are interesting to experts, who works on behavior low- value technological errors in shaped charges.
Magnetohydrodynamic Augmented Propulsion Experiment
NASA Technical Reports Server (NTRS)
Litchford, Ron J.; Cole, John; Lineberry, John; Chapman, Jim; Schmidt, Harold; Cook, Stephen (Technical Monitor)
2002-01-01
A fundamental obstacle to routine space access is the specific energy limitations associated with chemical fuels. In the case of vertical take-off, the high thrust needed for vertical liftoff and acceleration to orbit translates into power levels in the 10 GW range. Furthermore, useful payload mass fractions are possible only if the exhaust particle energy (i.e., exhaust velocity) is much greater than that available with traditional chemical propulsion. The electronic binding energy released by the best chemical reactions (e.g., LOX/LH2 for example, is less than 2 eV per product molecule (approx. 1.8 eV per H2O molecule), which translates into particle velocities less than 5 km/s. Useful payload fractions, however, will require exhaust velocities exceeding 15 km/s (i.e., particle energies greater than 20 eV). As an added challenge, the envisioned hypothetical RLV (reusable launch vehicle) should accomplish these amazing performance feats while providing relatively low acceleration levels to orbit (2-3g maximum). From such fundamental considerations, it is painfully obvious that planned and current RLV solutions based on chemical fuels alone represent only a temporary solution and can only result in minor gains, at best. What is truly needed is a revolutionary approach that will dramatically reduce the amount of fuel and size of the launch vehicle. This implies the need for new compact high-power energy sources as well as advanced accelerator technologies for increasing engine exhaust velocity. Electromagnetic acceleration techniques are of immense interest since they can be used to circumvent the thermal limits associated with conventional propulsion systems. This paper describes the Magnetohydrodynamic Augmented Propulsion Experiment (MAPX) being undertaken at NASA Marshall Space Flight Center (MSFC). In this experiment, a 1-MW arc heater is being used as a feeder for a 1-MW magnetohydrodynamic (MHD) accelerator. The purpose of the experiment is to demonstrate
Filamentary magnetohydrodynamic plasmas
Kinney, R.; Tajima, T.; Petviashvili, N.; McWilliams, J.C.
1993-05-01
A filamentary construct of magnetohydrodynamical plasma dynamics, based on the Elsasser variables was developed. This approach is modeled after discrete vortex models of hydrodynamical turbulence, which cannot be expected in general to produce results identical to ones based on a Fourier decomposition of the fields. In a highly intermittent plasma, the induction force is small compared to the convective motion, and when this force is neglected. the plasma vortex system is described by a Hamiltonian. For a system with many such vortices we present a statistical treatment of a collection of discrete current-vorticity concentrations. Canonical and microcanonical statistical calculations show that both the vorticity and the current spectra are peaked at long wavelengths, and the expected states revert to known hydrodynamical states as the magnetic field vanishes. These results differ from previous Fourier-based statistical theories. but it is found that when the filament calculation is expanded to include the inductive force, the results approach the Fourier equilibria in the low-temperature limit, and the previous Hamiltonian plasma vortex results in the high-temperature limit. Numerical simulations of a large number of filaments are carried out and support the theory. A three-dimensional vortex model is outlined as well, which is also Hamiltonian when the inductive force is neglected.
Generalized reduced magnetohydrodynamic equations
Kruger, S.E.
1999-02-01
A new derivation of reduced magnetohydrodynamic (MHD) equations is presented. A multiple-time-scale expansion is employed. It has the advantage of clearly separating the three time scales of the problem associated with (1) MHD equilibrium, (2) fluctuations whose wave vector is aligned perpendicular to the magnetic field, and (3) those aligned parallel to the magnetic field. The derivation is carried out without relying on a large aspect ratio assumption; therefore this model can be applied to any general configuration. By accounting for the MHD equilibrium and constraints to eliminate the fast perpendicular waves, equations are derived to evolve scalar potential quantities on a time scale associated with the parallel wave vector (shear-Alfven wave time scale), which is the time scale of interest for MHD instability studies. Careful attention is given in the derivation to satisfy energy conservation and to have manifestly divergence-free magnetic fields to all orders in the expansion parameter. Additionally, neoclassical closures and equilibrium shear flow effects are easily accounted for in this model. Equations for the inner resistive layer are derived which reproduce the linear ideal and resistive stability criterion of Glasser, Greene, and Johnson. The equations have been programmed into a spectral initial value code and run with shear flow that is consistent with the equilibrium input into the code. Linear results of tearing modes with shear flow are presented which differentiate the effects of shear flow gradients in the layer with the effects of the shear flow decoupling multiple harmonics.
Multi-symplectic magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Webb, G. M.; McKenzie, J. F.; Zank, G. P.; Zank
2014-10-01
A multi-symplectic formulation of ideal magnetohydrodynamics (MHD) is developed based on the Clebsch variable variational principle in which the Lagrangian consists of the kinetic minus the potential energy of the MHD fluid modified by constraints using Lagrange multipliers that ensure mass conservation, entropy advection with the flow, the Lin constraint, and Faraday's equation (i.e. the magnetic flux is Lie dragged with the flow). The analysis is also carried out using the magnetic vector potential Ã where α=Ã. d x is Lie dragged with the flow, and B=∇×Ã. The multi-symplectic conservation laws give rise to the Eulerian momentum and energy conservation laws. The symplecticity or structural conservation laws for the multi-symplectic system corresponds to the conservation of phase space. It corresponds to taking derivatives of the momentum and energy conservation laws and combining them to produce n(n-1)/2 extra conservation laws, where n is the number of independent variables. Noether's theorem for the multi-symplectic MHD system is derived, including the case of non-Cartesian space coordinates, where the metric plays a role in the equations.
Filamentary magnetohydrodynamic plasmas
Kinney, R.; Tajima, T. ); McWilliams, J.C. ); Petviashvili, N. )
1994-02-01
A filamentary construct of magnetohydrodynamical plasma dynamics based on the Elsaesser variables is developed. This approach is modeled after discrete vortex models of hydrodynamical turbulence, which cannot be expected in general to produce results identical to those based on a Fourier decomposition of the fields. In a highly intermittent plasma, the induction force is small compared to the convective motion, and when this force is neglected, the plasma vortex system is described by a Hamiltonian. A statistical treatment of a collection of discrete current-vorticity concentrations is given. Canonical and microcanonical statistical calculations show that both the vorticity and the current spectra are peaked at long wavelengths, and the expected states revert to known hydrodynamical states as the magnetic field vanishes. These results differ from previous Fourier-based statistical theories, but it is found that when the filament calculation is expanded to include the inductive force, the results approach the Fourier equilibria in the low-temperature limit, and the previous Hamiltonian plasma vortex results in the high-temperature limit. Numerical simulations of a large number of filaments are carried out and support the theory. A three-dimensional vortex model is presented as well, which is also Hamiltonian when the inductive force is neglected. A statistical calculation in the canonical ensemble and numerical simulations show that a nonzero large-scale magnetic field is statistically favored, and that the preferred shape of this field is a long, thin tube of flux. Possible applications to a variety of physical phenomena are suggested.
Filamentary magnetohydrodynamic plasmas
NASA Astrophysics Data System (ADS)
Kinney, R.; Tajima, T.; Petviashvili, N.; McWilliams, J. C.
1993-05-01
A filamentary construct of magnetohydrodynamical plasma dynamics, based on the Elsasser variables was developed. This approach is modeled after discrete vortex models of hydrodynamical turbulence, which cannot be expected in general to produce results identical to ones based on a Fourier decomposition of the fields. In a highly intermittent plasma, the induction force is small compared to the convective motion, and when this force is neglected, the plasma vortex system is described by a Hamiltonian. For a system with many such vortices we present a statistical treatment of a collection of discrete current-vorticity concentrations. Canonical and microcanonical statistical calculations show that both the vorticity and the current spectra are peaked at long wavelengths, and the expected states revert to known hydrodynamical states as the magnetic field vanishes. These results differ from previous Fourier-based statistical theories, but it is found that when the filament calculation is expanded to include the inductive force, the results approach the Fourier equilibria in the low-temperature limit, and the previous Hamiltonian plasma vortex results in the high-temperature limit. Numerical simulations of a large number of filaments are carried out and support the theory. A three-dimensional vortex model is outlined as well, which is also Hamiltonian when the inductive force is neglected.
Solar Flares: Magnetohydrodynamic Processes
NASA Astrophysics Data System (ADS)
Shibata, Kazunari; Magara, Tetsuya
2011-12-01
This paper outlines the current understanding of solar flares, mainly focused on magnetohydrodynamic (MHD) processes responsible for producing a flare. Observations show that flares are one of the most explosive phenomena in the atmosphere of the Sun, releasing a huge amount of energy up to about 10^32 erg on the timescale of hours. Flares involve the heating of plasma, mass ejection, and particle acceleration that generates high-energy particles. The key physical processes for producing a flare are: the emergence of magnetic field from the solar interior to the solar atmosphere (flux emergence), local enhancement of electric current in the corona (formation of a current sheet), and rapid dissipation of electric current (magnetic reconnection) that causes shock heating, mass ejection, and particle acceleration. The evolution toward the onset of a flare is rather quasi-static when free energy is accumulated in the form of coronal electric current (field-aligned current, more precisely), while the dissipation of coronal current proceeds rapidly, producing various dynamic events that affect lower atmospheres such as the chromosphere and photosphere. Flares manifest such rapid dissipation of coronal current, and their theoretical modeling has been developed in accordance with observations, in which numerical simulations proved to be a strong tool reproducing the time-dependent, nonlinear evolution of a flare. We review the models proposed to explain the physical mechanism of flares, giving an comprehensive explanation of the key processes mentioned above. We start with basic properties of flares, then go into the details of energy build-up, release and transport in flares where magnetic reconnection works as the central engine to produce a flare.
Compressible magnetohydrodynamic sawtooth crash
NASA Astrophysics Data System (ADS)
Sugiyama, Linda E.
2014-02-01
In a toroidal magnetically confined plasma at low resistivity, compressible magnetohydrodynamic (MHD) predicts that an m = 1/n = 1 sawtooth has a fast, explosive crash phase with abrupt onset, rate nearly independent of resistivity, and localized temperature redistribution similar to experimental observations. Large scale numerical simulations show that the 1/1 MHD internal kink grows exponentially at a resistive rate until a critical amplitude, when the plasma motion accelerates rapidly, culminating in fast loss of the temperature and magnetic structure inside q < 1, with somewhat slower density redistribution. Nonlinearly, for small effective growth rate the perpendicular momentum rate of change remains small compared to its individual terms ∇p and J × B until the fast crash, so that the compressible growth rate is determined by higher order terms in a large aspect ratio expansion, as in the linear eigenmode. Reduced MHD fails completely to describe the toroidal mode; no Sweet-Parker-like reconnection layer develops. Important differences result from toroidal mode coupling effects. A set of large aspect ratio compressible MHD equations shows that the large aspect ratio expansion also breaks down in typical tokamaks with rq =1/Ro≃1/10 and a /Ro≃1/3. In the large aspect ratio limit, failure extends down to much smaller inverse aspect ratio, at growth rate scalings γ =O(ɛ2). Higher order aspect ratio terms, including B˜ϕ, become important. Nonlinearly, higher toroidal harmonics develop faster and to a greater degree than for large aspect ratio and help to accelerate the fast crash. The perpendicular momentum property applies to other transverse MHD instabilities, including m ≥ 2 magnetic islands and the plasma edge.
Filamentary magnetohydrodynamic plasmas
NASA Astrophysics Data System (ADS)
Kinney, R.; Tajima, T.; McWilliams, J. C.; Petviashvili, N.
1994-02-01
A filamentary construct of magnetohydrodynamical plasma dynamics based on the Elsässer variables is developed. This approach is modeled after discrete vortex models of hydrodynamical turbulence, which cannot be expected in general to produce results identical to those based on a Fourier decomposition of the fields. In a highly intermittent plasma, the induction force is small compared to the convective motion, and when this force is neglected, the plasma vortex system is described by a Hamiltonian. A statistical treatment of a collection of discrete current-vorticity concentrations is given. Canonical and microcanonical statistical calculations show that both the vorticity and the current spectra are peaked at long wavelengths, and the expected states revert to known hydrodynamical states as the magnetic field vanishes. These results differ from previous Fourier-based statistical theories, but it is found that when the filament calculation is expanded to include the inductive force, the results approach the Fourier equilibria in the low-temperature limit, and the previous Hamiltonian plasma vortex results in the high-temperature limit. Numerical simulations of a large number of filaments are carried out and support the theory. A three-dimensional vortex model is presented as well, which is also Hamiltonian when the inductive force is neglected. A statistical calculation in the canonical ensemble and numerical simulations show that a nonzero large-scale magnetic field is statistically favored, and that the preferred shape of this field is a long, thin tube of flux. Possible applications to a variety of physical phenomena are suggested.
NASA Astrophysics Data System (ADS)
Owocki, Stan
A "stellar wind" is the continuous, supersonic outflow of matter from the surface layers of a star. Our sun has a solar wind, driven by the gas-pressure expansion of the hot (T > 106 K) solar corona. It can be studied through direct in situ measurement by interplanetary spacecraft; but analogous coronal winds in more distant solar-type stars are so tenuous and transparent that that they are difficult to detect directly. Many more luminous stars have winds that are dense enough to be opaque at certain wavelengths of the star's radiation, making it possible to study their wind outflows remotely through careful interpretation of the observed stellar spectra. Red giant stars show slow, dense winds that may be driven by the pressure from magnetohydrodyanmic waves. As stars with initial mass up to 8 M ⊙ evolve toward the Asymptotic Giant Branch (AGB), a combination of stellar pulsations and radiative scattering off dust can culminate in "superwinds" that strip away the entire stellar envelope, leaving behind a hot white dwarf stellar core with less than the Chandrasekhar mass of ˜ 1. 4M ⊙. The winds of hot, luminous, massive stars are driven by line-scattering of stellar radiation, but such massive stars can also exhibit superwind episodes, either as Red Supergiants or Luminous Blue Variable stars. The combined wind and superwind mass loss can strip the star's hydrogen envelope, leaving behind a Wolf-Rayet star composed of the products of earlier nuclear burning via the CNO cycle. In addition to such direct effects on a star's own evolution, stellar winds can be a substantial source of mass, momentum, and energy to the interstellar medium, blowing open large cavities or "bubbles" in this ISM, seeding it with nuclear processed material, and even helping trigger the formation of new stars, and influencing their eventual fate as white dwarves or core-collapse supernovae. This chapter reviews the properties of such stellar winds, with an emphasis on the various
An approximate single fluid 3-dimensional magnetohydrodynamic equilibrium model with toroidal flow
NASA Astrophysics Data System (ADS)
Cooper, W. A.; Hirshman, S. P.; Chapman, I. T.; Brunetti, D.; Faustin, J. M.; Graves, J. P.; Pfefferlé, D.; Raghunathan, M.; Sauter, O.; Tran, T. M.; Aiba, N.
2014-09-01
An approximate model for a single fluid three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium with pure isothermal toroidal flow with imposed nested magnetic flux surfaces is proposed. It recovers the rigorous toroidal rotation equilibrium description in the axisymmetric limit. The approximation is valid under conditions of nearly rigid or vanishing toroidal rotation in regions with significant 3D deformation of the equilibrium flux surfaces. Bifurcated helical core equilibrium simulations of long-lived modes in the MAST device demonstrate that the magnetic structure is only weakly affected by the flow but that the 3D pressure distortion is important. The pressure is displaced away from the major axis and therefore is not as noticeably helically deformed as the toroidal magnetic flux under the subsonic flow conditions measured in the experiment. The model invoked fails to predict any significant screening by toroidal plasma rotation of resonant magnetic perturbations in MAST free boundary computations.
NASA Astrophysics Data System (ADS)
Cooper, W. A.; Brunetti, D.; Duval, B. P.; Faustin, J. M.; Graves, J. P.; Kleiner, A.; Patten, H.; Pfefferlé, D.; Porte, L.; Raghunathan, M.; Reimerdes, H.; Sauter, O.; Tran, T. M.
2016-04-01
Free boundary magnetohydrodynamic equilibrium states with spontaneous three dimensional deformations of the plasma-vacuum interface are computed for the first time. The structures obtained have the appearance of saturated ideal external kink/peeling modes. High edge pressure gradients yield toroidal mode number n = 1 corrugations for a high edge bootstrap current and larger n distortions when this current is small. Deformations in the plasma boundary region induce a nonaxisymmetric Pfirsch-Schlüter current driving a field-aligned current ribbon consistent with reported experimental observations. A variation in the 3D equilibrium confirms that the n = 1 mode is a kink/peeling structure. We surmise that our calculated equilibrium structures constitute a viable model for the edge harmonic oscillations and outer modes associated with a quiescent H-mode operation in shaped tokamak plasmas.
NASA Technical Reports Server (NTRS)
Carpenter, Kenneth G.; Schrijver, Carolus J.; Karovska, Margarita
2006-01-01
The ultra-sharp images of the Stellar Imager (SI) will revolutionize our view of many dynamic astrophysical processes: The 0.1 milliarcsec resolution of this deep-space telescope will transform point sources into extended sources, and simple snapshots into spellbinding evolving views. SI s science focuses on the role of magnetism in the Universe, particularly on magnetic activity on the surfaces of stars like the Sun. SI s prime goal is to enable long-term forecasting of solar activity and the space weather that it drives in support of the Living With a Star program in the Exploration Era by imaging a sample of magnetically active stars with enough resolution to map their evolving dynamo patterns and their internal flows. By exploring the Universe at ultra-high resolution, SI will also revolutionize our understanding of the formation of planetary systems, of the habitability and climatology of distant planets, and of many magnetohydrodynamically controlled structures and processes in the Universe.
[Nonlinear magnetohydrodynamics]. Final report
Montgomery, D.C.
1998-11-01
This is a final report on the research activities carried out under the above grant at Dartmouth. During the period considered, the grant was identified as being for nonlinear magnetohydrodynamics, considered as the most tractable theoretical framework in which the plasma problems associated with magnetic confinement of fusion plasmas could be studied. During the first part of the grant`s lifetime, the author was associated with Los Alamos National Laboratory as a consultant and the work was motivated by the reversed-field pinch. Later, when that program was killed at Los Alamos, the problems became ones that could be motivated by their relation to tokamaks. Throughout the work, the interest was always on questions that were as fundamental as possible, compatible with those motivations. The intent was always to contribute to plasma physics as a science, as well as to the understanding of mission-oriented confined fusion plasmas. Twelve Ph.D. theses were supervised during this period and a comparable number of postdoctoral research associates were temporarily supported. Many of these have gone on to distinguished careers, though few have done so in the context of the controlled fusion program. Their work was a combination of theory and numerical computation, in gradually less and less idealized settings, moving from rectangular periodic boundary conditions in two dimensions, through periodic straight cylinders and eventually, before the grant was withdrawn, to toroids, with a gradually more prominent role for electrical and mechanical boundary conditions. The author never had access to a situation where he could initiate experiments and relate directly to the laboratory data he wanted. Computers were the laboratory. Most of the work was reported in referred publications in the open literature, copies of which were transmitted one by one to DOE at the time they appeared. The Appendix to this report is a bibliography of published work which was carried out under the
Pulsar Magnetohydrodynamic Winds
NASA Astrophysics Data System (ADS)
Okamoto, Isao; Sigalo, Friday B.
2006-12-01
The acceleration and collimation/decollimation of relativistic magnetocentrifugal winds are discussed concerning a cold plasma from a strongly magnetized, rapidly rotating neutron star in a steady axisymmetric state based on ideal magnetohydrodynamics. There exist unipolar inductors associated with the field line angular frequency, α, at the magnetospheric base surface, SB, with a huge potential difference between the poles and the equator, which drive electric current through the pulsar magnetosphere. Any ``current line'' must emanate from one terminal of the unipolar inductor and return to the other, converting the Poynting flux to the kinetic flux of the wind at finite distances. In a plausible field structure satisfying the transfield force-balance equation, the fast surface, SF, must exist somewhere between the subasymptotic and asymptotic domains, i.e., at the innermost point along each field line of the asymptotic domain of \\varpaA2/\\varpi2 ≪ 1, where \\varpiA is the Alfvénic axial distance. The criticality condition at SF yields the Lorentz factor, γF = μ\\varepsilon1/3, and the angular momentum flux, β, as the eigenvalues in terms of the field line angular velocity, α, the mass flux per unit flux tube, η, and one of the Bernoulli integrals, μδ, which are assumed to be specifiable as the boundary conditions at SB. The other Bernoulli integral, μɛ, is related to μδ as μɛ = μδ[1-(α2\\varpiA2/c2)]-1, and both μɛ and \\varpiA2 are eigenvalues to be determined by the criticality condition at SF. Ongoing MHD acceleration is possible in the superfast domain. This fact may be helpful in resolving a discrepancy between the wind theory and the Crab-nebula model. It is argued that the ``anti-collimation theorem'' holds for relativistic winds, based on the curvature of field streamlines determined by the transfield force balance. The ``theorem'' combines with the ``current-closure condition'' as a global condition in the wind zone to produce a
NASA Technical Reports Server (NTRS)
Parker, E. N.
1986-01-01
The present state of development of the theory of coronal heating is summarized. Coronal heating is the general cause of stellar X-ray emission, and it is also the cause of stellar mass loss in most stars. Hence a quantitive theory of coronal heating is an essential part of X-ray astronomy, and the development of a correct theory of coronal heating should be a primary concern of X-ray astronomers. The magnetohydrodynamical effects involved in coronal heating are not without interest in their own right, representing phenomena largely unknown in the terrestrial laboratory. Until these effects can be evaluated and assembled into a comprehensive theory of coronal heating for at least one star, the interpretation of the X-ray emissions of all stars is a phenomenological study at best, based on arbitrary organization and display of X-ray luminosity against bolometric luminosity, rotation rate, etc. The sun provides the one opportunity to pursue the exotic physical effects that combine to heat a stellar corona.
Reply to "Comment on 'Casimir force in the O(n→∞) model with free boundary conditions' ".
Dantchev, Daniel; Bergknoff, Jonathan; Rudnick, Joseph
2015-02-01
The preceding Comment raises a few points concerning our paper [Phys. Rev. E 89, 042116 (2014)]. In this Reply we stress that although Diehl et al. [Europhys. Lett. 100, 10004 (2012) and Phys. Rev. E 89, 062123 (2014)] use three different models to study the Casimir force for the O(n→∞) model with free boundary conditions we study a single model over the entire range of temperatures from above the bulk critical temperature T(c) to absolute temperatures down to T=0. The use of a single model renders more transparent the crossover from effects dominated by critical fluctuations in the vicinity of the bulk transition temperature to effects controlled by Goldstone modes at low temperatures. Contrary to the assertion in the Comment, we make no claim for the superiority of our model over any of those considered by Diehl et al. [Europhys. Lett. 100, 10004 (2012) and Phys. Rev. E 89, 062123 (2014)]. We also present additional evidence supporting our conclusion in Dantchev et al. [Phys. Rev. E 89, 042116 (2014)] that the temperature range in which our low-temperature analytical expansion for the Casimir force increases as L grows and remains accurate for values of the ratio T/T(c) that become closer and closer to unity, whereas T remains well outside of the critical region. PMID:25768643
COUNTER-ROTATION IN RELATIVISTIC MAGNETOHYDRODYNAMIC JETS
Cayatte, V.; Sauty, C.; Vlahakis, N.; Tsinganos, K.; Matsakos, T.; Lima, J. J. G.
2014-06-10
Young stellar object observations suggest that some jets rotate in the opposite direction with respect to their disk. In a recent study, Sauty et al. showed that this does not contradict the magnetocentrifugal mechanism that is believed to launch such outflows. Motion signatures that are transverse to the jet axis, in two opposite directions, have recently been measured in M87. One possible interpretation of this motion is that of counter-rotating knots. Here, we extend our previous analytical derivation of counter-rotation to relativistic jets, demonstrating that counter-rotation can indeed take place under rather general conditions. We show that both the magnetic field and a non-negligible enthalpy are necessary at the origin of counter-rotating outflows, and that the effect is associated with a transfer of energy flux from the matter to the electromagnetic field. This can be realized in three cases: if a decreasing enthalpy causes an increase of the Poynting flux, if the flow decelerates, or if strong gradients of the magnetic field are present. An illustration of the involved mechanism is given by an example of a relativistic magnetohydrodynamic jet simulation.
Realistic magnetohydrodynamical simulation of solar local supergranulation
NASA Astrophysics Data System (ADS)
Ustyugov, Sergey D.
2010-12-01
Three-dimensional numerical simulations of solar surface magnetoconvection using realistic model physics are conducted. The thermal structure of convective motions into the upper radiative layers of the photosphere, the main scales of convective cells and the penetration depths of convection are investigated. We take part of the solar photosphere with a size of 60×60 Mm2 in the horizontal direction and of depth 20 Mm from the level of the visible solar surface. We use a realistic initial model of the sun and apply the equation of state and opacities of stellar matter. The equations of fully compressible radiation magnetohydrodynamics (MHD) with dynamical viscosity and gravity are solved. We apply (i) the conservative total variation diminishing (TVD) difference scheme for MHD, (ii) the diffusion approximation for radiative transfer and (iii) dynamical viscosity from subgrid-scale modeling. In simulation, we take a uniform two-dimensional grid in the horizontal plane and a nonuniform grid in the vertical direction with the number of cells being 600×600×204. We use 512 processors with distributed memory multiprocessors on the supercomputer MVS-100k at the Joint Computational Centre of the Russian Academy of Sciences.
Nucleosynthesis in Magnetohydrodynamical Jets from Collapsars
Ono, M.; Hashimoto, M.; Fujimoto, S.; Kotake, K.
2011-10-28
We investigate the heavy-element nucleosynthesis of a massive star whose mass in the main sequence stage is M{sub ms} = 70 M{sub {center_dot}}. Detailed calculations of the nucleosynthesis are performed during the hydrostatic stellar evolution until the core composed of iron-group nuclei begins to collapse. As a supernova explosion model, a collapsar model is constructed whose jets are driven by magnetohydrodynamical effects of a differentially rotating core. The heavy-element nucleosynthesis inside the jet of a collapsar model is followed along the trajectories of stream lines of the jet. We combine the results of both hydrostatic and heavy-element nucleosyntheses to compare with the solar abundances. We find that neutron-rich elements of 70140.
Magnetohydrodynamic flow at microelectrodes
NASA Astrophysics Data System (ADS)
Ragsdale, Steven Ronald
1998-12-01
Voltammetric reduction of nitrobenzene (NB) at a 12.5 μm-radius Pt microdisk electrode in acetonitrile solutions containing 0.001/le x NB/le 0.999 is reported (x NB is the mole fraction of NB). The voltammetric response displays a reversible, sigmoidalshape wave, corresponding to the one-electron reduction of NB. The maximum limiting current occurs in solutions containing intermediate redox concentrations, x NB/le0.2. Voltammetric currents are analyzed using the Cullinan-Vignes model to describe the interdiffusion of the redox species and solvent. Mutual diffusivities are corrected for activity effects using isothermal liquid-vapor equilibrium data. Application of the activity-corrected diffusivities in the Cullinan- Vignes model yields reasonably accurate predictions of the dependence of the voltammetric current on solution composition. The influence of an external magnetic field (0-1 Tesla) on the voltammetric response of Pt and Au microdisk electrodes (0.1, 6.4, 12.5 and 25 μm radius) is described. Magnetohydrodynamic (MHD) flow within a microscopic volume element adjacent to the microdisk surface results from the magnetic force generated by the flux of electrogenerated ions through the magnetic field. An analytic expression is presented for the magnetic force generated during steady-state voltammetry at a hemispherical microelectrode immersed in a uniform magnetic field. The magnetic volume force, F/bf mag (N/m3), is shown to decrease as r-2 (where r is the distance from the center of the electrode). The dependence of F/bf mag on r-2 confines the MHD flow to small volumes very close to the electrode surface (e.g., ~2×10-9 L for a 12.5 μm-radius electrode). Scanning electrochemical microscopy (SECM) is used to map MHD flows at a 25 μm-radius Pt microdisk electrode during the one-electron reduction of NB. Unidirectional lateral flow is observed when the magnetic field is aligned parallel to the electrode surface; rotational or cyclotron flow is observed when
Reply to "Comment on `Casimir force in the O (n →∞ ) model with free boundary conditions' "
NASA Astrophysics Data System (ADS)
Dantchev, Daniel; Bergknoff, Jonathan; Rudnick, Joseph
2015-02-01
The preceding Comment raises a few points concerning our paper [Phys. Rev. E 89, 042116 (2014), 10.1103/PhysRevE.89.042116]. In this Reply we stress that although Diehl et al. [Europhys. Lett. 100, 10004 (2012), 10.1209/0295-5075/100/10004 and Phys. Rev. E 89, 062123 (2014), 10.1103/PhysRevE.89.062123] use three different models to study the Casimir force for the O (n →∞ ) model with free boundary conditions we study a single model over the entire range of temperatures from above the bulk critical temperature Tc to absolute temperatures down to T =0 . The use of a single model renders more transparent the crossover from effects dominated by critical fluctuations in the vicinity of the bulk transition temperature to effects controlled by Goldstone modes at low temperatures. Contrary to the assertion in the Comment, we make no claim for the superiority of our model over any of those considered by Diehl et al. [Europhys. Lett. 100, 10004 (2012), 10.1209/0295-5075/100/10004 and Phys. Rev. E 89, 062123 (2014), 10.1103/PhysRevE.89.062123]. We also present additional evidence supporting our conclusion in Dantchev et al. [Phys. Rev. E 89, 042116 (2014), 10.1103/PhysRevE.89.042116] that the temperature range in which our low-temperature analytical expansion for the Casimir force increases as L grows and remains accurate for values of the ratio T /Tc that become closer and closer to unity, whereas T remains well outside of the critical region.
Representation of Ideal Magnetohydrodynamic Modes
Roscoe B. White
2013-01-15
One of the most fundamental properties of ideal magnetohydrodynamics is the condition that plasma motion cannot change magnetic topology. The conventional representation of ideal magnetohydrodynamic modes by perturbing a toroidal equilibrium field through δ Β = ∇ X (xi X B) ensures that δ B • ∇ ψ = 0 at a resonance, with ψ labelling an equilibrium flux surface. Also useful for the analysis of guiding center orbits in a perturbed field is the representation δ Β = ∇ X αB. These two representations are equivalent, but the vanishing of δ B • ∇ψ at a resonance is necessary but not sufficient for the preservation of field line topology, and a indiscriminate use of either perturbation in fact destroys the original equilibrium flux topology. It is necessary to find the perturbed field to all orders in xi to conserve the original topology. The effect of using linearized perturbations on stability and growth rate calculations is discussed
BOOK REVIEW: Nonlinear Magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Shafranov, V.
1998-08-01
Nonlinear magnetohydrodynamics by Dieter Biskamp is a thorough introduction to the physics of the most impressive non-linear phenomena that occur in conducting magnetoplasmas. The basic systems, in which non-trivial dynamic processes are observed, accompanied by changes of geometry of the magnetic field and the effects of energy transformation (magnetic energy into kinetic energy or the opposite effect in magnetic dynamos), are the plasma magnetic confinement systems for nuclear fusion and space plasmas, mainly the solar plasma. A significant number of the examples of the dynamic processes considered are taken from laboratory plasmas, for which an experimental check of the theory is possible. Therefore, though the book is intended for researchers and students interested in both laboratory, including nuclear fusion, and astrophysical plasmas, it is most probably closer to the first category of reader. In the Introduction the author notes that unlike the hydrodynamics of non-conducting fluids, where the phenomena caused by rapid fluid motions are the most interesting, for plasmas in a strong magnetic field the quasi-static configurations inside which the local dynamic processes occur are often the most important. Therefore, the reader will also find in this book rather traditional material on the theory of plasma equilibrium and stability in magnetic fields. In addition, it is notable that, as opposed to a linear theory, the non-linear theory, as a rule, cannot give quite definite explanations or predictions of phenomena, and consequently there are in the book many results obtained by consideration of numerical models with the use of supercomputers. The treatment of non-linear dynamics is preceded by Chapters 2 to 4, in which the basics of MHD theory are presented with an emphasis on the role of integral invariants of the magnetic helicity type, a derivation of the reduced MHD equations is given, together with examples of the exact solutions of the equilibrium
MHD stability of the MHH2 stellarator
Garabedian, P.R.
1998-12-31
The NSTAB code provides a computer implementation of the variational principle of magnetohydrodynamics. Excellent resolution is obtained by combining a spectral representation in the toroidal and poloidal angles with a low order, but exceptionally accurate, finite difference scheme in the radial direction. Conservation form of the magnetostatics equations is used to capture islands and current sheets effectively on crude grids. This method enables one to discuss global stability by analyzing bifurcated solutions of the equilibrium problem. The author applies it to investigate the physics of the MHH2 stellarator, whose magnetic structure has a remarkable property of quasi-axial symmetry.
NASA Technical Reports Server (NTRS)
Tsiveriotis, K.; Brown, R. A.
1993-01-01
A new method is presented for the solution of free-boundary problems using Lagrangian finite element approximations defined on locally refined grids. The formulation allows for direct transition from coarse to fine grids without introducing non-conforming basis functions. The calculation of elemental stiffness matrices and residual vectors are unaffected by changes in the refinement level, which are accounted for in the loading of elemental data to the global stiffness matrix and residual vector. This technique for local mesh refinement is combined with recently developed mapping methods and Newton's method to form an efficient algorithm for the solution of free-boundary problems, as demonstrated here by sample calculations of cellular interfacial microstructure during directional solidification of a binary alloy.
JET FORMATION FROM MASSIVE YOUNG STARS: MAGNETOHYDRODYNAMICS VERSUS RADIATION PRESSURE
Vaidya, Bhargav; Porth, Oliver; Fendt, Christian; Beuther, Henrik E-mail: fendt@mpia.de
2011-11-20
Observations indicate that outflows from massive young stars are more collimated during their early evolution compared to later stages. Our paper investigates various physical processes that impact the outflow dynamics, i.e., its acceleration and collimation. We perform axisymmetric magnetohydrodynamic (MHD) simulations particularly considering the radiation pressure exerted by the star and the disk. We have modified the PLUTO code to include radiative forces in the line-driving approximation. We launch the outflow from the innermost disk region (r < 50 AU) by magnetocentrifugal acceleration. In order to disentangle MHD effects from radiative forces, we start the simulation in pure MHD and later switch on the radiation force. We perform a parameter study considering different stellar masses (thus luminosity), magnetic flux, and line-force strength. For our reference simulation-assuming a 30 M{sub Sun} star-we find substantial de-collimation of 35% due to radiation forces. The opening angle increases from 20 Degree-Sign to 32 Degree-Sign for stellar masses from 20 M{sub Sun} to 60 M{sub Sun }. A small change in the line-force parameter {alpha} from 0.60 to 0.55 changes the opening angle by {approx}8 Degree-Sign . We find that it is mainly the stellar radiation that affects the jet dynamics. Unless the disk extends very close to the star, its force is too small to have much impact. Essentially, our parameter runs with different stellar masses can be understood as a proxy for the time evolution of the star-outflow system. Thus, we have shown that when the stellar mass (thus luminosity) increases with age, the outflows become less collimated.
Compact stellarators with modular coils
Garabedian, P. R.
2000-01-01
Compact stellarator designs with modular coils and only two or three field periods are now available; these designs have both good stability and quasiaxial symmetry providing adequate transport for a magnetic fusion reactor. If the bootstrap current assumes theoretically predicted values a three field period configuration is optimal, but if that net current turns out to be lower, a device with two periods and just 12 modular coils might be better. There are also attractive designs with quasihelical symmetry and four or five periods whose properties depend less on the bootstrap current. Good performance requires that there be a satisfactory magnetic well in the vacuum field, which is a property lacking in a stellarator-tokamak hybrid that has been proposed for a proof of principle experiment. In this paper, we present an analysis of stability for these configurations that is based on a mountain pass theorem asserting that, if two solutions of the problem of magnetohydrodynamic equilibrium can be found, then there has to be an unstable solution. We compare results of our theory of equilibrium, stability, and transport with recently announced measurements from the large LHD experiment in Japan. PMID:10899993
Stellar Ablation of Planetary Atmospheres
NASA Technical Reports Server (NTRS)
Moore, Thomas E.; Horwitz, J. L.
2007-01-01
We review observations and theories of the solar ablation of planetary atmospheres, focusing on the terrestrial case where a large magnetosphere holds off the solar wind, so that there is little direct atmospheric impact, but also couples the solar wind electromagnetically to the auroral zones. We consider the photothermal escape flows known as the polar wind or refilling flows, the enhanced mass flux escape flows that result from localized solar wind energy dissipation in the auroral zones, and the resultant enhanced neutral atom escape flows. We term these latter two escape flows the "auroral wind." We review observations and theories of the heating and acceleration of auroral winds, including energy inputs from precipitating particles, electromagnetic energy flux at magnetohydrodynamic and plasma wave frequencies, and acceleration by parallel electric fields and by convection pickup processes also known as "centrifugal acceleration." We consider also the global circulation of ionospheric plasmas within the magnetosphere, their participation in magnetospheric disturbances as absorbers of momentum and energy, and their ultimate loss from the magnetosphere into the downstream solar wind, loading reconnection processes that occur at high altitudes near the magnetospheric boundaries. We consider the role of planetary magnetization and the accumulating evidence of stellar ablation of extrasolar planetary atmospheres. Finally, we suggest and discuss future needs for both the theory and observation of the planetary ionospheres and their role in solar wind interactions, to achieve the generality required for a predictive science of the coupling of stellar and planetary atmospheres over the full range of possible conditions.
Post, R.F.
1962-09-01
A method and means are described for injecting energetic neutral atoms or molecular ions into dense magnetically collimated plasma columns of stellarators and the like in such a manner that the atoms or ions are able to significantly penetrate the column before being ionized by collision with the plasma constituent particles. Penetration of the plasma column by the neutral atoms or molecular ions is facilitated by superposition of two closely spaced magnetic mirrors on the plasma confinement field. The mirrors are moved apart to magnetically sweep plasma from a region between the mirrors and establish a relatively low plasma density therein. By virture of the low density, neutral atoms or molecular ions injected into the region significantly penetrate the plasma column before being ionized. Thereafter, the mirrors are diminished to permit the injected material to admix with the plasma in the remainder of the column. (AEC)
Introduction to stellar evolution
NASA Astrophysics Data System (ADS)
Scilla, Degl’Innocenti
2016-04-01
This contribution is meant as a first brief introduction to stellar physics. First I shortly describe the main physical processes active in stellar structures then I summarize the most important features during the stellar life-cycle.
Stellar and galactic jets - Theoretical issues
NASA Technical Reports Server (NTRS)
Konigl, A.
1986-01-01
Theoretical issues pertaining to the modelling of jets in young stellar objects and in active galactic nuclei are reviewed. The strong morphological similarities between these two types of sources are emphasized, and observational constraints on the basic physical mechanisms that may be responsible for the jet phenomenon are outlined. Particular attention is given to the 'momentum-discharge problem' in molecular-cloud outflows and to its possible resolution in terms of a centrifugally driven magnetohydrodynamic wind from an accretion disk. In addition, various propagation effects are discussed, and the relevance to stellar jets of the de Laval collimation mechanism and of the accelerated-clump model for emission knots is assessed. The review concludes with a brief list of potentially useful observational tests.
Direct comparisons of compressible magnetohydrodynamics and reduced magnetohydrodynamics turbulence
NASA Astrophysics Data System (ADS)
Dmitruk, Pablo; Matthaeus, William H.; Oughton, Sean
2005-11-01
Direct numerical simulations of low Mach number compressible three-dimensional magnetohydrodynamic (CMHD3D) turbulence in the presence of a strong mean magnetic field are compared with simulations of reduced magnetohydrodynamics (RMHD). Periodic boundary conditions in the three spatial coordinates are considered. Different sets of initial conditions are chosen to explore the applicability of RMHD and to study how close the solution remains to the full compressible MHD solution as both freely evolve in time. In a first set, the initial state is prepared to satisfy the conditions assumed in the derivation of RMHD, namely, a strong mean magnetic field and plane-polarized fluctuations, varying weakly along the mean magnetic field. In those circumstances, simulations show that RMHD and CMHD3D evolve almost indistinguishably from one another. When some of the conditions are relaxed the agreement worsens but RMHD remains fairly close to CMHD3D, especially when the mean magnetic field is large enough. Moreover, the well-known spectral anisotropy effect promotes the dynamical attainment of the conditions for RMHD applicability. Global quantities (mean energies, mean-square current, and vorticity) and energy spectra from the two solutions are compared and point-to-point separation estimations are computed. The specific results shown here give support to the use of RMHD as a valid approximation of compressible MHD with a mean magnetic field under certain but quite practical conditions.
Method for manufacturing magnetohydrodynamic electrodes
Killpatrick, D.H.; Thresh, H.R.
1980-06-24
A method of manufacturing electrodes for use in a magnetohydrodynamic (MHD) generator is described comprising the steps of preparing a billet having a core of a first metal, a tubular sleeve of a second metal, and an outer sheath of an extrusile metal; evacuating the space between the parts of the assembled billet; extruding the billet; and removing the outer jacket. The extruded bar may be made into electrodes by cutting and bending to the shape required for an MHD channel frame. The method forms a bond between the first metal of the core and the second metal of the sleeve strong enough to withstand a hot and corrosive environment.
Adaptive wavelets and relativistic magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Hirschmann, Eric; Neilsen, David; Anderson, Matthe; Debuhr, Jackson; Zhang, Bo
2016-03-01
We present a method for integrating the relativistic magnetohydrodynamics equations using iterated interpolating wavelets. Such provide an adaptive implementation for simulations in multidimensions. A measure of the local approximation error for the solution is provided by the wavelet coefficients. They place collocation points in locations naturally adapted to the flow while providing expected conservation. We present demanding 1D and 2D tests includingthe Kelvin-Helmholtz instability and the Rayleigh-Taylor instability. Finally, we consider an outgoing blast wave that models a GRB outflow.
Magneto-Hydrodynamics Based Microfluidics
Qian, Shizhi; Bau, Haim H.
2009-01-01
In microfluidic devices, it is necessary to propel samples and reagents from one part of the device to another, stir fluids, and detect the presence of chemical and biological targets. Given the small size of these devices, the above tasks are far from trivial. Magnetohydrodynamics (MHD) offers an elegant means to control fluid flow in microdevices without a need for mechanical components. In this paper, we review the theory of MHD for low conductivity fluids and describe various applications of MHD such as fluid pumping, flow control in fluidic networks, fluid stirring and mixing, circular liquid chromatography, thermal reactors, and microcoolers. PMID:20046890
Non-axisymmetric equilibrium reconstruction for stellarators, reversed field pinches and tokamaks
NASA Astrophysics Data System (ADS)
Hanson, J. D.; Anderson, D. T.; Cianciosa, M.; Franz, P.; Harris, J. H.; Hartwell, G. H.; Hirshman, S. P.; Knowlton, S. F.; Lao, L. L.; Lazarus, E. A.; Marrelli, L.; Maurer, D. A.; Schmitt, J. C.; Sontag, A. C.; Stevenson, B. A.; Terranova, D.
2013-08-01
Axisymmetric equilibrium reconstruction using magnetohydrodynamic equilibrium solutions to the Grad-Shafranov equation has long been an important tool for interpreting tokamak experiments. This paper describes recent results in non-axisymmetric (three-dimensional) equilibrium reconstruction of nominally axisymmetric plasmas (tokamaks and reversed field pinches (RFPs)), and fully non-axisymmetric plasmas (stellarators). Results from applying the V3FIT code to CTH and HSX stellarator plasmas, RFX-mod RFP plasmas and the DIII-D tokamak are presented.
On the nature of incompressible magnetohydrodynamic turbulence
Gogoberidze, G.
2007-02-15
A novel model of incompressible magnetohydrodynamic turbulence in the presence of a strong external magnetic field is proposed for the explanation of recent numerical results. According to the proposed model, in the presence of the strong external magnetic field, incompressible magnetohydrodynamic turbulence becomes nonlocal in the sense that low-frequency modes cause decorrelation of interacting high-frequency modes from the inertial interval. It is shown that the obtained nonlocal spectrum of the inertial range of incompressible magnetohydrodynamic turbulence represents an anisotropic analogue of Kraichnan's nonlocal spectrum of hydrodynamic turbulence. Based on the analysis performed in the framework of the weak-coupling approximation, which represents one of the equivalent formulations of the direct interaction approximation, it is shown that incompressible magnetohydrodynamic turbulence could be both local and nonlocal, and therefore anisotropic analogues of both the Kolmogorov and Kraichnan spectra are realizable in incompressible magnetohydrodynamic turbulence.
NASA Astrophysics Data System (ADS)
Charbonneau, Paul
This chapter steps finally away from the sun and towards the stars, the idea being to apply the physical insight gained so far to see how much of stellar magnetism can be understood in terms of dynamo action. Dynamo action in the convective core of massive main-sequence stars is first considered and shown viable. For intermediate-mass main-sequence stars the fossil field hypothesis will carry the day, although possible dynamo alternatives are also briefly discussed. The extension of the solar dynamo models investigated in Chap. 3 (10.1007/978-3-642-32093-4_3) to other solar-type stars will first take us through an important detour in first having to understand rotational evolution in response to angular momentum loss in a magnetized wind. Dynamo action in fully convective stars comes next, and the chapter closes with an overview of the situation for pre- and post-main-sequence stars and compact objects, leading finally to the magnetic fields of galaxies and beyond.
Variational integrators for reduced magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Kraus, Michael; Tassi, Emanuele; Grasso, Daniela
2016-09-01
Reduced magnetohydrodynamics is a simplified set of magnetohydrodynamics equations with applications to both fusion and astrophysical plasmas, possessing a noncanonical Hamiltonian structure and consequently a number of conserved functionals. We propose a new discretisation strategy for these equations based on a discrete variational principle applied to a formal Lagrangian. The resulting integrator preserves important quantities like the total energy, magnetic helicity and cross helicity exactly (up to machine precision). As the integrator is free of numerical resistivity, spurious reconnection along current sheets is absent in the ideal case. If effects of electron inertia are added, reconnection of magnetic field lines is allowed, although the resulting model still possesses a noncanonical Hamiltonian structure. After reviewing the conservation laws of the model equations, the adopted variational principle with the related conservation laws is described both at the continuous and discrete level. We verify the favourable properties of the variational integrator in particular with respect to the preservation of the invariants of the models under consideration and compare with results from the literature and those of a pseudo-spectral code.
Representation of ideal magnetohydrodynamic modes
White, R. B.
2013-02-15
One of the most fundamental properties of ideal magnetohydrodynamics is the condition that plasma motion cannot change magnetic topology. The conventional representation of ideal magnetohydrodynamic modes by perturbing a toroidal equilibrium field through {delta}B(vector sign)={nabla} Multiplication-Sign ({xi}(vector sign) Multiplication-Sign B(vector sign)) ensures that {delta}B(vector sign){center_dot}{nabla}{psi}=0 at a resonance, with {psi} labelling an equilibrium flux surface. Also useful for the analysis of guiding center orbits in a perturbed field is the representation {delta}B(vector sign)={nabla} Multiplication-Sign {alpha}B(vector sign). These two representations are equivalent, but the vanishing of {delta}B(vector sign){center_dot}{nabla}{psi} at a resonance is necessary but not sufficient for the preservation of field line topology, and a indiscriminate use of either perturbation in fact destroys the original equilibrium flux topology. It is necessary to find the perturbed field to all orders in {xi}(vector sign) to conserve the original topology. The effect of using linearized perturbations on stability and growth rate calculations is discussed.
NASA Technical Reports Server (NTRS)
2002-01-01
[TOP LEFT AND RIGHT] The Hubble Space Telescope's Wide Field and Planetary Camera 2 has captured images of the birth of two planetary nebulae as they emerge from wrappings of gas and dust, like butterflies breaking out of their cocoons. These images highlight a fleeting phase in the stellar burnout process, occurring just before dying stars are transformed into planetary nebulae. The left-hand image is the Cotton Candy nebula, IRAS 17150-3224; the right-hand image, the Silkworm nebula, IRAS 17441-2411. Called proto-planetary nebulae, these dying stars have been caught in a transition phase between a red giant and a planetary nebula. This phase is only about 1,000 years long, very short in comparison to the 1 billion-year lifetime of a star. These images provide the earliest snapshots of the transition process. Studying images of proto-planetary nebulae is important to understanding the process of star death. A star begins to die when it has exhausted its thermonuclear fuel - hydrogen and helium. The star then becomes bright and cool (red giant phase) and swells to several tens of times its normal size. It begins puffing thin shells of gas off into space. These shells become the star's cocoon. In the Hubble images, the shells are the concentric rings seen around each nebula. But the images also reveal the nebulae breaking out from those shells. The butterfly-like wings of gas and dust are a common shape of planetary nebulae. Such butterfly shapes are created by the 'interacting winds' process, in which a more recent 'fast wind' - material propelled by radiation from the hot central star - punches a hole in the cocoon, allowing the nebula to emerge. (This 'interacting wind' theory was first proposed by Dr. Sun Kwok to explain the origin of planetary nebulae, and has been subsequently proven successful in explaining their shapes.) The nebulae are being illuminated by light from the invisible central star, which is then reflected toward us. We are viewing the nebulae
Magnetohydrodynamic turbulence: Observation and experiment
Brown, M. R.; Schaffner, D. A.; Weck, P. J.
2015-05-15
We provide a tutorial on the paradigms and tools of magnetohydrodynamic (MHD) turbulence. The principal paradigm is that of a turbulent cascade from large scales to small, resulting in power law behavior for the frequency power spectrum for magnetic fluctuations E{sub B}(f). We will describe five useful statistical tools for MHD turbulence in the time domain: the temporal autocorrelation function, the frequency power spectrum, the probability distribution function of temporal increments, the temporal structure function, and the permutation entropy. Each of these tools will be illustrated with an example taken from MHD fluctuations in the solar wind. A single dataset from the Wind satellite will be used to illustrate all five temporal statistical tools.
Method for manufacturing magnetohydrodynamic electrodes
Killpatrick, Don H.; Thresh, Henry R.
1982-01-01
A method of manufacturing electrodes for use in a magnetohydrodynamic (MHD) generator comprising the steps of preparing a billet having a core 10 of a first metal, a tubular sleeve 12 of a second metal, and an outer sheath 14, 16, 18 of an extrusile metal; evacuating the space between the parts of the assembled billet; extruding the billet; and removing the outer jacket 14. The extruded bar may be made into electrodes by cutting and bending to the shape required for an MDH channel frame. The method forms a bond between the first metal of the core 10 and the second metal of the sleeve 12 strong enough to withstand a hot and corrosive environment.
Magnetohydrodynamic Turbulence and the Geodynamo
NASA Technical Reports Server (NTRS)
Shebalin, John V.
2014-01-01
The ARES Directorate at JSC has researched the physical processes that create planetary magnetic fields through dynamo action since 2007. The "dynamo problem" has existed since 1600, when William Gilbert, physician to Queen Elizabeth I, recognized that the Earth was a giant magnet. In 1919, Joseph Larmor proposed that solar (and by implication, planetary) magnetism was due to magnetohydrodynamics (MHD), but full acceptance did not occur until Glatzmaier and Roberts solved the MHD equations numerically and simulated a geomagnetic reversal in 1995. JSC research produced a unique theoretical model in 2012 that provided a novel explanation of these physical observations and computational results as an essential manifestation of broken ergodicity in MHD turbulence. Research is ongoing, and future work is aimed at understanding quantitative details of magnetic dipole alignment in the Earth as well as in Mercury, Jupiter and its moon Ganymede, Saturn, Uranus, Neptune, and the Sun and other stars.
Electron magnetohydrodynamic turbulence: universal features
NASA Astrophysics Data System (ADS)
Shivamoggi, Bhimsen K.
2015-02-01
The energy cascade of electron magnetohydrodynamic (EMHD) turbulence is considered. Fractal and multi-fractal models for the energy dissipation field are used to determine the spatial intermittency corrections to the scaling behavior in the high-wavenumber (electron hydrodynamic limit) and low-wavenumber (magnetization limit) asymptotic regimes of the inertial range. Extrapolation of the multi-fractal scaling down to the dissipative microscales confirms in these asymptotic regimes a dissipative anomaly previously indicated by the numerical simulations of EMHD turbulence. Several basic features of the EMHD turbulent system are found to be universal which seem to transcend the existence of the characteristic length scale d e (which is the electron skin depth) in the EMHD problem: equipartition spectrum; Reynolds-number scaling of the dissipative microscales; scaling of the probability distribution function (PDF) of the electron-flow velocity (or magnetic field) gradient (even with intermittency corrections); dissipative anomaly; and critical exponent scaling.
Magnetohydrodynamics of chiral relativistic fluids
NASA Astrophysics Data System (ADS)
Boyarsky, Alexey; Fröhlich, Jürg; Ruchayskiy, Oleg
2015-08-01
We study the dynamics of a plasma of charged relativistic fermions at very high temperature T ≫m , where m is the fermion mass, coupled to the electromagnetic field. In particular, we derive a magnetohydrodynamical description of the evolution of such a plasma. We show that, compared to conventional magnetohydronamics (MHD) for a plasma of nonrelativistic particles, the hydrodynamical description of the relativistic plasma involves new degrees of freedom described by a pseudoscalar field originating in a local asymmetry in the densities of left-handed and right-handed fermions. This field can be interpreted as an effective axion field. Taking into account the chiral anomaly we present dynamical equations for the evolution of this field, as well as of other fields appearing in the MHD description of the plasma. Due to its nonlinear coupling to helical magnetic fields, the axion field significantly affects the dynamics of a magnetized plasma and can give rise to a novel type of inverse cascade.
Scaling laws in magnetohydrodynamic turbulence
Campanelli, Leonardo
2004-10-15
We analyze the decay laws of the kinetic and magnetic energies and the evolution of correlation lengths in freely decaying incompressible magnetohydrodynamic (MHD) turbulence. Scale invariance of MHD equations assures that, in the case of constant dissipation parameters (i.e., kinematic viscosity and resistivity) and null magnetic helicity, the kinetic and magnetic energies decay in time as E{approx}t{sup -1}, and the correlation lengths evolve as {xi}{approx}t{sup 1/2}. In the helical case, assuming that the magnetic field evolves towards a force-free state, we show that (in the limit of large magnetic Reynolds number) the magnetic helicity remains constant, and the kinetic and magnetic energies decay as E{sub v}{approx}t{sup -1} and E{sub B}{approx}t{sup -1/2} respectively, while both the kinetic and magnetic correlation lengths grow as {xi}{approx}t{sup 1/2}.
Anisotropic scaling of magnetohydrodynamic turbulence.
Horbury, Timothy S; Forman, Miriam; Oughton, Sean
2008-10-24
We present a quantitative estimate of the anisotropic power and scaling of magnetic field fluctuations in inertial range magnetohydrodynamic turbulence, using a novel wavelet technique applied to spacecraft measurements in the solar wind. We show for the first time that, when the local magnetic field direction is parallel to the flow, the spacecraft-frame spectrum has a spectral index near 2. This can be interpreted as the signature of a population of fluctuations in field-parallel wave numbers with a k(-2)_(||) spectrum but is also consistent with the presence of a "critical balance" style turbulent cascade. We also find, in common with previous studies, that most of the power is contained in wave vectors at large angles to the local magnetic field and that this component of the turbulence has a spectral index of 5/3. PMID:18999759
Micromachined magnetohydrodynamic actuators and sensors
Lee, Abraham P.; Lemoff, Asuncion V.
2000-01-01
A magnetohydrodynamic (MHD) micropump and microsensor which utilizes micromachining to integrate the electrodes with microchannels and includes a magnet for producing magnetic fields perpendicular to both the electrical current direction and the fluid flow direction. The magnet can also be micromachined and integrated with the micropump using existing technology. The MHD micropump, for example, can generate continuous, reversible flow, with readily controllable flow rates. The flow can be reversed by either reversing the electrical current flow or reversing the magnetic field. By mismatching the electrodes, a swirling vortex flow can be generated for potential mixing applications. No moving parts are necessary and the dead volume is minimal. The micropumps can be placed at any position in a fluidic circuit and a combination of micropumps can generate fluidic plugs and valves.
ANISOTROPIC INTERMITTENCY OF MAGNETOHYDRODYNAMIC TURBULENCE
Osman, K. T.; Kiyani, K. H.; Chapman, S. C.; Hnat, B.
2014-03-10
A higher-order multiscale analysis of spatial anisotropy in inertial range magnetohydrodynamic turbulence is presented using measurements from the STEREO spacecraft in fast ambient solar wind. We show for the first time that, when measuring parallel to the local magnetic field direction, the full statistical signature of the magnetic and Elsässer field fluctuations is that of a non-Gaussian globally scale-invariant process. This is distinct from the classic multiexponent statistics observed when the local magnetic field is perpendicular to the flow direction. These observations are interpreted as evidence for the weakness, or absence, of a parallel magnetofluid turbulence energy cascade. As such, these results present strong observational constraints on the statistical nature of intermittency in turbulent plasmas.
'Reduced' magnetohydrodynamics and minimum dissipation rates
NASA Technical Reports Server (NTRS)
Montgomery, David
1992-01-01
It is demonstrated that all solutions of the equations of 'reduced' magnetohydrodynamics approach a uniform-current, zero-flow state for long times, given a constant wall electric field, uniform scalar viscosity and resistivity, and uniform mass density. This state is the state of minimum energy dissipation rate for these boundary conditions. No steady-state turbulence is possible. The result contrasts sharply with results for full three-dimensional magnetohydrodynamics before the reduction occurs.
History of Stellar Interferometry
NASA Technical Reports Server (NTRS)
Lawson, Peter R.
2004-01-01
This viewgraph presentation reviews the history of stellar interferometry from the suggestion of Fizeau that stellar interferometry was possible,to the use of the Mark I, II and III for astrometry. Photographs, and parts of original articles are presented.
Computational Methods for Ideal Magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Kercher, Andrew D.
Numerical schemes for the ideal magnetohydrodynamics (MHD) are widely used for modeling space weather and astrophysical flows. They are designed to resolve the different waves that propagate through a magnetohydro fluid, namely, the fast, Alfven, slow, and entropy waves. Numerical schemes for ideal magnetohydrodynamics that are based on the standard finite volume (FV) discretization exhibit pseudo-convergence in which non-regular waves no longer exist only after heavy grid refinement. A method is described for obtaining solutions for coplanar and near coplanar cases that consist of only regular waves, independent of grid refinement. The method, referred to as Compound Wave Modification (CWM), involves removing the flux associated with non-regular structures and can be used for simulations in two- and three-dimensions because it does not require explicitly tracking an Alfven wave. For a near coplanar case, and for grids with 213 points or less, we find root-mean-square-errors (RMSEs) that are as much as 6 times smaller. For the coplanar case, in which non-regular structures will exist at all levels of grid refinement for standard FV schemes, the RMSE is as much as 25 times smaller. A multidimensional ideal MHD code has been implemented for simulations on graphics processing units (GPUs). Performance measurements were conducted for both the NVIDIA GeForce GTX Titan and Intel Xeon E5645 processor. The GPU is shown to perform one to two orders of magnitude greater than the CPU when using a single core, and two to three times greater than when run in parallel with OpenMP. Performance comparisons are made for two methods of storing data on the GPU. The first approach stores data as an Array of Structures (AoS), e.g., a point coordinate array of size 3 x n is iterated over. The second approach stores data as a Structure of Arrays (SoA), e.g. three separate arrays of size n are iterated over simultaneously. For an AoS, coalescing does not occur, reducing memory efficiency
The Stellar Imager (SI) - A Mission to Resolve Stellar Surfaces, Interiors, and Magnetic Activity
NASA Technical Reports Server (NTRS)
Carpenter, Kenneth; Schrijver, Carolus J.; Karovska, Margarita
2007-01-01
The Stellar Imager (SI) is a UV/Optical, Space-Based Interferometer designed to enable 0.1 milli-arcsecond (mas) spectral imaging of stellar surfaces and, via asteroseismology, stellar interiors and of the Universe in general. The ultra-sharp images of the Stellar Imager will revolutionize our view of many dynamic astrophysical processes by transforming point sources into extended sources, and snapshots into evolving views. SI's science focuses on the role of magnetism in the Universe, particularly on magnetic activity on the surfaces of stars like the Sun. SI's prime goal is to enable long-term forecasting of solar activity and the space weather that it drives. SI will also revolutionize our understanding of the formation of planetary systems, of the habitability and climatology of distant planets, and of many magneto-hydrodynamically controlled processes in the Universe. SI is included as a 'Flagship and Landmark Discovery Mission' in the 2005 NASA Sun Solar System Connection (SSSC) Roadmap and as a candidate for a 'Pathways to Life Observatory' in the NASA Exploration of the Universe Division (EUD) Roadmap (May, 2005). In this paper we discuss the science goals and technology needs of, and the baseline design for, the SI Mission (http://hires.gsfc.nasa.gov/si/) its ability to image the 'Biggest, Baddest, Coolest Stars'.
The Stellar Imager (SI) - A Mission to Resolve Stellar Surfaces, Interiors, and Magnetic Activity
NASA Astrophysics Data System (ADS)
Christensen-Dalsgaard, Jørgen; Carpenter, Kenneth G.; Schrijver, Carolus J.; Karovska, Margarita; Si Team
2011-01-01
The Stellar Imager (SI) is a space-based, UV/Optical Interferometer (UVOI) designed to enable 0.1 milli-arcsecond (mas) spectral imaging of stellar surfaces and of the Universe in general. It will also probe via asteroseismology flows and structures in stellar interiors. SI will enable the development and testing of a predictive dynamo model for the Sun, by observing patterns of surface activity and imaging of the structure and differential rotation of stellar interiors in a population study of Sun-like stars to determine the dependence of dynamo action on mass, internal structure and flows, and time. SI's science focuses on the role of magnetism in the Universe and will revolutionize our understanding of the formation of planetary systems, of the habitability and climatology of distant planets, and of many magneto-hydrodynamically controlled processes in the Universe. SI is a "Landmark/Discovery Mission" in the 2005 Heliophysics Roadmap, an implementation of the UVOI in the 2006 Astrophysics Strategic Plan, and a NASA Vision Mission ("NASA Space Science Vision Missions" (2008), ed. M. Allen). We present here the science goals of the SI Mission, a mission architecture that could meet those goals, and the technology development needed to enable this mission. Additional information on SI can be found at: http://hires.gsfc.nasa.gov/si/.
Buoyancy-driven Magnetohydrodynamic Waves
NASA Astrophysics Data System (ADS)
Hague, A.; Erdélyi, R.
2016-09-01
Turbulent motions close to the visible solar surface may generate low-frequency internal gravity waves (IGWs) that propagate through the lower solar atmosphere. Magnetic activity is ubiquitous throughout the solar atmosphere, so it is expected that the behavior of IGWs is to be affected. In this article we investigate the role of an equilibrium magnetic field on propagating and standing buoyancy oscillations in a gravitationally stratified medium. We assume that this background magnetic field is parallel to the direction of gravitational stratification. It is known that when the equilibrium magnetic field is weak and the background is isothermal, the frequencies of standing IGWs are sensitive to the presence of magnetism. Here, we generalize this result to the case of a slowly varying temperature. To do this, we make use of the Boussinesq approximation. A comparison between the hydrodynamic and magnetohydrodynamic cases allows us to deduce the effects due to a magnetic field. It is shown that the frequency of IGWs may depart significantly from the Brunt–Väisälä frequency, even for a weak magnetic field. The mathematical techniques applied here give a clearer picture of the wave mode identification, which has previously been misinterpreted. An observational test is urged to validate the theoretical findings.
Magnetohydrodynamic Simulations of Barred Galaxies
NASA Astrophysics Data System (ADS)
Kim, W.-T.
2013-04-01
Magnetic fields are pervasive in barred galaxies, especially in gaseous substructures such as dust lanes and nuclear rings. To explore the effects of magnetic fields on the formation of the substructures as well as on the mass inflow rates to the galaxy center, we run two-dimensional, ideal magnetohydrodynamic simulations. We use a modified version of the Athena code whose numerical magnetic diffusivity is shown to be of third order in space. In the bar regions, magnetic fields are compressed and abruptly bent around the dust-lane shocks. The associated magnetic stress not only reduces the peak density of the dust-lane shocks but also removes angular momentum further from the gas that is moving radially in. Nuclear rings that form at the location of centrifugal barrier rather than resonance with the bar are smaller and more radially distributed, and the mass flow rate to the galaxy center is correspondingly larger in models with stronger magnetic fields. Outside the bar regions, the bar potential and strong shear conspire to amplify the field strength near the corotation resonance. The amplified fields transport angular momentum outward, producing trailing magnetic arms with strong fields and low density. The base of the magnetic arms are found to be unstable to a tearing-mode instability of magnetic reconnection. This produces numerous magnetic islands that eventually make the outer regions highly chaotic.
Magnetohydrodynamic (MHD) driven droplet mixer
Lee, Abraham P.; Lemoff, Asuncion V.; Miles, Robin R.
2004-05-11
A magnetohydrodynamic fluidic system mixes a first substance and a second substance. A first substrate section includes a first flow channel and a first plurality of pairs of spaced electrodes operatively connected to the first flow channel. A second substrate section includes a second flow channel and a second plurality of pairs of spaced electrodes operatively connected to the second flow channel. A third substrate section includes a third flow channel and a third plurality of pairs of spaced electrodes operatively connected to the third flow channel. A magnetic section and a control section are operatively connected to the spaced electrodes. The first substrate section, the second substrate section, the third substrate section, the first plurality of pairs of spaced electrodes, the second plurality of pairs of spaced electrodes, the third plurality of pairs of spaced electrodes, the magnetic section, and the control section are operated to move the first substance through the first flow channel, the second substance through the second flow channel, and both the first substance and the second substance into the third flow channel where they are mixed.
Relativistic magnetohydrodynamics in one dimension.
Lyutikov, Maxim; Hadden, Samuel
2012-02-01
We derive a number of solutions for one-dimensional dynamics of relativistic magnetized plasma that can be used as benchmark estimates in relativistic hydrodynamic and magnetohydrodynamic numerical codes. First, we analyze the properties of simple waves of fast modes propagating orthogonally to the magnetic field in relativistically hot plasma. The magnetic and kinetic pressures obey different equations of state, so that the system behaves as a mixture of gases with different polytropic indices. We find the self-similar solutions for the expansion of hot strongly magnetized plasma into vacuum. Second, we derive linear hodograph and Darboux equations for the relativistic Khalatnikov potential, which describe arbitrary one-dimensional isentropic relativistic motion of cold magnetized plasma and find their general and particular solutions. The obtained hodograph and Darboux equations are very powerful: A system of highly nonlinear, relativistic, time-dependent equations describing arbitrary (not necessarily self-similar) dynamics of highly magnetized plasma reduces to a single linear differential equation. PMID:22463331
Magnetohydrodynamic Propulsion for the Classroom
NASA Astrophysics Data System (ADS)
Font, Gabriel I.; Dudley, Scott C.
2004-10-01
The cinema industry can sometimes prove to be an ally when searching for material with which to motivate students to learn physics. Consider, for example, the electromagnetic force on a current in the presence of a magnetic field. This phenomenon is at the heart of magnetohydrodynamic (MHD) propulsion systems. A submarine employing this type of propulsion was immortalized in the movie Hunt for Red October. While mentioning this to students certainly gets their attention, it often elicits comments that it is only fiction and not physically possible. Imagine their surprise when a working system is demonstrated! It is neither difficult nor expensive to construct a working system that can be demonstrated in the front of a classroom.2 In addition, all aspects of the engineering hurdles that must be surmounted and myths concerning this "silent propulsion" system are borne out in a simple apparatus. This paper details how to construct an inexpensive MHD propulsion boat that can be demonstrated for students in the classroom.
Two-dimensional magnetohydrodynamic equilibria with flow and studies of equilibrium fluctuations
Agim, Y.Z.
1989-01-01
A set of reduced ideal MHD (magnetohydrodynamic) equations is derived to investigate equilibria of plasmas with mass flow in general two-dimensional geometry. These equations provide a means of investigating the effects of flow on self-consistent equilibria in a number of new two-dimensional configurations such as helically symmetric configurations with helical axis, which are relevant to stellarators, as well as axisymmetric configurations. In the second part, magnetic fluctuations due to the thermally excited MHD waves are investigated using fluid and kinetic models to describe a stable, uniform, compressible plasma in the range above the drift wave frequency and below the ion cyclotron frequency.
NASA Astrophysics Data System (ADS)
Chabi Orou, Jean Bio; Pomalégni, Gisèle
2015-11-01
We investigate the combined effects of rotation , magnetic field and helical force on the onset of stationary and oscillatory convection in a horizontal electrically conducting fluid layer heated from below with free-free boundary conditions. For this investigation the linear stability analysis studied by Chandrasekhar (1961) is used. We obtain the condition for the formation of a single large scale structure. In (Pomalégni et al., 2014) it was shown the existence of a critical value Scr of the intensity of the helical force for which the apparition of two cells at marginal stability for the oscillatory convection is obtained. Then, we have shown here how the increasing of parameter Ta influences this critical value of the helical force intensity.
NASA Astrophysics Data System (ADS)
Fable, E.; Angioni, C.; Casson, F. J.; Told, D.; Ivanov, A. A.; Jenko, F.; McDermott, R. M.; Medvedev, S. Yu; Pereverzev, G. V.; Ryter, F.; Treutterer, W.; Viezzer, E.; the ASDEX Upgrade Team
2013-12-01
Tokamak scenario development requires an understanding of the properties that determine the kinetic profiles in non-steady plasma phases and of the self-consistent evolution of the magnetic equilibrium. Current ramps are of particular interest since many transport-relevant parameters explore a large range of values and their impact on transport mechanisms has to be assessed. To this purpose, a novel full-discharge modelling tool has been developed, which couples the transport code ASTRA (Pereverzev et al 1991 IPP Report 5/42) and the free boundary equilibrium code SPIDER (Ivanov et al 2005 32nd EPS Conf. on Plasma Physics vol 29C (ECA) P-5.063 and http://epsppd.epfl.ch/Tarragona/pdf/P5_063.pdf), utilizing a specifically designed coupling scheme. The current ramp-up phase can be accurately and reliably simulated using this scheme, where a plasma shape, position and current controller is applied, which mimics the one of ASDEX Upgrade. Transport of energy is provided by theory-based models (e.g. TGLF (Staebler et al 2007 Phys. Plasmas 14 055909)). A recipe based on edge-relevant parameters (Scott 2000 Phys. Plasmas 7 1845) is proposed to resolve the low current phase of the current ramps, where the impact of the safety factor on micro-instabilities could make quasi-linear approaches questionable in the plasma outer region. Current ramp scenarios, selected from ASDEX Upgrade discharges, are then simulated to validate both the coupling with the free-boundary evolution and the prediction of profiles. Analysis of the underlying transport mechanisms is presented, to clarify the possible physics origin of the observed L-mode empirical energy confinement scaling. The role of toroidal micro-instabilities (ITG, TEM) and of non-linear effects is discussed.
From the Einstein-Szilard Patent to Modern Magnetohydrodynamics.
ERIC Educational Resources Information Center
Povh, I. L.; Barinberg, A. D.
1979-01-01
Examines present-day and future prospects of the applications of modern magnetohydrodynamics in a number of countries. Explains how the electromagnetic pump, which was invented by Einstein and Leo Szilard, led to the development of applied magnetohydrodynamics. (HM)
Electron magnetohydrodynamics: dynamics and turbulence.
Lyutikov, Maxim
2013-11-01
We consider dynamics and turbulent interaction of whistler modes within the framework of inertialess electron magnetohydrodynamics (EMHD). We argue that there is no energy principle in EMHD: any stationary closed configuration is neutrally stable. On the other hand, the relaxation principle, the long term evolution of a weakly dissipative system towards Taylor-Beltrami state, remains valid in EMHD. We consider the turbulent cascade of whistler modes. We show that (i) harmonic whistlers are exact nonlinear solutions; (ii) collinear whistlers do not interact (including counterpropagating); (iii) waves with the same value of the wave vector k(1)=k(2) do not interact; (iv) whistler modes have a dispersion that allows a three-wave decay, including into a zero frequency mode; (v) the three-wave interaction effectively couples modes with highly different wave numbers and propagation angles. In addition, linear interaction of a whistler with a single zero mode can lead to spatially divergent structures via parametric instability. All these properties are drastically different from MHD, so that the qualitative properties of the Alfvén turbulence can not be transferred to the EMHD turbulence. We derive the Hamiltonian formulation of EMHD, and using Bogoliubov transformation reduce it to the canonical form; we calculate the matrix elements for the three-wave interaction of whistlers. We solve numerically the kinetic equation and show that, generally, the EMHD cascade develops within a broad range of angles, while transiently it may show anisotropic, nearly two-dimensional structures. Development of a cascade depends on the forcing (nonuniversal) and often fails to reach a steady state. Analytical estimates predict the spectrum of magnetic fluctuations for the quasi-isotropic cascade [proportionality]k(-2). The cascade remains weak (not critically balanced). The cascade is UV local, while the infrared locality is weakly (logarithmically) violated. PMID:24329368
Accurate, meshless methods for magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Hopkins, Philip F.; Raives, Matthias J.
2016-01-01
Recently, we explored new meshless finite-volume Lagrangian methods for hydrodynamics: the `meshless finite mass' (MFM) and `meshless finite volume' (MFV) methods; these capture advantages of both smoothed particle hydrodynamics (SPH) and adaptive mesh refinement (AMR) schemes. We extend these to include ideal magnetohydrodynamics (MHD). The MHD equations are second-order consistent and conservative. We augment these with a divergence-cleaning scheme, which maintains nabla \\cdot B≈ 0. We implement these in the code GIZMO, together with state-of-the-art SPH MHD. We consider a large test suite, and show that on all problems the new methods are competitive with AMR using constrained transport (CT) to ensure nabla \\cdot B=0. They correctly capture the growth/structure of the magnetorotational instability, MHD turbulence, and launching of magnetic jets, in some cases converging more rapidly than state-of-the-art AMR. Compared to SPH, the MFM/MFV methods exhibit convergence at fixed neighbour number, sharp shock-capturing, and dramatically reduced noise, divergence errors, and diffusion. Still, `modern' SPH can handle most test problems, at the cost of larger kernels and `by hand' adjustment of artificial diffusion. Compared to non-moving meshes, the new methods exhibit enhanced `grid noise' but reduced advection errors and diffusion, easily include self-gravity, and feature velocity-independent errors and superior angular momentum conservation. They converge more slowly on some problems (smooth, slow-moving flows), but more rapidly on others (involving advection/rotation). In all cases, we show divergence control beyond the Powell 8-wave approach is necessary, or all methods can converge to unphysical answers even at high resolution.
Electron magnetohydrodynamics: Dynamics and turbulence
NASA Astrophysics Data System (ADS)
Lyutikov, Maxim
2013-11-01
We consider dynamics and turbulent interaction of whistler modes within the framework of inertialess electron magnetohydrodynamics (EMHD). We argue that there is no energy principle in EMHD: any stationary closed configuration is neutrally stable. On the other hand, the relaxation principle, the long term evolution of a weakly dissipative system towards Taylor-Beltrami state, remains valid in EMHD. We consider the turbulent cascade of whistler modes. We show that (i) harmonic whistlers are exact nonlinear solutions; (ii) collinear whistlers do not interact (including counterpropagating); (iii) waves with the same value of the wave vector k1=k2 do not interact; (iv) whistler modes have a dispersion that allows a three-wave decay, including into a zero frequency mode; (v) the three-wave interaction effectively couples modes with highly different wave numbers and propagation angles. In addition, linear interaction of a whistler with a single zero mode can lead to spatially divergent structures via parametric instability. All these properties are drastically different from MHD, so that the qualitative properties of the Alfvén turbulence can not be transferred to the EMHD turbulence. We derive the Hamiltonian formulation of EMHD, and using Bogoliubov transformation reduce it to the canonical form; we calculate the matrix elements for the three-wave interaction of whistlers. We solve numerically the kinetic equation and show that, generally, the EMHD cascade develops within a broad range of angles, while transiently it may show anisotropic, nearly two-dimensional structures. Development of a cascade depends on the forcing (nonuniversal) and often fails to reach a steady state. Analytical estimates predict the spectrum of magnetic fluctuations for the quasi-isotropic cascade ∝k-2. The cascade remains weak (not critically balanced). The cascade is UV local, while the infrared locality is weakly (logarithmically) violated.
Saturation of Stellar Winds from Young Suns
NASA Astrophysics Data System (ADS)
Suzuki, Takeru K.; Imada, Shinsuke; Kataoka, Ryuho; Kato, Yoshiaki; Matsumoto, Takuma; Miyahara, Hiroko; Tsuneta, Saku
2013-10-01
We investigated mass losses via stellar winds from Sun-like main-sequence stars with a wide range of activity levels. We performed forward-type magnetohydrodynamical numerical experiments for Alfvén wave-driven stellar winds with a wide range of input Poynting flux from the photosphere. Increasing the magnetic field strength and the turbulent velocity at the stellar photosphere from the current solar level, the mass-loss rate rapidly at first increases, owing to suppression of the reflection of the Alfvén waves. The surface materials are lifted up by the magnetic pressure associated with the Alfvén waves, and the cool dense chromosphere is intermittently extended to 10%#8211;20% of the stellar radius. The dense atmospheres enhance the radiative losses, and eventually most of the input Poynting energy from the stellar surface escapes by radiation. As a result, there is no more sufficient energy remaining for the kinetic energy of the wind; the stellar wind saturates in very active stars, as observed in Wood et al. (2002, ApJ, 574, 412; 2005, ApJ, 628, L143). The saturation level is positively correlated with Br,0 f0, where Br,0 and f0 are the magnetic field strength and the filling factor of open flux tubes at the photosphere. If Br,0 f0 is relatively large gtrsim 5 G, the mass-loss rate could be as high as 1000 times. If such a strong mass loss lasts for ˜ 1 billion years, the stellar mass itself would be affected, which could be a solution to the faint young Sun paradox. We derived a Reimers-type scaling relation that estimates the mass-loss rate from an energetics consideration of our simulations. Finally, we derived the evolution of the mass-loss rates, dot;{M} ∝ t-1.23, of our simulations, combining with an observed time evolution of X-ray flux from Sun-like stars, which are shallower than dot;{M} ∝ t-2.33±0.55 in Wood et al. (2005).
Magnetohydrodynamically generated velocities in confined plasma
Morales, Jorge A. Bos, Wouter J. T.; Schneider, Kai; Montgomery, David C.
2015-04-15
We investigate by numerical simulation the rotational flows in a toroid confining a conducting magnetofluid in which a current is driven by the application of externally supported electric and magnetic fields. The computation involves no microscopic instabilities and is purely magnetohydrodynamic (MHD). We show how the properties and intensity of the rotations are regulated by dimensionless numbers (Lundquist and viscous Lundquist) that contain the resistivity and viscosity of the magnetofluid. At the magnetohydrodynamic level (uniform mass density and incompressible magnetofluids), rotational flows appear in toroidal, driven MHD. The evolution of these flows with the transport coefficients, geometry, and safety factor are described.
Generalized magnetofluid connections in relativistic magnetohydrodynamics.
Asenjo, Felipe A; Comisso, Luca
2015-03-20
The concept of magnetic connections is extended to nonideal relativistic magnetohydrodynamical plasmas. Adopting a general set of equations for relativistic magnetohydrodynamics including thermal-inertial, thermal electromotive, Hall, and current-inertia effects, we derive a new covariant connection equation showing the existence of generalized magnetofluid connections that are preserved during the dissipationless plasma dynamics. These connections are intimately linked to a general antisymmetric tensor that unifies the electromagnetic and fluid fields, allowing the extension of the magnetic connection notion to a much broader concept. PMID:25839284
New approach to nonrelativistic ideal magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Banerjee, Rabin; Kumar, Kuldeep
2016-07-01
We provide a novel action principle for nonrelativistic ideal magnetohydrodynamics in the Eulerian scheme exploiting a Clebsch-type parametrisation. Both Lagrangian and Hamiltonian formulations have been considered. Within the Hamiltonian framework, two complementary approaches have been discussed using Dirac's constraint analysis. In one case the Hamiltonian is canonical involving only physical variables but the brackets have a noncanonical structure, while the other retains the canonical structure of brackets by enlarging the phase space. The special case of incompressible magnetohydrodynamics is also considered where, again, both the approaches are discussed in the Hamiltonian framework. The conservation of the stress tensor reveals interesting aspects of the theory.
MHD computations for stellarators
Johnson, J.L.
1985-12-01
Considerable progress has been made in the development of computational techniques for studying the magnetohydrodynamic equilibrium and stability properties of three-dimensional configurations. Several different approaches have evolved to the point where comparison of results determined with different techniques shows good agreement. 55 refs., 7 figs.
Direct UV/Optical Imaging of Stellar Surfaces: The Stellar Imager (SI) Vision Mission
NASA Technical Reports Server (NTRS)
Carpenter, Kenneth G.; Lyon, Richard G.; Schrijver, Carolus; Karovska, Margarita; Mozurkewich, David
2007-01-01
The Stellar Imager (SI) is a UV/optical, space-based interferometer designed to enable 0.1 milli-arcsecond (mas) spectral imaging of stellar surfaces and, via asteroseismology, stellar interiors and of the Universe in general. SI's science focuses on the role of magnetism in the Universe, particularly on magnetic activity on the surfaces of stars like the Sun. SI's prime goal is to enable long-term forecasting of solar activity and the space weather that it drives, in support of the Living with a Star program in the Exploration Era. SI will also revolutionize our understanding of the formation of planetary systems, of the habitability and climatology of distant planets, and of many magneto-hydrodynamically controlled processes in thc Universe. SI is a "Flagship and Landmark Discovery Mission" in the 2005 Sun Solar System Connection (SSSC) Roadmap and a candidate for a "Pathways to Life Observatory" in the Exploration of the Universe Division (EUD) Roadmap. We discuss herein the science goals of the SI Mission, a mission architecture that could meet those goals, and the technologies needed to enable this mission. Additional information on SI can be found at: http://hires.gsfc.nasa.gov/si/.
Global Magnetohydrodynamic Modeling of the Solar Corona
NASA Technical Reports Server (NTRS)
Linker, Jon A.
2001-01-01
This report describes the progress made in the investigation of the solar corona using magnetohydrodynamic (MHD) simulations. Coronal mass ejections (CME) are believed to be the primary cause of nonrecurrent geomagnetic storms and these have been investigated through the use of three-dimensional computer simulation.
On energy conservation in extended magnetohydrodynamics
Kimura, Keiji; Morrison, P. J.
2014-08-15
A systematic study of energy conservation for extended magnetohydrodynamic models that include Hall terms and electron inertia is performed. It is observed that commonly used models do not conserve energy in the ideal limit, i.e., when viscosity and resistivity are neglected. In particular, a term in the momentum equation that is often neglected is seen to be needed for conservation of energy.
Solar-driven liquid metal magnetohydrodynamic generator
NASA Technical Reports Server (NTRS)
Lee, J. H.; Hohl, F.
1981-01-01
A solar oven heated by concentrated solar radiation as the heat source of a liquid metal magnetohydrodynamic (LMMHD) power generation system is proposed. The design allows the production of electric power in space, as well as on Earth, at high rates of efficiency. Two types of the solar oven suitable for the system are discussed.
Solar-driven liquid metal magnetohydrodynamic generator
NASA Astrophysics Data System (ADS)
Lee, J. H.; Hohl, F.
1981-05-01
A solar oven heated by concentrated solar radiation as the heat source of a liquid metal magnetohydrodynamic (LMMHD) power generation system is proposed. The design allows the production of electric power in space, as well as on Earth, at high rates of efficiency. Two types of the solar oven suitable for the system are discussed.
Validation of Magnetospheric Magnetohydrodynamic Models
NASA Astrophysics Data System (ADS)
Curtis, Brian
Magnetospheric magnetohydrodynamic (MHD) models are commonly used for both prediction and modeling of Earth's magnetosphere. To date, very little validation has been performed to determine their limits, uncertainties, and differences. In this work, we performed a comprehensive analysis using several commonly used validation techniques in the atmospheric sciences to MHD-based models of Earth's magnetosphere for the first time. The validation techniques of parameter variability/sensitivity analysis and comparison to other models were used on the OpenGGCM, BATS-R-US, and SWMF magnetospheric MHD models to answer several questions about how these models compare. The questions include: (1) the difference between the model's predictions prior to and following to a reversal of Bz in the upstream interplanetary field (IMF) from positive to negative, (2) the influence of the preconditioning duration, and (3) the differences between models under extreme solar wind conditions. A differencing visualization tool was developed and used to address these three questions. We find: (1) For a reversal in IMF Bz from positive to negative, the OpenGGCM magnetopause is closest to Earth as it has the weakest magnetic pressure near-Earth. The differences in magnetopause positions between BATS-R-US and SWMF are explained by the influence of the ring current, which is included in SWMF. Densities are highest for SWMF and lowest for OpenGGCM. The OpenGGCM tail currents differ significantly from BATS-R-US and SWMF; (2) A longer preconditioning time allowed the magnetosphere to relax more, giving different positions for the magnetopause with all three models before the IMF Bz reversal. There were differences greater than 100% for all three models before the IMF Bz reversal. The differences in the current sheet region for the OpenGGCM were small after the IMF Bz reversal. The BATS-R-US and SWMF differences decreased after the IMF Bz reversal to near zero; (3) For extreme conditions in the solar
NASA Astrophysics Data System (ADS)
Diehl, Hans Werner; Rutkevich, Sergei
The O (2) ϕ4 model on a 3D film of thickness L with free boundaries is relevant for the explanation of the thinning of wetting layers of 4He caused by critical Casimir forces near and below the λ-transition. Just as its O (n) analog, the model has long-range order below the bulk critical temperature Tc if L = ∞ , but remains disordered for all T > 0 when L < ∞ . A proper analysis of its scaling behavior near Tc is challenging: it involves a nontrivial dimensional crossover in addition to bulk, boundary, and finite-size critical behaviors. The n --> ∞ limit of the model can be solved exactly in terms of the eigenvalues and eigenfunctions of a self-consistent Schrödinger equation whose potential v (z) becomes singular at the boundary planes. Complementing recent numerically exact results, we derive various exact analytical results for series expansion coefficients of v (z) , its L = ∞ scattering data for all values m ⋛ 0 of the termperature scaling field, and the low-temperature asymptotic behavior of the residual free energy and the Casimir force using a combination of boundary-operator and short-distance expansions, proper extensions of inverse scattering theory, new trace formulae, and semi-classical expansions.
NASA Technical Reports Server (NTRS)
Musielak, Z. E.; Rosner, R.
1988-01-01
Magnetohydrodynamic (MHD) wave energy fluxes for late-type stars are calculated, using previously obtained formulae for the source functions for the generation of MHD waves in a stratified, but otherwise uniform, turbulent atmosphere; the magnetic fields in the wave generation region are assumed to be homogeneous. In contradiction to previous results, it is shown that in this uniform magnetic field case there is no significant increase in the efficiency of MHD wave generation, at least within the theory's limits of applicability. The major results are that the MHD energy fluxes calculated for late-type stars are less than those obtained for compressible modes in the magnetic field-free case, and that these MHD energy fluxes do not vary enough for a given spectral type to explain the observed range of UV and X-ray fluxes from such stars. It is therefore concluded that MHD waves in stellar atmospheres with homogeneous magnetic fields in the wave generation region cannot explain the observed stellar coronal emissions; if such MHD waves are responsible for a significant component of stellar coronal heating, then nonuniform fields within the generation region must be appealed to.
Stellar feedback efficiencies: supernovae versus stellar winds
NASA Astrophysics Data System (ADS)
Fierlinger, Katharina M.; Burkert, Andreas; Ntormousi, Evangelia; Fierlinger, Peter; Schartmann, Marc; Ballone, Alessandro; Krause, Martin G. H.; Diehl, Roland
2016-02-01
Stellar winds and supernova (SN) explosions of massive stars (`stellar feedback') create bubbles in the interstellar medium (ISM) and insert newly produced heavy elements and kinetic energy into their surroundings, possibly driving turbulence. Most of this energy is thermalized and immediately removed from the ISM by radiative cooling. The rest is available for driving ISM dynamics. In this work we estimate the amount of feedback energy retained as kinetic energy when the bubble walls have decelerated to the sound speed of the ambient medium. We show that the feedback of the most massive star outweighs the feedback from less massive stars. For a giant molecular cloud (GMC) mass of 105 M⊙ (as e.g. found in the Orion GMCs) and a star formation efficiency of 8 per cent the initial mass function predicts a most massive star of approximately 60 M⊙. For this stellar evolution model we test the dependence of the retained kinetic energy of the cold GMC gas on the inclusion of stellar winds. In our model winds insert 2.34 times the energy of an SN and create stellar wind bubbles serving as pressure reservoirs. We find that during the pressure-driven phases of the bubble evolution radiative losses peak near the contact discontinuity (CD), and thus the retained energy depends critically on the scales of the mixing processes across the CD. Taking into account the winds of massive stars increases the amount of kinetic energy deposited in the cold ISM from 0.1 per cent to a few per cent of the feedback energy.
Frontiers of stellar evolution
NASA Technical Reports Server (NTRS)
Lambert, David L. (Editor)
1991-01-01
The present conference discusses theoretical and observational views of star formation, spectroscopic constraints on the evolution of massive stars, very low mass stars and brown dwarfs, asteroseismology, globular clusters as tests of stellar evolution, observational tests of stellar evolution, and mass loss from cool evolved giant stars. Also discussed are white dwarfs and hot subdwarfs, neutron stars and black holes, supernovae from single stars, close binaries with evolved components, accretion disks in interacting binaries, supernovae in binary systems, stellar evolution and galactic chemical evolution, and interacting binaries containing compact components.
Moroz, P.E.
1997-09-01
A new stellarator configuration, the Double-Helix Stellarator (DHS), is introduced. This novel configuration features a double-helix center post as the only helical element of the stellarator coil system. The DHS configuration has many unique characteristics. One of them is the extreme low plasma aspect ratio, A {approx} 1--1.2. Other advantages include a high enclosed volume, appreciable rotational transform, and a possibility of extreme-high-{beta} MHD equilibria. Moreover, the DHS features improved transport characteristics caused by the absence of the magnetic field ripple on the outboard of the torus. Compactness, simplicity and modularity of the coil system add to the DHS advantages for fusion applications.
Moroz, P.E.
1997-03-01
A novel concept for magnetic plasma confinement, Stellarator-Spheromak (SSP), is proposed. Numerical analysis with the classical-stellarator-type outboard stellarator windings demonstrates a number of potential advantages of SSP for controlled nuclear fusion. Among the main ones are: simple and compact magnet coil configuration, absence of material structures (e.g. magnet coils or conducting walls) in the center of the torus, high rotational transform, and a possibility of MHD equilibria with very high {beta} (pressure/magnetic pressure) of the confined plasma.
Multi-region relaxed magnetohydrodynamics with flow
Dennis, G. R. Dewar, R. L.; Hole, M. J.; Hudson, S. R.
2014-04-15
We present an extension of the multi-region relaxed magnetohydrodynamics (MRxMHD) equilibrium model that includes plasma flow. This new model is a generalization of Woltjer's model of relaxed magnetohydrodynamics equilibria with flow. We prove that as the number of plasma regions becomes infinite, our extension of MRxMHD reduces to ideal MHD with flow. We also prove that some solutions to MRxMHD with flow are not time-independent in the laboratory frame, and instead have 3D structure which rotates in the toroidal direction with fixed angular velocity. This capability gives MRxMHD potential application to describing rotating 3D MHD structures such as 'snakes' and long-lived modes.
Geomagnetic main field modeling using magnetohydrodynamic constraints
NASA Technical Reports Server (NTRS)
Estes, R. H.
1985-01-01
The influence of physical constraints are investigated which may be approximately satisfied by the Earth's liquid core on models of the geomagnetic main field and its secular variation. A previous report describes the methodology used to incorporate nonlinear equations of constraint into the main field model. The application of that methodology to the GSFC 12/83 field model to test the frozen-flux hypothesis and the usefulness of incorporating magnetohydrodynamic constraints for obtaining improved geomagnetic field models is described.
Coal-burning magnetohydrodynamic power generation
Kessler, R.; Hals, F. )
1992-01-01
In this paper, coal-burning magnetohydrodynamic (MHD) electric power generation technology is described, and its economic and environmental advantages are discussed. advanced MHD/steam plants can achieve efficiencies of 55%-60% with less environmental intrusion than form conventional coal-burning steam plants. The national program for development of MHD power generation is outlined and the development status of individual components and subsystems is presented.
Lattice Boltzmann model for simulation of magnetohydrodynamics
NASA Technical Reports Server (NTRS)
Chen, Shiyi; Chen, Hudong; Martinez, Daniel; Matthaeus, William
1991-01-01
A numerical method, based on a discrete Boltzmann equation, is presented for solving the equations of magnetohydrodynamics (MHD). The algorithm provides advantages similar to the cellular automaton method in that it is local and easily adapted to parallel computing environments. Because of much lower noise levels and less stringent requirements on lattice size, the method appears to be more competitive with traditional solution methods. Examples show that the model accurately reproduces both linear and nonlinear MHD phenomena.
Potential vorticity formulation of compressible magnetohydrodynamics.
Arter, Wayne
2013-01-01
Compressible ideal magnetohydrodynamics is formulated in terms of the time evolution of potential vorticity and magnetic flux per unit mass using a compact Lie bracket notation. It is demonstrated that this simplifies analytic solution in at least one very important situation relevant to magnetic fusion experiments. Potentially important implications for analytic and numerical modelling of both laboratory and astrophysical plasmas are also discussed. PMID:23383802
Magnetohydrodynamic equilibria with incompressible flows: Symmetry approach
Cicogna, G.; Pegoraro, F.
2015-02-15
We identify and discuss a family of azimuthally symmetric, incompressible, magnetohydrodynamic plasma equilibria with poloidal and toroidal flows in terms of solutions of the Generalized Grad Shafranov (GGS) equation. These solutions are derived by exploiting the incompressibility assumption, in order to rewrite the GGS equation in terms of a different dependent variable, and the continuous Lie symmetry properties of the resulting equation and, in particular, a special type of “weak” symmetries.
Stellar atmospheric structural patterns
NASA Technical Reports Server (NTRS)
Thomas, R. N.
1983-01-01
The thermodynamics of stellar atmospheres is discussed. Particular attention is given to the relation between theoretical modeling and empirical evidence. The characteristics of distinctive atmospheric regions and their radical structures are discussed.
Guiding center equations for ideal magnetohydrodynamic modes
White, R. B.
2013-04-15
Guiding center simulations are routinely used for the discovery of mode-particle resonances in tokamaks, for both resistive and ideal instabilities and to find modifications of particle distributions caused by a given spectrum of modes, including large scale avalanches during events with a number of large amplitude modes. One of the most fundamental properties of ideal magnetohydrodynamics is the condition that plasma motion cannot change magnetic topology. The conventional representation of ideal magnetohydrodynamic modes by perturbing a toroidal equilibrium field through {delta}B-vector={nabla} Multiplication-Sign ({xi}-vector Multiplication-Sign B-vector), however, perturbs the magnetic topology, introducing extraneous magnetic islands in the field. A proper treatment of an ideal perturbation involves a full Lagrangian displacement of the field due to the perturbation and conserves magnetic topology as it should. In order to examine the effect of ideal magnetohydrodynamic modes on particle trajectories, the guiding center equations should include a correct Lagrangian treatment. Guiding center equations for an ideal displacement {xi}-vector are derived which preserve the magnetic topology and are used to examine mode particle resonances in toroidal confinement devices. These simulations are compared to others which are identical in all respects except that they use the linear representation for the field. Unlike the case for the magnetic field, the use of the linear field perturbation in the guiding center equations does not result in extraneous mode particle resonances.
Guiding Center Equations for Ideal Magnetohydrodynamic Modes
Roscoe B. White
2013-02-21
Guiding center simulations are routinely used for the discovery of mode-particle resonances in tokamaks, for both resistive and ideal instabilities and to find modifications of particle distributions caused by a given spectrum of modes, including large scale avalanches during events with a number of large amplitude modes. One of the most fundamental properties of ideal magnetohydrodynamics is the condition that plasma motion cannot change magnetic topology. The conventional representation of ideal magnetohydrodynamic modes by perturbing a toroidal equilibrium field through δ~B = ∇ X (ξ X B) however perturbs the magnetic topology, introducing extraneous magnetic islands in the field. A proper treatment of an ideal perturbation involves a full Lagrangian displacement of the field due to the perturbation and conserves magnetic topology as it should. In order to examine the effect of ideal magnetohydrodynamic modes on particle trajectories the guiding center equations should include a correct Lagrangian treatment. Guiding center equations for an ideal displacement ξ are derived which perserve the magnetic topology and are used to examine mode particle resonances in toroidal confinement devices. These simulations are compared to others which are identical in all respects except that they use the linear representation for the field. Unlike the case for the magnetic field, the use of the linear field perturbation in the guiding center equations does not result in extraneous mode particle resonances.
NASA Astrophysics Data System (ADS)
de Souza, R. A.; de Avellar, M. G. B.; Horvath, J. E.
2015-11-01
An appraisal of the behavior of stellar entropy along stellar evolution is made. It is shown that the entropy per baryon of a star of a fixed baryon number decreases monotonically with increasing compactness of the star. The same entropy per baryon increases only whenever an irreversible collapse of the star happens. The recent proposals for a gravitational entropy related to curvature may justify the huge increase of the entropy in the ultimate collapse to a black hole.
NASA Astrophysics Data System (ADS)
Zhang, Feng-Hui; Li, Li-Fang; Han, Zhan-Wen
In this paper we review the study of stellar population synthesis. So far there exist three methods in the study of the integrated light of stellar population-trail-and error, automated, and evolutionary population synthesis (EPS). We have discussed advantages and disadvantages for these methods. Among the three methods the EPS is the most direct approach to model galaxies. In this scheme, the model builder starts with knowledge of stellar evolution and attempts to build a model galaxy with physical input parameters such as star formation rate (SFR) and the initial mass function (IMF) slope. Therefore we have discussed emphatically the EPS method. First we have described and given the often used grids of several key ingredients in the EPS studies: (1) the library of evolutionary tracks used to calculate isochrones in the color-magnitude diagram (CMD), (2) the libraries of spectra adopted, which include empirical and theoretical stellar spectral libraries, star cluster library, active galactic nuclear (AGN) library and galaxy library, to derive the integrated spectral energy distributions (ISED) or magnitudes and colors in the suitable passbands, (3) the IMF used to evaluate the relative proportions of stars in the various evolutionary phases, and (4) the assumption for the underlying star formation rate (SFR) and chemical enrichment. Then we have listed several population synthesis criterions, i.e. broadband color indices, the integrated spectral energy distribution (ISED) and narrow band color indices, given the basic method of calculating broadband colors and flux-distribution for a simple stellar population (SSP). At last we have discussed simply the existed limitations, which are caused by some uncertainties in its two principal building blocks: stellar evolution models and spectral libraries in the studies of the EPS. Stellar evolution models are often subject to limitations in the following areas: the atomic data (radiative opacities, heavy element mixture
BOOK REVIEW: Stellarator and Heliotron Devices
NASA Astrophysics Data System (ADS)
Johnson, John L.
1999-02-01
ground for students without detracting from the usefulness of the book for knowledgeable fusion physicists. After a short, somewhat historical, introduction, Chapter 2 contains a good treatment of the basic properties of a toroidal magnetic configuration (the concepts of magnetic surfaces, rotational transform, shear and magnetic wells), averaging techniques which can often be used to simplify the calculations, helically invariant configurations, magnetic islands and line tracing techniques. Derivations and discussions of the basic tools of plasma theory, including the Vlasov equation, magnetohydrodynamic equations and their reduced form for low-β, large aspect ratio systems, properties of MHD waves, the drift kinetic equation and transport equations, are given in Chapter 3. Chapter 4 contains a good treatment of MHD equilibria, including a derivation of the three dimensional Grad-Shafranov equation, a discussion of the calculation of equilibria with a planar magnetic axis with both averaged equations and a variational approach, a comparison of the results of the two techniques, a formulation for stellarators with a helical magnetic axis and a good discussion of the Pfirsch-Schlüter current. The treatment of MHD instabilities in Chapter 5 is also excellent. It starts with a good derivation and discussion of the energy principle, gives a detailed treatment of ballooning modes where the wavelengths of the perturbation perpendicular to the field are short while those along B are long and derives the Mercier criterion from the ballooning mode equation. I personally prefer to obtain this criterion by making the low mode number assumption that dξ/dΨ>>dξ/dθ approx dξ/dζ, since non-ideal effects such as finite gyration radius corrections may provide less stabilization to these modes. A careful treatment of the resistive interchange mode is followed by a discussion of the role of localized stability criteria in the analysis of experiment and design studies, a study of
PREFACE: A Stellar Journey A Stellar Journey
NASA Astrophysics Data System (ADS)
Asplund, M.
2008-10-01
The conference A Stellar Journey was held in Uppsala, Sweden, 23 27June 2008, in honour of Professor Bengt Gustafsson's 65th birthday. The choice of Uppsala as the location for this event was obvious given Bengt's long-standing association with the city stemming back to his school days. With the exception of a two-year postdoc stint in Copenhagen, five years as professor at Stockholm University and two years as director of the Sigtuna foundation, Bengt has forged his illustrious professional career at Uppsala University. The symposium venue was Museum Gustavianum, once the main building of the oldest university in Scandinavia. The title of the symposium is a paraphrasing of Bengt's popular astronomy book Kosmisk Resa (in English: Cosmic Journey) written in the early eighties. I think this aptly symbolizes his career that has been an astronomical voyage from near to far, from the distant past to the present. The original book title was modified slightly to reflect that most of his work to date has dealt with stars in one way or another. In addition it also gives credit to Bengt's important role as a guiding light for a very large number of students, colleagues and collaborators, indeed for several generations of astronomers. For me personally, the book Kosmisk Resa bears particular significance as it has shaped my life rather profoundly. Although I had already decided to become an astronomer, when I first read the book as a 14-year-old I made up my mind then and there that I would study under Bengt Gustafsson and work on stars. Indeed I have remained true to this somewhat audacious resolution. I suspect that a great number of us have similar stories how Bengt has had a major influence on our lives, whether on the professional or personal level. Perhaps Bengt's most outstanding characteristic is his enthralling enthusiasm. This is equally true whether he is pondering some scientific conundrum, supervising students or performing in front of an audience, be it an
NASA Astrophysics Data System (ADS)
Chilingarian, Igor; Zolotukhin, Ivan; Beletsky, Yuri; Worthey, Guy
2015-08-01
Stellar libraries are fundamental tools required to understand stellar populations in star clusters and galaxies as well as properties of individual stars. Comprehensive libraries exist in the optical domain, but the near-infrared (NIR) domain stays a couple of decades behind. Here we present the Las Campanas Stellar Library project aiming at obtaining high signal-to-noise intermediate-resolution (R=8000) NIR spectra (0.83<λ<2.5μm) for a sample of 1200 stars in the Southern sky using the Folded-port InfraRed Echelette spectrograph at the 6.5-m Magellan Baade telescope. We developed a dedicated observing strategy and customized the telescope control software in order to achieve the highest possible level of data homogeniety. As of 2015, we observed about 600 stars of all spectral types and luminosity classes making our library the largest homogeneous collection of stellar spectra covering the entire NIR domain. We also re-calibrated in flux and wavelength the two existing optical stellar libraries, INDO-US and UVES-POP and followed up about 400 non-variable stars in the NIR in order to get complete optical-NIR coverage. Worth mentioning that our current sample includes about 80 AGB stars and a few dozens of bulge/LMC/SMC stars.
NASA Astrophysics Data System (ADS)
González Delgado, Rosa M.
2004-11-01
LLAGN that include low-ionization nuclear emission-line regions (LINERs), and transition-type objects (TOs) represent the most common type of nuclear activity. Here, a study of the central stellar population of LLAGN is presented. Our goal is to search for spectroscopic signatures of young and intermediate age stars, and to investigate their relationship with the ionization mechanism in LLAGN. The method used is based on the stellar population synthesis of the UV-optical continuum of the innermost (20-100 pc) regions in these galaxies. Half of the LINERs and TOs of the Palomar catalogue are analysed. It is found that weak-[OI] ([OI]/Hα≤0.25) LLAGN have an intermediate age stellar population that dominates the optical light. But young stellar clusters dominate the UV continuum in these objects. These clusters can co-exist with a black-hole in spatial scales of a few pc. Most of the strong-[OI] LLAGN have a predominantly old stellar population. These results suggest that young and intermediate age stars do not play a significant role in the ionization of LLAGN with strong [OI].
The infinite interface limit of multiple-region relaxed magnetohydrodynamics
Dennis, G. R.; Dewar, R. L.; Hole, M. J.; Hudson, S. R.
2013-03-15
We show the stepped-pressure equilibria that are obtained from a generalization of Taylor relaxation known as multi-region, relaxed magnetohydrodynamics (MRXMHD) are also generalizations of ideal magnetohydrodynamics (ideal MHD). We show this by proving that as the number of plasma regions becomes infinite, MRXMHD reduces to ideal MHD. Numerical convergence studies illustrating this limit are presented.
Magnetohydrodynamic modelling of exploding foil initiators
NASA Astrophysics Data System (ADS)
Neal, William
2015-06-01
Magnetohydrodynamic (MHD) codes are currently being developed, and used, to predict the behaviour of electrically-driven flyer-plates. These codes are of particular interest to the design of exploding foil initiator (EFI) detonators but there is a distinct lack of comparison with high-fidelity experimental data. This study aims to compare a MHD code with a collection of temporally and spatially resolved diagnostics including PDV, dual-axis imaging and streak imaging. The results show the code's excellent representation of the flyer-plate launch and highlight features within the experiment that the model fails to capture.
Lithium-sulfur hexafluoride magnetohydrodynamic power system
Dobran, F.
1987-02-24
A method is described to operate a two-phase flow magnetohydrodynamic electric power generation system with liquid lithium and gaseous sulfur-hexafluoride flowing through a diverging channel, with side electrodes to remove the electric current generated in the flowing liquid lithium, across the applied magnetic field that is perpendicular to both the flow velocity and electrodes. Sulfur-hexafluoride is dispersed in the form of small bubbles and reacts with liquid lithium that forms a continuous phase to conduct the current between the electrodes so as to produce a near isothermal two-phase flow mixture and provides for an expansion of lithium across the magnetic field in the generator.
Thermoelectric Magnetohydrodynamic Stirring of Liquid Metals
NASA Astrophysics Data System (ADS)
Jaworski, M. A.; Gray, T. K.; Antonelli, M.; Kim, J. J.; Lau, C. Y.; Lee, M. B.; Neumann, M. J.; Xu, W.; Ruzic, D. N.
2010-03-01
The direct observation of a thermoelectric magnetohydrodynamic (TEMHD) flow has been achieved and is reported here. The origin of the flow is identified based on a series of qualitative tests and corresponds, quantitatively, with a swirling flow TEMHD model. A theory for determining the dominant driver of a free-surface flow, TEMHD or thermocapillary (TC), is found to be consistent with the experimental results. The use of the analytical form for an open geometry develops a new dimensionless parameter describing the ratio of TEMHD to TC generated flows.
Magnetohydrodynamic effects in liquid metal batteries
NASA Astrophysics Data System (ADS)
Stefani, F.; Galindo, V.; Kasprzyk, C.; Landgraf, S.; Seilmayer, M.; Starace, M.; Weber, N.; Weier, T.
2016-07-01
Liquid metal batteries (LMBs) consist of two liquid metal electrodes and a molten salt ionic conductor sandwiched between them. The density ratios allow for a stable stratification of the three layers. LMBs were already considered as part of energy conversion systems in the 1960s and have recently received renewed interest for economical large-scale energy storage. In this paper, we concentrate on the magnetohydrodynamic aspects of this cell type with special focus on electro-vortex flows and possible effects of the Tayler instability.
Intrinsic rotation of toroidally confined magnetohydrodynamics.
Morales, Jorge A; Bos, Wouter J T; Schneider, Kai; Montgomery, David C
2012-10-26
The spatiotemporal self-organization of viscoresistive magnetohydrodynamics in a toroidal geometry is studied. Curl-free toroidal magnetic and electric fields are imposed. It is observed in our simulations that a flow is generated, which evolves from dominantly poloidal to toroidal when the Lundquist numbers are increased. It is shown that this toroidal organization of the flow is consistent with the tendency of the velocity field to align with the magnetic field. Up-down asymmetry of the geometry causes the generation of a nonzero toroidal angular momentum. PMID:23215195
Magnetohydrodynamic turbulence in the solar wind
NASA Technical Reports Server (NTRS)
Goldstein, Melvyn L.
1995-01-01
The fluctuations in magnetic field and plasma velocity in solar wind, which possess many features of fully developed magnetohydrodynamic (MHD) turbulence, are discussed. Direct spacecraft observations from 0.3 to over 20 AU, remote sensing radio scintillation observations, numerical simulations, and various models provide complementary methods that show that the fluctuations in the wind parameters undergo significant dynamical evolution independent of whatever turbulence might exist in the solar photosphere and corona. The Cluster mission, with high time resolution particle and field measurements and its variable separation strategies, should be able to provide data for answering many questions on MHD turbulence.
Magnetohydrodynamic thermochemotherapy and MRI of mouse tumors
NASA Astrophysics Data System (ADS)
Brusentsov, Nikolay A.; Brusentsova, Tatiana N.; Filinova, Elena Yu.; Jurchenko, Nikolay Y.; Kupriyanov, Dmitry A.; Pirogov, Yuri A.; Dubina, Andry I.; Shumskikh, Maxim N.; Shumakov, Leonid I.; Anashkina, Ekaterina N.; Shevelev, Alexandr A.; Uchevatkin, Andry A.
2007-04-01
A dextran-ferrite magnetic fluid was successfully tested as magnetic resonance imaging (MRI) contrast agent. The same magnetic fluid was then combined with Melphalan, a chemotherapeutic drug, and used for magnetohydrodynamic thermochemotherapy of different tumors. The placement of the tumors in an AC magnetic field led to hyperthermia at 46 °C for 30 min. In combination with tumor slime aspiration, a 30% regression of ˜130 mm 3 non-metastatic P388 tumors in BDF 1 mice was reached, together with a life span increase of 290%. The same procedure associated with cyclophosphamide treatment of ˜500 mm 3 metastases tumor increased the animal's life span by 180%.
Broken symmetry in ideal magnetohydrodynamic turbulence
NASA Technical Reports Server (NTRS)
Shebalin, John V.
1993-01-01
A numerical study of the long-time evolution of a number of cases of inviscid, isotropic, incompressible, three-dimensional fluid, and magneto-fluid turbulence has been completed. The results confirm that ideal magnetohydrodynamic turbulence is non-ergodic if there is no external magnetic field present. This is due essentially to a canonical symmetry being broken in an arbitrary dynamical representation. The broken symmetry manifests itself as a coherent structure, i.e., a non-zero time-averaged part of the turbulent magnetic field. The coherent structure is observed, in one case, to contain about eighteen percent of the total energy.
Stellar Astrophysical Fluid Dynamics
NASA Astrophysics Data System (ADS)
Thompson, Michael J.; Christensen-Dalsgaard, Jørgen
2003-05-01
Preface; 1. A selective overview Jørgen Christensen-Dalsgaard and Michael J. Thompson; Part I. Stellar Convection and Oscillations: 2. On the diversity of stellar pulsations Wojciech A. Dziembowski; 3. Acoustic radiation and mode excitation by turbulent convection Günter Houdek; 4. Understanding roAp stars Margarida S. Cunha; 5. Waves in the magnetised solar atmosphere Colin S. Rosenthal; Part II. Stellar Rotation and Magnetic Fields: 6. Stellar rotation: a historical survey Leon Mestel; 7. The oscillations of rapidly rotating stars Michel Rieutord; 8. Solar tachocline dynamics: eddy viscosity, anti-friction, or something in between? Michael E. McIntyre; 9. Dynamics of the solar tachocline Pascale Garaud; 10. Dynamo processes: the interaction of turbulence and magnetic fields Michael Proctor; 11. Dynamos in planets Chris Jones; Part III. Physics and Structure of Stellar Interiors: 12. Solar constraints on the equation of state Werner Däppen; 13. 3He transport and the solar neutrino problem Chris Jordinson; 14. Mixing in stellar radiation zones Jean-Paul Zahn; 15. Element settling and rotation-induced mixing in slowly rotating stars Sylvie Vauclair; Part IV. Helio- and Asteroseismology: 16. Solar structure and the neutrino problem Hiromoto Shibahashi; 17. Helioseismic data analysis Jesper Schou; 18. Seismology of solar rotation Takashi Sekii; 19. Telechronohelioseismology Alexander Kosovichev; Part V. Large-Scale Numerical Experiments: 20. Bridges between helioseismology and models of convection zone dynamics Juri Toomre; 21. Numerical simulations of the solar convection zone Julian R. Elliott; 22. Modelling solar and stellar magnetoconvection Nigel Weiss; 23. Nonlinear magnetoconvection in the presence of a strong oblique field Keith Julien, Edgar Knobloch and Steven M. Tobias; 24. Simulations of astrophysical fluids Marcus Brüggen; Part VI. Dynamics: 25. A magic electromagnetic field Donald Lynden-Bell; 26. Continuum equations for stellar dynamics Edward A
Stellar Astrophysical Fluid Dynamics
NASA Astrophysics Data System (ADS)
Thompson, Michael J.; Christensen-Dalsgaard, Jørgen
2008-02-01
Preface; 1. A selective overview Jørgen Christensen-Dalsgaard and Michael J. Thompson; Part I. Stellar Convection and Oscillations: 2. On the diversity of stellar pulsations Wojciech A. Dziembowski; 3. Acoustic radiation and mode excitation by turbulent convection Günter Houdek; 4. Understanding roAp stars Margarida S. Cunha; 5. Waves in the magnetised solar atmosphere Colin S. Rosenthal; Part II. Stellar Rotation and Magnetic Fields: 6. Stellar rotation: a historical survey Leon Mestel; 7. The oscillations of rapidly rotating stars Michel Rieutord; 8. Solar tachocline dynamics: eddy viscosity, anti-friction, or something in between? Michael E. McIntyre; 9. Dynamics of the solar tachocline Pascale Garaud; 10. Dynamo processes: the interaction of turbulence and magnetic fields Michael Proctor; 11. Dynamos in planets Chris Jones; Part III. Physics and Structure of Stellar Interiors: 12. Solar constraints on the equation of state Werner Däppen; 13. 3He transport and the solar neutrino problem Chris Jordinson; 14. Mixing in stellar radiation zones Jean-Paul Zahn; 15. Element settling and rotation-induced mixing in slowly rotating stars Sylvie Vauclair; Part IV. Helio- and Asteroseismology: 16. Solar structure and the neutrino problem Hiromoto Shibahashi; 17. Helioseismic data analysis Jesper Schou; 18. Seismology of solar rotation Takashi Sekii; 19. Telechronohelioseismology Alexander Kosovichev; Part V. Large-Scale Numerical Experiments: 20. Bridges between helioseismology and models of convection zone dynamics Juri Toomre; 21. Numerical simulations of the solar convection zone Julian R. Elliott; 22. Modelling solar and stellar magnetoconvection Nigel Weiss; 23. Nonlinear magnetoconvection in the presence of a strong oblique field Keith Julien, Edgar Knobloch and Steven M. Tobias; 24. Simulations of astrophysical fluids Marcus Brüggen; Part VI. Dynamics: 25. A magic electromagnetic field Donald Lynden-Bell; 26. Continuum equations for stellar dynamics Edward A
PREFACE: A Stellar Journey A Stellar Journey
NASA Astrophysics Data System (ADS)
Asplund, M.
2008-10-01
The conference A Stellar Journey was held in Uppsala, Sweden, 23 27June 2008, in honour of Professor Bengt Gustafsson's 65th birthday. The choice of Uppsala as the location for this event was obvious given Bengt's long-standing association with the city stemming back to his school days. With the exception of a two-year postdoc stint in Copenhagen, five years as professor at Stockholm University and two years as director of the Sigtuna foundation, Bengt has forged his illustrious professional career at Uppsala University. The symposium venue was Museum Gustavianum, once the main building of the oldest university in Scandinavia. The title of the symposium is a paraphrasing of Bengt's popular astronomy book Kosmisk Resa (in English: Cosmic Journey) written in the early eighties. I think this aptly symbolizes his career that has been an astronomical voyage from near to far, from the distant past to the present. The original book title was modified slightly to reflect that most of his work to date has dealt with stars in one way or another. In addition it also gives credit to Bengt's important role as a guiding light for a very large number of students, colleagues and collaborators, indeed for several generations of astronomers. For me personally, the book Kosmisk Resa bears particular significance as it has shaped my life rather profoundly. Although I had already decided to become an astronomer, when I first read the book as a 14-year-old I made up my mind then and there that I would study under Bengt Gustafsson and work on stars. Indeed I have remained true to this somewhat audacious resolution. I suspect that a great number of us have similar stories how Bengt has had a major influence on our lives, whether on the professional or personal level. Perhaps Bengt's most outstanding characteristic is his enthralling enthusiasm. This is equally true whether he is pondering some scientific conundrum, supervising students or performing in front of an audience, be it an
Global magnetohydrodynamic simulations on multiple GPUs
NASA Astrophysics Data System (ADS)
Wong, Un-Hong; Wong, Hon-Cheng; Ma, Yonghui
2014-01-01
Global magnetohydrodynamic (MHD) models play the major role in investigating the solar wind-magnetosphere interaction. However, the huge computation requirement in global MHD simulations is also the main problem that needs to be solved. With the recent development of modern graphics processing units (GPUs) and the Compute Unified Device Architecture (CUDA), it is possible to perform global MHD simulations in a more efficient manner. In this paper, we present a global magnetohydrodynamic (MHD) simulator on multiple GPUs using CUDA 4.0 with GPUDirect 2.0. Our implementation is based on the modified leapfrog scheme, which is a combination of the leapfrog scheme and the two-step Lax-Wendroff scheme. GPUDirect 2.0 is used in our implementation to drive multiple GPUs. All data transferring and kernel processing are managed with CUDA 4.0 API instead of using MPI or OpenMP. Performance measurements are made on a multi-GPU system with eight NVIDIA Tesla M2050 (Fermi architecture) graphics cards. These measurements show that our multi-GPU implementation achieves a peak performance of 97.36 GFLOPS in double precision.
Global Magnetohydrodynamic Modeling of the Solar Corona
NASA Technical Reports Server (NTRS)
Linker, Jon A.
1997-01-01
Under this contract, we have continued our investigations of the large scale structure of the solar corona and inner heliosphere using global magnetohydrodynamic (MHD) simulations. These computations have also formed the basis for studies of coronal mass ejections (CMES) using realistic coronal configurations. We have developed a technique for computing realistic magnetohydrodynamic (MHD) computations of the solar corona and inner heliosphere. To perform computations that can be compared with specific observations, it is necessary to incorporate solar observations into the boundary conditions. We have used the Wilcox Solar Observatory synoptic maps (collected during a solar rotation by daily measurements of the line-of-sight magnetic field at central meridian) to specify the radial magnetic field (B,) at the photosphere. For the initial condition, we use a potential magnetic field consistent with the specified distribution of B, at the lower boundary, and a wind solution consistent with the specified plasma density and temperature at the solar surface. Together this initial condition forms a (non-equilibrium) approximation of the state of the solar corona for the time-dependent MHD computation. The MHD equations are then integrated in time to steady state. Here we describe solutions relevant to a recent solar eclipse, as well as Ulysses observations. We have also developed a model configuration of solar minimum, useful for studying CME initiation and propagation.
Double-duct liquid metal magnetohydrodynamic engine
Haaland, Carsten M.
1995-01-01
An internal combustion, liquid metal (LM) magnetohydrodynamic (MHD) engine and an alternating current (AC) magnetohydrodynamic generator, are used in combination to provide useful AC electric energy output. The engine design has-four pistons and a double duct configuration, with each duct containing sodium potassium liquid metal confined between free pistons located at either end of the duct. The liquid metal is forced to flow back and forth in the duct by the movement of the pistons, which are alternatively driven by an internal combustion process. In the MHD generator, the two LM-MHD ducts pass in close proximity through a Hartmann duct with output transformer. AC power is produced by operating the engine with the liquid metal in the two generator ducts always flowing in counter directions. The amount of liquid metal maintained in the ducts may be varied. This provides a variable stroke length for the pistons. The engine/generator provides variable AC power at variable frequencies that correspond to the power demands of the vehicular propulsion. Also the engine should maintain nearly constant efficiency throughout the range of power usage. Automobiles and trucks could be powered by the invention, with no transmission or power converter devices being required.
Double-duct liquid metal magnetohydrodynamic engine
Haaland, Carsten M.
1997-01-01
An internal combustion, liquid metal (LM) magnetohydrodynamic (MHD) engine and an alternating current (AC) magnetohydrodynamic generator, are used in combination to provide useful AC electric energy output. The engine design has four pistons and a double duct configuration, with each duct containing sodium potassium liquid metal confined between free pistons located at either end of the duct. The liquid metal is forced to flow back and forth in the duct by the movement of the pistons, which are alternatively driven by an internal combustion process. In the MHD generator, the two LM-MHD ducts pass in close proximity through a Hartmann duct with output transformer. AC power is produced by operating the engine with the liquid metal in the two generator ducts always flowing in counter directions. The amount of liquid metal maintained in the ducts may be varied. This provides a variable stroke length for the pistons. The engine/generator provides variable AC power at variable frequencies that correspond to the power demands of the vehicular propulsion. Also the engine should maintain nearly constant efficiency throughout the range of power usage. Automobiles and trucks could be powered by the invention, with no transmission or power converter devices being required.
Magnetohydrodynamic boundary conditions for global models
NASA Technical Reports Server (NTRS)
Forbes, T. G.
1988-01-01
Boundary conditions in the ionosphere and the upstream solar wind are important in determining the dynamics of global magnetohydrodynamic models of the magnetosphere. It is generally recognized that the orientation of the magnetic field in the upstream solar wind strongly modulates the rate of energy input into the magnetosphere by magnetic reconnection. However, other aspects of the upstream boundary conditions may determine whether the reconnection occurs in a patchy manner, as in flux transfer events, or in a global manner, as in the Paschmann et al. (1979) events. Ionospheric boundary conditions should also affect the reconnection process. For example, ionospheric line-tying can cause x-line motion in the outer magnetosphere. If it is assumed that auroras occur on field lines mapping to x-lines, then auroral motions are different than the local convective motion of the plasma in which they occur. Global magnetohydrodynamic models which incorporate both magnetospheric reconnection and ionospheric convection could be used to investigate the effect of reconnection and convection upon dayside and nightside auroral motions during the course of a magnetic substorm.
Analytical study of magnetohydrodynamic propulsion stability
NASA Astrophysics Data System (ADS)
Abdollahzadeh Jamalabadi, M. Y.
2014-09-01
In this paper an analytical solution for the stability of the fully developed flow drive in a magneto-hydro-dynamic pump with pulsating transverse Eletro-magnetic fields is presented. To do this, a theoretical model of the flow is developed and the analytical results are obtained for both the cylindrical and Cartesian configurations that are proper to use in the propulsion of marine vessels. The governing parabolic momentum PDEs are transformed into an ordinary differential equation using approximate velocity distribution. The numerical results are obtained and asymptotic analyses are built to discover the mathematical behavior of the solutions. The maximum velocity in a magneto-hydro-dynamic pump versus time for various values of the Stuart number, electro-magnetic interaction number, Reynolds number, aspect ratio, as well as the magnetic and electrical angular frequency and the shift of the phase angle is presented. Results show that for a high Stuart number there is a frequency limit for stability of the fluid flow in a certain direction of the flow. This stability frequency is dependent on the geometric parameters of a channel.
Efficient acceleration of relativistic magnetohydrodynamic jets
NASA Astrophysics Data System (ADS)
Toma, Kenji; Takahara, Fumio
2013-08-01
Relativistic jets in active galactic nuclei, galactic microquasars, and gamma-ray bursts are widely considered to be magnetohydrodynamically driven by black hole accretion systems, although the conversion mechanism from the Poynting into the particle kinetic energy flux is still open. Recent detailed numerical and analytical studies of global structures of steady, axisymmetric magnetohydrodynamic (MHD) flows with specific boundary conditions have not reproduced as rapid an energy conversion as required by observations. In order to find more suitable boundary conditions, we focus on the flow along a poloidal magnetic field line just inside the external boundary, without treating the transfield force balance in detail. We find some examples of the poloidal field structure and corresponding external pressure profile for an efficient and rapid energy conversion as required by observations, and that the rapid acceleration requires a rapid decrease of the external pressure above the accretion disk. We also clarify the differences between the fast magnetosonic point of the MHD flow and the sonic point of the de Laval nozzle.
Numerical Simulations of Radiative Shocked Clumps: Application to the Young Stellar Outflow HH 29
NASA Astrophysics Data System (ADS)
Palmer, Jacob Nordin
Herbig-Haro objects are the nebulous emission regions arising from shocks associated with the bipolar outflows of young stellar objects. Observations of the emission lines of HH objects can reveal properties of these stellar jets, which are associated with star formation, and one tool to study stellar jets and HH objects is computer modeling. AstroBEAR is a 3D magnetohydrodynamic code with adaptive mesh refinement and parallelization capabilities. The goals of this research were (1) to improve the cooling in AstroBEAR and add the ability to create valid emission-line maps for HH objects, and (2) to model a jet overrunning and shocking clumps, and to explore a parameter space that may lead to morphologies like those of HH 29. The results generate emission-line maps with the correct spatial morphology in Halpha and [SII], and begin to map out the regimes (shock velocities, resolution, physical scales) in which the code remains valid.
Rogers, F J
1998-09-17
New efforts to calculate opacity have produced significant improvements in the quality of stellar models. The most dramatic effect has been large opacity enhancements for stars subject to large amplitude pulsations. Significant improvement in helioseismic modeling has also been obtained. A description and comparisons of the new opacity efforts are give
Progress Toward Attractive Stellarators
Neilson, G H; Brown, T G; Gates, D A; Lu, K P; Zarnstorff, M C; Boozer, A H; Harris, J H; Meneghini, O; Mynick, H E; Pomphrey, N; Reiman, A H; Xanthopoulos, P
2011-01-05
The quasi-axisymmetric stellarator (QAS) concept offers a promising path to a more compact stellarator reactor, closer in linear dimensions to tokamak reactors than previous stellarator designs. Concept improvements are needed, however, to make it more maintainable and more compatible with high plant availability. Using the ARIES-CS design as a starting point, compact stellarator designs with improved maintenance characteristics have been developed. While the ARIES-CS features a through-the-port maintenance scheme, we have investigated configuration changes to enable a sector-maintenance approach, as envisioned for example in ARIES AT. Three approaches are reported. The first is to make tradeoffs within the QAS design space, giving greater emphasis to maintainability criteria. The second approach is to improve the optimization tools to more accurately and efficiently target the physics properties of importance. The third is to employ a hybrid coil topology, so that the plasma shaping functions of the main coils are shared more optimally, either with passive conductors made of high-temperature superconductor or with local compensation coils, allowing the main coils to become simpler. Optimization tools are being improved to test these approaches.
ERIC Educational Resources Information Center
Ros, Rosa M.
2009-01-01
The main purpose of the stellar demonstrator is to help explain the movement of stars. In particular, students have difficulties understanding why, if they are living in the Northern Hemisphere, they may observe starts in the Southern Hemisphere, or why circumpolar stars are not the same in different parts of Europe. Using the demonstrator, these…
Stellar Ontogeny: From Dust...
ERIC Educational Resources Information Center
MOSAIC, 1978
1978-01-01
Discusses the process of star formation. Infrared and radio astronomy, particularly microwave astronomy is used to provide information on different stages of stellar formation. The role of dust and gas which swirl through the interstellar regions of a galaxy and the collapse of a cloud in star formation are also presented. (HM)
Trends of stellar entropy along stellar evolution
NASA Astrophysics Data System (ADS)
de Avellar, Guilherme Bronzato, Marcio; Alvares de Souza, Rodrigo; Horvath, Jorge Ernesto
2016-02-01
This paper is devoted to discussing the difference in the thermodynamic entropy budget per baryon in each type of stellar object found in the Universe. We track and discuss the actual decrease of the stored baryonic thermodynamic entropy from the most primitive molecular cloud up to the final fate of matter in black holes, passing through evolved states of matter as found in white dwarfs and neutron stars. We then discuss the case of actual stars with different masses throughout their evolution, clarifying the role of the virial equilibrium condition for the decrease in entropy and related issues. Finally, we discuss the role of gravity in driving the composition and the structural changes of stars with different Main Sequence masses during their evolution up to the final product. Particularly, we discuss the entropy of a black hole in this context arguing that the dramatic increase in its entropy, differently from the other cases, is due to the gravitational field itself.
NASA Astrophysics Data System (ADS)
Bastian, T. S.; Dulk, G. A.; Bookbinder, J. A.
The dynamic spectrum, a three dimensional record of the radio intensity as a function both of time and frequency, has long been used as a probe of plasma processes in the solar corona. Beginning with the work of Wild and McCready (1950) dynamic spectroscopy has been used to distinguish between the multitude of radio wave emitting phenomena which occur in the solar corona and to infer the physical mechanisms responsible. Stellar dynamic spectroscopy has always been a tantalizing prospect. The vast body of experience with solar dynamic spectroscopy would prove invaluable in interpreting stellar dynamic spectra. Further, the new parameter regimes presented by stellar coronas would allow further insight to be gained in the physical processes at work in stellar coronas. Recently, Bastian and Bookbinder (1987) used the Very Large Array
Stellar Pulsations and Stellar Evolution: Conflict, Cohabitation, or Symbiosis?
NASA Astrophysics Data System (ADS)
Weiss, Achim
While the analysis of stellar pulsations allows the determination of current properties of a star, stellar evolution models connect it with its previous history. In many cases results from both methods do not agree. In this review some classical and current cases of disagreement are presented. In some cases these conflicts led to an improvement of the theory of stellar evolution, while in others they still remain unsolved. Some well-known problems of stellar physics are pointed out as well, for which it is hoped that seismology—or in general the analysis of stellar pulsations—will help to resolve them. The limits of this symbiosis will be discussed as well.
Quantifying the impact of stellar feedback on molecular clouds
NASA Astrophysics Data System (ADS)
Boyden, Ryan; Koch, Eric; Offner, Stella
2016-01-01
Although the origin of molecular cloud turbulence remains debated, one possibility is that stellar feedback injects enough energy to drive observed motions on parsec scales. To investigate this possibility, we use magnetohydrodynamic simulations where we vary the stellar mass-loss rates and magnetic field strength. We generate synthetic 12CO(1-0) maps assuming that the simulation is at the distance of the nearby Perseus molecular cloud. By comparing different initial conditions and evolutionary times, we are able to identify differences in our synthetic observations. Using a variety of statistical techniques proposed in the literature, we quantify the differences by calculating the first, second, and higher order moment maps of the data cubes, analyzing the power spectrum, and convolving the data with Gaussian wavelets. We find that many turbulent statistics, such as the spectral correlation function, principal component analysis, and delta-variance, are sensitive to changes in mass-loss rates and/or turbulent structure. This demonstrates that stellar feedback influences molecular cloud turbulence and may be characterized using certain statistical metrics.
STELLAR DYNAMOS AND CYCLES FROM NUMERICAL SIMULATIONS OF CONVECTION
Dubé, Caroline; Charbonneau, Paul E-mail: paulchar@astro.umontreal.ca
2013-09-20
We present a series of kinematic axisymmetric mean-field αΩ dynamo models applicable to solar-type stars, for 20 distinct combinations of rotation rates and luminosities. The internal differential rotation and kinetic helicity profiles required to calculate source terms in these dynamo models are extracted from a corresponding series of global three-dimensional hydrodynamical simulations of solar/stellar convection, so that the resulting dynamo models end up involving only one free parameter, namely, the turbulent magnetic diffusivity in the convecting layers. Even though the αΩ dynamo solutions exhibit a broad range of morphologies, and sometimes even double cycles, these models manage to reproduce relatively well the observationally inferred relationship between cycle period and rotation rate. On the other hand, they fail in capturing the observed increase of magnetic activity levels with rotation rate. This failure is due to our use of a simple algebraic α-quenching formula as the sole amplitude-limiting nonlinearity. This suggests that α-quenching is not the primary mechanism setting the amplitude of stellar magnetic cycles, with magnetic reaction on large-scale flows emerging as the more likely candidate. This inference is coherent with analyses of various recent global magnetohydrodynamical simulations of solar/stellar convection.
A Multi-period Oscillation in a Stellar Superflare
NASA Astrophysics Data System (ADS)
Pugh, C. E.; Nakariakov, V. M.; Broomhall, A.-M.
2015-11-01
Flares that are orders of magnitude larger than the most energetic solar flares are routinely observed on Sun-like stars, raising the question of whether the same physical processes are responsible for both solar and stellar flares. In this Letter, we present a white-light stellar superflare on the star KIC 9655129, observed by NASA’s Kepler mission, with a rare multi-period quasi-periodic pulsation (QPP) pattern. Two significant periodic processes were detected using the wavelet and autocorrelation techniques, with periods of 78 ± 12 minutes and 32 ± 2 minutes. By comparing the phases and decay times of the two periodicities, the QPP signal was found to most likely be linear, suggesting that the two periodicities are independent, possibly corresponding either to different magnetohydrodynamic (MHD) modes of the flaring region or different spatial harmonics of the same mode. The presence of multiple periodicities is a good indication that the QPPs were caused by MHD oscillations and suggests that the physical processes in operation during stellar flares could be the same as those in solar flares.
Stellar radio emission (Review)
NASA Astrophysics Data System (ADS)
Zhelezniakov, V. V.
The current understanding of the radio-emission characteristics of 'ordinary' main sequence stars as well as giants and supergiants is examined. Particular consideration is given to radio emission from supergiants, Young T Tauri stars, magnetic Ap stars, flare stars of UV Ceti type, Alpha Sco, and RS CVn objects. It is noted that the study of stellar radio emission is in its initial stage. Further progress in this area depends on successes in finding new radio sources, associated, for example, with magnetic stars, and on an intensified investigation of the frequency spectra and polarization of already-discovered radio stars. It is also noted that, although the current knowledge of solar physics can help in understanding stellar radio emission, models and ideas developed for solar conditions should not be mechanically transferred to other stars by a simple change in scale.
Exploring Astrophysical Magnetohydrodynamics in the Laboratory
NASA Astrophysics Data System (ADS)
Manuel, Mario
2014-10-01
Plasma evolution in many astrophysical systems is dominated by magnetohydrodynamics. Specifically of interest to this talk are collimated outflows from accretion systems. Away from the central object, the Euler equations can represent the plasma dynamics well and may be scaled to a laboratory system. We have performed experiments to investigate the effects of a background magnetic field on an otherwise hydrodynamically collimated plasma. Laser-irradiated, cone targets produce hydrodynamically collimated plasma jets and a pulse-powered solenoid provides a constant background magnetic field. The application of this field is shown to completely disrupt the original flow and a new magnetically-collimated, hollow envelope is produced. Results from these experiments and potential implications for their astrophysical analogs will be discussed.
Dynamo Effect in the Kraichnan Magnetohydrodynamic Turbulence
NASA Astrophysics Data System (ADS)
Arponen, Heikki; Horvai, Peter
2007-10-01
The existence of a dynamo effect in a simplified magnetohydrodynamic model of turbulence is considered when the magnetic Prandtl number approaches zero or infinity. The magnetic field is interacting with an incompressible Kraichnan-Kazantsev model velocity field which incorporates also a viscous cutoff scale. An approximate system of equations in the different scaling ranges can be formulated and solved, so that the solution tends to the exact one when the viscous and magnetic-diffusive cutoffs approach zero. In this approximation we are able to determine analytically the conditions for the existence of a dynamo effect and give an estimate of the dynamo growth rate. Among other things we show that in the large magnetic Prandtl number case the dynamo effect is always present. Our analytical estimates are in good agreement with previous numerical studies of the Kraichnan-Kazantsev dynamo by Vincenzi (J. Stat. Phys. 106:1073-1091, 2002).
Ideal magnetohydrodynamic interchanges in low density plasmas
Huang Yimin; Goel, Deepak; Hassam, A.B.
2005-03-01
The ideal magnetohydrodynamic equations are usually derived under the assumption V{sub A}<
Magnetohydrodynamic stability of a compound liquid jet
NASA Astrophysics Data System (ADS)
Radwan, Ahmed E.
1989-10-01
The magnetohydrodynamics (MHD) stability of a compound nonmiscible fluid jet is discussed. A general eigenvalue relation, for that model which involves the fluid inertia, capillarity and electromagnetic forces, is derived. The model is capillary unstable only for small axisymmetric disturbances and stable for the rest. The magnetic fields interior and exterior to the gas-mantle jet have always a stabilizing influence. The radii ratio of the concentric jets plays an important role in the (instability) stability states and are (decreasing) increasing with increasing magnetic field intensity as the exterior radius is much larger than the interior radius; under some restrictions of the radii ratio and above a certain value of the magnetic field the capillary instability is omitted and completely suppressed and then stability sets in. The latter result is verified analytically and confirmed numerically in the case in which the cylindrical surface of the outer jet is sited at infinity.
Rarefaction wave in relativistic steady magnetohydrodynamic flows
Sapountzis, Konstantinos Vlahakis, Nektarios
2014-07-15
We construct and analyze a model of the relativistic steady-state magnetohydrodynamic rarefaction that is induced when a planar symmetric flow (with one ignorable Cartesian coordinate) propagates under a steep drop of the external pressure profile. Using the method of self-similarity, we derive a system of ordinary differential equations that describe the flow dynamics. In the specific limit of an initially homogeneous flow, we also provide analytical results and accurate scaling laws. We consider that limit as a generalization of the previous Newtonian and hydrodynamic solutions already present in the literature. The model includes magnetic field and bulk flow speed having all components, whose role is explored with a parametric study.
Numerical Methods for Radiation Magnetohydrodynamics in Astrophysics
Klein, R I; Stone, J M
2007-11-20
We describe numerical methods for solving the equations of radiation magnetohydrodynamics (MHD) for astrophysical fluid flow. Such methods are essential for the investigation of the time-dependent and multidimensional dynamics of a variety of astrophysical systems, although our particular interest is motivated by problems in star formation. Over the past few years, the authors have been members of two parallel code development efforts, and this review reflects that organization. In particular, we discuss numerical methods for MHD as implemented in the Athena code, and numerical methods for radiation hydrodynamics as implemented in the Orion code. We discuss the challenges introduced by the use of adaptive mesh refinement in both codes, as well as the most promising directions for future developments.
Entropy generation analysis of magnetohydrodynamic induction devices
NASA Astrophysics Data System (ADS)
Salas, Hugo; Cuevas, Sergio; López de Haro, Mariano
1999-10-01
Magnetohydrodynamic (MHD) induction devices such as electromagnetic pumps or electric generators are analysed within the approach of entropy generation. The flow of an electrically-conducting incompressible fluid in an MHD induction machine is described through the well known Hartmann model. Irreversibilities in the system due to ohmic dissipation, flow friction and heat flow are included in the entropy-generation rate. This quantity is used to define an overall efficiency for the induction machine that considers the total loss caused by process irreversibility. For an MHD generator working at maximum power output with walls at constant temperature, an optimum magnetic field strength (i.e. Hartmann number) is found based on the maximum overall efficiency.
Nonideal magnetohydrodynamic instabilities and toroidal magnetic confinement
Furth, H.P.
1985-05-01
The marked divergence of experimentally observed plasma instability phenomena from the predictions of ideal magnetohydrodynamics led in the early 1960s to the formulations of finite-resistivity stability theory. Beginning in the 1970s, advanced plasma diagnostics have served to establish a detailed correspondence between the predictions of the finite-resistivity theory and experimental plasma behavior - particularly in the case of the resistive kink mode and the tokamak plasma. Nonlinear resistive-kink phenomena have been found to govern the transport of magnetic flux and plasma energy in the reversed-field pinch. The other predicted finite-resistivity instability modes have been more difficult to identify directly and their implications for toroidal magnetic confinement are still unresolved.
Multidimensional numerical scheme for resistive relativistic magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Komissarov, Serguei S.
2007-12-01
The paper describes a new upwind conservative numerical scheme for special relativistic resistive magnetohydrodynamics with scalar resistivity. The magnetic field is kept approximately divergence free and the divergence of the electric field is kept consistent with the electric charge distribution via the method of Generalized Lagrange Multiplier. The hyperbolic fluxes are computed using the Harten-Lax-van Leer (HLL) prescription and the source terms are accounted via the time-splitting technique. The results of test simulations show that the scheme can handle equally well both resistive current sheets and shock waves, and thus can be a useful tool for studying phenomena of relativistic astrophysics that involve both colliding supersonic flows and magnetic reconnection.
Extended inertial range phenomenology of magnetohydrodynamic turbulence
NASA Technical Reports Server (NTRS)
Matthaeus, William H.; Zhou, YE
1989-01-01
A phenomenological treatment of the inertial range of isotropic statistically steady magnetohydrodynamic turbulence is presented, extending the theory of Kraichnan (1965). The role of Alfven wave propagation is treated on equal footing with nonlinear convection, leading to a simple generalization of the relations between the times characteristic of wave propagation, convection, energy transfer, and decay of triple correlations. The theory leads to a closed-form steady inertial range spectral law that reduces to the Kraichnan and Kolmogorov laws in appropriate limits. The Kraichnan constant is found to be related in a simple way to the Kolmogorov constant; for typical values of the latter constant, the former has values in the range 1.22-1.87. Estimates of the time scale associated with spectral transfer of energy also emerge from the new approach, generalizing previously presented 'golden rules' for relating the spectral transfer time scale to the Alfven and eddy-turnover time scales.
Resistive Magnetohydrodynamic Simulations of Relativistic Magnetic Reconnection
NASA Technical Reports Server (NTRS)
Zenitani, Seiji; Hesse, Michael; Klimas, Alex
2010-01-01
Resistive relativistic magnetohydrodynamic (RRMHD) simulations are applied to investigate the system evolution of relativistic magnetic reconnection. A time-split Harten-Lan-van Leer method is employed. Under a localized resistivity, the system exhibits a fast reconnection jet with an Alfv enic Lorentz factor inside a narrow Petschek-type exhaust. Various shock structures are resolved in and around the plasmoid such as the post-plasmoid vertical shocks and the "diamond-chain" structure due to multiple shock reflections. Under a uniform resistivity, Sweet-Parker-type reconnection slowly evolves. Under a current-dependent resistivity, plasmoids are repeatedly formed in an elongated current sheet. It is concluded that the resistivity model is of critical importance for RRMHD modeling of relativistic magnetic reconnection.
General Relativistic Magnetohydrodynamic Simulations of Collapsars
NASA Technical Reports Server (NTRS)
Mizuno, Yosuke; Yamada, S.; Koider, S.; Shipata, K.
2005-01-01
We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of collapsars including a rotating black hole. Initially, we assume that the core collapse has failed in this star. A rotating black hole of a few solar masses is inserted by hand into the calculation. The simulation results show the formation of a disklike structure and the generation of a jetlike outflow near the central black hole. The jetlike outflow propagates and accelerated mainly by the magnetic field. The total jet velocity is approximately 0.3c. When the rotation of the black hole is faster, the magnetic field is twisted strongly owing to the frame-dragging effect. The magnetic energy stored by the twisting magnetic field is directly converted to kinetic energy of the jet rather than propagating as an Alfven wave. Thus, as the rotation of the black hole becomes faster, the poloidal velocity of the jet becomes faster.
Lattice kinetic simulation of nonisothermal magnetohydrodynamics.
Chatterjee, Dipankar; Amiroudine, Sakir
2010-06-01
In this paper, a lattice kinetic algorithm is presented to simulate nonisothermal magnetohydrodynamics in the low-Mach number incompressible limit. The flow and thermal fields are described by two separate distribution functions through respective scalar kinetic equations and the magnetic field is governed by a vector distribution function through a vector kinetic equation. The distribution functions are only coupled via the macroscopic density, momentum, magnetic field, and temperature computed at the lattice points. The novelty of the work is the computation of the thermal field in conjunction with the hydromagnetic fields in the lattice Boltzmann framework. A 9-bit two-dimensional (2D) lattice scheme is used for the numerical computation of the hydrodynamic and thermal fields, whereas the magnetic field is simulated in a 5-bit 2D lattice. Simulation of Hartmann flow in a channel provides excellent agreement with corresponding analytical results. PMID:20866540
Numerical models for high beta magnetohydrodynamic flow
Brackbill, J.U.
1987-01-01
The fundamentals of numerical magnetohydrodynamics for highly conducting, high-beta plasmas are outlined. The discussions emphasize the physical properties of the flow, and how elementary concepts in numerical analysis can be applied to the construction of finite difference approximations that capture these features. The linear and nonlinear stability of explicit and implicit differencing in time is examined, the origin and effect of numerical diffusion in the calculation of convective transport is described, and a technique for maintaining solenoidality in the magnetic field is developed. Many of the points are illustrated by numerical examples. The techniques described are applicable to the time-dependent, high-beta flows normally encountered in magnetically confined plasmas, plasma switches, and space and astrophysical plasmas. 40 refs.
Hall magnetohydrodynamics: Conservation laws and Lyapunov stability
NASA Astrophysics Data System (ADS)
Holm, Darryl D.
1987-05-01
Hall electric fields produce circulating mass flow in confined ideal-fluid plasmas. The conservation laws, Hamiltonian structure, equilibrium state relations, and Lyapunov stability conditions are presented here for ideal Hall magnetohydrodynamics (HMHD) in two and three dimensions. The approach here is to use the remarkable array of nonlinear conservation laws for HMHD that follow from its Hamiltonian structure in order to construct explicit Lyapunov functionals for the HMHD equilibrium states. In this way, the Lyapunov stability analysis provides classes of HMHD equilibria that are stable and whose linearized initial-value problems are well posed (in the sense of possessing continuous dependence on initial conditions). Several examples are discussed in both two and three dimensions.
Magnetohydrodynamic Modeling of the Jovian Magnetosphere
NASA Technical Reports Server (NTRS)
Walker, Raymond
2005-01-01
Under this grant we have undertaken a series of magnetohydrodynamic (MHD) simulation and data analysis studies to help better understand the configuration and dynamics of Jupiter's magnetosphere. We approached our studies of Jupiter's magnetosphere in two ways. First we carried out a number of studies using our existing MHD code. We carried out simulation studies of Jupiter s magnetospheric boundaries and their dependence on solar wind parameters, we studied the current systems which give the Jovian magnetosphere its unique configuration and we modeled the dynamics of Jupiter s magnetosphere following a northward turning of the interplanetary magnetic field (IMF). Second we worked to develop a new simulation code for studies of outer planet magnetospheres.
Nuclear magnetohydrodynamic EMP, solar storms, and substorms
Rabinowitz, M. ); Meliopoulous, A.P.S.; Glytsis, E.N. . School of Electrical Engineering); Cokkinides, G.J. )
1992-10-20
In addition to a fast electromagnetic pulse (EMP), a high altitude nuclear burst produces a relatively slow magnetohydrodynamic EMP (MHD EMP), whose effects are like those from solar storm geomagnetically induced currents (SS-GIC). The MHD EMP electric field E [approx lt] 10[sup [minus] 1] V/m and lasts [approx lt] 10[sup 2] sec, whereas for solar storms E [approx gt] 10[sup [minus] 2] V/m and lasts [approx gt] 10[sup 3] sec. Although the solar storm electric field is lower than MHD EMP, the solar storm effects are generally greater due to their much longer duration. Substorms produce much smaller effects than SS-GIC, but occur much more frequently. This paper describes the physics of such geomagnetic disturbances and analyzes their effects.
Remarkable connections between extended magnetohydrodynamics models
Lingam, M. Morrison, P. J. Miloshevich, G.
2015-07-15
Through the use of suitable variable transformations, the commonality of all extended magnetohydrodynamics (MHD) models is established. Remarkable correspondences between the Poisson brackets of inertialess Hall MHD and inertial MHD (which has electron inertia, but not the Hall drift) and extended MHD (which has both effects) are established. The helicities (two in all) for each of these models are obtained through these correspondences. The commonality of all the extended MHD models is traced to the existence of two Lie-dragged 2-forms, which are closely associated with the canonical momenta of the two underlying species. The Lie-dragging of these 2-forms by suitable velocities also leads to the correct equations of motion. The Hall MHD Poisson bracket is analyzed in detail, the Jacobi identity is verified through a detailed proof, and this proof ensures the Jacobi identity for the Poisson brackets of all the models.
Symmetry transforms for ideal magnetohydrodynamics equilibria.
Bogoyavlenskij, Oleg I
2002-11-01
A method for constructing ideal magnetohydrodynamics (MHD) equilibria is introduced. The method consists of the application of symmetry transforms to any known MHD equilibrium [ O. I. Bogoyavlenskij, Phys. Rev. E. 62, 8616, (2000)]. The transforms break the geometrical symmetries of the field-aligned solutions and produce continuous families of the nonsymmetric MHD equilibria. The method of symmetry transforms also allows to obtain MHD equilibria with current sheets and exact solutions with noncollinear vector fields B and V. A model of the nonsymmetric astrophysical jets outside of their accretion disks is developed. The total magnetic and kinetic energy of the jet is finite in any layer c(1)
A photolithographic fabrication technique for magnetohydrodynamic micropumps
NASA Astrophysics Data System (ADS)
Kuenstner, Stephen; Baylor, Martha-Elizabeth
2014-03-01
Magnetohydrodynamic (MHD) devices use perpendicular electric and magnetic fields to exert a Lorentz body force on a conducting fluid. Miniaturized MHD devices have been used to create pumps, stirrers, heat exchangers, and microfluidic networks. Compared to mechanical micropumps, MHD micropumps are appealing because they require no moving parts, which simplifies fabrication, and because they are amenable to electronic control. This abstract reports the fabrication and testing of a centimeter-scale MHD pump using a thiol-ene/methacrylate-based photopolymer and mask-based photolithographic technique. Pumps like this one could simplify the fabrication of sophisticated optofluidic devices, including liquid-core, liquid cladding (L2) waveguides, which are usually created with PDMS using stamps, or etched into silicon wafers. The photolithographic technique demonstrated here requires only one masking step to create fluid channels with complex geometries.
MAGNETOHYDRODYNAMICS OF THE WEAKLY IONIZED SOLAR PHOTOSPHERE
Cheung, Mark C. M.; Cameron, Robert H.
2012-05-01
We investigate the importance of ambipolar diffusion and Hall currents for high-resolution comprehensive ({sup r}ealistic{sup )} photospheric simulations. To do so, we extended the radiative magnetohydrodynamics code MURaM to use the generalized Ohm's law under the assumption of local thermodynamic equilibrium. We present test cases comparing analytical solutions with numerical simulations for validation of the code. Furthermore, we carried out a number of numerical experiments to investigate the impact of these neutral-ion effects in the photosphere. We find that, at the spatial resolutions currently used (5-20 km per grid point), the Hall currents and ambipolar diffusion begin to become significant-with flows of 100 m s{sup -1} in sunspot light bridges, and changes of a few percent in the thermodynamic structure of quiet-Sun magnetic features. The magnitude of the effects is expected to increase rapidly as smaller-scale variations are resolved by the simulations.
Action principles for extended magnetohydrodynamic models
Keramidas Charidakos, I.; Lingam, M.; Morrison, P. J.; White, R. L.; Wurm, A.
2014-09-15
The general, non-dissipative, two-fluid model in plasma physics is Hamiltonian, but this property is sometimes lost or obscured in the process of deriving simplified (or reduced) two-fluid or one-fluid models from the two-fluid equations of motion. To ensure that the reduced models are Hamiltonian, we start with the general two-fluid action functional, and make all the approximations, changes of variables, and expansions directly within the action context. The resulting equations are then mapped to the Eulerian fluid variables using a novel nonlocal Lagrange-Euler map. Using this method, we recover Lüst's general two-fluid model, extended magnetohydrodynamic (MHD), Hall MHD, and electron MHD from a unified framework. The variational formulation allows us to use Noether's theorem to derive conserved quantities for each symmetry of the action.
Scale-locality of magnetohydrodynamic turbulence
Aluie, Hussein; Eyink, Gregory L
2009-01-01
We investigate the scale-locality of cascades of conserved invariants at high kinetic and magnetic Reynolds numbers in the 'inertial-inductive range' of magnetohydrodynamic (MHD) turbulence, where velocity and magnetic field increments exhibit suitable power-law scaling. We prove that fluxes of total energy and cross-helicity - or, equivalently, fluxes of Elsaesser energies - are dominated by the contributions of local triads. Corresponding spectral transfers are also scale-local when defined using octave wavenumber bands. Flux and transfer of magnetic helicity may be dominated by nonlocal triads. The magnetic stretching term also may be dominated by non-local triads but we prove that it can convert energy only between velocity and magnetic modes at comparable scales. We explain the disagreement with numerical studies that have claimed conversion non locally between disparate scales. We present supporting data from a 1024{sup 3} simulation of forced MHD turbulence.