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.
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.
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.
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.
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)
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).
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.
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.
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.
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.
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.
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.
3-D Relativistic MHD Simulations
NASA Astrophysics Data System (ADS)
Nishikawa, K.-I.; Frank, J.; Koide, S.; Sakai, J.-I.; Christodoulou, D. M.; Sol, H.; Mutel, R. L.
1998-12-01
We present 3-D numerical simulations of moderately hot, supersonic jets propagating initially along or obliquely to the field lines of a denser magnetized background medium with Lorentz factors of W = 4.56 and evolving in a four-dimensional spacetime. 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 in the simulations. This 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.
3-D Relativistic MHD Simulations
NASA Astrophysics Data System (ADS)
Nishikaw, K.-I.; Frank, J.; Christodoulou, D. M.; Koide, S.; Sakai, J.-I.; Sol, H.; Mutel, R. L.
1998-12-01
We present 3-D numerical simulations of moderately hot, supersonic jets propagating initially along or obliquely to the field lines of a denser magnetized background medium with Lorentz factors of W=4.56 and evolving in a four-dimensional spacetime. 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 in the simulations. This 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. We also simulate jets with the more realistic initial conditions for injecting jets for helical mangetic field, perturbed density, velocity, and internal energy, which are supposed to be caused in the process of jet generation. Three possible explanations for the observed variability are (i) tidal disruption of a star falling into the black hole, (ii) instabilities in the relativistic accretion disk, and (iii) jet-related PRocesses. New results will be reported at the meeting.
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.
Realistic and efficient 2D crack simulation
NASA Astrophysics Data System (ADS)
Yadegar, Jacob; Liu, Xiaoqing; Singh, Abhishek
2010-04-01
Although numerical algorithms for 2D crack simulation have been studied in Modeling and Simulation (M&S) and computer graphics for decades, realism and computational efficiency are still major challenges. In this paper, we introduce a high-fidelity, scalable, adaptive and efficient/runtime 2D crack/fracture simulation system by applying the mathematically elegant Peano-Cesaro triangular meshing/remeshing technique to model the generation of shards/fragments. The recursive fractal sweep associated with the Peano-Cesaro triangulation provides efficient local multi-resolution refinement to any level-of-detail. The generated binary decomposition tree also provides efficient neighbor retrieval mechanism used for mesh element splitting and merging with minimal memory requirements essential for realistic 2D fragment formation. Upon load impact/contact/penetration, a number of factors including impact angle, impact energy, and material properties are all taken into account to produce the criteria of crack initialization, propagation, and termination leading to realistic fractal-like rubble/fragments formation. The aforementioned parameters are used as variables of probabilistic models of cracks/shards formation, making the proposed solution highly adaptive by allowing machine learning mechanisms learn the optimal values for the variables/parameters based on prior benchmark data generated by off-line physics based simulation solutions that produce accurate fractures/shards though at highly non-real time paste. Crack/fracture simulation has been conducted on various load impacts with different initial locations at various impulse scales. The simulation results demonstrate that the proposed system has the capability to realistically and efficiently simulate 2D crack phenomena (such as window shattering and shards generation) with diverse potentials in military and civil M&S applications such as training and mission planning.
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.
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.
MHD simulation of the Bastille day event
NASA Astrophysics Data System (ADS)
Linker, Jon; Torok, Tibor; Downs, Cooper; Lionello, Roberto; Titov, Viacheslav; Caplan, Ronald M.; Mikić, Zoran; Riley, Pete
2016-03-01
We describe a time-dependent, thermodynamic, three-dimensional MHD simulation of the July 14, 2000 coronal mass ejection (CME) and flare. The simulation starts with a background corona developed using an MDI-derived magnetic map for the boundary condition. Flux ropes using the modified Titov-Demoulin (TDm) model are used to energize the pre-event active region, which is then destabilized by photospheric flows that cancel flux near the polarity inversion line. More than 1033 ergs are impulsively released in the simulated eruption, driving a CME at 1500 km/s, close to the observed speed of 1700km/s. The post-flare emission in the simulation is morphologically similar to the observed post-flare loops. The resulting flux rope that propagates to 1 AU is similar in character to the flux rope observed at 1 AU, but the simulated ICME center passes 15° north of Earth.
Simulation of Yeast Cooperation in 2D.
Wang, M; Huang, Y; Wu, Z
2016-03-01
Evolution of cooperation has been an active research area in evolutionary biology in decades. An important type of cooperation is developed from group selection, when individuals form spatial groups to prevent them from foreign invasions. In this paper, we study the evolution of cooperation in a mixed population of cooperating and cheating yeast strains in 2D with the interactions among the yeast cells restricted to their small neighborhoods. We conduct a computer simulation based on a game theoretic model and show that cooperation is increased when the interactions are spatially restricted, whether the game is of a prisoner's dilemma, snow drifting, or mutual benefit type. We study the evolution of homogeneous groups of cooperators or cheaters and describe the conditions for them to sustain or expand in an opponent population. We show that under certain spatial restrictions, cooperator groups are able to sustain and expand as group sizes become large, while cheater groups fail to expand and keep them from collapse. PMID:26988702
Shocked Magnetotail: ARTEMIS Observations and MHD Simulations
NASA Astrophysics Data System (ADS)
Zhou, Xiaoyan
2015-04-01
Interplanetary shocks can cause magnetospheric disturbances on various scales including kinetic and MHD processes. In this paper we study a shock event using ARTEMIS in situ observations and OpenGGCM MHD simulations, which shows how significant effect of interplanetary shocks could be on the magnetotail. The two ARTEMIS spacecraft were located near the tail current sheet and lobe center at (-60, 1, -5Re_GSM) when the shock arrived and recorded an abrupt tail compression leading to significant enhancements in the plasma density, temperature, magnetic field strength, and cross-tail current density, as well as to tailward flows and current sheet crossings. About 10 min later, the spacecraft entered the sheath solar wind unexpectedly. Two hypotheses are considered: either the tail was cut off by the high solar wind ram pressure (~25-30 nPa), or the compressed tail was pushed aside by the appreciable dawnward solar wind flow imposed by the shock. OpenGGMC simulation results confirmed the second hypothesis and revealed that during this 10 min interval, the lobe center moved dawnward by ~12 Re and the tail width in Y was reduced from ~40 to 26 Re, which eventually exposed ARTEMIS to the sheath solar wind. Comparisons of plasma and magnetic parameters between ARTEMIS in situ observations and simulations showed a satisfied consistence.
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].
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.
Global MHD simulations of plasmaspheric plumes
NASA Astrophysics Data System (ADS)
Lyon, J.; Ouellette, J.; Merkin, V. G.
2015-12-01
The plasmasphere represents a separate population from the rest of themagnetosphere, generally high density but cold. When the solar windturns strongly southward this plasma is convected toward the daysidemagnetopause and affects the interaction of the solar wind with themagnetosphere. We have used multi-fluid simulations using the LFMglobal MHD code to model this interaction. The plasmasphere isinitialized as a cold (~1eV) hydrogen plasma in a quiet northward IMFstate with a density distribution appropriate for K_p = 1. Thecorotation potential from the ionosphere spins up the plasmasphereinto rough corotation. After a initialization period of hours, asouthward IMF is introduced and the enhanced convection initiates asurge of plasmaspheric density to the dayside. We discuss two aspectsof this interaction, the effects on dayside reconnection and on theKelvin-Helmholtz instability (KHI). We find that the mass loading ofmagnetospheric flux tubes slows local reconnection rates, though notas much as predicted by Borovsky et al. [2013]. We findthat the total reconnection rate is reduced, although not as much aswould be predicted by just the sub-solar reconnection rate. The KHIis somewhat reduced by the plasmaspheric loading of density in the lowlatitude boundary layer. It has been suggested that the presence ofthe plasmasphere may lead to enhanced ULF wave power in the interiorof the magnetosphere from the KHI waves. We find only a minimal effect during northward IMF. For southward IMF, the situation is complicated by the interaction of KHI with non-steady reconnection.
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)
Griton, Léa; Pantellini, Filippo; Moncuquet, Michel
2016-04-01
We present 3D simulations of the interaction of the solar wind with Mercury's magnetosphere using the magnetohydrodynamic code AMRVAC. A procedure for the identification of standing MHD modes has been applied to these simulations showing that large scale standing slow mode structures may exist in Mercury's magnetosheath. The identification is mostly based on relatively simple approximate analytical solutions to the old problem of determining the family of all standing linear plane MHD waves in a flowing plasma. The question of the identification of standing slow mode structures using in situ measurements such as the future BepiColombo MMO mission to Mercury will be discussed as well.
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.
MHD simulations of supernova driven ISM turbulence
NASA Astrophysics Data System (ADS)
Gressel, Oliver; Ziegler, Udo
The dynamic evolution of the (stratified) turbulent interstellar medium (ISM) is simulated utilizing a three-dimensional MHD model including various physical effects. The computational domain covers a box of 0.5x0.5x2.0 kpc at a resolution of typically 128x128x1024 grid cells. The model includes (constant kinematic) viscosity and magnetic diffusivity. The adiabatic equation of state is supplemented by a parameterized heating- and cooling-function allowing for thermal instability (TI). The update due to heating and cooling is implemented implicitly using a Patankar-type discretization. Turbulence is driven by supernova explosions which are modelled as local injections of thermal energy, smeared over three standard-deviations of a Gaussian support with FWHM of 20pc. Supernova rates are adopted for typical cited values. Within our model we make a distinction between Type I and Type II SNe. Latter are statistically clustered by the (artificial) constraint that the density at the explosion site be above average (with respect to a horizontal slab) - former are spatially uncorrelated. The dual-energy feature of the conservative NIRVANA-code is used to tackle the extreme ratio of kinetic to internal energy that arises from the violent energy input. We stress the importance of using a conservative scheme to properly transfer the injected energy to kinetic motion. The model also includes a differentially rotating background (with shearing boundary conditions in radial direction) as well as vertical stratification. The initial density and pressure profiles are in hydrostatic equilibrium with respect to the equation of state given by the radiative equilibrium. Including z-dependent heating rates this leads to a considerable deviation from usual isothermal initial models. The primary focus of this work is on the galactic dynamo and the generation of large-scale magnetic fields. As a secondary target we are also interested in general properties of the ISM that are of importance
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.
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.
Three Dimensional Simulations of Compressible Hall MHD Plasmas
Shaikh, Dastgeer; Shukla, P. K.
2008-10-15
We have developed three dimensional, time dependent, compressible, non-adiabatic, driven and massively parallelized Hall magnetohydrodynamic (MHD) simulations to investigate turbulent spectral cascades in a regime where characteristic lengthscales associated with plasma fluctuations are smaller than ion gyro radii. Such regime is ubiquitously present in solar wind and many other collisionless space plasmas. Particularly in the solar wind, the high time resolution databases identify a spectral break at the end of MHD inertial range spectrum that corresponds to a high frequency regime. In the regime, turbulent cascades cannot be explained by the usual MHD models. With the help of our 3D Hall MHD code, we find that characteristic turbulent interactions in the high frequency regime evolve typically on kinetic Alfven time scales. The turbulent fluctuation associated with kinetic Alfven interactions are compressive and anisotropic and possess equipartition of kinetic and magnetic energies.
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).
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.
Simulating MEMS Chevron Actuator for Strain Engineering 2D Materials
NASA Astrophysics Data System (ADS)
Vutukuru, Mounika; Christopher, Jason; Bishop, David; Swan, Anna
2D materials pose an exciting paradigm shift in the world of electronics. These crystalline materials have demonstrated high electric and thermal conductivities and tensile strength, showing great potential as the new building blocks of basic electronic circuits. However, strain engineering 2D materials for novel devices remains a difficult experimental feat. We propose the integration of 2D materials with MEMS devices to investigate the strain dependence on material properties such as electrical and thermal conductivity, refractive index, mechanical elasticity, and band gap. MEMS Chevron actuators, provides the most accessible framework to study strain in 2D materials due to their high output force displacements for low input power. Here, we simulate Chevron actuators on COMSOL to optimize actuator design parameters and accurately capture the behavior of the devices while under the external force of a 2D material. Through stationary state analysis, we analyze the response of the device through IV characteristics, displacement and temperature curves. We conclude that the simulation precisely models the real-world device through experimental confirmation, proving that the integration of 2D materials with MEMS is a viable option for constructing novel strain engineered devices. The authors acknowledge support from NSF DMR1411008.
NASA Astrophysics Data System (ADS)
Daldorff, L. K. S.; Toth, G.; Borovikov, D.; Gombosi, T. I.; Lapenta, G.
2014-12-01
With the new modeling capability in the Space Weather Modeling Framework (SWMF) of embedding an implicit Particle-in-Cell (PIC) model iPIC3D into the BATS-R-US magnetohydrodynamics model (Daldorff et al. 2014, JCP, 268, 236) we are ready to locally handle the full physics of the reconnection and its implications on the full system where globally, away from the reconnection region, a magnetohydrodynamic description is satisfactory. As magnetic reconnection is one of the main drivers in magnetospheric and heliospheric plasma dynamics, the self-consistent description of the electron dynamics in the coupled MHD-EPIC model is well suited for investigating the nature of these systems. We will compare the new embedded MHD-EPIC model with pure MHD and Hall MHD simulations of the Earth's magnetosphere.
Global and Kinetic MHD Simulation by the Gpic-MHD Code
NASA Astrophysics Data System (ADS)
Hiroshi, Naitou; Yusuke, Yamada; Kenji, Kajiwara; Wei-li, Lee; Shinji, Tokuda; Masatoshi, Yagi
2011-10-01
In order to implement large-scale and high-beta tokamak simulation, a new algorithm of the electromagnetic gyrokinetic PIC (particle-in-cell) code was proposed and installed on the Gpic-MHD code [Gyrokinetic PIC code for magnetohydrodynamic (MHD) simulation]. In the new algorithm, the vorticity equation and the generalized Ohm's law along the magnetic field are derived from the basic equations of the gyrokinetic Vlasov, Poisson, and Ampere system and are used to describe the spatio-temporal evolution of the field quantities of the electrostatic potential varphi and the longitudinal component of the vector potential Az. The basic algorithm is equivalent to solving the reduced-MHD-type equations with kinetic corrections, in which MHD physics related to Alfven modes are well described. The estimation of perturbed electron pressure from particle dynamics is dominant, while the effects of other moments are negligible. Another advantage of the algorithm is that the longitudinal induced electric field, ETz = -∂Az/∂t, is explicitly estimated by the generalized Ohm's law and used in the equations of motion. Furthermore, the particle velocities along the magnetic field are used (vz-formulation) instead of generalized momentums (pz-formulation), hence there is no problem of ‘cancellation', which would otherwise appear when Az is estimated from the Ampere's law in the pz-formulation. The successful simulation of the collisionless internal kink mode by the new Gpic-MHD with realistic values of the large-scale and high-beta tokamaks revealed the usefulness of the new algorithm.
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.
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
A discrete simulation of 2-D fluid flow on TERASYS
Mullins, P.G.; Krolak, P.D.
1995-12-01
A discrete simulation of two-dimensional (2-D) fluid flow, on a recently designed novel architecture called TERASYS is presented. The simulation uses a cellular automaton approach, implemented in a new language called data-parallel bit C (dbC). A performance comparison between our implementation on TERASYS and an implementation on the Connection Machine is discussed. We comment briefly on the suitability of the TERASYS system for modeling fluid flow using cellular automata.
Initial simulation of MHD instabilites in a high speed plasma accelerator
NASA Astrophysics Data System (ADS)
Kim, Jin-Soo; Hughes, Tom; Thio, Francis
2005-10-01
High density, high Mach number plasma jets are under development for a variety of critical fusion applications. These applications include fueling, rotation driving, and disruption mitigation in magnetic fusion devices. They also include a range of innovative approaches to high energy density plasmas. FAR-TECH, Inc. has begun 3D MHD simulations using the LSP code [1] to examine such high speed plasma jets. An initial study to benchmark the code is currently underway. The blow-by instability will be simulated in a coaxial plasma accelerator using the 3D LSP code and compared with the 2D MACH2 code results. [1] LSP-Manual-MRC-ABQ-R-1942.pdf
Simulation of 2D Fields of Raindrop Size Distributions
NASA Astrophysics Data System (ADS)
Berne, A.; Schleiss, M.; Uijlenhoet, R.
2008-12-01
The raindrop size distribution (DSD hereafter) is of primary importance for quantitative applications of weather radar measurements. The radar reflectivity~Z (directly measured by radar) is related to the power backscattered by the ensemble of hydrometeors within the radar sampling volume. However, the rain rate~R (the flux of water to the surface) is the variable of interest for many applications (hydrology, weather forecasting, air traffic for example). Usually, radar reflectivity is converted into rain rate using a power law such as Z=aRb. The coefficients a and b of the Z-R relationship depend on the DSD. The variability of the DSD in space and time has to be taken into account to improve radar rain rate estimates. Therefore, the ability to generate a large number of 2D fields of DSD which are statistically homogeneous provides a very useful simulation framework that nicely complements experimental approaches based on DSD data, in order to investigate radar beam propagation through rain as well as radar retrieval techniques. The proposed approach is based on geostatistics for structural analysis and stochastic simulation. First, the DSD is assumed to follow a gamma distribution. Hence a 2D field of DSDs can be adequately described as a 2D field of a multivariate random function consisting of the three DSD parameters. Such fields are simulated by combining a Gaussian anamorphosis and a multivariate Gaussian random field simulation algorithm. Using the (cross-)variogram models fitted on data guaranties that the spatial structure of the simulated fields is consistent with the observed one. To assess its validity, the proposed method is applied to data collected during intense Mediterranean rainfall. As only time series are available, Taylor's hypothesis is assumed to convert time series in 1D range profile. Moreover, DSD fields are assumed to be isotropic so that the 1D structure can be used to simulate 2D fields. A large number of 2D fields of DSD parameters are
Hall MHD Simulations of Comet 67P/Churyumov-Gerasimenko
NASA Astrophysics Data System (ADS)
Shou, Y.; Combi, M. R.; Rubin, M.; Hansen, K. C.; Toth, G.; Gombosi, T. I.
2012-12-01
Comets have highly eccentric orbits and a wide range of gas production rates and thus they are ideal subjects to study the interaction between the solar wind and nonmagnetized bodies. Hansen et al. (2007, Space Sci. Rev. 128, 133) used a fluid-based MHD model and a semi-kinetic hybrid particle model to study the plasma environment of comet 67P/Churyumov-Gerasimenko (CG), the Rosetta mission target comet, at different heliocentric distances. They showed that for such a weak comet at a large heliocentric distance, the length scales of the cometosheath and the bow shock are comparable to or smaller than the ion gyroradius, which violates the underlying assumption for a valid fluid description of the plasma. As a result, the classical ideal MHD model is not able to always give physical results, while the hybrid model, which accounts for the kinetic effects of ions with both cometary and solar wind origin, is more reliable. However, hybrid models are computationally expensive and the results can be noisy. A compromise approach is Hall MHD [Toth et al., 2008], which includes the Hall term in the MHD equations and allows for the decoupling of the ion and electron fluids. We use a single ion species Hall MHD model to simulate the plasma environment of comet 67P/CG and compare the results with the two models mentioned above. We find that the Hall effect is capable of reproducing some features of the hybrid model and thus extends the applicability of MHD. In addition, this study helps to identify the conditions and regions in the cometary plasma where the Hall effect is not negligible. This work is supported by NSF Planetary Astronomy grant AST0707283 and JPL subcontract 1266313 under NASA grant NMO710889.
Evidence of Active MHD Instability in EULAG-MHD Simulations of Solar Convection
NASA Astrophysics Data System (ADS)
Lawson, Nicolas; Strugarek, Antoine; Charbonneau, Paul
2015-11-01
We investigate the possible development of magnetohydrodynamical instabilities in the EULAG-MHD “millennium simulation” of Passos & Charbonneau. This simulation sustains a large-scale magnetic cycle characterized by solar-like polarity reversals taking place on a regular multidecadal cadence, and in which zonally oriented bands of strong magnetic fields accumulate below the convective layers, in response to turbulent pumping from above in successive magnetic half-cycles. Key aspects of this simulation include low numerical dissipation and a strongly sub-adiabatic fluid layer underlying the convectively unstable layers corresponding to the modeled solar convection zone. These properties are conducive to the growth and development of two-dimensional instabilities that are otherwise suppressed by stronger dissipation. We find evidence for the action of a non-axisymmetric magnetoshear instability operating in the upper portions of the stably stratified fluid layers. We also investigate the possibility that the Tayler instability may be contributing to the destabilization of the large-scale axisymmetric magnetic component at high latitudes. On the basis of our analyses, we propose a global dynamo scenario whereby the magnetic cycle is driven primarily by turbulent dynamo action in the convecting layers, but MHD instabilities accelerate the dissipation of the magnetic field pumped down into the overshoot and stable layers, thus perhaps significantly influencing the magnetic cycle period. Support for this scenario is found in the distinct global dynamo behaviors observed in an otherwise identical EULAG-MHD simulations, using a different degree of sub-adiabaticity in the stable fluid layers underlying the convection zone.
Simulation of bootstrap current in 2D and 3D ideal magnetic fields in tokamaks
NASA Astrophysics Data System (ADS)
Raghunathan, M.; Graves, J. P.; Cooper, W. A.; Pedro, M.; Sauter, O.
2016-09-01
We aim to simulate the bootstrap current for a MAST-like spherical tokamak using two approaches for magnetic equilibria including externally caused 3D effects such as resonant magnetic perturbations (RMPs), the effect of toroidal ripple, and intrinsic 3D effects such as non-resonant internal kink modes. The first approach relies on known neoclassical coefficients in ideal MHD equilibria, using the Sauter (Sauter et al 1999 Phys. Plasmas 6 2834) expression valid for all collisionalities in axisymmetry, and the second approach being the quasi-analytic Shaing–Callen (Shaing and Callen 1983 Phys. Fluids 26 3315) model in the collisionless regime for 3D. Using the ideal free-boundary magnetohydrodynamic code VMEC, we compute the flux-surface averaged bootstrap current density, with the Sauter and Shaing–Callen expressions for 2D and 3D ideal MHD equilibria including an edge pressure barrier with the application of resonant magnetic perturbations, and equilibria possessing a saturated non-resonant 1/1 internal kink mode with a weak internal pressure barrier. We compare the applicability of the self-consistent iterative model on the 3D applications and discuss the limitations and advantages of each bootstrap current model for each type of equilibrium.
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.
Coupling MHD Simulations of CMEs to SEP Models
NASA Astrophysics Data System (ADS)
Torok, T.; Gorby, M.; Linker, J.; Schwadron, N.
2015-12-01
Large Solar Energetic Particle events (SEPs) are a main space weather hazard and extremely dangerous to astronauts and electronic equipmentin space. They are typically associated with fast Coronal Mass Ejections (CMEs). Recent results indicate that SEPs can be generated already inthe early phase of CME expansion low in the corona, but the underlyingphysical mechanisms are not yet well understood. State-of-the-artmagnetohydrodynamic (MHD) simulations of CME initiation and evolution,combined with numerical models of particle acceleration and propagation,provide a powerful tool to investigate these mechanisms. In this talk, we present recent developments in the coupling of CORHEL/MAS thermodynamicMHD simulations of fast CMEs to the EPREM particle code, and we discuss the insights that can be gained from such a combined modeling approach.
Quick Time-dependent Ionization Calculations Depending on MHD Simulations
NASA Astrophysics Data System (ADS)
Shen, Chengcai; Raymond, John C.; Murphy, Nicholas Arnold
2014-06-01
Time-dependent ionization is important in astrophysical environments where the thermodynamic time scale is shorter than ionization time scale. In this work, we report a FORTRAN program that performs fast non-equilibrium ionization calculations based on parallel computing. Using MHD simulation results, we trace the movements of plasma in a Lagrangian framework, and obtain evolutionary history of temperature and electron density. Then the time-dependent ionization equations are solved using the eigenvalue method. For any complex temperature and density histories, we introduce a advanced time-step strategy to improve the computational efficiency. Our tests show that this program has advantages of high numerical stability and high accuracy. In addition, it is also easy to integrate this solver with the other MHD routines.
Spectral Methods in General Relativistic MHD Simulations
NASA Astrophysics Data System (ADS)
Garrison, David
2012-03-01
In this talk I discuss the use of spectral methods in improving the accuracy of a General Relativistic Magnetohydrodynamic (GRMHD) computer code. I introduce SpecCosmo, a GRMHD code developed as a Cactus arrangement at UHCL, and show simulation results using both Fourier spectral methods and finite differencing. This work demonstrates the use of spectral methods with the FFTW 3.3 Fast Fourier Transform package integrated with the Cactus Framework to perform spectral differencing using MPI.
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.
Numerical simulation of rock cutting using 2D AUTODYN
NASA Astrophysics Data System (ADS)
Woldemichael, D. E.; Rani, A. M. Abdul; Lemma, T. A.; Altaf, K.
2015-12-01
In a drilling process for oil and gas exploration, understanding of the interaction between the cutting tool and the rock is important for optimization of the drilling process using polycrystalline diamond compact (PDC) cutters. In this study the finite element method in ANSYS AUTODYN-2D is used to simulate the dynamics of cutter rock interaction, rock failure, and fragmentation. A two-dimensional single PDC cutter and rock model were used to simulate the orthogonal cutting process and to investigate the effect of different parameters such as depth of cut, and back rake angle on two types of rocks (sandstone and limestone). In the simulation, the cutting tool was dragged against stationary rock at predetermined linear velocity and the depth of cut (1,2, and 3 mm) and the back rake angles(-10°, 0°, and +10°) were varied. The simulation result shows that the +10° back rake angle results in higher rate of penetration (ROP). Increasing depth of cut leads to higher ROP at the cost of higher cutting force.
3D MHD Simulations of Spheromak Compression
NASA Astrophysics Data System (ADS)
Stuber, James E.; Woodruff, Simon; O'Bryan, John; Romero-Talamas, Carlos A.; Darpa Spheromak Team
2015-11-01
The adiabatic compression of compact tori could lead to a compact and hence low cost fusion energy system. The critical scientific issues in spheromak compression relate both to confinement properties and to the stability of the configuration undergoing compression. We present results from the NIMROD code modified with the addition of magnetic field coils that allow us to examine the role of rotation on the stability and confinement of the spheromak (extending prior work for the FRC). We present results from a scan in initial rotation, from 0 to 100km/s. We show that strong rotational shear (10km/s over 1cm) occurs. We compare the simulation results with analytic scaling relations for adiabatic compression. Work performed under DARPA grant N66001-14-1-4044.
MHD Simulations of Plasma Dynamics with Non-Axisymmetric Boundaries
NASA Astrophysics Data System (ADS)
Hansen, Chris; Levesque, Jeffrey; Morgan, Kyle; Jarboe, Thomas
2015-11-01
The arbitrary geometry, 3D extended MHD code PSI-TET is applied to linear and non-linear simulations of MCF plasmas with non-axisymmetric boundaries. Progress and results from simulations on two experiments will be presented: 1) Detailed validation studies of the HIT-SI experiment with self-consistent modeling of plasma dynamics in the helicity injectors. Results will be compared to experimental data and NIMROD simulations that model the effect of the helicity injectors through boundary conditions on an axisymmetric domain. 2) Linear studies of HBT-EP with different wall configurations focusing on toroidal asymmetries in the adjustable conducting wall. HBT-EP studies the effect of active/passive stabilization with an adjustable ferritic wall. Results from linear verification and benchmark studies of ideal mode growth with and without toroidal asymmetries will be presented and compared to DCON predictions. Simulations of detailed experimental geometries are enabled by use of the PSI-TET code, which employs a high order finite element method on unstructured tetrahedral grids that are generated directly from CAD models. Further development of PSI-TET will also be presented including work to support resistive wall regions within extended MHD simulations. Work supported by DoE.
Relativistic MHD simulations of extragalactic jets
NASA Astrophysics Data System (ADS)
Leismann, T.; Antón, L.; Aloy, M. A.; Müller, E.; Martí, J. M.; Miralles, J. A.; Ibáñez, J. M.
2005-06-01
We have performed a comprehensive parameter study of the morphology and dynamics of axisymmetric, magnetized, relativistic jets by means of numerical simulations. The simulations have been performed with an upgraded version of the GENESIS code which is based on a second-order accurate finite volume method involving an approximate Riemann solver suitable for relativistic ideal magnetohydrodynamic flows, and a method of lines. Starting from pure hydrodynamic models we consider the effect of a magnetic field of increasing strength (up to β ≡ |b|2/2p ≈ 3.3 times the equipartition value) and different topology (purely toroidal or poloidal). We computed several series of models investigating the dependence of the dynamics on the magnetic field in jets of different beam Lorentz factor and adiabatic index. We find that the inclusion of the magnetic field leads to diverse effects which contrary to Newtonian magnetohydrodynamics models do not always scale linearly with the (relative) strength of the magnetic field. The relativistic models show, however, some clear trends. Axisymmetric jets with toroidal magnetic fields produce a cavity which consists of two parts: an inner one surrounding the beam which is compressed by magnetic forces, and an adjacent outer part which is inflated due to the action of the magnetic field. The outer border of the outer part of the cavity is given by the bow-shock where its interaction with the external medium takes place. Toroidal magnetic fields well below equipartition (β = 0.05) combined with a value of the adiabatic index of 4/3 yield extremely smooth jet cavities and stable beams. Prominent nose cones form when jets are confined by toroidal fields and carry a high Poynting flux (σ≡ |b|2/ρ>0.01 and β≥ 1). In contrast, none of our models possessing a poloidal field develops such a nose cone. The size of the nose cone is correlated with the propagation speed of the Mach disc (the smaller the speed the larger is the size). If two
Multiscale simulation of 2D elastic wave propagation
NASA Astrophysics Data System (ADS)
Zhang, Wensheng; Zheng, Hui
2016-06-01
In this paper, we develop the multiscale method for simulation of elastic wave propagation. Based on the first-order velocity-stress hyperbolic form of 2D elastic wave equation, the particle velocities are solved first ona coarse grid by the finite volume method. Then the stress tensor is solved by using the multiscale basis functions which can represent the fine-scale variation of the wavefield on the coarse grid. The basis functions are computed by solving a local problem with the finite element method. The theoretical formulae and description of the multiscale method for elastic wave equation are given in more detail. The numerical computations for an inhomogeneous model with random scatter are completed. The results show the effectiveness of the multiscale method.
Analysis and gyrokinetic simulation of MHD Alfven wave interactions
NASA Astrophysics Data System (ADS)
Nielson, Kevin Derek
The study of low-frequency turbulence in magnetized plasmas is a difficult problem due to both the enormous range of scales involved and the variety of physics encompassed over this range. Much of the progress that has been made in turbulence theory is based upon a result from incompressible magnetohydrodynamics (MHD), in which energy is only transferred from large scales to small via the collision of Alfven waves propagating oppositely along the mean magnetic field. Improvements in laboratory devices and satellite measurements have demonstrated that, while theories based on this premise are useful over inertial ranges, describing turbulence at scales that approach particle gyroscales requires new theory. In this thesis, we examine the limits of incompressible MHD theory in describing collisions between pairs of Alfven waves. This interaction represents the fundamental unit of plasma turbulence. To study this interaction, we develop an analytic theory describing the nonlinear evolution of interacting Alfven waves and compare this theory to simulations performed using the gyrokinetic code AstroGK. Gyrokinetics captures a much richer set of physics than that described by incompressible MHD, and is well-suited to describing Alfvenic turbulence around the ion gyroscale. We demonstrate that AstroGK is well suited to the study of physical Alfven waves by reproducing laboratory Alfven dispersion data collected using the LAPD. Additionally, we have developed an initialization alogrithm for use with AstroGK that allows exact Alfven eigenmodes to be initialized with user specified amplitudes and phases. We demonstrate that our analytic theory based upon incompressible MHD gives excellent agreement with gyrokinetic simulations for weakly turbulent collisions in the limit that k⊥rho i << 1. In this limit, agreement is observed in the time evolution of nonlinear products, and in the strength of nonlinear interaction with respect to polarization and scale. We also examine the
Quantum Simulation with 2D Arrays of Trapped Ions
NASA Astrophysics Data System (ADS)
Richerme, Philip
2016-05-01
The computational difficulty of solving fully quantum many-body spin problems is a significant obstacle to understanding the behavior of strongly correlated quantum matter. This work proposes the design and construction of a 2D quantum spin simulator to investigate the physics of frustrated materials, highly entangled states, mechanisms potentially underpinning high-temperature superconductivity, and other topics inaccessible to current 1D systems. The effective quantum spins will be encoded within the well-isolated electronic levels of trapped ions, confined in a two-dimensional planar geometry, and made to interact using phonon-mediated optical dipole forces. The system will be scalable to 100+ quantum particles, far beyond the realm of classical intractability, while maintaining individual-ion control, long quantum coherence times, and site-resolved projective spin measurements. Once constructed, the two-dimensional quantum simulator will implement a broad range of spin models on a variety of reconfigurable lattices and characterize their behavior through measurements of spin-spin correlations and entanglement. This versatile tool will serve as an important experimental resource for exploring difficult quantum many-body problems in a regime where classical methods fail.
MHD simulation studies of z-pinch shear flow stabilization
NASA Astrophysics Data System (ADS)
Paraschiv, I.; Bauer, B. S.; Sotnikov, V. I.; Makhin, V.; Siemon, R. E.
2003-10-01
The development of the m=0 instability in a z-pinch in the presence of sheared plasma flows is investigated with the aid of a two-dimensional magnetohydrodynamic (MHD) simulation code (MHRDR). The linear growth rates are compared to the results obtained by solving the ideal MHD linearized equations [1] and to the results obtained using a 3D hybrid simulation code [2]. The instability development is followed into the nonlinear regime where its growth and saturation are examined. [1] V.I. Sotnikov, I. Paraschiv, V. Makhin, B.S. Bauer, J.-N. Leboeuf, and J.M. Dawson, "Linear analysis of sheared flow stabilization of global magnetohydrodynamic instabilities based on the Hall fluid mode", Phys. Plasmas 9, 913 (2002). [2] V.I. Sotnikov, V. Makhin, B.S. Bauer, P. Hellinger, P. Travnicek, V. Fiala, J.-N. Leboeuf, "Hybrid Simulations of Current-Carrying Instabilities in Z-pinch Plasmas with Sheared Axial Flow", AIP Conference Proceedings, Volume 651, Dense Z-Pinches: 5th International Conference on Dense Z-Pinches, edited by J. Davis et al., page 396, June 2002.
Relativistic MHD simulations of stellar core collapse and magnetars
NASA Astrophysics Data System (ADS)
Font, José A.; Cerdá-Durán, Pablo; Gabler, Michael; Müller, Ewald; Stergioulas, Nikolaos
2011-02-01
We present results from simulations of magneto-rotational stellar core collapse along with Alfvén oscillations in magnetars. These simulations are performed with the CoCoA/CoCoNuT code, which is able to handle ideal MHD flows in dynamical spacetimes in general relativity. Our core collapse simulations highlight the importance of genuine magnetic effects, like the magneto-rotational instability, for the dynamics of the flow. For the modelling of magnetars we use the anelastic approximation to general relativistic MHD, which allows for an effective suppression of fluid modes and an accurate description of Alfvén waves. We further compute Alfvén oscillation frequencies along individual magnetic field lines with a semi-analytic approach. Our work confirms previous results based on perturbative approaches regarding the existence of two families of quasi-periodic oscillations (QPOs), with harmonics at integer multiples of the fundamental frequency. Additional material is presented in the accompanying contribution by Gabler et al (2010b) in these proceedings.
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.
Understanding Accretion Disks through Three Dimensional Radiation MHD Simulations
NASA Astrophysics Data System (ADS)
Jiang, Yan-Fei
I study the structures and thermal properties of black hole accretion disks in the radiation pressure dominated regime. Angular momentum transfer in the disk is provided by the turbulence generated by the magneto-rotational instability (MRI), which is calculated self-consistently with a recently developed 3D radiation magneto-hydrodynamics (MHD) code based on Athena. This code, developed by my collaborators and myself, couples both the radiation momentum and energy source terms with the ideal MHD equations by modifying the standard Godunov method to handle the stiff radiation source terms. We solve the two momentum equations of the radiation transfer equations with a variable Eddington tensor (VET), which is calculated with a time independent short characteristic module. This code is well tested and accurate in both optically thin and optically thick regimes. It is also accurate for both radiation pressure and gas pressure dominated flows. With this code, I find that when photon viscosity becomes significant, the ratio between Maxwell stress and Reynolds stress from the MRI turbulence can increase significantly with radiation pressure. The thermal instability of the radiation pressure dominated disk is then studied with vertically stratified shearing box simulations. Unlike the previous results claiming that the radiation pressure dominated disk with MRI turbulence can reach a steady state without showing any unstable behavior, I find that the radiation pressure dominated disks always either collapse or expand until we have to stop the simulations. During the thermal runaway, the heating and cooling rates from the simulations are consistent with the general criterion of thermal instability. However, details of the thermal runaway are different from the predictions of the standard alpha disk model, as many assumptions in that model are not satisfied in the simulations. We also identify the key reasons why previous simulations do not find the instability. The thermal
Global MHD Simulation of Mesoscale Structures at the Magnetospheric Boundary
NASA Technical Reports Server (NTRS)
Berchem, Jean
1998-01-01
The research carried out for this protocol was focused on the study of mesoscales structures at the magnetospheric boundary. We investigated three areas: (1) the structure of the magnetospheric boundary for steady solar wind conditions; (2) the dynamics of the dayside magnetospheric boundary and (3) the dynamics of the distant tail magnetospheric boundary. Our approach was to use high resolution three-dimensional global magnetohydrodynamic (MHD) simulations of the interaction of the solar wind with the Earth's magnetosphere. We first considered simple variations of the interplanetary conditions to obtain generic cases that helped us in establishing the basic cause and effect relationships for steady solar wind conditions. Subsequently, we used actual solar wind plasma and magnetic field parameters measured by an upstream spacecraft as input to the simulations and compared the simulation results with sequences of events observed by another or several other spacecraft located downstream the bow shock. In particular we compared results with observations made when spacecraft crossed the magnetospheric boundary.
A 2D simulation model for urban flood management
NASA Astrophysics Data System (ADS)
Price, Roland; van der Wielen, Jonathan; Velickov, Slavco; Galvao, Diogo
2014-05-01
The European Floods Directive, which came into force on 26 November 2007, requires member states to assess all their water courses and coast lines for risk of flooding, to map flood extents and assets and humans at risk, and to take adequate and coordinated measures to reduce the flood risk in consultation with the public. Flood Risk Management Plans are to be in place by 2015. There are a number of reasons for the promotion of this Directive, not least because there has been much urban and other infrastructural development in flood plains, which puts many at risk of flooding along with vital societal assets. In addition there is growing awareness that the changing climate appears to be inducing more frequent extremes of rainfall with a consequent increases in the frequency of flooding. Thirdly, the growing urban populations in Europe, and especially in the developing countries, means that more people are being put at risk from a greater frequency of urban flooding in particular. There are urgent needs therefore to assess flood risk accurately and consistently, to reduce this risk where it is important to do so or where the benefit is greater than the damage cost, to improve flood forecasting and warning, to provide where necessary (and possible) flood insurance cover, and to involve all stakeholders in decision making affecting flood protection and flood risk management plans. Key data for assessing risk are water levels achieved or forecasted during a flood. Such levels should of course be monitored, but they also need to be predicted, whether for design or simulation. A 2D simulation model (PriceXD) solving the shallow water wave equations is presented specifically for determining flood risk, assessing flood defense schemes and generating flood forecasts and warnings. The simulation model is required to have a number of important properties: -Solve the full shallow water wave equations using a range of possible solutions; -Automatically adjust the time step and
MHD simulations of ram pressure stripping of a disk galaxy
NASA Astrophysics Data System (ADS)
Ramos, Mariana; Gomez, Gilberto
2015-08-01
The removal of the ISM of disk galaxies through ram pressure stripping (RPS) has been extensively studied in numerous simulations. These models show that this process has a significant impact on galaxy evolution (the truncation of the ISM will lead to a decrease in the star formation and the galaxy will become redder).Nevertheless, the role of the magnetic fields (MFs) on the dynamics of the gas in this process has been hardly studied, although the influence of magnetic fields on the large scale disk structure is well established. The presence of MFs produce a less compressible gas, thus increasing the scale height of the gas in the galaxy, that is, gas can be found farther away from the galactic potential well, which may lead to an easier removal of gas. We test this idea by performing a 3D MHD simulation of a disk galaxy that experiences RPS under the wind-tunnel approximation.
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.
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.
Final Report: "Large-Eddy Simulation of Anisotropic MHD Turbulence"
Zikanov, Oleg
2008-06-23
To acquire better understanding of turbulence in flows of liquid metals and other electrically conducting fluids in the presence of steady magnetic fields and to develop an accurate and physically adequate LES (large-eddy simulation) model for such flows. The scientific objectives formulated in the project proposal have been fully completed. Several new directions were initiated and advanced in the course of work. Particular achievements include a detailed study of transformation of turbulence caused by the imposed magnetic field, development of an LES model that accurately reproduces this transformation, and solution of several fundamental questions of the interaction between the magnetic field and fluid flows. Eight papers have been published in respected peer-reviewed journals, with two more papers currently undergoing review, and one in preparation for submission. A post-doctoral researcher and a graduate student have been trained in the areas of MHD, turbulence research, and computational methods. Close collaboration ties have been established with the MHD research centers in Germany and Belgium.
Complexities of a 3-D plasmoid flux rope as shown by an MHD simulation
NASA Astrophysics Data System (ADS)
Farr, N. L.; Baker, D. N.; Wiltberger, M.
2008-12-01
The results of a global magnetohydrodynamic (MHD) simulation of a pair of magnetospheric substorms on 11 August 2002 are presented. Comparisons of data with simulation results reveal a good agreement regarding the sequence of events during substorm development. We give particular emphasis to results in the simulation of a flux rope formed during the second substorm. Unlike standard 2-D depictions of reconnection and plasmoid release during the substorm sequence, the simulation shows a highly complex structure that has considerable winding of both closed and open field lines. Additionally, the simulated flux rope does not move tailward uniformly, but rather it has asymmetric motion in which the dawn flank portion moves tailward prior to the dusk portion of the flux rope. This results in a skewed flux rope structure that runs almost parallel to the tail axis instead of perpendicular to it. The simulation compares well with both prior flux rope simulations as well as satellite observations of flux ropes. We use the global simulation to map flux tube properties to the ionosphere, which allows the complexity of the mapping of the magnetic field structure from the tail to the ionosphere to be seen in a novel manner.
A 2D simulation model for urban flood management
NASA Astrophysics Data System (ADS)
Price, Roland; van der Wielen, Jonathan; Velickov, Slavco; Galvao, Diogo
2014-05-01
The European Floods Directive, which came into force on 26 November 2007, requires member states to assess all their water courses and coast lines for risk of flooding, to map flood extents and assets and humans at risk, and to take adequate and coordinated measures to reduce the flood risk in consultation with the public. Flood Risk Management Plans are to be in place by 2015. There are a number of reasons for the promotion of this Directive, not least because there has been much urban and other infrastructural development in flood plains, which puts many at risk of flooding along with vital societal assets. In addition there is growing awareness that the changing climate appears to be inducing more frequent extremes of rainfall with a consequent increases in the frequency of flooding. Thirdly, the growing urban populations in Europe, and especially in the developing countries, means that more people are being put at risk from a greater frequency of urban flooding in particular. There are urgent needs therefore to assess flood risk accurately and consistently, to reduce this risk where it is important to do so or where the benefit is greater than the damage cost, to improve flood forecasting and warning, to provide where necessary (and possible) flood insurance cover, and to involve all stakeholders in decision making affecting flood protection and flood risk management plans. Key data for assessing risk are water levels achieved or forecasted during a flood. Such levels should of course be monitored, but they also need to be predicted, whether for design or simulation. A 2D simulation model (PriceXD) solving the shallow water wave equations is presented specifically for determining flood risk, assessing flood defense schemes and generating flood forecasts and warnings. The simulation model is required to have a number of important properties: -Solve the full shallow water wave equations using a range of possible solutions; -Automatically adjust the time step and
MHD Simulation of the Inverse Pinch Plasma Discharge
Esaulov, A; Bauer, B; Lindemuth, I; Makhin, V; Presura, R; Ryutov, D
2004-07-01
A wall confined plasma in an inverse pinch configuration holds potential as a plasma target for Magnetized Target Fusion (MTF) as well as the simple geometry to study wall-confined plasma. An experiment is planned to study the inverse pinch configuration using the Nevada Terawatt Facility (NTF) at the University of Nevada, Reno (UNR). The dynamics of the discharge formation have been analyzed using analytic models and numerical methods. Strong heating occurs by thermalization of directed energy when an outward moving current sheet (the inverse pinch effect) collides with the outer wall of the experimental chamber. Two dimensional MHD simulations show Rayleigh-Taylor and Richtmyer-Meshkov -like modes of instability, as expected because of the shock acceleration during plasma formation phase. The instabilities are not disruptive, but give rise to a mild level of turbulence. The conclusion from this work is that an interesting experiment relevant to wall confinement for MTF could be done using existing equipment at UNR.
Relative timing of substorm features in MHD simulations
NASA Technical Reports Server (NTRS)
Hesse, Michael; Birn, Joachim
1992-01-01
An investigation of the temporal sequence of substorm phenomena based on three dimensional MHD (magnetohydrodynamic) simulations of magnetic reconnection and plasmoid formation is presented. The investigation utilizes a spatially localized resistivity model which leads to a significantly faster evolution than found in previous investigations. The analysis of the results concentrates on substorm features that have received considerable attention in the past. The formation of magnetic neutral lines, the occurrence of fast flows directed both earthward and tailward, and the magnetic field changes leading to the formation of the substorm current wedge, and to the depolarization of the magnetic field earthward of the reconnection region and its dependence on the spatial distribution of resistivity, are discussed. These phenomena are seen as an integral part of the nonlinear evolution of the three dimensional tearing instability.
NASA Astrophysics Data System (ADS)
Ju, Wenhua; Stone, James M.; Zhu, Zhaohuan
2016-06-01
We present results from the first global 3D MHD simulations of accretion disks in cataclysmic variable (CV) systems in order to investigate the relative importance of angular momentum transport via turbulence driven by the magnetorotational instability (MRI) compared with that driven by spiral shock waves. Remarkably, we find that even with vigorous MRI turbulence, spiral shocks are an important component of the overall angular momentum budget, at least when temperatures in the disk are high (so that Mach numbers are low). In order to understand the excitation, propagation, and damping of spiral density waves in our simulations more carefully, we perform a series of 2D global hydrodynamical simulations with various equation of states, both with and without mass inflow via the Lagrangian point (L1). Compared with previous similar studies, we find the following new results. (1) The linear wave dispersion relation fits the pitch angles of spiral density waves very well. (2) We demonstrate explicitly that mass accretion is driven by the deposition of negative angular momentum carried by the waves when they dissipate in shocks. (3) Using Reynolds stress scaled by gas pressure to represent the effective angular momentum transport rate {α }{eff} is not accurate when mass accretion is driven by non-axisymmetric shocks. (4) Using the mass accretion rate measured in our simulations to directly measure α defined in standard thin-disk theory, we find 0.02≲ {α }{eff}≲ 0.05 for CV disks, consistent with observed values in quiescent states of dwarf novae. In this regime, the disk may be too cool and neutral for the MRI to operate and spiral shocks are a possible accretion mechanism. However, we caution that our simulations use unrealistically low Mach numbers in this regime and, therefore, future models with more realistic thermodynamics and non-ideal MHD are warranted.
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
MHD Simulations of the Initiation of Coronal Mass Ejections
NASA Astrophysics Data System (ADS)
Fan, Yuhong; Chatterjee, Piyali
Using three-dimensional MHD simulations, we model the quasi-static evolution and the onset of eruption of twisted magnetic flux ropes in the solar corona. We present simulations where the eruption is triggered by either the onset of the torus instability or the helical kink instability of the line-tied coronal flux rope. The simulations show that S (or inverse S) shaped current sheets develop along topological structures identified as Quasi Separatrix Layers (QSLs), during the quasi-static phase before the eruption. Reconnections in the current sheets effectively add twisted flux to the flux rope and thus allow it to rise quasi-statically to the critical height for the onset of the torus instability. We examine the thermal features produced by the current sheet formation and the associated reconnections and found that they can explain some of the observed features in coronal prominence cavities as well as in pre-eruption active regions. We also present simulations of the development of a homologous sequence of CMEs caused by the repeated formation and partial eruption of kink unstable flux ropes as a result of continued flux emergence. It is found that such homologous CMEs tend to be cannibalistic, leading to the formation of more energetic, highly twisted ejecta.
MHD Simulations of the Plasma Flow in the Magnetic Nozzle
NASA Technical Reports Server (NTRS)
Smith, T. E. R.; Keidar, M.; Sankaran, K.; olzin, K. A.
2013-01-01
The magnetohydrodynamic (MHD) flow of plasma through a magnetic nozzle is simulated by solving the governing equations for the plasma flow in the presence of an static magnetic field representing the applied nozzle. This work will numerically investigate the flow and behavior of the plasma as the inlet plasma conditions and magnetic nozzle field strength are varied. The MHD simulations are useful for addressing issues such as plasma detachment and to can be used to gain insight into the physical processes present in plasma flows found in thrusters that use magnetic nozzles. In the model, the MHD equations for a plasma, with separate temperatures calculated for the electrons and ions, are integrated over a finite cell volume with flux through each face computed for each of the conserved variables (mass, momentum, magnetic flux, energy) [1]. Stokes theorem is used to convert the area integrals over the faces of each cell into line integrals around the boundaries of each face. The state of the plasma is described using models of the ionization level, ratio of specific heats, thermal conductivity, and plasma resistivity. Anisotropies in current conduction due to Hall effect are included, and the system is closed using a real-gas equation of state to describe the relationship between the plasma density, temperature, and pressure.A separate magnetostatic solver is used to calculate the applied magnetic field, which is assumed constant for these calculations. The total magnetic field is obtained through superposition of the solution for the applied magnetic field and the self-consistently computed induced magnetic fields that arise as the flowing plasma reacts to the presence of the applied field. A solution for the applied magnetic field is represented in Fig. 1 (from Ref. [2]), exhibiting the classic converging-diverging field pattern. Previous research was able to demonstrate effects such as back-emf at a super-Alfvenic flow, which significantly alters the shape of the
High Resolution Simulations of Relativistic Hydrodynamic and MHD Turbulence
NASA Astrophysics Data System (ADS)
Zrake, Jonathan; MacFadyen, A.
2013-01-01
We present a program of simulations designed to investigate the basic properties of relativistic hydrodynamic and magnetohydrodynamic (MHD) turbulence. We employ a well-tested 5th-order accurate numerical scheme at resolutions of up to 2048^3 zones for hydrodynamic turbulence, and a minimally diffusive 2nd-order scheme at resolutions of up to 1024^3 in the case of relativistic MHD. For the hydrodynamic case, we simulate a relativistically hot gas in a cubic periodic domain continuously driven at large scales with Lorentz factor of about 3. We find that relativistic turbulent velocity fluctuations with Γ β > 1 persist from the driving scale down to scales an order of magnitude smaller, demonstrating the existence of a sustained relativistic turbulent cascade. The power spectrum of the fluid 4-velocity is broadly Kolmogorov-like, roughly obeying a power law with 5/3 index between scales 1/10 and 1/100 of the domain. Departures from 5/3 scaling are larger for the power spectrum of 3-velocity. We find that throughout the inertial interval, 25% of power is in dilatational modes, which obey strict power law scaling between 1/2 and 1/100 of the domain with an index of 1.88. Our program also explores turbulent amplification of magnetic fields in the conditions of merging neutron stars, using a realistic equation of state for dense nuclear matter (ρ ˜ 10^13 g/cm^3). We find that very robustly, seed fields are amplified to magnetar strength (≥ 4 * 10^16 Gauss) within ˜1 micro-second for fluid volumes near the size of the NS crust thickness <10 meters. We present power spectra of the kinetic and magnetic energy taken long into the fully stationary evolution of the highest resolution models, finding the magnetic energy to be in super-equipartition (4 times larger) with the kinetic energy through the inertial range. We believe that current global simulations of merging NS binaries are insufficiently resolved for studying field amplification via turbulent processes
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.
Preliminary analysis of the dynamic heliosphere by MHD simulations
Washimi, H.; Zank, G. P.; Tanaka, T.
2006-09-26
A preliminary analysis of the dynamic heliosphere to estimate the termination shock (TS) distance from the sun around the time when Voyager 1 passed the termination shock at December 16, 2004 is performed by using MHD simulations. For input to this simulation, we use the Voyager 2 solar-wind data. We first find a stationary solution of the 3-D outer heliosphere by assigning a set of LISM parameters as our outer boundary conditions and then the dynamical analysis is performed. The model TS crossing is within 6 months of the observed date. The TS is pushed outward every time a high ram-pressure solar wind pulse arrives. After the end of the high ram-pressure wind, the TS shock shrinks inward. When the last Halloween event passed through the TS at DOY 250, 2004, the TS began to shrink inward very quickly and the TS crossed V1. The highest inward speed of the TS is over 400 km/s. The high ram-pressure solar wind transmitted through the TS becomes a high thermal-pressure plasma in the heliosheath, acting to push the TS inward. This suggests that the position of the TS is determined not only by the steady-state pressure balance condition between the solar wind ram-pressure and the LISM pressure, but by the dynamical ram pressure too. The period when the high ram-pressure solar wind arrives at the TS shock seems to correspond to the period of the TS particle event (Stone et al, 2005, Decker et al., 2005). The TS crossing date will be revised in future simulations using a more appropriate set of parameters for the LISM. This will enable us to undertake a detailed comparison of the simulation results with the TS particle events.
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.
Complexities of a 3-D flux rope as shown by MHD simulation
NASA Astrophysics Data System (ADS)
Farr, N.; Baker, D. N.; Wiltberger, M.
2007-12-01
This paper presents the results of a global magnetohydrodynamic (MHD) simulation of a pair of substorms on August 11, 2002. Comparisons of data with simulation results reveal an agreement regarding the sequence of events in the magnetosphere. We then present the results in the simulation of a flux rope formed during the second substorm. Unlike standard 2-D depictions of reconnection and plasmoid release during a substorm, the simulation shows a highly complex structure that has considerable winding of both closed and open field lines. Additionally the flux rope does not move tailward uniformly, but rather has a assymetric motion where the dawn flank moves tailward prior to the dusk end of the flux rope, resulting in a a skewed flux rope that runs almost downtail instead of crosstail. These features can add considerably complexity to satellites observing a flux rope structure in-situ. A single spacecraft could observe particle populations that go through a sequence of alternating open and closed field lines and spacecraft separated by small spatial distances could observe quite different populations as well.
Observations and MHD Simulations for a Shocked Magnetotail
NASA Astrophysics Data System (ADS)
Zhou, X.; Zhou, X. Z.; Angelopoulos, V.; Raeder, J.; Oliveira, D.; Shi, Q.
2014-12-01
Recent studies disclosed that interplanetary shocks not only raise global auroral activity, but also cause significant tail disturbances, ranging from ULF wave excitation to abrupt cross-tail current sheet thinning and current density increase, generation of burst bulk flows and dipolarization fronts, and to magnetic reconnection enhancement. In addition, shocks can also cause significant deformation of the magnetotail at ~60 Re and beyond. In this paper we study a shock event using ARTEMIS in situ observations and OpenGGCM MHD simulations. The two ARTEMIS spacecraft were located near the tail current sheet and lobe center at (-60, 1, -5Re_GSM) when the shock arrived and recorded an abrupt tail compression leading to significant enhancements in the plasma density, temperature, magnetic field strength, and cross-tail current density, as well as to tailward flows. However, ~10 min later, the spacecraft entered the sheath solar wind unexpectedly. Two hypotheses are considered: either the tail was cut off by the high solar wind ram pressure (~25-30 nPa), or the compressed tail was pushed aside by the appreciable Vy solar wind flow component imposed by the shock. OpenGGMC simulation results confirmed the second hypothesis and disclosed that for this event the magnetic pressure played a dominant role at X=-60 Re for the compression. In addition to the shock normal direction and shock compression, the anisotropic (transverse) magnetic pressure also contributed to the significant reduction of the lobe Y dimension. Therefore, during this 10 min interval, the lobe center moved dawnward by ~12 Re and the tail width in Y was reduced from 40 to 26 Re, which eventually exposed ARTEMIS to the sheath solar wind. Comparisons of plasma and magnetic parameters between ARTEMIS in situ observations and simulations showed a satisfied consistence.
MHD simulation of RF current drive in MST
Hendries, E. R.; Anderson, J. K.; Forest, C. B.; Reusch, J. A.; Seltzman, A. H.; Sovinec, C. R.; Diem, S.; Harvey, R. W.
2014-02-12
Auxiliary heating and current drive using RF waves such as the electron Bernstein wave (EBW) promises to advance the performance of the reversed field pinch (RFP). In previous computational work [1], a hypothetical edge-localized current drive is shown to suppress the tearing activity which governs the macroscopic transport properties of the RFP. The ideal conditions for tearing stabilization include a reduced toroidal induction, and precise width and radial position of the Gaussian-shaped external current drive. In support of the EBW experiment on the Madison Symmetric Torus, an integrated modeling scheme now incorporates ray tracing and Fokker-Plank predictions of auxiliary current into single fluid MHD. Simulations at low Lundquist number (S ∼ 10{sup 4}) generally agree with the previous work; significantly more burdensome simulations at MST-like Lundquist number (S ∼ 3×10{sup 6}) show unexpected results. The effect on nonlinearly saturated current profile by a particular RF-driven external force decreases in magnitude and widens considerably as the Lundquist number increases toward experimental values. Simulations reproduce the periodic current profile relaxation events observed in experiment (sawteeth) in the absence of current profile control. Reduction of the tearing mode amplitudes is still observable; however, reduction is limited to periods between the large bursts of magnetic activity at each sawtooth. The sawtoothing pattern persists with up to 10 MW of externally applied RF power. Periods with prolonged low tearing amplitude are predicted with a combination of external current drive and a reduced toroidal loop voltage, consistent with previous conclusions. Finally, the resistivity profile is observed to have a strong effect on the optimal externally driven current profile for mode stabilization.
On the propagation of blobs in the magnetotail: MHD simulations
NASA Astrophysics Data System (ADS)
Birn, J.; Nakamura, R.; Hesse, M.
2013-09-01
Using three-dimensional magnetohydrodynamic (MHD) simulations of the magnetotail, we investigate the fate of entropy-enhanced localized magnetic flux tubes ("blobs"). Such flux tubes may be the result of a slippage process that also generates entropy-depleted flux tubes ("bubbles") or of a rapid localized energy increase, for instance, from wave absorption. We confirm the expectation that the entropy enhancement leads to a tailward motion and that the speed and distance traveled into the tail increase with the entropy enhancement, even though the blobs tend to break up into pieces. The vorticity on the outside of the blobs twists the magnetic field and generates field-aligned currents predominantly of region-2 sense (earthward on the dusk side and tailward on the dawn side), which might provide a possibility for remote identification from the ground. The breakup, however, leads to more turbulent flow patterns, associated with opposite vorticity and the generation of region-1 sense field-aligned currents of lower intensity but approximately equal integrated magnitude.
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.
A Real-time 3D Visualization of Global MHD Simulation for Space Weather Forecasting
NASA Astrophysics Data System (ADS)
Murata, K.; Matsuoka, D.; Kubo, T.; Shimazu, H.; Tanaka, T.; Fujita, S.; Watari, S.; Miyachi, H.; Yamamoto, K.; Kimura, E.; Ishikura, S.
2006-12-01
Recently, many satellites for communication networks and scientific observation are launched in the vicinity of the Earth (geo-space). The electromagnetic (EM) environments around the spacecraft are always influenced by the solar wind blowing from the Sun and induced electromagnetic fields. They occasionally cause various troubles or damages, such as electrification and interference, to the spacecraft. It is important to forecast the geo-space EM environment as well as the ground weather forecasting. Owing to the recent remarkable progresses of super-computer technologies, numerical simulations have become powerful research methods in the solar-terrestrial physics. For the necessity of space weather forecasting, NICT (National Institute of Information and Communications Technology) has developed a real-time global MHD simulation system of solar wind-magnetosphere-ionosphere couplings, which has been performed on a super-computer SX-6. The real-time solar wind parameters from the ACE spacecraft at every one minute are adopted as boundary conditions for the simulation. Simulation results (2-D plots) are updated every 1 minute on a NICT website. However, 3D visualization of simulation results is indispensable to forecast space weather more accurately. In the present study, we develop a real-time 3D webcite for the global MHD simulations. The 3-D visualization results of simulation results are updated every 20 minutes in the following three formats: (1)Streamlines of magnetic field lines, (2)Isosurface of temperature in the magnetosphere and (3)Isoline of conductivity and orthogonal plane of potential in the ionosphere. For the present study, we developed a 3-D viewer application working on Internet Explorer browser (ActiveX) is implemented, which was developed on the AVS/Express. Numerical data are saved in the HDF5 format data files every 1 minute. Users can easily search, retrieve and plot past simulation results (3D visualization data and numerical data) by using
NASA Astrophysics Data System (ADS)
Farr, Nathan; Baker, Daniel N.; Wiltberger, Michael
We present the results of a global magnetohydrodynamic (MHD) simulation of a pair of substorms on August 11, 2002. Comparisons of data with simulation results reveal an agreement regarding the timing and sequence of events in the magnetosphere. We then present the results in the simulation of a plasmoid flux rope formed during the second substorm. Unlike standard 2-D depictions of reconnection and plasmoid release during a substorm, the simulation shows a highly complex structure that has considerable winding of both closed and open field lines. Additionally, the flux rope does not move tailward uniformly, but rather has a asymmetric motion in which the dawn flank portion moves tailward prior to the dusk end of the flux rope, resulting in a a skewed flux rope that runs almost downtail instead of crosstail. The flux rope structure begins on closed field lines, then progresses through open field lines and finally onto IMF field lines. Using the global simulation we can connect the large scale magnetotail reconfiguration with features on the auroral boundary of the model. The results shown here are in agreement with existing MHD tail simulations as well as flux rope observations from satellite data.
COYOTE: A computer program for 2-D reactive flow simulations
Cloutman, L.D.
1990-04-01
We describe the numerical algorithm used in the COYOTE two- dimensional, transient, Eulerian hydrodynamics program for reactive flows. The program has a variety of options that provide capabilities for a wide range of applications, and it is designed to be robust and relatively easy to use while maintaining adequate accuracy and efficiency to solve realistic problems. It is based on the ICE method, and it includes a general species and chemical reaction network for simulating reactive flows. It also includes swirl, turbulence transport models, and a nonuniform mesh capability. We describe several applications of the program. 33 refs., 4 figs.
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.
NASA Astrophysics Data System (ADS)
Olmi, B.; Del Zanna, L.; Amato, E.; Bucciantini, N.; Bandiera, R.
2015-09-01
Pulsar wind nebulae are among the most powerful particle accelerators in the Galaxy with acceleration efficiencies that reach up to 30% and maximum particle energies in the PeV range. In recent years relativistic axisymmetric MHD models have proven to be excellent tools for describing the physics of such objects, and particularly successful at explaining their high energy morphology, down to very fine details. Nevertheless, some important aspects of the physics of PWNe are still obscure: the mechanism(s) responsible for the acceleration of particles of all energies is (are) still unclear, and the origin of the lowest energy (radio emitting) particles is most mysterious. The correct interpretation of the origin of radio emitting particles is of fundamental importance, as this holds information about the amount of pair production in the pulsar magnetosphere, and hence on the role of pulsars as antimatter factories. On the other hand, the long lifetimes of these particles against synchrotron losses, allows them to travel far from their injection location, making their acceleration site difficult to constrain. As far as the highest energy (X and gamma-ray emitting) particles are concerned, their acceleration is commonly believed to occur at the pulsar wind termination shock. But since the upstream flow is thought to have non-uniform properties along the shock surface, important constraints on the acceleration mechanism(s) could come from exact knowledge of the location and flow properties where particles are being accelerated. We investigate in detail both topics by means of 2D numerical MHD simulations. Different assumptions on the origin of radio particles and more generally on the injection sites of all particles are considered, and the corresponding emission properties are computed. We discuss the physical constraints that can be inferred from comparison of the synthetic emission properties against multiwavelength observations of the PWN class prototype, the Crab
An MHD simulation of plasmoid instability in the dayside ionosphere of an unmagnetized planet
NASA Astrophysics Data System (ADS)
Hitoshi, S.; Terada, N.; Kasaba, Y.
2015-12-01
A numerical simulation of magnetic reconnection in the dayside ionosphere of an unmagnetized planet and a comparison of the size distribution of flux ropes obtained from simulation with that from observation will be reported. Flux ropes have been frequently observed in the dayside ionospheres of Venus and Mars[Russell and Elphic, 1979; Cloutier et al.,1999] and their radius has been found to be between 6 to 12 km near the subsolar location of Venus[Russell et al., 1990]. Dreher et al. [1995] suggested using an MHD simulation that reconnection caused by an IMF rotation can generate flux ropes at the Venus ionopause. However, Dreher et al. [1995] examined only the linear stage of reconnection, so the nonlinear stage that takes into consideration the vertical convection of the reconnection site along the intrinsic convection in the Venus ionosphere has yet to be investigated. In this study, using a 2-D multi-species MHD simulation, the spatiotemporal evolution of reconnection in the ionosphere of Venus is examined. The size distribution of flux ropes is also examined and the validity of the generation mechanism of flux ropes is discussed by comparing the rope size distribution with the observed one. In the ionosphere of Venus, our simulation result shows that plasmoid instability [Loureiro et al., 2007] occurs in a Sweet-Parker (SP) current sheet above the altitude where Lundquist number exceeds 106, and consequently many plasmoids are generated. In the nonlinear stage, secondary reconnections occur in the current sheets, which exist between adjacent pairs of plasmoids, and thus smaller flux ropes are created. It has been suggested that the smaller side of the size distribution increases as a result of hierarchical reconnections in the SP current sheet [Shibata et al., 2001]. The observational size distribution [Vignes et al., 2003] shows that the population of small flux ropes is larger than that of large ones and it is consistent with the simulation result. Through
Relativistic MHD Simulations of Poynting Flux-driven Jets
NASA Astrophysics Data System (ADS)
Guan, Xiaoyue; Li, Hui; Li, Shengtai
2014-01-01
Relativistic, magnetized jets are observed to propagate to very large distances in many active galactic nuclei (AGNs). We use three-dimensional relativistic MHD simulations to study the propagation of Poynting flux-driven jets in AGNs. These jets are already assumed to be being launched from the vicinity (~103 gravitational radii) of supermassive black holes. Jet injections are characterized by a model described in Li et al., and we follow the propagation of these jets to ~parsec scales. We find that these current-carrying jets are always collimated and mildly relativistic. When α, the ratio of toroidal-to-poloidal magnetic flux injection, is large the jet is subject to nonaxisymmetric current-driven instabilities (CDI) which lead to substantial dissipation and reduced jet speed. However, even with the presence of instabilities, the jet is not disrupted and will continue to propagate to large distances. We suggest that the relatively weak impact by the instability is due to the nature of the instability being convective and the fact that the jet magnetic fields are rapidly evolving on Alfvénic time scales. We present the detailed jet properties and show that far from the jet launching region, a substantial amount of magnetic energy has been transformed into kinetic energy and thermal energy, producing a jet magnetization number σ < 1. In addition, we have also studied the effects of a gas pressure supported "disk" surrounding the injection region, and qualitatively similar global jet behaviors were observed. We stress that jet collimation, CDIs, and the subsequent energy transitions are intrinsic features of current-carrying jets.
Simulation of subgrid orographic precipitation with an embedded 2-D cloud-resolving model
NASA Astrophysics Data System (ADS)
Jung, Joon-Hee; Arakawa, Akio
2016-03-01
By explicitly resolving cloud-scale processes with embedded two-dimensional (2-D) cloud-resolving models (CRMs), superparameterized global atmospheric models have successfully simulated various atmospheric events over a wide range of time scales. Up to now, however, such models have not included the effects of topography on the CRM grid scale. We have used both 3-D and 2-D CRMs to simulate the effects of topography with prescribed "large-scale" winds. The 3-D CRM is used as a benchmark. The results show that the mean precipitation can be simulated reasonably well by using a 2-D representation of topography as long as the statistics of the topography such as the mean and standard deviation are closely represented. It is also shown that the use of a set of two perpendicular 2-D grids can significantly reduce the error due to a 2-D representation of topography.
Attempts to Simulate Anisotropies of Solar Wind Fluctuations Using MHD with a Turning Magnetic Field
NASA Technical Reports Server (NTRS)
Ghosh, Sanjoy; Roberts, D. Aaron
2010-01-01
We examine a "two-component" model of the solar wind to see if any of the observed anisotropies of the fields can be explained in light of the need for various quantities, such as the magnetic minimum variance direction, to turn along with the Parker spiral. Previous results used a 3-D MHD spectral code to show that neither Q2D nor slab-wave components will turn their wave vectors in a turning Parker-like field, and that nonlinear interactions between the components are required to reproduce observations. In these new simulations we use higher resolution in both decaying and driven cases, and with and without a turning background field, to see what, if any, conditions lead to variance anisotropies similar to observations. We focus especially on the middle spectral range, and not the energy-containing scales, of the simulation for comparison with the solar wind. Preliminary results have shown that it is very difficult to produce the required variances with a turbulent cascade.
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
Comparison of Iridium Determined Field-Aligned Current Patterns with MHD Simulations
NASA Astrophysics Data System (ADS)
Korth, H.; Anderson, B. J.; Goodrich, C. C.; Waters, C. L.; Merkine, V. G.
2002-05-01
The engineering magnetometers aboard the 70+ Iridium satellites arranged in six equally spaced polar orbital planes provide a unique database for determination of global field-aligned currents [Waters et al., 2001]. In this study we compare these field-aligned currents with MHD simulation results to quantitatively evaluate the MHD results in a global way. We report analysis for three events of steady interplanetary magnetic field (IMF) orientation, stable to within 25o of the average direction. The start times of these intervals are August~11, 1999 (22:36), November~23, 1999 (07:15), and August~10, 2000 (22:11), and the events extend between eight and ten hours in duration. The IMF clock angles for the events are -124o, 125o, and 160o, respectively, and the IMF cone angles for all three intervals are within 25o of 90o. The solar wind flow speeds for the events averages 430, 453, and 386~km/s, and the mean solar wind densities are 3.7, 3.6, and 12.0 {cm}-3, respectively. The field aligned current densities in the MHD simulations are evaluated at the inner simulation boundary (2~Re) and mapped on dipole field lines to ionospheric altitudes. Preliminary results show a reasonably good agreement in the morphology of the Region-1 currents, although the field-aligned currents of the MHD simulations are displaced somewhat poleward with respect to the Iridium patterns. DMSP particle source identifications are used to compare source regions of Region 1 in the observed FAC maps with those in the MHD simulations. The Region-2 currents show expectedly larger differences since ring current drift physics necessary to drive these currents in the magnetosphere is not implemented in the MHD evaluations. The ratio between Region~1 and Region~2 is used to measure the relative deficit of Region-2 currents in the MHD simulation results.
NASA Astrophysics Data System (ADS)
Toth, G.; Daldorff, L. K. S.; Jia, X.; Gombosi, T. I.; Lapenta, G.
2014-12-01
We have recently developed a new modeling capability to embed theimplicit Particle-in-Cell (PIC) model iPIC3D into the BATS-R-USmagnetohydrodynamic model. The PIC domain can cover the regions wherekinetic effects are most important, such as reconnection sites. TheBATS-R-US code, on the other hand, can efficiently handle the rest ofthe computational domain where the MHD or Hall MHD description issufficient. As one of the very first applications of the MHD-EPICalgorithm (Daldorff et al. 2014, JCP, 268, 236) we simulate theinteraction between Jupiter's magnetospheric plasma with Ganymede'smagnetosphere, where the separation of kinetic and global scalesappears less severe than for the Earth's magnetosphere. Because theexternal Jovian magnetic field remains in an anti-parallel orientationwith respect to Ganymede's intrinsic magnetic field, magneticreconnection is believed to be the major process that couples the twomagnetospheres. As the PIC model is able to describe self-consistentlythe electron behavior, our coupled MHD-EPIC model is well suited forinvestigating the nature of magnetic reconnection in thisreconnection-driven mini-magnetosphere. We will compare the MHD-EPICsimulations with pure Hall MHD simulations and compare both modelresults with Galileo plasma and magnetic field measurements to assess therelative importance of ion and electron kinetics in controlling theconfiguration and dynamics of Ganymede's magnetosphere.
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.
2-D simulation of a waveguide free electron laser having a helical undulator
Kim, S.K.; Lee, B.C.; Jeong, Y.U.
1995-12-31
We have developed a 2-D simulation code for the calculation of output power from an FEL oscillator having a helical undulator and a cylindrical waveguide. In the simulation, the current and the energy of the electron beam is 2 A and 400 keV, respectively. The parameters of the permanent-magnet helical undulator are : period = 32 mm, number of periods = 20, magnetic field = 1.3 kG. The gain per pass is 10 and the output power is calculated to be higher than 10 kW The results of the 2-D simulation are compared with those of 1-D simulation.
Daily Coronal MHD Simulation Using HMI Near-Real-Time Magnetograms
NASA Astrophysics Data System (ADS)
Hayashi, Keiji; HMI Team
2012-05-01
SDO/HMI is making full-disk line-of-sight magnetogram measurements with a cadence of 45 seconds. The HMI analysis pipeline regularly generates two types of synoptic map of the solar surface magnetic field. Definitive calibrated data maps are created every Carrington Rotation, about every 27 days and a preliminary synoptic map is updated on a near-real-time basis. As an application of the near-real-time data, we have been running a daily MHD simulation of the global solar corona using the photospheric map as the boundary condition ( http://hmi.stanford.edu/MHD ). The daily MHD model assumes a polytropic gas with the specific heat ratio of 1.05, and the simulation is conducted in a 4-pi spherical grid system with latitudinal and longitudinal grid sizes of pi/64. The output available at hmi.stanford.edu/MHD includes the three-dimensional volume data, the shape of the open-field regions corresponding to the coronal holes, and the LoS-integration of the coronal density mimicking coronagraph observations. For validation, we compare the results of the low-resolution daily MHD simulation and the high-resolution PFSS calculation with SDO/AIA and SOHO/C2 and C3 image data. In the future the simulation region will be extended to 1 AU, and models of coronal heating and acceleration will be applied to allow a timely prediction of solar wind at the Earth for space weather purposes.
Accuracy of MHD simulations: Effects of simulation initialization in GUMICS-4
NASA Astrophysics Data System (ADS)
Lakka, Antti; Pulkkinen, Tuija; Dimmock, Andrew; Osmane, Adnane; Palmroth, Minna; Honkonen, Ilja
2016-04-01
We conducted a study aimed at revealing how different global magnetohydrodynamic (MHD) simulation initialization methods affect the dynamics in different parts of the Earth's magnetosphere-ionosphere system. While such magnetosphere-ionosphere coupling codes have been used for more than two decades, their testing still requires significant work to identify the optimal numerical representation of the physical processes. We used the Grand Unified Magnetosphere-Ionosphere Coupling Simulation (GUMICS-4), the only European global MHD simulation being developed by the Finnish Meteorological Institute. GUMICS-4 was put to a test that included two stages: 1) a 10 day Omni data interval was simulated and the results were validated by comparing both the bow shock and the magnetopause spatial positions predicted by the simulation to actual measurements and 2) the validated 10 day simulation run was used as a reference in a comparison of five 3 + 12 hour (3 hour synthetic initialisation + 12 hour actual simulation) simulation runs. The 12 hour input was not only identical in each simulation case but it also represented a subset of the 10 day input thus enabling quantifying the effects of different synthetic initialisations on the magnetosphere-ionosphere system. The used synthetic initialisation data sets were created using stepwise, linear and sinusoidal functions. Switching the used input from the synthetic to real Omni data was immediate. The results show that the magnetosphere forms in each case within an hour after the switch to real data. However, local dissimilarities are found in the magnetospheric dynamics after formation depending on the used initialisation method. This is evident especially in the inner parts of the lobe.
Multidimensional MHD Simulations Of DSA Using AstroBEAR
NASA Astrophysics Data System (ADS)
Edmon, Paul; Jones, T.; Mitran, S.; Cunningham, A.; Frank, A.
2009-05-01
We present a modification to the AstroBEAR (Astronomical Boundary Embedded Adaptive Refinement) MHD code (Cunningham et. al. 2007) that allows it to treat time dependent Diffusive Shock Acceleration (DSA) of cosmic rays in multiple dimensions including dynamical feedback from the cosmic rays. Utilizing the power of Adaptive Mesh Refinement (AMR) in tandem with efficient methods for cosmic ray diffusion and advection, this allows us for the first time to explore the evolution of modified MHD shocks in more than one spatial dimension. Among the early applications of the code will be investigations of colliding and clumpy stellar winds, type II supernova remnants and cosmic ray driven instabilities. This work is supported at the University of Minnesota by NSF, NASA and the Minnesota Supercomputing Institute.
MHD simulations for investigating interaction processes between a CME and ambient solar wind
NASA Astrophysics Data System (ADS)
An, Junmo; Magara, Tetsuya
2016-05-01
The interaction between coronal mass ejections (CMEs) and ambient solar winds is one of the important issues of space weather because it affects the trajectory of a flying CME, which determines whether the CME hits the Earth and produces geomagnetic disturbances or not. In this study, two-step 3D magnetohydrodynamic (MHD) simulations including a spheromak-type CME and an ambient solar wind are performed to investigate their interaction processes such as deflection and rotation of a CME. We perform the 1st-step MHD simulation using averaged surface magnetic field data to construct a steady state with an ambient solar wind. A spheromak-type CME is then injected through the solar surface, and subsequent evolution is reproduced by performing the 2nd-step MHD simulation. We discuss key parameters that characterize interaction processes between a CME and ambient solar wind.
Forced Reconnection in the Near Magnetotail: Onset and Energy Conversion in PIC and MHD Simulations
NASA Technical Reports Server (NTRS)
Birn, J.; Hesse, Michael
2014-01-01
Using two-dimensional particle-in-cell (PIC) together with magnetohydrodynamic (MHD) Q1 simulations of magnetotail dynamics, we investigate the evolution toward onset of reconnection and the subsequent energy transfer and conversion. In either case, reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. The boundary deformation leads to the formation of thin embedded current sheets, which are bifurcated in the near tail, converging to a single sheet farther out in the MHD simulations. The thin current sheets in the PIC simulation are carried by electrons and are associated with a strong perpendicular electrostatic field, which may provide a connection to parallel potentials and auroral arcs and an ionospheric signal even prior to the onset of reconnection. The PIC simulation very well satisfies integral entropy conservation (intrinsic to ideal MHD) during this phase, supporting ideal ballooning stability. Eventually, the current intensification leads to the onset of reconnection, the formation and ejection of a plasmoid, and a collapse of the inner tail. The earthward flow shows the characteristics of a dipolarization front: enhancement of Bz, associated with a thin vertical electron current sheet in the PIC simulation. Both MHD and PIC simulations show a dominance of energy conversion from incoming Poynting flux to outgoing enthalpy flux, resulting in heating of the inner tail. Localized Joule dissipation plays only a minor role.
Phase transition-like behavior of magnetospheric substorms: Global MHD simulation results
NASA Astrophysics Data System (ADS)
Shao, X.; Sitnov, M. I.; Sharma, S. A.; Papadopoulos, K.; Goodrich, C. C.; Guzdar, P. N.; Milikh, G. M.; Wiltberger, M. J.; Lyon, J. G.
2003-01-01
Using nonlinear dynamical techniques, we statistically investigate whether the simulated substorms from global magnetohydrodynamic (MHD) models have a combination of global and multiscale features, revealed in substorm dynamics by [2000] and featured the phase transition-like behavior. We simulate seven intervals of total duration of 280 hours from the data set used in the above works [, 1985]. We analyze the input-output (vBs-pseudo AL index) system obtained from the global MHD model and compare the results to those inferred from the original set (vBs-observed AL index). The analysis of the coupled vBs-pseudo AL index system shows the first-order phase transition map, which is consistent with the map obtained for the vBs-observed AL index system. Although the comparison between observations and global MHD simulations for individual events may vary, the overall global transition pattern during the substorm cycle revealed by singular spectrum analysis (SSA) is statistically consistent between simulations and observations. The coupled vBs-pseudo AL index system also shows multiscale behavior (scale-invariant power law dependence) in SSA power spectrum. Besides, we find the critical exponent of the nonequilibrium transitions in the magnetosphere, which reflects the multiscale aspect of the substorm activity, different from power law frequency of autonomous systems. The exponent relates input and output parameters of the magnetosphere. We also discuss the limitations of the global MHD model in reproducing the multiscale behavior when compared to the real system.
The flare position obtained from MHD simulation and comparison with X-ray observations
NASA Astrophysics Data System (ADS)
Podgorny, Alexander; Podgorny, Igor
It was for the first time shown that the position of the current sheet, obtained by numerical MHD simulation, coincides with the position of the thermal X-ray source. In our 3D MHD simulation we do not use any hypotheses about the flare mechanism. Several mechanisms of solar flare production are considered by different authors. Usually the initial conditions at numerical simulation are artificially set such a way that it is required for development of the proposed mechanism. In this approach, the unstable configuration of the magnetic field is set as the initial conditions, and the possibility of forming such an unstable system at the real evolution of the active region before the flare is not considered. Here the flare mechanism is obtained from the numerical MHD simulations in which all the conditions are taken from observations in the active region. It is shown that flare energy accumulation occurs in the current sheet magnetic field created by disturbances focusing in the vicinity of an X-type singular line. According to the developed solar flare electrodynamical model the thermal X-ray emission source appears in a current sheet, where plasma is heated due to magnetic field dissipation. Using 3D MHD numerical simulation the position of source of thermal X-ray emission are found for the flare occurred May 27, 2003 at 02:53. To find positions of sources of thermal X-ray radiation in the corona from MHD simulation results the graphical system is developed. The comparison with RHESSI X-ray observations show the coincidence of current sheet and observed the thermal X-ray emission source.
Numerical Simulation of MHD Effect in Liquid Metal Blankets with Flow Channel Insert
NASA Astrophysics Data System (ADS)
Mao, J.; Pan, H. C.
2011-09-01
The magnetohydrodynamic effect in liquid metal blankets with flow channel insert and pressure equalization slot for fusion liquid metal blanket is studied by numerical simulation based on two dimensional fully developed flow model. The code is verified by comparing analytical solution and numerical solution of Hunt Case II. The velocity field and MHD pressure drop varying with electric conductivity of the FCI is analyzed. The result shows that the average velocity in central area of the cross section decreases with the increase of the electric conductivity of FCI. While the average velocity in gap zone is reverse. Comparing with MHD duct flow without FCI, MHD pressure drop is reduced significantly when the FCI material is electrically insulating.
The simulation of 3D microcalcification clusters in 2D digital mammography and breast tomosynthesis
Shaheen, Eman; Van Ongeval, Chantal; Zanca, Federica; Cockmartin, Lesley; Marshall, Nicholas; Jacobs, Jurgen; Young, Kenneth C.; Dance, David R.; Bosmans, Hilde
2011-12-15
Purpose: This work proposes a new method of building 3D models of microcalcification clusters and describes the validation of their realistic appearance when simulated into 2D digital mammograms and into breast tomosynthesis images. Methods: A micro-CT unit was used to scan 23 breast biopsy specimens of microcalcification clusters with malignant and benign characteristics and their 3D reconstructed datasets were segmented to obtain 3D models of microcalcification clusters. These models were then adjusted for the x-ray spectrum used and for the system resolution and simulated into 2D projection images to obtain mammograms after image processing and into tomographic sequences of projection images, which were then reconstructed to form 3D tomosynthesis datasets. Six radiologists were asked to distinguish between 40 real and 40 simulated clusters of microcalcifications in two separate studies on 2D mammography and tomosynthesis datasets. Receiver operating characteristic (ROC) analysis was used to test the ability of each observer to distinguish between simulated and real microcalcification clusters. The kappa statistic was applied to assess how often the individual simulated and real microcalcification clusters had received similar scores (''agreement'') on their realistic appearance in both modalities. This analysis was performed for all readers and for the real and the simulated group of microcalcification clusters separately. ''Poor'' agreement would reflect radiologists' confusion between simulated and real clusters, i.e., lesions not systematically evaluated in both modalities as either simulated or real, and would therefore be interpreted as a success of the present models. Results: The area under the ROC curve, averaged over the observers, was 0.55 (95% confidence interval [0.44, 0.66]) for the 2D study, and 0.46 (95% confidence interval [0.29, 0.64]) for the tomosynthesis study, indicating no statistically significant difference between real and simulated
Phase Transition-like Behavior of Magnetospheric Substorms: Global MHD Simulation Results
NASA Astrophysics Data System (ADS)
Shao, X.; Sitnov, M.; Sharma, A. S.; Papadopoulos, K.; Guzdar, P. N.; Goodrich, C. C.; Milikh, G. M.; Wiltberger, M. J.; Lyon, J. G.
2001-12-01
Because of their relevance to massive global energy loading and unloading, lots of observations and studies have been made for magnetic substorm events. Using nonlinear dynamical techniques, we investigate whether the simulated substorms from global MHD models have the non-equilibrium phase transition-like features revealed by \\markcite{Sitnov et al. [2000]}. We simulated 6 intervals of total duration of 240 hours from the same data set used in Sitnov et al. [2000]. We analyzed the input-output (vBs--pseudo-AL index) system obtained from the global MHD model and compared the results to those in \\markcite{Sitnov et al. [2000, 2001]}. The analysis of the coupled vBs--pseudo-AL index system shows the first-order phase transition map, which is consistent with the map obtained for the vBs--observed-AL index system from Sitnov et al. [2000]. The explanation lies in the cusp catastrophe model proposed by Lewis [1991]. Although, the comparison between observation and individual global MHD simulations may vary, the overall global transition pattern during the substorm cycle revealed by Singular Spectrum Analysis (SSA) is consistent between simulations and observations. This is an important validation of the global MHD simulations of the magnetosphere. The coupled vBs--pseudo-AL index system shows multi-scale behavior (scale-invarianet power-law dependence) in singular power spectrum. We found critical exponents of the non-equilibrium transitions in the magnetosphere, which reflect the multi-scale aspect of the substorm activity, different from power-law frequency of autonomous systems. The exponents relate input and output parameters of the magnetosphere and distinguish the second order phase transition model from the self-organized criticality model. We also discuss the limitations of the global MHD model in reproducing the multi-scale behavior when compared to the real system.
Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations
NASA Technical Reports Server (NTRS)
Frank, Adam; Jones, T. W.; Ryu, Dongsu
1991-01-01
We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.
3D multiple-point statistics simulation using 2D training images
NASA Astrophysics Data System (ADS)
Comunian, A.; Renard, P.; Straubhaar, J.
2012-03-01
One of the main issues in the application of multiple-point statistics (MPS) to the simulation of three-dimensional (3D) blocks is the lack of a suitable 3D training image. In this work, we compare three methods of overcoming this issue using information coming from bidimensional (2D) training images. One approach is based on the aggregation of probabilities. The other approaches are novel. One relies on merging the lists obtained using the impala algorithm from diverse 2D training images, creating a list of compatible data events that is then used for the MPS simulation. The other (s2Dcd) is based on sequential simulations of 2D slices constrained by the conditioning data computed at the previous simulation steps. These three methods are tested on the reproduction of two 3D images that are used as references, and on a real case study where two training images of sedimentary structures are considered. The tests show that it is possible to obtain 3D MPS simulations with at least two 2D training images. The simulations obtained, in particular those obtained with the s2Dcd method, are close to the references, according to a number of comparison criteria. The CPU time required to simulate with the method s2Dcd is from two to four orders of magnitude smaller than the one required by a MPS simulation performed using a 3D training image, while the results obtained are comparable. This computational efficiency and the possibility of using MPS for 3D simulation without the need for a 3D training image facilitates the inclusion of MPS in Monte Carlo, uncertainty evaluation, and stochastic inverse problems frameworks.
Numerical simulation of ( T 2, T 1) 2D NMR and fluid responses
NASA Astrophysics Data System (ADS)
Tan, Mao-Jin; Zou, You-Long; Zhang, Jin-Yan; Zhao, Xin
2012-12-01
One-dimensional nuclear magnetic resonance (1D NMR) logging technology is limited for fluid typing, while two-dimensional nuclear magnetic resonance (2D NMR) logging can provide more parameters including longitudinal relaxation time ( T 1) and transverse relaxation time ( T 2) relative to fluid types in porous media. Based on the 2D NMR relaxation mechanism in a gradient magnetic field, echo train simulation and 2D NMR inversion are discussed in detail. For 2D NMR inversion, a hybrid inversion method is proposed based on the damping least squares method (LSQR) and an improved truncated singular value decomposition (TSVD) algorithm. A series of spin echoes are first simulated with multiple waiting times ( T W s) in a gradient magnetic field for given fluid models and these synthesized echo trains are inverted by the hybrid method. The inversion results are consistent with given models. Moreover, the numerical simulation of various fluid models such as the gas-water, light oil-water, and vicious oil-water models were carried out with different echo spacings ( T E s) and T W s by this hybrid method. Finally, the influences of different signal-to-noise ratios (SNRs) on inversion results in various fluid models are studied. The numerical simulations show that the hybrid method and optimized observation parameters are applicable to fluid typing of gas-water and oil-water models.
Using Two-Ribbon Flare Observations and MHD Simulations to Constrain Flare Properties
NASA Astrophysics Data System (ADS)
Kazachenko, Maria D.; Lynch, Benjamin J.; Welsch, Brian
2016-05-01
Flare ribbons are emission structures that are frequently observed during flares in transition-region and chromospheric radiation. These typically straddle a polarity inversion line (PIL) of the radial magnetic field at the photosphere, and move apart as the flare progresses. The ribbon flux - the amount of unsigned photospheric magnetic flux swept out by flare ribbons - is thought to be related to the amount coronal magnetic reconnection, and hence provides a key diagnostic tool for understanding the physical processes at work in flares and CMEs. Previous measurements of the magnetic flux swept out by flare ribbons required time-consuming co-alignment between magnetograph and intensity data from different instruments, explaining why those studies only analyzed, at most, a few events. The launch of the Helioseismic and Magnetic Imager (HMI) and the Atmospheric Imaging Assembly (AIA), both aboard the Solar Dynamics Observatory (SDO), presented a rare opportunity to compile a much larger sample of flare-ribbon events than could readily be assembled before. We created a dataset of 363 events of both flare ribbon positions and fluxes, as a function of time, for all C9.-class and greater flares within 45 degrees of disk center observed by SDO from June 2010 till April 2015. For this purpose, we used vector magnetograms (2D magnetic field maps) from HMI and UV images from AIA. A critical problem with using unprocessed AIA data is the existence of spurious intensities in AIA data associated with strong flare emission, most notably "blooming" (spurious smearing of saturated signal into neighboring pixels, often in streaks). To overcome this difficulty, we have developed an algorithmic procedure that effectively excludes artifacts like blooming. We present our database and compare statistical properties of flare ribbons, e.g. evolutions of ribbon reconnection fluxes, reconnection flux rates and vertical currents with the properties from MHD simulations.
Comparison of empirical magnetic field models and global MHD simulations: The near-tail currents
NASA Technical Reports Server (NTRS)
Pulkkinen, T. I.; Baker, D. N.; Walker, R. J.; Raeder, J.; Ashour-Abdalla, M.
1995-01-01
The tail currents predicted by empirical magnetic field models and global MHD simulations are compared. It is shown that the near-Earth currents obtained from the MHD simulations are much weaker than the currents predicted by the Tsyganenko models, primarily because the ring current is not properly represented in the simulations. On the other hand, in the mid-tail and distant tail the lobe field strength predicted by the simulations is comparable to what is observed at about 50 R(sub E) distance, significantly larger than the very low lobe field values predicted by the Tsyganenko models at that distance. Ways to improve these complementary approaches to model the actual magnetospheric configuration are discussed.
Comparison of empirical magnetic field models and global MHD simulations: The near-tail currents
Pulkkinen, T.I.; Baker, D.N.; Walker, R.J.
1995-03-15
The tail currents predicted by empirical magnetic field models and global MHD simulations are compared. It is shown that the near-Earth currents obtained from the MHD simulations are much weaker than the currents predicted by the Tsyganenko models, primarily because the ring current is not properly represented in the simulations. On the other hand, in the mid-tail and distant tail the lobe field strength predicted by the simulations is comparable to what is observed at about 50R{sub E} distance, significantly larger than the very low lobe field values predicted by the Tsyganenko models at that distance. Ways to improve these complementary approaches to model the actual magnetospheric configuration are discussed. 11 refs., 3 figs.
NASA Astrophysics Data System (ADS)
Huang, Z.; Jia, X.; Rubin, M.; Fougere, N.; Gombosi, T. I.; Tenishev, V.; Combi, M. R.; Bieler, A. M.; Toth, G.; Hansen, K. C.; Shou, Y.
2014-12-01
We study the plasma environment of the comet Churyumov-Gerasimenko, which is the target of the Rosetta mission, by performing large scale numerical simulations. Our model is based on BATS-R-US within the Space Weather Modeling Framework that solves the governing multifluid MHD equations, which describe the behavior of the cometary heavy ions, the solar wind protons, and electrons. The model includes various mass loading processes, including ionization, charge exchange, dissociative ion-electron recombination, as well as collisional interactions between different fluids. The neutral background used in our MHD simulations is provided by a kinetic Direct Simulation Monte Carlo (DSMC) model. We will simulate how the cometary plasma environment changes at different heliocentric distances.
2D and 3D Mass Transfer Simulations in β Lyrae System
NASA Astrophysics Data System (ADS)
Nazarenko, V. V.; Glazunova, L. V.; Karetnikov, V. G.
2001-12-01
2D and 3D mass transfer simulations of the mass transfer in β Lyrae binary system. We have received that from a point L3 40 per cent of mass transfer from L1-point is lost.The structure of a gas envelope, around system is calculated.3-D mass transfer simulations has shown presence the spiral shock in the disk around primary star's and a jet-like structures (a mass flow in vertical direction) over a stream.
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.
The simulation of 3D mass models in 2D digital mammography and breast tomosynthesis
Shaheen, Eman De Keyzer, Frederik; Bosmans, Hilde; Ongeval, Chantal Van; Dance, David R.; Young, Kenneth C.
2014-08-15
Purpose: This work proposes a new method of building 3D breast mass models with different morphological shapes and describes the validation of the realism of their appearance after simulation into 2D digital mammograms and breast tomosynthesis images. Methods: Twenty-five contrast enhanced MRI breast lesions were collected and each mass was manually segmented in the three orthogonal views: sagittal, coronal, and transversal. The segmented models were combined, resampled to have isotropic voxel sizes, triangularly meshed, and scaled to different sizes. These masses were referred to as nonspiculated masses and were then used as nuclei onto which spicules were grown with an iterative branching algorithm forming a total of 30 spiculated masses. These 55 mass models were projected into 2D projection images to obtain mammograms after image processing and into tomographic sequences of projection images, which were then reconstructed to form 3D tomosynthesis datasets. The realism of the appearance of these mass models was assessed by five radiologists via receiver operating characteristic (ROC) analysis when compared to 54 real masses. All lesions were also given a breast imaging reporting and data system (BIRADS) score. The data sets of 2D mammography and tomosynthesis were read separately. The Kendall's coefficient of concordance was used for the interrater observer agreement assessment for the BIRADS scores per modality. Further paired analysis, using the Wilcoxon signed rank test, of the BIRADS assessment between 2D and tomosynthesis was separately performed for the real masses and for the simulated masses. Results: The area under the ROC curves, averaged over all observers, was 0.54 (95% confidence interval [0.50, 0.66]) for the 2D study, and 0.67 (95% confidence interval [0.55, 0.79]) for the tomosynthesis study. According to the BIRADS scores, the nonspiculated and the spiculated masses varied in their degrees of malignancy from normal (BIRADS 1) to highly
2D-simulation of wet steam flow in a steam turbine with spontaneous condensation
NASA Astrophysics Data System (ADS)
Sun, Lan-Xin; Zheng, Qun; Liu, Shun-Long
2007-06-01
Removal of condensates from wet steam flow in the last stages of steam turbines significantly promotes stage efficiency and prevents erosion of rotors. In this paper, homogeneous spontaneous condensation in transonic steam flow in the 2-D rotor-tip section of a stage turbine is investigated. Calculated results agree with experimental data reasonably well. On the basis of the above work, a 2-D numerical simulation of wet steam flow in adjacent root sections of a complex steam turbine stage was carried out. Computational results were analyzed and provide insights into effective removal of humidity.
Global Hall-MHD simulations of magnetorotational instability in a plasma Couette flow experiment
Ebrahimi, F.; Lefebvre, B.; Bhattacharjee, A.; Forest, C. B.
2011-06-15
Global MHD and Hall-MHD numerical simulations relevant to the Madison plasma Couette flow experiment (MPCX) have been performed using the extended MHD code NIMROD. The MPCX has been constructed to study the magnetorotational instability (MRI) in a plasma. The two-fluid Hall effect, which is relevant to some astrophysical situations such as protostellar disks, is also expected to be important in the MPCX. Here, we first derive the local Hall dispersion relation including viscosity, extending earlier work by Balbus and Terquem [Astrophys. J. 552, 235 (2001)]. The predictions of the local analysis are then compared with nonlocal calculations of linear stability of the MRI for a parameter range relevant to the MPCX. It is found that the MHD stability limit and mode structure are altered by the Hall term, and nonlocal analysis is necessary to obtain quantitatively reliable predictions for MPCX. Two-fluid physics also significantly changes the nonlinear evolution and saturation of the axisymmetric MRI. Both the Reynolds and Maxwell stresses contribute significantly to momentum transport. In the Hall regime, when the magnetic field is parallel to the rotation axis, the Maxwell stress is larger than the Reynolds stress (similar to the MHD regime). However, when the magnetic field is antiparallel to the rotation axis in the Hall regime, the Reynolds stress is much larger than the Maxwell stress. To further study the role of non-axisymmetric modes, we have also carried out fully nonlinear MHD computations. Non-axisymmetric modes play an increasingly important role as the magnetic Reynolds number increases and grow to large amplitudes in a saturated turbulent state.
2D-3D hybrid stabilized finite element method for tsunami runup simulations
NASA Astrophysics Data System (ADS)
Takase, S.; Moriguchi, S.; Terada, K.; Kato, J.; Kyoya, T.; Kashiyama, K.; Kotani, T.
2016-09-01
This paper presents a two-dimensional (2D)-three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier-Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.
2D-3D hybrid stabilized finite element method for tsunami runup simulations
NASA Astrophysics Data System (ADS)
Takase, S.; Moriguchi, S.; Terada, K.; Kato, J.; Kyoya, T.; Kashiyama, K.; Kotani, T.
2016-05-01
This paper presents a two-dimensional (2D)-three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier-Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.
Modeling anisotropic MHD turbulence in simulations of liquid metal flows
NASA Astrophysics Data System (ADS)
Widlund, O.
2001-06-01
The dynamical properties of the MHD turbulence model proposed by Widlund etal. are examined for the case of homogeneous decaying turbulence. The model is a Reynolds stress closure, extended with a transport equation for a dimensional anisotropy variable, α, which carries information about length scale anisotropy. The analysis suggests that the model term originally proposed for the nonlinear energy transfer in the α equation should be modified. A unique set of model coefficients could be determined, which makes the model consistent with theory and experiments for interaction parameters N ranging from zero to infinity. The model coincides with the standard K-eps model when there is no magnetic field. In the linear regime of large N, it produces the K˜ t^{-1/2} energy decay predicted by linear theory. When nonlinear effects are important, the model predicts K˜ t^{-1.7} and L_∥ ˜ t^{0.65}, in agreement with the classical experiments by Alemany etal. Figs 5, Refs 11.
FRANC2D: A two-dimensional crack propagation simulator. Version 2.7: User's guide
NASA Technical Reports Server (NTRS)
Wawrzynek, Paul; Ingraffea, Anthony
1994-01-01
FRANC 2D (FRacture ANalysis Code, 2 Dimensions) is a menu driven, interactive finite element computer code that performs fracture mechanics analyses of 2-D structures. The code has an automatic mesh generator for triangular and quadrilateral elements. FRANC2D calculates the stress intensity factor using linear elastic fracture mechanics and evaluates crack extension using several methods that may be selected by the user. The code features a mesh refinement and adaptive mesh generation capability that is automatically developed according to the predicted crack extension direction and length. The code also has unique features that permit the analysis of layered structure with load transfer through simulated mechanical fasteners or bonded joints. The code was written for UNIX workstations with X-windows graphics and may be executed on the following computers: DEC DecStation 3000 and 5000 series, IBM RS/6000 series, Hewlitt-Packard 9000/700 series, SUN Sparc stations, and most Silicon Graphics models.
The Contribution of Jets to Coronal and Solar Wind Energetics: MHD Simulations
NASA Astrophysics Data System (ADS)
Lionello, Roberto; Torok, Tibor; Titov, Viacheslav; Linker, Jon A.; Mikic, Zoran; Leake, James E.; Linton, Mark
2016-05-01
Transient collimated plasma eruptions in the corona, commonly known as coronal jets, are among the most interesting manifestations of solar activity.We use the 3D MHD model with thermodynamics developed at PSI to investigate the origin, dynamics, and plasma properties of coronal jets.Our model is coupled with 3D MHD flux emergence simulations, i.e, we use boundary conditions provided by such simulations to drive a time-dependent coronal evolution. It includes parametric coronal heating, radiative losses, and thermal conduction in the energy equations.This enables us to simulate the energy transfer in coronal jets in a more realistic manner than done so far and to study the amount of energy and mass transported by these phenomena into the higher corona and inner heliosphere. We discuss our results and compare them with previous estimations obtained from observations.
NASA Astrophysics Data System (ADS)
Ando, Tsutomu; Ueno, Kazuyuki; Sawada, Keisuke
Numerical simulation at the same condition as an experiment is carried out under the magnetic Stokes approximation for small shielding parameter. Results of the simulation compensate for the information of molten metal flow that we could not directly obtain in the experiment. In this paper, we study the molten metal flow at a starting condition and quasi-steady state. Besides, the energy conversion in the MHD pump is discussed. The simulation result shows that the proposed MHD pump causes the spiral induced current in a molten gallium and produces an axial flow with swirl. At quasi-steady state, it is confirmed that the centrifugal force by the excessive swirl flow produces high pressure at a duct wall and low pressure around the central axis. Since the excessive swirl flow results in large viscous dissipation, the mechanical power output of the pump uses only about 1% of the mechanical energy production in the molten gallium.
Simulation of Cardiac Arrhythmias Using a 2D Heterogeneous Whole Heart Model
Balakrishnan, Minimol; Chakravarthy, V. Srinivasa; Guhathakurta, Soma
2015-01-01
Simulation studies of cardiac arrhythmias at the whole heart level with electrocardiogram (ECG) gives an understanding of how the underlying cell and tissue level changes manifest as rhythm disturbances in the ECG. We present a 2D whole heart model (WHM2D) which can accommodate variations at the cellular level and can generate the ECG waveform. It is shown that, by varying cellular-level parameters like the gap junction conductance (GJC), excitability, action potential duration (APD) and frequency of oscillations of the auto-rhythmic cell in WHM2D a large variety of cardiac arrhythmias can be generated including sinus tachycardia, sinus bradycardia, sinus arrhythmia, sinus pause, junctional rhythm, Wolf Parkinson White syndrome and all types of AV conduction blocks. WHM2D includes key components of the electrical conduction system of the heart like the SA (Sino atrial) node cells, fast conducting intranodal pathways, slow conducting atriovenctricular (AV) node, bundle of His cells, Purkinje network, atrial, and ventricular myocardial cells. SA nodal cells, AV nodal cells, bundle of His cells, and Purkinje cells are represented by the Fitzhugh-Nagumo (FN) model which is a reduced model of the Hodgkin-Huxley neuron model. The atrial and ventricular myocardial cells are modeled by the Aliev-Panfilov (AP) two-variable model proposed for cardiac excitation. WHM2D can prove to be a valuable clinical tool for understanding cardiac arrhythmias. PMID:26733873
Simulation of Cardiac Arrhythmias Using a 2D Heterogeneous Whole Heart Model.
Balakrishnan, Minimol; Chakravarthy, V Srinivasa; Guhathakurta, Soma
2015-01-01
Simulation studies of cardiac arrhythmias at the whole heart level with electrocardiogram (ECG) gives an understanding of how the underlying cell and tissue level changes manifest as rhythm disturbances in the ECG. We present a 2D whole heart model (WHM2D) which can accommodate variations at the cellular level and can generate the ECG waveform. It is shown that, by varying cellular-level parameters like the gap junction conductance (GJC), excitability, action potential duration (APD) and frequency of oscillations of the auto-rhythmic cell in WHM2D a large variety of cardiac arrhythmias can be generated including sinus tachycardia, sinus bradycardia, sinus arrhythmia, sinus pause, junctional rhythm, Wolf Parkinson White syndrome and all types of AV conduction blocks. WHM2D includes key components of the electrical conduction system of the heart like the SA (Sino atrial) node cells, fast conducting intranodal pathways, slow conducting atriovenctricular (AV) node, bundle of His cells, Purkinje network, atrial, and ventricular myocardial cells. SA nodal cells, AV nodal cells, bundle of His cells, and Purkinje cells are represented by the Fitzhugh-Nagumo (FN) model which is a reduced model of the Hodgkin-Huxley neuron model. The atrial and ventricular myocardial cells are modeled by the Aliev-Panfilov (AP) two-variable model proposed for cardiac excitation. WHM2D can prove to be a valuable clinical tool for understanding cardiac arrhythmias. PMID:26733873
MHD simulation of the solar wind interaction with the magnetosphere of Mercury
NASA Astrophysics Data System (ADS)
Varela, Jacobo; Pantellini, Filippo; Moncuquet, Michel
2014-05-01
We show MHD simulations of the solar wind interaction with the magnetosphere of Mercury. We use the open source codes Pluto and MPI-AMRVAC in 3 dimensional spherical geometry. In order to appreciate the limits of the MHD approach in the context of Mercury's environment we do first compare our simulations with hybrid simulation (e.g. Trávníček et al, Icarus, 209, pp 11-22, 2010). We do also compare magnetic field profiles from the magnetometer on Messenger with profiles sampled along the corresponding spacecraft trajectory in the simulations. These comparisons show that despite the lack of kinetic effects, MHD simulation provide a more than fair description of the interaction of the solar wind with Mercury at low computational cost making it a useful tool to help decrypt data from current and future exploratory missions in the hermean magnetosphere (e.g. Bepi Colombo-MMO). The research leading to these results has received funding from the European Commission's Seventh Framework Programme (FP7/2007-2013) under the grant agreement SHOCK (project number 284515).
Test-particle Orbit Simulations in Fields from a Realistic 3D MHD Simulation
NASA Astrophysics Data System (ADS)
Decker, R. B.; Opher, M.; Hill, M. E.
2007-05-01
Models designed to explore the global structure of the heliosphere have become increasing sophisticated. Incentives to increase and to further explore the predictive capabilities of such models include the entry of the Voyager spacecraft into the foreshock region of the termination shock (TS), Voyager 1 in mid-2002 and Voyager 2 in late 2004, and the crossing of the TS and passage into the heliosheath (HSH) of Voyager 1 in 2004 day 351. Using the electric and magnetic fields generated by a MHD model of a 3D, asymmetric heliosphere [Opher et al., Ap. J. L., 640, 2006], we have developed full-particle and adiabatic-orbit codes to simulate the motion of test particles in the solar wind, TS, and HSH environments. The full-particle orbits are necessary to investigate energetic ion (e.g., anomalous and galactic cosmic ray) motion at the TS and within the heliospheric current sheet that is included in the MHD model. Adiabatic orbits are used to study particle motion in the much larger volume of the HSH where the non-homogeneous model fields produce complex guiding center motions, including mirroring in local field compressions. We will present results from these orbit computations, which are intended to provide an initial, albeit simplified, look at the propagation of high-energy charged particles, in the scatter-free limit, in the best model of the TS/HSH field configurations currently available. We will also display drift paths of high-energy ions in the HSH fields using the guiding center drift equations that are applicable in the limit of diffusive propagation.
NASA Technical Reports Server (NTRS)
Tang, H. T.; Hofmann, R.; Yee, G.; Vaughan, D. K.
1980-01-01
Transient, nonlinear soil-structure interaction simulations of an Electric Power Research Institute, SIMQUAKE experiment were performed using the large strain, time domain STEALTH 2D code and a cyclic, kinematically hardening cap soil model. Results from the STEALTH simulations were compared to identical simulations performed with the TRANAL code and indicate relatively good agreement between all the STEALTH and TRANAL calculations. The differences that are seen can probably be attributed to: (1) large (STEALTH) vs. small (TRANAL) strain formulation and/or (2) grid discretization differences.
Effects of the driving mechanism in MHD simulations of coronal mass ejections
NASA Astrophysics Data System (ADS)
Linker, J. A.; van Hoven, G.; Schnack, D. D.
Results of time-dependent MHD simulations of mass ejections in the solar coronal are presented. Previous authors have shown that results from simulations using a thermal driving mechanism are consistent with the observations only if an elaborate model of the initial corona is used. The first simulation effort, using a simple model of a plasmoid as the driving mechanism and a simple model of the initial corona, produces results that are also consistent with many observational features, suggesting that the nature of the driving mechanism plays an important role in determining the subsequent evolution of mass ejections. First simulations are based on the assumption that mass ejections are driven by magnetic forces.
Effects of the driving mechanism in MHD simulations of coronal mass ejections
NASA Technical Reports Server (NTRS)
Linker, J. A.; Van Hoven, G.; Schnack, D. D.
1990-01-01
Results of time-dependent MHD simulations of mass ejections in the solar coronal are presented. Previous authors have shown that results from simulations using a thermal driving mechanism are consistent with the observations only if an elaborate model of the initial corona is used. The first simulation effort, using a simple model of a plasmoid as the driving mechanism and a simple model of the initial corona, produces results that are also consistent with many observational features, suggesting that the nature of the driving mechanism plays an important role in determining the subsequent evolution of mass ejections. First simulations are based on the assumption that mass ejections are driven by magnetic forces.
Momentum Transport: 2D and 3D Cloud Resolving Model Simulations
NASA Technical Reports Server (NTRS)
Tao, Wei-Kuo
2001-01-01
The major objective of this study is to investigate the momentum budgets associated with several convective systems that developed during the TOGA COARE IOP (west Pacific warm pool region) and GATE (east Atlantic region). The tool for this study is the improved Goddard Cumulas Ensemble (GCE) model which includes a 3-class ice-phase microphysical scheme, explicit cloud radiative interactive processes and air-sea interactive surface processes. The model domain contains 256 x 256 grid points (with 2 km resolution) in the horizontal and 38 grid points (to a depth of 22 km) in the vertical. The 2D domain has 1024 grid points. The simulations were performed over a 7-day time period (December 19-26, 1992, for TOGA COARE and September 1-7, 1994 for GATE). Cyclic literal boundary conditions are required for this type of long-term integration. Two well organized squall systems (TOGA, COARE February 22, 1993, and GATE September 12, 1994) were also simulated using the 3D GCE model. Only 9 h simulations were required to cover the life time of the squall systems. the lateral boundary conditions were open for these two squall systems simulations. the following will be examined: (1) the momentum budgets in the convective and stratiform regions, (2) the relationship between momentum transport and cloud organization (i.e., well organized squall lines versus less organized convective), (3) the differences and similarities in momentum transport between 2D and 3D simulated convective systems, and (4) the differences and similarities in momentum budgets between cloud systems simulated with open and cyclic lateral boundary conditions. Preliminary results indicate that there are only small differences between 2D and 3D simulated momentum budgets. Major differences occur, however, between momentum budgets associated with squall systems simulated using different lateral boundary conditions.
Numerical Simulation of 3-D Supersonic Viscous Flow in an Experimental MHD Channel
NASA Technical Reports Server (NTRS)
Kato, Hiromasa; Tannehill, John C.; Gupta, Sumeet; Mehta, Unmeel B.
2004-01-01
The 3-D supersonic viscous flow in an experimental MHD channel has been numerically simulated. The experimental MHD channel is currently in operation at NASA Ames Research Center. The channel contains a nozzle section, a center section, and an accelerator section where magnetic and electric fields can be imposed on the flow. In recent tests, velocity increases of up to 40% have been achieved in the accelerator section. The flow in the channel is numerically computed using a new 3-D parabolized Navier-Stokes (PNS) algorithm that has been developed to efficiently compute MHD flows in the low magnetic Reynolds number regime. The MHD effects are modeled by introducing source terms into the PNS equations which can then be solved in a very e5uent manner. To account for upstream (elliptic) effects, the flowfield can be computed using multiple streamwise sweeps with an iterated PNS algorithm. The new algorithm has been used to compute two test cases that match the experimental conditions. In both cases, magnetic and electric fields are applied to the flow. The computed results are in good agreement with the available experimental data.
Constrained-transport Hall-MHD simulations using CWENO reconstruction with libMRC
NASA Astrophysics Data System (ADS)
Lin, Liwei; Germaschewski, Kai; Abbott, Stephen; Maynard, Kris; Raeder, Jimmy
2013-10-01
We present a new CWENO (Centrally-Weighted Essentially Non-Oscillatory) reconstruction based extended MHD (XMHD) solver that has been built for libMRC. libMRC is a library for creating efficient parallel PDE solvers on structured grids, which is used in the MRC (Magnetic Reconnection Code), OpenGGCM (Open Global Geospace Circulation Model) and PSC (Plasma Simulation Code) codes. The use of libMRC gives us access to its core functionality of providing an automated code generation framework which takes a user provided PDE right hand side in symbolic form to generate an efficient, computer-architecture specific, parallel code. libMRC also supports block-structured adaptive mesh refinement, and implicit-time stepping through integration with the PETSc library. We demonstrate validation of the new CWENO MHD solver against existing solvers both in standard test problems as well as in 3D global magnetosphere simulations.
Direct simulation of multi-phase MHD flows on an unstructured Cartesian adaptive system
NASA Astrophysics Data System (ADS)
Zhang, Jie; Ni, Ming-Jiu
2014-08-01
An approach for direct simulation of the multi-phase magnetohydrodynamics (MHD) flows has been developed in the present study on an unstructured Cartesian adaptive system. The approach is based on the volume-of-fluid (VOF) method for capturing the interface with the adaptive mesh refinement (AMR) technique used to well resolve the interface and the boundary layer. The Lorentz force is calculated using the consistent and conservative scheme, which is specially designed on a Cartesian adaptive mesh to conserve the physical conservation laws. The continuous-surface-tension (CSF) formulation is adopted for surface tension calculation. Moreover, the interfacial flows driven by thermal Marangoni effects at multifluid interfaces are also studied with a special numerical treatment presented. The method is able to simulate bubble motion in liquid metal under magnetic field irrespective of high density ratio and electric conductivity ratio. The proposed scheme for multi-phase MHD flows is validated by experimental results as well as analytical solutions.
Substorm effects in MHD and test particle simulations of magnetotail dynamics
Birn, J.; Hesse, M.
1998-12-31
Recent magnetohydrodynamic simulations demonstrate that a global tail instability, initiated by localized breakdown of MHD, can cause plasmoid formation and ejection as well as dipolarization and the current diversion of the substorm current wedge. The connection between the reconnection process and the current wedge signatures is provided by earthward flow from the reconnection site. Its braking and diversion in the inner magnetosphere causes dipolarization and the magnetic field distortions of the current wedge. The authors demonstrate the characteristic properties of this process and the current systems involved. The strong localized electric field associated with the flow burst and the dipolarization is also the cause of particle acceleration and energetic particle injections. Test particle simulations of orbits in the MHD fields yield results that are quite consistent with observed injection signatures.
3D simulations of fluctuation spectra in the hall-MHD plasma.
Shaikh, Dastgeer; Shukla, P K
2009-01-30
Turbulent spectral cascades are investigated by means of fully three-dimensional (3D) simulations of a compressible Hall-magnetohydrodynamic (H-MHD) plasma in order to understand the observed spectral break in the solar wind turbulence spectra in the regime where the characteristic length scales associated with electromagnetic fluctuations are smaller than the ion gyroradius. In this regime, the results of our 3D simulations exhibit that turbulent spectral cascades in the presence of a mean magnetic field follow an omnidirectional anisotropic inertial-range spectrum close to k(-7/3). The latter is associated with the Hall current arising from nonequal electron and ion fluid velocities in our 3D H-MHD plasma model. PMID:19257431
Solar Wind Turbulence from MHD to Sub-ion Scales: High-resolution Hybrid Simulations
NASA Astrophysics Data System (ADS)
Franci, Luca; Verdini, Andrea; Matteini, Lorenzo; Landi, Simone; Hellinger, Petr
2015-05-01
We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wavenumbers. The simulation results simultaneously exhibit several properties of the observed solar wind fluctuations: spectral indices of the magnetic, kinetic, and residual energy spectra in the magnetohydrodynamic (MHD) inertial range along with a flattening of the electric field spectrum, an increase in magnetic compressibility, and a strong coupling of the cascade with the density and the parallel component of the magnetic fluctuations at sub-proton scales. Our findings support the interpretation that in the solar wind, large-scale MHD fluctuations naturally evolve beyond proton scales into a turbulent regime that is governed by the generalized Ohm’s law.
Gas Core Reactor Numerical Simulation Using a Coupled MHD-MCNP Model
NASA Technical Reports Server (NTRS)
Kazeminezhad, F.; Anghaie, S.
2008-01-01
Analysis is provided in this report of using two head-on magnetohydrodynamic (MHD) shocks to achieve supercritical nuclear fission in an axially elongated cylinder filled with UF4 gas as an energy source for deep space missions. The motivation for each aspect of the design is explained and supported by theory and numerical simulations. A subsequent report will provide detail on relevant experimental work to validate the concept. Here the focus is on the theory of and simulations for the proposed gas core reactor conceptual design from the onset of shock generations to the supercritical state achieved when the shocks collide. The MHD model is coupled to a standard nuclear code (MCNP) to observe the neutron flux and fission power attributed to the supercritical state brought about by the shock collisions. Throughout the modeling, realistic parameters are used for the initial ambient gaseous state and currents to ensure a resulting supercritical state upon shock collisions.
Simulation of surface tension in 2D and 3D with smoothed particle hydrodynamics method
NASA Astrophysics Data System (ADS)
Zhang, Mingyu
2010-09-01
The methods for simulating surface tension with smoothed particle hydrodynamics (SPH) method in two dimensions and three dimensions are developed. In 2D surface tension model, the SPH particle on the boundary in 2D is detected dynamically according to the algorithm developed by Dilts [G.A. Dilts, Moving least-squares particle hydrodynamics II: conservation and boundaries, International Journal for Numerical Methods in Engineering 48 (2000) 1503-1524]. The boundary curve in 2D is reconstructed locally with Lagrangian interpolation polynomial. In 3D surface tension model, the SPH particle on the boundary in 3D is detected dynamically according to the algorithm developed by Haque and Dilts [A. Haque, G.A. Dilts, Three-dimensional boundary detection for particle methods, Journal of Computational Physics 226 (2007) 1710-1730]. The boundary surface in 3D is reconstructed locally with moving least squares (MLS) method. By transforming the coordinate system, it is guaranteed that the interface function is one-valued in the local coordinate system. The normal vector and curvature of the boundary surface are calculated according to the reconstructed boundary surface and then surface tension force can be calculated. Surface tension force acts only on the boundary particle. Density correction is applied to the boundary particle in order to remove the boundary inconsistency. The surface tension models in 2D and 3D have been applied to benchmark tests for surface tension. The ability of the current method applying to the simulation of surface tension in 2D and 3D is proved.
Modeling and 2-D discrete simulation of dislocation dynamics for plastic deformation of metal
NASA Astrophysics Data System (ADS)
Liu, Juan; Cui, Zhenshan; Ou, Hengan; Ruan, Liqun
2013-05-01
Two methods are employed in this paper to investigate the dislocation evolution during plastic deformation of metal. One method is dislocation dynamic simulation of two-dimensional discrete dislocation dynamics (2D-DDD), and the other is dislocation dynamics modeling by means of nonlinear analysis. As screw dislocation is prone to disappear by cross-slip, only edge dislocation is taken into account in simulation. First, an approach of 2D-DDD is used to graphically simulate and exhibit the collective motion of a large number of discrete dislocations. In the beginning, initial grains are generated in the simulation cells according to the mechanism of grain growth and the initial dislocation is randomly distributed in grains and relaxed under the internal stress. During the simulation process, the externally imposed stress, the long range stress contribution of all dislocations and the short range stress caused by the grain boundaries are calculated. Under the action of these forces, dislocations begin to glide, climb, multiply, annihilate and react with each other. Besides, thermal activation process is included. Through the simulation, the distribution of dislocation and the stress-strain curves can be obtained. On the other hand, based on the classic dislocation theory, the variation of the dislocation density with time is described by nonlinear differential equations. Finite difference method (FDM) is used to solve the built differential equations. The dislocation evolution at a constant strain rate is taken as an example to verify the rationality of the model.
The influence of slope profile extraction techniques and DEM resolution on 2D rockfall simulation
NASA Astrophysics Data System (ADS)
Wang, X.; Frattini, P.; Agliardi, F.; Crosta, G. B.
2012-04-01
The development of advanced 3D rockfall modelling algorithms and tools during the last decade has allowed to gain insights in the topographic controls on the quality and reliability of rockfall simulation results. These controls include DEM resolution and roughness, and depend on the adopted rockfall simulation approach and DEM generation techniques. Despite the development of 3D simulations, the 2D modelling approach still remains suitable and convenient in some cases. Therefore, the accuracy of high-quality 3D descriptions of topography must be preserved when extracting slope profiles for 2D simulations. In this perspective, this study compares and evaluates three different techniques commonly used to extract slope profiles from DEM, in order to assess their suitability and effects on rockfall simulation results. These methods include: (A) an "interpolated shape" method (ESRI 3D Analyst), (B) a raw raster sampling method (EZ Profiler), and (C) a vector TIN sampling method (ESRI 3D Analyst). The raster DEMs used in the study were all derived from the same TIN DEM used for method C. For raster DEM, the "interpolated shape" method (A) extracts the profile by bi-linear interpolating the elevation among the four neighbouring cells at each sampling location along the profile trace. The EZ Profiler extension (B) extracts the profile by sampling elevation values directly from the DEM raster grid at each sampling location. These methods have been compared to the extraction of profiles from TIN DEM (C), where slope profile elevations are directly obtained by sampling the TIN triangular facets. 2D rockfall simulations performed using a widely used commercial software (RocfallTM) with the different profiles show that: (1) method A and C provide similar results; (2) runout simulated using profiles obtained by method A is usually shorter than method C; (3) method B presents abrupt horizontal steps in the profiles, resulting in unrealistic runout. To study the influence of DEM
Thermodynamic MHD Simulation of the 2000 July 14 "Bastille Day" Eruption
NASA Astrophysics Data System (ADS)
Torok, Tibor; Downs, Cooper; Lionello, Roberto; Linker, Jon A.; Titov, Viacheslav S.; Mikic, Zoran; Riley, Pete
2015-04-01
The "Bastille Day" event that occurred on 2000 July 14 is one of the most extensively studied solar eruptions. It originated in a complex active region close to disk center and produced an X5.7 flare, a fast halo CME, and an intense geomagnetic storm. Accurate numerical simulations of such events, in particular the matching of parameters relevant for space weather such as the CME velocity and magnetic orientation, require a realistic model of the large-scale magnetic field and plasma environment into which the eruption propagates and interacts, as well as a modeling of the pre-eruptive configuration and eruption initiation that are as realistic as possible. Here we present an MHD simulation of the Bastille Day event that complies with these requirements. We first produce a steady-state MHD solution of the background corona that incorporates realistic energy transport ("thermodynamic MHD"), photospheric magnetic field measurements, and the solar wind. In order to model the pre-eruptive magnetic field, we then insert a stable, elongated flux rope that resides above the highly curved polarity inversion line of the active region. Finally, we produce the eruption by imposing photospheric flows that slowly converge towards the polarity inversion line. In this presentation we describe our method, compare the simulation results with the observations, and discuss the challenges and limitations involved in modeling such complex and powerful eruptions.
A three-dimensional MHD simulation analysis of the origin of the slow solar wind
NASA Astrophysics Data System (ADS)
Washimi, H.; Zank, G. P.; Hu, Q.; Nakamizo, A.; Tanaka, T.; Kojima, M.; Kubo, Y.
2012-12-01
We have developed a 3D MHD simulation model for the study of the solar-wind acceleration mechanism and for reproducing a realistic configuration of solar wind plasma by using observed photospheric magnetic field at each Carrington rotation cycle. Using an unstructured mesh coordinate system on spherical surface with fine spacing in radial direction, we aim to reproduce a wide range of solar-wind plasma configuration from the photosphere to 1AU. We have incorporated external source terms into the momentum and energy equations in our MHD simulation. The energy source term consists of two volumetric heating functions: one is a new term, as a new development from our original model (Nakamizo et al. JGR 114, A07109, 2009), for the heating in a very narrow region around the transition region. The other one is an additional heating source which probably comes from some nonlinear wave phenomena which are effective over a radial distance of an order of the solar radius in the corona. The Spitzer-type thermal conduction term is also taken into account. The momentum source term is given in a form similar to that of the nonlinear wave heating function noted above. Using this MHD simulation system, we will study the origin of the slow solar wind from nearby regions of some isolated active regions during CR1900-CR1913 in some details.
Quantum simulation of 2D topological physics in a 1D array of optical cavities
Luo, Xi-Wang; Zhou, Xingxiang; Li, Chuan-Feng; Xu, Jin-Shi; Guo, Guang-Can; Zhou, Zheng-Wei
2015-01-01
Orbital angular momentum of light is a fundamental optical degree of freedom characterized by unlimited number of available angular momentum states. Although this unique property has proved invaluable in diverse recent studies ranging from optical communication to quantum information, it has not been considered useful or even relevant for simulating nontrivial physics problems such as topological phenomena. Contrary to this misconception, we demonstrate the incredible value of orbital angular momentum of light for quantum simulation by showing theoretically how it allows to study a variety of important 2D topological physics in a 1D array of optical cavities. This application for orbital angular momentum of light not only reduces required physical resources but also increases feasible scale of simulation, and thus makes it possible to investigate important topics such as edge-state transport and topological phase transition in a small simulator ready for immediate experimental exploration. PMID:26145177
A faster method for 3D/2D medical image registration--a simulation study.
Birkfellner, Wolfgang; Wirth, Joachim; Burgstaller, Wolfgang; Baumann, Bernard; Staedele, Harald; Hammer, Beat; Gellrich, Niels Claudius; Jacob, Augustinus Ludwig; Regazzoni, Pietro; Messmer, Peter
2003-08-21
3D/2D patient-to-computed-tomography (CT) registration is a method to determine a transformation that maps two coordinate systems by comparing a projection image rendered from CT to a real projection image. Iterative variation of the CT's position between rendering steps finally leads to exact registration. Applications include exact patient positioning in radiation therapy, calibration of surgical robots, and pose estimation in computer-aided surgery. One of the problems associated with 3D/2D registration is the fact that finding a registration includes solving a minimization problem in six degrees of freedom (dof) in motion. This results in considerable time requirements since for each iteration step at least one volume rendering has to be computed. We show that by choosing an appropriate world coordinate system and by applying a 2D/2D registration method in each iteration step, the number of iterations can be grossly reduced from n6 to n5. Here, n is the number of discrete variations around a given coordinate. Depending on the configuration of the optimization algorithm, this reduces the total number of iterations necessary to at least 1/3 of it's original value. The method was implemented and extensively tested on simulated x-ray images of a tibia, a pelvis and a skull base. When using one projective image and a discrete full parameter space search for solving the optimization problem, average accuracy was found to be 1.0 +/- 0.6(degrees) and 4.1 +/- 1.9 (mm) for a registration in six parameters, and 1.0 +/- 0.7(degrees) and 4.2 +/- 1.6 (mm) when using the 5 + 1 dof method described in this paper. Time requirements were reduced by a factor 3.1. We conclude that this hardware-independent optimization of 3D/2D registration is a step towards increasing the acceptance of this promising method for a wide number of clinical applications. PMID:12974581
Broadband Electron Precipitation in Global MHD Simulation and its Effect on the Ionosphere
NASA Astrophysics Data System (ADS)
Zhang, B.; Lotko, W.; Brambles, O. J.; Wiltberger, M. J.
2010-12-01
A broadband electron (BBE) precipitation model is implemented and analyzed in the MI coupling module of the Lyon-Fedder-Mobarry MHD simulation. Both number flux and energy flux of precipitating BBEs are regulated by MHD variables calculated near the low-altitude boundary of the LFM simulation. An empirical relation deduced from results of Keiling et al. (2003) is used to relate the AC Poynting flux to the energy flux precipitating BBEs in the simulation. We are investigating two different ways of regulating the number flux of BBE precipitation, one using an empirical relation between AC Poynting flux and number flux (Strangeway, unpublished) and another by constraining the intensity and cut-off energy of a fixed-pitch angle distribution of BBEs in terms of MHD simulation variables. The contributions to ionospheric conductance from BBE precipitation are evaluated using empirical relations derived by Robinson et al. (1987). The BBE-induced-conductance is added to the “standard” auroral contribution to conductance derived from monoenergetic and diffuse electron precipitation in the existing LFM precipitation model. The simulation is driven by ideal SW/IMF conditions with Vsw=400 km/s, Nsw=5/cc and Bz=-5 nT. The simulated time-average AC Poynting flux pattern resembles statistical patterns from Polar data (Keiling et al. 2003), and the simulated statistical pattern of BBE number flux resembles the statistical maps derived from DMSP data (Newell et al. 2009) on the nightside with a similar dawn-dusk asymmetry. The ionospheric Pedersen and Hall conductances are enhanced about 20% by the BBE precipitation. The number flux produced by BBEs is the same order of magnitude as that of monoenergetic and diffuse electrons. We thus expect BBE precipitation to have a moderate effect on the E-region ionosphere and a more significant influence on the density distribution of the F-region ionosphere.
MHD-PIC interlocked simulation model in space plasma
NASA Astrophysics Data System (ADS)
Sugiyama, T.; Kuasano, K.
2008-12-01
We have developed a new type of simulation technique by directly interlocking a traditional Ion-Particle Hybrid simulation model (Hybrid) and an Energetic-Particle Hybrid simulation (EP-HYB) model. In the traditional Hybrid model, all ions are kinetically treated as particles. In the EP-HYB model, non-thermal energetic ions are kinetically treated, and the thermal component is calculated as a fluid. The interlocked model is applied to a two-dimensional collisionless shock problem. The domain for the Hybrid model is embedded in a part of the system, and the bounded data are exchanged to each other to keep the consistency between both models. It can handle the full ion kinetics to investigate the injection problem at the shock transition region, as well as the wave-particle interactions in even far upstream region. We have carried out the long-term simulation of the shock acceleration process using this interlocked model, and successfully reproduced the power-law distribution function, which is consistent with the diffusive acceleration theory. Since the calculation cost of the EP-HYB model is much smaller than that of the Hybrid model, we can considerably reduce the computational demand.
Comparison between 2D and 3D Numerical Modelling of a hot forging simulative test
Croin, M.; Ghiotti, A.; Bruschi, S.
2007-04-07
The paper presents the comparative analysis between 2D and 3D modelling of a simulative experiment, performed in laboratory environment, in which operating conditions approximate hot forging of a turbine aerofoil section. The plane strain deformation was chosen as an ideal case to analyze the process because of the thickness variations in the final section and the consequent distributions of contact pressure and sliding velocity at the interface that are closed to the conditions of the real industrial process. In order to compare the performances of 2D and 3D approaches, two different analyses were performed and compared with the experiments in terms of loads and temperatures peaks at the interface between the dies and the workpiece.
Simulation of the flow and mass transfer for KDP crystals undergoing 2D translation during growth
NASA Astrophysics Data System (ADS)
Zhou, Chuan; Li, Mingwei; Hu, Zhitao; Yin, Huawei; Wang, Bangguo; Cui, Qidong
2016-09-01
In this study, a novel motion mode for crystals during growth, i.e., 2D translation, is proposed. Numerical simulations of flow and mass transfer are conducted for the growth of large-scale potassium dihydrogen phosphate (KDP) crystals subjected to the new motion mode. Surface supersaturation and shear stress are obtained as functions of the translational velocity, distance, size, orientation of crystals. The dependence of these two parameters on the flow fields around the crystals is also discussed. The thicknesses of the solute boundary layer varied with translational velocity are described. The characteristics of solution flow and surface supersaturation distribution are summarized, where it suggests that the morphological stability of a crystal surface can be enhanced if the proposed 2D translation is applied to crystal growth.
NASA Astrophysics Data System (ADS)
Wu, S. T.; Zheng, Huinan; Wang, S.; Thompson, B. J.; Plunkett, S. P.; Zhao, X. P.; Dryer, M.
2001-11-01
We investigate the global large amplitude waves propagating across the solar disk as observed by the SOHO/Extreme Ultraviolet Imaging Telescope (EIT). These waves appear to be similar to those observed in Hα in the chromosphere and which are known as ``Moreton waves,'' associated with large solar flares [Moreton, 1960, 1964]. Uchida [1968] interpreted these Moreton waves as the propagation of a hydromagnetics disturbance in the corona with its wavefront intersecting the chromosphere to produce the Moreton wave as observed in movie sequences of Hα images. To search for an understanding of the physical characteristics of these newly observed EIT waves, we constructed a three-dimensional, time-dependent, numerical magnetohydrodynamic (MHD) model. Measured global magnetic fields, obtained from the Wilcox Solar Observatory (WSO) at Stanford University, are used as the initial magnetic field to investigate hydromagnetics wave propagation in a three-dimensional spherical geometry. Using magnetohydrodynamic wave theory together with simulation, we are able to identify these observed EIT waves as fast mode MHD waves dominated by the acoustic mode, called magnetosonic waves. The results to be presented include the following: (1) comparison of observed and simulated morphology projected on the disk and the distance-time curves on the solar disk; (2) three-dimensional evolution of the disturbed magnetic field lines at various viewing angles; (3) evolution of the plasma density profile at a specific location as a function of latitude; and (4) computed Friedrich's diagrams to identify the MHD wave characteristics.
Two-dimensional magnetohydrodynamic simulations of poloidal flows in tokamaks and MHD pedestal
Guazzotto, L.; Betti, R.
2011-09-15
Poloidal rotation is routinely observed in present-day tokamak experiments, in particular near the plasma edge and in the high-confinement mode of operation. According to the magnetohydrodynamic (MHD) equilibrium theory [R. Betti and J. P. Freidberg, Phys. Plasmas 7, 2439 (2000)], radial discontinuities form when the poloidal velocity exceeds the poloidal sound speed (or rather, more correctly, the poloidal magneto-slow speed). Two-dimensional compressible magnetohydrodynamic simulations show that the transonic discontinuities develop on a time scale of a plasma poloidal revolution to form an edge density pedestal and a localized velocity shear layer at the pedestal location. While such an MHD pedestal surrounds the entire core, the outboard side of the pedestal is driven by the transonic discontinuity while the inboard side is caused by a poloidal redistribution of the mass. The MHD simulations use a smooth momentum source to drive the poloidal flow. Soon after the flow exceeds the poloidal sound speed, the density pedestal and the velocity shear layer form and persist into a quasi steady state. These results may be relevant to the L-H transition, the early stages of the pedestal and edge transport barrier formation.
PROPERTIES OF UMBRAL DOTS AS MEASURED FROM THE NEW SOLAR TELESCOPE DATA AND MHD SIMULATIONS
Kilcik, A.; Yurchyshyn, V. B.; Abramenko, V.; Goode, P. R.; Cao, W.; Rempel, M.; Kitai, R.; Watanabe, H.
2012-02-01
We studied bright umbral dots (UDs) detected in a moderate size sunspot and compared their statistical properties to recent MHD models. The study is based on high-resolution data recorded by the New Solar Telescope at the Big Bear Solar Observatory and three-dimensional (3D) MHD simulations of sunspots. Observed UDs, living longer than 150 s, were detected and tracked in a 46 minute long data set, using an automatic detection code. A total of 1553 (620) UDs were detected in the photospheric (low chromospheric) data. Our main findings are (1) none of the analyzed UDs is precisely circular, (2) the diameter-intensity relationship only holds in bright umbral areas, and (3) UD velocities are inversely related to their lifetime. While nearly all photospheric UDs can be identified in the low chromospheric images, some small closely spaced UDs appear in the low chromosphere as a single cluster. Slow-moving and long-living UDs seem to exist in both the low chromosphere and photosphere, while fast-moving and short-living UDs are mainly detected in the photospheric images. Comparison to the 3D MHD simulations showed that both types of UDs display, on average, very similar statistical characteristics. However, (1) the average number of observed UDs per unit area is smaller than that of the model UDs, and (2) on average, the diameter of model UDs is slightly larger than that of observed ones.
Neoclassical viscous stress tensor for non-linear MHD simulations with XTOR-2F
NASA Astrophysics Data System (ADS)
Mellet, N.; Maget, P.; Lütjens, H.; Meshcheriakov, D.; the Tore Supra Team
2013-04-01
The neoclassical viscous stress tensor is implemented in the non-linear MHD code XTOR-2F (Lütjens and Luciani 2010 J. Comput. Phys. 229 8130-43), allowing consistent bi-fluid simulations of MHD modes, including the metastable branch of neoclassical tearing modes (NTMs) (Carrera et al 1986 Phys. Fluids 29 899-902). Equilibrium flows and bootstrap current from the neoclassical theory are formally recovered in this Chew-Goldberger-Low formulation. The non-linear behaviour of the new model is verified on a test case coming from a Tore Supra non-inductive discharge. A NTM threshold that is larger than with the previous model is obtained. This is due to the fact that the velocity is now part of the bootstrap current and that it differs from the theoretical neoclassical value.
MHD simulations of large conducting bodies moving through a planetary magnetosphere
NASA Astrophysics Data System (ADS)
Kopp, A.; Schröer, A.
1998-01-01
The interaction between a large body with finite conductivity moving relatively to a planetary magnetic field and the magnetospheric plasma is studied by means of numerical simulations in the framework of resistive MHD. Applications discussed are spacecrafts (tethered satellites) in the Earth's ionosphere and the interaction of Jupiter with its satellite Io. The interaction excites MHD waves in which consequence an electrical current flow through the conductor is generated which extends into the surrounding plasma and propagates along the Alfvén characteristics. Thus, kinetic energy is converted into electrical energy which may be utilized for satellite projects or explain Jupiter's decametric radiation. The linear and nonlinear evolution of the plasma environment and the current system is investigated with special care taken of the principal phenomena occuring in the vicinity of the conducting body.
Explosive Turbulent Magnetic Reconnection: A New Approach of MHD-Turbulent Simulation
NASA Astrophysics Data System (ADS)
Hoshino, Masahiro; Yokoi, Nobumitsu; Higashimori, Katsuaki
2013-04-01
Turbulent flows are often observed in association with magnetic reconnection in space and astrophysical plasmas, and it is often hypothesized that the turbulence can contribute to the fast magnetic reconnection through the enhancement of magnetic dissipation. In this presentation, we demonstrate that an explosive turbulent reconnection can happen by using a new turbulent MHD simulation, in which the evolution of the turbulent transport coefficients are self-consistently solved together with the standard MHD equations. In our model, the turbulent electromotive force defined by the correlation of turbulent fluctuations between v and B is added to the Ohm's law. We discuss that the level of turbulent can control the topology of reconnection, namely the transition from the Sweet-Parker reconnection to the Petscheck reconnection occurs when the level of fluctuations becomes of order of the ambient physical quantities, and show that the growth of the turbulent Petscheck reconnection becomes much faster than the conventional one.
Simulation of 2D NMR Spectra of Carbohydrates Using GODESS Software.
Kapaev, Roman R; Toukach, Philip V
2016-06-27
Glycan Optimized Dual Empirical Spectrum Simulation (GODESS) is a web service, which has been recently shown to be one of the most accurate tools for simulation of (1)H and (13)C 1D NMR spectra of natural carbohydrates and their derivatives. The new version of GODESS supports visualization of the simulated (1)H and (13)C chemical shifts in the form of most 2D spin correlation spectra commonly used in carbohydrate research, such as (1)H-(1)H TOCSY, COSY/COSY-DQF/COSY-RCT, and (1)H-(13)C edHSQC, HSQC-COSY, HSQC-TOCSY, and HMBC. Peaks in the simulated 2D spectra are color-coded and labeled according to the signal assignment and can be exported in JCAMP-DX format. Peak widths are estimated empirically from the structural features. GODESS is available free of charge via the Internet at the platform of the Carbohydrate Structure Database project ( http://csdb.glycoscience.ru ). PMID:27227420
Application of 2-D simulations to hollow z-pinch implosions
Peterson, D.L.; Bowers, R.L.; Brownell, J.H.
1997-12-01
The application of simulations of z-pinch implosions should have at least two goals: first, to properly model the most important physical processes occurring in the pinch allowing for a better understanding of the experiments and second, provide a design capability for future experiments. Beginning with experiments fielded at Los Alamos on the Pegasus 1 and Pegasus 2 capacitor banks, the authors have developed a methodology for simulating hollow z-pinches in two dimensions which has reproduced important features of the measured experimental current drive, spectrum, radiation pulse shape, peak power and total radiated energy. This methodology employs essentially one free parameter, the initial level of the random density perturbations imposed at the beginning of the 2-D simulation, but in general no adjustments to other parameters are required. Currently the authors are applying this capability to the analysis of recent Saturn and PBFA-Z experiments. The code results provide insight into the nature of the pinch plasma prior to arrival on-axis, during thermalization and development after peak pinch time. Among other things, the simulation results provide an explanation for the production of larger amounts of radiated energy than would be expected from a simple slug-model kinetic energy analysis and the appearance of multiple peaks in the radiation power. The 2-D modeling has also been applied to the analysis of Saturn dynamic hohlraum experiments and is being used in the design of this and other Z-Pinch applications on PBFA-Z.
MHD simulations of coronal mass ejections - Importance of the driving mechanism
NASA Astrophysics Data System (ADS)
Linker, J. A.; van Hoven, G.; Schnack, D. D.
1990-04-01
The importance of the form of the driving mechanism in MHD simulations of coronal mass ejections is investigated. A model simulation problem is devised, and it is found that the use of a simple form for the initial corona, with an upward moving parcel of cold, dense plasma as the driving mechanism, can produce results that are consistent with many of the features observed by coronagraphs. The results imply that the nature of the driving mechanism may play an important role in determining the dynamical evolution of mass ejections.
3D MHD simulations of planet migration in turbulent stratified disks
NASA Astrophysics Data System (ADS)
Uribe, Ana; Klahr, Hubert; Flock, Mario; Henning, Thomas
2011-11-01
We performed 3D MHD numerical simulations of planet migration in stratified disks using the Godunov code PLUTO (Mignone et al. 2007). The disk is invaded by turbulence generated by the magnetorotational instability (MRI). We study the migration for planets with different mass to primary mass ratio. The migration of the low-mass planet (q=Mp/Ms=10-5) is dominated by random fluctuations in the torque and there is no defined direction of migration on timescales of 100 orbits. The intermediate-mass planet (q=Mp/Ms=10-4) can experience systematic outwards migration that was sustained for the times we were able to simulate.
THE SUBMILLIMETER BUMP IN Sgr A* FROM RELATIVISTIC MHD SIMULATIONS
Dexter, Jason; Agol, Eric; Fragile, P. Chris; McKinney, Jonathan C.
2010-07-10
Recent high resolution observations of the Galactic center black hole allow for direct comparison with accretion disk simulations. We compare two-temperature synchrotron emission models from three-dimensional, general relativistic magnetohydrodynamic simulations to millimeter observations of Sgr A*. Fits to very long baseline interferometry and spectral index measurements disfavor the monochromatic face-on black hole shadow models from our previous work. Inclination angles {<=}20{sup 0} are ruled out to 3{sigma}. We estimate the inclination and position angles of the black hole, as well as the electron temperature of the accretion flow and the accretion rate, to be i=50{sup o+35o}{sub -15}{sup o}, {xi}=-23{sup o+97o}{sub -22}{sup o}, T{sub e} = (5.4 {+-} 3.0) x 10{sup 10} K, and M-dot =5{sup +15}{sub -2}x10{sup -9} M{sub sun} yr{sup -1}, respectively, with 90% confidence. The black hole shadow is unobscured in all best-fit models, and may be detected by observations on baselines between Chile and California, Arizona, or Mexico at 1.3 mm or .87 mm either through direct sampling of the visibility amplitude or using closure phase information. Millimeter flaring behavior consistent with the observations is present in all viable models and is caused by magnetic turbulence in the inner radii of the accretion flow. The variability at optically thin frequencies is strongly correlated with that in the accretion rate. The simulations provide a universal picture of the 1.3 mm emission region as a small region near the midplane in the inner radii of the accretion flow, which is roughly isothermal and has {nu}/{nu} {sub c} {approx} 1-20, where {nu} {sub c} is the critical frequency for thermal synchrotron emission.
The Submillimeter Bump in Sgr A* from Relativistic MHD Simulations
NASA Astrophysics Data System (ADS)
Dexter, Jason; Agol, Eric; Fragile, P. Chris; McKinney, Jonathan C.
2010-07-01
Recent high resolution observations of the Galactic center black hole allow for direct comparison with accretion disk simulations. We compare two-temperature synchrotron emission models from three-dimensional, general relativistic magnetohydrodynamic simulations to millimeter observations of Sgr A*. Fits to very long baseline interferometry and spectral index measurements disfavor the monochromatic face-on black hole shadow models from our previous work. Inclination angles <=20° are ruled out to 3σ. We estimate the inclination and position angles of the black hole, as well as the electron temperature of the accretion flow and the accretion rate, to be i={50°}^{+35°}_{-15°}, ξ ={-23°}^{+97°}_{-22°}, Te = (5.4 ± 3.0) × 1010 K, and \\dot{M}=5^{+15}_{-2}× 10^{-9} M_⊙ yr^{-1}, respectively, with 90% confidence. The black hole shadow is unobscured in all best-fit models, and may be detected by observations on baselines between Chile and California, Arizona, or Mexico at 1.3 mm or .87 mm either through direct sampling of the visibility amplitude or using closure phase information. Millimeter flaring behavior consistent with the observations is present in all viable models and is caused by magnetic turbulence in the inner radii of the accretion flow. The variability at optically thin frequencies is strongly correlated with that in the accretion rate. The simulations provide a universal picture of the 1.3 mm emission region as a small region near the midplane in the inner radii of the accretion flow, which is roughly isothermal and has ν/ν c ~ 1-20, where ν c is the critical frequency for thermal synchrotron emission.
NASA Astrophysics Data System (ADS)
Hatori, Tomoharu; Ito, Atsushi M.; Nunami, Masanori; Usui, Hideyuki; Miura, Hideaki
2016-08-01
We propose a numerical method to determine the artificial viscosity in magnetohydrodynamics (MHD) simulations with adaptive mesh refinement (AMR) method, where the artificial viscosity is adaptively changed due to the resolution level of the AMR hierarchy. Although the suitable value of the artificial viscosity depends on the governing equations and the model of target problem, it can be determined by von Neumann stability analysis. By means of the new method, "level-by-level artificial viscosity method," MHD simulations of Rayleigh-Taylor instability (RTI) are carried out with the AMR method. The validity of the level-by-level artificial viscosity method is confirmed by the comparison of the linear growth rates of RTI between the AMR simulations and the simple simulations with uniform grid and uniform artificial viscosity whose resolution is the same as that in the highest level of the AMR simulation. Moreover, in the nonlinear phase of RTI, the secondary instability is clearly observed where the hierarchical data structure of AMR calculation is visualized as high resolution region floats up like terraced fields. In the applications of the method to general fluid simulations, the growth of small structures can be sufficiently reproduced, while the divergence of numerical solutions can be suppressed.
3D MHD simulation of polarized emission in SN 1006
NASA Astrophysics Data System (ADS)
Schneiter, E. M.; Velázquez, P. F.; Reynoso, E. M.; Esquivel, A.; De Colle, F.
2015-05-01
We use three-dimensional magnetohydrodynamic simulations to model the supernova remnant SN 1006. From our numerical results, we have carried out a polarization study, obtaining synthetic maps of the polarized intensity, the Stokes parameter Q, and the polar-referenced angle, which can be compared with observational results. Synthetic maps were computed considering two possible particle acceleration mechanisms: quasi-parallel and quasi-perpendicular. The comparison of synthetic maps of the Stokes parameter Q maps with observations proves to be a valuable tool to discern unambiguously which mechanism is taking place in the remnant of SN 1006, giving strong support to the quasi-parallel model.
MHD simulations of homologous and cannibalistic coronal mass ejections
NASA Astrophysics Data System (ADS)
Fan, Yuhong; Chatterjee, Piyali
2014-06-01
We present magneto-hydrodynamic simulations of the development of a homologous sequence of coronal mass ejections (CMEs) and demonstrate their so-called cannibalistic behavior. These CMEs originate from the repeated formations and partial eruptions of kink unstable flux ropes as a result of the continued emergence of a twisted flux rope across the lower boundary into a pre-existing coronal potential arcade field. The simulations show that a CME erupting into the open magnetic field created by a preceding CME has a higher speed, and therefore tends to be cannibalistic, catching up and merging with the preceding one into a single fast CME. All the CMEs attained speeds of about 1000 km/s as they exit the domain. The reformation of a twisted flux rope after each CME eruption during the sustained flux emergence can naturally explain the X-ray observations of repeated reformations of sigmoids and “sigmoid-under-cusp” configurations at a low-coronal source of homologous CMEs.
Constrained Transport vs. Divergence Cleanser Options in Astrophysical MHD Simulations
NASA Astrophysics Data System (ADS)
Lindner, Christopher C.; Fragile, P.
2009-01-01
In previous work, we presented results from global numerical simulations of the evolution of black hole accretion disks using the Cosmos++ GRMHD code. In those simulations we solved the magnetic induction equation using an advection-split form, which is known not to satisfy the divergence-free constraint. To minimize the build-up of divergence error, we used a hyperbolic cleanser function that simultaneously damped the error and propagated it off the grid. We have since found that this method produces qualitatively and quantitatively different behavior in high magnetic field regions than results published by other research groups, particularly in the evacuated funnels of black-hole accretion disks where Poynting-flux jets are reported to form. The main difference between our earlier work and that of our competitors is their use of constrained-transport schemes to preserve a divergence-free magnetic field. Therefore, to study these differences directly, we have implemented a constrained transport scheme into Cosmos++. Because Cosmos++ uses a zone-centered, finite-volume method, we can not use the traditional staggered-mesh constrained transport scheme of Evans & Hawley. Instead we must implement a more general scheme; we chose the Flux-CT scheme as described by Toth. Here we present comparisons of results using the divergence-cleanser and constrained transport options in Cosmos++.
Fast Acceleration of 2D Wave Propagation Simulations Using Modern Computational Accelerators
Wang, Wei; Xu, Lifan; Cavazos, John; Huang, Howie H.; Kay, Matthew
2014-01-01
Recent developments in modern computational accelerators like Graphics Processing Units (GPUs) and coprocessors provide great opportunities for making scientific applications run faster than ever before. However, efficient parallelization of scientific code using new programming tools like CUDA requires a high level of expertise that is not available to many scientists. This, plus the fact that parallelized code is usually not portable to different architectures, creates major challenges for exploiting the full capabilities of modern computational accelerators. In this work, we sought to overcome these challenges by studying how to achieve both automated parallelization using OpenACC and enhanced portability using OpenCL. We applied our parallelization schemes using GPUs as well as Intel Many Integrated Core (MIC) coprocessor to reduce the run time of wave propagation simulations. We used a well-established 2D cardiac action potential model as a specific case-study. To the best of our knowledge, we are the first to study auto-parallelization of 2D cardiac wave propagation simulations using OpenACC. Our results identify several approaches that provide substantial speedups. The OpenACC-generated GPU code achieved more than speedup above the sequential implementation and required the addition of only a few OpenACC pragmas to the code. An OpenCL implementation provided speedups on GPUs of at least faster than the sequential implementation and faster than a parallelized OpenMP implementation. An implementation of OpenMP on Intel MIC coprocessor provided speedups of with only a few code changes to the sequential implementation. We highlight that OpenACC provides an automatic, efficient, and portable approach to achieve parallelization of 2D cardiac wave simulations on GPUs. Our approach of using OpenACC, OpenCL, and OpenMP to parallelize this particular model on modern computational accelerators should be applicable to other computational models of wave propagation in
Global MHD modeling of resonant ULF waves: Simulations with and without a plasmasphere
NASA Astrophysics Data System (ADS)
Claudepierre, S. G.; Toffoletto, F. R.; Wiltberger, M.
2016-01-01
We investigate the plasmaspheric influence on the resonant mode coupling of magnetospheric ultralow frequency (ULF) waves using the Lyon-Fedder-Mobarry (LFM) global magnetohydrodynamic (MHD) model. We present results from two different versions of the model, both driven by the same solar wind conditions: one version that contains a plasmasphere (the LFM coupled to the Rice Convection Model, where the Gallagher plasmasphere model is also included) and another that does not (the stand-alone LFM). We find that the inclusion of a cold, dense plasmasphere has a significant impact on the nature of the simulated ULF waves. For example, the inclusion of a plasmasphere leads to a deeper (more earthward) penetration of the compressional (azimuthal) electric field fluctuations, due to a shift in the location of the wave turning points. Consequently, the locations where the compressional electric field oscillations resonantly couple their energy into local toroidal mode field line resonances also shift earthward. We also find, in both simulations, that higher-frequency compressional (azimuthal) electric field oscillations penetrate deeper than lower frequency oscillations. In addition, the compressional wave mode structure in the simulations is consistent with a radial standing wave oscillation pattern, characteristic of a resonant waveguide. The incorporation of a plasmasphere into the LFM global MHD model represents an advance in the state of the art in regard to ULF wave modeling with such simulations. We offer a brief discussion of the implications for radiation belt modeling techniques that use the electric and magnetic field outputs from global MHD simulations to drive particle dynamics.
NASA Astrophysics Data System (ADS)
Shao, X.; Sitnov, M. I.; Sharma, A. S.; Papadopoulos, K.; Goodrich, C. C.; Guzdar, P. N.; Milikh, G. M.; Wiltberger, M. J.; Lyon, J. G.
2002-05-01
Studies of the magnetosphere during substorms based on the observational data of the solar wind and the geomagnetic indices have shown clear features of phase transition-like behavior [Sitnov et al., 2000]. The global MHD simulations of the events in the Bargatze et al. [1985] database are used to study the non-equilibrium phase transition-like features of substorms. We simulated 7 intervals of total duration of 280 hours from the same data set used in Sitnov et al. [2000]. From the simulations the AL index is computed from the maximum of the westward Hall current and is referred to as the pseudo-AL index. We analyzed the input-output (vBs-pseudo-AL index) system obtained from the global MHD model and compare the results to those in Sitnov et al. [2000, 2001]. The analysis of the coupled vBs-pseudo-AL index system shows the first-order phase transition characterizing global beahavior, similar to the case of vBs-observed-AL index [Sitnov et al., 2000]. Although, the comparison between observations and global MHD simulations for individual events may vary, the overall global transition pattern during the substorm cycle revealed by singular spectrum snalysis is statistically consistent between simulations and observations. The coupled vBs-pseudo-AL index system shows multi-scale behavior (scale-invariant power-law dependence) in singular power spectrum. We find critical exponents of the non-equilibrium transitions in the magnetosphere, which reflect the multi-scale aspect of the substorm activity, different from power-law frequency of autonomous systems. The exponents relate input and output parameters of the magnetosphere.
2D PIC/MC simulations of electrical asymmetry effect in capacitive coupled plasma
NASA Astrophysics Data System (ADS)
Zhang, Quan-Zhi; Jiang, Wei; Wang, You-Nian
2011-10-01
Recently a so-called electrical asymmetry effect (EAE), which could achieve high-degree separate control of ion flux and energy in dual-frequency capacitively coupled plasmas, was discovered theoretically by Heil et al. and was confirmed by experiments and theory/numerical simulations later on. However, since there always is a bigger grounded surface area for experiment devices, which reduces the geometrical symmetry, and all the simulations were limited to 1D before, it is, thus, worth studying the EAE when coupling the electrically and geometrically asymmetric discharges theoretically. Here, we perform 2D PIC/MC simulations, which can include both electrically and geometrically asymmetric factors. The EAE on plasma parameters, such as dc self-bias voltage, density profiles, ion energy distribution and power absorption of electron have been examined for different pressures and geometry conditions. Recently a so-called electrical asymmetry effect (EAE), which could achieve high-degree separate control of ion flux and energy in dual-frequency capacitively coupled plasmas, was discovered theoretically by Heil et al. and was confirmed by experiments and theory/numerical simulations later on. However, since there always is a bigger grounded surface area for experiment devices, which reduces the geometrical symmetry, and all the simulations were limited to 1D before, it is, thus, worth studying the EAE when coupling the electrically and geometrically asymmetric discharges theoretically. Here, we perform 2D PIC/MC simulations, which can include both electrically and geometrically asymmetric factors. The EAE on plasma parameters, such as dc self-bias voltage, density profiles, ion energy distribution and power absorption of electron have been examined for different pressures and geometry conditions. This work was supported by the National Natural Science Foundation of China (Grant No 10635010) and the Important National Science & Technology Specific Project (Grant No
Phase Transitions in Quasi-2D Plasma-Dust Systems: Simulations and Experiments
NASA Astrophysics Data System (ADS)
Petrov, Oleg; Vasiliev, Mikhail; Statsenko, Konstantin; Koss, Xeniya; Vasilieva, Elena; Myasnikov, Maxim; Lisin, Evgeny
2015-11-01
A nature of phase transition in quasi-2D dusty plasma structures was studied and the influence of the quasi-2D cluster size (a number of particles in it) on the features of the phase transition was investigated. Experiments and numerical simulation was conducted for the systems consisting of small (~ 10) and large (~ 103) number of particles. To investigate the phase state of the system with 7, 18 and 100 particles observed in numerical and laboratory experiments, we used the method based on analysis of dynamic entropy. Numerical modeling of small systems was conducted by the Langevin molecular dynamic method with the Langevin force, responsible for the stochastic nature of the motion of particles with a given kinetic temperature. Phase state of systems with the number of elements in the order of 103, was studied using the methods of statistical thermodynamics. Here we present new results of an experimental study of the change of translational and orientational order and topological defects, and the pair interactions at 2D melting of dust cluster in rf discharge plasma. The experimental results have revealed the existence of hexatic phase as well as solid-to-hexatic phase and hexatic-to-liquid transitions. This work was supported by the Russian Science Foundation (O.F. Petrov, M.M.Vasiliev, K.B. Stacenko, X.G. Koss, E.V. Vasilieva, M.I.Myasnikov and E.?.Lisin) through Grant No. 14-12-01440).
2-D/3-D ECE imaging data for validation of turbulence simulations
NASA Astrophysics Data System (ADS)
Choi, Minjun; Lee, Jaehyun; Yun, Gunsu; Lee, Woochang; Park, Hyeon K.; Park, Young-Seok; Sabbagh, Steve A.; Wang, Weixing; Luhmann, Neville C., Jr.
2015-11-01
The 2-D/3-D KSTAR ECEI diagnostic can provide a local 2-D/3-D measurement of ECE intensity. Application of spectral analysis techniques to the ECEI data allows local estimation of frequency spectra S (f) , wavenumber spectra S (k) , wavernumber and frequency spectra S (k , f) , and bispectra b (f1 ,f2) of ECE intensity over the 2-D/3-D space, which can be used to validate turbulence simulations. However, the minimum detectable fluctuation amplitude and the maximum detectable wavenumber are limited by the temporal and spatial resolutions of the diagnostic system, respectively. Also, the finite measurement area of the diagnostic channel could introduce uncertainty in the spectra estimation. The limitations and accuracy of the ECEI estimated spectra have been tested by a synthetic ECEI diagnostic with the model and/or fluctuations calculated by GTS. Supported by the NRF of Korea under Contract No. NRF-2014M1A7A1A03029881 and NRF-2014M1A7A1A03029865 and by U.S. DOE grant DE-FG02-99ER54524.
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.
Multipacting Simulation Study for 56 MHz Quarter Wave Resonator using 2D Code
Naik,D.; Ben-Zvi, I.
2009-01-02
A beam excited 56 MHz Radio Frequency (RF) Niobium Quarter Wave Resonator (QWR) has been proposed to enhance RHIC beam luminosity and bunching. Being a RF cavity, multipacting is expected; therefore an extensive study was carried out with the Multipac 2.1 2D simulation code. The study revealed that multipacting occurs in various bands up to peak surface electric field 50 kV/m and is concentrated mostly above the beam gap and on the outer conductor. To suppress multipacting, a ripple structure was introduced to the outer conductor and the phenomenon was successfully eliminated from the cavity.
NASA Astrophysics Data System (ADS)
Tzeferacos, Petros; Daley, Christopher; Fatenejad, Milad; Flocke, Norbert; Graziani, Carlo; Lamb, Donald Q.; Lee, Dongwook; Scopatz, Anthony; Weide, Klaus; Doyle, Hugo; Gregori, Gianluca; Meinecke, Jena; Reville, Brian; Miniati, Francesco
2013-10-01
The process of generation and amplification of Biermann battery magnetic fields is closely linked to the development of turbulence. In an astrophysical environment, a small seed field can be formed in asymmetric supernova remnant blast waves due to misaligned pressure and density gradients. Inhomogeneities in the density distribution can cause the flow to become turbulent and the B-field can be amplified via dynamo action. In this context, the COSMOLAB team will perform experiments using the Omega EP laser at LLE, that represent a scaled-down model of the astrophysical process in a controlled environment. The experiments involve the illumination of a slab-like target, which produces a plasma flow and a Biermann battery field. The flow then propagates through a grid that creates turbulence and amplifies the field. In this study we describe 2D and 3D radiative MHD simulations of the experimental setup, carried out using the FLASH code on Mira (BG/Q) at ALCF. The objective of these simulations is to explore the morphology and strength of the B-fields generated by ablation of target material by the laser, and their amplification due to the grid. This work was supported by DOE NNSA ASC.
Energy storage and dissipation in the magnetotail during substorms. 2. MHD simulations
Steinolfson, R.S. ); Winglee, R.M. )
1993-05-01
The authors present a global MHD simulation of the magnetotail in an effort to study magnetic storm development. They address the question of energy storage in the current sheet in the early phases of storm growth, which previous simulations have not shown. They address this problem by dealing with the variation of the resistivity throughout the magnetosphere. They argue that MHD theory should provide a suitable representation to this problem on a global scale, even if it does not handle all details adequately. For their simulation they use three different forms for the resistivity. First is a uniform and constant resistivity. Second is a resistivity proportional to the current density, which is related to argument that resistivity is driven by wave-particle interactions which should be strongest in regions where the current is the greatest. Thirdly is a model where the resistivity varies with the magnetic field strength, which was suggested by previous results from particle simulations of the same problem. The simulation then gives approximately the same response of the magnetosphere for all three of the models. Each results in the formation and ejection of plasmoids, but the energy stored in the magnetotail, the timing of substorm onset in relation to the appearance of a southward interplanetary magnetic field, and the speed of ejection of the plasmoids formed differ with the resistivity models.
NASA Technical Reports Server (NTRS)
Klimas, A. J.; Uritsky, V.; Vassiliadis, D.; Baker, D. N.
2005-01-01
Loading and consequent unloading of magnetic flux is an essential element of the substorm cycle in Earth's magnetotail. We are unaware of an available global MHD magnetospheric simulation model that includes a loading- unloading cycle in its behavior. Given the central role that MHD models presently play in the development of our understanding of magnetospheric dynamics, and given the present plans for the central role that these models will play in ongoing space weather prediction programs, it is clear that this failure must be corrected. A 2-dimensional numerical driven current-sheet model has been developed that incorporates an idealized current- driven instability with a resistive MHD system. Under steady loading, the model exhibits a global loading- unloading cycle. The specific mechanism for producing the loading-unloading cycle will be discussed. It will be shown that scale-free avalanching of electromagnetic energy through the model, from loading to unloading, is carried by repetitive bursts of localized reconnection. Each burst leads, somewhat later, to a field configuration that is capable of exciting a reconnection burst again. This process repeats itself in an intermittent manner while the total field energy in the system falls. At the end of an unloading interval, the total field energy is reduced to well below that necessary to initiate the next unloading event and, thus, a loading-unloading cycle results. It will be shown that, in this model, it is the topology of bursty localized reconnection that is responsible for the appearance of the loading-unloading cycle.
Simulation of 3-D Nonequilibrium Seeded Air Flow in the NASA-Ames MHD Channel
NASA Technical Reports Server (NTRS)
Gupta, Sumeet; Tannehill, John C.; Mehta, Unmeel B.
2004-01-01
The 3-D nonequilibrium seeded air flow in the NASA-Ames experimental MHD channel has been numerically simulated. The channel contains a nozzle section, a center section, and an accelerator section where magnetic and electric fields can be imposed on the flow. In recent tests, velocity increases of up to 40% have been achieved in the accelerator section. The flow in the channel is numerically computed us ing a 3-D parabolized Navier-Stokes (PNS) algorithm that has been developed to efficiently compute MHD flows in the low magnetic Reynolds number regime: The MHD effects are modeled by introducing source terms into the PNS equations which can then be solved in a very efficient manner. The algorithm has been extended in the present study to account for nonequilibrium seeded air flows. The electrical conductivity of the flow is determined using the program of Park. The new algorithm has been used to compute two test cases that match the experimental conditions. In both cases, magnetic and electric fields are applied to the seeded flow. The computed results are in good agreement with the experimental data.
Analysis and statistics of discontinuities as obtained from 3D simulation of MHD turbulence
NASA Astrophysics Data System (ADS)
Zhang, Lei; He, Jian-Sen; Tu, Chuan-Yi; Yang, Li-Ping; Wang, Xin; Marsch, Eckart; Wang, Ling-Hua
2016-03-01
The turbulent solar wind abounds with MHD discontinuities, and such discontinuities are often found in close connection with turbulence intermittency, constituting a possible main contributor to the turbulence dissipation and solar wind heating. Among the discontinuities, tangential (TD) and rotational (RD) ones are two most important types. Recently, the connection between turbulence intermittency and proton thermodynamics has been being intensively investigated. Such connections are founded to be involved with MHD instablilities, but the difference of TDs an RDs in this process has not yet been covered. Herewith we define new methods for identifying TDs and RDs obtained from a three-dimensional MHD simulation with pressure anisotropy. Especially, we define the Total Variance of Increments (TVI) as a new measure of magnetic field changes. Based on the identified cases, we compare their occurrence rates and heating effects. More specifically, we find that the thermal states embedding TDs, compared with their RD counterparts, tend to be more associated with extreme plasma parameters or instabilites. Some other possible applications of TVI-like norms are also herewith discussed.
Global MHD simulations of cosmic ray driven galactic winds
NASA Astrophysics Data System (ADS)
Ruszkowski, Mateusz; Yang, Hsiang-Yi Karen; Gould Zweibel, Ellen
2016-04-01
Galactic outflows play an important role in galactic evolution. Despite their importance, a detailed understanding of the physical mechanisms responsible for the driving of these winds is lacking. In an effort to gain more insight into the nature of these flows, we perform global three-dimensional magneto-hydrodynamical simulations of an isolated starbursting galaxy. We focus on the dynamical role of cosmic rays injected by supernovae, and specifically on the impact of the streaming and anisotropic diffusion of cosmic rays along the magnetic fields. We find that these microphysical effects can have a significant effect on the wind launching and mass loading factors depending on the details of the plasma physics. Cosmic rays stream away from the densest regions near the galactic disk along partially ordered magnetic fields and, in the process, accelerate more tenuous gas away from the galaxy. For cosmic ray acceleration efficiencies broadly consistent with the observational constraints, cosmic rays are likely to have a notable impact on the wind launching.
Coronal extension of the MURaM radiative MHD code: From quiet sun to flare simulations
NASA Astrophysics Data System (ADS)
Rempel, Matthias D.; Cheung, Mark
2016-05-01
We present a new version of the MURaM radiative MHD code, which includes a treatment of the solar corona in terms of MHD, optically thin radiative loss and field-aligned heat conduction. In order to relax the severe time-step constraints imposed by large Alfven velocities and heat conduction we use a combination of semi-relativistic MHD with reduced speed of light ("Boris correction") and a hyperbolic formulation of heat conduction. We apply the numerical setup to 4 different setups including a mixed polarity quiet sun, an open flux region, an arcade solution and an active region setup and find all cases an amount of coronal heating sufficient to maintain a corona with temperatures from 1 MK (quiet sun) to 2 MK (active region, arcade). In all our setups the Poynting flux is self-consistently created by photospheric and sub-photospheric magneto-convection in the lower part of our simulation domain. Varying the maximum allowed Alfven velocity ("reduced speed of light") leads to only minor changes in the coronal structure as long as the limited Alfven velocity remains larger than the speed of sound and about 1.5-3 times larger than the peak advection velocity. We also found that varying details of the numerical diffusivities that govern the resistive and viscous energy dissipation do not strongly affect the overall coronal heating, but the ratio of resistive and viscous energy dependence is strongly dependent on the effective numerical magnetic Prandtl number. We use our active region setup in order to simulate a flare triggered by the emergence of a twisted flux rope into a pre-existing bipolar active region. Our simulation yields a series of flares, with the strongest one reaching GOES M1 class. The simulation reproduces many observed properties of eruptions such as flare ribbons, post flare loops and a sunquake.
Double Dynamo Signatures in a Global MHD Simulation and Mean-field Dynamos
NASA Astrophysics Data System (ADS)
Beaudoin, Patrice; Simard, Corinne; Cossette, Jean-François; Charbonneau, Paul
2016-08-01
The 11 year solar activity cycle is the most prominent periodic manifestation of the magnetohydrodynamical (MHD) large-scale dynamo operating in the solar interior, yet longer and shorter (quasi-) periodicities are also present. The so-called “quasi-biennial” signal appearing in many proxies of solar activity has been gaining increasing attention since its detection in p-mode frequency shifts, which suggests a subphotospheric origin. A number of candidate mechanisms have been proposed, including beating between co-existing global dynamo modes, dual dynamos operating in spatially separated regions of the solar interior, and Rossby waves driving short-period oscillations in the large-scale solar magnetic field produced by the 11 year activity cycle. In this article, we analyze a global MHD simulation of solar convection producing regular large-scale magnetic cycles, and detect and characterize shorter periodicities developing therein. By constructing kinematic mean-field α 2Ω dynamo models incorporating the turbulent electromotive force (emf) extracted from that same simulation, we find that dual-dynamo behavior materializes in fairly wide regions of the model’s parameters space. This suggests that the origin of the similar behavior detected in the MHD simulation lies with the joint complexity of the turbulent emf and differential rotation profile, rather that with dynamical interactions such as those mediated by Rossby waves. Analysis of the simulation also reveals that the dual dynamo operating therein leaves a double-period signature in the temperature field, consistent with a dual-period helioseismic signature. Order-of-magnitude estimates for the magnitude of the expected frequency shifts are commensurate with helioseismic measurements. Taken together, our results support the hypothesis that the solar quasi-biennial oscillations are associated with a secondary dynamo process operating in the outer reaches of the solar convection zone.
NASA Astrophysics Data System (ADS)
Borrero, J. M.; Lites, B. W.; Lagg, A.; Rezaei, R.; Rempel, M.
2014-12-01
Milne-Eddington (M-E) inversion codes for the radiative transfer equation are the most widely used tools to infer the magnetic field from observations of the polarization signals in photospheric and chromospheric spectral lines. Unfortunately, a comprehensive comparison between the different M-E codes available to the solar physics community is still missing, and so is a physical interpretation of their inferences. In this contribution we offer a comparison between three of those codes (VFISV, ASP/HAO, and HeLIx+). These codes are used to invert synthetic Stokes profiles that were previously obtained from realistic non-grey three-dimensional magnetohydrodynamical (3D MHD) simulations. The results of the inversion are compared with each other and with those from the MHD simulations. In the first case, the M-E codes retrieve values for the magnetic field strength, inclination and line-of-sight velocity that agree with each other within σB ≤ 35 (Gauss), σγ ≤ 1.2°, and σv ≤ 10 m s-1, respectively. Additionally, M-E inversion codes agree with the numerical simulations, when compared at a fixed optical depth, within σB ≤ 130 (Gauss), σγ ≤ 5°, and σv ≤ 320 m s-1. Finally, we show that employing generalized response functions to determine the height at which M-E codes measure physical parameters is more meaningful than comparing at a fixed geometrical height or optical depth. In this case the differences between M-E inferences and the 3D MHD simulations decrease to σB ≤ 90 (Gauss), σγ ≤ 3°, and σv ≤ 90 m s-1.
The Substorm Current Wedge: Further Insights from MHD Simulations
NASA Technical Reports Server (NTRS)
Birn, J.; Hesse, M.
2015-01-01
Using a recent magnetohydrodynamic simulation of magnetotail dynamics, we further investigate the buildup and evolution of the substorm current wedge (SCW), resulting from flow bursts generated by near-tail reconnection. Each flow burst generates an individual current wedge, which includes the reduction of cross-tail current and the diversion to region 1 (R1)-type field-aligned currents (earthward on the dawn and tailward on the duskside), connecting the tail with the ionosphere. Multiple flow bursts generate initially multiple SCW patterns, which at later times combine to a wider single SCW pattern. The standard SCWmodel is modified by the addition of several current loops, related to particular magnetic field changes: the increase of Bz in a local equatorial region (dipolarization), the decrease of |Bx| away from the equator (current disruption), and increases in |By| resulting from azimuthally deflected flows. The associated loop currents are found to be of similar magnitude, 0.1-0.3 MA. The combined effect requires the addition of region 2 (R2)-type currents closing in the near tail through dawnward currents but also connecting radially with the R1 currents. The current closure at the inner boundary, taken as a crude proxy of an idealized ionosphere, demonstrates westward currents as postulated in the original SCW picture as well as North-South currents connecting R1- and R2-type currents, which were larger than the westward currents by a factor of almost 2. However, this result should be applied with caution to the ionosphere because of our neglect of finite resistance and Hall effects.
NASA Astrophysics Data System (ADS)
Krause, M.; Camenzind, M.
2001-12-01
In the present paper, we examine the convergence behavior and inter-code reliability of astrophysical jet simulations in axial symmetry. We consider both pure hydrodynamic jets and jets with a dynamically significant magnetic field. The setups were chosen to match the setups of two other publications, and recomputed with the MHD code NIRVANA. We show that NIRVANA and the two other codes give comparable, but not identical results. We explain the differences by the different application of artificial viscosity in the three codes and numerical details, which can be summarized in a resolution effect, in the case without magnetic field: NIRVANA turns out to be a fair code of medium efficiency. It needs approximately twice the resolution as the code by Lind (Lind et al. 1989) and half the resolution as the code by Kössl (Kössl & Müller 1988). We find that some global properties of a hydrodynamical jet simulation, like e.g. the bow shock velocity, converge at 100 points per beam radius (ppb) with NIRVANA. The situation is quite different after switching on the toroidal magnetic field: in this case, global properties converge even at 10 ppb. In both cases, details of the inner jet structure and especially the terminal shock region are still insufficiently resolved, even at our highest resolution of 70 ppb in the magnetized case and 400 ppb for the pure hydrodynamic jet. The magnetized jet even suffers from a fatal retreat of the Mach disk towards the inflow boundary, which indicates that this simulation does not converge, in the end. This is also in definite disagreement with earlier simulations, and challenges further studies of the problem with other codes. In the case of our highest resolution simulation, we can report two new features: first, small scale Kelvin-Helmholtz instabilities are excited at the contact discontinuity next to the jet head. This slows down the development of the long wavelength Kelvin-Helmholtz instability and its turbulent cascade to smaller
Ion acoustic wave collapse via two-ion wave decay: 2D Vlasov simulation and theory
NASA Astrophysics Data System (ADS)
Chapman, Thomas; Berger, Richard; Banks, Jeffrey; Brunner, Stephan
2015-11-01
The decay of ion acoustic waves (IAWs) via two-ion wave decay may transfer energy from the electric field of the IAWs to the particles, resulting in a significant heating of resonant particles. This process has previously been shown in numerical simulations to decrease the plasma reflectivity due to stimulated Brillouin scattering. Two-ion wave decay is a fundamental property of ion acoustic waves that occurs over most if not all of the parameter space of relevance to inertial confinement fusion experiments, and can lead to a sudden collapse of IAWs. The treatment of all species kinetically, and in particular the electrons, is required to describe the decay process correctly. We present fully kinetic 2D+2V Vlasov simulations of IAWs undergoing decay to a highly nonlinear turbulent state using the code LOKI. The scaling of the decay rate with characteristic plasma parameters and wave amplitude is shown. A new theory describing two-ion wave decay in 2D, that incorporates key kinetic properties of the electrons, is presented and used to explain quantitatively for the first time the observed decay of IAWs. Work performed under auspices of U.S. DoE by LLNL, Contract DE-AC52-07NA2734. Funded by LDRD 15-ERD-038 and supported by LLNL Grand Challenge allocation.
Superclusters of galaxies from the 2dF redshift survey. 2. Comparison with simulations
Einasto, Jaan; Einasto, M.; Saar, E.; Tago, E.; Liivamagi, L.J.; Joeveer, M.J; Suhhonenko, I.; Hutsi, G.; Jaaniste, J.; Heinamaki, P.; Muller, V.; Knebe, A.; Tucker, D.; /Fermilab
2006-04-01
We investigate properties of superclusters of galaxies found on the basis of the 2dF Galaxy Redshift Survey, and compare them with properties of superclusters from the Millennium Simulation.We study the dependence of various characteristics of superclusters on their distance from the observer, on their total luminosity, and on their multiplicity. The multiplicity is defined by the number of Density Field (DF) clusters in superclusters. Using the multiplicity we divide superclusters into four richness classes: poor, medium, rich and extremely rich.We show that superclusters are asymmetrical and have multi-branching filamentary structure, with the degree of asymmetry and filamentarity being higher for the more luminous and richer superclusters. The comparison of real superclusters with Millennium superclusters shows that most properties of simulated superclusters agree very well with real data, the main differences being in the luminosity and multiplicity distributions.
Calibration and simulation of ASM2d at different temperatures in a phosphorus removal pilot plant.
García-Usach, F; Ferrer, J; Bouzas, A; Seco, A
2006-01-01
In this work, an organic and nutrient removal pilot plant was used to study the temperature influence on phosphorus accumulating organisms. Three experiments were carried out at 13, 20 and 24.5 degrees C, achieving a high phosphorus removal percentage in all cases. The ASM2d model was calibrated at 13 and 20 degrees C and the Arrhenius equation constant was obtained for phosphorus removal processes showing that the temperature influences on the biological phosphorus removal subprocesses in a different degree. The 24.5 degrees C experiment was simulated using the model parameters obtained by means of the Arrhenius equation. The simulation results for the three experiments showed good correspondence with the experimental data, demonstrating that the model and the calibrated parameters were able to predict the pilot plant behaviour. PMID:16889256
Well-posedness and generalized plane waves simulations of a 2D mode conversion model
NASA Astrophysics Data System (ADS)
Imbert-Gérard, Lise-Marie
2015-12-01
Certain types of electro-magnetic waves propagating in a plasma can undergo a mode conversion process. In magnetic confinement fusion, this phenomenon is very useful to heat the plasma, since it permits to transfer the heat at or near the plasma center. This work focuses on a mathematical model of wave propagation around the mode conversion region, from both theoretical and numerical points of view. It aims at developing, for a well-posed equation, specific basis functions to study a wave mode conversion process. These basis functions, called generalized plane waves, are intrinsically based on variable coefficients. As such, they are particularly adapted to the mode conversion problem. The design of generalized plane waves for the proposed model is described in detail. Their implementation within a discontinuous Galerkin method then provides numerical simulations of the process. These first 2D simulations for this model agree with qualitative aspects studied in previous works.
NASA Astrophysics Data System (ADS)
Bezzeccheri, E.; Colasanti, S.; Falco, A.; Liguori, R.; Rubino, A.; Lugli, P.
2016-05-01
Vertical Organic Transistors and Phototransistors have been proven to be promising technologies due to the advantages of reduced channel length and larger sensitive area with respect to planar devices. Nevertheless, a real improvement of their performance is subordinate to the quantitative description of their operation mechanisms. In this work, we present a comparative study on the modeling of vertical and planar Organic Phototransistor (OPT) structures. Computer-based simulations of the devices have been carried out with Synopsys Sentaurus TCAD in a 2D Drift-Diffusion framework. The photoactive semiconductor material has been modeled using the virtual semiconductor approach as the archetypal P3HT:PC61BM bulk heterojunction. It has been found that both simulated devices have comparable electrical and optical characteristics, accordingly to recent experimental reports on the subject.
Extended MHD Simulations of Tearing Instabilities and the Dynamo Effect in the Reversed-Field Pinch
NASA Astrophysics Data System (ADS)
Germaschewski, K.; Dearborn, J.; Bhattacharjee, A.
2009-11-01
Observations on MST indicate the importance of the Hall current in sawtooth crashes and the dynamo effect in a RFP. We employ our Magnetic Reconnection Code (MRC) to perform fully 3D extended MHD simulations in the RFP, including the Hall current and electron pressure gradient in a generalized Ohm's law. The MRC is an MPI-parallelized finite-volume based simulation code that integrates the extended MHD equations. It supports arbitrary curvilinear coordinate mappings, allowing it to be adapted to cylindrical and toroidal geometries. In order to overcome restrictive time-step limits, it uses implicit time integration. We have benchmarked the code for linear tearing instabilities, and performed fully nonlinear simulations. Due to the presence of the Hall current, novel vortical flows are seen in the vicinity of rational surfaces, akin to those seen in recent sawtooth studies in tokamaks, when the peak of the current density separates from the stagnation point of the flow. We calculate the dynamo field by averaging, and compare simulations with observations.
Highly-resolved 2D HYDRA simulations of Double-Shell Ignition Designs
Milovich, J L; Amendt, P; Hamza, A; Marinak, M; Robey, H
2006-06-30
Double-shell (DS) targets (Amendt, P. A. et al., 2002) offer a complementary approach to the cryogenic baseline design (Lindl, J. et al., 2004) for achieving ignition on the National Ignition Facility (NIF). Among the expected benefits are the ease of room temperature preparation and fielding, the potential for lower laser backscatter and the reduced need for careful shock timing. These benefits are offset, however, by demanding fabrication tolerances, e.g., shell concentricity and shell surface smoothness. In particular, the latter is of paramount importance since DS targets are susceptible to the growth of interface perturbations from impulsive and time-dependent accelerations. Previous work (Milovich, J. L. et al., 2004) has indicated that the growth of perturbations on the outer surface of the inner shell is potentially disruptive. To control this instability new designs have been proposed requiring bimetallic inner shells and material-matching mid-Z nanoporous foam. The challenges in manufacturing such exotic foams have led to a further evaluation of the densities and pore sizes needed to reduce the seeding of perturbations on the outer surface of the inner shell, thereby guiding the ongoing material science research efforts. Highly-resolved 2D simulations of porous foams have been performed to establish an upper limit on the allowable pore sizes for instability growth. Simulations indicate that foams with higher densities than previously thought are now possible. Moreover, while at the present time we are only able to simulate foams with average pore sizes larger than 1 micron (due to computational limitations), we can conclude that these pore sizes are potentially problematic. Furthermore, the effect of low-order hohlraum radiation asymmetries on the growth of intrinsic surface perturbations is also addressed. Highly-resolved 2D simulations indicate that the transverse flows that are set up by these low-order mode features (which can excite Kelvin
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
Propagation of Pi2 pulsations through the braking region in global MHD simulations
NASA Astrophysics Data System (ADS)
Ream, J. B.; Walker, R. J.; Ashour-Abdalla, M.; El-Alaoui, M.; Wiltberger, M.; Kivelson, M. G.; Goldstein, M. L.
2015-12-01
We investigate the propagation of Pi2 period pulsations from their origin in the plasma sheet through the braking region, the region where the fast flows are slowed as they approach the inner edge of the plasma sheet. Our approach is to use both the University of California, Los Angeles (UCLA) and Lyon-Fedder-Mobarry (LFM) global magnetohydrodynamic (MHD) computer codes to simulate the Earth's magnetosphere during a substorm that occurred on 14 September 2004 when Pi2 pulsations were observed. We use two different MHD models in order to test the robustness of our conclusions about Pi2. The simulations are then compared with ground-based and satellite data. We find that the propagation of the pulsations in the simulations, especially through the braking region, depends strongly on the ionospheric models used at the inner boundary of the MHD models. With respect to typical observed values, the modeled conductances are high in the UCLA model and low in the LFM model. The different conductances affect the flows, producing stronger line tying that slows the flow in the braking region more in the UCLA model than in the LFM model. Therefore, perturbations are able to propagate much more freely into the inner magnetosphere in the LFM results. However, in both models Pi2 period perturbations travel with the dipolarization front (DF) that forms at the earthward edge of the flow channel, but as the DF slows in the braking region, -8≤x≤-6 RE, the Pi2 period perturbations begin to travel ahead of it into the inner magnetosphere. This indicates that the flow channels generate compressional waves with periods that fall within the Pi2 range and that, as the flows themselves are stopped in the braking region, the compressional wave continues to propagate into the inner magnetosphere.
Reynolds-Averaged Navier-Stokes Simulation of a 2D Circulation Control Wind Tunnel Experiment
NASA Technical Reports Server (NTRS)
Allan, Brian G.; Jones, Greg; Lin, John C.
2011-01-01
Numerical simulations are performed using a Reynolds-averaged Navier-Stokes (RANS) flow solver for a circulation control airfoil. 2D and 3D simulation results are compared to a circulation control wind tunnel test conducted at the NASA Langley Basic Aerodynamics Research Tunnel (BART). The RANS simulations are compared to a low blowing case with a jet momentum coefficient, C(sub u), of 0:047 and a higher blowing case of 0.115. Three dimensional simulations of the model and tunnel walls show wall effects on the lift and airfoil surface pressures. These wall effects include a 4% decrease of the midspan sectional lift for the C(sub u) 0.115 blowing condition. Simulations comparing the performance of the Spalart Allmaras (SA) and Shear Stress Transport (SST) turbulence models are also made, showing the SST model compares best to the experimental data. A Rotational/Curvature Correction (RCC) to the turbulence model is also evaluated demonstrating an improvement in the CFD predictions.
Spot size variation FCS in simulations of the 2D Ising model
NASA Astrophysics Data System (ADS)
Burns, Margaret C.; Nouri, Mariam; Veatch, Sarah L.
2016-06-01
Spot variation fluorescence correlation spectroscopy (svFCS) was developed to study the movement and organization of single molecules in plasma membranes. This experimental technique varies the size of an illumination area while measuring correlations in time using standard fluorescence correlation methods. Frequently, this data is interpreted using the assumption that correlation measurements reflect the dynamics of single molecule motions, and not motions of the average composition. Here, we explore how svFCS measurements report on the dynamics of components diffusing within simulations of a 2D Ising model with a conserved order parameter. Simulated correlation functions report on both the fast dynamics of single component mobility and the slower dynamics of the average composition. Over a range of simulation conditions, a conventional svFCS analysis suggests the presence of anomalous diffusion even though single molecule motions are nearly Brownian in these simulations. This misinterpretation is most significant when the surface density of the fluorescent label is elevated, therefore we suggest future measurements be made over a range of tracer densities. Some simulation conditions reproduce qualitative features of published svFCS experimental data. Overall, this work emphasizes the need to probe membranes using multiple complimentary experimental methodologies in order to draw conclusions regarding the nature of spatial and dynamical heterogeneity in these systems.
Isolated Coronal Mass Ejections and Associated Phenomena: MHD Simulations and STEREO Observations
NASA Astrophysics Data System (ADS)
Lugaz, Noé; Roussev, Ilia; Sokolov, Igor; Shibata, Kazunari; Schwadron, Nathan; Downs, Cooper
2012-07-01
Coronal Mass Ejections (CMEs), the most energetic events occurring in our solar system, are associated with a range of other phenomena such as EUV waves, dimming regions and solar energetic particles (SEPs). With the advancement of global numerical simulations and the new fleet of spacecraft observing the Sun and the heliosphere, it is possible to combine simulations with data analyses to gain new insight into the complex Sun-Earth system. In this talk, I will discuss magneto-hydrodynamic (MHD) simulations of CMEs and associated phenomena. In particular, I will focus on the changes in magnetic connectivity and the opening of previously closed field lines during and after an eruption and discuss the consequences for the acceleration and transport of energetic particles. I will also present combined numerical simulations and observations by SDO and STEREO/SECCHI of CMEs and EUV waves, which have led to new insights into CME deflection, expansion and rotation.
Relativistic modeling capabilities in PERSEUS extended MHD simulation code for HED plasmas
Hamlin, Nathaniel D.; Seyler, Charles E.
2014-12-15
We discuss the incorporation of relativistic modeling capabilities into the PERSEUS extended MHD simulation code for high-energy-density (HED) plasmas, and present the latest hybrid X-pinch simulation results. The use of fully relativistic equations enables the model to remain self-consistent in simulations of such relativistic phenomena as X-pinches and laser-plasma interactions. By suitable formulation of the relativistic generalized Ohm’s law as an evolution equation, we have reduced the recovery of primitive variables, a major technical challenge in relativistic codes, to a straightforward algebraic computation. Our code recovers expected results in the non-relativistic limit, and reveals new physics in the modeling of electron beam acceleration following an X-pinch. Through the use of a relaxation scheme, relativistic PERSEUS is able to handle nine orders of magnitude in density variation, making it the first fluid code, to our knowledge, that can simulate relativistic HED plasmas.
NASA Astrophysics Data System (ADS)
Gao, Donghong
Interest in utilizing liquid metal film flows to protect the plasma-facing solid structures places increasing demand on understanding the magnetohydrodynamics (MHD) of such flows in a magnetic field with spatial variation. The field gradient effect is studied by a two-dimensional (2D) model in Cartesian coordinates. The thin film flow down an inclined plane in spanwise (z-direction) magnetic field with constant streamwise gradient and applied current is analyzed. The solution to the equilibrium flow shows forcefully the M-shaped velocity profile and dependence of side layer thickness on Ha-1/2 whose definition is based on field gradient. The major part of the dissertation is the numerical simulation of free surface film flows and understanding the results. The VOF method is employed to track the free surface, and the CSF model is combined with VOF method to account for surface dynamics condition. The code is validated with respect to Navier-Stokes solver and MHD implementation by computations of ordinary wavy films, MHD flat films and a colleague proposed film flow. The comparisons are performed against respective experimental, theoretical or numerical solutions, and the results are well matched with them. It is found for the ordinary water falling films, at low frequency and high flowrate, the small forcing disturbance at inlet flowrate develops into big roll waves preceded by small capillary bow waves; at high frequency and low Re, it develops into nearly sinusoidal waves with small amplitude and without fore-running capillary waves. The MHD surface instability is investigated for two kinds of film flows in constant streamwise field gradient: one with spatial disturbance and without surface tension, the other with inlet forcing disturbance and with surface tension. At no surface tension condition, the finite amplitude disturbance is rapidly amplified and degrades to irregular shape. With surface tension to maintain smooth interface, finite amplitude regular waves
NASA Astrophysics Data System (ADS)
Tarditi, Alfonso G.; Shebalin, John V.
2002-11-01
A simulation study with the NIMROD code [1] is being carried on to investigate the efficiency of the thrust generation process and the properties of the plasma detachment in a magnetic nozzle. In the simulation, hot plasma is injected in the magnetic nozzle, modeled as a 2D, axi-symmetric domain. NIMROD has two-fluid, 3D capabilities but the present runs are being conducted within the MHD, 2D approximation. As the plasma travels through the magnetic field, part of its thermal energy is converted into longitudinal kinetic energy, along the axis of the nozzle. The plasma eventually detaches from the magnetic field at a certain distance from the nozzle throat where the kinetic energy becomes larger than the magnetic energy. Preliminary NIMROD 2D runs have been benchmarked with a particle trajectory code showing satisfactory results [2]. Further testing is here reported with the emphasis on the analysis of the diffusion rate across the field lines and of the overall nozzle efficiency. These simulation runs are specifically designed for obtaining comparisons with laboratory measurements of the VASIMR experiment, by looking at the evolution of the radial plasma density and temperature profiles in the nozzle. VASIMR (Variable Specific Impulse Magnetoplasma Rocket, [3]) is an advanced space propulsion concept currently under experimental development at the Advanced Space Propulsion Laboratory, NASA Johnson Space Center. A plasma (typically ionized Hydrogen or Helium) is generated by a RF (Helicon) discharge and heated by an Ion Cyclotron Resonance Heating antenna. The heated plasma is then guided into a magnetic nozzle to convert the thermal plasma energy into effective thrust. The VASIMR system has no electrodes and a solenoidal magnetic field produced by an asymmetric mirror configuration ensures magnetic insulation of the plasma from the material surfaces. By powering the plasma source and the heating antenna at different levels it is possible to vary smoothly of the
Rise characteristics of gas bubbles in a 2D rectangular column: VOF simulations vs experiments
Krishna, R.; Baten, J.M. van
1999-10-01
About five centuries ago, Leonardo da Vinci described the sinuous motion of gas bubbles rising in water. The authors have attempted to simulate the rise trajectories of bubbles of 4, 5, 7, 8, 9, 12, and 20 mm in diameter rising in a 2D rectangular column filled with water. The simulations were carried out using the volume-of-fluid (VOF) technique developed by Hirt and Nichols (J. Computational Physics, 39, 201--225 (1981)). To solve the Navier-Stokes equations of motion the authors used a commercial solver, CFX 4.1c of AEA Technology, UK. They developed their own bubble-tracking algorithm to capture sinuous bubble motions. The 4 and 5 mm bubbles show large lateral motions observed by Da Vinci. The 7, 8 and 9 mm bubble behave like jellyfish. The 12 mm bubble flaps its wings like a bird. The extent of lateral motion of the bubbles decreases with increasing bubble size. Bubbles larger than 20 mm in size assume a spherical cap form and simulations of the rise characteristics match experiments exactly. VOF simulations are powerful tools for a priori determination of the morphology and rise characteristics of bubbles rising in a liquid. Bubble-bubble interactions are also properly modeled by the VOF technique.
Numerical simulation of 2D buoyant jets in ice-covered and temperature-stratified water
NASA Astrophysics Data System (ADS)
Gu, Ruochuan
A two-dimensional (2D) unsteady simulation model is applied to the problem of a submerged warm water discharge into a stratified lake or reservoir with an ice cover. Numerical simulations and analyses are conducted to gain insight into large-scale convective recirculation and flow processes in a cold waterbody induced by a buoyant jet. Jet behaviors under various discharge temperatures are captured by directly modeling flow and thermal fields. Flow structures and processes are described by the simulated spatial and temporal distributions of velocity and temperature in various regions: deflection, recirculation, attachment, and impingement. Some peculiar hydrothermal and dynamic features, e.g. reversal of buoyancy due to the dilution of a warm jet by entraining cold ambient water, are identified and examined. Simulation results show that buoyancy is the most important factor controlling jet behavior and mixing processes. The inflow boundary is treated as a liquid wall from which the jet is offset. Similarity and difference in effects of boundaries perpendicular and parallel to flow, and of buoyancy on jet attachment and impingement, are discussed. Symmetric flow configuration is used to de-emphasize the Coanda effect caused by offset.
Incorporating a Turbulence Transport Model into 2-D Hybrid Hall Thruster Simulations
NASA Astrophysics Data System (ADS)
Cha, Eunsun; Cappelli, Mark A.; Fernandez, Eduardo
2014-10-01
2-D hybrid simulations of Hall plasma thrusters that do not resolve cross-field transport-generating fluctuations require a model to capture how electrons migrate across the magnetic field. We describe the results of integrating a turbulent electron transport model into simulations of plasma behavior in a plane spanned by the E and B field vectors. The simulations treat the electrons as a fluid and the heavy species (ions/neutrals) as discrete particles. The transport model assumes that the turbulent eddy cascade in the electron fluid to smaller scales is the primary means of electron energy dissipation. Using this model, we compare simulations to experimental measurements made on a laboratory Hall discharge over a range of discharge voltage. Both the current-voltage trends as well as the plasma properties such as plasma temperature, electron density, and ion velocities seem agree favorably with experiments, where a simple Bohm transport model tends to perform poorly in capturing much of the discharge behavior.
NASA Technical Reports Server (NTRS)
Ding, D. Q.; Denton, . E.; Hudson, M. K.; Lysak, R. L.
1995-01-01
The poloidal mode field line resonance in the Earth's dipole magnetic field is investigated using cold plasma ideal MHD simulations in dipole geometry. In order to excite the poloidal mode resonance, we use either an initial or a continuous velocity perturbation to drive the system. The perturbation is localized at magnetic shell L = 7 with plasma flow in the radial direction (electric field component in the azimuthal direction). It is found that with the initial perturbation alone, no polodial mode resonance can be obtained and the initially localized perturbation spreads out across all magnetic L shells. With the continuous perturbation, oscillating near the poloidal resonance frequency, a global-scale poloidal cavity mode can be obtained. For the first time, a localized guided poloidal mode resonance is obtained when a radial component of electric field is added to the initial perturbation such that the curl of the electric field is everywhere perpendicular to the background dipole magnetic field. During the localized poloidal resonance, plasma vortices parallel/antiparallel to the background dipole magnetic field B(sub 0). This circular flow, elongated radially, results in twisting of magnetic field flux tubes, which, in turn, leads to the slowdown of the circular plasma flow and reversal of the plasma vortices. The energy associated with the localized poloidal resonance is conserved as it shifts back and forth between the oscillating plasma vortices and the alternately twisted magnetic flux tubes. In the simulations the eigenfunctions associated with the localized poloidal resonance are grid-scale singular functions. This result indicates that ideal MHD is inadequate to describe the underlying problem and nonideal MHD effects are needed for mode broadening.
Jupiter Magnetotail Interaction with a Variable Solar Wind: A 3D MHD Simulation
NASA Astrophysics Data System (ADS)
Ranquist, D. A.; Bagenal, F.; Delamere, P. A.; Ma, X.
2015-12-01
Jupiter's magnetosphere is the largest object within the heliosphere. Voyager 2 detected its influence at Saturn's orbit, 4.3 AU away. It takes considerable time, therefore, for the solar wind to propagate such lengths down the tail. This propagation time is much greater than typical periods between changes in direction of the interplanetary magnetic field (IMF). We expect these variable magnetic fields to create a jumbled structure in Jupiter's magnetotail, resulting in magnetic reconnection and other magnetic processes. We simulate the global interaction of the solar wind with Jupiter's magnetosphere using a 3D magnetohydrodynamics (MHD) code. Delamere & Bagenal (2010) argue that the interaction is largely viscous, so we simulate the jovian magnetosphere as a region where the momentum equation has an added loss term. We also use in situ data gathered by the Ulysses spacecraft near Jupiter's orbit for solar wind input. Here, we report on the simulated dynamics in Jupiter's tail region.
NASA Astrophysics Data System (ADS)
Li, Jinghe; Song, Linping; Liu, Qing Huo
2016-02-01
A simultaneous multiple frequency contrast source inversion (CSI) method is applied to reconstructing hydrocarbon reservoir targets in a complex multilayered medium in two dimensions. It simulates the effects of a salt dome sedimentary formation in the context of reservoir monitoring. In this method, the stabilized biconjugate-gradient fast Fourier transform (BCGS-FFT) algorithm is applied as a fast solver for the 2D volume integral equation for the forward computation. The inversion technique with CSI combines the efficient FFT algorithm to speed up the matrix-vector multiplication and the stable convergence of the simultaneous multiple frequency CSI in the iteration process. As a result, this method is capable of making quantitative conductivity image reconstruction effectively for large-scale electromagnetic oil exploration problems, including the vertical electromagnetic profiling (VEP) survey investigated here. A number of numerical examples have been demonstrated to validate the effectiveness and capacity of the simultaneous multiple frequency CSI method for a limited array view in VEP.
Relaxation of ferroelectric states in 2D distributions of quantum dots: EELS simulation
NASA Astrophysics Data System (ADS)
Cortés, C. M.; Meza-Montes, L.; Moctezuma, R. E.; Carrillo, J. L.
2016-06-01
The relaxation time of collective electronic states in a 2D distribution of quantum dots is investigated theoretically by simulating EELS experiments. From the numerical calculation of the probability of energy loss of an electron beam, traveling parallel to the distribution, it is possible to estimate the damping time of ferroelectric-like states. We generate this collective response of the distribution by introducing a mean field interaction among the quantum dots, and then, the model is extended incorporating effects of long-range correlations through a Bragg–Williams approximation. The behavior of the dielectric function, the energy loss function, and the relaxation time of ferroelectric-like states is then investigated as a function of the temperature of the distribution and the damping constant of the electronic states in the single quantum dots. The robustness of the trends and tendencies of our results indicate that this scheme of analysis can guide experimentalists to develop tailored quantum dots distributions for specific applications.
A new model for two-dimensional numerical simulation of pseudo-2D gas-solids fluidized beds
Li, Tingwen; Zhang, Yongmin
2013-10-11
Pseudo-two dimensional (pseudo-2D) fluidized beds, for which the thickness of the system is much smaller than the other two dimensions, is widely used to perform fundamental studies on bubble behavior, solids mixing, or clustering phenomenon in different gas-solids fluidization systems. The abundant data from such experimental systems are very useful for numerical model development and validation. However, it has been reported that two-dimensional (2D) computational fluid dynamic (CFD) simulations of pseudo-2D gas-solids fluidized beds usually predict poor quantitative agreement with the experimental data, especially for the solids velocity field. In this paper, a new model is proposed to improve the 2D numerical simulations of pseudo-2D gas-solids fluidized beds by properly accounting for the frictional effect of the front and back walls. Two previously reported pseudo-2D experimental systems were simulated with this model. Compared to the traditional 2D simulations, significant improvements in the numerical predictions have been observed and the predicted results are in better agreement with the available experimental data.
Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions
NASA Technical Reports Server (NTRS)
Mikellides, Ioannis G.; Katz, Ira; Hofer, Richard H.; Goebel, Dan M.
2009-01-01
A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.
Nonlinear MHD simulation of DC helicity injection in the Pegasus spherical tokamak
NASA Astrophysics Data System (ADS)
Bayliss, Adam; Sovinec, Carl
2006-10-01
DC helicity injection has been successfully employed in spherical tokamaks (ST's) to produce a tokamak-like plasma with either a poloidal-gap voltage known as coaxial helicity injection [HIT-II, NSTX] or a biased cathode gun configuration [CDX, PEGASUS]. In PEGASUS, the tokamak-like plasma which is subsequently ohmically driven is the product of a reversal of vacuum poloidal flux and a merger of gun-injected current filaments. A 3D nonlinear MHD computation using the NIMROD code [Sovinec et al. JCP 195, 355 (2004)] simulates the formation, merger, and relaxation of the gun-injected current filaments to the tokamak-like plasma. The reversal of poloidal flux due to the field induced by the helicity drive is reproduced and the MHD processes leading to the merger and relaxation of the current filaments are described. Over the lifetime of a helically-driven experimental shot (approximately 10ms), the extent to which the merged plasma exhibits amplication of poloidal flux and the injected current in the relaxed state, reported in PEGASUS, is explored. The results are compared with simulations of current drive in NSTX via coaxial helicity injection which exhibit an n=1 open field-line kink [Tang and Boozer, Phys. Plasmas 11, 2679 (2004)].
Chatterjee, Dipankar; Amiroudine, Sakir
2011-02-01
A comprehensive non-isothermal Lattice Boltzmann (LB) algorithm is proposed in this article to simulate the thermofluidic transport phenomena encountered in a direct-current (DC) magnetohydrodynamic (MHD) micropump. Inside the pump, an electrically conducting fluid is transported through the microchannel by the action of an electromagnetic Lorentz force evolved out as a consequence of the interaction between applied electric and magnetic fields. The fluid flow and thermal characteristics of the MHD micropump depend on several factors such as the channel geometry, electromagnetic field strength and electrical property of the conducting fluid. An involved analysis is carried out following the LB technique to understand the significant influences of the aforementioned controlling parameters on the overall transport phenomena. In the LB framework, the hydrodynamics is simulated by a distribution function, which obeys a single scalar kinetic equation associated with an externally imposed electromagnetic force field. The thermal history is monitored by a separate temperature distribution function through another scalar kinetic equation incorporating the Joule heating effect. Agreement with analytical, experimental and other available numerical results is found to be quantitative. PMID:21053082
Numerical simulation of HTPB combustion in a 2D hybrid slab combustor
NASA Astrophysics Data System (ADS)
Gariani, Gabriela; Maggi, Filippo; Galfetti, Luciano
2011-09-01
A code for the numerical simulation of combustion processes in hybrid rockets, developed at the Space Propulsion Laboratory of Politecnico di Milano (SPLab), is presented. The code deals with Navier-Stokes equations solved with RANS approach, blowing effect, combustion kinetics and radiation. The equations are closed with k-epsilon turbulence model and well stirred reactor model. The P1 model, a simplification of the PN radiation model, is adopted. Specific simulation tools were developed using OpenFOAM®open source technology. The computational domain is 2D and split in two subdomains, simulating the reacting gas mixture on one side and the solid fuel grain on the other. The interface between the two regions plays a key role as the solid grain pyrolysis comes from a straight solution of the model without shortcuts. A propellant combination with polybutadiene and gaseous oxygen has been chosen and a reduced kinetic model for combustion of butadiene, considered as the major gaseous constituent coming from polybutadiene pyrolysis, has been developed for reactions occurring in oxygen atmosphere. The computational domain tries to replicate the real experimental setup and is split into three areas: pre-chamber, slab zone and post-chamber. High speed camera visualizations of the combustion processes allow to compare the flame height, obtained by the code and by experimental tests, along the grain for given boundary conditions.
Simulation of abrasive flow machining process for 2D and 3D mixture models
NASA Astrophysics Data System (ADS)
Dash, Rupalika; Maity, Kalipada
2015-12-01
Improvement of surface finish and material removal has been quite a challenge in a finishing operation such as abrasive flow machining (AFM). Factors that affect the surface finish and material removal are media viscosity, extrusion pressure, piston velocity, and particle size in abrasive flow machining process. Performing experiments for all the parameters and accurately obtaining an optimized parameter in a short time are difficult to accomplish because the operation requires a precise finish. Computational fluid dynamics (CFD) simulation was employed to accurately determine optimum parameters. In the current work, a 2D model was designed, and the flow analysis, force calculation, and material removal prediction were performed and compared with the available experimental data. Another 3D model for a swaging die finishing using AFM was simulated at different viscosities of the media to study the effects on the controlling parameters. A CFD simulation was performed by using commercially available ANSYS FLUENT. Two phases were considered for the flow analysis, and multiphase mixture model was taken into account. The fluid was considered to be a
What Can We Learn about Magnetotail Reconnection from 2D PIC Harris-Sheet Simulations?
NASA Astrophysics Data System (ADS)
Goldman, M. V.; Newman, D. L.; Lapenta, G.
2016-03-01
The Magnetosphere Multiscale Mission (MMS) will provide the first opportunity to probe electron-scale physics during magnetic reconnection in Earth's magnetopause and magnetotail. This article will address only tail reconnection—as a non-steady-state process in which the first reconnected field lines advance away from the x-point in flux pile-up fronts directed Earthward and anti-Earthward. An up-to-date microscopic physical picture of electron and ion-scale collisionless tail reconnection processes is presented based on 2-D Particle-In-Cell (PIC) simulations initiated from a Harris current sheet and on Cluster and Themis measurements of tail reconnection. The successes and limitations of simulations when compared to measured reconnection are addressed in detail. The main focus is on particle and field diffusion region signatures in the tail reconnection geometry. The interpretation of these signatures is vital to enable spacecraft to identify physically significant reconnection events, to trigger meaningful data transfer from MMS to Earth and to construct a useful overall physical picture of tail reconnection. New simulation results and theoretical interpretations are presented for energy transport of particles and fields, for the size and shape of electron and ion diffusion regions, for processes occurring near the fronts and for the j × B (Hall) electric field.
Lattice Boltzmann simulations of 2D laminar flows past two tandem cylinders
NASA Astrophysics Data System (ADS)
Mussa, Alberto; Asinari, Pietro; Luo, Li-Shi
2009-03-01
We apply the lattice Boltzmann equation (LBE) with multiple-relaxation-time (MRT) collision model to simulate laminar flows in two-dimensions (2D). In order to simulate flows in an unbounded domain with the LBE method, we need to address two issues: stretched non-uniform mesh and inflow and outflow boundary conditions. We use the interpolated grid stretching method to address the need of non-uniform mesh. We demonstrate that various inflow and outflow boundary conditions can be easily and consistently realized with the MRT-LBE. The MRT-LBE with non-uniform stretched grids is first validated with a number of test cases: the Poiseuille flow, the flow past a cylinder asymmetrically placed in a channel, and the flow past a cylinder in an unbounded domain. We use the LBE method to simulate the flow past two tandem cylinders in an unbounded domain with Re = 100. Our results agree well with existing ones. Through this work we demonstrate the effectiveness of the MRT-LBE method with grid stretching.
Application of 2-D simulations to hollow Z-pinch implosions
Peterson, D. L.; Bowers, R. L.; Brownell, J. H.; Lund, C.; Matuska, W.; McLenithan, K.; Oona, H.; Deeney, C.; Derzon, M.; Spielman, R. B.; Nash, T. J.; Chandler, G.; Mock, R. C.; Sanford, T. W. L.; Matzen, M. K.; Roderick, N. F.
1997-05-05
The application of simulations of z-pinch implosions should have at least two goals: first, to properly model the most important physical processes occurring in the pinch allowing for a better understanding of the experiments and second, provide a design capability for future experiments. Beginning with experiments fielded at Los Alamos on the Pegasus I and Pegasus II capacitor banks, we have developed a methodology for simulating hollow z-pinches in two dimensions which has reproduced important features of the measured experimental current drive, spectrum, radiation pulse shape, peak power and total radiated energy (1,2,3). This methodology employs essentially one free parameter, the initial level of the random density perturbations imposed at the beginning of the 2-D simulation, but in general no adjustments to other parameters (such as the resistivity) are required (1). Limitations in the use of this approach include the use of the 3-T, gray diffusion treatment of radiation and the fact that the initial perturbation conditions are not known a priori. Nonetheless, the approach has been successful in reproducing important experimental features of such implosions over a wide variety of timescales (tens of nanoseconds to microseconds), current drives (3 to 16 MA), masses (submilligram to tens of milligrams), initial radii (<1 cm to 5 cm), materials (Al and W) and initial configurations (thin foils and wire arrays with 40 to 240 wires). Currently we are applying this capability to the analysis of recent Saturn and PBFA-Z experiments (4,5). The code results provide insight into the nature of the pinch plasma prior to arrival on-axis, during thermalization and development after peak pinch time. Among other things, the simulation results provide an explanation for the production of larger amounts of radiated energy than would be expected from a simple slug-model kinetic energy analysis and the appearance of multiple peaks in the radiation power. The 2-D modeling has
Application of 2-D simulations to hollow Z-pinch implosions
Peterson, D.L.; Bowers, R.L.; Brownell, J.H.; Lund, C.; Matuska, W.; McLenithan, K.; Oona, H.; Deeney, C.; Derzon, M.; Spielman, R.B.; Nash, T.J.; Chandler, G.; Mock, R.C.; Sanford, T.W.; Matzen, M.K.; Roderick, N.F.
1997-05-01
The application of simulations of z-pinch implosions should have at least two goals: first, to properly model the most important physical processes occurring in the pinch allowing for a better understanding of the experiments and second, provide a design capability for future experiments. Beginning with experiments fielded at Los Alamos on the Pegasus I and Pegasus II capacitor banks, we have developed a methodology for simulating hollow z-pinches in two dimensions which has reproduced important features of the measured experimental current drive, spectrum, radiation pulse shape, peak power and total radiated energy (1,2,3). This methodology employs essentially one free parameter, the initial level of the random density perturbations imposed at the beginning of the 2-D simulation, but in general no adjustments to other parameters (such as the resistivity) are required (1). Limitations in the use of this approach include the use of the 3-T, gray diffusion treatment of radiation and the fact that the initial perturbation conditions are not known {ital a priori}. Nonetheless, the approach has been successful in reproducing important experimental features of such implosions over a wide variety of timescales (tens of nanoseconds to microseconds), current drives (3 to 16 MA), masses (submilligram to tens of milligrams), initial radii ({lt}1cm to 5 cm), materials (Al and W) and initial configurations (thin foils and wire arrays with 40 to 240 wires). Currently we are applying this capability to the analysis of recent Saturn and PBFA-Z experiments (4,5). The code results provide insight into the nature of the pinch plasma prior to arrival on-axis, during thermalization and development after peak pinch time. Among other things, the simulation results provide an explanation for the production of larger amounts of radiated energy than would be expected from a simple slug-model kinetic energy analysis and the appearance of multiple peaks in the radiation power. The 2-D
One year in the Earth's magnetosphere: A global MHD simulation and spacecraft measurements
NASA Astrophysics Data System (ADS)
Facskó, G.; Honkonen, I.; Živković, T.; Palin, L.; Kallio, E.; Ã gren, K.; Opgenoorth, H.; Tanskanen, E. I.; Milan, S.
2016-05-01
The response of the Earth's magnetosphere to changing solar wind conditions is studied with a 3-D Magnetohydrodynamic (MHD) model. One full year (155 Cluster orbits) of the Earth's magnetosphere is simulated using Grand Unified Magnetosphere Ionosphere Coupling simulation (GUMICS-4) magnetohydrodynamic code. Real solar wind measurements are given to the code as input to create the longest lasting global magnetohydrodynamics simulation to date. The applicability of the results of the simulation depends critically on the input parameters used in the model. Therefore, the validity and the variance of the OMNIWeb data are first investigated thoroughly using Cluster measurement close to the bow shock. The OMNIWeb and the Cluster data were found to correlate very well before the bow shock. The solar wind magnetic field and plasma parameters are not changed significantly from the L1 Lagrange point to the foreshock; therefore, the OMNIWeb data are appropriate input to the GUMICS-4. The Cluster SC3 footprints are determined by magnetic field mapping from the simulation results and the Tsyganenko (T96) model in order to compare two methods. The determined footprints are in rather good agreement with the T96. However, it was found that the footprints agree better in the Northern Hemisphere than the Southern one during quiet conditions. If the By is not zero, the agreement of the GUMICS-4 and T96 footprint is worse in longitude in the Southern Hemisphere. Overall, the study implies that a 3-D MHD model can increase our insight of the response of the magnetosphere to solar wind conditions.
Nonlinear MHD simulations of Quiescent H-mode plasmas in DIII-D
NASA Astrophysics Data System (ADS)
Liu, F.; Huijsmans, G. T. A.; Loarte, A.; Garofalo, A. M.; Solomon, W. M.; Snyder, P. B.; Hoelzl, M.; Zeng, L.
2015-09-01
In the Quiescent H-mode (QH-mode) regime, the edge harmonic oscillation (EHO), thought to be a saturated kink-peeling mode (KPM) driven unstable by current and rotation, is found in experiment to provide sufficient stationary edge particle transport to avoid the periodic expulsion of particles and energy by edge localized modes (ELMs). In this paper, both linear and nonlinear MHD modelling of QH-mode plasmas from the DIII-D tokamak have been investigated to understand the mechanism leading to the appearance of the EHO in QH-mode plasmas. For the first time nonlinear MHD simulations with low-n modes both with ideal wall and resistive wall boundary conditions have been carried out with the 3D non-linear MHD code JOREK. The results show, in agreement with the original conjectures, that in the non-linear phase, kink peeling modes are the main unstable modes in QH-mode plasmas of DIII-D and that the kink-peeling modes saturate non-linearly leading to a 3D stationary state. The characteristics of the kink-peeling modes, in terms of mode structure and associated decrease of the edge plasma density associated with them, are in good agreement with experimental measurements of the EHO in DIII-D. The effect of plasma resistivity, the role of plasma parallel rotation as well as the effect of the conductivity of the vacuum vessel wall on the destabilization and saturation of kink-peeling modes have been evaluated for experimental QH-mode plasma conditions in DIII-D.
NASA Astrophysics Data System (ADS)
Xi, S.; Lotko, W.; Zhang, B.; Brambles, O.; Wiltberger, M. J.; Lyon, J.; Merkin, V. G.
2010-12-01
In global modeling, magnetosphere-ionosphere (MI) coupling physically connects a global magnetospheric (GM) model and a global ionospheric-thermospheric (GIT) model. The field-aligned current from the GM model and the conductance distributions from the GIT model are used in a Poisson equation derived from the ionospheric Ohm's law combined with current continuity to determine the electrostatic potential in the ionosphere. In current GM models, this electrostatic potential is mapped to the inner boundary of the GM simulation to determine electrostatic boundary conditions on the electric field and MHD velocity there. Inductive effects and the finite Alfven transit time between the low-altitude GM boundary and the high-altitude GIT boundary (MI gap region) are neglected in this formulation of MI coupling. Using fields and currents derived from Lyon-Fedder-Mobarry GM simulations, and conductance distributions derived from its standalone empirical conductance model in the MI coupling Poisson equation, we have computed the fast Fourier transform of the electrostatic field at the low-altitude LFM simulation boundary as described above, and the FFT of the inductive electric field at the boundary under the assumption that μ 0 Σ P vA ≤ 1, where Σ P is the ionospheric Pedersen conductance and vA is the smallest value of the Alfven speed in the MI gap region. In this regime, the complete electric field at the low-altitude simulation boundary includes the usual mapped electrostatic field with an inductive addition for which the finite Alfven transit time and the diversion of field-aligned into polarization currents in the gap region are negligible (Lotko, 2004). By comparing the boundary-averaged spectra of the electrostatic and so-determined inductive fields, we confirm that the purely electrostatic formulation of MI coupling is valid when the MHD state varies on times scales exceeding about 200 s. For faster MHD time variations, the inductive electric field is shown to
NASA Technical Reports Server (NTRS)
Benyo, Theresa L.
2011-01-01
Flow matching has been successfully achieved for an MHD energy bypass system on a supersonic turbojet engine. The Numerical Propulsion System Simulation (NPSS) environment helped perform a thermodynamic cycle analysis to properly match the flows from an inlet employing a MHD energy bypass system (consisting of an MHD generator and MHD accelerator) on a supersonic turbojet engine. Working with various operating conditions (such as the applied magnetic field, MHD generator length and flow conductivity), interfacing studies were conducted between the MHD generator, the turbojet engine, and the MHD accelerator. This paper briefly describes the NPSS environment used in this analysis. This paper further describes the analysis of a supersonic turbojet engine with an MHD generator/accelerator energy bypass system. Results from this study have shown that using MHD energy bypass in the flow path of a supersonic turbojet engine increases the useful Mach number operating range from 0 to 3.0 Mach (not using MHD) to a range of 0 to 7.0 Mach with specific net thrust range of 740 N-s/kg (at ambient Mach = 3.25) to 70 N-s/kg (at ambient Mach = 7). These results were achieved with an applied magnetic field of 2.5 Tesla and conductivity levels in a range from 2 mhos/m (ambient Mach = 7) to 5.5 mhos/m (ambient Mach = 3.5) for an MHD generator length of 3 m.
NASA Astrophysics Data System (ADS)
Yamada, Susumu; Kitamura, Akihiro; Kurikami, Hiroshi; Machida, Masahiko
2015-04-01
Fukushima Daiichi Nuclear Power Plant (FDNPP) accident on March 2011 released significant quantities of radionuclides to atmosphere. The most significant nuclide is radioactive cesium isotopes. Therefore, the movement of the cesium is one of the critical issues for the environmental assessment. Since the cesium is strongly sorbed by soil particles, the cesium transport can be regarded as the sediment transport which is mainly brought about by the aquatic system such as a river and a lake. In this research, our target is the sediment transport on Ogaki dam reservoir which is located in about 16 km northwest from FDNPP. The reservoir is one of the principal irrigation dam reservoirs in Fukushima Prefecture and its upstream river basin was heavily contaminated by radioactivity. We simulate the sediment transport on the reservoir using 2-D river simulation code named Nays2D originally developed by Shimizu et al. (The latest version of Nays2D is available as a code included in iRIC (http://i-ric.org/en/), which is a river flow and riverbed variation analysis software package). In general, a 2-D simulation code requires a huge amount of calculation time. Therefore, we parallelize the code and execute it on a parallel computer. We examine the relationship between the behavior of the sediment transport and the height of the reservoir exit. The simulation result shows that almost all the sand that enter into the reservoir deposit close to the entrance of the reservoir for any height of the exit. The amounts of silt depositing within the reservoir slightly increase by raising the height of the exit. However, that of the clay dramatically increases. Especially, more than half of the clay deposits, if the exit is sufficiently high. These results demonstrate that the water level of the reservoir has a strong influence on the amount of the clay discharged from the reservoir. As a result, we conclude that the tuning of the water level has a possibility for controlling the
The PLX- α project: Radiation-MHD Simulations of Imploding Plasma Liners Using USim
NASA Astrophysics Data System (ADS)
Beckwith, Kristian; Stoltz, Peter; Kundrapu, Madhusudhan; Hsu, Scott; PLX-α Team
2015-11-01
USim is a tool for modeling high energy density plasmas using multi-fluid models coupled to electromagnetics using fully-implicit iterative solvers, combined with finite volume discretizations on unstructured meshes. Prior work has demonstrated application of USim models and algorithms to simulation of supersonic plasma jets relevant to the Plasma Liner Experiment (PLX) and compared synthetic interferometry to that gathered from the experiment. Here, we give an overview of the models and algorithms included in USim; review results from prior modeling campaigns for the PLX; and describe plans for radiation magnetohydrodynamic (MHD) simulation efforts focusing on integrated plasma-liner implosion and target compression in a fusion-relevant regime using USim for the PLX- α project. Supported by ARPA-E's ALPHA program. Original PLX construction supported by OFES. USim development supported in part by Air Force Office of Scientific Research.
Interpreting Irradiance Distributions Using High-Resolution 3D MHD Simulations
NASA Astrophysics Data System (ADS)
Peck, Courtney; Rast, Mark; Criscuoli, Serena; Uitenbroek, Han; Rempel, Matthias D.
2016-05-01
We present initial results of studies aimed at understanding the impact of the unresolved magnetic field distribution on solar spectral irradiance. Using high-resolution 3D MHD simulations (from MURaM code) and spectral synthesis (with the RH code), we examine the emergent spectra of two atmospheres with similar mean field strengths but differing imposed-field conditions at wavelengths spanning from visible to infrared. Comparing the contrast against the magnetic field strength for the two magnetic simulations, we find differences in the distributions of contrasts versus field strength. We repeat the analysis after convolving the images with the PSF of a typical solar telescope (1-meter) and discuss the potential implications for irradiance modeling and future steps.
The role of condensation and heat conduction in the formation of prominences - An MHD simulation
NASA Technical Reports Server (NTRS)
Wu, S. T.; Bao, J. J.; An, C. H.; Tandberg-Hanssen, E.
1990-01-01
The effects of condensation and thermal conduction on the formation of Kippenhahn-Schlueter (K-S) type prominences in quiet regions (QP) due to symmetric mass injection are investigated. To implement this investigation a self-consistent, two-dimensional, nonplanar, time-dependent MHD simulation model is developed. In the model, various values of the injection velocity, density, and magnetic field strength are used to determine the most favorable conditions for the QP formation. Based on these simulation results, it is found that the formation of a K-S type field configuration should be considered as a dynamic process which needs both condensation amd mass injection to supply enough mass to maintain such a configuration to complete the quiescent prominence formation process.
1D and 2D simulations of seismic wave propagation in fractured media
NASA Astrophysics Data System (ADS)
Möller, Thomas; Friederich, Wolfgang
2016-04-01
Fractures and cracks have a significant influence on the propagation of seismic waves. Their presence causes reflections and scattering and makes the medium effectively anisotropic. We present a numerical approach to simulation of seismic waves in fractured media that does not require direct modelling of the fracture itself, but uses the concept of linear slip interfaces developed by Schoenberg (1980). This condition states that at an interface between two imperfectly bonded elastic media, stress is continuous across the interface while displacement is discontinuous. It is assumed that the jump of displacement is proportional to stress which implies a jump in particle velocity at the interface. We use this condition as a boundary condition to the elastic wave equation and solve this equation in the framework of a Nodal Discontinuous Galerkin scheme using a velocity-stress formulation. We use meshes with tetrahedral elements to discretise the medium. Each individual element face may be declared as a slip interface. Numerical fluxes have been derived by solving the 1D Riemann problem for slip interfaces with elastic and viscoelastic rheology. Viscoelasticity is realised either by a Kelvin-Voigt body or a Standard Linear Solid. These fluxes are not limited to 1D and can - with little modification - be used for simulations in higher dimensions as well. The Nodal Discontinuous Galerkin code "neXd" developed by Lambrecht (2013) is used as a basis for the numerical implementation of this concept. We present examples of simulations in 1D and 2D that illustrate the influence of fractures on the seismic wavefield. We demonstrate the accuracy of the simulation through comparison to an analytical solution in 1D.
Simulation and analysis of solute transport in 2D fracture/pipe networks: The SOLFRAC program
NASA Astrophysics Data System (ADS)
Bodin, Jacques; Porel, Gilles; Delay, Fred; Ubertosi, Fabrice; Bernard, Stéphane; de Dreuzy, Jean-Raynald
2007-01-01
The Time Domain Random Walk (TDRW) method has been recently developed by Delay and Bodin [Delay, F. and Bodin, J., 2001. Time domain random walk method to simulate transport by advection-dispersion and matrix diffusion in fracture networks. Geophys. Res. Lett., 28(21): 4051-4054.] and Bodin et al. [Bodin, J., Porel, G. and Delay, F., 2003c. Simulation of solute transport in discrete fracture networks using the time domain random walk method. Earth Planet. Sci. Lett., 6566: 1-8.] for simulating solute transport in discrete fracture networks. It is assumed that the fracture network can reasonably be represented by a network of interconnected one-dimensional pipes (i.e. flow channels). Processes accounted for are: (1) advection and hydrodynamic dispersion in the channels, (2) matrix diffusion, (3) diffusion into stagnant zones within the fracture planes, (4) sorption reactions onto the fracture walls and in the matrix, (5) linear decay, and (6) mass sharing at fracture intersections. The TDRW method is handy and very efficient in terms of computation costs since it allows for the one-step calculation of the particle residence time in each bond of the network. This method has been programmed in C++, and efforts have been made to develop an efficient and user-friendly software, called SOLFRAC. This program is freely downloadable at the URL http://labo.univ-poitiers.fr/hydrasa/intranet/telechargement.htm. It calculates solute transport into 2D pipe networks, while considering different types of injections and different concepts of local dispersion within each flow channel. Post-simulation analyses are also available, such as the mean velocity or the macroscopic dispersion at the scale of the entire network. The program may be used to evaluate how a given transport mechanism influences the macroscopic transport behaviour of fracture networks. It may also be used, as is the case, e.g., with analytical solutions, to interpret laboratory or field tracer test experiments
2-D Three Fluid Simulation of Upstreaming Ions Above Auroral Precipitation
NASA Astrophysics Data System (ADS)
Danielides, M. A.; Lummerzheim, D.; Otto, A.; Stevens, R. J.
2006-12-01
The ionosphere is a rich reservoir of charged particles from which a variable fraction is transported to the magnetosphere. An important transport phenomena is the formation of upward ion flow above auroral structure. A primary region of the outflow is not known, but contributions come from polar cap, dayside cusp/cleft region, auroral oval, or even from mid-latitudes. In the past global magnetospheric models and fluid codes were used to simulate large scale ion outflow above, e.g., the polar-cap aurora. However, satellites orbiting at low- altitudes have repeatingly detected localized ion outflow above the auroral oval. Ionosphere-magnetosphere coupling simulations gave first insides into the small-scale dynamics of aurora. The aim of this study is the investigation of coupled plasma and neutral dynamics in smaller scale aurora to explain the generation, structure, and dynamics of vertical ion upstream. We consider auroral electron precipitation at ionospheric heights in a 2-D three fluid ionospheric-magnetospheric coupling code (Otto and Zhu, 2003). Specially we examine the effects of the electron precipitation, heat conduction and heating in field- aligned current through coulomb collisions or turbulence causing: i) electron heating, ii) electron pressure gradients, and iii) upstreaming of ions through a resulting ambipolar electric field. Our first case studies are performed for different boundary conditions and for different auroral electron precipitation parameters (variation in characteristic auroral energy, auroral energy flux and horizontal scale). The results shall clarify how auroral precipitation can drive ions upwards. Finally we discuss the effect of ion drag and the interaction of the upstreaming ions with a stable neutral constituent. Otto, O. and H. Zhu, Fluid plasma simulation of coupled systems: Ionosphere and magnetosphere, Space Plasma Simulation. Edited by J. Buechner, C. Dum, and M. Scholer., Lecture Notes in Physics, vol. 615, p.193
2D properties of core turbulence on DIII-D and comparison to gyrokinetic simulations
Shafer, Morgan W; Fonck, R. J.; McKee, G. R.; Holland, Chris; White, A. E.; Schlossberg, D J
2012-01-01
Quantitative 2D characteristics of localized density fluctuations are presented over the range of 0.3 < r/a < 0.9 in L-mode plasmas on DIII-D [J. L. Luxon, Nucl. Fusion 42, 614 (2002)]. Broadband density fluctuations increase in amplitude from (n) over tilde/n < 0.5% in the deep core to (n) over tilde/n similar to 2.5% near the outer region. The observed Doppler-shift due to the E x B velocity matches well with the measured turbulence group and phase velocities (in toroidally rotating neutral beam heated plasmas). Turbulence decorrelation rates are found to be similar to 200 kHz at the edge and to decrease toward the core (0.45 < r/a < 0.9) where they approach the E x B shearing rate (similar to 50 kHz). Radial and poloidal correlation lengths are found to scale with the ion gyroradius and exhibit an asymmetric poloidally elongated eddy structure. The ensemble-averaged turbulent eddy structure changes its tilt with respect to the radial-poloidal coordinates in the core, consistent with an E x B shear mechanism. The 2D spatial correlation and wavenumber spectra [S(k(r); k(theta))] are presented and compared to nonlinear flux-tube GYRO simulations at two radii, r/a = 0.5 and r/a = 0.75, showing reasonable overall agreement, but the GYRO spectrum exhibits a peak at finite kr for r/a = 0.75 that is not observed experimentally; E x B shear may cause this discrepancy. (C) 2012 American Institute of Physics.
2D IR spectra of cyanide in water investigated by molecular dynamics simulations
Lee, Myung Won; Carr, Joshua K.; Göllner, Michael; Hamm, Peter; Meuwly, Markus
2013-01-01
Using classical molecular dynamics simulations, the 2D infrared (IR) spectroscopy of CN− solvated in D2O is investigated. Depending on the force field parametrizations, most of which are based on multipolar interactions for the CN− molecule, the frequency-frequency correlation function and observables computed from it differ. Most notably, models based on multipoles for CN− and TIP3P for water yield quantitatively correct results when compared with experiments. Furthermore, the recent finding that T 1 times are sensitive to the van der Waals ranges on the CN− is confirmed in the present study. For the linear IR spectrum, the best model reproduces the full widths at half maximum almost quantitatively (13.0 cm−1 vs. 14.9 cm−1) if the rotational contribution to the linewidth is included. Without the rotational contribution, the lines are too narrow by about a factor of two, which agrees with Raman and IR experiments. The computed and experimental tilt angles (or nodal slopes) α as a function of the 2D IR waiting time compare favorably with the measured ones and the frequency fluctuation correlation function is invariably found to contain three time scales: a sub-ps, 1 ps, and one on the 10-ps time scale. These time scales are discussed in terms of the structural dynamics of the surrounding solvent and it is found that the longest time scale (≈10 ps) most likely corresponds to solvent exchange between the first and second solvation shell, in agreement with interpretations from nuclear magnetic resonance measurements.
2D IR spectra of cyanide in water investigated by molecular dynamics simulations.
Lee, Myung Won; Carr, Joshua K; Göllner, Michael; Hamm, Peter; Meuwly, Markus
2013-08-01
Using classical molecular dynamics simulations, the 2D infrared (IR) spectroscopy of CN(-) solvated in D2O is investigated. Depending on the force field parametrizations, most of which are based on multipolar interactions for the CN(-) molecule, the frequency-frequency correlation function and observables computed from it differ. Most notably, models based on multipoles for CN(-) and TIP3P for water yield quantitatively correct results when compared with experiments. Furthermore, the recent finding that T1 times are sensitive to the van der Waals ranges on the CN(-) is confirmed in the present study. For the linear IR spectrum, the best model reproduces the full widths at half maximum almost quantitatively (13.0 cm(-1) vs. 14.9 cm(-1)) if the rotational contribution to the linewidth is included. Without the rotational contribution, the lines are too narrow by about a factor of two, which agrees with Raman and IR experiments. The computed and experimental tilt angles (or nodal slopes) α as a function of the 2D IR waiting time compare favorably with the measured ones and the frequency fluctuation correlation function is invariably found to contain three time scales: a sub-ps, 1 ps, and one on the 10-ps time scale. These time scales are discussed in terms of the structural dynamics of the surrounding solvent and it is found that the longest time scale (≈10 ps) most likely corresponds to solvent exchange between the first and second solvation shell, in agreement with interpretations from nuclear magnetic resonance measurements. PMID:23927269
2D IR spectra of cyanide in water investigated by molecular dynamics simulations
NASA Astrophysics Data System (ADS)
Lee, Myung Won; Carr, Joshua K.; Göllner, Michael; Hamm, Peter; Meuwly, Markus
2013-08-01
Using classical molecular dynamics simulations, the 2D infrared (IR) spectroscopy of CN- solvated in D2O is investigated. Depending on the force field parametrizations, most of which are based on multipolar interactions for the CN- molecule, the frequency-frequency correlation function and observables computed from it differ. Most notably, models based on multipoles for CN- and TIP3P for water yield quantitatively correct results when compared with experiments. Furthermore, the recent finding that T1 times are sensitive to the van der Waals ranges on the CN- is confirmed in the present study. For the linear IR spectrum, the best model reproduces the full widths at half maximum almost quantitatively (13.0 cm-1 vs. 14.9 cm-1) if the rotational contribution to the linewidth is included. Without the rotational contribution, the lines are too narrow by about a factor of two, which agrees with Raman and IR experiments. The computed and experimental tilt angles (or nodal slopes) α as a function of the 2D IR waiting time compare favorably with the measured ones and the frequency fluctuation correlation function is invariably found to contain three time scales: a sub-ps, 1 ps, and one on the 10-ps time scale. These time scales are discussed in terms of the structural dynamics of the surrounding solvent and it is found that the longest time scale (≈10 ps) most likely corresponds to solvent exchange between the first and second solvation shell, in agreement with interpretations from nuclear magnetic resonance measurements.
2D simulation of transport and degradation in the River Rhine.
Teichmann, L; Reuschenbach, P; Müller, B; Horn, H
2002-01-01
A simple 2D model has been developed for the simulation of mass transport and degradation of substances in the river Rhine. The model describes mass transport in the flow direction with a convective and a dispersive term. Transversal transport is described by segmenting the river and formulating a transversal exchange coefficient between the segments. Degradation can be formulated with any kinetics from first order to complex enzyme kinetics. The model was verified with monitoring data from the river Rhine. The hydrodynamic parameters such as dispersion coefficients and exchange coefficients were fitted to the conductivity, which was assumed to be non-degradable. The degradation term was fitted to ammonia values. The model was used to simulate measured concentrations of a readily (Aniline) and a poorly biodegradable substance (1,4-Dioxan) 10 m from the left river bank. It was the objective of this research program to develop a model which allows a realistic estimation of the locally and regionally predicted environmental concentration of chemical substances in the EU risk assessment scheme. PMID:12380980
Simulations of SH wave scattering due to cracks by the 2-D finite difference method
NASA Astrophysics Data System (ADS)
Suzuki, Y.; Kawahara, J.; Okamoto, T.; Miyashita, K.
2006-05-01
We simulate SH wave scattering by 2-D parallel cracks using the finite difference method (FDM), instead of the popularly used boundary integral equation method (BIEM). Here special emphasis is put on simplicity; we apply a standard FDM (fourth-order velocity-stress scheme with a staggered grid) to media in cluding traction-freecracks, which are expressed by arrays of grid points with zero traction. Two types of accuracy tests based oncomparison with a reliable BIEM, suggest that the present method gives practically sufficient accuracy, except for the wavefields in the vicinity of cracks, which can be well handled if the second-order FDM is used instead. As an application of this method, we also simulate wave propagation in media with randomly distributed cracks of the same length. We experimentally determine the attenuation and velocity dispersion induced by scattering from the synthetic seismograms, using a waveform averaging technique. It is shown that the results are well explained by a theory based on the Foldy approximation for crack densities of up to about 01. The presence of a free surface does not affect the validity of the theory. A preliminary experiment also suggests that the validity will not change even for multi-scale cracks.
2D Mesoscale Simulation of Shock Response of Dry Sand in Plate Impact Experiments
NASA Astrophysics Data System (ADS)
Pei, L.; Teeter, R. D.; Dwivedi, S. K.; Gupta, Y. M.
2007-06-01
The one-dimensional approach with a homogenized continuum model used in the literature to derive the shock Hugoniot of sand from plate impact experimental data neglects heterogeneous deformation and cannot incorporate mesoscale phenomena. We present a 2D mesoscale simulation approach to probe the shock response of dry sand with the main objectives to identify important mesoscale phenomena and the role of inter granular friction. The in-house code ISP-SAND was used to generate sand with desired grain size distribution and porosity. The explicit finite element code ISP-TROTP was used to simulate plate impact experiments of assumed configurations. The deformation of individual sand grains was modeled by non-linear mean stress volume compression relation with an assumed mean stress dependent yield strength. The results show heterogeneous deformation with finite lateral velocity and regions of stress concentrations in the sand sample. The effects of grain size distribution, porosity and friction between grains are discussed by comparing the particle velocity profiles at the window interface. Work supported by DOE and AFOSR.
OpenGGCM Simulation of Ballooning and Axial MHD Mode at Substorm Onset (Invited)
NASA Astrophysics Data System (ADS)
Raeder, J.; Zhu, P.; Ge, Y.; Siscoe, G. L.
2010-12-01
It is generally accepted that magnetic reconnection is the main mechanism that dissipates power during a substorm. It is less clear, however, whether the beginning of magnetic reconnection in the magnetotail also signifies the onset of the substorm expansion phase itself, i.e., whether the "outside-in" scenario applies, or if a different process happens first closer to Earth that triggers the reconnection onset in the magnetotail, i.e., the "inside-out" scenario. Global MHD simulations have generally supported the "outside-in" scenario. However, ideal MHD instabilities that could possibly trigger tail reconnection may have been missed due to coarse numerical resolution or due to other numerical effects. Here, we present results from an OpenGGCM simulation of the March 23, 2007 substorm that clearly shows growth of the ballooning mode as suggested by our earlier analysis (Zhu et al., 2009), as well as growth of an ideal-like instability that is purely axial and was previously reported by Siscoe et al. (2009). Both instabilities occur simultaneously and are immediately followed by reconnection onset. The simulations results are in accordance with recent Geotail observations of ballooning with a wavelength of approximately 0.5 RE, and the time scales agree with that of the explosive growth phase. The exact relation between the three instabilities, i.e., ballooning, the axial mode, and tearing, is not entirely clear yet; however, having demonstrated that the OpenGGCM reproduces all of them, they can now be analyzed in more detail. Furthermore, we present the expected auroral signature of these processes, which will allow for a more detailed comparison of the simulation results with ground based imagers.
NASA Astrophysics Data System (ADS)
Jung, J.; Arakawa, A.
2015-12-01
Through explicitly resolved cloud-scale processes by embedded 2-D cloud-resolving models (CRMs), the Multiscale Modeling Framework (MMF) known as the superparameterization has been reasonably successful to simulate various atmospheric events over a wide range of time scales. One thing to be justified is, however, if the influence of complex 3-D topography can be adequately represented by the embedded 2-D CRMs. In this study, simulations are performed in the presence of a variety of topography with embedded 3-D and 2-D CRMs in a single-column inactive GCM. Through the comparison between these simulations, it is demonstrated that the 2-D representation of topography is able to simulate the statistics of precipitation due to 3-D topography reasonably well as long as the topographic characteristics, such as the mean and standard deviation, are closely recognized. It is also shown that the use of two perpendicular sets of 2-D representations tends to reduce the error due to a 2-D representation.
MHD SIMULATIONS OF ACCRETION ONTO Sgr A*: QUIESCENT FLUCTUATIONS, OUTBURSTS, AND QUASIPERIODICITY
Chan Chikwan; Liu Siming; Fryer, Christopher L.; Psaltis, Dimitrios; Oezel, Feryal; Melia, Fulvio; Rockefeller, Gabriel
2009-08-10
High-resolution observations of Sgr A* have revealed a wide variety of phenomena, ranging from intense rapid flares to quasi-periodic oscillations (QPOs), making this object an ideal system to study the properties of low luminosity accreting black holes. In this paper, we use a pseudospectral algorithm to construct and evolve a three-dimensional magnetohydrodynamic (MHD) model of the accretion disk in Sgr A*. Assuming a hybrid thermal-nonthermal emission scheme and calibrating the parameters by observations, we show that the MHD turbulence in the environment of Sgr A* can by itself only produce factor two fluctuations in luminosity. These fluctuations cannot explain the magnitude of flares observed in this system. However, we also demonstrate that external forcing of the accretion disk, which may be generated by the 'clumpy material' raining down onto the disk from the large-scale flow, do produce outbursts qualitatively similar to those observed by XMM-Newton in X-rays and by ground-based facilities in the near infrared. Strong, but short-term QPOs emerge naturally in the simulated light curves. We attribute these to nonaxisymmetric density perturbations that emerge as the disk evolves back toward its quiescent state.
Sub-Alfvenic Non-Ideal MHD Turbulence Simulations with Ambipolar Diffusion: I. Turbulence Statistics
Klein, R I; Li, P S; McKee, C F; Fisher, R
2008-04-10
Most numerical investigations on the role of magnetic fields in turbulent molecular clouds (MCs) are based on ideal magneto-hydrodynamics (MHD). However, MCs are weakly ionized, so that the time scale required for the magnetic field to diffuse through the neutral component of the plasma by ambipolar diffusion (AD) can be comparable to the dynamical time scale. We have performed a series of 256{sup 3} and 512{sup 3} simulations on supersonic but sub-Alfvenic turbulent systems with AD using the Heavy-Ion Approximation developed in Li et al. (2006). Our calculations are based on the assumption that the number of ions is conserved, but we show that these results approximately apply to the case of time-dependent ionization in molecular clouds as well. Convergence studies allow us to determine the optimal value of the ionization mass fraction when using the heavy-ion approximation for low Mach number, sub-Alfvenic turbulent systems. We find that ambipolar diffusion steepens the velocity and magnetic power spectra compared to the ideal MHD case. Changes in the density PDF, total magnetic energy, and ionization fraction are determined as a function of the AD Reynolds number. The power spectra for the neutral gas properties of a strongly magnetized medium with a low AD Reynolds number are similar to those for a weakly magnetized medium; in particular, the power spectrum of the neutral velocity is close to that for Burgers turbulence.
Two-Species, 3D, MHD Simulation of Europa's Interaction with Jupiter's Magnetosphere
NASA Technical Reports Server (NTRS)
Liu, Yifan; Nagy, Andrew F.; Kabin, Konstantin; Combi, Michael R.; DeZeeuw, Darren L.; Gombosi, Tamas I.; Powell, Kenneth G.
2000-01-01
The interaction of Europa with the Jovian a magnetosphere has been studied by using a two species in ideal magnetohydrodynamic (MHD) numerical model. This model considers the upstream plasma in the Jovian magnetosphere and the molecular oxygen ions in the ionosphere of Europa, separately. We present results a from simulation studies, which take into account impact ionization, recombination, and the effect of a possible induced dipole magnetic field of Europa. The total mass loading of the magnetospheric flow and the ionization frequency used in the model are consistent with the estimates of Europa's ionosphere and atmosphere. The multi-species MHD equations are solved by using a finite volume, high-order, Godunov-type method on an adoptively refined unstructured grid, which allows a detailed modeling of the region near Europa's surface, while still resolving both the upstream region and the satellite's wake. We have paid special attention to the wake of Europa, in order to be able to make comparisons with the Galileo's E4 flyby observations, as well as other model calculations. The calculated escape flux of a O2+ down the tail was found to be about 5.6 x 10(exp 25) s(sup -1).
On Europa's Magnetospheric Interaction: A MHD Simulation of the E4 Flyby
NASA Technical Reports Server (NTRS)
Kabin, K.; Combi, M. R.; Gombosi, T. I.; Nagy, A. F.; DeZeeuw, D. L.; Powell, K. G.
1999-01-01
The global three-dimensional interaction of Europa with the Jovian magnetosphere is modeled by using a complete set of ideal magnetohydrodynamic (MHD) equations. The model accounts for exospheric mass loading, ion-neutral charge exchange, recombination, and a possible intrinsic dipole magnetic field of Europa. The single-fluid MHD equations are solved by using a modem, finite volume, higher-order, Godunov-type method on an adoptively refined unstructured grid, which allows detailed modeling of the region near Europa while still resolving both the upstream region and the satellite's wake. The magnetic field and plasma density measured during Galileo's EGA flyby of December 19, 1996, are reproduced reasonably well in the simulation. We find the agreement between the data and our model particularly convincing if we assume that the plasma velocity during the EGA flyby deviated from the nominal corotation direction by approximately 20 deg. Evidence from the Galileo energetic particle detector also supports this assumption. In this case, we can fit the data using a dipole with orientation close to that of an induced dipole arising from the interaction of a hypothetical conducting subsurface layer on Europa with the periodically changing magnetic field of Jupiter. However, the magnitude of the dipole in our model is somewhat smaller (70%) than that suggested by Khurana et al. The total mass loading and ion-neutral charge exchange rates are consistent with the estimates of Europa's atmosphere and ionosphere.
NASA Astrophysics Data System (ADS)
Pantellini, Filippo; Griton, Léa; Varela, Jacobo
2015-07-01
We show that slow mode compressional fronts form upstream of the day side magnetopause in MHD simulations of Mercury's magnetosphere. The strongest compressional fronts are located upstream of the magnetopause with strong magnetic shear. Compressional fronts are crossed by magnetic field lines connecting the interplanetary magnetic field and the planet's intrinsic field, their role is to bend the magnetic field in the magnetosheath towards the magnetopause. Besides these compressional fronts, already observed in space and theoretically discussed by various authors for the case of the Earth, we observe the formation of a slow mode standing rarefaction wave spatially growing over a substantial fraction of the distance between the bow shock and the magnetopause. The slow mode source region for the rarefaction waves is located in the magnetosheath, near the bow shock's nose. The generated standing rarefaction waves, however, form even at large distances from the source region along the magnetospheric flanks. They fine-tune the magnetic field line draping and plasma flow around the magnetopause. In ideal MHD the magnetospheres of Mercury, the Earth and the giant planets do closely resemble each other, we therefore expect the mentioned slow mode structures not to be specific to Mercury.
NASA Astrophysics Data System (ADS)
Kanki, Takashi; Nagata, Masayoshi; Kagei, Yasuhiro
2011-10-01
The dynamics of structures of magnetic field, current density, and plasma flow generated during multi-pulsed coaxial helicity injection in spherical torus is investigated by 3-D nonlinear MHD simulations. During the driven phase, the flux and current amplifications occur due to the merging and magnetic reconnection between the preexisting plasma in the confinement region and the ejected plasma from the gun region involving the n = 1 helical kink distortion of the central open flux column (COFC). Interestingly, the diamagnetic poloidal flow which tends toward the gun region is then observed due to the steep pressure gradients of the COFC generated by ohmic heating through an injection current winding around the inboard field lines, resulting in the formation of the strong poloidal flow shear at the interface between the COFC and the core region. This result is consistent with the flow shear observed in the HIST. During the decay phase, the configuration approaches the axisymmetric MHD equilibrium state without flow because of the dissipation of magnetic fluctuation energy to increase the closed flux surfaces, suggesting the generation of ordered magnetic field structure. The parallel current density λ concentrated in the COFC then diffuses to the core region so as to reduce the gradient in λ, relaxing in the direction of the Taylor state.
Using high resolution bathymetric lidar data for a Telemac2D simulation
NASA Astrophysics Data System (ADS)
Dobler, Wolfgang; Baran, Ramona; Steinbacher, Frank; Ritter, Marcel; Aufleger, Markus
2014-05-01
Knowledge about the hydraulic situation in a mountain torrent is relevant to quantify flood risks, to study sediment transport and to assess the waterbodies' ecology. To conduct reliable calculations, high-quality terrain data of riverbeds, riverbanks and floodplains are required. Typically, digital terrain models (DTMs) of floodplains are derived from classical airborne laserscanning (red wavelength) together with terrestrial surveys along riverbeds and riverbanks. Usually, these are restricted to a limited number of cross sections. Terrestrial surveys are required since laser measurement systems cannot penetrate the water column of the observed waterbodies. Consequently, data describing the geometry of riverbeds and bank structures are hardly available at high spatial resolutions and extents, comparable to the airborne-laser scanning derived data for river floodplains. In this study, a newly available, water-penetrating airborne laser system (green wavelength, FFG research project between the University of Innsbruck and Riegl LMS) was used to survey a mountain torrent. Detailed and extensive data (~30 points/m² on topo-bathy side) of the riverbed and the riverbanks were acquired with this single sensor. In order to construct a 2D-Telemac simulation, the point cloud was down-sampled to an appropriate resolution required for the simulation. The creation of the mesh was carried out with the Software HydroVish and imported into Blue Kenue for further boundary treatment. On one hand the calibration of the numerical model was based on a known water discharge-rate and on the other on abundant data points of the water surface. The green laser system demonstrates its great potential for such an analysis. The final results of the numerical simulation show clearly the supremacy of using such a high resolution data basis in contrast to the traditional way of terrestrial surveying of cross sections along riverbeds.
Lagrangian MHD Particle-in-Cell simulations of coronal interplanetary shocks driven by observations
NASA Astrophysics Data System (ADS)
Lapenta, Giovanni; Bacchini, Fabio; Bemporad, Alessandro; Susino, Roberto; Olshevskyi, Vyacheslav
2016-04-01
In this work, we compare the spatial distribution of the plasma parameters along the June 11, 1999 CME-driven shock front with the results obtained from a CME-like event simulated with the FLIPMHD3D code, based on the FLIP-MHD Particle-in-Cell (PiC) method. The observational data are retrieved from the combination of white-light (WL) coronagraphic data (for the upstream values) and the application of the Rankine-Hugoniot (RH) equations (for the downstream values). The comparison shows a higher compression ratio X and Alfvénic Mach number MA at the shock nose, and a stronger magnetic field deflection d towards the flanks, in agreement with observations. Then, we compare the spatial distribution of MA with the profiles obtained from the solutions of the shock adiabatic equation relating MA, X, and the angle between the upstream magnetic field and the shock front normal for the special cases of parallel and perpendicular shock, and with a semi-empirical expression for a generically oblique shock. The semi-empirical curve approximates the actual values of MA very well, if the effects of a non-negligible shock thickness and plasma-to magnetic pressure ratio are taken into account throughout the computation. Moreover, the simulated shock turns out to be supercritical at the nose and sub-critical at the flanks. Finally, we develop a new 1D Lagrangian ideal MHD method based on the GrAALE code, to simulate the ion-electron temperature decoupling due to the shock transit. Two models are used, a simple solar wind model and a variable-gamma model. Both produce results in agreement with observations, the second one being capable of introducing the physics responsible for the additional electron heating due to secondary effects (collisions, Alfvén waves, etc.). Work supported by the European Commission under the SWIFF project (swiff.eu)
NASA Astrophysics Data System (ADS)
Ge, Y.; Raeder, J.; Angelopoulos, V.; Gilson, M. L.; Runov, A.
2010-12-01
A global MHD simulation has been performed to investigate the THEMIS substorm on February 27, 2009. During this substorm the conjugated observations from the space and on the THEMIS ground observatories are available. The location and time of this substorm onset can be determined based on these observations. The initial auroral brightening is found at around 07:49 UT in the field of view of Fort Smith station (FSMI), with a pre-existing auroral arc located equatorward. A couple minutes later, the in situ observations recorded a sharp dipolarization front sunward passing through THEMIS spacecraft, which travels almost 10 RE in the magnetotail. In this study our global MHD model, i.e., OpenGGCM, driven by the real-time solar wind/IMF conditions, is able to reproduce the key features of these substorm signatures, including the auroral breakup at FSMI with the same onset time as the observations, and a strong earthward Bursty Bulk Flow (BBF) and dipolarization fronts that cause the substorm onset signatures. It is found in the simulation that the auroral breakup is caused by the strong flow shear and the flow vortices which form as the BBF moves earthward. Investigation of the tail BBF and its dipolarization front (DF) reveals that the bipolar change of the Bz component ahead of the DF can be produced by the interaction between two distinct plasmas from separate X lines: the anti-sunward moving southward flux tubes in the tailward flows emanating from an inner magnetic reconnection region, and the sunward traveling dipolarized tubes within the front of a strong earthward BBF that originates in a mid-tail reconnection region. The rebound and oscillations of the intruding BBF reported by the recent THEMIS observations are also seen in the simulation when the BBF encounters the high-pressure inner magnetosphere.
Dayside Proton Aurora: Comparisons between Global MHD Simulations and Image Observations
NASA Technical Reports Server (NTRS)
Berchem, J.; Fuselier, S. A.; Petrinec, S.; Frey, H. U.; Burch, J. L.
2003-01-01
The IMAGE mission provides a unique opportunity to evaluate the accuracy of current global models of the solar wind interaction with the Earth's magnetosphere. In particular, images of proton auroras from the Far Ultraviolet Instrument (FUV) onboard the IMAGE spacecraft are well suited to support investigations of the response of the Earth's magnetosphere to interplanetary disturbances. Accordingly, we have modeled two events that occurred on June 8 and July 28, 2000, using plasma and magnetic field parameters measured upstream of the bow shock as input to three-dimensional magnetohydrodynamic (MHD) simulations. This paper begins with a discussion of images of proton auroras from the FUV SI-12 instrument in comparison with the simulation results. The comparison showed a very good agreement between intensifications in the auroral emissions measured by FUV SI-12 and the enhancement of plasma flows into the dayside ionosphere predicted by the global simulations. Subsequently, the IMAGE observations are analyzed in the context of the dayside magnetosphere's topological changes in magnetic field and plasma flows inferred from the simulation results. Finding include that the global dynamics of the auroral proton precipitation patterns observed by IMAGE are consistent with magnetic field reconnection occurring as a continuous process while the iMF changes in direction and the solar wind dynamic pressure varies. The global simulations also indicate that some of the transient patterns observed by IMAGE are consistent with sporadic reconnection processes. Global merging patterns found in the simulations agree with the antiparallel merging model. though locally component merging might broaden the merging region, especially in the region where shocked solar wind discontinuities first reach the magnetopause. Finally, the simulations predict the accretion of plasma near the bow shock in the regions threaded by newly open field lines on which plasma flows into the dayside
NASA Technical Reports Server (NTRS)
Kuznetsova, M. M.; Sibeck, D. G.; Hesse, M.; Wang, Y.; Rastaetter, L.; Toth, G.; Ridley, A.
2009-01-01
We use the global magnetohydrodynamic (MHD) code BATS-R-US to model multipoint observations of Flux Transfer Event (FTE) signatures. Simulations with high spatial and temporal resolution predict that cavities of weak magnetic field strength protruding into the magnetosphere trail FTEs. These predictions are consistent with recently reported multi-point Cluster observations of traveling magnetopause erosion regions (TMERs).