NASA Astrophysics Data System (ADS)
Zhai, Cuili; Zhang, Ting
2016-09-01
In this article, we consider the global existence and uniqueness of the solution to the 2D incompressible non-resistive MHD system with non-equilibrium background magnetic field. Our result implies that a strong enough non-equilibrium background magnetic field will guarantee the stability of the nonlinear MHD system. Beside the classical energy method, the interpolation inequalities and the algebraic structure of the equations coming from the incompressibility of the fluid are crucial in our arguments.
On Regularity Criteria for the 2D Generalized MHD System
NASA Astrophysics Data System (ADS)
Jiang, Zaihong; Wang, Yanan; Zhou, Yong
2016-06-01
This paper deals with the problem of regularity criteria for the 2D generalized MHD system with fractional dissipative terms {-Λ^{2α}u} for the velocity field and {-Λ^{2β}b} for the magnetic field respectively. Various regularity criteria are established to guarantee smoothness of solutions. It turns out that our regularity criteria imply previous global existence results naturally.
2-D skin-current toroidal-MHD-equilibrium code
Feinberg, B.; Niland, R.A.; Coonrod, J.; Levine, M.A.
1982-09-01
A two-dimensional, toroidal, ideal MHD skin-current equilibrium computer code is described. The code is suitable for interactive implementation on a minicomptuer. Some examples of the use of the code for design and interpretation of toroidal cusp experiments are presented.
Atmospheric Outflows from Hot Jupiters: 2D MHD Simulations
NASA Astrophysics Data System (ADS)
Uribe, A.; Matsakos, T.; Konigl, A.
2015-01-01
Recent observations of stellar hydrogen Ly-α line absorption during transits of some hot Jupiter exoplanets suggest the presence of a dense, fast wind that is blowing from planetary atmosphere tep{2003Natur.422..143V,2007ApJ...671L..61B}. Modeling efforts include 1D hydrodynamic models tep{2009ApJ...693...23M,2004Icar..170..167Y,2007P&SS...55.1426G} and 2D isothermal magnetized wind models tep{2014arXiv1404.5817T}, among others. In this work, we model the 2D structure of the irradiated upper atmosphere of a hot Jupiter planet and its interaction with the planetary magnetic field. We calculate self consistently the heating by stellar UV radiation and the cooling of the atmosphere by Ly-α emission. We solve for the ionization structure assuming a 100% hydrogen atmosphere, accounting for the radiative ionization, recombination and advection of the gas. We show the effect of stellar tides and planetary magnetic field on the planet outflow and calculate the Ly-α transmission spectra of the resulting atmosphere.
Global small solutions of 2-D incompressible MHD system
NASA Astrophysics Data System (ADS)
Lin, Fanghua; Xu, Li; Zhang, Ping
2015-11-01
In this paper, we consider the global wellposedness of 2-D incompressible magneto-hydrodynamical system with smooth initial data which is close to some non-trivial steady state. It is a coupled system between the Navier-Stokes equations and a free transport equation with a universal nonlinear coupling structure. The main difficulty of the proof lies in exploring the dissipative mechanism of the system. To achieve this and to avoid the difficulty of propagating anisotropic regularity for the free transport equation, we first reformulate our system (1.1) in the Lagrangian coordinates (2.19). Then we employ anisotropic Littlewood-Paley analysis to establish the key a prioriL1 (R+ ; Lip (R2)) estimate for the Lagrangian velocity field Yt. With this estimate, we can prove the global wellposedness of (2.19) with smooth and small initial data by using the energy method. We emphasize that the algebraic structure of (2.19) is crucial for the proofs to work. The global wellposedness of the original system (1.1) then follows by a suitable change of variables.
Problems in nonlinear resistive MHD
Turnbull, A.D.; Strait, E.J.; La Haye, R.J.; Chu, M.S.; Miller, R.L.
1998-12-31
Two experimentally relevant problems can relatively easily be tackled by nonlinear MHD codes. Both problems require plasma rotation in addition to the nonlinear mode coupling and full geometry already incorporated into the codes, but no additional physics seems to be crucial. These problems discussed here are: (1) nonlinear coupling and interaction of multiple MHD modes near the B limit and (2) nonlinear coupling of the m/n = 1/1 sawtooth mode with higher n gongs and development of seed islands outside q = 1.
Implicit adaptive mesh refinement for 2D reduced resistive magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Philip, Bobby; Chacón, Luis; Pernice, Michael
2008-10-01
An implicit structured adaptive mesh refinement (SAMR) solver for 2D reduced magnetohydrodynamics (MHD) is described. The time-implicit discretization is able to step over fast normal modes, while the spatial adaptivity resolves thin, dynamically evolving features. A Jacobian-free Newton-Krylov method is used for the nonlinear solver engine. For preconditioning, we have extended the optimal "physics-based" approach developed in [L. Chacón, D.A. Knoll, J.M. Finn, An implicit, nonlinear reduced resistive MHD solver, J. Comput. Phys. 178 (2002) 15-36] (which employed multigrid solver technology in the preconditioner for scalability) to SAMR grids using the well-known Fast Adaptive Composite grid (FAC) method [S. McCormick, Multilevel Adaptive Methods for Partial Differential Equations, SIAM, Philadelphia, PA, 1989]. A grid convergence study demonstrates that the solver performance is independent of the number of grid levels and only depends on the finest resolution considered, and that it scales well with grid refinement. The study of error generation and propagation in our SAMR implementation demonstrates that high-order (cubic) interpolation during regridding, combined with a robustly damping second-order temporal scheme such as BDF2, is required to minimize impact of grid errors at coarse-fine interfaces on the overall error of the computation for this MHD application. We also demonstrate that our implementation features the desired property that the overall numerical error is dependent only on the finest resolution level considered, and not on the base-grid resolution or on the number of refinement levels present during the simulation. We demonstrate the effectiveness of the tool on several challenging problems.
A pressure-based high resolution numerical method for resistive MHD
NASA Astrophysics Data System (ADS)
Xisto, Carlos M.; Páscoa, José C.; Oliveira, Paulo J.
2014-10-01
In the paper we describe in detail a numerical method for the resistive magnetohydrodynamic (MHD) equations involving viscous flow and report the results of application to a number of typical MHD test cases. The method is of the finite volume type but mixes aspects of pressure-correction and density based solvers; the algorithm arrangement is patterned on the well-known PISO algorithm, which is a pressure method, while the flux computation makes use of the AUSM-MHD scheme, which originates from density based methods. Five groups of test cases are addressed to verify and validate the method. We start with two resistive MHD cases, namely the Shercliff and Hunt flow problems, which are intended to validate the method for low-speed resistive MHD flows. The remaining three test cases, namely the cloud-shock interaction, the MHD rotor and the MHD blast wave, are standard 2D ideal MHD problems that serve to validate the method under high-speed flow and complex interaction of MHD shocks. Finally, we demonstrate the method with a more complex application problem, and discuss results of simulation for a quasi-bi-dimensional self-field magnetoplasmadynamic (MPD) thruster, for which we study the effect of cathode length upon the electromagnetic nozzle performance.
OpenMHD: Godunov-type code for ideal/resistive magnetohydrodynamics (MHD)
NASA Astrophysics Data System (ADS)
Zenitani, Seiji
2016-04-01
OpenMHD is a Godunov-type finite-volume code for ideal/resistive magnetohydrodynamics (MHD). It is written in Fortran 90 and is parallelized by using MPI-2 and OpenMP. The code was originally developed for studying magnetic reconnection problems and has been made publicly available in the hope that others may find it useful.
First results from ideal 2-D MHD reconstruction: magnetopause reconnection event seen by Cluster
NASA Astrophysics Data System (ADS)
Teh, W.-L.; Ã-. Sonnerup, B. U.
2008-09-01
We have applied a new reconstruction method (Sonnerup and Teh, 2008), based on the ideal single-fluid MHD equations in a steady-state, two-dimensional geometry, to a reconnection event observed by the Cluster-3 (C3) spacecraft on 5 July 2001, 06:23 UT, at the dawn-side Northern-Hemisphere magnetopause. The event has been previously studied by use of Grad-Shafranov (GS) reconstruction, performed in the deHoffmann-Teller frame, and using the assumption that the flow effects were either negligible or the flow was aligned with the magnetic field. Our new method allows the reconstruction to be performed in the frame of reference moving with the reconnection site (the X-line). In the event studied, this motion is tailward/equatorward at 140 km/s. The principal result of the study is that the new method functions well, generating a magnetic field map that is qualitatively similar to those obtained in the earlier GS-based reconstructions but now includes the reconnection site itself. In comparison with the earlier map by Hasegawa et al. (2004), our new map has a slightly improved ability (cc=0.979 versus cc=0.975) to predict the fields measured by the other three Cluster spacecraft, at distances from C3 ranging from 2132 km (C1) to 2646 km (C4). The new field map indicates the presence of a magnetic X-point, located some 5300 km tailward/equatorward of C3 at the time of its traversal of the magnetopause. In the immediate vicinity of the X-point, the ideal-MHD assumption breaks down, i.e. resistive and/or other effects should be included. We have circumvented this problem by an ad-hoc procedure in which we allow the axial part of convection electric field to be non-constant near the reconnection site. The new reconstruction method also provides a map of the velocity field, in which the inflow into the wedge of reconnected field lines and the plasma jet within it can be seen, and maps of the electric potential and of the electric current distribution. Even though the
Spectrum of resistive MHD modes in cylindrical plasmas
Ryu, C.M.; Grimm, R.C.
1983-07-01
A numerical study of the normal modes of a compressible resistive MHD fluid in cylindrical geometry is presented. Resistivity resolves the shear Alfven and slow magnetosonic continua of ideal MHD into discrete spectra and gives rise to heavily damped modes whose frequencies lie on specific lines in the complex plane. Fast magnetosonic waves are less affected but are also damped. Overstable modes arise from the shear Alfven spectrum. The stabilizing effect of favorable average curvature is shown. Eigenfunctions illustrating the nature of typical normal modes are displayed.
Towards a Scalable Fully-Implicit Fully-coupled Resistive MHD Formulation with Stabilized FE Methods
Shadid, J N; Pawlowski, R P; Banks, J W; Chacon, L; Lin, P T; Tuminaro, R S
2009-06-03
This paper presents an initial study that is intended to explore the development of a scalable fully-implicit stabilized unstructured finite element (FE) capability for low-Mach-number resistive MHD. The discussion considers the development of the stabilized FE formulation and the underlying fully-coupled preconditioned Newton-Krylov nonlinear iterative solver. To enable robust, scalable and efficient solution of the large-scale sparse linear systems generated by the Newton linearization, fully-coupled algebraic multilevel preconditioners are employed. Verification results demonstrate the expected order-of-acuracy for the stabilized FE discretization of a 2D vector potential form for the steady and transient solution of the resistive MHD system. In addition, this study puts forth a set of challenging prototype problems that include the solution of an MHD Faraday conduction pump, a hydromagnetic Rayleigh-Bernard linear stability calculation, and a magnetic island coalescence problem. Initial results that explore the scaling of the solution methods are presented on up to 4096 processors for problems with up to 64M unknowns on a CrayXT3/4. Additionally, a large-scale proof-of-capability calculation for 1 billion unknowns for the MHD Faraday pump problem on 24,000 cores is presented.
Energy structure of MHD flow coupling with outer resistance circuit
NASA Astrophysics Data System (ADS)
Huang, Z. Y.; Liu, Y. J.; Chen, Y. Q.; Peng, Z. L.
2015-08-01
Energy structure of MHD flow coupling with outer resistance circuit is studied to illuminate qualitatively and quantitatively the energy relation of this basic MHD flow system with energy input and output. Energy structure are analytically derived based on the Navier-Stocks equations for two-dimensional fully-developed flow and generalized Ohm's Law. The influences of applied magnetic field, Hall parameter and conductivity on energy structure are discussed based on the analytical results. Associated energies in MHD flow are deduced and validated by energy conservation. These results reveal that energy structure consists of two sub structures: electrical energy structure and internal energy structure. Energy structure and its sub structures provide an integrated theoretical energy path of the MHD system. Applied magnetic field and conductivity decrease the input energy, dissipation by fluid viscosity and internal energy but increase the ratio of electrical energy to input energy, while Hall parameter has the opposite effects. These are caused by their different effects on Bulk velocity, velocity profiles, voltage and current in outer circuit. Understanding energy structure helps MHD application designers to actively adjust the allocation of different parts of energy so that it is more reasonable and desirable.
2D Radiation MHD K-shell Modeling of Single Wire Array Stainless Steel Experiments on the Z Machine
Thornhill, J. W.; Giuliani, J. L.; Apruzese, J. P.; Chong, Y. K.; Davis, J.; Dasgupta, A.; Whitney, K. G.; Clark, R. W.; Jones, B.; Coverdale, C. A.; Ampleford, D. J.; Cuneo, M. E.; Deeney, C.
2009-01-21
Many physical effects can produce unstable plasma behavior that affect K-shell emission from arrays. Such effects include: asymmetry in the initial density profile, asymmetry in power flow, thermal conduction at the boundaries, and non-uniform wire ablation. Here we consider how asymmetry in the radiation field also contributes to the generation of multidimensional plasma behavior that affects K-shell power and yield. To model this radiation asymmetry, we have incorporated into the MACH2 r-z MHD code a self-consistent calculation of the non-LTE population kinetics based on radiation transport using multi-dimensional ray tracing. Such methodology is necessary for modeling the enhanced radiative cooling that occurs at the anode and cathode ends of the pinch during the run-in phase of the implosion. This enhanced radiative cooling is due to reduced optical depth at these locations producing an asymmetric flow of radiative energy that leads to substantial disruption of large initial diameter (>5 cm) pinches and drives 1D into 2D fluid (i.e., Rayleigh-Taylor like) flows. The impact of this 2D behavior on K-shell power and yield is investigated by comparing 1D and 2D model results with data obtained from a series of single wire array stainless steel experiments performed on the Z generator.
Evolution of generalized two-dimensional magnetotail equilibria in ideal and resistive MHD
NASA Astrophysics Data System (ADS)
Merkin, V. G.; Sitnov, M. I.; Lyon, J. G.
2015-03-01
We present results of two-dimensional (2-D) magnetohydrodynamic (MHD) simulations of the terrestrial magnetotail. A regional adaptation of the Lyon-Fedder-Mobarry global MHD model is used. As initial conditions, we employ a class of asymptotic magnetotail equilibria with and without an accumulation of magnetic flux at the tailward end (a Bz hump). The former have been recently shown by full particle simulations to be unstable to a kinetic mode with formal properties of ion tearing. Thus, our goal here is to investigate the evolution of the same equilibria in the MHD approximation and assist in the physical interpretation of the kinetic simulations. This is additionally motivated by the energy principle considerations which suggest that if the system is unstable kinetically, it may also be unstable ideally. To seek dynamical MHD regimes similar to those observed in kinetic simulations, we implement two sets of boundary conditions (velocity balanced, VB, and momentum balanced, MB), one allowing plasma flows through the boundaries and the other inhibiting such flows. The use of more reflecting MB boundary conditions results in suppression of any significant dynamics, and we see no substantial changes beyond initial equilibrium relaxation. On the other hand, VB boundary conditions allow a more efficient relaxation of initial equilibrium and absorb subsequently generated plasma flows. With these boundary conditions we find the equilibrium without a flux accumulation (i.e., with constant magnetic field component normal to the current sheet) to develop an apparently resistive mode accompanied by tailward plasma flows. At the same time, the equilibria with a Bz hump of sufficiently large amplitude develop a different, ideal, mode characterized by spontaneous generation of earthward plasma flows and an exponential growth of the corresponding electric field. This growth is qualitatively similar to the corresponding fully kinetic simulations although no explosive growth of
Global well-posedness for the 2D MHD equations without magnetic diffusion in a strip domain
NASA Astrophysics Data System (ADS)
Ren, Xiaoxia; Xiang, Zhaoyin; Zhang, Zhifei
2016-04-01
We study the initial boundary value problem of two dimensional MHD equations without magnetic diffusion in a strip domain. It was proved that the MHD equations have a unique global strong solution around the equilibrium state ≤ft(0,{{\\mathbf{e}}1}\\right) for both the non-slip boundary condition and Navier slip boundary condition on the velocity.
Kinetic Effects of Energetic Particles on Resistive MHD Stability
Takahashi, R.; Brennan, D. P.; Kim, C. C.
2009-04-03
We show that the kinetic effects of energetic particles can play a crucial role in the stability of the m/n=2/1 tearing mode in tokamaks (e.g., JET, JT-60U, and DIII-D), where the fraction of energetic particle {beta}{sub frac} is high. Using model equilibria based on DIII-D experimental reconstructions, the nonideal MHD linear stability of cases unstable to the 2/1 mode is investigated including a {delta}f particle-in-cell model for the energetic particles coupled to the nonlinear 3D resistive MHD code NIMROD[C. C. Kim et al., Phys. Plasmas 15, 072507 (2008)]. It is observed that energetic particles have significant damping and stabilizing effects at experimentally relevant {beta}, {beta}{sub frac}, and S, and excite a real frequency of the 2/1 mode. Extrapolation of the results is discussed for implications to JET and ITER, where the effects are projected to be significant.
The ideal tearing mode: theory and resistive MHD simulations
NASA Astrophysics Data System (ADS)
Del Zanna, L.; Landi, S.; Papini, E.; Pucci, F.; Velli, M.
2016-05-01
Classical MHD reconnection theories, both the stationary Sweet-Parker model and the tearing instability, are known to provide rates which are too slow to explain the observations. However, a recent analysis has shown that there exists a critical threshold on current sheet's thickness, namely a/L ∼ S -1/3, beyond which the tearing modes evolve on fast macroscopic Alfvénic timescales, provided the Lunquist number S is high enough, as invariably found in solar and astrophysical plasmas. Therefore, the classical Sweet-Parker scenario, for which the diffusive region scales as a/L ∼ S -1/2 and thus can be up to ∼ 100 times thinner than the critical value, is likely to be never realized in nature, as the current sheet itself disrupts in the elongation process. We present here two-dimensional, compressible, resistive MHD simulations, with S ranging from 105 to 107, that fully confirm the linear analysis. Moreover, we show that a secondary plasmoid instability always occurs when the same critical scaling is reached on the local, smaller scale, leading to a cascading explosive process, reminiscent of the flaring activity.
Resistivity inversion in 2-D anisotropic media: numerical experiments
NASA Astrophysics Data System (ADS)
Wiese, Timothy; Greenhalgh, Stewart; Zhou, Bing; Greenhalgh, Mark; Marescot, Laurent
2015-04-01
Many rocks and layered/fractured sequences have a clearly expressed electrical anisotropy although it is rare in practice to incorporate anisotropy into resistivity inversion. In this contribution, we present a series of 2.5-D synthetic inversion experiments for various electrode configurations and 2-D anisotropic models. We examine and compare the image reconstructions obtained using the correct anisotropic inversion code with those obtained using the false but widely used isotropic assumption. Superior reconstruction in terms of reduced data misfit, true anomaly shape and position, and anisotropic background parameters were obtained when the correct anisotropic assumption was employed for medium to high coefficients of anisotropy. However, for low coefficient values the isotropic assumption produced better-quality results. When an erroneous isotropic inversion is performed on medium to high level anisotropic data, the images are dominated by patterns of banded artefacts and high data misfits. Various pole-pole, pole-dipole and dipole-dipole data sets were investigated and evaluated for the accuracy of the inversion result. The eigenvalue spectra of the pseudo-Hessian matrix and the formal resolution matrix were also computed to determine the information content and goodness of the results. We also present a data selection strategy based on high sensitivity measurements which drastically reduces the number of data to be inverted but still produces comparable results to that of the comprehensive data set. Inversion was carried out using transversely isotropic model parameters described in two different co-ordinate frames for the conductivity tensor, namely Cartesian versus natural or eigenframe. The Cartesian frame provided a more stable inversion product. This can be simply explained from inspection of the eigenspectra of the pseudo-Hessian matrix for the two model descriptions.
MHD modeling of magnetotail instability for localized resistivity
NASA Technical Reports Server (NTRS)
Hesse, Michael; Birn, Joachim
1994-01-01
We present results of a three-dimensional magnetohydrodynamic (MHD) simulation of magnetotail evolution initiated by a sudden occurrence or increase of spatially localized resistivity as the major expected concequence of some localized microinstability. Because of the absence of a quantitative model, possible variations of resistivity levels with current density, or the reduction thereof, are not incorporated in the present investigation. The emphasis of the study is on an investigation of the changes to the overall evolution brought about by this localization, in particular, on the disruption and diversion of the cross-tail current and the nonlinear evolution of the magnetotail instability. The immediate consequences of the occurrence of the localized resistance and the resulting electric field are a reduction and diversion of the electric current around the region of high resistivity, associated with an increase of B(sub z) ('dipolarization') at the earthward edge and a decrease of B(sub z) at the tailward edge of this region. These effects, however, are localized and do not involve a reduction of the total cross-tail current and hence do not lead to the global development of a 'substorm current wedge,' which includes not only the reduction of the cross-tail current but also the buildup of a global field-aligned current system of 'regional 1' type (toward the Earth on the dawnside and away on the duskside of the tail). Such signatures develop at a later time, as consequences of a three-dimensional tearing instability, which is triggered by the occurrence of the resistivity. These features are found in combination with plasmoid formation and ejection, quite similar to results of earlier simulations with uniform resistivity. Differences are found in the timescale of the evolution, which tends to be shorter for localized resistivity, and in the propagation of the dipolarization effects in the equatorial plane. Whereas for uniform resistivity the temporal increase in
ELRIS2D: A MATLAB Package for the 2D Inversion of DC Resistivity/IP Data
NASA Astrophysics Data System (ADS)
Akca, Irfan
2016-04-01
ELRIS2D is an open source code written in MATLAB for the two-dimensional inversion of direct current resistivity (DCR) and time domain induced polarization (IP) data. The user interface of the program is designed for functionality and ease of use. All available settings of the program can be reached from the main window. The subsurface is discretized using a hybrid mesh generated by the combination of structured and unstructured meshes, which reduces the computational cost of the whole inversion procedure. The inversion routine is based on the smoothness constrained least squares method. In order to verify the program, responses of two test models and field data sets were inverted. The models inverted from the synthetic data sets are consistent with the original test models in both DC resistivity and IP cases. A field data set acquired in an archaeological site is also used for the verification of outcomes of the program in comparison with the excavation results.
2D numerical simulation of the resistive reconnection layer
D. A. Uzdensky; R. M. Kulsrud
2000-07-21
In this paper the authors present a two-dimensional numerical simulation of a reconnection current layer in incompressible resistive magnetohydrodynamics with uniform resistivity in the limit of very large Lundquist numbers. They use realistic boundary conditions derived consistently from the outside magnetic field, and they also take into account the effect of the backpressure from flow into the separatrix region. They find that within a few Alfven times the system reaches a steady state consistent with the Sweet-Parker model, even if the initial state is Petschek-like.
2D Numerical Simulation of the Resistive Reconnection Layer
Kulsrud, R.M.; Uzdensky, D.A.
1999-03-01
In this paper we present a two-dimensional numerical simulation of a reconnection current layer in incompressible resistive magnetohydrodynamics with uniform resistivity in the limit of very large Lundquist numbers. We use realistic boundary conditions derived consistently from the outside magnetic field, and we also take into account the effect of the back pressure from flow into the separatrix region. We find that within a few Alfvén times the system reaches a steady state consistent with the Sweet-Parker model, even if the initial state is Petschek-like.
2D Resistive Magnetohydrodynamics Calculations with an Arbitrary Lagrange Eulerian Code
NASA Astrophysics Data System (ADS)
Rousculp, C. L.; Gianakon, T. A.; Lipnikov, K. N.; Nelson, E. M.
2015-11-01
Single fluid resistive MHD is useful for modeling Z-pinch configurations in cylindrical geometry. One such example is thin walled liners for shock physics or HEDP experiments driven by capacitor banks such as the LANL's PHELIX or Sandia-Z. MHD is also useful for modeling high-explosive-driven flux compression generators (FCGs) and their high-current switches. The resistive MHD in our arbitrary Lagrange Eulerian (ALE) code operates in one and two dimensions in both Cartesian and cylindrical geometry. It is implemented as a time-step split operator, which consists of, ideal MHD connected to the explicit hydro momentum and energy equations and a second order mimetic discretization solver for implicit solution of the magnetic diffusion equation. In a staggered grid scheme, a single-component of cell-centered magnetic flux is conserved in the Lagrangian frame exactly, while magnetic forces are accumulated at the nodes. Total energy is conserved to round off. Total flux is conserved under the ALE relaxation and remap. The diffusion solver consistently computes Ohmic heating. Both Neumann and Dirichlet boundary conditions are available with coupling to external circuit models. Example calculations will be shown.
NASA Astrophysics Data System (ADS)
Wu, C.; Chang, T.
2010-12-01
A new method in describing the multifractal characteristics of intermittent events was introduced by Cheng and Wu [Chang T. and Wu C.C., Physical Rev, E77, 045401(R), 2008]. The procedure provides a natural connection between the rank-ordered spectrum and the idea of one-parameter scaling for monofractals. This technique has been demonstrated using results obtained from a 2D MHD simulation. It has also been successfully applied to in-situ solar wind observations [Chang T., Wu, C.C. and Podesta, J., AIP Conf Proc. 1039, 75, 2008], and the broadband electric field oscillations from the auroral zone [Tam, S.W.Y. et al., Physical Rev, E81, 036414, 2010]. We take the next step in this procedure. By using the ROMA spectra and the scaled probability distribution functions (PDFs), raw PDFs can be calculated, which can be compared directly with PDFs from observations or simulation results. In addition to 2D MHD simulation results and in-situ solar wind observation, we show clearly using the ROMA analysis the multifractal character of the 3D fluid simulation data obtained from the JHU turbulence database cluster at http://turbulence.pha.jhu.edu. In particular, we show the scaling of the non-symmetrical PDF for the parallel-velocity fluctuations of this 3D fluid data.
A Numerical Study of Resistivity and Hall Effects for a Compressible MHD Model
NASA Technical Reports Server (NTRS)
Yee, H. C.; Sjogreen, B.
2005-01-01
The effect of resistive, Hall, and viscous terms on the flow structure compared with compressible ideal MHD is studied numerically for a one-fluid non-ideal MHD model. The goal of the present study is to shed some light on the emerging area of non-ideal MHD modeling and simulation. Numerical experiments are performed on a hypersonic blunt body flow with future application to plasma aerodynamics flow control in reentry vehicles. Numerical experiments are also performed on a magnetized time-developing mixing layer with possible application to magnetic/turbulence mixing.
NASA Astrophysics Data System (ADS)
Chuang, Hsun-Jen; Chamlagain, Bhim; Koehler, Michael; Perera, Meeghage Madusanka; Yan, Jiaqiang; Mandrus, David; Tománek, David; Zhou, Zhixian
2016-03-01
We report a new strategy for fabricating 2D/2D low-resistance ohmic contacts for a variety of transition metal dichalcogenides (TMDs) using van der Waals assembly of substitutionally doped TMDs as drain/source contacts and TMDs with no intentional doping as channel materials. We demonstrate that few-layer WSe2 field-effect transistors (FETs) with 2D/2D contacts exhibit low contact resistances of ~ 0.3 k ohm.um, high on/off ratios up to > 109, and high drive currents exceeding 320 uA um-1. These favorable characteristics are combined with a two-terminal field-effect hole mobility ~ 2x102 cm2 V-1 s-1 at room temperature, which increases to >2x103 cm2 V-1 s-1 at cryogenic temperatures. We observe a similar performance also in MoS2 and MoSe2 FETs with 2D/2D drain and source contacts. The 2D/2D low-resistance ohmic contacts presented here represent a new device paradigm that overcomes a significant bottleneck in the performance of TMDs and a wide variety of other 2D materials as the channel materials in post-silicon electronics.
Chuang, Hsun-Jen; Chamlagain, Bhim; Koehler, Michael; Perera, Meeghage Madusanka; Yan, Jiaqiang; Mandrus, David; Tománek, David; Zhou, Zhixian
2016-03-01
We report a new strategy for fabricating 2D/2D low-resistance ohmic contacts for a variety of transition metal dichalcogenides (TMDs) using van der Waals assembly of substitutionally doped TMDs as drain/source contacts and TMDs with no intentional doping as channel materials. We demonstrate that few-layer WSe2 field-effect transistors (FETs) with 2D/2D contacts exhibit low contact resistances of ∼0.3 kΩ μm, high on/off ratios up to >10(9), and high drive currents exceeding 320 μA μm(-1). These favorable characteristics are combined with a two-terminal field-effect hole mobility μFE ≈ 2 × 10(2) cm(2) V(-1) s(-1) at room temperature, which increases to >2 × 10(3) cm(2) V(-1) s(-1) at cryogenic temperatures. We observe a similar performance also in MoS2 and MoSe2 FETs with 2D/2D drain and source contacts. The 2D/2D low-resistance ohmic contacts presented here represent a new device paradigm that overcomes a significant bottleneck in the performance of TMDs and a wide variety of other 2D materials as the channel materials in postsilicon electronics. PMID:26844954
NASA Astrophysics Data System (ADS)
Choi, M. J.; Park, H. K.; Yun, G. S.; Lee, W.; Luhmann, N. C., Jr.; Lee, K. D.; Ko, W.-H.; Park, Y.-S.; Park, B. H.; In, Y.
2016-06-01
Minor and major disruptions by explosive MHD instabilities were observed with the novel quasi 3D electron cyclotron emission imaging (ECEI) system in the KSTAR plasma. The fine electron temperature (T e) fluctuation images revealed two types of minor disruptions: a small minor disruption is a q∼ 2 localized fast transport event due to a single m/n = 2/1 magnetic island growth, while a large minor disruption is partial collapse of the q≤slant 2 region with two successive fast heat transport events by the correlated m/n = 2/1 and m/n = 1/1 instabilities. The m/n = 2/1 magnetic island growth during the minor disruption is normally limited below the saturation width. However, as the additional interchange-like perturbation grows near the inner separatrix of the 2/1 island, the 2/1 island can expand beyond the limit through coupling with the cold bubble formed by the interchange-like perturbation.
NASA Astrophysics Data System (ADS)
Westerhof, E.; de Blank, H. J.; Pratt, J.
2016-03-01
Two dimensional reduced MHD simulations of neoclassical tearing mode growth and suppression by ECCD are performed. The perturbation of the bootstrap current density and the EC drive current density perturbation are assumed to be functions of the perturbed flux surfaces. In the case of ECCD, this implies that the applied power is flux surface averaged to obtain the EC driven current density distribution. The results are consistent with predictions from the generalized Rutherford equation using common expressions for Δ \\text{bs}\\prime and Δ \\text{ECCD}\\prime . These expressions are commonly perceived to describe only the effect on the tearing mode growth of the helical component of the respective current perturbation acting through the modification of Ohm’s law. Our results show that they describe in addition the effect of the poloidally averaged current density perturbation which acts through modification of the tearing mode stability index. Except for modulated ECCD, the largest contribution to the mode growth comes from this poloidally averaged current density perturbation.
Resistive MHD studies of high-beta Tokamak plasmas
NASA Astrophysics Data System (ADS)
Lynch, V. E.; Hicks, H. R.; Holmes, J. A.; Carreras, B. A.; Garcia, L.
1982-02-01
The magnetohydrodynamic (MHD) activity in high beta Tokamaks such as ISX-B was calculated. These initial value calculations are built on earlier low beta techniques, but the beta effects create several new numerical issues. In addition to time stepping modules, the system of computer codes includes equilibrium solvers (used to provide an initial condition) and output modules, such as a magnetic field line follower and an X-ray diagnostic code. The transition from current driven modes a low beta to predominantly pressure driven modes at high beta is described. The nonlinear studies yield X-ray emissivity plots which are compared with experiment.
Resistive MHD studies of high-. beta. -tokamak plasmas
Lynch, V.E.; Carreras, B.A.; Hicks, H.R.; Holmes, J.A.; Garcia, L.
1981-01-01
Numerical calculations have been performed to study the MHD activity in high-..beta.. tokamaks such as ISX-B. These initial value calculations built on earlier low ..beta.. techniques, but the ..beta.. effects create several new numerical issues. These issues are discussed and resolved. In addition to time-stepping modules, our system of computer codes includes equilibrium solvers (used to provide an initial condition) and output modules, such as a magnetic field line follower and an X-ray diagnostic code. The transition from current driven modes at low ..beta.. to predominantly pressure driven modes at high ..beta.. is described. The nonlinear studies yield X-ray emissivity plots which are compared with experiment.
High-resistance liquid-crystal lens array for rotatable 2D/3D autostereoscopic display.
Chang, Yu-Cheng; Jen, Tai-Hsiang; Ting, Chih-Hung; Huang, Yi-Pai
2014-02-10
A 2D/3D switchable and rotatable autostereoscopic display using a high-resistance liquid-crystal (Hi-R LC) lens array is investigated in this paper. Using high-resistance layers in an LC cell, a gradient electric-field distribution can be formed, which can provide a better lens-like shape of the refractive-index distribution. The advantages of the Hi-R LC lens array are its 2D/3D switchability, rotatability (in the horizontal and vertical directions), low driving voltage (~2 volts) and fast response (~0.6 second). In addition, the Hi-R LC lens array requires only a very simple fabrication process. PMID:24663563
2D resistivity inversion using conjugate gradients for a finite element discretization
NASA Astrophysics Data System (ADS)
Bortolozo, C. A.; Santos, F. M.; Porsani, J. L.
2014-12-01
In this work we present a DC 2D inversion algorithm using conjugate gradients relaxation to solve the maximum likelihood inverse equations. We apply, according to Zhang (1995), the maximum likelihood inverse theory developed by Tarantola and Valette (1982) to our 2D resistivity inversion. This algorithm was chosen to this research because it doesn't need to calculate the field's derivatives. Since conjugate gradient techniques only need the results of the sensitivity matrix Ã or its transpose ÃT multiplying a vector, the actual computation of the sensitivity matrix are not performed, according to the methodology described in Zhang (1995). In Zhang (1995), the terms Ãx and ÃTy, are dependent of the stiffness matrix K and its partial derivative ∂K⁄∂ρ. The inversion methodology described in Zhang (1995) is for the case of 3D electrical resistivity by finite differences discretization. So it was necessary to make a series of adjustments to obtain a satisfactory result for 2D electrical inversion using finite element method. The difference between the modeling of 3D resistivity with finite difference and the 2D finite element method are in the integration variable, used in the 2D case. In the 2D case the electrical potential are initially calculated in the transformed domain, including the stiffness matrix, and only in the end is transformed in Cartesian domain. In the case of 3D, described by Zhang (1995) this is done differently, the calculation is done directly in the Cartesian domain. In the literature was not found any work describing how to deal with this problem. Because the calculations of Ãx and ÃTy must be done without having the real stiffness matrix, the adaptation consist in calculate the stiffness matrix and its partial derivative using a set of integration variables. We transform those matrix in the same form has in the potential case, but with different sets of variables. The results will be presented and are very promising.
Linear Analysis of the m=0 Instability for a Visco-resistive Hall MHD Plasma
Oliver, B.V.; Genoni, T.C.; Mehlhorn, T.A.
2006-01-05
We present a comprehensive analysis of the linear dispersion relation for the axisymmetric (m=0) compressible interchange instability of Bennett equilibria in a visco-resistive, Hall MHD plasma. The full anisotropic stress tensor with Braginskii viscous coefficients is considered. The eigenvalues are obtained numerically. For small axial mode number, Hall currents enhance the growth rates, whereas in the limit of high mode number the growth rates are suppressed, eventually resulting in wave cutoff, even in the ideal limit. For the visco-resistive plasma the unstable spectra are weakly dependent on the off-diagonal elements of the stress tensor.
Interelectrode resistance and performance of small and large scale MHD generators
Doss, E.D.; Picologlou, B.F.
1983-09-01
The effect of reduced interelectrode resistance in MHD generators on the generator power output is investigated. The analytical model used in the investigation allows for the solution for the electric field and current density distributions in the cross plane of the generator. The power output, expressed as a fraction of the power output of a perfectly insulated generator, is found to be a function of the wall temperature, the ratio of boundary layer thickness to channel transverse dimension, and the product of interelectrode resistance and channel cross-sectional area. The interelectrode resistance is assumed to be inversely proportional to the channel transverse dimension and the variation of power output ratio with channel size is calculated. It is found that deterioration of performance of NHD generators, resulting from reduced interelectrode resistance, diminishes with generator size and is negligible for large-scale generators, provided that the interelectrode resistance remains larger than an order of one-tenth ohm.
NASA Astrophysics Data System (ADS)
Den, M.; Horiuchi, R.; Fujita, S.; Tanaka, T.
2011-12-01
Magnetic reconnection is considered to play an important role in space phenomena such as substorm in the Earth's magnetosphere. Tanaka and Fujita reproduced substorm evolution process by numerical simulation with the global MHD code [1]. In the MHD framework, the dissipation model is introduced for modeling of the kinetic effects. They found that the normalized reconnection viscosity, one of the dissipation model employed there, gave a large effect for the dipolarization, central phenomenon in the substorm development process, though that viscosity was assumed to be a constant parameter. It is well known that magnetic reconnection is controlled by microscopic kinetic mechanism. Frozen-in condition is broken due to particle kinetic effects and collisionless reconnection is triggered when current sheet is compressed as thin as ion kinetic scales under the influence of external driving flow [2, 3]. Horiuchi and his collaborators showed that reconnection electric field generated by microscopic physics evolves inside ion meandering scale so as to balance the flux inflow rate at the inflow boundary, which is controlled by macroscopic physics [2]. That is, effective resistivity generated through this process can be expressed by balance equation between micro and macro physics. In this paper, we perform substorm simulation by using the global MHD code developed by Tanaka [3] with this effective resistivity instead of the empirical resistivity model. We obtain the AE indices from simulation data, in which substorm onset can be seen clearly, and investigate the relationship between the substorm development and the effective resistivity model. [1] T. Tanaka, A, Nakamizo, A. Yoshikawa, S. Fujita, H. Shinagawa, H. Shimazu, T. Kikuchi, and K. K. Hashimoto, J. Geophys. Res. 115 (2010) A05220,doi:10.1029/2009JA014676. [2] W. Pei, R. Horiuchi, and T. Sato, Physics of Plasmas,Vol. 8 (2001), pp. 3251-3257. [3] A. Ishizawa, and R. Horiuchi, Phys. Rev. Lett., Vol. 95, 045003 (2005). [4
Minimally implicit Runge-Kutta methods for Resistive Relativistic MHD
NASA Astrophysics Data System (ADS)
Aloy, Miguel-Á.; Cordero-Carrión, Isabel
2016-05-01
The Relativistic Resistive Magnetohydrodynamic (RRMHD) equations are a hyperbolic system of partial differential equations used to describe the dynamics of relativistic magnetized fluids with a finite conductivity. Close to the ideal magnetohydrodynamic regime, the source term proportional to the conductivity becomes potentially stiff and cannot be handled with standard explicit time integration methods. We propose a new class of methods to deal with the stiffness fo the system, which we name Minimally Implicit Runge-Kutta methods. These methods avoid the development of numerical instabilities without increasing the computational costs in comparison with explicit methods, need no iterative extra loop in order to recover the primitive (physical) variables, the analytical inversion of the implicit operator is trivial and the several stages can actually be viewed as stages of explicit Runge-Kutta methods with an effective time-step. We test these methods with two different one-dimensional test beds in varied conductivity regimes, and show that our second-order schemes satisfy the theoretical expectations.
Impact of resistive MHD plasma response on perturbation field sidebands
NASA Astrophysics Data System (ADS)
Orlov, D. M.; Evans, T. E.; Moyer, R. A.; Lyons, B. C.; Ferraro, N. M.; Park, G.-Y.
2016-07-01
Single fluid linear simulations of a KSTAR RMP ELM suppressed discharge with the M3D-C1 resistive magnetohydrodynamic code have been performed for the first time. The simulations show that the application of the n = 1 perturbation using the KSTAR in-vessel control coils (IVCC), which apply modest levels of n = 3 sidebands (~20% of the n = 1), leads to levels of n = 3 sideband that are comparable to the n = 1 when plasma response is included. This is due to the reduced level of screening of the rational-surface-resonant n = 3 component relative to the rational-surface-resonant n = 1 component. The n = 3 sidebands could play a similar role in ELM suppression on KSTAR as the toroidal sidebands (n = 1, 2, 4) in DIII-D n = 3 ELM suppression with missing I-coil segments (Paz Soldan et al 2014 Nucl. Fusion 54 073013). This result may help to explain the uniqueness of ELM suppression with n = 1 perturbations in KSTAR since the effective perturbation is a mixed n = 1/n = 3 perturbation similar to n = 3 ELM suppression in DIII-D.
NASA Astrophysics Data System (ADS)
Ali, Nisa'; Saad, Rosli; Muztaza, Nordiana M.; Ismail, Noer E. H.
2013-05-01
The purpose of this study was to locate the geological contact using 2D resistivity method for Light Rail Transit (LRT) track alignment. The resistivity method was conducted on eight survey lines with the length of line 1 was 600m. The length of line 2, 3, 4, 5, 6, and 7 were 200m each while line 8 is 115m. All the survey used minimum electrode spacing of 5m and using Pole-dipole array with minimum current is 2mA and maximum was 20mA. The result obtained from the pseudosection showed that the area generally divided into three main zones, fill materials/residual soil with a resistivity value of <500 Ωm, saturated zone with a resistivity value of 30-100 Ωm and bedrock with a resistivity value of >2000 Ωm. Three fractured zones were detected along line L1 and a lot of boulders were detected at L1, L3, L4, L5 and L6. The geological contact was between the residual soil and granite bedrock.
A Novel Crosstalk Suppression Method of the 2-D Networked Resistive Sensor Array
Wu, Jianfeng; Wang, Lei; Li, Jianqing; Song, Aiguo
2014-01-01
The 2-D resistive sensor array in the row–column fashion suffered from the crosstalk problem for parasitic parallel paths. Firstly, we proposed an Improved Isolated Drive Feedback Circuit with Compensation (IIDFCC) based on the voltage feedback method to suppress the crosstalk. In this method, a compensated resistor was specially used to reduce the crosstalk caused by the column multiplexer resistors and the adjacent row elements. Then, a mathematical equivalent resistance expression of the element being tested (EBT) of this circuit was analytically derived and verified by the circuit simulations. The simulation results show that the measurement method can greatly reduce the influence on the EBT caused by parasitic parallel paths for the multiplexers' channel resistor and the adjacent elements. PMID:25046011
NASA Astrophysics Data System (ADS)
Yan, Bo; Li, Yuguo; Liu, Ying
2016-07-01
In this paper, we present an adaptive finite element (FE) algorithm for direct current (DC) resistivity modeling in 2-D generally anisotropic conductivity structures. Our algorithm is implemented on an unstructured triangular mesh that readily accommodates complex structures such as topography and dipping layers and so on. We implement a self-adaptive, goal-oriented grid refinement algorithm in which the finite element analysis is performed on a sequence of refined grids. The grid refinement process is guided by an a posteriori error estimator. The problem is formulated in terms of total potentials where mixed boundary conditions are incorporated. This type of boundary condition is superior to the Dirichlet type of conditions and improves numerical accuracy considerably according to model calculations. We have verified the adaptive finite element algorithm using a two-layered earth with azimuthal anisotropy. The FE algorithm with incorporation of mixed boundary conditions achieves high accuracy. The relative error between the numerical and analytical solutions is less than 1% except in the vicinity of the current source location, where the relative error is up to 2.4%. A 2-D anisotropic model is used to demonstrate the effects of anisotropy upon the apparent resistivity in DC soundings.
MHD computation of feedback of resistive-shell instabilities in the reversed field pinch
Zita, E.J.; Prager, S.C. . Plasma Physics Research); Ho, Y.L.; Schnack, D.D )
1992-05-01
MHD computation demonstrates that feedback can sustain reversal and reduce loop voltage in resistive-shell reversed field pinch (RFP) plasmas. Edge feedback on {approximately}2R/a tearing modes resonant near axis is found to restore plasma parameters to nearly their levels with a close-fitting conducting shell. When original dynamo modes are stabilized, neighboring tearing modes grow to maintain the RFP dynamo more efficiently. This suggests that experimentally observed limits on RFP pulselengths to the order of the shell time can be overcome by applying feedback to a few helical modes.
Shadid, J. N.; Pawlowski, R. P.; Cyr, E. C.; Tuminaro, R. S.; Chacon, L.; Weber, P. D.
2016-02-10
Here, we discuss that the computational solution of the governing balance equations for mass, momentum, heat transfer and magnetic induction for resistive magnetohydrodynamics (MHD) systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena, as well as the significant range of time- and length-scales that the interactions of these physical mechanisms produce. This paper explores the development of a scalable, fully-implicit stabilized unstructured finite element (FE) capability for 3D incompressible resistive MHD. The discussion considers the development of a stabilized FE formulation in context of the variational multiscale (VMS) method,more » and describes the scalable implicit time integration and direct-to-steady-state solution capability. The nonlinear solver strategy employs Newton–Krylov methods, which are preconditioned using fully-coupled algebraic multilevel preconditioners. These preconditioners are shown to enable a robust, scalable and efficient solution approach for the large-scale sparse linear systems generated by the Newton linearization. Verification results demonstrate the expected order-of-accuracy for the stabilized FE discretization. The approach is tested on a variety of prototype problems, that include MHD duct flows, an unstable hydromagnetic Kelvin–Helmholtz shear layer, and a 3D island coalescence problem used to model magnetic reconnection. Initial results that explore the scaling of the solution methods are also presented on up to 128K processors for problems with up to 1.8B unknowns on a CrayXK7.« less
NASA Astrophysics Data System (ADS)
Kiflu, H. G.; Kruse, S. E.; Harro, D.; Loke, M. H.; Wilkinson, P. B.
2013-12-01
Electrical resistivity tomography is commonly used to identify geologic features associated with sinkhole formation. In covered karst terrain, however, it can be difficult to resolve the depth to top of limestone with this method. This is due to the fact that array lengths, and hence depth of resolution, are often limited by residential or commercial lot dimensions in urban environments. Furthermore, the sediments mantling the limestone are often clay-rich and highly conductive. The resistivity method has limited sensitivity to resistive zones beneath conductive zones. This sensitivity can be improved significantly with electrodes implanted at depth in the cover sediments near the top of limestone. An array of deep electrodes is installed with direct push technology in the karst cover. When combined with a surface array in which each surface electrode is underlain by a deep electrode, the array geometry is similar to a borehole array turned on its side. This method, called the Multi-Electrode Resistivity Implant Technique (MERIT), offers the promise of significantly improved resolution of epikarst and cover collapse development zones in the overlying sediment, the limestone or at the sediment-bedrock interface in heterogeneous karst environments. With a non-traditional array design, the question of optimal array geometries arises. Optimizing array geometries is complicated by the fact that many plausible 4-electrode readings will produce negative apparent resistivity values, even in homogeneous terrain. Negative apparent resistivities cannot be used in inversions based on the logarithm of the apparent resistivity. New algorithms for seeking optimal array geometries have been developed by modifying the 'Compare R' method of Wilkinson and Loke. The optimized arrays show significantly improved resolution over basic arrays adapted from traditional 2D surface geometries. Several MERIT case study surveys have been conducted in covered karst in west-central Florida, with
Electrical resistivity tomography applied to a complex lava dome: 2D and 3D models comparison
NASA Astrophysics Data System (ADS)
Portal, Angélie; Fargier, Yannick; Lénat, Jean-François; Labazuy, Philippe
2015-04-01
The study of volcanic domes growth (e.g. St. Helens, Unzen, Montserrat) shows that it is often characterized by a succession of extrusion phases, dome explosions and collapse events. Lava dome eruptive activity may last from days to decades. Therefore, their internal structure, at the end of the eruption, is complex and includes massive extrusions and lava lobes, talus and pyroclastic deposits as well as hydrothermal alteration. The electrical resistivity tomography (ERT) method, initially developed for environmental and engineering exploration, is now commonly used for volcano structure imaging. Because a large range of resistivity values is often observed in volcanic environments, the method is well suited to study the internal structure of volcanic edifices. We performed an ERT survey on an 11ka years old trachytic lava dome, the Puy de Dôme volcano (French Massif Central). The analysis of a recent high resolution DEM (LiDAR 0.5 m), as well as other geophysical data, strongly suggest that the Puy de Dôme is a composite dome. 11 ERT profiles have been carried out, both at the scale of the entire dome (base diameter of ~2 km and height of 400 m) on the one hand, and at a smaller scale on the summit part on the other hand. Each profile is composed of 64 electrodes. Three different electrode spacing have been used depending on the study area (35 m for the entire dome, 10 m and 5 m for its summit part). Some profiles were performed with half-length roll-along acquisitions, in order to keep a good trade-off between depth of investigation and resolution. Both Wenner-alpha and Wenner-Schlumberger protocols were used. 2-D models of the electrical resistivity distribution were computed using RES2DINV software. In order to constrain inversion models interpretation, the depth of investigation (DOI) method was applied to those results. It aims to compute a sensitivity index on inversion results, illustrating how the data influence the model and constraining models
Broken Ergodicity in MHD Turbulence in a Spherical Domain
NASA Technical Reports Server (NTRS)
Shebalin, John V.; wang, Yifan
2011-01-01
Broken ergodicity (BE) occurs in Fourier method numerical simulations of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence. Although naive statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that low-wave-number coefficients have non-zero mean values that can be very large compared to the associated standard deviation. In other words, large-scale coherent structure (i.e., broken ergodicity) in homogeneous MHD turbulence can spontaneously grow out of random initial conditions. Eigenanalysis of the modal covariance matrices in the probability density functions of ideal statistical theory leads to a theoretical explanation of observed BE in homogeneous MHD turbulence. Since dissipation is minimal at the largest scales, BE is also relevant for resistive magnetofluids, as evidenced in numerical simulations. Here, we move beyond model magnetofluids confined by periodic boxes to examine BE in rotating magnetofluids in spherical domains using spherical harmonic expansions along with suitable boundary conditions. We present theoretical results for 3-D and 2-D spherical models and also present computational results from dynamical simulations of 2-D MHD turbulence on a rotating spherical surface. MHD turbulence on a 2-D sphere is affected by Coriolus forces, while MHD turbulence on a 2-D plane is not, so that 2-D spherical models are a useful (and simpler) intermediate stage on the path to understanding the much more complex 3-D spherical case.
Plasmoid dynamics in 3D resistive MHD simulations of magnetic reconnection
NASA Astrophysics Data System (ADS)
Samtaney, R.; Loureiro, N. F.; Uzdensky, D. A.; Schekochihin, A. A.
2012-04-01
Magnetic reconnection is a well known plasma process believed to lie at the heart of a variety of phenomena such as sub-storms in the Earth's magnetosphere, solar/stellar and accretion-disk flares, sawteeth activity in fusion devices, etc. During reconnection, the global magnetic field topology changes rapidly, leading to the violent release of magnetic energy. Over the past few years, the basic understanding of this fundamental process has undergone profound changes. The validity of the most basic, and widely accepted, reconnection paradigm - the famous Sweet-Parker (SP) model, which predicts that, in MHD, reconnection is extremely slow, its rate scaling as S-1/2, where S is the Lundquist number of the system - has been called into question as it was analytically demonstrated that, for S ≫ 1, SP-like current sheets are violently unstable to the formation of a large number of secondary islands, or plasmoids. Subsequent numerical simulations in 2D have confirmed the validity of the linear theory, and shown that plasmoids quickly grow to become wider than the thickness of the original SP current sheet, thus effectively changing the underlying reconnection geometry. Ensuing numerical work has revealed that the process of plasmoid formation, coalescence and ejection from the sheet drastically modifies the steady state picture assumed by Sweet and Parker, and leads to the unexpected result that MHD reconnection is independent of S. In this talk, we review these recent developments and present results from three-dimensional simulations of high-Lundquist number reconnection in the presence of a guide field. A parametric study varying the strength of the guide field is presented. Plasmoid flux and width distribution functions are quantified and compared with corresponding two dimensional simulations.
Altered NK Cell Development and Enhanced NK Cell-Mediated Resistance to MCMV in NKG2D-Deficient Mice
Zafirova, Biljana; Mandarić, Sanja; Antulov, Ronald; Krmpotić, Astrid; Jonsson, Helena; Yokoyama, Wayne M.; Jonjić, Stipan; Polić, Bojan
2009-01-01
Summary NKG2D is a potent activating receptor on NK cells which acts as a molecular sensor for stressed cells expressing NKG2D ligands such as infected or tumor transformed cells. Although NKG2D is expressed on NK cell precursors, its role in NK cell development is still not known. We have generated NKG2D-deficient mice by targeting the Klrk1 locus. Here we provide evidence for an important regulatory role of NKG2D in the development of NK cells. The absence of NKG2D causes faster division of NK cells, perturbation in size of some NK cell subpopulations and their augmented sensitivity to apoptosis. As expected, NKG2D−/− NK cells are less responsive to tumor targets expressing NKG2D ligands. NKG2D−/− mice, however, show an enhanced NK cell-mediated resistance to MCMV infection as a consequence of NK cell dysregulation. Altogether, these findings provide evidence for yet unknown regulatory function of NKG2D in NK cell physiology. PMID:19631564
Residual resistance of 2D and 3D structures and Joule heat release.
Gurevich, V L; Kozub, V I
2011-06-22
We consider a residual resistance and Joule heat release in 2D nanostructures as well as in ordinary 3D conductors. We assume that elastic scattering of conduction electrons by lattice defects is predominant. Within a rather intricate situation in such systems we discuss in detail two cases. (1) The elastic scattering alone (i.e. without regard of inelastic mechanisms of scattering) leads to a transition of the mechanical energy (stored by the electrons under the action of an electric field) into heat in a traditional way. This process can be described by the Boltzmann equation where it is possible to do the configuration averaging over defect positions in the electron-impurity collision term. The corresponding conditions are usually met in metals. (2) The elastic scattering can be considered with the help of the standard electron-impurity collision integral only in combination with some additional averaging procedure (possibly including inelastic scattering or some mechanisms of electron wavefunction phase destruction). This situation is typical for degenerate semiconductors with a high concentration of dopants and conduction electrons. Quite often, heat release can be observed via transfer of heat to the lattice, i.e. via inelastic processes of electron-phonon collisions and can take place at distances much larger than the size of the device. However, a direct heating of the electron system can be registered too by, for instance, local measurements of the current noise or direct measurement of an electron distribution function. PMID:21628783
2D-resistivity surveys of deteriorating historic stonework in Oxford, UK
NASA Astrophysics Data System (ADS)
Sass, O.; Viles, H. A.
2009-04-01
Historic stonework deteriorates in often very complex ways and despite many years of research on the topic, we are still a long way from being able to predict its occurrence and severity. As most deterioration processes involve water, techniques which provide a better picture of the moisture contents and fluctuations within stonework are very valuable in attempts to improve understanding. 2D resistivity methods can provide useful information about moisture distributions within porous historic stonework. We report on a series of experiments on historic walls within the centre of Oxford, UK, which illustrate varying degrees of deterioration including catastrophic decay. Using medical electrodes we have been able to carry out non-invasive and non-destructive 2D resistivity surveys to study the distribution and amount of water stored in deteriorating limestone walls. Fifteen vertical profiles, each 2-2.5 m in length, have been monitored at five sites. Furthermore, simulated driving rain experiments have been carried out at two sites. The data indicate the diversity and complexity of moisture distributions within these walls. Replacement stone patches show consistently higher moisture conditions than the surrounding stone. Some profiles show wetter sections towards the base of the wall, usually where a plinth is absent. Conversely, hard stone plinths obviously reduce capillary rise from ground water. However, at several sites we noticed a wetter zone immediately above the top of the plinth which often correlates with the occurrence of catastrophic decay - indicating that the plinth may encourage concentration of decay. Most profiles indicate the presence of wetter patches 5-10cm behind the wall face under blackened crusts. Such patches of heightened absolute moisture contents could play a very important role in encouraging catastrophic decay. Severely decayed sections of profiles often exhibit wetter near-surface conditions than surrounding stonework, whilst areas with
Computational modeling of neoclassical and resistive MHD tearing modes in tokamaks
Gianakon, T.A.
1996-02-01
Numerical studies of the linear and nonlinear evolution of magnetic tearing type modes in three-dimensional toroidal geometry are presented. In addition to traditional resistive MHD effects, where the parameter {Delta}{prime} determines the stability properties, neoclassical effects have been included for the first time in such models. The inclusion of neoclassical physics introduces and additional free-energy source for the nonlinear formation of magnetic islands through the effects of a bootstrap current in Ohm`s law. The neoclassical tearing mode is demonstrated to be destabilized in plasmas which are otherwise {Delta}{prime} stable, albeit once an island width threshold is exceeded. The simulations are based on a set of neoclassical reduced magnetohydrodynamic (MHD) equations in three-dimensional toroidal geometry derived from the two-fluid equations in the limit of small inverse aspect ratio {epsilon} and low plasma pressure {beta} with neoclassical closures for the viscous force {del} {center_dot} {leftrightarrow}{pi}. The poloidal magnetic flux {psi}, the toroidal vorticity {omega}{sup {zeta}}, and the plasma pressure p are time advanced using the parallel projection of Ohm`s law, the toroidal projection of the curl of the momentum equation, and a pressure evolution equation with anisotropic pressure transport parallel to and across magnetic field lines. The equations are implemented in an initial value code which Fourier decomposes equilibrium and perturbation quantities in the poloidal and toroidal directions, and finite differences them radially based on a equilibrium straight magnetic field line representation. An implicit algorithm is used to advance the linear terms; the nonlinear terms are advanced explicitly. The simulations are benchmarked linearly and nonlinearly against single and multiple helicity {Delta}{prime} tearing modes in toroidal geometry in the absence of neo-classical effects.
2-D resistivity distribution in the seismo-genic zone of the 2011 Tohoku earthquake (M9.0)
NASA Astrophysics Data System (ADS)
Ichihara, H.; Kasaya, T.; Baba, K.
2012-12-01
We discussed 2-D resistivity distribution around the rupture zone of the 2011 Tohoku earthquake (M9.0) in order to clarify structural properties that enabled such a large fault slip. Electro-magnetic time series were collected using ocean bottom electro-magnetometers (OBEMs) along a survey line in the latitude 38N between 2009 and 2011. The MT responses show the out of the quadrant phases in TE mode. The out of quadrant phases were explained by coastal effect inducing strong distortion especially in the magnetic field [Key and Constable, 2011] and indicate strong resistivity contrast between seawater and crustal area. Then we carried out 2-D resistivity inversion after the correction of 3-D bathymetric effect. The inverted resistivity model entirely consists of conductive surface layer (< 3 ohm-m) and subsequence resistive area (> 300 ohm-m). While the main rupture zone shows high resistivity (> 3000 ohm-m), resistivity near the trench and the deeper area along the plate boundary is relatively low. The former conductive zone probably reflects thick sediments containing large amount of fluid. The lower conductive zone possibly indicates connected pore fluids due to dehydration. Thus the main rupture zone seems to be located on a resistive area sandwiched by conductive zones. It implies that heterogeneous distribution of fluid controlled fault activity of the 2011 earthquake.
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)
Pérez-Corona, M.; García, J. A.; Taller, G.; Polgár, D.; Bustos, E.; Plank, Z.
2016-02-01
The purpose of geophysical electrical surveys is to determine the subsurface resistivity distribution by making measurements on the ground surface. From these measurements, the true resistivity of the subsurface can be estimated. The ground resistivity is related to various geological parameters, such as the mineral and fluid content, porosity and degree of water saturation in the rock. Electrical resistivity surveys have been used for many decades in hydrogeological, mining and geotechnical investigations. More recently, they have been used for environmental surveys. To obtain a more accurate subsurface model than is possible with a simple 1-D model, a more complex model must be used. In a 2-D model, the resistivity values are allowed to vary in one horizontal direction (usually referred to as the x direction) but are assumed to be constant in the other horizontal (the y) direction. A more realistic model would be a fully 3-D model where the resistivity values are allowed to change in all three directions. In this research, a simulation of the cone penetration test and 2D imaging resistivity are used as tools to simulate the distribution of hydrocarbons in soil.
NASA Astrophysics Data System (ADS)
Kumar, Dewashish; Rao, V. Ananda; Sarma, V. S.
2014-04-01
Electrical resistivity method is a versatile and economical technique for groundwater prospecting in different geological settings due to wide spectrum of resistivity compared to other geophysical parameters. Exploration and exploitation of groundwater, a vital and precious resource, is a challenging task in hard rock, which exhibits inherent heterogeneity. In the present study, two-dimensional Electrical Resistivity Tomography (2D-ERT) technique using two different arrays, viz., pole-dipole and pole-pole, were deployed to look into high signal strength data in a tectonically disturbed hard rock ridge region for groundwater. Four selected sites were investigated. 2D subsurface resistivity tomography data were collected using Syscal Pro Switch-10 channel system and covered a 2 km long profile in a tough terrain. The hydrogeological interpretation based on resistivity models reveal the water horizons trap within the clayey sand and weathered/fractured quartzite formations. Aquifer resistivity lies between ˜3-35 and 100-200 Ωm. The results of the resistivity models decipher potential aquifer lying between 40 and 88 m depth, nevertheless, it corroborates with the static water level measurements in the area of study. The advantage of using pole-pole in conjunction with the pole-dipole array is well appreciated and proved worth which gives clear insight of the aquifer extent, variability and their dimension from shallow to deeper strata from the hydrogeological perspective in the present geological context.
2D Electrical Resistivity Tomography surveys optimisation of the solutes transports in porous media.
NASA Astrophysics Data System (ADS)
Lekmine, G.; Pessel, M.; Auradou, H.
2009-04-01
Electrical resistivity tomography applied in borehole or cross-borehole is a method often used to follow the invasion process of pollutant [Daily, 1991]. The aim of this work is to test experimentally the electrode arrays and inversion process used to obtain a spatial representation of tracer propagation in porous media. Experiments were conducted in a plexiglas container with glass beads of 166 microns in diameter. The height of the container is 275 mm, its width 85 mm and its thickness 10 mm. 21 electrodes, equally spaced, are placed along each of the lateral sides of the porous medium : these electrodes are used to perform the electrical measurements. The porous medium is lightened from behind and a video camera records the propagation of the fluids. The fluid containing the tracer (i.e the pollutant) is a water solution containing a small amount of dye together with NaCl (0.5g/l up to 2.0g/l). The medium is first saturated by a water solution containing a slight concentration of NaCl so that its density is smaller than the injected fluid. An upward flow is first established, then the denser fluid is injected at the bottom and over the full width of the medium. In this way, the flow is stabilized by gravity avoiding the development of unstable fingers. Still, the fluids are miscible and a mixing front develops during the flow: in the present study, both the determination by optical and electrical imaging of the mean position of the front and its width are of interest. The comparison of the two techniques allows to study the ability of the inversion process to quantify the solute transport. The electrical measurements are acquired by a standard multi electrode system (IRIS Instruments) and the data are inverted with the Res2Dinv software which models the 2D distribution of conductivity contrasts. The obtained bulk conductivity can be related through Archie's law to fluid conductivity by the porosity and the cementation factor which have been experimentally
NASA Astrophysics Data System (ADS)
Azwin, I. N.; Saad, Rosli; Saidin, Mokhtar; Nordiana, M. M.; Anderson Bery, Andy; Hidayah, I. N. E.
2015-01-01
Interest in studying impact crater on earth has increased tremendously due to its importance in geologic events, earth inhabitant history as well as economic value. The existences of few shock metamorphism and crater morphology evidences are discovered in Bukit Bunuh, Malaysia thus detailed studies are performed using geophysical and geotechnical methods to verify the type of the crater and characteristics accordingly. This paper presents the combined analysis of 2-D electrical resistivity, seismic refraction, geotechnical SPT N value, moisture content and RQD within the study area. Three stages of data acquisition are made starting with regional study followed by detailed study on West side and East side. Bulk resistivity and p-wave seismic velocity were digitized from 2-D resistivity and seismic sections at specific distance and depth for corresponding boreholes and samples taken. Generally, Bukit Bunuh shows the complex crater characteristics. Standard table of bulk resistivity and p-wave seismic velocity against SPT N value, moisture content and RQD are produce according to geological classifications of impact crater; inside crater, rim/slumped terrace and outside crater.
Baker, Theresa; Nerle, Sujata; Pritchard, Justin; Zhao, Boyang; Rivera, Victor M.
2015-01-01
Although targeted therapies have revolutionized cancer treatment, overcoming acquired resistance remains a major clinical challenge. EZH2 inhibitors (EZH2i), EPZ-6438 and GSK126, are currently in the early stages of clinical evaluation and the first encouraging signs of efficacy have recently emerged in the clinic. To anticipate mechanisms of resistance to EZH2i, we used a forward genetic platform combining a mutagenesis screen with next generation sequencing technology and identified a hotspot of secondary mutations in the EZH2 D1 domain (Y111 and I109). Y111D mutation within the WT or A677G EZH2 allele conferred robust resistance to both EPZ-6438 and GSK126, but it only drove a partial resistance within the Y641F allele. EZH2 mutants required histone methyltransferase (HMT) catalytic activity and the polycomb repressive complex 2 (PRC2) components, SUZ12 and EED, to drive drug resistance. Furthermore, D1 domain mutations not only blocked the ability of EZH2i to bind to WT and A677G mutant, but also abrogated drug binding to the Y641F mutant. These data provide the first cellular validation of the mechanistic model underpinning the oncogenic function of WT and mutant EZH2. Importantly, our findings suggest that acquired-resistance to EZH2i may arise in WT and mutant EZH2 patients through a single mutation that remains targetable by second generation EZH2i. PMID:26360609
IP4DI: A software for time-lapse 2D/3D DC-resistivity and induced polarization tomography
NASA Astrophysics Data System (ADS)
Karaoulis, M.; Revil, A.; Tsourlos, P.; Werkema, D. D.; Minsley, B. J.
2013-04-01
We propose a 2D/3D forward modelling and inversion package to invert direct current (DC)-resistivity, time-domain induced polarization (TDIP), and frequency-domain induced polarization (FDIP) data. Each cell used for the discretization of the 2D/3D problems is characterized by a DC-resistivity value and a chargeability or complex conductivity for TDIP/FDIP problems, respectively. The governing elliptic partial differential equations are solved with the finite element method, which can be applied for both real and complex numbers. The inversion can be performed either for a single snapshot of data or for a sequence of snapshots in order to monitor a dynamic process such as a salt tracer test. For the time-lapse inversion, we have developed an active time constrained (ATC) approach that is very efficient in filtering out noise in the data that is not correlated over time. The forward algorithm is benchmarked with simple analytical solutions. The inversion package IP4DI is benchmarked with three tests, two including simple geometries. The last one corresponds to a time-lapse resistivity problem for cross-well tomography during enhanced oil recovery. The algorithms are based on MATLAB® code package and a graphical user interface (GUI).
Body edge delineation in 2D DC resistivity imaging using differential method
NASA Astrophysics Data System (ADS)
Susanto, Kusnahadi; Fitrah Bahari, Mohammad
2016-01-01
DC resistivity is widely used to identify the kind of rock and the lithology contact. However, the image resulting from resistivity processing is shown in a contour image. There is be a problem to interpret where the edge of body location is. This study uses differential method to delineate the edge of body in DC resistivity contour. This method was applied to the boundary between gravel and underlying clay layer. The first and the second order differential method is applied to the delineation of lithology contact. The profiling curve has to be sliced and extracted from the resistivity contour before the differential method can be used. The spectral analysis shows the frequency and wavenumber of the profiling curve used to make gridding. The slicing process was conducted horizontally and vertically in order to get the mesh size which will be used in the differential method. The second order differential, the Laplace operator, is able to show the edge of body more clearly than the first order differential and shows the contact between gravel and clay.
Drug Abuse Resistance Education Act of 1990. House of Representatives, 101st Congress, 2d Session.
ERIC Educational Resources Information Center
Congress of the U.S., Washington, DC. House Committee on Education and Labor.
This document is a congressional report on the Drug Abuse Resistance Education Act of 1990, House of Representatives (H.R.) 5064 which requires the Secretary of Education to reserve $15 million from the appropriations made for the Drug Free Schools and Communities Act for grants to consortia of local educational agencies and other agencies to…
Resistively detected high-order magnetoplasmons in a high-quality 2D electron gas
NASA Astrophysics Data System (ADS)
Zudov, M. A.; Shi, Q.; Pfeiffer, L. N.; West, K. W.; Watson, J. D.; Manfra, M. J.
We report on high-order magnetoplasmon resonances detected in photoresistance in high-mobility GaAs quantum wells. These resonances manifest themselves as a series of resistance extrema in the regime of Shubnikov-de Haas oscillations. Extending to orders above 20, the extrema exhibit alternating strength, being less (more) pronounced at even (odd) order magnetoplasmon modes. The lower magnetoplasmon modes reveal the importance of retardation effects.
Tuning and simulating a 193-nm resist for 2D applications
NASA Astrophysics Data System (ADS)
Howard, William B.; Wiaux, Vincent; Ercken, Monique; Bui, Bang; Byers, Jeff D.; Pochkowski, Mike
2002-07-01
For some applications, the usefulness of lithography simulation results depends strongly on the matching between experimental conditions and the simulation input parameters. If this matching is optimized and other sources of error are minimized, then the lithography model can be used to explain printed wafer experimental results. Further, simulation can be useful in predicting the results or in choosing the correct set of experiments. In this paper, PROLITH and ProDATA AutoTune were used to systematically vary simulation input parameters to match measured results on printed wafers used in a 193 nm process. The validity of the simulation parameters was then checked using 3D simulation compared to 2D top-down SEM images. The quality of matching was evaluated using the 1D metrics of average gate CD and Line End Shortening (LES). To ensure the most accurate simulation, a new approach was taken to create a compound mask from GDSII contextual information surrounding an accurate SEM image of the reticle region of interest. Corrections were made to account for all metrology offsets.
NASA Astrophysics Data System (ADS)
Kamaruddin, Nur Aminuda; Saad, Rosli; Nordiana, M. M.; Azwin, I. N.
2015-04-01
The Great Sumatra Fault system was split into two sub-parallel lines or segments at the Northern Sumatra. This event is one of the impacts of powerful earthquakes that hit Sumatra Island especially one that occurred in 2004. These two sub-parallel segments known as Aceh and Seulimeum fault. The study is focused on the Seulimeum fault and two geophysical methods chosen aimed to compare and verified the result obtained respectively. 2-D resistivity method is a common geophysical method used in determination of near surface structures such as faults, cavities, voids and sinkholes. Meanwhile, the magnetic method often chosen to delineate subsurface structures, determine depth of magnetic source bodies and possibly sediment thickness. Three survey lines of resistivity method and randomly magnetic stations were carried out covering Krueng district. The resistivity data processed using Res2Dinv and result presented using Surfer software. The fault identified by the contrast of low and high resistivity value. Meanwhile, the magnetic data were presented in magnetic residual contour map and the extended fault system is suspected represent by the contrast value of the magnetic anomalies. Within suspected fault zone, the results of resistivity are tally with magnetic result.
Magnetic turbulence and resistive MHD instabilities in a 0. 6 < q < 3 poloidal divertor tokamak
Agim, Y.Z.; Callen, J.D.; Chang, Z.; Dexter, R.N.; Goetz, J.A.; Graessle, D.E.; Haines, E.; Kortbawi, D.; LaPointe, M.A.; Moyer, R.A.
1988-09-01
Detailed statistical properties of internal magnetic turbulence, and internal disruptions in magnetically- and materially-limited discharges, are studied in the Tokapole II poloidal divertor tokamak over the safety factor range 0.6 < q{sub a} < 3. A nonlinear MHD code treats tearing modes in the divertor geometry. 9 refs., 2 figs.
A 2D vector map watermarking algorithm resistant to simplication attack
NASA Astrophysics Data System (ADS)
Wang, Chuanjian; Liang, Bin; Zhao, Qingzhan; Qiu, Zuqi; Peng, Yuwei; Yu, Liang
2009-12-01
Vector maps are valuable asset of data producers. How to protect copyright of vector maps effectively using digital watermarking is a hot research issue. In this paper, we propose a new robust and blind watermarking algorithm resilient to simplification attack. We proof that spatial topological relation between map objects bears an important property of approximate simplification invariance. We choose spatial topological relations as watermark feature domain and embed watermarks by slightly modifying spatial topological relation between map objects. Experiment shows that our algorithm has good performance to resist simplification attack and tradeoff of the robustness and data fidelity is acquired.
Optimized arrays for 2-D resistivity survey lines with a large number of electrodes
NASA Astrophysics Data System (ADS)
Loke, M. H.; Wilkinson, P. B.; Chambers, J. E.; Uhlemann, S. S.; Sorensen, J. P. R.
2015-01-01
Previous studies show that optimized arrays generated using the 'Compare R' method have significantly better resolution than conventional arrays. This method determines the optimum set of arrays by selecting those that give the maximum model resolution. The number of possible arrays (the comprehensive data set) increases with the fourth power of the number of electrodes. The optimization method faces practical limitations for 2-D survey lines with more than 60 electrodes where the number of possible arrays exceeds a million. Several techniques are proposed to reduce the calculation time for such survey lines. A single-precision version of the 'Compare R' algorithm using a new ranking function reduces the calculation time by two to eight times while providing results similar to the double-precision version. Recent improvements in computer GPU technology can reduce the calculation time by about seven times. The calculation time is reduced by half by using the fact that arrays that are symmetrical about the center of the line produce identical changes in the model resolution values. It is further reduced by more than thirty times by calculating the Sherman-Morrison update for all the possible two-electrode combinations, which are then used to calculate the model resolution values for the four-electrode arrays. The calculation time is reduced by more then ten times by using a subset of the comprehensive data set consisting of only symmetrical arrays. Tests with a synthetic model and field data set show that optimized arrays derived from this subset produce inversion models with differences of less than 10% from those derived using the full comprehensive data set. The optimized data sets produced models that are more accurate than the Wenner-Schlumberger array data sets in all the tests.
2D Time-lapse Resistivity Monitoring of an Organic Produced Gas Plume in a Landfill using ERT.
NASA Astrophysics Data System (ADS)
Amaral, N. D.; Mendonça, C. A.; Doherty, R.
2014-12-01
This project has the objective to study a landfill located on the margins of Tietê River, in São Paulo, Brazil, using the electroresistivity tomography method (ERT). Due to huge organic matter concentrations in the São Paulo Basin quaternary sediments, there is subsurface depth related biogas accumulation (CH4 and CO2), induced by anaerobic degradation of the organic matter. 2D resistivity sections were obtained from a test area since March 2012, a total of 7 databases, being the last one dated from October 2013. The studied line has the length of 56m, the electrode interval is of 2m. In addition, there are two boreholes along the line (one with 3 electrodes and the other one with 2) in order to improve data quality and precision. The boreholes also have a multi-level sampling system that indicates the fluid (gas or water) presence in relation to depth. With our results it was possible to map the gas plume position and its area of extension in the sections as it is a positive resistivity anomaly, with the gas level having approximately 5m depth. With the time-lapse analysis (Matlab script) between the obtained 2D resistivity sections from the site, it was possible to map how the biogas volume and position change in the landfill in relation to time. Our preliminary results show a preferential gas pathway through the subsurface studied area. A consistent relation between the gas depth and obtained microbiological data from archea and bacteria population was also observed.
Comparison of 1D, 2D and 2.5D Constrained Inversion of Electrical Resistivity Data
NASA Astrophysics Data System (ADS)
Catt, L. M.; West, J.; Clark, R. A.
2007-05-01
Clay-rich till plains cover much of the UK. Such sites are attractive locations for landfills, since the till cover lowers the risk of landfill leachate entering groundwater. However, such tills often contain discrete sand and gravel bodies that can act as leachate flow routes. Such bodies may not be detected by conventional site investigation techniques such as drilling boreholes and trial pitting. A combined geoelectrical survey was carried out at a study site typical of such till plains and close to cliff exposures, which allowed direct mapping of sand bodies. Electrical resistivity tomography (ERT), resistivity cone penetrometry (RCPT) and frequency-domain electromagnetic data were collected. In a previous study, the electromagnetic and RCPT data were used to construct reference models for 2D inversion of the ERT data. The use of these reference models improved the solution models produced by inversion. We showed that the best solution model produced by inversion with a range of reference models could be determined without a priori knowledge of the true geoelectrical structure. This was done by using the area-weighted L2 norm between the solution models and associated reference models as a proxy for the misfit between the solution models and the true geoelectrical structure of the ground. In order to assess the most suitable method for combining invasive and non-invasive measurements, we compare both constrained and unconstrained 1D, 2D and 2.5D inversions of resistivity data collected at the study site. Preliminary results suggest that for 2.5D inversion the true 3D geoelectrical structure of the ground at the field study site is not sufficiently well known for comparison between the solution models and the true geoelectrical structure of the ground to be made. The results of work in progress evaluating layer-depth- constrained 1D inversion will be presented at the meeting.
Gállego, Isaac; Oncins, Gerard; Sisquella, Xavier; Fernàndez-Busquets, Xavier; Daban, Joan-Ramon
2010-01-01
In a previous study, we found that metaphase chromosomes are formed by thin plates, and here we have applied atomic force microscopy (AFM) and friction force measurements at the nanoscale (nanotribology) to analyze the properties of these planar structures in aqueous media at room temperature. Our results show that high concentrations of NaCl and EDTA and extensive digestion with protease and nuclease enzymes cause plate denaturation. Nanotribology studies show that native plates under structuring conditions (5 mM Mg2+) have a relatively high friction coefficient (μ ≈ 0.3), which is markedly reduced when high concentrations of NaCl or EDTA are added (μ ≈ 0.1). This lubricant effect can be interpreted considering the electrostatic repulsion between DNA phosphate groups and the AFM tip. Protease digestion increases the friction coefficient (μ ≈ 0.5), but the highest friction is observed when DNA is cleaved by micrococcal nuclease (μ ≈ 0.9), indicating that DNA is the main structural element of plates. Whereas nuclease-digested plates are irreversibly damaged after the friction measurement, native plates can absorb kinetic energy from the AFM tip without suffering any damage. These results suggest that plates are formed by a flexible and mechanically resistant two-dimensional network which allows the safe storage of DNA during mitosis. PMID:21156137
Computational modeling of neoclassical and resistive MHD tearing modes in tokamaks
Gianakon, T.A.; Hegna, C.C.; Callen, J.D.
1996-06-24
Numerical studies of the nonlinear evolution of MHD-type tearing modes in three-dimensional toroidal geometry with neoclassical effects are presented. The inclusion of neoclassical physics introduces an additional free-energy source for the nonlinear formation of magnetic islands through the effects of a bootstrap current in Ohm`s law. The neoclassical tearing mode is demonstrated to be destabilized in plasmas which are otherwise {Delta}` stable, albeit once an island width threshold is exceeded. The plasma pressure dynamics and neoclassical tearing growth is shown to be sensitive to the choice of the ratio of the parallel to perpendicular diffusivity ({Chi}{parallel}/{Chi}{perpendicular}). The study is completed with a demonstration and theoretical comparison of the threshold for single helicity neoclassical MHD tearing modes, which is described based on parameter scans of the local pressure gradient, the ratio of perpendicular to parallel pressure diffusivities {Chi}{perpendicular}/{Chi}{parallel}, and the magnitude of an initial seed magnetic perturbation.
NASA Astrophysics Data System (ADS)
Jinmin, M.; Saad, R.; Saidin, M.; Ismail, N. A.
2015-03-01
A 2-D resistivity imaging (2-DERI) study was conducted at Bukit Bunuh, Lenggong, Perak. Archaeological Global Research Centre, Universiti Sains Malaysia shows the field evidence of shock metamorphisms (suevite breccia) and crater morphology at Bukit Bunuh. A regional 2-DERI study focusing at Bukit Bunuh to identify the features of subsurface and detail study was then executed to verify boundary of the crater with the rebound effects at Bukit Bunuh which covered approximately 132.25 km2. 2-DERI survey used resistivity equipment by ABEM SAS4000 Terrameter and ES10-64C electrode slector with pole-dipole array. The survey lines were carried out using `roll-along' technique. The data were processed and analysed using RES2DINV, Excel and Surfer software to obtain resistivity results for qualitative interpretations. Bedrock depths were digitized from section by sections obtained. 2-DERI results gives both regional and detail study shows that the study area was divided into two main zones, overburden consists of alluvium mix with boulders embedded with resistivity value of 10-800 Ωm and granitic bedrock with resistivity value of >1500 Ωm and depth 5-50 m. The low level bedrock was circulated by high level bedrock (crater rim) was formed at the same area with few spots of high level bedrock which appeared at the centre of the rim which suspected as rebound zones (R). Assimilations of 2-DERI with boreholes are successful give valid and reliable results. The results of the study indicates geophysical method are capable to retrieve evidence of meteorite impact subsurface of the studied area.
3D MHD Simulations of Tokamak Disruptions
NASA Astrophysics Data System (ADS)
Woodruff, Simon; Stuber, James
2014-10-01
Two disruption scenarios are modeled numerically by use of the CORSICA 2D equilibrium and NIMROD 3D MHD codes. The work follows the simulations of pressure-driven modes in DIII-D and VDEs in ITER. The aim of the work is to provide starting points for simulation of tokamak disruption mitigation techniques currently in the CDR phase for ITER. Pressure-driven instability growth rates previously observed in simulations of DIIID are verified; Halo and Hiro currents produced during vertical displacements are observed in simulations of ITER with implementation of resistive walls in NIMROD. We discuss plans to exercise new code capabilities and validation.
NASA Astrophysics Data System (ADS)
Schaa, R.; Gross, L.; du Plessis, J.
2016-04-01
We present a general finite-element solver, escript, tailored to solve geophysical forward and inverse modeling problems in terms of partial differential equations (PDEs) with suitable boundary conditions. Escript’s abstract interface allows geoscientists to focus on solving the actual problem without being experts in numerical modeling. General-purpose finite element solvers have found wide use especially in engineering fields and find increasing application in the geophysical disciplines as these offer a single interface to tackle different geophysical problems. These solvers are useful for data interpretation and for research, but can also be a useful tool in educational settings. This paper serves as an introduction into PDE-based modeling with escript where we demonstrate in detail how escript is used to solve two different forward modeling problems from applied geophysics (3D DC resistivity and 2D magnetotellurics). Based on these two different cases, other geophysical modeling work can easily be realized. The escript package is implemented as a Python library and allows the solution of coupled, linear or non-linear, time-dependent PDEs. Parallel execution for both shared and distributed memory architectures is supported and can be used without modifications to the scripts.
Proceedings of the workshop on nonlinear MHD and extended MHD
1998-12-01
Nonlinear MHD simulations have proven their value in interpreting experimental results over the years. As magnetic fusion experiments reach higher performance regimes, more sophisticated experimental diagnostics coupled with ever expanding computer capabilities have increased both the need for and the feasibility of nonlinear global simulations using models more realistic than regular ideal and resistive MHD. Such extended-MHD nonlinear simulations have already begun to produce useful results. These studies are expected to lead to ever more comprehensive simulation models in the future and to play a vital role in fully understanding fusion plasmas. Topics include the following: (1) current state of nonlinear MHD and extended-MHD simulations; (2) comparisons to experimental data; (3) discussions between experimentalists and theorists; (4) /equations for extended-MHD models, kinetic-based closures; and (5) paths toward more comprehensive simulation models, etc. Selected papers have been indexed separately for inclusion in the Energy Science and Technology Database.
NASA Astrophysics Data System (ADS)
Freidberg, Jeffrey P.
2014-06-01
1. Introduction; 2. The ideal MHD model; 3. General properties of ideal MHD; 5. Equilibrium: one-dimensional configurations; 6. Equilibrium: two-dimensional configurations; 7. Equilibrium: three-dimensional configurations; 8. Stability: general considerations; 9. Alternate MHD models; 10. MHD stability comparison theorems; 11. Stability: one-dimensional configurations; 12. Stability: multi-dimensional configurations; Appendix A. Heuristic derivation of the kinetic equation; Appendix B. The Braginskii transport coefficients; Appendix C. Time derivatives in moving plasmas; Appendix D. The curvature vector; Appendix E. Overlap limit of the high b and Greene-Johnson stellarator models; Appendix F. General form for q(y); Appendix G. Natural boundary conditions; Appendix H. Upper and lower bounds on dQKIN.
Ding, Lin; Song, An; Dai, Meng; Xu, Min; Sun, Wanwan; Xu, Baihui; Sun, Jichao; Wang, Tiange; Xu, Yu; Lu, Jieli; Wang, Weiqing; Bi, Yufang; Ning, Guang
2015-01-01
Lipoprotein (a) [Lp(a)], an LDL-like particle, has been proposed as a causal risk factor for CVD among general populations. Meanwhile, both serum Lp(a) and diabetes increase the risk of CVD. However, the relationship between serum Lp(a) and T2D is poorly characterized, especially in the Asian population. Therefore, we conducted a cross-sectional study in 10,122 participants aged 40 years or older in Jiading District, Shanghai, China. Our study found that the prevalence of T2D was decreased from 20.9% to 15.0% from the lowest quartile to the highest quartile of serum Lp(a) concentrations (P for trend <0.0001). Logistic regression analyses showed that the odds ratios and 95% confidence intervals of prevalent T2D for quartiles 2–4 versus quartile 1 were 0.86 (0.73–1.01), 0.88 (0.75–1.04), and 0.76 (0.64–0.90) (P for trend = 0.0002), after adjustment for traditional confounding factors. Moreover, the risks for prevalent prediabetes, insulin resistance, and hyperinsulinemia were also decreased from the lowest to the top quartile. This inverse association between serum Lp(a) and T2D was not appreciably changed after we adjusted hypoglycemic medications or excluded the subjects with hypoglycemic and/or lipid-lowering agents and/or a history of self-reported CVD. PMID:25649924
NASA Astrophysics Data System (ADS)
Yi, Myeong-Jong; Sasaki, Yutaka
2015-11-01
Frequency-domain loop-loop electromagnetic (EM) methods are sensitive to the magnetic susceptibility of the Earth as well as its resistivity. Thus, inversion techniques have been used to simultaneously reconstruct both resistivity and susceptibility models from EM data. However, to take full advantage of inversion methods, calibration errors must be assessed and removed because ignoring them can result in misleading models. We present a multidimensional inversion method that jointly inverts EM and direct current (DC) resistivity data to derive offset errors as well as resistivity and susceptibility models, assuming that calibration errors can be represented by in-phase and quadrature offsets at each frequency. Addition of independent data such as DC data is effective for more accurately estimating the offsets, resulting in more reliable subsurface models. Synthetic examples involving small-loop EM data show that simultaneous inversion for resistivity and susceptibility is not stable, because of strong correlations between in-phase offset parameters and background susceptibility, but that the offsets are well determined when the data misfit is reduced rapidly in the early iteration step. Improvements achieved by joint inversion are mainly on the resistivity model. For airborne electromagnetic (AEM) data, the inversion process is stable, because AEM data are acquired using more loop-loop geometries and a wider range of frequencies. As a result, both the resistivity and susceptibility models are significantly improved by joint inversion.
Generalized reduced MHD equations
Kruger, S.E.; Hegna, C.C.; Callen, J.D.
1998-07-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 toroidal 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.
NASA Astrophysics Data System (ADS)
Downs, C. M.; Nowicki, R. S.; Rains, M. C.; Kruse, S.
2015-12-01
In west-central Florida, wetland and lake distribution is strongly controlled by karst landforms. Sandhill wetlands and lakes are sand-filled upland basins whose water levels are groundwater driven. Lake dimensions only reach wetland edges during extreme precipitation events. Current wetland classification schemes are inappropriate for identifying sandhill wetlands due to their unique hydrologic regime and ecologic expression. As a result, it is difficult to determine whether or not a wetland is impacted by groundwater pumping, development, and climate change. A better understanding of subsurface structures and how they control the hydrologic regime is necessary for development of an identification and monitoring protocol. Long-term studies record vegetation diversity and distribution, shallow ground water levels and surface water levels. The overall goals are to determine the hydrologic controls (groundwater, seepage, surface water inputs). Most recently a series of geophysical surveys was conducted at select sites in Hernando and Pasco County, Florida. Electrical resistivity and ground penetrating radar were employed to image sand-filled basins and the top of the limestone bedrock and stratigraphy of wetland slopes, respectively. The deepest extent of these sand-filled basins is generally reflected in topography as shallow depressions. Resistivity along inundated wetlands suggests the pools are surface expressions of the surficial aquifer. However, possible breaches in confining clay layers beneath topographic highs between depressions are seen in resistivity profiles as conductive anomalies and in GPR as interruptions in otherwise continuous horizons. These data occur at sites where unconfined and confined water levels are in agreement, suggesting communication between shallow and deep groundwater. Wetland plants are observed outside the historic wetland boundary at many sites, GPR profiles show near-surface layers dipping towards the wetlands at a shallower
NASA Astrophysics Data System (ADS)
Hirshman, S. P.; Shafer, M. W.; Seal, S. K.; Canik, J. M.
2016-04-01
> The SIESTA magnetohydrodynamic (MHD) equilibrium code has been used to compute a sequence of ideally stable equilibria resulting from numerical variation of the helical resonant magnetic perturbation (RMP) applied to an axisymmetric DIII-D plasma equilibrium. Increasing the perturbation strength at the dominant , resonant surface leads to lower MHD energies and increases in the equilibrium island widths at the (and sidebands) surfaces, in agreement with theoretical expectations. Island overlap at large perturbation strengths leads to stochastic magnetic fields which correlate well with the experimentally inferred field structure. The magnitude and spatial phase (around the dominant rational surfaces) of the resonant (shielding) component of the parallel current are shown to change qualitatively with the magnetic island topology.
NASA Astrophysics Data System (ADS)
Chang, S.; Chang, P.; Wu, Y.; Chen, J.; Huang, C.; Wang, Y.; Chen, W.
2011-12-01
In the study we try to characteristize the hydraulic boundary that separates the Pleistocene Toukeshan Formation from the Holocene sediments with Electrical Resistivity Imaging (ERI) method near the Bizetou Pass of the Choushui River. Comparing 2D Electrical resistivity imaging (ERI) with water level logs and core data from observation wells, we attempted to map the distribution of the shallow groundwater surfaces and the composition variations of the shallow unconfined aquifer from upper to lower fan across the Bizetou Pass. We found that the shallow groundwater surface drops sudenly from about 110m to 70 m after passing the Bizetou Pass from observation wells at the east and west side of the Pass and with about 3 km apart. The inverted resistivity images also show that the hydraulic gradient estimated from the resistivity images is getting larger to about 7% near the Bizetou Pass and then becoming less than 3% in the west the Pass. In addition, we found a significant change from about 1500 ohm-m to 500 ohm-m in the unconfined aquifer after passing the Bizetou Pass from the upstream side. The high resistivity value (about 1500 ohm-m) in the east of Bizetou Pass may represent the compacted Pleistocene Tokeshan sand and gravel layers, and the relatively low resistivity value (about 500 ohm-m) in the west of Bizetou Pass was the loose Holocene alluvium sediments. Since the previous study shows that the Changhua fault is not outcropped at the ground surface, our findings imply that the hydraulic boundary may be due to the less permeable Toukoshan Formation. To confirm this, the future work will extend the study area in North-South direction between Changhua County and Yunlin County.
NASA Astrophysics Data System (ADS)
Kellerer-Pirklbauer, Andreas
2015-04-01
Mountain permafrost covers some 2000 km² of the Austrian Alps which is less than 2.5% of the national territory. Delineating the altitudinal lower limit of permafrost in the mountains of Austria is difficult due the complex topography, the rather sparseness of field verification data and the lack of long-term permafrost monitoring data. Such monitoring data should cover different slope aspects, different elevations, different substrates and different mountain regions of Austria. In this study it was attempted to delineate the lower limit of permafrost at two study sites in the Tauern Range, Austria, applying two-dimensional geoelectrical resistivity tomography (ERT). In addition, multi-annual ground temperature data collected by miniature temperature datalogger (MDT) were used to validate the results. At the study site Hochreichart (maximum elevation 2416 m asl), located in the Seckauer Tauern Range, 14 ERT profiles (lengths 48-196 m; electrode spacing 2, 2.5 or 4 m) were measured at elevations between 1805 and 2416 m asl. Measurements were carried out at two cirques (Reichart, Schöneben) and at the summit plateau of Hochreichart. Results at this site indicate that permafrost lenses are detectable at elevations down to c.1900 m asl at radiation-sheltered sites. Furthermore, at the summit plateau permafrost only occurs as rather small lenses. The ERT-based permafrost pattern is generally confirmed by the MTD data with negative mean annual ground temperature values at only a few monitoring sites. However, the possibility of air-filled cavities causing higher resistive zones faking permafrost existence cannot be excluded because coarse-grained sediments (i.e. relict rock glaciers and autochthonous block fields) are widespread at this study site. At the second study site Kögele Cirque (maximum elevation 3030 m asl) located in the Schober Mountains 12 ERT profiles (lengths 48 m; electrode spacing 2 m) were measured at elevations between 2631 and 2740 m asl. Spatially
Kura, Nura Umar; Ramli, Mohammad Firuz; Ibrahim, Shaharin; Sulaiman, Wan Nor Azmin; Zaudi, Muhammad Amar; Aris, Ahmad Zaharin
2014-01-01
The existing knowledge regarding seawater intrusion and particularly upconing, in which both problems are linked to pumping, entirely relies on theoretical assumptions. Therefore, in this paper, an attempt is made to capture the effects of pumping on seawater intrusion and upconing using 2D resistivity measurement. For this work, two positions, one perpendicular and the other parallel to the sea, were chosen as profile line for resistivity measurement in the coastal area near the pumping wells of Kapas Island, Malaysia. Subsequently, water was pumped out of two pumping wells simultaneously for about five straight hours. Then, immediately after the pumping stopped, resistivity measurements were taken along the two stationed profile lines. This was followed by additional measurements after four and eight hours. The results showed an upconing with low resistivity of about 1-10 Ωm just beneath the pumping well along the first profile line that was taken just after the pumping stopped. The resistivity image also shows an intrusion of saline water (water enriched with diluted salt) from the sea coming towards the pumping well with resistivity values ranging between 10 and 25 Ωm. The subsequent measurements show the recovery of freshwater in the aquifer and how the saline water is gradually diluted or pushed out of the aquifer. Similarly the line parallel to the sea (L2) reveals almost the same result as the first line. However, in the second and third measurements, there were some significant variations which were contrary to the expectation that the freshwater may completely flush out the saline water from the aquifer. These two time series lines show that as the areas with the lowest resistivity (1 Ωm) shrink with time, the low resistivity (10 Ωm) tends to take over almost the entire area implying that the freshwater-saltwater equilibrium zone has already been altered. These results have clearly enhanced our current understanding and add more scientific
Ramli, Mohammad Firuz; Ibrahim, Shaharin; Sulaiman, Wan Nor Azmin; Aris, Ahmad Zaharin
2014-01-01
The existing knowledge regarding seawater intrusion and particularly upconing, in which both problems are linked to pumping, entirely relies on theoretical assumptions. Therefore, in this paper, an attempt is made to capture the effects of pumping on seawater intrusion and upconing using 2D resistivity measurement. For this work, two positions, one perpendicular and the other parallel to the sea, were chosen as profile line for resistivity measurement in the coastal area near the pumping wells of Kapas Island, Malaysia. Subsequently, water was pumped out of two pumping wells simultaneously for about five straight hours. Then, immediately after the pumping stopped, resistivity measurements were taken along the two stationed profile lines. This was followed by additional measurements after four and eight hours. The results showed an upconing with low resistivity of about 1–10 Ωm just beneath the pumping well along the first profile line that was taken just after the pumping stopped. The resistivity image also shows an intrusion of saline water (water enriched with diluted salt) from the sea coming towards the pumping well with resistivity values ranging between 10 and 25 Ωm. The subsequent measurements show the recovery of freshwater in the aquifer and how the saline water is gradually diluted or pushed out of the aquifer. Similarly the line parallel to the sea (L2) reveals almost the same result as the first line. However, in the second and third measurements, there were some significant variations which were contrary to the expectation that the freshwater may completely flush out the saline water from the aquifer. These two time series lines show that as the areas with the lowest resistivity (1 Ωm) shrink with time, the low resistivity (10 Ωm) tends to take over almost the entire area implying that the freshwater-saltwater equilibrium zone has already been altered. These results have clearly enhanced our current understanding and add more scientific
Resistive reduced MHD modeling of multi-edge-localized-mode cycles in Tokamak X-point plasmas.
Orain, F; Bécoulet, M; Huijsmans, G T A; Dif-Pradalier, G; Hoelzl, M; Morales, J; Garbet, X; Nardon, E; Pamela, S; Passeron, C; Latu, G; Fil, A; Cahyna, P
2015-01-23
The full dynamics of a multi-edge-localized-mode (ELM) cycle is modeled for the first time in realistic tokamak X-point geometry with the nonlinear reduced MHD code jorek. The diamagnetic rotation is found to be instrumental to stabilize the plasma after an ELM crash and to model the cyclic reconstruction and collapse of the plasma pressure profile. ELM relaxations are cyclically initiated each time the pedestal gradient crosses a triggering threshold. Diamagnetic drifts are also found to yield a near-symmetric ELM power deposition on the inner and outer divertor target plates, consistent with experimental measurements. PMID:25659004
Turco, F. Hanson, J. M.; Navratil, G. A.; Turnbull, A. D.
2015-02-15
Experiments conducted at DIII-D investigate the role of drift kinetic damping and fast neutral beam injection (NBI)-ions in the approach to the no-wall β{sub N} limit. Modelling results show that the drift kinetic effects are significant and necessary to reproduce the measured plasma response at the ideal no-wall limit. Fast neutral-beam ions and rotation play important roles and are crucial to quantitatively match the experiment. In this paper, we report on the model validation of a series of plasmas with increasing β{sub N}, where the plasma stability is probed by active magnetohydrodynamic (MHD) spectroscopy. The response of the plasma to an externally applied field is used to probe the stable side of the resistive wall mode and obtain an indication of the proximity of the equilibrium to an instability limit. We describe the comparison between the measured plasma response and that calculated by means of the drift kinetic MARS-K code [Liu et al., Phys. Plasmas 15, 112503 (2008)], which includes the toroidal rotation, the electron and ion drift-kinetic resonances, and the presence of fast particles for the modelled plasmas. The inclusion of kinetic effects allows the code to reproduce the experimental results within ∼13% for both the amplitude and phase of the plasma response, which is a significant improvement with respect to the undamped MHD-only model. The presence of fast NBI-generated ions is necessary to obtain the low response at the highest β{sub N} levels (∼90% of the ideal no-wall limit). The toroidal rotation has an impact on the results, and a sensitivity study shows that a large variation in the predicted response is caused by the details of the rotation profiles at high β{sub N}.
FLIP MHD - A particle-in-cell method for magnetohydrodynamics
NASA Technical Reports Server (NTRS)
Brackbill, J. U.
1991-01-01
The fluid-implicit-particle, or 'FLIP' method presently extended to 2D and 3D MHD flow incorporates a Lagrangian field representation and yields a grid magnetic Reynolds number of up to 16 while preserving contact continuities that retain the Galilean invariance of the MHD flow equations. Analytical arguments and numerical examples demonstrate the conservation of mass, momentum, magnetic flux, and energy; 2D calculation results for the illustrative cases of contact discontinuity convection, Rayleigh-Taylor unstable flow.
Heeter, R F; Fasoli, A; Testa, D; Sharapov, S; Berk, H L; Breizman, B; Gondhalekar, A; Mantsinen, M
2004-03-23
Experiments are conducted on the JET tokamak to assess the diagnostic potential of MHD active and passive spectroscopy, for the plasma bulk and its suprathermal components, using Alfv{acute e}n Eigenmodes (AEs) excited by external antennas and by energetic particles. The measurements of AE frequencies and mode numbers give information on the bulk plasma. Improved equilibrium reconstruction, in particular in terms of radial profiles of density and safety factor, is possible from the comparison between the antenna driven spectrum and that calculated theoretically. Details of the time evolution of the non-monotonic safety factor profile in advanced scenarios can be reconstructed from the frequency of ICRH-driven energetic particle modes. The plasma effective mass can be inferred from the resonant frequency of externally driven AEs in discharges with similar equilibrium profiles. The stability thresholds and the nonlinear development of the instabilities can give clues on energy and spatial distribution of the fast particle population. The presence of unstable AEs provides lower limits in the energy of ICRH generated fast ion tails. Fast ion pressure gradients and their evolution can be inferred from the stability of AEs at different plasma radial positions. Finally, the details of the AE spectrum in the nonlinear stage can be used to obtain information about the fast particle velocity space diffusion.
NASA Astrophysics Data System (ADS)
Vlahakis, Nektarios
2010-03-01
Outflows emanating from the environment of stellar or galactic objects are a widespread phenomenon in astrophysics. Their morphology ranges from nearly spherically symmetric winds to highly collimated jets. In some cases, e.g., in jets associated with young stellar objects, the bulk outflow speeds are nonrelativistic, while in others, e.g., in jets associated with active galactic nuclei or gamma-ray bursts, it can even be highly relativistic. The main driving mechanism of collimated outflows is likely related to magnetic fields. These fields are able to tap the rotational energy of the compact object or disk, accelerate, and collimate matter ejecta. To zeroth order these outflows can be described by the highly intractable theory of magnetohydrodynamics (MHD). Even in systems where the assumptions of zero resistivity (ideal MHD), steady state, axisymmetry, one fluid description, and polytropic equation of state are applicable, the problem remains difficult. In this case the problem reduces to only two equations, corresponding to the two components of the momentum equation along the flow and in the direction perpendicular to the magnetic field (transfield direction). The latter equation is the most difficult to solve, but also the most important. It answers the question on the degree of the collimation, but also crucially affects the solution of the first, the acceleration efficiency and the bulk velocity of the flow. The first and second parts of this chapter refer to nonrelativistic and relativistic flows, respectively. These Parts can be read independently. In each one, the governing equations are presented and discussed, focusing on the case of flows that are magnetically dominated near the central source. The general characteristics of the solutions in relation to the acceleration and collimation mechanisms are analyzed. As specific examples of exact solutions of the full system of the MHD equations that satisfy all the analyzed general characteristics, self
A 2D model to design MHD induction pumps
NASA Astrophysics Data System (ADS)
Stieglitz, R.; Zeininger, J.
2006-09-01
Technical liquid metal systems accompanied by a thermal transfer of energy such as reactor systems, metallurgical processes, metal refinement, casting, etc., require a forced convection of the fluid. The increased temperatures and more often the environmental conditions as, e.g., in a nuclear environment, pumping principles are required, in which rotating parts are absent. Additionally, in many applications a controlled atmosphere is indispensable, in order to ensure the structural integrity of the duct walls. An interesting option to overcome the sealing problem of a mechanical pump towards the surrounding is offered by induction systems. Although their efficiency compared to that of turbo machines is quite low, they have several advantages, which are attractive to the specific requirements in liquid metal applications such as: - low maintenance costs due to the absence of sealings, bearings and moving parts; - low degradation rate of the structural material; - simple replacement of the inductor without cut of the piping system; - fine regulation of flow rate by different inductor connections; - change of pump characteristics without change of the mechanical set-up. Within the article, general design requirements of electromagnetic pumps (EMP) are elaborated. The design of two annular linear induction pumps operating with sodium and lead-bismuth are presented and the calculated pump characteristics and experimentally obtained data are compared. In this context, physical effects leading to deviations between the model and the real data are addressed. Finally, the main results are summarized. Tables 4, Figs 4, Refs 12.
Wang, Chong; Liu, Chang-Ming; Wei, Li-Liang; Shi, Li-Ying; Pan, Zhi-Fen; Mao, Lian-Gen; Wan, Xiao-Chen; Ping, Ze-Peng; Jiang, Ting-Ting; Chen, Zhong-Liang; Li, Zhong-Jie; Li, Ji-Cheng
2016-01-01
The epidemic of pulmonary tuberculosis (TB), especially multidrug-resistance tuberculosis (MDR-TB) presented a major challenge for TB treatment today. We performed iTRAQ labeling coupled with two-dimensional liquid chromatography-tandem mass spectrometry (2D LC-MS/MS) and Solexa sequencing among MDR-TB patients, drug-sensitive tuberculosis (DS-TB) patients, and healthy controls. A total of 50 differentially expressed proteins and 43 differentially expressed miRNAs (fold change >1.50 or <0.60, P<0.05) were identified in the MDR-TB patients compared to both DS-TB patients and healthy controls. We found that 22.00% of differentially expressed proteins and 32.56% of differentially expressed miRNAs were related, and could construct a network mainly in complement and coagulation cascades. Significant differences in CD44 antigen (CD44), coagulation factor XI (F11), kininogen-1 (KNG1), miR-4433b-5p, miR-424-5p, and miR-199b-5p were found among MDR-TB patients, DS-TB patients and healthy controls (P<0.05) by enzyme-linked immunosorbent assay (ELISA) and SYBR green qRT-PCR validation. A strong negative correlation, consistent with the target gene prediction, was found between miR-199b-5p and KNG1 (r=-0.232, P=0.017). Moreover, we established the MDR-TB diagnostic model based on five biomarkers (CD44, KNG1, miR-4433b-5p, miR-424-5p, and miR-199b-5p). Our study proposes potential biomarkers for MDR-TB diagnosis, and also provides a new experimental basis to understand the pathogenesis of MDR-TB. PMID:26884721
ERIC Educational Resources Information Center
Kantrowitz, Arthur; Rosa, Richard J.
1975-01-01
Explains the operation of the Magnetohydrodynamic (MHD) generator and advantages of the system over coal, oil or nuclear powered generators. Details the development of MHD generators in the United States and Soviet Union. (CP)
2d PDE Linear Symmetric Matrix Solver
Energy Science and Technology Software Center (ESTSC)
1983-10-01
ICCG2 (Incomplete Cholesky factorized Conjugate Gradient algorithm for 2d symmetric problems) was developed to solve a linear symmetric matrix system arising from a 9-point discretization of two-dimensional elliptic and parabolic partial differential equations found in plasma physics applications, such as resistive MHD, spatial diffusive transport, and phase space transport (Fokker-Planck equation) problems. These problems share the common feature of being stiff and requiring implicit solution techniques. When these parabolic or elliptic PDE''s are discretized withmore » finite-difference or finite-element methods,the resulting matrix system is frequently of block-tridiagonal form. To use ICCG2, the discretization of the two-dimensional partial differential equation and its boundary conditions must result in a block-tridiagonal supermatrix composed of elementary tridiagonal matrices. The incomplete Cholesky conjugate gradient algorithm is used to solve the linear symmetric matrix equation. Loops are arranged to vectorize on the Cray1 with the CFT compiler, wherever possible. Recursive loops, which cannot be vectorized, are written for optimum scalar speed. For matrices lacking symmetry, ILUCG2 should be used. Similar methods in three dimensions are available in ICCG3 and ILUCG3. A general source containing extensions and macros, which must be processed by a pre-compiler to obtain the standard FORTRAN source, is provided along with the standard FORTRAN source because it is believed to be more readable. The pre-compiler is not included, but pre-compilation may be performed by a text editor as described in the UCRL-88746 Preprint.« less
NASA Astrophysics Data System (ADS)
Marchand, P.; Masson, J.; Chabrier, G.; Hennebelle, P.; Commerçon, B.; Vaytet, N.
2016-07-01
We develop a detailed chemical network relevant to calculate the conditions that are characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of potassium, sodium, and hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to nH = 1012 cm-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process. The multi-dimensional multi-species equilibrium abundance table and a copy of the code are only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (http://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/592/A18
2-D Imaging of Electron Temperature in Tokamak Plasmas
T. Munsat; E. Mazzucato; H. Park; C.W. Domier; M. Johnson; N.C. Luhmann Jr.; J. Wang; Z. Xia; I.G.J. Classen; A.J.H. Donne; M.J. van de Pol
2004-07-08
By taking advantage of recent developments in millimeter wave imaging technology, an Electron Cyclotron Emission Imaging (ECEI) instrument, capable of simultaneously measuring 128 channels of localized electron temperature over a 2-D map in the poloidal plane, has been developed for the TEXTOR tokamak. Data from the new instrument, detailing the MHD activity associated with a sawtooth crash, is presented.
He, Weiqi; Kuang, Yongqin; Xing, Xuemin; Simpson, Richard J; Huang, Haidong; Yang, Tao; Chen, Jingmin; Yang, Libin; Liu, Enyu; He, Weifeng; Gu, Jianwen
2014-05-01
Three-dimensional cell culture techniques can better reflect the in vivo characteristics of tumor cells compared with traditional monolayer cultures. Compared with their 2D counterparts, 3D-cultured tumor cells showed enhanced resistance to the cytotoxic T cell-mediated immune response. However, it remains unclear whether 3D-cultured tumor cells have an enhanced resistance to NK cell cytotoxicity. In this study, a total of 363 differentially expressed proteins were identified between the 2D- and 3D-cultured U251 cells by comparative proteomics, and an immune-associated protein-protein interaction (PPI) network based on these differential proteins was constructed by bioinformatics. Within the network, HLA-E, as a molecule for inhibiting NK cell activation, was significantly up-regulated in the 3D-cultured tumor cells. Then, we found that the 3D-cultured U251 cells exhibited potent resistance to NK cell cytotoxicity in vitro and were prone to tumor formation in vivo. The resistance of the 3D-cultured tumor cells to NK cell lysis was mediated by the HLA-E/NKG2A interaction because the administration of antibodies that block either HLA-E or NKG2A completely eliminated this resistance and significantly decreased tumor formation. Taken together, our findings indicate that HLA-E up-regulation in 3D-cultured cells may result in enhanced tumor resistance to NK cell-mediated immune response. PMID:24742303
Technology Transfer Automated Retrieval System (TEKTRAN)
The wheat stem rust resistance gene Sr6, present in several wheat cultivars, confers a high level of resistance against a wide range of races of Puccinia graminis f. sp. tritici. The inheritance and expression of the Sr6 gene is complex. Resistance conferred by Sr6 is influenced by temperature, lig...
Benchmarking a hybrid MHD/kinetic code with C-2 experimental data
NASA Astrophysics Data System (ADS)
Magee, Richard; Clary, Ryan; Dettrick, Sean; Korepanov, Sergey; Onofri, Marco; Smirnov, Artem; TAE Team
2013-10-01
The C-2 device creates field-reversed configuration (FRC) plasmas via the dynamic merging of two compact toroids and heated with neutral beams. Simulations of these plasmas are performed with Q2D - a hybrid MHD/Monte Carlo code that evolves the plasma according to the resistive MHD equations and treats the neutral beam injected fast ions as a minority kinetic species. Recent Q2D runs have resulted in testable predictions, namely that the axial profile of the fast ions is double-peaked, and charge-exchange neutrals are localized in pitch-angle. In some simulations, the fast particle population can induce magnetic fluctuations. These fluctuations are largest in the radial component, have a characteristic frequency approximately equal to the fast ion bounce frequency (f ~ 150 kHz), and a broad k spectrum. These fluctuations have the beneficial effect of smoothing out the double-peaked axial fast ion density profile, resulting in an increased fast ion density at the mid-plane. We will present results from a benchmarking study to quantitatively compare the results of Q2D runs to existing C-2 experimental data.
Numerical MHD codes for modeling astrophysical flows
NASA Astrophysics Data System (ADS)
Koldoba, A. V.; Ustyugova, G. V.; Lii, P. S.; Comins, M. L.; Dyda, S.; Romanova, M. M.; Lovelace, R. V. E.
2016-05-01
We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar coordinates, and fully 3D cylindrical coordinates. Viscosity and diffusivity are implemented in the code to control the accretion rate in the disk and the rate of penetration of the disk matter through the magnetic field lines. The code has been utilized for the numerical investigations of a number of different astrophysical problems, several examples of which are shown.
Magnetic reconnection in a compressible MHD plasma
Hesse, Michael; Zenitani, Seiji; Birn, Joachim
2011-04-15
Using steady-state resistive MHD, magnetic reconnection is reinvestigated for conditions of high resistivity/low magnetic Reynolds number, when the thickness of the diffusion region is no longer small compared to its length. Implicit expressions for the reconnection rate and other reconnection parameters are derived based on the requirements of mass, momentum, and energy conservation. These expressions are solved via simple iterative procedures. Implications specifically for low Reynolds number/high resistivity are being discussed.
MHD simulations on an unstructured mesh
Strauss, H.R.; Park, W.; Belova, E.; Fu, G.Y.; Longcope, D.W.; Sugiyama, L.E.
1998-12-31
Two reasons for using an unstructured computational mesh are adaptivity, and alignment with arbitrarily shaped boundaries. Two codes which use finite element discretization on an unstructured mesh are described. FEM3D solves 2D and 3D RMHD using an adaptive grid. MH3D++, which incorporates methods of FEM3D into the MH3D generalized MHD code, can be used with shaped boundaries, which might be 3D.
MHD Turbulence and Magnetic Dynamos
NASA Technical Reports Server (NTRS)
Shebalin, John V
2014-01-01
Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much
FEMHD: An adaptive finite element method for MHD and edge modelling
Strauss, H.R.
1995-07-01
This paper describes the code FEMHD, an adaptive finite element MHD code, which is applied in a number of different manners to model MHD behavior and edge plasma phenomena on a diverted tokamak. The code uses an unstructured triangular mesh in 2D and wedge shaped mesh elements in 3D. The code has been adapted to look at neutral and charged particle dynamics in the plasma scrape off region, and into a full MHD-particle code.
Energy Science and Technology Software Center (ESTSC)
2005-07-01
Aniso2d is a two-dimensional seismic forward modeling code. The earth is parameterized by an X-Z plane in which the seismic properties Can have monoclinic with x-z plane symmetry. The program uses a user define time-domain wavelet to produce synthetic seismograms anrwhere within the two-dimensional media.
NASA Astrophysics Data System (ADS)
Takizuka, T.
2014-09-01
The tokamak hybrid scenario, a combination of the inductive current drive and the non-inductive current drive, aims at high neutron fluence and a long pulse length for the ITER engineering test towards DEMO. It has been reported experimentally that the concentration of current density, J, peculiar to inductive discharges, was much smaller than that of a simple estimation. We propose a new model to explain the above observations: the central core plasma column can be non-axisymmetric due to the m/n = 1/1 MHD-mode deformation (m and n are poloidal and toroidal mode numbers, respectively). The neo-classical resistivity is enhanced, and the diagnostics are affected. An internal kink mode is a candidate of this instability. Compared with the JT-60U experiment, the reconstructed J profile agrees well with the calculated J profile including the enhanced neo-classical resistivity.
Kang, Dongwei; Fang, Zengjun; Li, Zhenyu; Huang, Boshi; Zhang, Heng; Lu, Xueyi; Xu, Haoran; Zhou, Zhongxia; Ding, Xiao; Daelemans, Dirk; De Clercq, Erik; Pannecouque, Christophe; Zhan, Peng; Liu, Xinyong
2016-09-01
We designed and synthesized a series of human immunodeficiency virus type 1 (HIV-1) non-nucleoside reverse transcriptase inhibitors (NNRTIs) with a piperidine-substituted thiophene[3,2-d]pyrimidine scaffold, employing a strategy of structure-based molecular hybridization and substituent decorating. Most of the synthesized compounds exhibited broad-spectrum activity with low (single-digit) nanomolar EC50 values toward a panel of wild-type (WT), single-mutant, and double-mutant HIV-1 strains. Compound 27 was the most potent; compared with ETV, its antiviral efficacy was 3-fold greater against WT, 5-7-fold greater against Y181C, Y188L, E138K, and F227L+V106A, and nearly equipotent against L100I and K103N, though somewhat weaker against K103N+Y181C. Importantly, 27 has lower cytotoxicity (CC50 > 227 μM) and a huge selectivity index (SI) value (ratio of CC50/EC50) of >159101. 27 also showed favorable, drug-like pharmacokinetic and safety properties in rats in vivo. Molecular docking studies and the structure-activity relationships provide important clues for further molecular elaboration. PMID:27541578
Greg Flach, Frank Smith
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assigns an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.
Energy Science and Technology Software Center (ESTSC)
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assignsmore » an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.« less
NASA Astrophysics Data System (ADS)
Lotsch, Bettina V.
2015-07-01
Graphene's legacy has become an integral part of today's condensed matter science and has equipped a whole generation of scientists with an armory of concepts and techniques that open up new perspectives for the postgraphene area. In particular, the judicious combination of 2D building blocks into vertical heterostructures has recently been identified as a promising route to rationally engineer complex multilayer systems and artificial solids with intriguing properties. The present review highlights recent developments in the rapidly emerging field of 2D nanoarchitectonics from a materials chemistry perspective, with a focus on the types of heterostructures available, their assembly strategies, and their emerging properties. This overview is intended to bridge the gap between two major—yet largely disjunct—developments in 2D heterostructures, which are firmly rooted in solid-state chemistry or physics. Although the underlying types of heterostructures differ with respect to their dimensions, layer alignment, and interfacial quality, there is common ground, and future synergies between the various assembly strategies are to be expected.
Fully Parallel MHD Stability Analysis Tool
NASA Astrophysics Data System (ADS)
Svidzinski, Vladimir; Galkin, Sergei; Kim, Jin-Soo; Liu, Yueqiang
2014-10-01
Progress on full parallelization of the plasma stability code MARS will be reported. MARS calculates eigenmodes in 2D axisymmetric toroidal equilibria in MHD-kinetic plasma models. It is a powerful tool for studying MHD and MHD-kinetic instabilities and it is widely used by fusion community. Parallel version of MARS is intended for simulations on local parallel clusters. It will be an efficient tool for simulation of MHD instabilities with low, intermediate and high toroidal mode numbers within both fluid and kinetic plasma models, already implemented in MARS. Parallelization of the code includes parallelization of the construction of the matrix for the eigenvalue problem and parallelization of the inverse iterations algorithm, implemented in MARS for the solution of the formulated eigenvalue problem. Construction of the matrix is parallelized by distributing the load among processors assigned to different magnetic surfaces. Parallelization of the solution of the eigenvalue problem is made by repeating steps of the present MARS algorithm using parallel libraries and procedures. Initial results of the code parallelization will be reported. Work is supported by the U.S. DOE SBIR program.
Fully Parallel MHD Stability Analysis Tool
NASA Astrophysics Data System (ADS)
Svidzinski, Vladimir; Galkin, Sergei; Kim, Jin-Soo; Liu, Yueqiang
2013-10-01
Progress on full parallelization of the plasma stability code MARS will be reported. MARS calculates eigenmodes in 2D axisymmetric toroidal equilibria in MHD-kinetic plasma models. It is a powerful tool for studying MHD and MHD-kinetic instabilities and it is widely used by fusion community. Parallel version of MARS is intended for simulations on local parallel clusters. It will be an efficient tool for simulation of MHD instabilities with low, intermediate and high toroidal mode numbers within both fluid and kinetic plasma models, already implemented in MARS. Parallelization of the code includes parallelization of the construction of the matrix for the eigenvalue problem and parallelization of the inverse iterations algorithm, implemented in MARS for the solution of the formulated eigenvalue problem. Construction of the matrix is parallelized by distributing the load among processors assigned to different magnetic surfaces. Parallelization of the solution of the eigenvalue problem is made by repeating steps of the present MARS algorithm using parallel libraries and procedures. Preliminary results of the code parallelization will be reported. Work is supported by the U.S. DOE SBIR program.
Fully Parallel MHD Stability Analysis Tool
NASA Astrophysics Data System (ADS)
Svidzinski, Vladimir; Galkin, Sergei; Kim, Jin-Soo; Liu, Yueqiang
2015-11-01
Progress on full parallelization of the plasma stability code MARS will be reported. MARS calculates eigenmodes in 2D axisymmetric toroidal equilibria in MHD-kinetic plasma models. It is a powerful tool for studying MHD and MHD-kinetic instabilities and it is widely used by fusion community. Parallel version of MARS is intended for simulations on local parallel clusters. It will be an efficient tool for simulation of MHD instabilities with low, intermediate and high toroidal mode numbers within both fluid and kinetic plasma models, already implemented in MARS. Parallelization of the code includes parallelization of the construction of the matrix for the eigenvalue problem and parallelization of the inverse iterations algorithm, implemented in MARS for the solution of the formulated eigenvalue problem. Construction of the matrix is parallelized by distributing the load among processors assigned to different magnetic surfaces. Parallelization of the solution of the eigenvalue problem is made by repeating steps of the present MARS algorithm using parallel libraries and procedures. Results of MARS parallelization and of the development of a new fix boundary equilibrium code adapted for MARS input will be reported. Work is supported by the U.S. DOE SBIR program.
Dynamo theory and liquid metal MHD experiments
NASA Astrophysics Data System (ADS)
Lielausis, O.
1994-06-01
High values of magnetic Reynolds number Rm are characteristic not only to astrophysics, but also to other interesting objects, including liquid metal (LM) flows. LM experiments have been performed illustrating important predictions of the dynamo theory, for example, about the existence and features of the alpha effect. Consideration of so called 'laminar' dynamos provides a theoretical base for direct experimental realization and examination of the dynamo process. First step results, gathered a subcritical conditions, confirm the statement that self-excitation in LM experiments can be achieved practically today. In such devices as LM (sodium) cooled fast breeders Rm can reach values of up to 50 and specific MHD phenomena have been observed in operating fast reactors. Cautions against crisis like processes have been expressed. It is important for the dynamo theory to understand what kind of perturbed motion is able to coexist with the generated magnetic field. Fundamentally new ideas here are issuing from the theory of 2D MHD turbulence. LM MHD served for the first direct proves, confirming, that the predicted surprising features of 2D turbulence can be observed in reality. It is worth incorporating these already not new ideas in the dynamo theory. In such a way a field for new solutions could be established.
Wallin, Erin L.; Johnson, Timothy C.; Greenwood, William J.; Zachara, John M.
2013-03-29
The Hanford 300 Area is located adjacent to the Columbia River in south-central Washington State, USA, and was a former site for nuclear fuel processing operations. Waste disposal practices resulted in persistent unsaturated zone and groundwater contamination, the primary contaminant of concern being uranium. Uranium behavior at the site is intimately linked with river stage driven groundwater-river water exchange such that understanding the nature of river water intrusion into the 300 Area is critical for predicting uranium desorption and transport. In this paper we use time-lapse electrical resistivity tomography (ERT) to image the inland intrusion of river during high stage conditions. We demonstrate a modified time-lapse inversion approach, whereby the transient water table elevation is explicitly modeled by removing regularization constraints across the water table boundary. This implementation was critical for producing meaningful imaging results. We inverted approximately 1200 data sets (400 per line over 3 lines) using high performance computing resources to produce a time-lapse sequence of changes in bulk conductivity caused by river water intrusion during the 2011 spring runoff cycle over approximately 125 days. The resulting time series for each mesh element was then analyzed using common time series analysis to reveal the timing and location of river water intrusion beneath each line. The results reveal non-uniform flows characterized by preferred flow zones where river water enters and exits quickly with stage increase and decrease, and low permeability zones with broader bulk conductivity ‘break through’ curves and longer river water residence times. The time-lapse ERT inversion approach removes the deleterious effects of changing water table elevation and enables remote and spatial continuous groundwater-river water exchange monitoring using surface based ERT arrays under conditions where groundwater and river water conductivity are in contrast.
Survey of MHD plant applications
NASA Technical Reports Server (NTRS)
Lynch, J. J.; Seikel, G. R.; Cutting, J. C.
1979-01-01
Open-cycle MHD is one of the major R&D efforts in the Department of Energy's program to meet the national goal of reducing U.S. dependence on oil through increased utilization of coal. MHD offers an effective way to use coal to produce electric power at low cost in a highly efficient and environmentally acceptable manner. Open-cycle MHD plants are categorized by the MHD combustor oxidizer, its temperature and the method of preheat. The paper discusses MHD baseline plant design, open-cycle MHD plant in the Energy Conversion Alternatives Study (ECAS), early commercial MHD plants, conceptual studies of the engineering test facility, retrofit (addition of an MHD topping cycle to an existing steam plant), and other potential applications and concepts. Emphasis is placed on a survey of both completed and ongoing studies to define both commercial and pilot plant design, cost, and performance.
Extended MHD Simulations of Spheromaks
NASA Astrophysics Data System (ADS)
Howell, E. C.; Sovinec, C. R.
2012-10-01
Nonlinear extended MHD simulations of a spheromak in a cylindrical flux conserver are performed using the NIMROD code (JCP 195, 2004). An idealized series of simulations starting from a Grad-Shafranov equilibrium and small non-axisymmetric perturbations are performed to model the sustained decay phase. The resulting confinement leads to steep resistivity gradients. Strong current gradients develop, driving tearing modes that dominate the evolution of the spheromak. Absent in these simulations are the remains of n=1 fluctuations created during the formation process. A second series of simulations start from vacuum fields and model the full spheromak evolution, including the formation process where the n=1 fluctuations dominate. To understand the role of pressure driven instabilities in the evolution of the spheromak, a numerical diagnostic is developed to calculate the Mercier stability criterion from the axisymmetric fields.
Petrick, Michael; Pierson, Edward S.; Schreiner, Felix
1980-01-01
According to the present invention, coal combustion gas is the primary working fluid and copper or a copper alloy is the electrodynamic fluid in the MHD generator, thereby eliminating the heat exchangers between the combustor and the liquid-metal MHD working fluids, allowing the use of a conventional coalfired steam bottoming plant, and making the plant simpler, more efficient and cheaper. In operation, the gas and liquid are combined in a mixer and the resulting two-phase mixture enters the MHD generator. The MHD generator acts as a turbine and electric generator in one unit wherein the gas expands, drives the liquid across the magnetic field and thus generates electrical power. The gas and liquid are separated, and the available energy in the gas is recovered before the gas is exhausted to the atmosphere. Where the combustion gas contains sulfur, oxygen is bubbled through a side loop to remove sulfur therefrom as a concentrated stream of sulfur dioxide. The combustor is operated substoichiometrically to control the oxide level in the copper.
NASA Technical Reports Server (NTRS)
Retallick, F. D.
1980-01-01
Directly-fired, separately-fired, and oxygen-augmented MHD power plants incorporating a disk geometry for the MHD generator were studied. The base parameters defined for four near-optimum-performance MHD steam power systems of various types are presented. The finally selected systems consisted of (1) two directly fired cases, one at 1920 K (2996F) preheat and the other at 1650 K (2500 F) preheat, (2) a separately-fired case where the air is preheated to the same level as the higher temperature directly-fired cases, and (3) an oxygen augmented case with the same generator inlet temperature of 2839 (4650F) as the high temperature directly-fired and separately-fired cases. Supersonic Mach numbers at the generator inlet, gas inlet swirl, and constant Hall field operation were specified based on disk generator optimization. System pressures were based on optimization of MHD net power. Supercritical reheat stream plants were used in all cases. Open and closed cycle component costs are summarized and compared.
NASA Astrophysics Data System (ADS)
Silva, Pablo G.
2010-05-01
The ancient roman city of Baelo Claudia has been subject of several papers on earthquake environmental effects (EEE) and well as earthquake archaeological effects (EAE). During the field training course on archaeoseismology and palaeoseismology conducted in September 2009 (INQUA-IGCP567 Workshop) held at Baelo Claudia, four Electric Resistivity Tomography (ERT) profiles were carried out, by the teams of the Salamanca University (Spain), RWTH Aachen University (Germany) and the Geological Survey of Spain (IGME). ERT surveys were developed in the eastern side of the ancient roman Forum across the unexcavated sector of the archaeological site heading on the 1st Century AD Isis Temple. Each ERT profile was constituted by a 48 multielectrode array with spacing of 2 m resulting in a total length of investigation of around 384 m. ERT lines were separated 10 m each other resulting in a total research area of 3840 m2 to a mean investigation depth of 16 m. The selected survey configurations were Pole-Dipole and Wenner in order to get detailed information about lateral resistivity contrasts, but with a reasonable depth of investigation. The resulting 2D resistivity pseudosections clearly display deformations of the buried roman pavements which propagated in depth within the pre-roman clayey substratum of the Bolonia Bay area.. 3D modelling of the 2D pseudosections indicates that the observed deformations are related to near-surface landsliding, being possible to calculate the minimum volume of mobilized material. ERT 3D imaging allow to refine previous GPR surveys conducted at this same area and to get a subsurface picture of ground deformations caused by repeated earthquakes during the 1st and 3rd Centuries AD. Preliminary calculated volume for the mobilized materials affecting the foundations of the Isis Temple and Forum clearly points to a minimum ESI-07 VIII Intensity validating previous research in the zone. This study has been supported by the Spanish Research Projects
Extended-MHD modeling of diamagnetic-drift tearing instabilities
NASA Astrophysics Data System (ADS)
King, Jacob; Kruger, Scott
2012-03-01
We use analytics and computations with the NIMROD code to examine tearing stability in large-guide-field slab cases with a nonzero equilibrium pressure gradient. A well known result from drift-reduced MHD is the diamagnetic drift associated with the pressure gradient has a stabilizing influence were the dispersion relation becomes (γ+iφ*e)^3γ(γ+iφ*i)=γrMHD^5 [1]. Here φ*i and φ*e are the ion- and electron-diamagnetic frequencies and γrMHD is the tearing growth rate with a resistive-MHD model. Preliminary computational results with an unreduced extended-MHD model do not produce the expected drift-reduced result. For moderate values of φ*i (φ*i<=3γrMHD), the computations follow the dispersion relation that would result if the ∇pe term were not included in the drift-reduced parallel Ohm's law: (γ+iφ*e)^4(γ+iφ*i)=γrMHD^5. Analytics, guided by computational diagnostics, are used to examine the significant terms in the flux evolution equation and investigate the discrepancy with the drift-reduced result.[4pt] [1] For example Coppi, PoF 7, 1501 (1964); Biskamp, NF 18, 1059 (1978).
NASA Astrophysics Data System (ADS)
Schüssler, M.
1999-05-01
Two aspects of solar MHD are discussed in relation to the work of the MHD simulation group at KIS. Photospheric magneto-convection, the nonlinear interaction of magnetic field and convection in a strongly stratified, radiating fluid, is a key process of general astrophysical relevance. Comprehensive numerical simulations including radiative transfer have significantly improved our understanding of the processes and have become an important tool for the interpretation of observational data. Examples of field intensification in the solar photosphere ('convective collapse') are shown. The second line of research is concerned with the dynamics of flux tubes in the convection zone, which has far-reaching implications for our understanding of the solar dynamo. Simulations indicate that the field strength in the region where the flux is stored before erupting to form sunspot groups is of the order of 105 G, an order of magnitude larger than previous estimates based on equipartition with the kinetic energy of convective flows.
NASA Astrophysics Data System (ADS)
Morrison, P. J.; Abdelhamid, H. M.; Grasso, D.; Hazeltine, R. D.; Lingam, M.; Tassi, E.
2015-11-01
Over the years various reduced fluid models have been obtained for modeling plasmas, with the goal of capturing important physics while maintaining computability. Such models have included the physics contained in various generalizations of Ohm's law, including Hall drift and electron inertia. In a recent publication it was shown that full 3D extended MHD is a Hamiltonian system by finding its noncanonical Poisson bracket. Subsequently, this bracket was shown to be derivable from that for Hall MHD by a series of remarkable transformations, which greatly simplifies the proof of the Jacobi identity and allows one to immediately obtain generalizations of the helicity and cross helicity. In this poster we use this structure to obtain exact reduced fluid models with the effects of full two-fluid theory. Results of numerical computations of collisionless reconnection using an exact reduced 4-field model will be presented and analytical comparisons of mode structure of previous reduced models will be made.
NASA Astrophysics Data System (ADS)
Wang, Jin; Ma, Jianyong; Zhou, Changhe
2014-11-01
A 3×3 high divergent 2D-grating with period of 3.842μm at wavelength of 850nm under normal incidence is designed and fabricated in this paper. This high divergent 2D-grating is designed by the vector theory. The Rigorous Coupled Wave Analysis (RCWA) in association with the simulated annealing (SA) is adopted to calculate and optimize this 2D-grating.The properties of this grating are also investigated by the RCWA. The diffraction angles are more than 10 degrees in the whole wavelength band, which are bigger than the traditional 2D-grating. In addition, the small period of grating increases the difficulties of fabrication. So we fabricate the 2D-gratings by direct laser writing (DLW) instead of traditional manufacturing method. Then the method of ICP etching is used to obtain the high divergent 2D-grating.
Stochastic Inversion of 2D Magnetotelluric Data
Energy Science and Technology Software Center (ESTSC)
2010-07-01
The algorithm is developed to invert 2D magnetotelluric (MT) data based on sharp boundary parametrization using a Bayesian framework. Within the algorithm, we consider the locations and the resistivity of regions formed by the interfaces are as unknowns. We use a parallel, adaptive finite-element algorithm to forward simulate frequency-domain MT responses of 2D conductivity structure. Those unknown parameters are spatially correlated and are described by a geostatistical model. The joint posterior probability distribution function ismore » explored by Markov Chain Monte Carlo (MCMC) sampling methods. The developed stochastic model is effective for estimating the interface locations and resistivity. Most importantly, it provides details uncertainty information on each unknown parameter. Hardware requirements: PC, Supercomputer, Multi-platform, Workstation; Software requirements C and Fortan; Operation Systems/version is Linux/Unix or Windows« less
Stochastic Inversion of 2D Magnetotelluric Data
Chen, Jinsong
2010-07-01
The algorithm is developed to invert 2D magnetotelluric (MT) data based on sharp boundary parametrization using a Bayesian framework. Within the algorithm, we consider the locations and the resistivity of regions formed by the interfaces are as unknowns. We use a parallel, adaptive finite-element algorithm to forward simulate frequency-domain MT responses of 2D conductivity structure. Those unknown parameters are spatially correlated and are described by a geostatistical model. The joint posterior probability distribution function is explored by Markov Chain Monte Carlo (MCMC) sampling methods. The developed stochastic model is effective for estimating the interface locations and resistivity. Most importantly, it provides details uncertainty information on each unknown parameter. Hardware requirements: PC, Supercomputer, Multi-platform, Workstation; Software requirements C and Fortan; Operation Systems/version is Linux/Unix or Windows
3-D Relativistic MHD Simulations of Extragalactic Jets
NASA Astrophysics Data System (ADS)
Nishikawa, K.-I.; Koide, S.; Sakai, J.-I.; Frank, J.; Christodoulou, D. M.; Sol, H.; Mutel, R. L.
1997-12-01
We present the numerical simulations of relativistic jets propagating initially oblique to the field lines of a magnetized ambient medium. Our simulations incorporate relativistic MHD in a four-dimensional spacetime and clearly show that (a) relatively weak, oblique fields (at 1/16 of the equipartition value) have only a negligible influence on the propagating jet and they are passively pushed away by the relativistically moving head; (b) oblique fields in equipartition with the ambient plasma provide more resistance and cause bending at the jet head, but the magnitude of this deflection and the associated backflow are small compared to those identified by previous studies with a 2-D slab model. The new results are understood as follows: Relativistic simulations have consistently shown that these jets are effectively heavy and so they do not suffer substantial momentum losses and are not decelerated as efficiently as their nonrelativistic counterparts. In addition, the ambient magnetic field, however strong, can be pushed aside with relative ease by the beam, provided that the degrees of freedom associated with all three spatial dimensions are followed self-consistently during the simulations. The effect is analogous to pushing Japanese ``noren'' or vertical Venetian blinds out of the way while the slats are allowed to bend in 3-D space rather than as a 2-D slab structure. Applied to relativistic extragalactic jets from blazars, the new results are encouraging since superluminal outflows exhibit bending near their sources and their environments are profoundly magnetized---but observations do not provide support for irregular kinematics such as large-scale vortical motions and pronounced reverse flows near the points of origin.
NASA Astrophysics Data System (ADS)
Stone, James M.; Norman, Michael L.
1992-06-01
A detailed description of ZEUS-2D, a numerical code for the simulation of fluid dynamical flows including a self-consistent treatment of the effects of magnetic fields and radiation transfer is presented. Attention is given to the hydrodynamic (HD) algorithms which form the foundation for the more complex MHD and radiation HD algorithms. The effect of self-gravity on the flow dynamics is accounted for by an iterative solution of the sparse-banded matrix resulting from discretizing the Poisson equation in multidimensions. The results of an extensive series of HD test problems are presented. A detailed description of the MHD algorithms in ZEUS-2D is presented. A new method of computing the electromotive force is developed using the method of characteristics (MOC). It is demonstrated through the results of an extensive series of MHD test problems that the resulting hybrid MOC-constrained transport method provides for the accurate evolution of all modes of MHD wave families.
FLASH MHD simulations of experiments that study shock-generated magnetic fields
NASA Astrophysics Data System (ADS)
Tzeferacos, P.; Fatenejad, M.; Flocke, N.; Graziani, C.; Gregori, G.; Lamb, D. Q.; Lee, D.; Meinecke, J.; Scopatz, A.; Weide, K.
2015-12-01
We summarize recent additions and improvements to the high energy density physics capabilities in FLASH, highlighting new non-ideal magneto-hydrodynamic (MHD) capabilities. We then describe 3D Cartesian and 2D cylindrical FLASH MHD simulations that have helped to design and analyze experiments conducted at the Vulcan laser facility. In these experiments, a laser illuminates a carbon rod target placed in a gas-filled chamber. A magnetic field diagnostic (called a Bdot) employing three very small induction coils is used to measure all three components of the magnetic field at a chosen point in space. The simulations have revealed that many fascinating physical processes occur in the experiments. These include megagauss magnetic fields generated by the interaction of the laser with the target via the Biermann battery mechanism, which are advected outward by the vaporized target material but decrease in strength due to expansion and resistivity; magnetic fields generated by an outward expanding shock via the Biermann battery mechanism; and a breakout shock that overtakes the first wave, the contact discontinuity between the target material and the gas, and then the initial expanding shock. Finally, we discuss the validation and predictive science we have done for this experiment with FLASH.
Toroidal Theory of MHD Instabilities
Goedbloed, J.P.
2004-03-15
We continue with the adventures of the Alfven wave and its two magnetosonic companions as they travel in the curved space of magnetic surfaces and field lines (Sec. 2), find themselves trapped in singularities of an unprecedented richness (Sec. 3), decide to get themselves better maps of the landscape to do the required twisting while some of their youthful energy is leaking away (Sec. 4), cause trouble at the edge of a powerful empire (Sec. 5), and finally see the light in a distant future (Sec. 6). Needed on the trip are the evolution equations of both ideal and resistive MHD 'derived' in reference [1], the solutions to the toroidal equilibrium equations discussed in reference [2], the general background on spectral theory of inhomogeneous plasmas presented in reference [3], which is extended in the two directions of toroidal geometry and resistivity in this lecture [4]. This leads to such intricate dynamics that numerical techniques are virtually the only way to proceed. This aspect is further elaborated in reference [5] on numerical techniques.
MHD Modeling in Complex 3D Geometries: Towards Predictive Simulation of SIHI Current Drive
NASA Astrophysics Data System (ADS)
Hansen, Christopher James
The HIT-SI experiment studies Steady Inductive Helicity Injection (SIHI) for the purpose of forming and sustaining a spheromak plasma. A spheromak is formed in a nearly axisymmetric flux conserver, with a bow tie cross section, by means of two semi-toroidal injectors. The plasma-facing surfaces of the device, which are made of copper for its low resistivity, are covered in an insulating coating in order to operate in a purely inductive manner. Following formation, the spheromak flux and current are increased during a quiescent period marked by a decrease in the global mode activity. A proposed mechanism, Imposed Dynamo Current Drive (IDCD), is expected to be responsible for this phase of quiescent current drive. Due to the geometric complexity of the experiment, previous numerical modeling efforts have used a simplified geometry that excludes the injector volumes from the simulated domain. The effect of helicity injection is then modeled by boundary conditions on this reduced plasma volume. The work presented here has explored and developed more complete computational models of the HIT-SI device. This work is separated into 3 distinct but complementary areas: 1) Development of a 3D MHD equilibrium code that can incorporate the non-axisymmetric injector fields present in HIT-SI and investigation of equilibria of interest during spheromak sustainment. 2) A 2D axisymmetric MHD equilibrium code that was used to explore reduced order models for mean-field evolution using equations derived from IDCD theory including coupling to 3D equilibria. 3) A 3D time-dependent non-linear MHD code that is capable of modeling the entire plasma volume including dynamics within the injectors. Although HIT-SI was the motivation for, and experiment studied in this research, the tools and methods developed are general --- allowing their application to a broad range of magnetic confinement experiments. These tools constitute a significant advance for modeling plasma dynamics in devices with
Implicit Predictor-Corrector finite difference scheme for the ideal MHD simulations
NASA Astrophysics Data System (ADS)
Tsai, T.; Yu, H.; Lai, S.
2012-12-01
A innovative simulation code for ideal magnetohydrodynamics (MHD) is developed. We present a multiple-dimensional MHD code based on high-order implicit predictor-corrector finite difference scheme (high-order IPCFD scheme). High-order IPCFD scheme adopts high-order predictor-corrector scheme for the time integration and high-order central difference method as the spatial derivative solver. We use Elimination-of-the-Runoff-Errors (ERE) technology to avoid the numerical oscillations and numerical instability in the simulation results. In one-dimensional MHD problem, our simulation results show good agreement with the Brio & Wu MHD shock tube problem. The divergent B constraint remains fully satisfied, that is the divergent B equals to zero throughout the simulation. When solving the two-dimensional (2D) linear wave in MHD plasma, we clearly obtain the group-velocity Friedrichs diagrams of the MHD waves. Here we demonstrate 2D simulation results of rotor problem, Orszag-Tang vortex system, vortex type K-H instability, and kink type K-H instability by using our IPCFD MHD code and discuss the advantage of our simulation code.
Kinetic MHD simulation of large 'circ; tearing mode
NASA Astrophysics Data System (ADS)
Cheng, Jianhua; Chen, Yang; Parker, Scott; Uzdensky, Dmitri
2012-03-01
We have developed a second-order accurate semi-implicit δ method for kinetic MHD simulation with Lorentz force ions and fluid electrons. The model has been used to study the resistive tearing mode instability, which involves multiple spatial scales. In small 'circ; cases, the linear growth rate and eigenmode structure are consistent with resistive MHD analysis. The Rutherford stage and saturation are demonstrated, but the simulation exhibits different saturation island widths compared with previous MHD simulations. In large 'circ; cases, nonlinear simulations show multiple islands forming, followed by the islands coalescing at later times. The competition between these two processes strongly influences the reconnection rates and eventually leads to a steady state reconnection. We will present various parameter studies and show that our hybrid results agree with fluid analysis in certain limits (e.g., relatively large resisitivities).
Commercialization of MHD power technology
Aleman, D.J.; Jensen, A.D.; Probert, P.B.
1984-08-01
This paper presents an approach to the commercialization of Magnetohydrodynamics (MHD) technology from the perspective of an equipment manufacturer. It discusses and recommends actions to be taken in solving technical problems and mitigating risk for the first commercial MHD power plant.
2-D Magnetohydrodynamic Modeling of A Pulsed Plasma Thruster
NASA Technical Reports Server (NTRS)
Thio, Y. C. Francis; Cassibry, J. T.; Wu, S. T.; Rodgers, Stephen L. (Technical Monitor)
2002-01-01
Experiments are being performed on the NASA Marshall Space Flight Center (MSFC) MK-1 pulsed plasma thruster. Data produced from the experiments provide an opportunity to further understand the plasma dynamics in these thrusters via detailed computational modeling. The detailed and accurate understanding of the plasma dynamics in these devices holds the key towards extending their capabilities in a number of applications, including their applications as high power (greater than 1 MW) thrusters, and their use for producing high-velocity, uniform plasma jets for experimental purposes. For this study, the 2-D MHD modeling code, MACH2, is used to provide detailed interpretation of the experimental data. At the same time, a 0-D physics model of the plasma initial phase is developed to guide our 2-D modeling studies.
Fully implicit adaptive mesh refinement MHD algorithm
NASA Astrophysics Data System (ADS)
Philip, Bobby
2005-10-01
In the macroscopic simulation of plasmas, the numerical modeler is faced with the challenge of dealing with multiple time and length scales. The former results in stiffness due to the presence of very fast waves. The latter requires one to resolve the localized features that the system develops. Traditional approaches based on explicit time integration techniques and fixed meshes are not suitable for this challenge, as such approaches prevent the modeler from using realistic plasma parameters to keep the computation feasible. We propose here a novel approach, based on implicit methods and structured adaptive mesh refinement (SAMR). Our emphasis is on both accuracy and scalability with the number of degrees of freedom. To our knowledge, a scalable, fully implicit AMR algorithm has not been accomplished before for MHD. As a proof-of-principle, we focus on the reduced resistive MHD model as a basic MHD model paradigm, which is truly multiscale. The approach taken here is to adapt mature physics-based technologyootnotetextL. Chac'on et al., J. Comput. Phys. 178 (1), 15- 36 (2002) to AMR grids, and employ AMR-aware multilevel techniques (such as fast adaptive composite --FAC-- algorithms) for scalability. We will demonstrate that the concept is indeed feasible, featuring optimal scalability under grid refinement. Results of fully-implicit, dynamically-adaptive AMR simulations will be presented on a variety of problems.
Stability of a two-volume MRxMHD model in slab geometry
NASA Astrophysics Data System (ADS)
Tuen, Li Huey; Hole, Matthew; Dewar, Robert; Dennis, Graham
2015-11-01
Ideal MHD models are inadequate to describe various physical attributes of toroidal plasmas with no continuous symmetry, such as magnetic islands and stochastic regions. A new MHD formulation, MRxMHD, uses only a finite number of ideal-MHD flux surface interfaces, with relaxed plasma regions in between, thus allowing a stepped-pressure profile, magnetic islands, and stochastic regions. In toroidally asymmetric plasma, the existence of interfaces in MRxMHD is contingent on the field pitch, or rotational transform, of flux surfaces being irrational; a KAM argument shows that some good equilibrium flux surfaces continue to exist for small perturbations to an integrable system (foliated by flux surfaces), provided that the rotational transforms on both sides of each interface are sufficiently irrational. Building upon the MRxMHD stability model by Hole, we study the effects of irrationality of the rotational transform at interfaces in MRxMHD on plasma stability. Investigating the plasma stability of a two-volume MRxMHD periodic-slab, we find that the 2D system stability conditions are dependent on the interface and resonance layer magnetic field pitch at minimised energy states. While tearing instabilities exist at low order rational resonances, investigating instability of high-order rationals requires study of pressure-driven instabilities.
NASA Astrophysics Data System (ADS)
Schüssler, M.
1999-05-01
Two aspects of solar MHD are discussed in relation to the work of the MHD simulation group at KIS. Photospheric magneto-convection, the nonlinear interaction of magnetic field and convection in a strongly stratified, radiating fluid, is a key process of general astrophysical relevance. Comprehensive numerical simulations including radiative transfer have significantly improved our understanding of the processes and have become an important tool for the interpretation of observational data. Examples of field intensification in the solar photosphere ('convective collapse') are shown. The second line of research is concerned with the dynamics of flux tubes in the convection zone, which has far-reaching implications for our understanding of the solar dynamo. Simulations indicate that the field strength in the region where the flux is stored before erupting to form sunspot groups is of the order of 105 G, an order of magnitude larger than previous estimates based on equipartition with the kinetic energy of convective flows.Key words. Solar physics · astrophysics and astronomy (photosphere and chromosphere; stellar interiors and dynamo theory; numerical simulation studies).
NASA Technical Reports Server (NTRS)
Wegmann, R.; Schmidt, H. U.; Huebner, W. F.; Boice, D. C.
1987-01-01
An MHD and chemical comet-coma model was developed, applying the computer program of Huebner (1985) for the detailed chemical evolution of a spherically expanding coma and the program of Schmidt and Wegman (1982) and Wegman (1987) for the MHD flow of plasma and magnetic field in a comet to the Giotto-mission data on the ion abundances measured by the HIS ion mass spectrometer. The physics and chemistry of the coma are modeled in great detail, including photoprocesses, gas-phase chemical kinetics, energy balance with a separate electron temperature, multifluid hydrodynamics with a transition to free molecular flow, fast-streaming atomic and molecular hydrogen, counter and cross streaming of the ionized species relative to the neutral species in the coma-solar wind interaction region with momentum exchange by elastic collisions, mass-loading through ion pick-up, and Lorentz forces of the advected magnetic field. The results, both inside and outside of the contact surface, are discussed and compared with the relevant HIS ion mass spectra.
Energy Science and Technology Software Center (ESTSC)
2004-08-01
AnisWave2D is a 2D finite-difference code for a simulating seismic wave propagation in fully anisotropic materials. The code is implemented to run in parallel over multiple processors and is fully portable. A mesh refinement algorithm has been utilized to allow the grid-spacing to be tailored to the velocity model, avoiding the over-sampling of high-velocity materials that usually occurs in fixed-grid schemes.
Energy storage and dissipation in the magnetotail during substorms 2. MHD simulations
NASA Astrophysics Data System (ADS)
Steinolfson, R. S.; Winglee, R. M.
1993-05-01
The effects of temporal and spatial variations in the plasma resistivity on the evolution of the magnetosphere during substorms are examined with numerical solutions of the two-dimensional magnetohydrodynamic (MHD) equations. The global MHD simulations self-consistently consider the interaction of the solar wind with the dayside magnetosphere as well as the evolution of the tail region. These solutions are used to study how various solar wind states generate conditions in the tail, such as pressure gradients and cross-tail currents, that have the potential of leading to a substorm. Although the MHD formalism does provide information on the large-scale evolution, the essential mechanism for substorm development may involve microscopic or particle processes not present in an MHD approach. As a result, this MHD study is carried out in association with particle simulations (Winglee and Steinolfson, this issue). Since one connection between the MHD and particle approaches is through the resistivity, the effects of various resistivity distributions on the global MHD configuration are examined. The resistivity distributions considered here are (1) a temporally constant and spatially uniform resistivity, (2) a resistivity proportional to the square of the local current density, and (3) a resistivity proportional to the square of the local magnetic field strength. The latter distribution is suggested by the above particle simulations and represents effects produced by the increased magnetization of particles and the differential motion between electrons and ions. For all three cases a plasmoid is formed and ejected tailward. However, when the resistivity depends on the field strength, considerably more energy is stored in the tail prior to plasmoid formation, and plasmoid formation is delayed relative to the results for the other two resistivity distributions. Furthermore, when the plasmoid is eventually ejected, it moves down the tail with a higher speed. The MHD results
Numerical linearized MHD model of flapping oscillations
NASA Astrophysics Data System (ADS)
Korovinskiy, D. B.; Ivanov, I. B.; Semenov, V. S.; Erkaev, N. V.; Kiehas, S. A.
2016-06-01
Kink-like magnetotail flapping oscillations in a Harris-like current sheet with earthward growing normal magnetic field component Bz are studied by means of time-dependent 2D linearized MHD numerical simulations. The dispersion relation and two-dimensional eigenfunctions are obtained. The results are compared with analytical estimates of the double-gradient model, which are found to be reliable for configurations with small Bz up to values ˜ 0.05 of the lobe magnetic field. Coupled with previous results, present simulations confirm that the earthward/tailward growth direction of the Bz component acts as a switch between stable/unstable regimes of the flapping mode, while the mode dispersion curve is the same in both cases. It is confirmed that flapping oscillations may be triggered by a simple Gaussian initial perturbation of the Vz velocity.
Observational Tests of Recent MHD Turbulence Perspectives
NASA Astrophysics Data System (ADS)
Ghosh, Sanjoy
2001-06-01
This grant seeks to analyze the Heliospheric Missions data to test current theories on the angular dependence (with respect to mean magnetic field direction) of magnetohydrodynamic (MHD) turbulence in the solar wind. Solar wind turbulence may be composed of two or more dynamically independent components. Such components include magnetic pressure-balanced structures, velocity shears, quasi-2D turbulence, and slab (Alfven) waves. We use a method, developed during the first two years of this grant, for extracting the individual reduced spectra of up to three separate turbulence components from a single spacecraft time series. The method has been used on ISEE-3 data, Pioneer Venus Orbiter, Ulysses, and Voyager data samples. The correlation of fluctuations as a function of angle between flow direction and magnetic-field direction is the focus of study during the third year.
Observational Tests of Recent MHD Turbulence Perspectives
NASA Technical Reports Server (NTRS)
Ghosh, Sanjoy; Guhathakurta, M. (Technical Monitor)
2001-01-01
This grant seeks to analyze the Heliospheric Missions data to test current theories on the angular dependence (with respect to mean magnetic field direction) of magnetohydrodynamic (MHD) turbulence in the solar wind. Solar wind turbulence may be composed of two or more dynamically independent components. Such components include magnetic pressure-balanced structures, velocity shears, quasi-2D turbulence, and slab (Alfven) waves. We use a method, developed during the first two years of this grant, for extracting the individual reduced spectra of up to three separate turbulence components from a single spacecraft time series. The method has been used on ISEE-3 data, Pioneer Venus Orbiter, Ulysses, and Voyager data samples. The correlation of fluctuations as a function of angle between flow direction and magnetic-field direction is the focus of study during the third year.
MHD Simulation Heliospheric Magnetic Fields and Turbulence
NASA Technical Reports Server (NTRS)
Roberts, D. Aaron
2005-01-01
This talk will present a summary of our results on simulations of heliospheric structure and dynamics. We use a three-dimensional MHD code in spherical coordinates to produce a solar wind containing a rotating, tilted heliospheric current sheet, fast-slow stream and microstream shear layers, waves, 2-D turbulence, and pressure balanced structures that are input to the inner (superAlfvenic) boundary. The evolution of various combinations of these has led to a deeper understanding of sector structure, magnetic holes, fluctuation anisotropies, and general turbulent evolution. We show how the sectors are likely to be connected, how spiral fields can arise, and how field line diffusion can be caused by waves with transverse structure and microstream shears.
Not Available
1987-06-01
This report presents a study of the nuclear weapons magnetohydrodynamic (MHD) effects on submarine communications cables. The study consisted of the analysis and interpretation of currently available data on submarine cable systems TAT-4, TAT-6, and TAT-7. The primary result of the study is that decrease of the effective resistivity with frequency over the available experimental range, coupled with the model results, leads to quite small effective resistivities at the MHD characteristic frequencies, and hence small earth potential differences. Thus, it appears that submarine cable systems are less susceptible to an MHD threat than their land-based counter-parts.
MHD Energy Bypass Scramjet Engine
NASA Technical Reports Server (NTRS)
Mehta, Unmeel B.; Bogdanoff, David W.; Park, Chul; Arnold, Jim (Technical Monitor)
2001-01-01
Revolutionary rather than evolutionary changes in propulsion systems are most likely to decrease cost of space transportation and to provide a global range capability. Hypersonic air-breathing propulsion is a revolutionary propulsion system. The performance of scramjet engines can be improved by the AJAX energy management concept. A magneto-hydro-dynamics (MHD) generator controls the flow and extracts flow energy in the engine inlet and a MHD accelerator downstream of the combustor accelerates the nozzle flow. A progress report toward developing the MHD technology is presented herein. Recent theoretical efforts are reviewed and ongoing experimental efforts are discussed. The latter efforts also include an ongoing collaboration between NASA, the US Air Force Research Laboratory, US industry, and Russian scientific organizations. Two of the critical technologies, the ionization of the air and the MHD accelerator, are briefly discussed. Examples of limiting the combustor entrance Mach number to a low supersonic value with a MHD energy bypass scheme are presented, demonstrating an improvement in scramjet performance. The results for a simplified design of an aerospace plane show that the specific impulse of the MHD-bypass system is better than the non-MHD system and typical rocket over a narrow region of flight speeds and design parameters. Equilibrium ionization and non-equilibrium ionization are discussed. The thermodynamic condition of air at the entrance of the engine inlet determines the method of ionization. The required external power for non-equilibrium ionization is computed. There have been many experiments in which electrical power generation has successfully been achieved by magneto-hydrodynamic (MHD) means. However, relatively few experiments have been made to date for the reverse case of achieving gas acceleration by the MHD means. An experiment in a shock tunnel is described in which MHD acceleration is investigated experimentally. MHD has several
NASA Astrophysics Data System (ADS)
Sych, Robert
2016-02-01
The study of magnetohydrodynamic (MHD) waves and oscillations in the solar atmosphere is one of the fastest developing fields in solar physics, and lies in the mainstream of using solar instrumentation data. This chapter first addresses the spatial frequency morphology of sources of sunspot oscillations and waves, including their localization, size, oscillation periods, and height localization with the mechanism of cutoff frequency that forms the observed emission variability. Then, it presents a review dynamic of sunspot wave processes, provides the information about the structure of wave fronts and their time variations, and investigates the oscillation frequency transformation depending on the wave energy. The chapter also addresses the initializing solar flares caused by trigger agents like magnetoacoustic waves, accelerated particle beams, and shocks. Special attention is paid to the relation between the flare reconnection periodic initialization and the dynamics of sunspot slow magnetoacoustic waves.
Lacey, James J.; Kurtzrock, Roy C.; Bienstock, Daniel
1976-08-24
A hot gaseous fluid of low ash content, suitable for use in open-cycle MHD (magnetohydrodynamic) power generation, is produced by means of a three-stage process comprising (1) partial combustion of a fossil fuel to produce a hot gaseous product comprising CO.sub.2 CO, and H.sub.2 O, (2) reformation of the gaseous product from stage (1) by means of a fluidized char bed, whereby CO.sub.2 and H.sub.2 O are converted to CO and H.sub.2, and (3) combustion of CO and H.sub.2 from stage (2) to produce a low ash-content fluid (flue gas) comprising CO.sub.2 and H.sub.2 O and having a temperature of about 4000.degree. to 5000.degree.F.
MHD channel performance for potential early commercial MHD power plants
NASA Technical Reports Server (NTRS)
Swallom, D. W.
1981-01-01
The commercial viability of full and part load early commercial MHD power plants is examined. The load conditions comprise a mass flow of 472 kg/sec in the channel, Rosebud coal, 34% by volume oxygen in the oxidizer preheated to 922 K, and a one percent by mass seeding with K. The full load condition is discussed in terms of a combined cycle plant with optimized electrical output by the MHD channel. Various electrical load parameters, pressure ratios, and magnetic field profiles are considered for a baseload MHD generator, with a finding that a decelerating flow rate yields slightly higher electrical output than a constant flow rate. Nominal and part load conditions are explored, with a reduced gas mass flow rate and an enriched oxygen content. An enthalpy extraction of 24.6% and an isentropic efficiency of 74.2% is predicted for nominal operation of a 526 MWe MHD generator, with higher efficiencies for part load operation.
MHD properties of magnetosheath flow
NASA Astrophysics Data System (ADS)
Siscoe, G. L.; Crooker, N. U.; Erickson, G. M.; Sonnerup, B. U. Ö.; Maynard, N. C.; Schoendorf, J. A.; Siebert, K. D.; Weimer, D. R.; White, W. W.; Wilson, G. R.
2002-04-01
We discuss four aspects of magnetosheath flow that require MHD for their calculation and understanding. We illustrate these aspects with computations using a numerical MHD code that simulates the global magnetosphere and its magnetosheath. The four inherently MHD aspects of magnetosheath flow that we consider are the depletion layer, the magnetospheric sash, MHD flow deflections, and the magnetosheath's slow-mode expansion into the magnetotail. We introduce new details of these aspects or illustrate known details in a new way, including the dependence of the depletion layer on interplanetary magnetic filed clock angle; agreement between the locations of the antiparallel regions of Luhmann et al. (J. Geophys. Res. 89 (1984) 1739) and the magnetospheric sash, and deflections corresponding separately to a stagnation line and magnetic reconnection.
MHD technology in aluminum casting
Kalinichenko, I.
1984-08-01
The use of MHD technology in aluminum casting is discussed. Associates of the Latvian Academy of Sciences Institute of Physics developed magnetohydrodynamic units for the Siberian plant. A MHD unit made it possible to free five persons from heavy work at the plant. Labor productivity doubled in this section. With the aid of the magnetic field, the alloy silumin is obtained in only three hours. Specialists of the Irkutsk affiliate of the All-Union Scientific Research and Design Institute of the Aluminum, Magnesium and Electrode Industry are convinced that MHD technology has a bright future. However, this will necessitate the development of new MHD technology for different types of casting facilities, with their specific features taken into account.
Propagation and Dissipation of MHD Waves in Coronal Holes
NASA Astrophysics Data System (ADS)
Dwivedi, B. N.
2006-11-01
bholadwivedi@gmail.com In view of the landmark result on the solar wind outflow, starting between 5 Mm and 20 Mm above the photosphere in magnetic funnels, we investigate the propagation and dissipation of MHD waves in coronal holes. We underline the importance of Alfvén wave dissipation in the magnetic funnels through the viscous and resistive plasma. Our results show that Alfvén waves are one of the primary energy sources in the innermost part of coronal holes where the solar wind outflow starts. We also consider compressive viscosity and thermal conductivity to study the propagation and dissipation of long period slow longitudinal MHD waves in polar coronal holes. We discuss their likely role in the line profile narrowing, and in the energy budget for coronal holes and the solar wind. We compare the contribution of longitudinal MHD waves with high frequency Alfvén waves.
MHD Simulations of Thermal Plasma Jets in Coaxial Plasma Accelerators
NASA Astrophysics Data System (ADS)
Subramaniam, Vivek; Raja, Laxminarayan
2015-09-01
The development of a magneto-hydrodynamics (MHD) numerical tool to study high energy density thermal plasma in coaxial plasma accelerators is presented. The coaxial plasma accelerator is a device used simulate the conditions created at the confining wall of a thermonuclear fusion reactor during an edge localized mode (ELM) disruption event. This is achieved by creating magnetized thermal plasma in a coaxial volume which is then accelerated by the Lorentz force to form a high velocity plasma jet. The simulation tool developed solves the resistive MHD equation using a finite volume method (FVM) framework. The acceleration and subsequent demagnetization of the plasma as it travels down the length of the accelerator is simulated and shows good agreement with experiments. Additionally, a model to study the thermalization of the plasma at the inlet is being developed in order to give self-consistent initial conditions to the MHD solver.
Perturbed Stability Analysis of External Ideal MHD Modes
NASA Astrophysics Data System (ADS)
Comer, K. J.; Callen, J. D.; Hegna, C. C.; Garstka, G. D.; Turnbull, A. D.; Garofalo, A. M.; Cowley, S. C.
2002-11-01
Traditionally, numerical parameter scans are performed to study the effects of equilibrium shaping and profiles on long wavelength ideal MHD instabilities. Previously, we introduced a new perturbative technique to more efficiently explore these dependencies: changes in delta-W due to small equilibrium variations are found using a perturbation of the energy principle rather than with an eigenvalue-solver instability code. With this approach, the stability properties of similar equilibria can be efficiently explored without generating complete numerical results for every set of parameters (which is time-intensive for accurate representations of several configurations). Here, we apply this approach to toroidal geometry using GATO (an ideal MHD stability code) and experimental equilibria. In particular, we explore ideal MHD stability of external kink modes in the spherical tokamak Pegasus and resistive wall modes in DIII-D.
NASA Astrophysics Data System (ADS)
Fuentes-Fernández, J.; Parnell, C. E.; Hood, A. W.
2010-05-01
Context. For the last thirty years, most of the studies on the relaxation of stressed magnetic fields in the solar environment have only considered the Lorentz force, neglecting plasma contributions, and therefore, limiting every equilibrium to that of a force-free field. Aims: Here we begin a study of the non-resistive evolution of finite beta plasmas and their relaxation to magnetohydrostatic states, where magnetic forces are balanced by plasma-pressure gradients, by using a simple 2D scenario involving a hydromagnetic disturbance to a uniform magnetic field. The final equilibrium state is predicted as a function of the initial disturbances, with aims to demonstrate what happens to the plasma during the relaxation process and to see what effects it has on the final equilibrium state. Methods: A set of numerical experiments are run using a full MHD code, with the relaxation driven by magnetoacoustic waves damped by viscous effects. The numerical results are compared with analytical calculations made within the linear regime, in which the whole process must remain adiabatic. Particular attention is paid to the thermodynamic behaviour of the plasma during the relaxation. Results: The analytical predictions for the final non force-free equilibrium depend only on the initial perturbations and the total pressure of the system. It is found that these predictions hold surprisingly well even for amplitudes of the perturbation far outside the linear regime. Conclusions: Including the effects of a finite plasma beta in relaxation experiments leads to significant differences from the force-free case.
Mathematical modelling in MHD technology
Scheindlin, A.E.; Medin, S.A. )
1990-01-01
The technological scheme and the general parameters of the commercial scale pilot MHD power plant are described. The characteristics of the flow train components and the electrical equipment are discussed. The basic ideas of the mathematical modelling of the processes and the devices operation in MHD systems are considered. The application of different description levels in computer simulation is analyzed and the examples of typical solutions are presented.
NASA Astrophysics Data System (ADS)
Mayor, Louise
2016-05-01
Graphene might be the most famous example, but there are other 2D materials and compounds too. Louise Mayor explains how these atomically thin sheets can be layered together to create flexible “van der Waals heterostructures”, which could lead to a range of novel applications.
Dynamo action in dissipative, forced, rotating MHD turbulence
NASA Astrophysics Data System (ADS)
Shebalin, John V.
2016-06-01
Magnetohydrodynamic (MHD) turbulence is an inherent feature of large-scale, energetic astrophysical and geophysical magnetofluids. In general, these are rotating and are energized through buoyancy and shear, while viscosity and resistivity provide a means of dissipation of kinetic and magnetic energy. Studies of unforced, rotating, ideal (i.e., non-dissipative) MHD turbulence have produced interesting results, but it is important to determine how these results are affected by dissipation and forcing. Here, we extend our previous work and examine dissipative, forced, and rotating MHD turbulence. Incompressibility is assumed, and finite Fourier series represent turbulent velocity and magnetic field on a 643 grid. Forcing occurs at an intermediate wave number by a method that keeps total energy relatively constant and allows for injection of kinetic and magnetic helicity. We find that 3-D energy spectra are asymmetric when forcing is present. We also find that dynamo action occurs when forcing has either kinetic or magnetic helicity, with magnetic helicity injection being more important. In forced, dissipative MHD turbulence, the dynamo manifests itself as a large-scale coherent structure that is similar to that seen in the ideal case. These results imply that MHD turbulence, per se, may play a fundamental role in the creation and maintenance of large-scale (i.e., dipolar) stellar and planetary magnetic fields.
Extended MHD Stabiliy Calculations of Spheromak Equilibria
NASA Astrophysics Data System (ADS)
Howell, E. C.; Sovinec, C. R.
2013-10-01
Linear extended MHD calculations of spheromak equilibria in a cylindrical flux conserver are performed using the NIMROD code (Sovinec et al., JCP 195, 2004). A series of Grad-Sharfranov equilibria are generated with β ranging from 0 . 4 % to 4 . 2 % , corresponding to peak electron temperatures ranging 50 to 300 eV. These equilibria use a λ profile representative of SSPX shot 14590, which measured a peak electron temperature of 325 eV (McLean et al., POP 13, 2006). Resistive MHD calculations find that the β = 0 . 4 % case is unstable to resonant resistive interchange modes with γτA <= 2 . 3 % . These modes transition to ideal interchange as the equilibrium pressure is increased. Growth rates as large as γτA = 20 % are calculated for the 4 . 2 % β case. Calculations including ion-gyroviscosity show a minimal reduction of growth rate. Effects from including the Hall and Electron pressure terms in Ohm's Law and the cross-field diamagnetic heat flux are investigated. Results of related nonlinear simulations are also presented. Work Supported by US DOE.
Graphene suspensions for 2D printing
NASA Astrophysics Data System (ADS)
Soots, R. A.; Yakimchuk, E. A.; Nebogatikova, N. A.; Kotin, I. A.; Antonova, I. V.
2016-04-01
It is shown that, by processing a graphite suspension in ethanol or water by ultrasound and centrifuging, it is possible to obtain particles with thicknesses within 1-6 nm and, in the most interesting cases, 1-1.5 nm. Analogous treatment of a graphite suspension in organic solvent yields eventually thicker particles (up to 6-10 nm thick) even upon long-term treatment. Using the proposed ink based on graphene and aqueous ethanol with ethylcellulose and terpineol additives for 2D printing, thin (~5 nm thick) films with sheet resistance upon annealing ~30 MΩ/□ were obtained. With the ink based on aqueous graphene suspension, the sheet resistance was ~5-12 kΩ/□ for 6- to 15-nm-thick layers with a carrier mobility of ~30-50 cm2/(V s).
The Biermann catastrophe of numerical MHD
NASA Astrophysics Data System (ADS)
Graziani, C.; Tzeferacos, P.; Lee, D.; Lamb, D. Q.; Weide, K.; Fatenejad, M.; Miller, J.
2016-05-01
The Biermann Battery effect is frequently invoked in cosmic magnetogenesis and studied in High-Energy Density laboratory physics experiments. Unfortunately, direct implementation of the Biermann effect in MHD codes is known to produce unphysical magnetic fields at shocks whose value does not converge with resolution. We show that this convergence breakdown is due to naive discretization, which fails to account for the fact that discretized irrotational vector fields have spurious solenoidal components that grow without bound near a discontinuity. We show that careful consideration of the kinetics of ion viscous shocks leads to a formulation of the Biermann effect that gives rise to a convergent algorithm. We note a novel physical effect a resistive magnetic precursor in which Biermann-generated field in the shock “leaks” resistively upstream. The effect appears to be potentially observable in experiments at laser facilities.
Energy Science and Technology Software Center (ESTSC)
2001-01-31
This software reduces the data from two-dimensional kSA MOS program, k-Space Associates, Ann Arbor, MI. Initial MOS data is recorded without headers in 38 columns, with one row of data per acquisition per lase beam tracked. The final MOSS 2d data file is reduced, graphed, and saved in a tab-delimited column format with headers that can be plotted in any graphing software.
Spatially Resolved Synthetic Spectra from 2D Simulations of Stainless Steel Wire Array Implosions
Clark, R. W.; Giuliani, J. L.; Thornhill, J. W.; Chong, Y. K.; Dasgupta, A.; Davis, J.
2009-01-21
A 2D radiation MHD model has been developed to investigate stainless steel wire array implosion experiments on the Z and refurbished Z machines. This model incorporates within the Mach2 MHD code a self-consistent calculation of the non-LTE kinetics and ray trace based radiation transport. Such a method is necessary in order to account for opacity effects in conjunction with ionization kinetics of K-shell emitting plasmas. Here the model is used to investigate multi-dimensional effects of stainless steel wire implosions. In particular, we are developing techniques to produce non-LTE, axially and/or radially resolved synthetic spectra based upon snapshots of our 2D simulations. Comparisons between experimental spectra and these synthetic spectra will allow us to better determine the state of the experimental pinches.
An advanced implicit solver for MHD
NASA Astrophysics Data System (ADS)
Udrea, Bogdan
A new implicit algorithm has been developed for the solution of the time-dependent, viscous and resistive single fluid magnetohydrodynamic (MHD) equations. The algorithm is based on an approximate Riemann solver for the hyperbolic fluxes and central differencing applied on a staggered grid for the parabolic fluxes. The algorithm employs a locally aligned coordinate system that allows the solution to the Riemann problems to be solved in a natural direction, normal to cell interfaces. The result is an original scheme that is robust and reduces the complexity of the flux formulas. The evaluation of the parabolic fluxes is also implemented using a locally aligned coordinate system, this time on the staggered grid. The implicit formulation employed by WARP3 is a two level scheme that was applied for the first time to the single fluid MHD model. The flux Jacobians that appear in the implicit scheme are evaluated numerically. The linear system that results from the implicit discretization is solved using a robust symmetric Gauss-Seidel method. The code has an explicit mode capability so that implementation and test of new algorithms or new physics can be performed in this simpler mode. Last but not least the code was designed and written to run on parallel computers so that complex, high resolution runs can be per formed in hours rather than days. The code has been benchmarked against analytical and experimental gas dynamics and MHD results. The benchmarks consisted of one-dimensional Riemann problems and diffusion dominated problems, two-dimensional supersonic flow over a wedge, axisymmetric magnetoplasmadynamic (MPD) thruster simulation and three-dimensional supersonic flow over intersecting wedges and spheromak stability simulation. The code has been proven to be robust and the results of the simulations showed excellent agreement with analytical and experimental results. Parallel performance studies showed that the code performs as expected when run on parallel
Nonlinear Alfvén wave dynamics at a 2D magnetic null point: ponderomotive force
NASA Astrophysics Data System (ADS)
Thurgood, J. O.; McLaughlin, J. A.
2013-07-01
Context. In the linear, β = 0 MHD regime, the transient properties of magnetohydrodynamic (MHD) waves in the vicinity of 2D null points are well known. The waves are decoupled and accumulate at predictable parts of the magnetic topology: fast waves accumulate at the null point; whereas Alfvén waves cannot cross the separatricies. However, in nonlinear MHD mode conversion can occur at regions of inhomogeneous Alfvén speed, suggesting that the decoupled nature of waves may not extend to the nonlinear regime. Aims: We investigate the behaviour of low-amplitude Alfvén waves about a 2D magnetic null point in nonlinear, β = 0 MHD. Methods: We numerically simulate the introduction of low-amplitude Alfvén waves into the vicinity of a magnetic null point using the nonlinear LARE2D code. Results: Unlike in the linear regime, we find that the Alfvén wave sustains cospatial daughter disturbances, manifest in the transverse and longitudinal fluid velocity, owing to the action of nonlinear magnetic pressure gradients (viz. the ponderomotive force). These disturbances are dependent on the Alfvén wave and do not interact with the medium to excite magnetoacoustic waves, although the transverse daughter becomes focused at the null point. Additionally, an independently propagating fast magnetoacoustic wave is generated during the early stages, which transports some of the initial Alfvén wave energy towards the null point. Subsequently, despite undergoing dispersion and phase-mixing due to gradients in the Alfvén-speed profile (∇cA ≠ 0) there is no further nonlinear generation of fast waves. Conclusions: We find that Alfvén waves at 2D cold null points behave largely as in the linear regime, however they sustain transverse and longitudinal disturbances - effects absent in the linear regime - due to nonlinear magnetic pressure gradients.
2D superconductivity by ionic gating
NASA Astrophysics Data System (ADS)
Iwasa, Yoshi
2D superconductivity is attracting a renewed interest due to the discoveries of new highly crystalline 2D superconductors in the past decade. Superconductivity at the oxide interfaces triggered by LaAlO3/SrTiO3 has become one of the promising routes for creation of new 2D superconductors. Also, the MBE grown metallic monolayers including FeSe are also offering a new platform of 2D superconductors. In the last two years, there appear a variety of monolayer/bilayer superconductors fabricated by CVD or mechanical exfoliation. Among these, electric field induced superconductivity by electric double layer transistor (EDLT) is a unique platform of 2D superconductivity, because of its ability of high density charge accumulation, and also because of the versatility in terms of materials, stemming from oxides to organics and layered chalcogenides. In this presentation, the following issues of electric filed induced superconductivity will be addressed; (1) Tunable carrier density, (2) Weak pinning, (3) Absence of inversion symmetry. (1) Since the sheet carrier density is quasi-continuously tunable from 0 to the order of 1014 cm-2, one is able to establish an electronic phase diagram of superconductivity, which will be compared with that of bulk superconductors. (2) The thickness of superconductivity can be estimated as 2 - 10 nm, dependent on materials, and is much smaller than the in-plane coherence length. Such a thin but low resistance at normal state results in extremely weak pinning beyond the dirty Boson model in the amorphous metallic films. (3) Due to the electric filed, the inversion symmetry is inherently broken in EDLT. This feature appears in the enhancement of Pauli limit of the upper critical field for the in-plane magnetic fields. In transition metal dichalcogenide with a substantial spin-orbit interactions, we were able to confirm the stabilization of Cooper pair due to its spin-valley locking. This work has been supported by Grant-in-Aid for Specially
MHD Integrated Topping Cycle Project
Not Available
1992-03-01
The Magnetohydrodynamics (MHD) Integrated Topping Cycle (ITC) Project represents the culmination of the proof-of-concept (POC) development stage in the US Department of Energy (DOE) program to advance MHD technology to early commercial development stage utility power applications. The project is a joint effort, combining the skills of three topping cycle component developers: TRW, Avco/TDS, and Westinghouse. TRW, the prime contractor and system integrator, is responsible for the 50 thermal megawatt (50 MW{sub t}) slagging coal combustion subsystem. Avco/TDS is responsible for the MHD channel subsystem (nozzle, channel, diffuser, and power conditioning circuits), and Westinghouse is responsible for the current consolidation subsystem. The ITC Project will advance the state-of-the-art in MHD power systems with the design, construction, and integrated testing of 50 MW{sub t} power train components which are prototypical of the equipment that will be used in an early commercial scale MHD utility retrofit. Long duration testing of the integrated power train at the Component Development and Integration Facility (CDIF) in Butte, Montana will be performed, so that by the early 1990's, an engineering data base on the reliability, availability, maintainability and performance of the system will be available to allow scaleup of the prototypical designs to the next development level. This Sixteenth Quarterly Technical Progress Report covers the period May 1, 1991 to July 31, 1991.
Nanoimprint lithography: 2D or not 2D? A review
NASA Astrophysics Data System (ADS)
Schift, Helmut
2015-11-01
Nanoimprint lithography (NIL) is more than a planar high-end technology for the patterning of wafer-like substrates. It is essentially a 3D process, because it replicates various stamp topographies by 3D displacement of material and takes advantage of the bending of stamps while the mold cavities are filled. But at the same time, it keeps all assets of a 2D technique being able to pattern thin masking layers like in photon- and electron-based traditional lithography. This review reports about 20 years of development of replication techniques at Paul Scherrer Institut, with a focus on 3D aspects of molding, which enable NIL to stay 2D, but at the same time enable 3D applications which are "more than Moore." As an example, the manufacturing of a demonstrator for backlighting applications based on thermally activated selective topography equilibration will be presented. This technique allows generating almost arbitrary sloped, convex and concave profiles in the same polymer film with dimensions in micro- and nanometer scale.
Equilibrium and MHD Activity in PEGASUS Plasmas
NASA Astrophysics Data System (ADS)
Garstka, G. D.
2000-10-01
The primary research goals of the PEGASUS toroidal experiment involve the exploration of MHD stability boundaries at high beta and extremely low aspect ratio. A novel nonlinear least-squares fitting technique is used to reconstruct equilibria; the more conventional TokaMac equilibrium code is used as well. The continuous resistive vacuum vessel is modeled as a set of inductively coupled coils for the purpose of the reconstructions. A variety of MHD phenomena have been observed on PEGASUS ohmic discharges. Internal reconnection events (IREs) are often seen. The characteristics of these events match those observed on other STs: they are observed more frequently when the machine is dirty and when the plasma is overdriven into the central column. An n=1 mode with a frequency of 3-8 kHz is often present throughout the discharge. This mode is associated with fast current ramps (> 30 MA/s) and can limit the discharge evolution. There is also substantial evidence of double tearing modes during fast current ramps. Present work involves the exploration of the edge kink stability boundary at near-unity aspect ratio.
MHD stability of tokamak plasmas
Chance, M.S. Sun, Y.C.; Jardin, S.C.; Kessel, C.E.; Okabayashi, M.
1992-08-01
This paper will give an overview of the some of the methods which are used to simulate the ideal MHD properties of tokamak plasmas. A great deal of the research in this field is necessarily numerical and the substantial progress made during the past several years has roughly paralleled the continuing availability of more advanced supercomputers. These have become essential to accurately model the complex configurations necessary for achieving MHD stable reactor grade conditions. Appropriate tokamak MHD equilibria will be described. Then the stability properties is discussed in some detail, emphasizing the difficulties of obtaining stable high {beta} discharges in plasmas in which the current is mainly ohmically driven and thus demonstrating the need for tailoring the current and pressure profiles of the plasma away from the ohmic state. The outline of this paper will roughly follow the physics development to attain the second region of stability in the PBX-M device at The Princeton Plasmas Physics Laboratory.
Magnetohydrodynamic (MHD) channel corner seal
Spurrier, Francis R.
1980-01-01
A corner seal for an MHD duct includes a compressible portion which contacts the duct walls and an insulating portion which contacts the electrodes, sidewall bars and insulators. The compressible portion may be a pneumatic or hydraulic gasket or an open-cell foam rubber. The insulating portion is segmented into a plurality of pieces of the same thickness as the electrodes, insulators and sidewall bars and aligned therewith, the pieces aligned with the insulator being of a different size from the pieces aligned with the electrodes and sidewall bars to create a stepped configuration along the corners of the MHD channel.
NASA Astrophysics Data System (ADS)
Yokoyama, Takaaki
Temporal evolution of a current sheet with initial perturbations is studied by using the threedimensional resistive magnetohydrodynamic (MHD) simulations. The magnetic reconnection is considered to be the main engine of the energy rele ase in solar flares. The structure of the diffusion region is, however, not stil l understood under the circumstances with enormously large magnetic Reynolds num ber as the solar corona. In particular, the relationship between the flare's macroscopic physics and the microscopic ones are unclear. It is generally believed that the MHD turbulence s hould play a role in the intermediate scale. The initial current sheet is in an approximately hydromagnetic equilibrium with anti-parallel magnetic field in the y-direction. We imposed a finite-amplitude perturbations (=50ee what happens. Special attention is paid upon the evolution of a three-dimens ional structure in the direction along the initial electric current (z-direction ). Our preliminary results are as follows: (1) In the early phase of the evolut ion, high wavenumber modes in the z-direction are excited and grow. (2) Many "X "-type neutral points (lines) are generated along the magnetic neutral line (pla ne) in the current sheet. When they evolve into the non-linear phase, three-dime nsional structures in the z-direction also evolve. The spatial scale in the z-di rection seems to be almost comparable with that in the xy-plane. (3) The energy release rate is reduced in case of 3D simulations compared with 2D ones probably because of the reduction of the inflow cross sections by the formation of pattc hy structures in the current sheet.
On 2-D Boussinesq equations for MHD convection with stratification effects
NASA Astrophysics Data System (ADS)
Bian, Dongfen; Gui, Guilong
2016-08-01
This paper is concerned with the two-dimensional magnetohydrodynamics-Boussinesq system with the temperature-dependent viscosity, thermal diffusivity and electrical conductivity. The first progress on this topic was made independently by Chae and Hou-Li [8,26] where the Boussinesq system with partial constant viscosity is obtained. Recently, Wang-Zhang [45] considered the temperature-dependent viscosity and thermal diffusivity, and Li-Xu [16] generalized the Wang-Zhang's result to the inviscid case with temperature-dependent thermal diffusivity. In this paper, we include the stratification and magnetic effects and consider the full system, in the framework of low regularity. We prove that, without any smallness assumption on the initial data, the full system is globally well-posed. Moreover, by applying the uniformly bounded generalized Oseen operator, time decay estimate of the solution is obtained.
NASA Astrophysics Data System (ADS)
Stone, James M.; Norman, Michael L.
1992-06-01
In this, the second of a series of three papers, we continue a detailed description of ZEUS-2D, a numerical code for the simulation of fluid dynamical flows in astrophysics including a self-consistent treatment of the effects of magnetic fields and radiation transfer. In this paper, we give a detailed description of the magnetohydrodynamical (MHD) algorithms in ZEUS-2D. The recently developed constrained transport (CT) algorithm is implemented for the numerical evolution of the components of the magnetic field for MHD simulations. This formalism guarantees the numerically evolved field components will satisfy the divergence-free constraint at all times. We find, however, that the method used to compute the electromotive forces must be chosen carefully to propagate accurately all modes of MHD wave families (in particular shear Alfvén waves). A new method of computing the electromotive force is developed using the method of characteristics (MOC). It is demonstrated through the results of an extensive series of MHD test problems that the resulting hybrid MOC-CT method provides for the accurate evolution of all modes of MHD wave families.
Wall surface leakage effects on MHD power generator performance
Pian, C.C.P.; Schmitt, E.W.
1994-12-31
Internal surface leakage effects on the MHD generator performance were studied using a combined experimental and analytical approach. A method to determine the wall resistances and slag layer conductivities from seed shut-off test data is introduced. These measured resistance values are then utilized in generator performance analyses. Calculated results were compared with measured data from MHD generator tests to verify the modeling approach. Finally, these calculated results were used to investigate the distribution of internal leakage currents as a function of generator size, generator operating conditions, and iron oxide injection rates. An advantage of this analysis methodology is the ability to differentiate between wall leakage and apparent leakage effects in the measured test data.
Corrosion and arc erosion in MHD channels
Rosa, R.J. . Dept. of Mechanical Engineering); Pollina, R.J. . Dept. of Mechanical Engineering Avco-Everett Research Lab., Everett, MA )
1991-10-01
The objective of this task is to study the corrosion and arc erosion of MHD materials in a cooperative effort with, and to support, the MHD topping cycle program. Materials tested in the Avco Research Laboratory/Textron facility, or materials which have significant MHD importance, will be analyzed to document their physical deterioration. Conclusions shall be drawn about their wear mechanisms and lifetime in the MHD environment with respect to the following issues; sulfur corrosion, electrochemical corrosion, and arc erosion. The impact of any materials or slag conditions on the level of power output and on the level of leakage current in the MHD channel will also be noted, where appropriate.
Open cycle gas fired MHD power plants
Medin, S.A. ); Negrini, F. )
1991-01-01
In this paper, the main objectives for the present development of gas fired MHD power generation are considered. The state of the world-wide natural gas consumption and its utilization for electricity production is analyzed. The experimental efforts in gas-fired MHD studies are briefly described. The essential features of the two major world gas-fired MHD project - the Ryazan MHDES-580 (U-500) power plant and the Italian 230 MWt retrofit are presented. New suggestions for improving the efficiency of MHD systems and the theoretical and experimental aspects of MHD development are discussed.
MHD integrated topping cycle project
NASA Astrophysics Data System (ADS)
1992-01-01
The overall objective of the project is to design and construct prototypical hardware for an integrated MHD topping cycle, and conduct long duration proof-of-concept tests of integrated system at the US DOE Component Development and Integration Facility in Butte, Montana. The results of the long duration tests will augment the existing engineering design data base on MHD power train reliability, availability, maintainability, and performance, and will serve as a basis for scaling up the topping cycle design to the next level of development, an early commercial scale power plant retrofit. The components of the MHD power train to be designed, fabricated, and tested include: A slagging coal combustor with a rated capacity of 50 MW thermal input, capable of operation with an Eastern (Illinois no. 6) or Western (Montana Rosebud) coal, a segmented supersonic nozzle, a supersonic MHD channel capable of generating at least 1.5 MW of electrical power, a segmented supersonic diffuser section to interface the channel with existing facility quench and exhaust systems, a complete set of current control circuits for local diagonal current control along the channel, and a set of current consolidation circuits to interface the channel with the existing facility inverter.
MHD Integrated Topping Cycle Project
Not Available
1992-01-01
The overall objective of the project is to design and construct prototypical hardware for an integrated MHD topping cycle, and conduct long duration proof-of-concept tests of integrated system at the US DOE Component Development and Integration Facility in Butte, Montana. The results of the long duration tests will augment the existing engineering design data base on MHD power train reliability, availability, maintainability, and performance, and will serve as a basis for scaling up the topping cycle design to the next level of development, an early commercial scale power plant retrofit. The components of the MHD power train to be designed, fabricated, and tested include: A slagging coal combustor with a rated capacity of 50 MW thermal input, capable of operation with an Eastern (Illinois {number sign}6) or Western (Montana Rosebud) coal, a segmented supersonic nozzle, a supersonic MHD channel capable of generating at least 1.5 MW of electrical power, a segmented supersonic diffuser section to interface the channel with existing facility quench and exhaust systems, a complete set of current control circuits for local diagonal current control along the channel, and a set of current consolidation circuits to interface the channel with the existing facility inverter.
Variational Integration for Ideal MHD with Built-in Advection Equations
Zhou, Yao; Qin, Hong; Burby, J. W.; Bhattacharjee, A.
2014-08-05
Newcomb's Lagrangian for ideal MHD in Lagrangian labeling is discretized using discrete exterior calculus. Variational integrators for ideal MHD are derived thereafter. Besides being symplectic and momentum preserving, the schemes inherit built-in advection equations from Newcomb's formulation, and therefore avoid solving them and the accompanying error and dissipation. We implement the method in 2D and show that numerical reconnection does not take place when singular current sheets are present. We then apply it to studying the dynamics of the ideal coalescence instability with multiple islands. The relaxed equilibrium state with embedded current sheets is obtained numerically.
Metrology for graphene and 2D materials
NASA Astrophysics Data System (ADS)
Pollard, Andrew J.
2016-09-01
The application of graphene, a one atom-thick honeycomb lattice of carbon atoms with superlative properties, such as electrical conductivity, thermal conductivity and strength, has already shown that it can be used to benefit metrology itself as a new quantum standard for resistance. However, there are many application areas where graphene and other 2D materials, such as molybdenum disulphide (MoS2) and hexagonal boron nitride (h-BN), may be disruptive, areas such as flexible electronics, nanocomposites, sensing and energy storage. Applying metrology to the area of graphene is now critical to enable the new, emerging global graphene commercial world and bridge the gap between academia and industry. Measurement capabilities and expertise in a wide range of scientific areas are required to address this challenge. The combined and complementary approach of varied characterisation methods for structural, chemical, electrical and other properties, will allow the real-world issues of commercialising graphene and other 2D materials to be addressed. Here, examples of metrology challenges that have been overcome through a multi-technique or new approach are discussed. Firstly, the structural characterisation of defects in both graphene and MoS2 via Raman spectroscopy is described, and how nanoscale mapping of vacancy defects in graphene is also possible using tip-enhanced Raman spectroscopy (TERS). Furthermore, the chemical characterisation and removal of polymer residue on chemical vapour deposition (CVD) grown graphene via secondary ion mass spectrometry (SIMS) is detailed, as well as the chemical characterisation of iron films used to grow large domain single-layer h-BN through CVD growth, revealing how contamination of the substrate itself plays a role in the resulting h-BN layer. In addition, the role of international standardisation in this area is described, outlining the current work ongoing in both the International Organization of Standardization (ISO) and the
Broken Ergodicity in MHD Turbulence
NASA Technical Reports Server (NTRS)
Shebalin, John V.
2010-01-01
Ideal magnetohydrodynamic (MHD) turbulence may be represented by finite Fourier series, where the inherent periodic box serves as a surrogate for a bounded astrophysical plasma. Independent Fourier coefficients form a canonical ensemble described by a Gaussian probability density function containing a Hermitian covariance matrix with positive eigenvalues. The eigenvalues at lowest wave number can be very small, resulting in a large-scale coherent structure: a turbulent dynamo. This is seen in computations and a theoretical explanation in terms of 'broken ergodicity' contains Taylor s theory of force-free states. An important problem for future work is the case of real, i.e., dissipative flows. In real flows, broken ergodicity and coherent structure are still expected to occur in MHD turbulence at the largest scale, as suggested by low resolution simulations. One challenge is to incorporate coherent structure at the largest scale into the theory of turbulent fluctuations at smaller scales.
NASA Technical Reports Server (NTRS)
Wu, S. T.; Song, M. T.; Martens, P. C. H.; Dryer, M.
1991-01-01
A situation wherein a bipolar magnetic field embedded in a stratified solar atmosphere undergoes symmetrical shear motion at the footpoints is investigated via a 2D (nonplanar) MHD simulation. It was found that the vertical plasma flow velocities grow exponentially, leading to a new type of global MHD instability. The growth rate increases almost linearly until it reaches the same order of magnitude as the Alfven speed. Then a nonlinear MHD instability occurs beyond this point. It was found that the central loops are pinched by opposing Lorentz forces, and the outer closed loops stretch upward with the vertically-rising mass flow. The nonlinear dynamical shearing instability is illustrated by a numerical example that is given for three different values of the plasma beta that span several orders of magnitude.
Newtonian CAFE: a new ideal MHD code to study the solar atmosphere
NASA Astrophysics Data System (ADS)
González, J. J.; Guzmán, F.
2015-12-01
In this work we present a new independent code designed to solve the equations of classical ideal magnetohydrodynamics (MHD) in three dimensions, submitted to a constant gravitational field. The purpose of the code centers on the analysis of solar phenomena within the photosphere-corona region. In special the code is capable to simulate the propagation of impulsively generated linear and non-linear MHD waves in the non-isothermal solar atmosphere. We present 1D and 2D standard tests to demonstrate the quality of the numerical results obtained with our code. As 3D tests we present the propagation of MHD-gravity waves and vortices in the solar atmosphere. The code is based on high-resolution shock-capturing methods, uses the HLLE flux formula combined with Minmod, MC and WENO5 reconstructors. The divergence free magnetic field constraint is controlled using the Flux Constrained Transport method.
MHD Contractors' Review Meeting: Abstracts
NASA Astrophysics Data System (ADS)
The objectives of the Integrated Topping Cycle project are to design, construct, and deliver all prototypical hardware necessary to conduct long duration integrated MHD topping cycle proof-of-concept tests at the Component Development and Integration Facility (CDIF) in Butte, Montana. The results of the long duration tests will augment the existing engineering data base on MHD power train reliability, maintainability, durability, and performance, and will serve as a basis for scaling up to the early commercial scale power plant retrofit. The components of the MHD power train to be designed, fabricated, and tested include the following three systems: (1) a slagging coal combustion subsystem with a rated capacity of 50 MW thermal input, capable of operation with eastern (Illinois) or western (Montana Rosebud) coal; (2) a channel subsystem consisting of a segmented supersonic nozzle, channel (with current controls), and diffuser, capable of power output of 1.5 MW(sub e); and (3) a current consolidation subsystem to interface the channel with the existing facility inverter.
NASA Astrophysics Data System (ADS)
McLaughlin, J. A.; De Moortel, I.; Hood, A. W.; Brady, C. S.
2009-01-01
Context: This paper extends the models of Craig & McClymont (1991, ApJ, 371, L41) and McLaughlin & Hood (2004, A&A, 420, 1129) to include finite β and nonlinear effects. Aims: We investigate the nature of nonlinear fast magnetoacoustic waves about a 2D magnetic X-point. Methods: We solve the compressible and resistive MHD equations using a Lagrangian remap, shock capturing code (Arber et al. 2001, J. Comp. Phys., 171, 151) and consider an initial condition in {v}×{B} \\cdot {hat{z}} (a natural variable of the system). Results: We observe the formation of both fast and slow oblique magnetic shocks. The nonlinear wave deforms the X-point into a “cusp-like” point which in turn collapses to a current sheet. The system then evolves through a series of horizontal and vertical current sheets, with associated changes in connectivity, i.e. the system exhibits oscillatory reconnection. Our final state is non-potential (but in force balance) due to asymmetric heating from the shocks. Larger amplitudes in our initial condition correspond to larger values of the final current density left in the system. Conclusions: The inclusion of nonlinear terms introduces several new features to the system that were absent from the linear regime. A movie is available in electronic form at http://www.aanda.org
MHD simple waves and the divergence wave
Webb, G. M.; Pogorelov, N. V.; Zank, G. P.
2010-03-25
In this paper we investigate magnetohydrodynamic (MHD) simple divergence waves in MHD, for models in which nablacentre dotBnot =0. These models are related to the eight wave Riemann solvers in numerical MHD, in which the eighth wave is the divergence wave associated with nablacentre dotBnot =0. For simple wave solutions, all physical variables (the gas density, pressure, fluid velocity, entropy, and magnetic field induction in the MHD case) depend on a single phase function phi. We consider the form of the MHD equations used by both Powell et al. and Janhunen. It is shown that the Janhunen version of the equations possesses fully nonlinear, exact simple wave solutions for the divergence wave, but no physically meaningful simple divergence wave solution exists for the Powell et al. system. We suggest that the 1D simple, divergence wave solution for the Janhunen system, may be useful for the testing and validation of numerical MHD codes.
NKG2D ligands as therapeutic targets
Spear, Paul; Wu, Ming-Ru; Sentman, Marie-Louise; Sentman, Charles L.
2013-01-01
The Natural Killer Group 2D (NKG2D) receptor plays an important role in protecting the host from infections and cancer. By recognizing ligands induced on infected or tumor cells, NKG2D modulates lymphocyte activation and promotes immunity to eliminate ligand-expressing cells. Because these ligands are not widely expressed on healthy adult tissue, NKG2D ligands may present a useful target for immunotherapeutic approaches in cancer. Novel therapies targeting NKG2D ligands for the treatment of cancer have shown preclinical success and are poised to enter into clinical trials. In this review, the NKG2D receptor and its ligands are discussed in the context of cancer, infection, and autoimmunity. In addition, therapies targeting NKG2D ligands in cancer are also reviewed. PMID:23833565
Rapporteur report: MHD electric power plants
NASA Technical Reports Server (NTRS)
Seikel, G. R.
1980-01-01
Five US papers from the Proceedings of the Seventh International Conference on MHD Electrical Power Generation at the Massachusetts Institute of Technology are summarized. Results of the initial parametric phase of the US effort on the study of potential early commercial MHD plants are reported and aspects of the smaller commercial prototype plant termed the Engineering Test Facility are discussed. The alternative of using a disk geometry generator rather than a linear generator in baseload MHD plants is examined. Closed-cycle as well as open-cycle MHD plants are considered.
MHD activity in the ISX-B tokamak: experimental results and theoretical interpretation
Carreras, B.A.; Dunlap, J.L.; Bell, J.D.; Charlton, L.A.; Cooper, W.A.; Dory, R.A.; Hender, T.C.; Hicks, H.R.; Holmes, J.A.; Lynch, V.E.
1982-01-01
The observed spectrum of MHD fluctuations in the ISX-B tokamak is clearly dominated by the n=1 mode when the q=1 surface is in the plasma. This fact agrees well with theoretical predictions based on 3-D resistive MHD calculations. They show that the (m=1; n=1) mode is then the dominant instability. It drives other n=1 modes through toroidal coupling and n>1 modes through nonlinear couplings. These theoretically predicted mode structures have been compared in detail with the experimentally measured wave forms (using arrays of soft x-ray detectors). The agreement is excellent. More detailed comparisons between theory and experiment have required careful reconstructions of the ISX-B equilibria. The equilibria so constructed have permitted a precise evaluation of the ideal MHD stability properties of ISX-B. The present results indicate that the high ..beta.. ISX-B equilibria are marginally stable to finite eta ideal MHD modes. The resistive MHD calculations also show that at finite ..beta.. there are unstable resistive pressure driven modes.
Perspectives for spintronics in 2D materials
NASA Astrophysics Data System (ADS)
Han, Wei
2016-03-01
The past decade has been especially creative for spintronics since the (re)discovery of various two dimensional (2D) materials. Due to the unusual physical characteristics, 2D materials have provided new platforms to probe the spin interaction with other degrees of freedom for electrons, as well as to be used for novel spintronics applications. This review briefly presents the most important recent and ongoing research for spintronics in 2D materials.
Magnetotail dynamics: MHD simulations of driven and spontaneous dynamic changes
Birn, J.; Schindler, K.; Hesse, M.
1994-05-01
The dynamic evolution of the magnetotail during growth phase and expansion phase of a substorm is studied through threedimensional time-dependent MHD simulations. To model growth phase effects, an external electric field with an equatorward inflow is applied at the boundaries over a finite time period. This leads to the formation of a thin current sheet with greatly enhanced current density in the near tail, embedded in the wider plasma/current sheet, which becomes diminished in strength. A faster, spontaneous current sheet formation occurs when entropy conservation is released in an isobaric model, while the ideal MHD constraint persists. This may be a suitable model for the late, explosive part of the growth phase. The transition to the substorm expansive phase is modeled by an increase in anomalous resistivity, using either uniform resistivity or a current density dependent resistivity which is turned on when the current density exceeds a certain threshold. In both cases the violation of ideal MHD leads to resistive instability and the formation of a near-Earth neutral line, fast flow, and plasmoid ejection, together with the dipolarization and current reduction in the region further earthward. The spontaneous increase in total region 1 type field-aligned currents associated with the disruptions of the thin current sheets is less significant than that found in earlier simulations of the disruption of a wider current sheet, whereas the driven increase in the region 1 type current is substantial. The results demonstrate that the same dynamic process which appears spontaneous in the behavior of some quantities might be interpreted as entirely driven from the observation of others.
Kinematic MHD Models of Collapsing Magnetic Traps: Extension to 3D
Grady, Keith J.; Neukirch, Thomas
2009-02-16
We show how fully 3D kinematic MHD models of collapsing magnetic traps (CMTs) can be constructed, thus extending previous work on 2D trap models. CMTs are thought to form in the relaxing magnetic field lines in solar flares and it has been proposed that they play an important role in the acceleration of high-energy particles. This work is a first step to understanding the physics of CMTs better.
MHD Integrated Topping Cycle Project
Not Available
1992-07-01
This eighteenth quarterly technical progress report of the MHD Integrated Topping cycle Project presents the accomplishments during the period November 1, 1991 to January 31, 1992. The precombustor is fully assembled. Manufacturing of all slagging stage components has been completed. All cooling panels were welded in place and the panel/shell gap was filled with RTV. Final combustor assembly is in progress. The low pressure cooling subsystem (LPCS) was delivered to the CDIF. Second stage brazing issues were resolved. The construction of the two anode power cabinets was completed.
Ceramic component for MHD electrode
Marchant, David D.; Bates, Junior L.
1981-01-01
A ceramic component which exhibits electrical conductivity down to near room temperatures has the formula: Hf.sub.x In.sub.y A.sub.z O.sub.2 where x=0.1 to 0.4, y=0.3 to 0.6, z=0.1 to 0.4 and A is a lanthanide rare earth or yttrium. The component is suitable for use in the fabrication of MHD electrodes or as the current leadout portion of a composite electrode with other ceramic components.
Ceramic components for MHD electrode
Marchant, D.D.
A ceramic component which exhibits electrical conductivity down to near room temperatures has the formula: Hf/sub x/In/sub y/A/sub z/O/sub 2/ where x = 0.1 to 0.4, y = 0.3 to 0.6, z = 0.1 to 0.4 and A is a lanthanide rare earth or yttrium. The component is suitable for use in the fabrication of MHD electrodes or as the current leadout portion of a composite electrode with other ceramic components.
NASA Astrophysics Data System (ADS)
Kasprzak, Marek; Traczyk, Andrzej
2014-06-01
In urbanized areas, particularly in lowland terrains and floors of large river valleys, the natural land configuration is often hard to recognize due to a long history of human activity. Accordingly, archaeological works in cities, which supply knowledge on settlement conditions, are usually accompanied by geological and geomophological research. Lately, data from light detection and ranging (LiDAR) have become a valuable source of information on urban land configuration. Geophysical methods are also becoming increasingly popular in background studies. The paper presents a method of using and linking these sources of spatial information about landforms in such areas. The main aim is to identify to what extent these complementary sources of data and the proposed method can be used in such a specific environment to reconstruct natural, buried terrain morphology. The city of Wrocław in Central Europe serves as an example. To this end geomorphometric studies were conducted with the use of digital elevation models (DEMs) based on LiDAR scanning and derivated land-surface parameters—SAGA Wetness Index, Channel Network Base Level and Altitude above Channel Network. The study also involved determining morphological edges and measurements of the meanders of the Odra, as well as expanding information on the spatial distribution of alluvia and the structure of slope breaks. To this end, geophysical measurements were conducted using the Two-Dimensional Electrical Resistivity Tomography method. Additionally, five typical sequences of man-made ground present within the perimeter of the city were distinguished. As a result, a map of the main landforms of Wrocław is presented. Finally, we argue that although high resolution DEM and derivate land-surface parameters are very useful in terrain analysis, places with thick man-made ground or strongly levelled areas must be recognized by geoarchaeological excavations or geological bore holes. The geophysical survey is useful to
Numerical simulation of the operation of a MHD generator in transient regimes in MHD power stations
Bityurin, V.A.; Ivanov, P.P.; Koryagina, G.M.; Lyubimov, G.A.; Medin, S.A.; Morozov, G.N.; Prokop, A.S.
1982-09-01
Transient regimes of a MHD generator operating in combination with equipment in a MHD power station are analzyed with the help of a numerical model. The MHD generator, whose flow-through part consists of a nozzle, a channel, and a diffuser, is regulated by changing the flow rate and the load. Three types of MHD channels are studied: Faraday supersonic and subsonic, and diagonal supersonic. Their characteristics are presented and the efficiency of the MHD power station under nonrated regimes is determined. It is established that a MHD generator and the MHD power station as a whole admit quite efficient and deep regulation of the change in the flow rate of the working body.
Workshop on Feedback Stabilization of MHD Stabilities
McGuire, K.; Kugel, H.; La Haye, R.; Mauel, M.; Nevins, W.; Prager, S.
1996-12-31
The feedback stabilization of MHD instabilities is an area of research that is critical for improving the performance and economic attractiveness of magnetic confinement devices. A Workshop dedicated to feedback stabilization of MHD instabilities was held from December 11-13, 1996 at the Princeton Plasma Physics Laboratory, Princeton NJ, USA. The resulting presentations, conclusions, and recommendations are summarized.
MHD (Magnetohydrodynamics) Program Plan, FY 1989
NASA Astrophysics Data System (ADS)
1989-05-01
The essential elements of the current program, which is a continuation of the program outlined in the FY 1988 MHD Program Plan, are to: develop technical and environmental data for the integrated MHD topping cycle system through long-term (1000 hours) proof of concept (POC) testing; develop technical and environmental data for the integrated MHD bottoming cycle subsystem through long-term (4000 hours) POC testing; design and construct a seed regeneration system capable of independent operation, using spent seed materials from the MHD process; prepare a conceptual design for an MHD retrofit plant; and continue system studies and supporting research necessary for system testing. Results of the topping cycle POC tests at the Component Development and Integration Facility (CDIF), coupled with the bottoming cycle POC test results obtained at the Coal Fired Flow Facility (CFFF), and the seed regeneration POC effort will provide the critical engineering data base for the private sector's final decision on proceeding with the design, construction, and operation of an MHD retrofit. The development schedule, decision points, and resource requirements are discussed. As part of the MHD program, international activities of several nations are monitored and evaluated through contact with the international MHD scientific and technical community.
Annotated Bibliography of EDGE2D Use
J.D. Strachan and G. Corrigan
2005-06-24
This annotated bibliography is intended to help EDGE2D users, and particularly new users, find existing published literature that has used EDGE2D. Our idea is that a person can find existing studies which may relate to his intended use, as well as gain ideas about other possible applications by scanning the attached tables.
Staring 2-D hadamard transform spectral imager
Gentry, Stephen M.; Wehlburg, Christine M.; Wehlburg, Joseph C.; Smith, Mark W.; Smith, Jody L.
2006-02-07
A staring imaging system inputs a 2D spatial image containing multi-frequency spectral information. This image is encoded in one dimension of the image with a cyclic Hadamarid S-matrix. The resulting image is detecting with a spatial 2D detector; and a computer applies a Hadamard transform to recover the encoded image.
MHD Turbulence through the Heliosphere
NASA Astrophysics Data System (ADS)
Veltri, P.
Velocity and magnetic field fluctuations in a wide range of space and time scales have been directly detected in the interplanetary medium In the solar corona the presence of MHD turbulence is naturally generated by the mechanical and magnetic energy input from the photosphere and it could be related to coronal heating as well as to energy release events like micro and nanoflares A certain amount of fluctuations from the solar corona arrives in the solar wind mainly as Alfvénic turbulence i e strongly correlated velocity and magnetic field fluctuations with a very low level of compressible density magnetic field intensity temperature fluctuations The whole system formed by the solar corona and the solar wind represents a sort of wind tunnel extremely useful to study the MHD turbulence properties The presence of magnetic turbulence in the heliosphere is identified as the source of charged particle collisionless diffusion which according the values of parameters like the energy level on magnetic fluctuations or the turbulence correlation length can display both a normal gaussian random walk and an anomalous subdiffusive or super diffusive behavior The former case is obtained in a situation of global stochasticity high level of fluctuation energy while the latter in a situation of weak chaos low level of fluctuation energy The talk will discuss turbulence generation at photospheric level its propagation and its interaction with heliospheric structures and its effects on anomalous transport processes of charged
Feasibility of MHD submarine propulsion
Doss, E.D. ); Sikes, W.C. )
1992-09-01
This report describes the work performed during Phase 1 and Phase 2 of the collaborative research program established between Argonne National Laboratory (ANL) and Newport News Shipbuilding and Dry Dock Company (NNS). Phase I of the program focused on the development of computer models for Magnetohydrodynamic (MHD) propulsion. Phase 2 focused on the experimental validation of the thruster performance models and the identification, through testing, of any phenomena which may impact the attractiveness of this propulsion system for shipboard applications. The report discusses in detail the work performed in Phase 2 of the program. In Phase 2, a two Tesla test facility was designed, built, and operated. The facility test loop, its components, and their design are presented. The test matrix and its rationale are discussed. Representative experimental results of the test program are presented, and are compared to computer model predictions. In general, the results of the tests and their comparison with the predictions indicate that thephenomena affecting the performance of MHD seawater thrusters are well understood and can be accurately predicted with the developed thruster computer models.
Cusp geometry in MHD simulations
NASA Astrophysics Data System (ADS)
Siscoe, George; Crooker, Nancy; Siebert, Keith; Maynard, Nelson; Weimer, Daniel; White, Willard
2005-01-01
The MHD simulations described here show that the latitude of the high-altitude cusp decreases as the IMF swings from North to South, that there is a pronounced dawn dusk asymmetry at high-altitude associated with a dawn dusk component of the IMF, and that at the same time there is also a pronounced dawn dusk asymmetry at low-altitude. The simulations generate a feature that represents what has been called the cleft. It appears as a tail (when the IMF has a By component) attached to the cusp, extending either toward the dawn flank or the dusk flank depending on the dawn dusk orientation of the IMF. This one-sided cleft connects the cusp to the magnetospheric sash. We compare cusp geometry predicted by MHD simulations against published observations based on Hawkeye and DMSP data. Regarding the high-altitude predictions, the comparisons are not definitive, mainly because the observations are incomplete or mutually inconsistent. Regarding the low-altitude prediction of a strong dawn dusk asymmetry, the observations are unambiguous and are in good qualitative agreement with the prediction.
Light field morphing using 2D features.
Wang, Lifeng; Lin, Stephen; Lee, Seungyong; Guo, Baining; Shum, Heung-Yeung
2005-01-01
We present a 2D feature-based technique for morphing 3D objects represented by light fields. Existing light field morphing methods require the user to specify corresponding 3D feature elements to guide morph computation. Since slight errors in 3D specification can lead to significant morphing artifacts, we propose a scheme based on 2D feature elements that is less sensitive to imprecise marking of features. First, 2D features are specified by the user in a number of key views in the source and target light fields. Then the two light fields are warped view by view as guided by the corresponding 2D features. Finally, the two warped light fields are blended together to yield the desired light field morph. Two key issues in light field morphing are feature specification and warping of light field rays. For feature specification, we introduce a user interface for delineating 2D features in key views of a light field, which are automatically interpolated to other views. For ray warping, we describe a 2D technique that accounts for visibility changes and present a comparison to the ideal morphing of light fields. Light field morphing based on 2D features makes it simple to incorporate previous image morphing techniques such as nonuniform blending, as well as to morph between an image and a light field. PMID:15631126
2D materials for nanophotonic devices
NASA Astrophysics Data System (ADS)
Xu, Renjing; Yang, Jiong; Zhang, Shuang; Pei, Jiajie; Lu, Yuerui
2015-12-01
Two-dimensional (2D) materials have become very important building blocks for electronic, photonic, and phononic devices. The 2D material family has four key members, including the metallic graphene, transition metal dichalcogenide (TMD) layered semiconductors, semiconducting black phosphorous, and the insulating h-BN. Owing to the strong quantum confinements and defect-free surfaces, these atomically thin layers have offered us perfect platforms to investigate the interactions among photons, electrons and phonons. The unique interactions in these 2D materials are very important for both scientific research and application engineering. In this talk, I would like to briefly summarize and highlight the key findings, opportunities and challenges in this field. Next, I will introduce/highlight our recent achievements. We demonstrated atomically thin micro-lens and gratings using 2D MoS2, which is the thinnest optical component around the world. These devices are based on our discovery that the elastic light-matter interactions in highindex 2D materials is very strong. Also, I would like to introduce a new two-dimensional material phosphorene. Phosphorene has strongly anisotropic optical response, which creates 1D excitons in a 2D system. The strong confinement in phosphorene also enables the ultra-high trion (charged exciton) binding energies, which have been successfully measured in our experiments. Finally, I will briefly talk about the potential applications of 2D materials in energy harvesting.
Inertial solvation in femtosecond 2D spectra
NASA Astrophysics Data System (ADS)
Hybl, John; Albrecht Ferro, Allison; Farrow, Darcie; Jonas, David
2001-03-01
We have used 2D Fourier transform spectroscopy to investigate polar solvation. 2D spectroscopy can reveal molecular lineshapes beneath ensemble averaged spectra and freeze molecular motions to give an undistorted picture of the microscopic dynamics of polar solvation. The transition from "inhomogeneous" to "homogeneous" 2D spectra is governed by both vibrational relaxation and solvent motion. Therefore, the time dependence of the 2D spectrum directly reflects the total response of the solvent-solute system. IR144, a cyanine dye with a dipole moment change upon electronic excitation, was used to probe inertial solvation in methanol and propylene carbonate. Since the static Stokes' shift of IR144 in each of these solvents is similar, differences in the 2D spectra result from solvation dynamics. Initial results indicate that the larger propylene carbonate responds more slowly than methanol, but appear to be inconsistent with rotational estimates of the inertial response. To disentangle intra-molecular vibrations from solvent motion, the 2D spectra of IR144 will be compared to the time-dependent 2D spectra of the structurally related nonpolar cyanine dye HDITCP.
Internal Photoemission Spectroscopy of 2-D Materials
NASA Astrophysics Data System (ADS)
Nguyen, Nhan; Li, Mingda; Vishwanath, Suresh; Yan, Rusen; Xiao, Shudong; Xing, Huili; Cheng, Guangjun; Hight Walker, Angela; Zhang, Qin
Recent research has shown the great benefits of using 2-D materials in the tunnel field-effect transistor (TFET), which is considered a promising candidate for the beyond-CMOS technology. The on-state current of TFET can be enhanced by engineering the band alignment of different 2D-2D or 2D-3D heterostructures. Here we present the internal photoemission spectroscopy (IPE) approach to determine the band alignments of various 2-D materials, in particular SnSe2 and WSe2, which have been proposed for new TFET designs. The metal-oxide-2-D semiconductor test structures are fabricated and characterized by IPE, where the band offsets from the 2-D semiconductor to the oxide conduction band minimum are determined by the threshold of the cube root of IPE yields as a function of photon energy. In particular, we find that SnSe2 has a larger electron affinity than most semiconductors and can be combined with other semiconductors to form near broken-gap heterojunctions with low barrier heights which can produce a higher on-state current. The details of data analysis of IPE and the results from Raman spectroscopy and spectroscopic ellipsometry measurements will also be presented and discussed.
MHD waves and instabilities for gravitating, magnetized configurations in motion
NASA Astrophysics Data System (ADS)
Keppens, Rony; Goedbloed, Hans J. P.
Seismic probing of equilibrium configurations is of course well-known from geophysics, but has also been succesfully used to determine the internal structure of the Sun to an amazing accuracy. The results of helioseismology are quite impressive, although they only exploit an equilibrium structure where inward gravity is balanced by a pressure gradient in a 1D radial fashion. In principle, one can do the same for stationary, gravitating, magnetized plasma equilibria, as needed to perform MHD seismology in astrophysical jets or accretion disks. The introduction of (sheared) differential rotation does require the important switch from diagnosing static to stationary equilibrium configurations. The theory to describe all linear waves and instabilities in ideal MHD, given an exact stationary, gravitating, magnetized plasma equilibrium, in any dimensionality (1D, 2D, 3D) has been known since 1960, and is governed by the Frieman-Rotenberg equation. The full (mathematical) power of spectral theory governing physical eigenmode determination comes into play when using the Frieman-Rotenberg equation for moving equilibria, as applicable to astrophysical jets, accretion disks, but also solar flux ropes with stationary flow patterns. I will review exemplary seismic studies of flowing equilibrium configurations, covering solar to astrophysical configurations in motion. In that case, even essentially 1D configurations require quantification of the spectral web of eigenmodes, organizing the complex eigenfrequency plane.
Corrosion and arc erosion in MHD channels
NASA Astrophysics Data System (ADS)
Rosa, R. J.; Pollina, R. J.
1991-04-01
The objective of this task is to study the corrosion and arc erosion of MHD materials in a cooperative effort with, and to support, the MHD topping cycle program. Materials tested in the Avco Research Laboratory/Textron facility, or materials which have significant MHD importance, will be analyzed to document their physical deterioration. Conclusions shall be drawn about their wear mechanisms and lifetime in the MHD environment with respect to the following issues; sulfur corrosion, electrochemical corrosion, and arc erosion. The impact of any materials or slag conditions on the level of power output and on the level of leakage current in the MHD channel will also be noted, where appropriate. Two phenomena that can effect the analysis of slag leakage current have been investigated and found significant. These are: (1) transverse current along the slag layer in the insulator walls of an MHD duct, and (2) electrode surface voltage drops. Both tend to reduce the value inferred for average plasma conductivity and increase the value inferred for axial leakage current. These two effects in combination are potentially capable of explaining the high leakage inferred. Corrosion on the water side of metal MHD duct wall elements has been examined in CDIF and Mark 7 generators. It appears to be controllable by adjusting the pH of the water and/or by controlling the dissolved oxygen content.
2D radiation-magnetohydrodynamic simulations of SATURN imploding Z-pinches
Hammer, J.H.; Eddleman, J.L.; Springer, P.T.
1995-11-06
Z-pinch implosions driven by the SATURN device at Sandia National Laboratory are modeled with a 2D radiation magnetohydrodynamic (MHD) code, showing strong growth of magneto-Rayleigh Taylor (MRT) instability. Modeling of the linear and nonlinear development of MRT modes predicts growth of bubble-spike structures that increase the time span of stagnation and the resulting x-ray pulse width. Radiation is important in the pinch dynamics keeping the sheath relatively cool during the run-in and releasing most of the stagnation energy. The calculations give x-ray pulse widths and magnitudes in reasonable agreement with experiments, but predict a radiating region that is too dense and radially localized at stagnation. We also consider peaked initial density profiles with constant imploding sheath velocity that should reduce MRT instability and improve performance. 2D krypton simulations show an output x-ray power > 80 TW for the peaked profile.
MHD Technology Transfer, Integration and Review Committee
Not Available
1989-10-01
As part of the MHD Integrated Topping Cycle (ITC) project, TRW was given the responsibility to organize, charter and co-chair, with the Department of Energy (DOE), an MHD Technology Transfer, Integration and Review Committee (TTIRC). The Charter of the TTIRC, which was approved by the DOE in June 1988 and distributed to the committee members, is included as part of this Summary. As stated in the Charter, the purpose of this committee is to: (1) review all Proof-of-Concept (POC) projects and schedules in the national MHD program; to assess their compatibility with each other and the first commercial MHD retrofit plant; (2) establish and implement technology transfer formats for users of this technology; (3) identify interfaces, issues, and funding structures directly impacting the success of the commercial retrofit; (4) investigate and identify the manner in which, and by whom, the above should be resolved; and (5) investigate and assess other participation (foreign and domestic) in the US MHD Program. The DOE fiscal year 1989 MHD Program Plan Schedule is included at the end of this Summary. The MHD Technology Transfer, Integration and Review Committee's activities to date have focused primarily on the technology transfer'' aspects of its charter. It has provided a forum for the dissemination of technical and programmatic information among workers in the field of MHD and to the potential end users, the utilities, by holding semi-annual meetings. The committee publishes this semi-annual report, which presents in Sections 2 through 11 capsule summaries of technical progress for all DOE Proof-of-Concept MHD contracts and major test facilities.
MHD Simulations of Core Collapse Supernovae with Cosmos++
NASA Astrophysics Data System (ADS)
Akiyama, Shizuka; Salmonson, Jay
2010-10-01
We performed 2D, axisymmetric, MHD simulations with Cosmos++ in order to examine the growth of the magnetorotational instability (MRI) in core-collapse supernovae. We have initialized a non-rotating 15 Msolar progenitor, infused with differential rotation and poloidal magnetic fields. The collapse of the iron core is simulated with the Shen EOS, and the parametric Ye and entropy evolution. The wavelength of the unstable mode in the post-collapse environment is expected to be only ~200 m. In order to achieve the fine spatial resolution requirement, we employed remapping technique after the iron core has collapsed and bounced. The MRI unstable region appears near the equator and angular momentum and entropy are transported outward. Higher resolution remap run display more vigorous overturns and stronger transport of angular momentum and entropy. Our results are in agreement with the earlier work by Akiyama et al. [1] and Obergaulinger et al. [2].
MHD Code Optimizations and Jets in Dense Gaseous Halos
NASA Astrophysics Data System (ADS)
Gaibler, Volker; Vigelius, Matthias; Krause, Martin; Camenzind, Max
We have further optimized and extended the 3D-MHD-code NIRVANA. The magnetized part runs in parallel, reaching 19 Gflops per SX-6 node, and has a passively advected particle population. In addition, the code is MPI-parallel now - on top of the shared memory parallelization. On a 512^3 grid, we reach 561 Gflops with 32 nodes on the SX-8. Also, we have successfully used FLASH on the Opteron cluster. Scientific results are preliminary so far. We report one computation of highly resolved cocoon turbulence. While we find some similarities to earlier 2D work by us and others, we note a strange reluctancy of cold material to enter the low density cocoon, which has to be investigated further.
Application of Magnetohydrodynamics (MHD) and Recent Research Trend
NASA Astrophysics Data System (ADS)
Harada, Nobuhiro
As the applications of Magnetohydrodynamic (MHD) energy conversion, research and development for high-efficiency and low emission electric power generation system, MHD accelerations and/or MHD thrusters, and flow control around hypersonic and re-entry vehicles are introduced. For closed cycle MHD power generation, high-efficiency MHD single system is the most hopeful system and space power system using mixed inert gas (MIG) working medium is proposed. For open cycle MHD, high-efficiency coal fired MHD system with CO2 recovery has been proposed. As inverse process of MHD power generation, MHD accelerators/thrusters are expected as the next generation propulsion system. Heat flux reduction to protect re-entry vehicles is expected by an MHD process for safety return from space missions.
Alfven Wave Tomography for Cold MHD Plasmas
I.Y. Dodin; N.J. Fisch
2001-09-07
Alfven waves propagation in slightly nonuniform cold plasmas is studied by means of ideal magnetohydrodynamics (MHD) nonlinear equations. The evolution of the MHD spectrum is shown to be governed by a matrix linear differential equation with constant coefficients determined by the spectrum of quasi-static plasma density perturbations. The Alfven waves are shown not to affect the plasma density inhomogeneities, as they scatter off of them. The application of the MHD spectrum evolution equation to the inverse scattering problem allows tomographic measurements of the plasma density profile by scanning the plasma volume with Alfven radiation.
Local structures of homogeneous Hall MHD turbulence
NASA Astrophysics Data System (ADS)
Miura, H.; Araki, K.
2011-12-01
Local structures of decaying homogeneous and isotropic Hall MHD turbulence are studied by means of direct numerical simulations. Regions of strong vorticity and strong current density in Hall MHD turbulence are compared to those of single-fluid MHD turbulence. An analysis by the use of a low-pass filter reveals that the introduction of the Hall term can modify not only small-scale structures of the current density but also structures of the vorticity field, especially at the scales smaller than the ion skin depth.
NASA Astrophysics Data System (ADS)
Brown, Michael R.; Schaffner, David A.
2015-06-01
A new turbulent plasma source at the Swarthmore Spheromak Experiment (SSX) facility is described. The MHD wind tunnel configuration employs a magnetized plasma gun to inject high-beta plasma into a large, well-instrumented, vacuum drift region. This provides unique laboratory conditions approaching that in the solar wind: there is no applied background magnetic field in the drift region and has no net axial magnetic flux; the plasma flow speed is on the order of the local sound speed (M ~ 1), so flow energy density is comparable to thermal energy density; and the ratio of thermal to magnetic pressure is of order unity (plasma β ~ 1) so thermal energy density is also comparable to magnetic energy density. Results presented here and referenced within demonstrate the new capabilities and show how the new platform is proving useful for fundamental plasma turbulence studies.
MHD Integrated Topping Cycle Project
Not Available
1992-07-01
This seventeenth quarterly technical progress report of the MHD Integrated Topping Cycle Project presents the accomplishments during the period August 1, 1991 to October 31, 1991. Manufacturing of the prototypical combustor pressure shell has been completed including leak, proof, and assembly fit checking. Manufacturing of forty-five cooling panels was also completed including leak, proof, and flow testing. All precombustor internal components (combustion can baffle and swirl box) were received and checked, and integration of the components was initiated. A decision was made regarding the primary and backup designs for the 1A4 channel. The assembly of the channel related prototypical hardware continued. The cathode wall electrical wiring is now complete. The mechanical design of the diffuser has been completed.
Cosmological AMR MHD with Enzo
Xu, Hao; Li, Hui; Li, Shengtai
2009-01-01
In this work, we present EnzoMHD, the extension of the cosmological code Enzoto include magnetic fields. We use the hyperbolic solver of Li et al. (2008) for the computation of interface fluxes. We use constrained transport methods of Balsara & Spicer (1999) and Gardiner & Stone (2005) to advance the induction equation, the reconstruction technique of Balsara (2001) to extend the Adaptive Mesh Refinement of Berger & Colella (1989) already used in Enzo, though formulated in a slightly different way for ease of implementation. This combination of methods preserves the divergence of the magnetic field to machine precision. We use operator splitting to include gravity and cosmological expansion. We then present a series of cosmological and non cosmologjcal tests problems to demonstrate the quality of solution resulting from this combination of solvers.
MHD seed recovery and regeneration
NASA Astrophysics Data System (ADS)
1988-10-01
The TRW Econoseed MHD Seed Regeneration Process is based on the reaction of calcium formate with potassium sulfate spent seed from an MHD electric power generation plant. The process was tested at bench scale, design a proof of concept (POC) test plant, plan and cost a Phase 2 project for a POC plant evaluation and prepare a conceptual design of a 300 MW (t) commercial plant. The results of the project are as follows: (1) each of the unit operations is demonstrated, and (2) the data are incorporated into a POC plant design and project cost, as well as a 300 MW (t) commercial retrofit plant design and cost estimate. Specific results are as follows: (1) calcium formate can be produced at 100 percent yield in a total retention time of less than 5 minutes, (2) utilizing the calcium formate, spent seed can quantitatively be converted to potassium formate, potassium carbonate or mixtures of these with potassium sulfate as per the commercial design without measurable loss of potassium to insolubles at a total retention time under 20 minutes and ambient pressure, (3) the solid rejects form the process meet RCRA EP Toxicity requirements for safe disposal, and (4) filtration and evaporation data, as well as reaction data cited above, show that the Econoseed technology is ready for scale up to POC plant scale. Economics forecast studies show that the total cost per unit of potassium for seed regeneration by the Econoseed Process is in the range of $0.23 to $0.27/lb, a cost which is less than half the potassium cost of $0.63/lb for purchasing new potassium carbonate.
MHD Integrated Topping Cycle Project
Not Available
1992-02-01
A summary of the work is excerpted here. Final design of an MHD channel for the ITC program POC test has been completed. The channel was designed to be capable of 1.5 MW {sub e} power output and a lifetime of 2000 hours. Emphasis was placed upon durability and reliability. Hence, specific measures were taken to design against channel damage due to electric faults. The life-limiting issues associated with electrochemical corrosion and erosion of gas-side surfaces were addressed by the use of various materials with proven wear characteristics in a coal-fired MHD channel environment. Pitting of prototypical sidewall coupons was observed in the CDIF workhorse testing. The most likely cause of the observed pitting, water leaks resulting from cooling water tube braze failures, has been remedied. New brazing procedures and isolation of the sidebar gas-side material from water contact will prevent sidebar pitting in the prototypical channel. Water-side corrosion tests reported in this quarterly report include the latest results of tungsten-copper elements at controlled pH, heat flux and voltage levels. In the combustion subsystem, efforts continued to focus on understanding and improving the current levels of slag recovery and seed utilization achieved by the combustor. Analytical support was also provided in the areas of slag rejection system operation, precombustor operation, and oil burner design modification. Channel data analysis activities continued in support of prototypical coupon testing at the CDIF. Analyses are presented on channel wall slagging behavior and sidewall voltage distributions.
Brittle damage models in DYNA2D
Faux, D.R.
1997-09-01
DYNA2D is an explicit Lagrangian finite element code used to model dynamic events where stress wave interactions influence the overall response of the system. DYNA2D is often used to model penetration problems involving ductile-to-ductile impacts; however, with the advent of the use of ceramics in the armor-anti-armor community and the need to model damage to laser optics components, good brittle damage models are now needed in DYNA2D. This report will detail the implementation of four brittle damage models in DYNA2D, three scalar damage models and one tensor damage model. These new brittle damage models are then used to predict experimental results from three distinctly different glass damage problems.
Ginsparg, P.
1991-01-01
These are introductory lectures for a general audience that give an overview of the subject of matrix models and their application to random surfaces, 2d gravity, and string theory. They are intentionally 1.5 years out of date.
Ginsparg, P.
1991-12-31
These are introductory lectures for a general audience that give an overview of the subject of matrix models and their application to random surfaces, 2d gravity, and string theory. They are intentionally 1.5 years out of date.
2D electronic materials for army applications
NASA Astrophysics Data System (ADS)
O'Regan, Terrance; Perconti, Philip
2015-05-01
The record electronic properties achieved in monolayer graphene and related 2D materials such as molybdenum disulfide and hexagonal boron nitride show promise for revolutionary high-speed and low-power electronic devices. Heterogeneous 2D-stacked materials may create enabling technology for future communication and computation applications to meet soldier requirements. For instance, transparent, flexible and even wearable systems may become feasible. With soldier and squad level electronic power demands increasing, the Army is committed to developing and harnessing graphene-like 2D materials for compact low size-weight-and-power-cost (SWAP-C) systems. This paper will review developments in 2D electronic materials at the Army Research Laboratory over the last five years and discuss directions for future army applications.
2-d Finite Element Code Postprocessor
Energy Science and Technology Software Center (ESTSC)
1996-07-15
ORION is an interactive program that serves as a postprocessor for the analysis programs NIKE2D, DYNA2D, TOPAZ2D, and CHEMICAL TOPAZ2D. ORION reads binary plot files generated by the two-dimensional finite element codes currently used by the Methods Development Group at LLNL. Contour and color fringe plots of a large number of quantities may be displayed on meshes consisting of triangular and quadrilateral elements. ORION can compute strain measures, interface pressures along slide lines, reaction forcesmore » along constrained boundaries, and momentum. ORION has been applied to study the response of two-dimensional solids and structures undergoing finite deformations under a wide variety of large deformation transient dynamic and static problems and heat transfer analyses.« less
Distribution of temperature in the slot of an induction MHD pump
NASA Astrophysics Data System (ADS)
Krishbergs, R.
2007-09-01
At thermal calculations of the induction MHD machines, the mean temperatures of main parts of the machine are often defined. Yet, with the increase of specific electromagnetic loads, with the use of new heat-resistant insulating materials and with the intention to use electrically short windings (a small number of poles that causes asymmetry of phase currents) in the pumps of mean power determining are the local temperature peaks in the winding, though the winding mean temperature is far from limiting. The paper considers the temperature distribution over the slot height in the induction MHD pump. Figs 2, Refs 13.
Chemical Approaches to 2D Materials.
Samorì, Paolo; Palermo, Vincenzo; Feng, Xinliang
2016-08-01
Chemistry plays an ever-increasing role in the production, functionalization, processing and applications of graphene and other 2D materials. This special issue highlights a selection of enlightening chemical approaches to 2D materials, which nicely reflect the breadth of the field and convey the excitement of the individuals involved in it, who are trying to translate graphene and related materials from the laboratory into a real, high-impact technology. PMID:27478083
Extended 2D generalized dilaton gravity theories
NASA Astrophysics Data System (ADS)
de Mello, R. O.
2008-09-01
We show that an anomaly-free description of matter in (1+1) dimensions requires a deformation of the 2D relativity principle, which introduces a non-trivial centre in the 2D Poincaré algebra. Then we work out the reduced phase space of the anomaly-free 2D relativistic particle, in order to show that it lives in a noncommutative 2D Minkowski space. Moreover, we build a Gaussian wave packet to show that a Planck length is well defined in two dimensions. In order to provide a gravitational interpretation for this noncommutativity, we propose to extend the usual 2D generalized dilaton gravity models by a specific Maxwell component, which guages the extra symmetry associated with the centre of the 2D Poincaré algebra. In addition, we show that this extension is a high energy correction to the unextended dilaton theories that can affect the topology of spacetime. Further, we couple a test particle to the general extended dilaton models with the purpose of showing that they predict a noncommutativity in curved spacetime, which is locally described by a Moyal star product in the low energy limit. We also conjecture a probable generalization of this result, which provides strong evidence that the noncommutativity is described by a certain star product which is not of the Moyal type at high energies. Finally, we prove that the extended dilaton theories can be formulated as Poisson Sigma models based on a nonlinear deformation of the extended Poincaré algebra.
2D-2D tunneling field-effect transistors using WSe2/SnSe2 heterostructures
NASA Astrophysics Data System (ADS)
Roy, Tania; Tosun, Mahmut; Hettick, Mark; Ahn, Geun Ho; Hu, Chenming; Javey, Ali
2016-02-01
Two-dimensional materials present a versatile platform for developing steep transistors due to their uniform thickness and sharp band edges. We demonstrate 2D-2D tunneling in a WSe2/SnSe2 van der Waals vertical heterojunction device, where WSe2 is used as the gate controlled p-layer and SnSe2 is the degenerately n-type layer. The van der Waals gap facilitates the regulation of band alignment at the heterojunction, without the necessity of a tunneling barrier. ZrO2 is used as the gate dielectric, allowing the scaling of gate oxide to improve device subthreshold swing. Efficient gate control and clean interfaces yield a subthreshold swing of ˜100 mV/dec for >2 decades of drain current at room temperature, hitherto unobserved in 2D-2D tunneling devices. The subthreshold swing is independent of temperature, which is a clear signature of band-to-band tunneling at the heterojunction. A maximum switching ratio ION/IOFF of 107 is obtained. Negative differential resistance in the forward bias characteristics is observed at 77 K. This work bodes well for the possibilities of two-dimensional materials for the realization of energy-efficient future-generation electronics.
3D simulation studies of tokamak plasmas using MHD and extended-MHD models
Park, W.; Chang, Z.; Fredrickson, E.; Fu, G.Y.
1996-12-31
The M3D (Multi-level 3D) tokamak simulation project aims at the simulation of tokamak plasmas using a multi-level tokamak code package. Several current applications using MHD and Extended-MHD models are presented; high-{beta} disruption studies in reversed shear plasmas using the MHD level MH3D code, {omega}{sub *i} stabilization and nonlinear island saturation of TAE mode using the hybrid particle/MHD level MH3D-K code, and unstructured mesh MH3D{sup ++} code studies. In particular, three internal mode disruption mechanisms are identified from simulation results which agree which agree well with experimental data.
2D FEM Heat Transfer & E&M Field Code
Energy Science and Technology Software Center (ESTSC)
1992-04-02
TOPAZ and TOPAZ2D are two-dimensional implicit finite element computer codes for heat transfer analysis. TOPAZ2D can also be used to solve electrostatic and magnetostatic problems. The programs solve for the steady-state or transient temperature or electrostatic and magnetostatic potential field on two-dimensional planar or axisymmetric geometries. Material properties may be temperature or potential-dependent and either isotropic or orthotropic. A variety of time and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation.more » By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functional representation of boundary conditions and internal heat generation. The programs can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.« less
2D FEM Heat Transfer & E&M Field Code
1992-04-02
TOPAZ and TOPAZ2D are two-dimensional implicit finite element computer codes for heat transfer analysis. TOPAZ2D can also be used to solve electrostatic and magnetostatic problems. The programs solve for the steady-state or transient temperature or electrostatic and magnetostatic potential field on two-dimensional planar or axisymmetric geometries. Material properties may be temperature or potential-dependent and either isotropic or orthotropic. A variety of time and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation. By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functional representation of boundary conditions and internal heat generation. The programs can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.
Open Boundary Conditions for Dissipative MHD
Meier, E T
2011-11-10
In modeling magnetic confinement, astrophysics, and plasma propulsion, representing the entire physical domain is often difficult or impossible, and artificial, or 'open' boundaries are appropriate. A novel open boundary condition (BC) for dissipative MHD, called Lacuna-based open BC (LOBC), is presented. LOBC, based on the idea of lacuna-based truncation originally presented by V.S. Ryaben'kii and S.V. Tsynkov, provide truncation with low numerical noise and minimal reflections. For hyperbolic systems, characteristic-based BC (CBC) exist for separating the solution into outgoing and incoming parts. In the hyperbolic-parabolic dissipative MHD system, such separation is not possible, and CBC are numerically unstable. LOBC are applied in dissipative MHD test problems including a translating FRC, and coaxial-electrode plasma acceleration. Solution quality is compared to solutions using CBC and zero-normal derivative BC. LOBC are a promising new open BC option for dissipative MHD.
[Nonlinear magnetohydrodynamics]. [Threshold unstable MHD activity
Not Available
1992-01-01
Theoretical predictions were compared with available data from JET on the threshold unstable MHD activity in toroidal confinement devices. In particular, questions arising as to Hartmans number and the selection of a kinematic viscosity are discussed.
Solar driven liquid metal MHD power generator
NASA Astrophysics Data System (ADS)
Lee, J. H.; Hohl, F.
1983-06-01
A solar energy collector focuses solar energy onto a solar oven which is attached to a mixer which in turn is attached to the channel of a MHD generator. Gas enters the oven and a liquid metal enters the mixer. The gas/liquid metal mixture is heated by the collected solar energy and moves through the MHD generator thereby generating electrical power. The mixture is then separated and recycled.
Solar driven liquid metal MHD power generator
NASA Technical Reports Server (NTRS)
Lee, J. H.; Hohl, F. (Inventor)
1983-01-01
A solar energy collector focuses solar energy onto a solar oven which is attached to a mixer which in turn is attached to the channel of a MHD generator. Gas enters the oven and a liquid metal enters the mixer. The gas/liquid metal mixture is heated by the collected solar energy and moves through the MHD generator thereby generating electrical power. The mixture is then separated and recycled.
Results from the Pamir-3U pulsed portable MHD power system program
Swallom, D.W.; Goldfarb, V.M.; Gibbs, J.S.
1996-12-31
The Pamir-3U MHD Power System Acceptance Test Program was successfully performed in the United States in 1995. The Pamir-3U MHD system is a portable power system that can be transported to various operational locations. The power system is self-contained and does not require extensive support equipment to generate the design power. The work was performed by the IVTAN-Association of the Russian Academy of Sciences (Russia) and Textron Systems Division (US). Major subcontractors that participated in this effort are Nizhny Novgorod Machine Building Plant (Russia), Federal Center of Dual Technologies Soyuz (Russia) and Aerojet Corporation (US). A preliminary acceptance test program, consisting of five power tests and several preliminary tests, was conducted during August 1994 at Geodesiya Research and Development Institute, Krasnoarmejsk, Russia. During this test program, power levels as high as 15 MW(e) were obtained. Operation of the facility in various operating modes was demonstrated, and several tests were conducted where the resistance was varied during the hot-fire test run. For the final Pamir-3U MHD Power System Acceptance Test Program, eight hot-fire tests were performed. The performance levels of the Pamir-3U MHD Power System were confirmed during these tests. As a result of the experiments performed, the operational ability of the Pamir-3U MHD power system to operate in a variety of performance modes and under a variety of operating conditions was confirmed.
HYBRID AND HALL-MHD SIMULATIONS OF COLLISIONLESS RECONNECTION: EFFECTS OF PLASMA PRESSURE TENSOR
L. YIN; D. WINSKE; ET AL
2001-05-01
In this study we performed two-dimensional hybrid (particle ions, massless fluid electrons) and Hall-MHD simulations of collisionless reconnection in a thin current sheet. Both calculations include the full electron pressure tensor (instead of a localized resistivity) in the generalized Ohm's law to initiate reconnection, and in both an initial perturbation to the Harris equilibrium is applied. First, electron dynamics from the two calculations are compared, and we find overall agreement between the two calculations in both the reconnection rate and the global configuration. To address the issue of how kinetic treatment for the ions affects the reconnection dynamics, we compared the fluid-ion dynamics from the Hall-MHD calculation to the particle-ion dynamics obtained from the hybrid simulation. The comparison demonstrates that off-diagonal elements of the ion pressure tensor are important in correctly modeling the ion out-of-plane momentum transport from the X point. It is that these effects can be modeled efficiently using a particle Hall-MHD simulation method in which particle ions used in a predictor/corrector to implement the ion gyro-radius corrections. We also investigate the micro- macro-scale coupling in the magnetotail dynamics by using a new integrated approach in which particle Hall-MHD calculations are embedded inside a MHD simulation. Initial results of the simulation concerning current sheet thinning and reconnection dynamics are discussed.
MHD technology transfer, integration, and review committee
NASA Astrophysics Data System (ADS)
1990-05-01
As part of Task 8 of the magnetohydrodynamic (MHD) Integrated Topping Cycle (ITC) project, TRW was given the responsibility to organize, charter and co-chair, with the Department of Energy (DOE), an MHD Technology Transfer, Integration and Review Committee (TTIRC). The TTIRC consists of an Executive Committee (EC) which acts as the governing body, and a General Committee (GC), also referred to as the main or full committee, consisting of representatives from the various POC contractors, participating universities and national laboratories, utilities, equipment suppliers, and other potential MHD users or investors. The purpose of the TTIRC is to: (1) review all Proof-of-Concept (POC) projects and schedules in the national MHD program; to assess their compatibility with each other and the first commercial MHD retrofit plant; (2) establish and implement technology transfer formats for users of this technology; (3) identify interfaces, issues, and funding structures directly impacting the success of the commercial retrofit; (4) investigate and identify the manner in which, and by whom, the above should be resolved; and (5) investigate and assess other participation (foreign and domestic) in the U.S. MHD Program. There are seven sections: introduction; Executive Committee and General Committee activity; Committee activities related to technology transfer; ongoing POC integration activities being performed under the auspices of the Executive Committee; recommendations passed on to the DOE by the Executive Committee; Planned activities for the next six months.
Numerical simulations of resistive magnetohydrodynamic instabilities in a poloidal divertor tokamak
NASA Astrophysics Data System (ADS)
Uchimoto, E.
1988-03-01
A new 3-D resistive MHD initial value code RPD has been successfully developed from scratch to study the linear and nonlinear evolution of long wavelength resistive MHD instabilities in a square cross-section tokamak with or without a poloidal divertor. The code numerically advances the full set of compressible resistive MHD equations in a toroidal geometry, with an important option of permitting the divertor separatrix and the region outside it to be in the computational domain. A severe temporal step size restriction for numerical stability imposed by the fast compressional waves was removed by developing and implementing a new, efficient semi-implicit scheme extending one first proposed by Harned and Kerner. As a result, the code typically runs faster than that with a mostly explicit scheme by a factor of about the aspect ratio. The equilibrium input for RPD is generated by a new 2-D code EQPD that is based on the Chodura-Schluter method. The RPD code, as well as the new semi-implicit scheme, has passed very extensive numerical tests in both divertor and divertorless geometries. Linear and nonlinear simulations in a divertorless geometry have reproduced the standard, previously known results. In a geometry with a four-node divertor the m = 2, n = 1 (2/1) tearing mode tends to be linearly stabilized as the q = 2 surface approaches the divertor separatrix. However, the m = 1, n = 1 (1/1) resistive kink mode remains relatively unaffected by the nearness of the q = 1 surface to the divertor separatrix. When plasma current is added to the region outside the divertor separatrix, the 2/1 tearing mode is linearly stabilized not by this current, but by the profile modifications induced near the q = 2 surface and the divertor separatrix. A similar stabilization effect is seen for the 1/1 resistive kink mode, but to a lesser extent.
Optical modulators with 2D layered materials
NASA Astrophysics Data System (ADS)
Sun, Zhipei; Martinez, Amos; Wang, Feng
2016-04-01
Light modulation is an essential operation in photonics and optoelectronics. With existing and emerging technologies increasingly demanding compact, efficient, fast and broadband optical modulators, high-performance light modulation solutions are becoming indispensable. The recent realization that 2D layered materials could modulate light with superior performance has prompted intense research and significant advances, paving the way for realistic applications. In this Review, we cover the state of the art of optical modulators based on 2D materials, including graphene, transition metal dichalcogenides and black phosphorus. We discuss recent advances employing hybrid structures, such as 2D heterostructures, plasmonic structures, and silicon and fibre integrated structures. We also take a look at the future perspectives and discuss the potential of yet relatively unexplored mechanisms, such as magneto-optic and acousto-optic modulation.
Large Area Synthesis of 2D Materials
NASA Astrophysics Data System (ADS)
Vogel, Eric
Transition metal dichalcogenides (TMDs) have generated significant interest for numerous applications including sensors, flexible electronics, heterostructures and optoelectronics due to their interesting, thickness-dependent properties. Despite recent progress, the synthesis of high-quality and highly uniform TMDs on a large scale is still a challenge. In this talk, synthesis routes for WSe2 and MoS2 that achieve monolayer thickness uniformity across large area substrates with electrical properties equivalent to geological crystals will be described. Controlled doping of 2D semiconductors is also critically required. However, methods established for conventional semiconductors, such as ion implantation, are not easily applicable to 2D materials because of their atomically thin structure. Redox-active molecular dopants will be demonstrated which provide large changes in carrier density and workfunction through the choice of dopant, treatment time, and the solution concentration. Finally, several applications of these large-area, uniform 2D materials will be described including heterostructures, biosensors and strain sensors.
2D microwave imaging reflectometer electronics
Spear, A. G.; Domier, C. W. Hu, X.; Muscatello, C. M.; Ren, X.; Luhmann, N. C.; Tobias, B. J.
2014-11-15
A 2D microwave imaging reflectometer system has been developed to visualize electron density fluctuations on the DIII-D tokamak. Simultaneously illuminated at four probe frequencies, large aperture optics image reflections from four density-dependent cutoff surfaces in the plasma over an extended region of the DIII-D plasma. Localized density fluctuations in the vicinity of the plasma cutoff surfaces modulate the plasma reflections, yielding a 2D image of electron density fluctuations. Details are presented of the receiver down conversion electronics that generate the in-phase (I) and quadrature (Q) reflectometer signals from which 2D density fluctuation data are obtained. Also presented are details on the control system and backplane used to manage the electronics as well as an introduction to the computer based control program.
2D microwave imaging reflectometer electronics
NASA Astrophysics Data System (ADS)
Spear, A. G.; Domier, C. W.; Hu, X.; Muscatello, C. M.; Ren, X.; Tobias, B. J.; Luhmann, N. C.
2014-11-01
A 2D microwave imaging reflectometer system has been developed to visualize electron density fluctuations on the DIII-D tokamak. Simultaneously illuminated at four probe frequencies, large aperture optics image reflections from four density-dependent cutoff surfaces in the plasma over an extended region of the DIII-D plasma. Localized density fluctuations in the vicinity of the plasma cutoff surfaces modulate the plasma reflections, yielding a 2D image of electron density fluctuations. Details are presented of the receiver down conversion electronics that generate the in-phase (I) and quadrature (Q) reflectometer signals from which 2D density fluctuation data are obtained. Also presented are details on the control system and backplane used to manage the electronics as well as an introduction to the computer based control program.
2D microwave imaging reflectometer electronics.
Spear, A G; Domier, C W; Hu, X; Muscatello, C M; Ren, X; Tobias, B J; Luhmann, N C
2014-11-01
A 2D microwave imaging reflectometer system has been developed to visualize electron density fluctuations on the DIII-D tokamak. Simultaneously illuminated at four probe frequencies, large aperture optics image reflections from four density-dependent cutoff surfaces in the plasma over an extended region of the DIII-D plasma. Localized density fluctuations in the vicinity of the plasma cutoff surfaces modulate the plasma reflections, yielding a 2D image of electron density fluctuations. Details are presented of the receiver down conversion electronics that generate the in-phase (I) and quadrature (Q) reflectometer signals from which 2D density fluctuation data are obtained. Also presented are details on the control system and backplane used to manage the electronics as well as an introduction to the computer based control program. PMID:25430247