WorldWideScience

Sample records for 3-d global mhd

  1. Trapping Solids at the Inner Edge of the Dead Zone: 3-D Global MHD Simulations

    Dzyurkevich, Natalia; Flock, Mario; Turner, Neal J.; Klahr, Hubert; Henning, Thomas

    2010-01-01

    The poorly-ionized interior of the protoplanetary disk is the location where dust coagulation processes may be most efficient. However even here, planetesimal formation may be limited by the loss of solid material through radial drift, and by collisional fragmentation of the particles. Our aim is to investigate the possibility that solid particles are trapped at local pressure maxima in the dynamically evolving disk. We perform the first 3-D global non-ideal MHD calculations of the disk treat...

  2. Capabilities of a Global 3D MHD Model for Monitoring Extremely Fast CMEs

    Wu, C. C.; Plunkett, S. P.; Liou, K.; Socker, D. G.; Wu, S. T.; Wang, Y. M.

    2015-12-01

    Since the start of the space era, spacecraft have recorded many extremely fast coronal mass ejections (CMEs) which have resulted in severe geomagnetic storms. Accurate and timely forecasting of the space weather effects of these events is important for protecting expensive space assets and astronauts and avoiding communications interruptions. Here, we will introduce a newly developed global, three-dimensional (3D) magnetohydrodynamic (MHD) model (G3DMHD). The model takes the solar magnetic field maps at 2.5 solar radii (Rs) and intepolates the solar wind plasma and field out to 18 Rs using the algorithm of Wang and Sheeley (1990, JGR). The output is used as the inner boundary condition for a 3D MHD model. The G3DMHD model is capable of simulating (i) extremely fast CME events with propagation speeds faster than 2500 km/s; and (ii) multiple CME events in sequence or simultaneously. We will demonstrate the simulation results (and comparison with in-situ observation) for the fastest CME in record on 23 July 2012, the shortest transit time in March 1976, and the well-known historic Carrington 1859 event.

  3. Turbulence and Steady Flows in 3D Global Stratified MHD Simulations of Accretion Disks

    Flock, M; Klahr, H; Turner, N J; Henning, Th

    2011-01-01

    We present full 2 Pi global 3-D stratified MHD simulations of accretion disks. We interpret our results in the context of proto-planetary disks. We investigate the turbulence driven by the magneto-rotational instability (MRI) using the PLUTO Godunov code in spherical coordinates with the accurate and robust HLLD Riemann solver. We follow the turbulence for more than 1500 orbits at the innermost radius of the domain to measure the overall strength of turbulent motions and the detailed accretion flow pattern. We find that regions within two scale heights of the midplane have a turbulent Mach number of about 0.1 and a magnetic pressure two to three orders of magnitude less than the gas pressure, while outside three scale heights the magnetic pressure equals or exceeds the gas pressure and the turbulence is transonic, leading to large density fluctuations. The strongest large-scale density disturbances are spiral density waves, and the strongest of these waves has m=5. No clear meridional circulation appears in t...

  4. Trapping Solids at the Inner Edge of the Dead Zone: 3-D Global MHD Simulations

    Dzyurkevich, Natalia; Turner, Neal J; Klahr, Hubert; Henning, Thomas

    2010-01-01

    The poorly-ionized interior of the protoplanetary disk is the location where dust coagulation processes may be most efficient. However even here, planetesimal formation may be limited by the loss of solid material through radial drift, and by collisional fragmentation of the particles. Our aim is to investigate the possibility that solid particles are trapped at local pressure maxima in the dynamically evolving disk. We perform the first 3-D global non-ideal MHD calculations of the disk treating the turbulence driven by the magneto-rotational instability. The domain contains an inner MRI-active region near the young star and an outer midplane dead zone, with the transition between the two modeled by a sharp increase in the magnetic diffusivity. The azimuthal magnetic fields generated in the active zone oscillate over time, changing sign about every 150 years. We thus observe the radial structure of the `butterfly pattern' seen previously in local shearing-box simulations. The mean magnetic field diffuses from...

  5. Trapping solids at the inner edge of the dead zone: 3-D global MHD simulations

    Dzyurkevich, N.; Flock, M.; Turner, N. J.; Klahr, H.; Henning, Th.

    2010-06-01

    Context. The poorly-ionized interior of the protoplanetary disk or “dead zone” is the location where dust coagulation processes may be most efficient. However even here, planetesimal formation may be limited by the loss of solid material through radial drift, and by collisional fragmentation of the particles. Both depend on the turbulent properties of the gas. Aims: Our aim here is to investigate the possibility that solid particles are trapped at local pressure maxima in the dynamically evolving disk. We perform the first 3-D global non-ideal magnetohydrodynamical (MHD) calculations of a section of the disk treating the turbulence driven by the magneto-rotational instability (MRI). Methods: We use the ZeusMP code with a fixed Ohmic resistivity distribution. The domain contains an inner MRI-active region near the young star and an outer midplane dead zone, with the transition between the two modeled by a sharp increase in the magnetic diffusivity. Results: The azimuthal magnetic fields generated in the active zone oscillate over time, changing sign about every 150 years. We thus observe the radial structure of the “butterfly pattern” seen previously in local shearing-box simulations. The mean magnetic field diffuses from the active zone into the dead zone, where the Reynolds stress nevertheless dominates, giving a residual α between 10-4 and 10-3. The greater total accretion stress in the active zone leads to a net reduction in the surface density, so that after 800 years an approximate steady state is reached in which a local radial maximum in the midplane pressure lies near the transition radius. We also observe the formation of density ridges within the active zone. Conclusions: The dead zone in our models possesses a mean magnetic field, significant Reynolds stresses and a steady local pressure maximum at the inner edge, where the outward migration of planetary embryos and the efficient trapping of solid material are possible.

  6. 3D Dynamics of Magnetopause Reconnection Using Hall-MHD Global Simulations

    Maynard, K.; Germaschewski, K.; Raeder, J.; Bhattacharjee, A.

    2011-12-01

    Magnetic reconnection at Earth's magnetopause and in the magnetotail is of crucial importance for the dynamics of the global magnetosphere and space weather. Even though the plasma conditions in the magnetosphere are largely in the collisionless regime, most of the existing research using global computational models employ single-fluid magnetohydrodynamics (MHD) with artificial resistivity. Studies of reconnection in simplified, two-dimensional geometries have established that two-fluid and kinetic effects can dramatically alter dynamics and reconnection rates when compared with single-fluid models. These enhanced models also introduce particular signatures, for example a quadrupolar out-of-plane magnetic field component that has already been observed in space by satellite measurements. However, results from simplified geometries cannot be translated directly to the dynamics of three-dimensional magnetospheric reconnection. For instance, magnetic flux originating from the solar wind and arriving at the magnetopause can either reconnect or be advected around the magnetosphere. In this study, we use a new version of the OpenGGCM code that incorporates the Hall term in a Generalized Ohm's Law to study magnetopause reconnection under synthetic solar wind conditions and investigate how reconnection rates and dynamics of flux transfer events depend on the strength of the Hall term. The OpenGGCM, a global model of Earth's magnetosphere, has recently been ported to exploit modern computing architectures like the Cell processor and SIMD capabilities of conventional processors using an automatic code generator. These enhancements provide us with the performance needed to include the computationally expensive Hall physics.

  7. 3-D Relativistic MHD Simulations

    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.

  8. A high-order Godunov scheme for global 3D MHD accretion disks simulations. I. The linear growth regime of the magneto-rotational instability

    Flock, M; Klahr, H; Mignone, A

    2009-01-01

    We employ the PLUTO code for computational astrophysics to assess and compare the validity of different numerical algorithms on simulations of the magneto-rotational instability in 3D accretion disks. In particular we stress on the importance of using a consistent upwind reconstruction of the electro-motive force (EMF) when using the constrained transport (CT) method to avoid the onset of numerical instabilities. We show that the electro-motive force (EMF) reconstruction in the classical constrained transport (CT) method for Godunov schemes drives a numerical instability. The well-studied linear growth of magneto-rotational instability (MRI) is used as a benchmark for an inter-code comparison of PLUTO and ZeusMP. We reproduce the analytical results for linear MRI growth in 3D global MHD simulations and present a robust and accurate Godunov code which can be used for 3D accretion disk simulations in curvilinear coordinate systems.

  9. High-order Godunov schemes for global 3D MHD simulations of accretion disks. I. Testing the linear growth of the magneto-rotational instability

    Flock, M.; Dzyurkevich, N.; Klahr, H.; Mignone, A.

    2010-06-01

    We assess the suitability of various numerical MHD algorithms for astrophysical accretion disk simulations with the PLUTO code. The well-studied linear growth of the magneto-rotational instability is used as the benchmark test for a comparison between the implementations within PLUTO and against the ZeusMP code. The results demonstrate the importance of using an upwind reconstruction of the electro-motive force (EMF) in the context of a constrained transport scheme, which is consistent with plane-parallel, grid-aligned flows. In contrast, constructing the EMF from the simple average of the Godunov fluxes leads to a numerical instability and the unphysical growth of the magnetic energy. We compare the results from 3D global calculations using different MHD methods against the analytical solution for the linear growth of the MRI, and discuss the effect of numerical dissipation. The comparison identifies a robust and accurate code configuration that is vital for realistic modeling of accretion disk processes.

  10. 3D simulation studies of tokamak plasmas using MHD and extended-MHD models

    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-β disruption studies in reversed shear plasmas using the MHD level MH3D code, ω*i stabilization and nonlinear island saturation of TAE mode using the hybrid particle/MHD level MH3D-K code, and unstructured mesh MH3D++ code studies. In particular, three internal mode disruption mechanisms are identified from simulation results which agree which agree well with experimental data

  11. 3D simulation studies of tokamak plasmas using MHD and extended-MHD models

    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-β disruption studies in reversed shear plasmas using the MHD level MH3D code, ω*i stabilization and nonlinear island rotation studies using the two-fluid level MH3D-T code, studies of nonlinear saturation of TAE modes using the hybrid particle/MHD level MH3D-K code, and unstructured mesh MH3D++ code studies. In particular, three internal mode disruption mechanisms are identified from simulation results which agree well with experimental data

  12. 3D simulation studies of tokamak plasmas using MHD and extended-MHD models

    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-β disruption studies in reversed shear plasmas using the MHD level MH3D code; ω*i stabilization and nonlinear island rotation studies using the two-fluid level MH3D-T code; studies of nonlinear saturation of TAE modes using the hybrid particle/MHD level MH3D-K code; and unstructured mesh MH3D++ code studies. In particular, three internal mode disruption mechanisms are identified from simulation results which agree well with experimental data. (author). 18 refs, 5 figs

  13. 3D MHD Flux emergence experiments

    Hood, A.W.; Archontis, V.; Mactaggart, David

    2012-01-01

    This paper reviews some of the many 3D numerical experiments of the emergence of magnetic fields from the solar interior and the subsequent interaction with the pre-existing coronal magnetic field. The models described here are idealised, in the sense that the internal energy equation only involv...

  14. Suppression of Magnetic Flux Diffusion in Reduced 3D MHD

    Bayliss, A.; Ware, A. S.; Diamond, P. H.; Kim, E.-J.

    1999-11-01

    The important impact of small scale magnetic fields on self-organization (i.e., dynamo) in MHD turbulence was originally pin-pointed by the observation that magnetic flux diffusion (anomalous resistivity) is drastically reduced in 2D MHD turbulence. This reduction is a consequence of mean square magnetic potential in two dimensions. It is natural, then, to investigate magnetic flux diffusion in 3D reduced MHD; since in that system conservation is broken only by linear field line bending (symptomatic of Alfvén wave propagation along B_z), and resistive dissipation. In particular, the Ohm's Law nonlinearity conserves . Not surprisingly, it is possible to derive an exact constraint upon the spatial flux of magnetic potential from the condition of balance. This expression may then be used to simplify the calculation of the turbulent resistivity, which is found to be suppressed, as in 2D MHD, up to corrections resulting from hat z-direction Alfvén wave propagation effects. These corrections vanish in the limit of unity magnetic Prandtl number. Work on understanding the self-consistent alpha effect in reduced MHD is ongoing and will be discussed.

  15. PIXIE3D: An efficient, fully implicit, parallel, 3D extended MHD code for fusion plasma modeling

    PIXIE3D is a modern, parallel, state-of-the-art extended MHD code that employs fully implicit methods for efficiency and accuracy. It features a general geometry formulation, and is therefore suitable for the study of many magnetic fusion configurations of interest. PIXIE3D advances the state of the art in extended MHD modeling in two fundamental ways. Firstly, it employs a novel conservative finite volume scheme which is remarkably robust and stable, and demands very small physical and/or numerical dissipation. This is a fundamental requirement when one wants to study fusion plasmas with realistic conductivities. Secondly, PIXIE3D features fully-implicit time stepping, employing Newton-Krylov methods for inverting the associated nonlinear systems. These methods have been shown to be scalable and efficient when preconditioned properly. Novel preconditioned ideas (so-called physics based), which were prototypes in the context of reduced MHD, have been adapted for 3D primitive-variable resistive MHD in PIXIE3D, and are currently being extended to Hall MHD. PIXIE3D is fully parallel, employing PETSc for parallelism. PIXIE3D has been thoroughly benchmarked against linear theory and against other available extended MHD codes on nonlinear test problems (such as the GEM reconnection challenge). We are currently in the process of extending such comparisons to fusion-relevant problems in realistic geometries. In this talk, we will describe both the spatial discretization approach and the preconditioning strategy employed for extended MHD in PIXIE3D. We will report on recent benchmarking studies between PIXIE3D and other 3D extended MHD codes, and will demonstrate its usefulness in a variety of fusion-relevant configurations such as Tokamaks and Reversed Field Pinches. (Author)

  16. Testing Observational Techniques with 3D MHD Jets in Clusters

    Mendygral, Peter J; Jones, Tom W

    2009-01-01

    Observations of X-ray cavities formed by powerful jets from AGN in galaxy cluster cores are commonly used to estimate the mechanical luminosity of these sources. We test the reliability of observationally measuring this power with synthetic X-ray observations of 3-D MHD simulations of jets in a galaxy cluster environment. We address the role that factors such as jet intermittency and orientation of the jets on the sky have on the reliability of observational measurements of cavity enthalpy and age. An estimate of the errors in these quantities can be made by directly comparing ``observationally'' derived values with values from the simulations. In our tests, cavity enthalpy, age and mechanical luminosity derived from observations are within a factor of two of the simulation values.

  17. 3-D Relativistic MHD Simulations of Extragalactic Jets

    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.

  18. 3D MHD simulation of polarized emission in SN 1006

    Schneiter, E M; Reynoso, E M; Esquivel, A; De Colle, F

    2015-01-01

    We use three dimensional magnetohydrodynamic (MHD) 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.

  19. 3D MHD simulations of pellet injection and disruptions in tokamak plasmas

    Nonlinear MHD simulation results of pellet injection show that MHD forces can accelerate large pellets, injected on the high field side of a tokamak, to the plasma center. Magnetic reconnection can produce a reverse shear q profile. Ballooning instability caused by pellets is also reduced by high field side injection. Studies are also reported of the current quench phase of disruptions, which can cause 3D halo currents and runaway electrons. (author)

  20. FARGO3D: A new GPU-oriented MHD code

    Benítez-Llambay, Pablo

    2016-01-01

    We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on protoplanetary disks physics and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve the divergence-free property of the magnetic field to machine accuracy, coupled to a method of characteristics for the evaluation of electromotive forces and Lorentz forces. Orbital advection is implemented, and an N-body solver is included to simulate planets or stars interacting with the gas. We present our implementation in detail and present a number of widely known tests for comparison purposes. One strength of FARGO3D is that it can run on both "Graphical Processing Units" (GPUs) or "Central Processing unit" (CPUs), achieving large speed up with respect to CPU cores. We describe our implementation choices, whi...

  1. 3-D MHD Numerical Simulations of Cloud-Wind Interactions

    Gregori, G.; Miniati, Francesco; Ryu, Dongsu; Jones, T. W.

    2000-01-01

    We present results from three-dimensional (3-D) numerical simulations investigating the magnetohydrodynamics of cloud-wind interactions. The initial cloud is spherical while the magnetic field is uniform and transverse to the cloud motion. A simplified analytical model that describes the magnetic energy evolution in front of the cloud is developed and compared with simulation results. In addition, it is found the interaction of the cloud with a magnetized interstellar medium (ISM) results in ...

  2. Study of magnetic island using a 3D MHD equilibrium calculation code

    Coupling the magnetic diagnostics and a 3D MHD equilibrium calculation code, the magnetic island is studied in the Large Helical Device (LHD) experiment. In an experiment, the collapse in the plasma core was observed in a configuration, which has large magnetic island produced by external perturbation coils. At the collapse, the temperature profile was flattened. This suggests the magnetic island evolved. The magnetic island was observed by the magnetic diagnostics. The magnetic diagnostics also suggests evolving the magnetic island. A 3D MHD equilibrium is calculated by the 3D MHD equilibrium code then signals of the magnetic diagnostics are simulated. Since the comparison of observed and calculated signals is comparable, the magnetic island in calculated equilibrium is similar to one of the experiment. (author)

  3. 3D MHD modeling of twisted coronal loops

    Reale, F; Guarrasi, M; Mignone, A; Peres, G; Hood, A W; Priest, E R

    2016-01-01

    We perform MHD modeling of a single bright coronal loop to include the interaction with a non-uniform magnetic field. The field is stressed by random footpoint rotation in the central region and its energy is dissipated into heating by growing currents through anomalous magnetic diffusivity that switches on in the corona above a current density threshold. We model an entire single magnetic flux tube, in the solar atmosphere extending from the high-beta chromosphere to the low-beta corona through the steep transition region. The magnetic field expands from the chromosphere to the corona. The maximum resolution is ~30 km. We obtain an overall evolution typical of loop models and realistic loop emission in the EUV and X-ray bands. The plasma confined in the flux tube is heated to active region temperatures (~3 MK) after ~2/3 hr. Upflows from the chromosphere up to ~100 km/s fill the core of the flux tube to densities above 10^9 cm^-3. More heating is released in the low corona than the high corona and is finely ...

  4. FARGO3D: A New GPU-oriented MHD Code

    Benítez-Llambay, Pablo; Masset, Frédéric S.

    2016-03-01

    We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on the physics of protoplanetary disks and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite-difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve the divergence-free property of the magnetic field to machine accuracy, coupled to a method of characteristics for the evaluation of electromotive forces and Lorentz forces. Orbital advection is implemented, and an N-body solver is included to simulate planets or stars interacting with the gas. We present our implementation in detail and present a number of widely known tests for comparison purposes. One strength of FARGO3D is that it can run on either graphical processing units (GPUs) or central processing units (CPUs), achieving large speed-up with respect to CPU cores. We describe our implementation choices, which allow a user with no prior knowledge of GPU programming to develop new routines for CPUs, and have them translated automatically for GPUs.

  5. 3-D MHD Numerical Simulations of Cloud-Wind Interactions

    Gregori, G; Ryu, D; Jones, T W; Miniati, Francesco; Ryu, Dongsu

    2000-01-01

    We present results from three-dimensional (3-D) numerical simulations investigating the magnetohydrodynamics of cloud-wind interactions. The initial cloud is spherical while the magnetic field is uniform and transverse to the cloud motion. A simplified analytical model that describes the magnetic energy evolution in front of the cloud is developed and compared with simulation results. In addition, it is found the interaction of the cloud with a magnetized interstellar medium (ISM) results in the formation of a highly structured magnetotail. The magnetic flux in the wake of the cloud organizes into flux ropes and a reconnection, current sheet is developed, as field lines of opposite polarity are brought close together near the symmetry axis. At the same time, magnetic pressure is strongly enhanced at the leading edge of the cloud from the stretching of the field lines that occurs there. This has an important dynamical effect on the subsequent evolution of the cloud, since some unstable modes tend to be strongl...

  6. 3D-MHD simulations of an accretion disk with star-disk boundary layer

    Steinacker, A; Steinacker, Adriane; Papaloizou, John C.B.

    2002-01-01

    We present global 3D MHD simulations of geometrically thin but unstratified accretion disks in which a near Keplerian disk rotates between two bounding regions with initial rotation profiles that are stable to the MRI. The inner region models the boundary layer between the disk and an assumed more slowly rotating central, non magnetic star. We investigate the dynamical evolution of this system in response to initial vertical and toroidal fields imposed in a variety of domains contained within the near Keplerian disk. Cases with both non zero and zero net magnetic flux are considered and sustained dynamo activity found in runs for up to fifty orbital periods at the outer boundary of the near Keplerian disk. Simulations starting from fields with small radial scale and with zero net flux lead to the lowest levels of turbulence and smoothest variation of disk mean state variables. For our computational set up, average values of the Shakura & Sunyaev (1973) $\\alpha$ parameter in the Keplerian disk are typicall...

  7. 3D Tomography of MHD Fluctuations in the H-1NF Heliac

    Haskey, S R; Seiwald, B; Howard, J

    2014-01-01

    A 3D tomographic reconstruction technique which does not rely on a set of radial basis functions is described for inversion of a set of limited-angle high-resolution 2D visible light emission projections (extended in the vertical and toroidal directions) of global MHD eigenmodes in the H-1NF heliac. This paper deals with some of the features and challenges that will arise in the application of tomographic imaging systems to fusion reactors, especially the strong shaping of optimised stellarator/heliotron configurations, and limited access in all types. The fluctuations are represented as a finite sum of Fourier modes characterised by toroidal and poloidal mode numbers having fixed amplitude and phase in a set of nested cylindrical flux volumes in Boozer space. The amplitudes and phases are calculated using iterative tomographic inversion techniques such as ART, SIRT and standard linear least-squares methods. The tomography is applied to synchronous camera images of singly charged carbon impurity ion emission ...

  8. 3D hybrid and MHD/particle simulations of field-reversed configurations

    A nonlinear 3D code in cylindrical geometry is being developed for the stability studies of FRC. Two numerical schemes have been implemented: a hybrid scheme with particle ions and fluid electrons, and MHD/particle scheme in which the background plasma is described by MHD equations. And energetic ions are treated via particle simulations. The MHD equations are advanced on a finite-difference mesh in a cylindrical coordinate system, while particle pushing is done on 3D Cartesian grids. Full ion dynamics is retained in order to include large-orbit effects (with s∼1), which are important for the tilt mode stabilization in FRC. Also, in contrast to the previous work, δf method is utilized to reduce numerical noise in the simulations. The code has been benchmarked against previous MHD simulation of tilting instability in FRC. It was found that rigid rotation reduces the growth rate, but does not stabilize the mode even for rotation rates equal to the Alfven time. Sheared rotation is found to be destabilizing for the velocity profile considered. Simulations with a fast ion beam with 1 % of the bulk ion density and s∼3 did not show a reduction in growth rate of the tilting instability. (author)

  9. Implementation of a 3-D nonlinear MHD [magnetohydrodynamics] calculation on the Intel hypercube

    The optimization of numerical schemes and increasing computer capabilities in the last ten years have improved the efficiency of 3-D nonlinear resistive MHD calculations by about two to three orders of magnitude. However, we are still very limited in performing these types of calculations. Hypercubes have a large number of processors with only local memory and bidirectional links among neighbors. The Intel Hypercube at Oak Ridge has 64 processors with 0.5 megabytes of memory per processor. The multiplicity of processors opens new possibilities for the treatment of such computations. The constraint on time and resources favored the approach of using the existing RSF code which solves as an initial value problem the reduced set of MHD equations for a periodic cylindrical geometry. This code includes minimal physics and geometry, but contains the basic three dimensionality and nonlinear structure of the equations. The code solves the reduced set of MHD equations by Fourier expansion in two angular coordinates and finite differences in the radial one. Due to the continuing interest in these calculations and the likelihood that future supercomputers will take greater advantage of parallelism, the present study was initiated by the ORNL Exploratory Studies Committee and funded entirely by Laboratory Discretionary Funds. The objectives of the study were: to ascertain the suitability of MHD calculation for parallel computation, to design and implement a parallel algorithm to perform the computations, and to evaluate the hypercube, and in particular, ORNL's Intel iPSC, for use in MHD computations

  10. 3D MHD Jet in a Non-Uniform Magnetic Field

    Huang Hulin; Han Dong

    2005-01-01

    The purpose of this paper is to present a two-phase 3D magnetohydrodynamics (MHD) flow model that combines the volume of fluid (VOF) method with the technique derived from induced-magnetic-field equations for liquid metal free surface MHD-jet-flow. Analogy between the induced-magnetic-filed equation and the conventional computational fluid dynamics (CFD) equation is made, so that the equation can be conveniently accounted for by CFD. A penalty factor numerical method is introduced in order to force the local divergence-free condition of the magnetic fields and an extension of the void insulating calculation domain is applied to ensure that the induced-magnetic field at its boundaries is null. These simulation results for lithium liquid metal jets under magnetic field configurations of Magnetic Torus (Mtor) and National Spherical Torus Experiment (NSTX) outboard divertor have shown that three dimensional jet can not be annihilated by magnetic braking and its cross-section will deform in such a way that the momentum flux of the jet is conserved. 3D MHD effects from a magnetic field gradient cause return currents to interact with applied magnetic fields and produce unfavorable Lorentz forces.Under 3D applied non-uniform magnetic fields of the divertor, unfavorable Lorentz forces lead to a substantial change in flow pattern and a reduction in flow velocity, with the jet cross-section moving to one side of the jet space. These critical phenomena can not be revealed by 2D models.

  11. 3D MHD free surface fluid flow simulation based on magnetic-field induction equations

    The purpose of this paper is to present our recent efforts on 3D MHD model development and our results based on the technique derived from induced-magnetic-field equations. Two important features are utilized in our numerical method to obtain convergent solutions. First, a penalty factor is introduced in order to force the local divergence free condition of the magnetic fields. The second is that we extend the insulating wall thickness to ensure that the induced magnetic field at its boundaries is null. These simulation results for lithium film free surface flows under NSTX outboard mid-plane magnetic field configurations have shown that 3D MHD effects from a surface normal field gradient cause return currents to interact with surface normal fields and produce unfavorable MHD forces. This leads to a substantial change in flow pattern and a reduction in flow velocity, with most of the flow spilling over one side of the chute. These critical phenomena can not be revealed by 2D models. Additionally, a design which overcomes these undesired flow characteristics is obtained

  12. Roles of initial current carrier in the distribution of field-aligned current in 3-D Hall MHD simulations

    ZHANG XianGuo; PU ZuYin; MA ZhiWei; ZHOU XuZhi

    2008-01-01

    A three-dimensional (3-D) Hall MHD simulation is carried out to study the roles of initial current carrier in the topology of magnetic field,the generation and distribuering the contribution of ions to the initial current,the topology of the obtained magnetic field turns to be more complex. In some cases,it is found that not only the traditional By quadrupole structure but also a reversal By quadrupole structure appears in the simulation box. This can explain the observational features near the diffusion region,which are inconsistent with the Hall MHD theory with the total initial current carried by electrons. Several other interesting features are also emerged. First,motions of electrons and ions are decoupled from each other in the small plasma region (Hall effect region) with a scale less than or comparable with the ion inertial length or ion skin depth di=c/ωp. In the non-Hall effect region,the global magnetic structure is shifted in +y direction under the influence of ions with initial y directional motion. However,in the Hall effect region,magnetic field lines are bent in -y direction,mainly controlled by the motion of electrons,then By is generated. Second,FACs emerge as a result of the appearance of By. Compared with the prior Hall MHD simulation results,the generated FACs shift in +y direction,

  13. 3D Simulations of MHD Jet Propagation Through Uniform and Stratified External Environments

    O'Neill, S. M.; Tregillis, I. L.; Jones, T. W.; Ryu, Dongsu

    2005-01-01

    We present a set of high-resolution 3D MHD simulations of steady light, supersonic jets, exploring the influence of jet Mach number and the ambient medium on jet propagation and energy deposition over long distances. The results are compared to simple self-similar scaling relations for the morphological evolution of jet-driven structures and to previously published 2D simulations. For this study we simulated the propagation of light jets with internal Mach numbers 3 and 12 to lengths exceedin...

  14. Nanoflare statistics in an active region 3D MHD coronal model

    Bingert, Sven

    2012-01-01

    Context. We investigate the statistics of the spatial and temporal distribution of the coronal heating in a three-dimensional magneto- hydrodynamical (3D MHD) model. The model describes the temporal evolution of the corona above an observed active region. The model is driven by photospheric granular motions which braid the magnetic field lines. This induces currents and their dissipation heats the plasma. We evaluate the transient heating as subsequent heating events and analyze their statistics. The results are then interpreted in the context of observed flare statistics and coronal heating mechanisms. Methods. To conduct the numerical experiment we use a high order finite difference code which solves the partial differential equations for the conservation of mass, the momentum and energy balance, and the induction equation. The energy balance includes the Spitzer heat conduction and the optical thin radiative loss in the corona. Results. The temporal and spatial distribution of the Ohmic heating in the 3D M...

  15. Effect of magnetic perturbations on the 3D MHD self-organization of shaped tokamak plasmas

    Bonfiglio, D; Veranda, M; Chacón, L; Escande, D F

    2016-01-01

    The effect of magnetic perturbations (MPs) on the helical self-organization of shaped tokamak plasmas is discussed in the framework of the nonlinear 3D MHD model. Numerical simulations performed in toroidal geometry with the \\textsc{pixie3d} code [L. Chac\\'on, Phys. Plasmas {\\bf 15}, 056103 (2008)] show that $n=1$ MPs significantly affect the spontaneous quasi-periodic sawtoothing activity of such plasmas. In particular, the mitigation of sawtooth oscillations is induced by $m/n=1/1$ and $2/1$ MPs. These numerical findings provide a confirmation of previous circular tokamak simulations, and are in agreement with tokamak experiments in the RFX-mod and DIII-D devices. Sawtooth mitigation via MPs has also been observed in reversed-field pinch simulations and experiments. The effect of MPs on the stochastization of the edge magnetic field is also discussed.

  16. 3D Neutronic Analysis in MHD Calculations at ARIES-ST Fusion Reactors Systems

    Hançerliogulları, Aybaba; Cini, Mesut

    2013-10-01

    In this study, we developed new models for liquid wall (FW) state at ARIES-ST fusion reactor systems. ARIES-ST is a 1,000 MWe fusion reactor system based on a low aspect ratio ST plasma. In this article, we analyzed the characteristic properties of magnetohydrodynamics (MHD) and heat transfer conditions by using Monte-Carlo simulation methods (ARIES Team et al. in Fusion Eng Des 49-50:689-695, 2000; Tillack et al. in Fusion Eng Des 65:215-261, 2003) . In fusion applications, liquid metals are traditionally considered to be the best working fluids. The working liquid must be a lithium-containing medium in order to provide adequate tritium that the plasma is self-sustained and that the fusion is a renewable energy source. As for Flibe free surface flows, the MHD effects caused by interaction with the mean flow is negligible, while a fairly uniform flow of thick can be maintained throughout the reactor based on 3-D MHD calculations. In this study, neutronic parameters, that is to say, energy multiplication factor radiation, heat flux and fissile fuel breeding were researched for fusion reactor with various thorium and uranium molten salts. Sufficient tritium amount is needed for the reactor to work itself. In the tritium breeding ratio (TBR) >1.05 ARIES-ST fusion model TBR is >1.1 so that tritium self-sufficiency is maintained for DT fusion systems (Starke et al. in Fusion Energ Des 84:1794-1798, 2009; Najmabadi et al. in Fusion Energ Des 80:3-23, 2006).

  17. A novel code for numerical 3-D MHD studies of CME expansion

    J. Kleimann

    2009-03-01

    Full Text Available A recent third-order, essentially non-oscillatory central scheme to advance the equations of single-fluid magnetohydrodynamics (MHD in time has been implemented into a new numerical code. This code operates on a 3-D Cartesian, non-staggered grid, and is able to handle shock-like gradients without producing spurious oscillations.

    To demonstrate the suitability of our code for the simulation of coronal mass ejections (CMEs and similar heliospheric transients, we present selected results from test cases and perform studies of the solar wind expansion during phases of minimum solar activity. We can demonstrate convergence of the system into a stable Parker-like steady state for both hydrodynamic and MHD winds. The model is subsequently applied to expansion studies of CME-like plasma bubbles, and their evolution is monitored until a stationary state similar to the initial one is achieved. In spite of the model's (current simplicity, we can confirm the CME's nearly self-similar evolution close to the Sun, thus highlighting the importance of detailed modelling especially at small heliospheric radii.

    Additionally, alternative methods to implement boundary conditions at the coronal base, as well as strategies to ensure a solenoidal magnetic field, are discussed and evaluated.

  18. Radiative transfer with scattering for domain-decomposed 3D MHD simulations of cool stellar atmospheres

    Hayek, W; Carlsson, M; Trampedach, R; Collet, R; Gudiksen, B V; Hansteen, V H; Leenaarts, J

    2010-01-01

    We present the implementation of a radiative transfer solver with coherent scattering in the new BIFROST code for radiative magneto-hydrodynamical (MHD) simulations of stellar surface convection. The code is fully parallelized using MPI domain decomposition, which allows for large grid sizes and improved resolution of hydrodynamical structures. We apply the code to simulate the surface granulation in a solar-type star, ignoring magnetic fields, and investigate the importance of coherent scattering for the atmospheric structure. A scattering term is added to the radiative transfer equation, requiring an iterative computation of the radiation field. We use a short-characteristics-based Gauss-Seidel acceleration scheme to compute radiative flux divergences for the energy equation. The effects of coherent scattering are tested by comparing the temperature stratification of three 3D time-dependent hydrodynamical atmosphere models of a solar-type star: without scattering, with continuum scattering only, and with bo...

  19. Quasi 3D ECE imaging system for study of MHD instabilities in KSTAR.

    Yun, G S; Lee, W; Choi, M J; Lee, J; Kim, M; Leem, J; Nam, Y; Choe, G H; Park, H K; Park, H; Woo, D S; Kim, K W; Domier, C W; Luhmann, N C; Ito, N; Mase, A; Lee, S G

    2014-11-01

    A second electron cyclotron emission imaging (ECEI) system has been installed on the KSTAR tokamak, toroidally separated by 1/16th of the torus from the first ECEI system. For the first time, the dynamical evolutions of MHD instabilities from the plasma core to the edge have been visualized in quasi-3D for a wide range of the KSTAR operation (B0 = 1.7∼3.5 T). This flexible diagnostic capability has been realized by substantial improvements in large-aperture quasi-optical microwave components including the development of broad-band polarization rotators for imaging of the fundamental ordinary ECE as well as the usual 2nd harmonic extraordinary ECE. PMID:25430233

  20. Quasi 3D ECE imaging system for study of MHD instabilities in KSTAR

    Yun, G. S., E-mail: gunsu@postech.ac.kr; Choi, M. J.; Lee, J.; Kim, M.; Leem, J.; Nam, Y.; Choe, G. H. [Department of Physics, Pohang University of Science and Technology, Pohang 790-784 (Korea, Republic of); Lee, W.; Park, H. K. [Ulsan National Institute of Science and Technology, Ulsan 689-798 (Korea, Republic of); Park, H.; Woo, D. S.; Kim, K. W. [School of Electrical Engineering, Kyungpook National University, Daegu 702-701 (Korea, Republic of); Domier, C. W.; Luhmann, N. C. [Department of Electrical and Computer Engineering, University of California, Davis, California 95616 (United States); Ito, N. [KASTEC, Kyushu University, Kasuga-shi, Fukuoka 812-8581 (Japan); Mase, A. [Ube National College of Technology, Ube-shi, Yamaguchi 755-8555 (Japan); Lee, S. G. [National Fusion Research Institute, Daejeon 305-333 (Korea, Republic of)

    2014-11-15

    A second electron cyclotron emission imaging (ECEI) system has been installed on the KSTAR tokamak, toroidally separated by 1/16th of the torus from the first ECEI system. For the first time, the dynamical evolutions of MHD instabilities from the plasma core to the edge have been visualized in quasi-3D for a wide range of the KSTAR operation (B{sub 0} = 1.7∼3.5 T). This flexible diagnostic capability has been realized by substantial improvements in large-aperture quasi-optical microwave components including the development of broad-band polarization rotators for imaging of the fundamental ordinary ECE as well as the usual 2nd harmonic extraordinary ECE.

  1. Roles of initial current carrier in the distribution of field-aligned current in 3-D Hall MHD simulations

    2008-01-01

    A three-dimensional (3-D) Hall MHD simulation is carried out to study the roles of initial current carrier in the topology of magnetic field, the generation and distribu- tion of field aligned currents (FACs), and the appearance of Alfvén waves. Consid- ering the contribution of ions to the initial current, the topology of the obtained magnetic field turns to be more complex. In some cases, it is found that not only the traditional By quadrupole structure but also a reversal By quadrupole structure appears in the simulation box. This can explain the observational features near the diffusion region, which are inconsistent with the Hall MHD theory with the total ini- tial current carried by electrons. Several other interesting features are also emerged. First, motions of electrons and ions are decoupled from each other in the small plasma region (Hall effect region) with a scale less than or comparable with the ion inertial length or ion skin depth di=c/ωp. In the non-Hall effect region, the global magnetic structure is shifted in +y direction under the influence of ions with initial y directional motion. However, in the Hall effect region, magnetic field lines are bent in ?y direction, mainly controlled by the motion of electrons, then By is generated. Second, FACs emerge as a result of the appearance of By. Compared with the prior Hall MHD simulation results, the generated FACs shift in +y direction, and hence the dawn-dusk symmetry is broken. Third, the Walén relation in our simulations is consistent with the Walén relation in Hall plasma, thus the presence of Alfvén wave is confirmed.

  2. A global simulation of ICRF heating in a 3D magnetic configuration

    A global simulation code for the ICRF heating analysis in a three-dimensional (3D) magnetic configuration is developed combining two global simulation codes; a drift kinetic equation solver, GNET, and a wave field solver, TASK/WM. Both codes take into account 3D geometry using the numerically obtained 3D MHD equilibrium. The developed simulation code is applied to the LHD configuration as an example. Characteristics of energetic ion distributions in the phase space are clarified in LHD. The simulation results are also compared with experimental results by evaluating the count number of the neutral particle analyzer using the obtained energetic ion distribution, and a relatively good agreement is obtained. (author)

  3. Scaling laws of coronal loops compared to a 3D MHD model of an Active Region

    Bourdin, Philippe-A; Peter, Hardi

    2016-01-01

    Context. The structure and heating of coronal loops are investigated since decades. Established scaling laws relate fundamental quantities like the loop apex temperature, pressure, length, and the coronal heating. Aims. We test such scaling laws against a large-scale 3D MHD model of the Solar corona, which became feasible with nowadays high-performance computing. Methods. We drive an active region simulation a with photospheric observations and found strong similarities to the observed coronal loops in X-rays and EUV wavelength. A 3D reconstruction of stereoscopic observations showed that our model loops have a realistic spatial structure. We compare scaling laws to our model data extracted along an ensemble of field lines. Finally, we fit a new scaling law that represents well hot loops and also cooler structures, which was not possible before only based on observations. Results. Our model data gives some support for scaling laws that were established for hot and EUV-emissive coronal loops. For the RTV scali...

  4. 3D MHD Simulations of Planet Migration in Turbulent Stratified Disks

    Uribe, Ana; Flock, Mario; Henning, Thomas

    2011-01-01

    We performed 3D MHD simulations of planet migration in stratified disks using the Godunov code PLUTO, where the disk is turbulent due to the magnetorotational instability. We study the migration for planets with different planet-star mass ratios $q=M_{p}/M_{s}$. In agreement with previous studies, for the low-mass planet cases ($q=5\\times10^{-6}$ and $10^{-5}$), migration is dominated by random fluctuations in the torque. For a Jupiter-mass planet $(q=M_{p}/M_{s}=10^{-3}$ for $M_{s}=1M_{\\odot})$, we find a reduction of the magnetic stress inside the orbit of the planet and around the gap region. After an initial stage where the torque on the planet is positive, it reverses and we recover migration rates similar to those found in disks where the turbulent viscosity is modelled by an $\\alpha$ viscosity. For the intermediate-mass planets ($q=5\\times10^{-5}, 10^{-4}$ and $2\\times10^{-4}$) we find a new and so far unexpected behavior. In some cases they experience sustained and systematic outwards migration for th...

  5. Coronal energy input and dissipation in a solar active region 3D MHD model

    Bourdin, Philippe-A; Peter, Hardi

    2015-01-01

    Context. We have conducted a 3D MHD simulation of the solar corona above an active region in full scale and high resolution, which shows coronal loops, and plasma flows within them, similar to observations. Aims. We want to find the connection between the photospheric energy input by field-line braiding with the coronal energy conversion by Ohmic dissipation of induced currents. Methods. To this end we compare the coronal energy input and dissipation within our simulation domain above different fields of view, e.g. for a small loops system in the active region (AR) core. We also choose an ensemble of field lines to compare, e.g., the magnetic energy input to the heating per particle along these field lines. Results. We find an enhanced Ohmic dissipation of currents in the corona above areas that also have enhanced upwards-directed Poynting flux. These regions coincide with the regions where hot coronal loops within the AR core are observed. The coronal density plays a role in estimating the coronal temperatur...

  6. Radiative 3D MHD simulations of the spontaneous small-scale eruptions in the solar atmosphere

    Kitiashvili, Irina N.

    2015-08-01

    Studying non-linear turbulent dynamics of the solar atmosphere is important for understanding mechanism of the solar and stellar brightness variations. High-resolution observations of the quiet Sun reveal ubiquitous distributions of high-speed jets, which are transport mass and energy into the solar corona and feeding the solar wind. However, the origin of these eruption events is still unknown. Using 3D realistic MHD numerical simulations we find that small-scale eruptions are produced by ubiquitous magnetized vortex tubes generated by the Sun's turbulent convection in subsurface layers. The swirling vortex tubes (resembling tornadoes) penetrate into the solar atmosphere, capture and stretch background magnetic field, and push the surrounding material up, generating shocks. Our simulations reveal complicated high-speed flow patterns and thermodynamic and magnetic structure in the erupting vortex tubes and shows that the eruptions are initiated in the subsurface layers and are driven by high-pressure gradients in the subphotosphere and photosphere and by the Lorentz force in the higher atmosphere layers. I will discuss about properties of these eruptions, their effects on brightness and spectral variations and comparison with observations.

  7. Linearly perturbed MHD equilibria and 3D eddy current coupling via the control surface method

    Portone, A.; Villone, F.; Liu, Y.; Albanese, R.; Rubinacci, G.

    2008-08-01

    In this paper, a coupling strategy based on the control surface concept is used to self-consistently couple linear MHD solvers to 3D codes for the eddy current computation of eddy currents in the metallic structures surrounding the plasma. The coupling is performed by assuming that the plasma inertia (and, with it, all Alfven wave-like phenomena) can be neglected on the time scale of interest, which is dictated by the relevant electromagnetic time of the metallic structures. As is shown, plasma coupling with the metallic structures results in perturbations to the inductance matrix operator. In particular, by adopting the Fourier decomposition in poloidal and toroidal modes, it turns out that each toroidal mode can be associated with a matrix (additively) perturbing the inductance matrix that commonly describes the magnetic coupling of currents in vacuum. In this way, the treatment of resistive wall modes instabilities of various toroidal mode numbers and their possible cross-talk through the currents induced in the metallic structures can be easily studied.

  8. Global Value Chains from a 3D Printing Perspective

    Laplume, André O; Petersen, Bent; Pearce, Joshua M.

    2016-01-01

    This article outlines the evolution of additive manufacturing technology, culminating in 3D printing and presents a vision of how this evolution is affecting existing global value chains (GVCs) in production. In particular, we bring up questions about how this new technology can affect the...... whether in some industries diffusion of 3D printing technologies may change the role of multinational enterprises as coordinators of GVCs by inducing the engagement of a wider variety of firms, even households....

  9. 3D Multifluid MHD simulation for Uranus and Neptune: the seasonal variations of their magnetosphere

    Cao, X.; Paty, C. S.

    2015-12-01

    The interaction between Uranus' intrinsic magnetic field and the solar wind is quite different from the magnetospheric interactions of other planets. Uranus' large obliquity, coupled with the fact that its dipole moment is off-centered and highly tilted relative to the rotation axis, leads to unique and seasonally dependent interaction geometries with the solar wind. We present results from adapting a multifluid MHD simulation to examine these seasonally dependent geometries in terms of the global magnetospheric structure, magnetopause and bow shock location, and magnetotail configuration. The Voyager 2 spacecraft encountered Uranus near solstice, and was able to observe the magnetic field structure and plasma characteristics of a twisted magnetotail [Behannon et al., 1987]. We use such magnetometer and plasma observations as a basis for benchmarking our simulations for the solstice scenario. Auroral observations made by the Hubble Space Telescope during equinox [Lamy et al.,2012] give some indication of the magnetospheric interaction with the solar wind. We also demonstrate the structural difference of the magnetosphere between solstice and equinox seasons. The magnetosphere at equinox is quite distinct due to the orientation and rotation of the magnetic axis relative to the solar wind direction.

  10. 3-D resistive MHD calculations for tokamak plasmas: beyond the simple reduced set of equations

    Numerical studies of the resistive stability of tokamak plasmas in cylindrical geometry have been performed using: (1) the full set of resistive Magnetohydrodynamic (MHD) equations and (2) an extended version of the reduced set of resistive MHD equations including diamagnetic and electron temperature effects. In particular, the nonlinear interaction of tearing modes of many helicities has been investigated. The numerical results confirm many of the features uncovered previously using the simple reduced equations. (author)

  11. MODELING STATISTICAL PROPERTIES OF SOLAR ACTIVE REGIONS THROUGH DIRECT NUMERICAL SIMULATIONS OF 3D-MHD TURBULENCE

    Malapaka, Shiva Kumar; Mueller, Wolf-Christian [Max-Planck Institute for Plasma Physics, Boltzmannstrasse 2, D-85748 Garching bei Muenchen (Germany)

    2013-09-01

    Statistical properties of the Sun's photospheric turbulent magnetic field, especially those of the active regions (ARs), have been studied using the line-of-sight data from magnetograms taken by the Solar and Heliospheric Observatory and several other instruments. This includes structure functions and their exponents, flatness curves, and correlation functions. In these works, the dependence of structure function exponents ({zeta}{sub p}) of the order of the structure functions (p) was modeled using a non-intermittent K41 model. It is now well known that the ARs are highly turbulent and are associated with strong intermittent events. In this paper, we compare some of the observations from Abramenko et al. with the log-Poisson model used for modeling intermittent MHD turbulent flows. Next, we analyze the structure function data obtained from the direct numerical simulations (DNS) of homogeneous, incompressible 3D-MHD turbulence in three cases: sustained by forcing, freely decaying, and a flow initially driven and later allowed to decay (case 3). The respective DNS replicate the properties seen in the plots of {zeta}{sub p} against p of ARs. We also reproduce the trends and changes observed in intermittency in flatness and correlation functions of ARs. It is suggested from this analysis that an AR in the onset phase of a flare can be treated as a forced 3D-MHD turbulent system in its simplest form and that the flaring stage is representative of decaying 3D-MHD turbulence. It is also inferred that significant changes in intermittency from the initial onset phase of a flare to its final peak flaring phase are related to the time taken by the system to reach the initial onset phase.

  12. MODELING STATISTICAL PROPERTIES OF SOLAR ACTIVE REGIONS THROUGH DIRECT NUMERICAL SIMULATIONS OF 3D-MHD TURBULENCE

    Statistical properties of the Sun's photospheric turbulent magnetic field, especially those of the active regions (ARs), have been studied using the line-of-sight data from magnetograms taken by the Solar and Heliospheric Observatory and several other instruments. This includes structure functions and their exponents, flatness curves, and correlation functions. In these works, the dependence of structure function exponents (ζp) of the order of the structure functions (p) was modeled using a non-intermittent K41 model. It is now well known that the ARs are highly turbulent and are associated with strong intermittent events. In this paper, we compare some of the observations from Abramenko et al. with the log-Poisson model used for modeling intermittent MHD turbulent flows. Next, we analyze the structure function data obtained from the direct numerical simulations (DNS) of homogeneous, incompressible 3D-MHD turbulence in three cases: sustained by forcing, freely decaying, and a flow initially driven and later allowed to decay (case 3). The respective DNS replicate the properties seen in the plots of ζp against p of ARs. We also reproduce the trends and changes observed in intermittency in flatness and correlation functions of ARs. It is suggested from this analysis that an AR in the onset phase of a flare can be treated as a forced 3D-MHD turbulent system in its simplest form and that the flaring stage is representative of decaying 3D-MHD turbulence. It is also inferred that significant changes in intermittency from the initial onset phase of a flare to its final peak flaring phase are related to the time taken by the system to reach the initial onset phase

  13. Global smooth solutions of MHD equations with large data

    Lin, Yurui; Zhang, Huali; Zhou, Yi

    2016-07-01

    In this paper, we establish the global existence of smooth solutions of the three-dimensional MHD system for a class of large initial data. Both the initial velocity and magnetic field can be arbitrarily large in the critical norm.

  14. Plasma response measurements of external magnetic perturbations using electron cyclotron emission and comparisons to 3D ideal MHD equilibrium

    Willensdorfer, M; Strumberger, E; Suttrop, W; Vanovac, B; Brida, D; Cavedon, M; Classen, I; Dunne, M; Fietz, S; Fischer, R; Kirk, A; Laggner, F M; Liu, Y Q; Odstrcil, T; Ryan, D A; Viezzer, E; Zohm, H; Luhmann, I C

    2016-01-01

    The plasma response from an external n = 2 magnetic perturbation field in ASDEX Upgrade has been measured using mainly electron cyclotron emission (ECE) diagnostics and a rigid rotating field. To interpret ECE and ECE-imaging (ECE-I) measurements accurately, forward modeling of the radiation transport has been combined with ray tracing. The measured data is compared to synthetic ECE data generated from a 3D ideal magnetohydrodynamics (MHD) equilibrium calculated by VMEC. The measured amplitudes of the helical displacement in the midplane are in reasonable agreement with the one from the synthetic VMEC diagnostics. Both exceed the vacuum field calculations and indicate the presence of an amplified kink response at the edge. Although the calculated magnetic structure of this edge kink peaks at poloidal mode numbers larger than the resonant components |m| > |nq|, the displacement measured by ECE-I is almost resonant |m| ~ |nq|. This is expected from ideal MHD in the proximity of rational surfaces. VMEC and MARS-...

  15. Global MHD Mode Stabilization and Control for Tokamak Disruption Avoidance

    Sabbagh, S. A.; Berkery, J. W.; Bialek, J. M.; Hanson, J. M.; Park, Y. S.; Bell, R. E.; Gates, D. A.; Gerhardt, S. P.; Goumiri, I.; Grierson, B.; Holcomb, C.

    2015-11-01

    The near-complete elimination of plasma disruptions in fusion-producing tokamaks is the present ``grand challenge'' for stability research. Meeting this goal requires multiple approaches, important components of which are prediction, stabilization, and control of global MHD instabilities. Research on NSTX and its upgrade is synergizing these elements to make quantified progress on this challenge. Initial results from disruption characterization and prediction analyses describe physical disruption event chains in NSTX. Analysis of NSTX and DIII-D experiments show that stabilization of global modes is dominated by precession drift and bounce orbit resonances respectively. Stability therefore depends on the plasma rotation profile. A model-based rotation profile controller for NSTX-U using both neutral beams and neoclassical toroidal viscosity is shown in simulations to evolve profiles away from unstable states. Active RWM control is addressed using dual field component sensor feedback and a model-based RWM state-space controller. Comparison of measurements and synthetic diagnostics is examined for off-normal event handling. A planned 3D coil system upgrade can allow RWM control close to the ideal n = 1 with-wall limit. Supported by US DOE Contracts DE-FG02-99ER54524 and DE-AC02-09CH11466.

  16. Behaviour of magnetic islands in 3D MHD equilibria of helical devices

    Magnetic island formation in 3D finite-β equilibria in the H-1 Heliac is studied by using the HINT code. It is found that the size of a dangerous island should increase with β but that destruction of the equilibrium at low β is avoided because the rotational transform evolves to exclude the rational surface concerned. At higher β there is evidence of near-resonant flux surface deformations which may lead to an equilibrium limit. A reconnected equilibrium at still higher β exhibits a double island structure which is similar to homoclinic phase portraits which have been observed after separatrix reconnection in Hamiltonian systems. The physical mechanisms of island formation in finite-β helical equilibria have been investigated. The HINT code predicts that the global effect to the Pfirsch-Schlueter currents can lead to self-healing of magnetic islands independent of whether or not the plasma is stable to resistive interchange modes. This result has been compared with the predictions of a boundary-layer analysis which has been extended to consider configurations with islands in the vacuum magnetic field. (author). 5 refs, 1 fig

  17. A gravimetric 3D global inversion for cavity detection

    Camacho, A.G.; Vieira, R.; Montesinos, F.G. (Inst. de Astronomia y Geodesia, Madrid (Spain). Faculty de CC. Matematicas); Cuellar, V. (Lab. de Geotecnia del Centro de Estudios y Experimentacion de Obras Publicas, Madrid (Spain))

    1994-02-01

    A gravimetric survey, covering a site 200 m square, was carried out in order to locate karstic cavities. After eliminating the regional trend using a polynomial fit, the residual is modeled by least-squares prediction. Correlated signals for several wavelengths are detected. The inversion of these anomalies is performed by a global 3D adjustment using spherical bodies as models. The adjustment is repeated in order to obtain a stable configuration. The results show the probable presence of a system of cavities and galleries. Data collected from boreholes and the subsequent appearance of sink-holes are consistent with the results.

  18. 3-D Rigid Models from Partial Views - Global Factorization

    Aguiar, Pedro M Q; Gonçalves, Bruno B

    2010-01-01

    The so-called factorization methods recover 3-D rigid structure from motion by factorizing an observation matrix that collects 2-D projections of features. These methods became popular due to their robustness - they use a large number of views, which constrains adequately the solution - and computational simplicity - the large number of unknowns is computed through an SVD, avoiding non-linear optimization. However, they require that all the entries of the observation matrix are known. This is unlikely to happen in practice, due to self-occlusion and limited field of view. Also, when processing long videos, regions that become occluded often appear again later. Current factorization methods process these as new regions, leading to less accurate estimates of 3-D structure. In this paper, we propose a global factorization method that infers complete 3-D models directly from the 2-D projections in the entire set of available video frames. Our method decides whether a region that has become visible is a region tha...

  19. A Novel High-Order, Entropy Stable, 3D AMR MHD Solver with Guaranteed Positive Pressure

    Derigs, Dominik; Gassner, Gregor J; Walch, Stefanie

    2016-01-01

    We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum, and energy and is entropy stable. The method is designed to treat the divergence-free constraint on the magnetic field in a similar fashion to a hyperbolic divergence cleaning technique. The solver described herein is especially well-suited for flows involving strong discontinuities. Furthermore, we present a new formulation to guarantee positivity of the pressure. We present the underlying theory and implementation of the new solver into the multi-physics, multi-scale adaptive mesh refinement (AMR) simulation code $\\texttt{FLASH}$ (http://flash.uchicago.edu). The accuracy, robustness and computational efficiency is demonstrated with a number of tests, including comparisons to available MHD implementations in $\\texttt{FLASH}$.

  20. Virial theorem analysis of 3D numerical simulations of MHD self-gravitating turbulence

    Shadmehri, Mohsen; Vazquez-Semadeni, Enrique; Ballesteros-Paredes, Javier

    2001-01-01

    We discuss the virial balance of all members of a cloud ensemble in numerical simulations of self-gravitating MHD turbulence. We first discuss the choice of reference frame for evaluating the terms entering the virial theorem (VT), concluding that the balance of each cloud should be measured in its own reference frame. We then report preliminary results suggesting that a) the clouds are far from virial equilibrium, with the ``geometric'' (time derivative) terms dominating the VT. b) The surfa...

  1. A comparative study on 3-D solar wind structure observed by Ulysses and MHD simulation

    FENG Xueshang; XIANG Changqing; ZHONG Dingkun; FAN Quanlin

    2005-01-01

    During Ulysses' first rapid pole-to-pole transit from September 1994 to June 1995, its observations showed that middle- or high-speed solar winds covered all latitudes except those between -20° and +20° near the ecliptic plane,where the velocity was 300-450 km/s. At poleward 40°,however, only fast solar winds at the speed of 700-870 km/s were observed. In addition, the transitions from low-speed wind to high-speed wind or vice versa were abrupt. In this paper, the large-scale structure of solar wind observed by Ulysses near solar minimum is simulated by using the three-dimensional numerical MHD model. The model combines TVD Lax-Friedrich scheme and MacCormack Ⅱ scheme and decomposes the calculation region into two regions: one from 1 to 22 Rs and the other from 18 Rs to 1 AU.Based on the observations of the solar photospheric magnetic field and an addition of the volumetric heating to MHD equations, the large-scale solar wind structure mentioned above is reproduced by using the three-dimensional MHD model and the numerical results are roughly consistent with Ulysses' observations. Our simulation shows that the initial magnetic field topology and the addition of volume heating may govern the bimodal structure of solar wind observed by Ulysses and also demonstrates that the three-dimensional MHD numerical model used here is efficient in modeling the large-scale solar wind structure.

  2. An Arbitrary Lagrangian-Eulerian Discretization of MHD on 3D Unstructured Grids

    Rieben, R N; White, D A; Wallin, B K; Solberg, J M

    2006-06-12

    We present an arbitrary Lagrangian-Eulerian (ALE) discretization of the equations of resistive magnetohydrodynamics (MHD) on unstructured hexahedral grids. The method is formulated using an operator-split approach with three distinct phases: electromagnetic diffusion, Lagrangian motion, and Eulerian advection. The resistive magnetic dynamo equation is discretized using a compatible mixed finite element method with a 2nd order accurate implicit time differencing scheme which preserves the divergence-free nature of the magnetic field. At each discrete time step, electromagnetic force and heat terms are calculated and coupled to the hydrodynamic equations to compute the Lagrangian motion of the conducting materials. By virtue of the compatible discretization method used, the invariants of Lagrangian MHD motion are preserved in a discrete sense. When the Lagrangian motion of the mesh causes significant distortion, that distortion is corrected with a relaxation of the mesh, followed by a 2nd order monotonic remap of the electromagnetic state variables. The remap is equivalent to Eulerian advection of the magnetic flux density with a fictitious mesh relaxation velocity. The magnetic advection is performed using a novel variant of constrained transport (CT) that is valid for unstructured hexahedral grids with arbitrary mesh velocities. The advection method maintains the divergence free nature of the magnetic field and is second order accurate in regions where the solution is sufficiently smooth. For regions in which the magnetic field is discontinuous (e.g. MHD shocks) the method is limited using a novel variant of algebraic flux correction (AFC) which is local extremum diminishing (LED) and divergence preserving. Finally, we verify each stage of the discretization via a set of numerical experiments.

  3. Energetic particle effects on global MHD modes

    The effects of energetic particles on MHD type modes are studied by analytical theories and the nonvariational kinetic-MHD stability code (NOVA-K). In particular we address the problems of (1) the stabilization of ideal MHD internal kink modes and the excitation of resonant ''fishbone'' internal modes and (2) the alpha particle destabilization of toroidicity-induced Alfven eigenmodes (TAE) via transit resonances. Analytical theories are presented to help explain the NOVA-K results. For energetic trapped particles generated by neutral-beam injection (NBI) or ion cyclotron resonant heating (ICRH), a stability window for the n=1 internal kink mode in the hot particle beat space exists even in the absence of core ion finite Larmor radius effect (finite ω*i). On the other hand, the trapped alpha particles are found to resonantly excite instability of the n=1 internal mode and can lower the critical beta threshold. The circulating alpha particles can strongly destabilize TAE modes via inverse Landau damping associated with the spatial gradient of the alpha particle pressure. 23 refs., 5 figs

  4. First 3D MHD simulation of a massive-star magnetosphere with application to H\\alpha emission from \\theta^1 Ori C

    ud-Doula, Asif; Owocki, Stanley P; Petit, Veronique; Townsend, Richard H D

    2012-01-01

    We present the first fully 3D MHD simulation for magnetic channeling and confinement of a radiatively driven, massive-star wind. The specific parameters are chosen to represent the prototypical slowly rotating magnetic O star \\theta^1 Ori C, for which centrifugal and other dynamical effects of rotation are negligible. The computed global structure in latitude and radius resembles that found in previous 2D simulations, with unimpeded outflow along open field lines near the magnetic poles, and a complex equatorial belt of inner wind trapping by closed loops near the stellar surface, giving way to outflow above the Alfv\\'{e}n radius. In contrast to this previous 2D work, the 3D simulation described here now also shows how this complex structure fragments in azimuth, forming distinct clumps of closed loop infall within the Alfv\\'{e}n radius, transitioning in the outer wind to radial spokes of enhanced density with characteristic azimuthal separation of $15-20 \\degr$. Applying these results in a 3D code for line r...

  5. A novel high-order, entropy stable, 3D AMR MHD solver with guaranteed positive pressure

    Derigs, Dominik; Winters, Andrew R.; Gassner, Gregor J.; Walch, Stefanie

    2016-07-01

    We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum, and energy and is entropy stable. The method is designed to treat the divergence-free constraint on the magnetic field in a similar fashion to a hyperbolic divergence cleaning technique. The solver described herein is especially well-suited for flows involving strong discontinuities. Furthermore, we present a new formulation to guarantee positivity of the pressure. We present the underlying theory and implementation of the new solver into the multi-physics, multi-scale adaptive mesh refinement (AMR) simulation code FLASH (http://flash.uchicago.edu)

  6. A novel high-order, entropy stable, 3D AMR MHD solver with guaranteed positive pressure

    Derigs, Dominik; Winters, Andrew R.; Gassner, Gregor J.; Walch, Stefanie

    2016-07-01

    We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum, and energy and is entropy stable. The method is designed to treat the divergence-free constraint on the magnetic field in a similar fashion to a hyperbolic divergence cleaning technique. The solver described herein is especially well-suited for flows involving strong discontinuities. Furthermore, we present a new formulation to guarantee positivity of the pressure. We present the underlying theory and implementation of the new solver into the multi-physics, multi-scale adaptive mesh refinement (AMR) simulation code FLASH (http://flash.uchicago.edu).

  7. MHD and deep mixing in evolved stars. 1. 2D and 3D analytical models for the AGB

    Nucci, M C

    2014-01-01

    The advection of thermonuclear ashes by magnetized domains emerging from near the H-shell was suggested to explain AGB star abundances. Here we verify this idea quantitatively through exact MHD models. Starting with a simple 2D geometry and in an inertia frame, we study plasma equilibria avoiding the complications of numerical simulations. We show that, below the convective envelope of an AGB star, variable magnetic fields induce a natural expansion, permitted by the almost ideal MHD conditions, in which the radial velocity grows as the second power of the radius. We then study the convective envelope, where the complexity of macro-turbulence allows only for a schematic analytical treatment. Here the radial velocity depends on the square root of the radius. We then verify the robustness of our results with 3D calculations for the velocity, showing that, for both the studied regions, the solution previously found can be seen as a planar section of a more complex behavior, in which anyway the average radial vel...

  8. Non-twist map bifurcation of drift-lines and drift-island formation in saturated 3D MHD equilibria

    Pfefferle, David; Cooper, Wilfred A.; Graves, Jonathan P.

    2015-11-01

    Based on non-canonical perturbation theory, guiding-centre drift equations are identified as perturbed magnetic field-line equations. The topology of passing-particle orbits, called drift-lines, is completely determined by the magnetic configuration. In axisymmetric tokamak fields, drift-lines lie on shifted flux-surfaces, called drift-surfaces. Field-lines and drift-lines are subject to island structures at rational surfaces only when a non-axisymmetric component is added. The picture is different in the case of 3D saturated MHD equilibrium like the helical core associated with a non-resonant internal kink mode. In assuming nested flux-surfaces, these bifurcated states, expected for a reversed q-profile with qmin close yet above unity and conveniently obtained in VMEC, feature integrable field-lines. The helical drift-lines however become resonant with the axisymmetric component in the region of qmin and spontaneously generate drift-islands. Due to the locally reversed sheared q-profile, the drift-island structure follows the bifurcation/reconnection mechanism of non-twist maps. This result provides a theoretical interpretation of NBI fast ion helical hot-spots in Long-Lived Modes as well as snake-like impurity density accumulation in internal MHD activity.

  9. 3D relaxation MHD modeling with FOI-PERFECT code for electromagnetically driven HED systems

    Wang, Ganghua; Duan, Shuchao; Xie, Weiping; Kan, Mingxian; Institute of Fluid Physics Collaboration

    2015-11-01

    One of the challenges in numerical simulations of electromagnetically driven high energy density (HED) systems is the existence of vacuum region. The electromagnetic part of the conventional model adopts the magnetic diffusion approximation (magnetic induction model). The vacuum region is approximated by artificially increasing the resistivity. On one hand the phase/group velocity is superluminal and hence non-physical in the vacuum region, on the other hand a diffusion equation with large diffusion coefficient can only be solved by implicit scheme. Implicit method is usually difficult to parallelize and converge. A better alternative is to solve the full electromagnetic equations for the electromagnetic part. Maxwell's equations coupled with the constitutive equation, generalized Ohm's law, constitute a relaxation model. The dispersion relation is given to show its transition from electromagnetic propagation in vacuum to resistive MHD in plasma in a natural way. The phase and group velocities are finite for this system. A better time stepping is adopted to give a 3rd full order convergence in time domain without the stiff relaxation term restriction. Therefore it is convenient for explicit & parallel computations. Some numerical results of FOI-PERFECT code are also given. Project supported by the National Natural Science Foundation of China (Grant No. 11172277,11205145).

  10. Relativistic MHD simulations of core-collapse GRB jets: 3D instabilities and magnetic dissipation

    Bromberg, Omer

    2015-01-01

    Relativistic jets naturally occur in astrophysical systems that involve accretion onto compact objects, such as core collapse of massive stars in gamma-ray bursts (GRBs) and accretion onto supermassive black holes in active galactic nuclei (AGN). It is generally accepted that these jets are powered electromagnetically, by the magnetised rotation of a central compact object. However, how they produce the observed emission and survive the propagation for many orders of magnitude in distance without being disrupted by current-driven non-axisymmetric instabilities is the subject of active debate. We carry out time-dependent 3D relativistic magnetohydrodynamic simulations of relativistic, Poynting flux dominated jets. The jets are launched self-consistently by the rotation of a strongly magnetised central compact object. This determines the natural degree of azimuthal magnetic field winding, a crucial factor that controls jet stability. We find that the jets are susceptible to two types of instability: (i) a globa...

  11. The 3D MHD code GOEMHD3 for astrophysical plasmas with large Reynolds numbers. Code description, verification, and computational performance

    Skála, J.; Baruffa, F.; Büchner, J.; Rampp, M.

    2015-08-01

    Context. The numerical simulation of turbulence and flows in almost ideal astrophysical plasmas with large Reynolds numbers motivates the implementation of magnetohydrodynamical (MHD) computer codes with low resistivity. They need to be computationally efficient and scale well with large numbers of CPU cores, allow obtaining a high grid resolution over large simulation domains, and be easily and modularly extensible, for instance, to new initial and boundary conditions. Aims: Our aims are the implementation, optimization, and verification of a computationally efficient, highly scalable, and easily extensible low-dissipative MHD simulation code for the numerical investigation of the dynamics of astrophysical plasmas with large Reynolds numbers in three dimensions (3D). Methods: The new GOEMHD3 code discretizes the ideal part of the MHD equations using a fast and efficient leap-frog scheme that is second-order accurate in space and time and whose initial and boundary conditions can easily be modified. For the investigation of diffusive and dissipative processes the corresponding terms are discretized by a DuFort-Frankel scheme. To always fulfill the Courant-Friedrichs-Lewy stability criterion, the time step of the code is adapted dynamically. Numerically induced local oscillations are suppressed by explicit, externally controlled diffusion terms. Non-equidistant grids are implemented, which enhance the spatial resolution, where needed. GOEMHD3 is parallelized based on the hybrid MPI-OpenMP programing paradigm, adopting a standard two-dimensional domain-decomposition approach. Results: The ideal part of the equation solver is verified by performing numerical tests of the evolution of the well-understood Kelvin-Helmholtz instability and of Orszag-Tang vortices. The accuracy of solving the (resistive) induction equation is tested by simulating the decay of a cylindrical current column. Furthermore, we show that the computational performance of the code scales very

  12. A 3-D MHD equilibrium description of nonlinearly saturated ideal external kink/peeling structures in tokamaks

    Cooper, W. A.; Graves, J. P.; Duval, B. P.; Porte, L.; Reimerdes, H.; Sauter, O.; Tran, T.-M.

    2015-12-01

    > Novel free boundary magnetohydrodynamic equilibrium states with spontaneous three-dimensional (3-D) deformations of the plasma-vacuum interface are computed. The structures obtained look like saturated ideal external kink/peeling modes. Large edge pressure gradients yield toroidal mode number distortions when the edge bootstrap current is large and higher corrugations when this current is small. Linear ideal MHD stability analyses confirm the nonlinear saturated ideal kink equilibrium states produced and we can identify the Pfirsch-Schlüter current as the main linear instability driving mechanism when the edge pressure gradient is large. The dominant non-axisymmetric component of this Pfirsch-Schlüter current drives a near resonant helical parallel current density ribbon that aligns with the near vanishing magnetic shear region caused by the edge bootstrap current. This current ribbon is a manifestation of the outer mode previously found on JET (Solano 2010). We claim that the equilibrium corrugations describe structures that are commonly observed in quiescent H-mode tokamak discharges.

  13. Toward the analysis of waves in the solar atmosphere based on NLTE spectral synthesis from 3D MHD simulations

    Haberreiter, M; McIntosh, S; Wedemeyer-Boehm, S

    2010-01-01

    From the analysis of Dopplergrams in the K I 7699 A and Na I 5890 A spectral lines observed with the Magneto-Optical filter at Two Heights (MOTH) experiment during the austral summer in 2002-03 we find upward traveling waves in magnetic regions. Our analysis shows that the dispersion relation of these waves strongly depends on whether the wave is detected in the low-beta or high-beta regime. Moreover, the observed dispersion relation does not show the expected decrease of the acoustic cut-off frequency for the field guided slow magnetic wave. Instead, we detected an increase of the travel times below the acoustic cut-off frequency and at the same time a decrease of the travel time above it. To study the formation height of the spectral lines employed by MOTH in greater detail we are currently in the process of employing 3D MHD simulations carried out with CO5BOLD to perform NLTE spectral synthesis.

  14. Particle entry through "Sash" groove simulated by Global 3D Electromagnetic Particle code with duskward IMF By

    Yan, X.; Cai, D.; Nishikawa, K.; Lembege, B.

    2004-12-01

    We made our efforts to parallelize the global 3D HPF Electromagnetic particle model (EMPM) for several years and have also reported our meaningful simulation results that revealed the essential physics involved in interaction of the solar wind with the Earth's magnetosphere using this EMPM (Nishikawa et al., 1995; Nishikawa, 1997, 1998a, b, 2001, 2002) in our PC cluster and supercomputer(D.S. Cai et al., 2001, 2003). Sash patterns and related phenomena have been observed and reported in some satellite observations (Fujumoto et al. 1997; Maynard, 2001), and have motivated 3D MHD simulations (White and al., 1998). We also investigated it with our global 3D parallelized HPF EMPM with dawnward IMF By (K.-I. Nishikawa, 1998) and recently new simulation with dusk-ward IMF By was accomplished in the new VPP5000 supercomputer. In the new simulations performed on the new VPP5000 supercomputer of Tsukuba University, we used larger domain size, 305×205×205, smaller grid size (Δ ), 0.5R E(the radium of the Earth), more total particle number, 220,000,000 (about 8 pairs per cell). At first, we run this code until we get the so-called quasi-stationary status; After the quasi-stationary status was established, we applied a northward IMF (B z=0.2), and then wait until the IMF arrives around the magnetopuase. After the arrival of IMF, we begin to change the IMF from northward to duskward (IMF B y=-0.2). The results revealed that the groove structure at the day-side magnetopause, that causes particle entry into inner magnetosphere and the cross structure or S-structure at near magneto-tail are formed. Moreover, in contrast with MHD simulations, kinetic characteristic of this event is also analyzed self-consistently with this simulation. The new simulation provides new and more detailed insights for the observed sash event.

  15. On the characterization of magnetic reconnection in global MHD simulations

    T. V. Laitinen

    2006-11-01

    Full Text Available The conventional definition of reconnection rate as the electric field parallel to an x-line is problematic in global MHD simulations for several reasons: the x-line itself may be hard to find in a non-trivial geometry such as at the magnetopause, and the lack of realistic resistivity modelling leaves us without reliable non-convective electric field. In this article we describe reconnection characterization methods that avoid those problems and are practical to apply in global MHD simulations. We propose that the reconnection separator line can be identified as the region where magnetic field lines of different topological properties meet, rather than by local considerations. The global convection associated with reconnection is then quantified by calculating the transfer of mass, energy or magnetic field across the boundary of closed and open field line regions. The extent of the diffusion region is determined from the destruction of electromagnetic energy, given by the divergence of the Poynting vector. Integrals of this energy conversion provide a way to estimate the total reconnection efficiency.

  16. Global MHD Simulations of Accretion Disks in Cataclysmic Variables (CVs): I. The Importance of Spiral Shocks

    Ju, Wenhua; Zhu, Zhaohuan

    2016-01-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 to that driven by spiral shock waves. Remarkably, we find that even with vigorous MRI turbulence, spiral shocks are an important component to 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 and both with and without mass inflow via the Lagrangian point (L1). Compared with previous similar studies, we find the following new results. 1) Linear wave dispersion relation fits the pitch angles of spiral density waves very well. 2) We demonstrate explicitly that mass accreti...

  17. Global stability analysis of turbulent 3D wakes

    Rigas, Georgios; Sipp, Denis; Juniper, Matthew

    2015-11-01

    At low Reynolds numbers, corresponding to laminar and transitional regimes, hydrodynamic stability theory has aided the understanding of the dynamics of bluff body wake-flows and the application of effective control strategies. However, flows of fundamental importance to many industries, in particular the transport industry, involve high Reynolds numbers and turbulent wakes. Despite their turbulence, such wake flows exhibit organisation which is manifested as coherent structures. Recent work has shown that the turbulent coherent structures retain the shape of the symmetry-breaking laminar instabilities and only those manifest as large-scale structures in the near wake (Rigas et al., JFM vol. 750:R5 2014, JFM vol. 778:R2 2015). Based on the findings of the persistence of the laminar instabilities at high Reynolds numbers, we investigate the global stability characteristics of a turbulent wake generated behind a bluff three-dimensional axisymmetric body. We perform a linear global stability analysis on the experimentally obtained mean flow and we recover the dynamic characteristics and spatial structure of the coherent structures, which are linked to the transitional instabilities. A detailed comparison of the predictions with the experimental measurements will be provided.

  18. A global sustainability perspective on 3D printing technologies

    Three-dimensional printing (3DP) represents a relative novel technology in manufacturing which is associated with potentially strong stimuli for sustainable development. Until now, research has merely assessed case study-related potentials of 3DP and described specific aspects of 3DP. This study represents the first comprehensive assessment of 3DP from a global sustainability perspective. It contains a qualitative assessment of 3DP-induced sustainability implications and quantifies changes in life cycle costs, energy and CO2 emissions globally by 2025. 3DP is identified to cost-effectively lower manufacturing inputs and outputs in markets with low volume, customized and high-value production chains as aerospace and medical component manufacturing. This lowers energy use, resource demands and related CO2 emissions over the entire product life cycle, induces changes in labour structures and generates shifts towards more digital and localized supply chains. The model calculations show that 3DP contains the potential to reduce costs by 170–593 billion US $, the total primary energy supply by 2.54–9.30 EJ and CO2 emissions by 130.5–525.5 Mt by 2025. The great range within the saving potentials can be explained with the immature state of the technology and the associated uncertainties of predicting market and technology developments. The energy and CO2 emission intensities of industrial manufacturing are reducible by maximally 5% through 3DP by 2025, as 3DP remains a niche technology. If 3DP was applicable to larger production volumes in consumer products or automotive manufacturing, it contains the (theoretical) potential to absolutely decouple energy and CO2 emission from economic activity. - Highlights: • Global sustainability aspects of 3DP in manufacturing are assessed in two ways. • 3DP will strongly influence manufacturing in aerospace, medical components, tooling. • 3DP re-shifts production to consumer countries due to decreased labour costs.

  19. Extracting 3D layout from a single image using global image structures.

    Lou, Zhongyu; Gevers, Theo; Hu, Ninghang

    2015-10-01

    Extracting the pixel-level 3D layout from a single image is important for different applications, such as object localization, image, and video categorization. Traditionally, the 3D layout is derived by solving a pixel-level classification problem. However, the image-level 3D structure can be very beneficial for extracting pixel-level 3D layout since it implies the way how pixels in the image are organized. In this paper, we propose an approach that first predicts the global image structure, and then we use the global structure for fine-grained pixel-level 3D layout extraction. In particular, image features are extracted based on multiple layout templates. We then learn a discriminative model for classifying the global layout at the image-level. Using latent variables, we implicitly model the sublevel semantics of the image, which enrich the expressiveness of our model. After the image-level structure is obtained, it is used as the prior knowledge to infer pixel-wise 3D layout. Experiments show that the results of our model outperform the state-of-the-art methods by 11.7% for 3D structure classification. Moreover, we show that employing the 3D structure prior information yields accurate 3D scene layout segmentation. PMID:25966478

  20. Study of the 3D Coronal Magnetic Field of Active Region 11117 Around the Time of a Confined Flare Using a Data-Driven CESE--MHD Model

    Jiang, Chaowei; Wu, S T; Hu, Qiang

    2012-01-01

    We apply a data-driven MHD model to investigate the three-dimensional (3D) magnetic field of NOAA active region (AR) 11117 around the time of a C-class confined flare occurred on 2010 October 25. The MHD model, based on the spacetime conservation-element and solution-element (CESE) scheme, is designed to focus on the magnetic-field evolution and to consider a simplified solar atomsphere with finite plasma $\\beta$. Magnetic vector-field data derived from the observations at the photoshpere is inputted directly to constrain the model. Assuming that the dynamic evolution of the coronal magnetic field can be approximated by successive equilibria, we solve a time sequence of MHD equilibria basing on a set of vector magnetograms for AR 11117 taken by the Helioseismic and Magnetic Imager (HMI) on board the {\\it Solar Dynamic Observatory (SDO)} around the time of flare. The model qualitatively reproduces the basic structures of the 3D magnetic field, as supported by the visual similarity between the field lines and t...

  1. One year in the Earth's magnetosphere: A global MHD simulation and spacecraft measurements

    Facsko, G; Zivkovic, T; Palin, L; Kallio, E; Agren, K; Opgenoorth, H; Tanskanen, E I; Milan, S E

    2016-01-01

    The response of the Earth's magnetosphere to changing solar wind conditions are studied with a 3D 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 is 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 $L_1$ Lagrange point to the foreshock, therefore the OMNIWeb data is appropriate input to the GUMICS-4. The Cluster SC3 footprints are dete...

  2. Self-organisation in protoplanetary disks: global, non-stratified Hall-MHD simulations

    Béthune, William; Ferreira, Jonathan

    2016-01-01

    Recent observations revealed organised structures in protoplanetary disks, such as axisymmetric rings or horseshoe concen- trations evocative of large-scale vortices. These structures are often interpreted as the result of planet-disc interactions. However, these disks are also known to be unstable to the magneto-rotational instability (MRI) which is believed to be one of the dominant angular momentum transport mechanism in these objects. It is therefore natural to ask if the MRI itself could produce these structures without invoking planets. The nonlinear evolution of the MRI is strongly affected by the low ionisation fraction in protoplanetary disks. The Hall effect in particular, which is dominant in dense and weakly ionised parts of these objects, has been shown to spontaneously drive self- organising flows in shearing box simulations. Here, we investigate the behaviour of global MRI-unstable disc models dominated by the Hall effect and characterise their dynamics. We perform 3D unstratified Hall-MHD simu...

  3. MHD Simulations of Global Accretion Disks with Vertical Magnetic Fields

    Suzuki, Takeru K

    2013-01-01

    (Abridged) We report results of three dimensional MHD simulations of global accretion disks threaded with weak vertical magnetic fields. We perform the simulations in the spherical coordinates with different temperature profiles and accordingly different rotation profiles. In the cases with a spatially constant temperature, because the rotation frequency is vertically constant in the equilibrium condition, general properties of the turbulence are quantitatively similar to those obtained in local shearing box simulations. On the other hand, in the cases with a radially variable temperature profile, the vertical differential rotation, which is inevitable in the equilibrium condition, winds up the magnetic field lines, in addition to the usual radial differential rotation. As a result, the coherent wound magnetic fields contribute to the Maxwell stress in the surface regions. We obtain nondimensional density and velocity fluctuations ~0.1-0.2 at the midplane. The azimuthal power spectra of the magnetic fields sh...

  4. Tela global e oceano azul: cinema 3D, o caso Brasil

    Andre Gatti

    2012-12-01

    Full Text Available A exibição de filmes de longa-metragem em 3D é a nova onda da indústria cinematográfica global. O que está acontecendo no Brasil? Quais problemas da implantação do 3D no país? Tudo ainda é muito recente, mas os resultados desta onda já podem ser percebidos.

  5. One year in the Earth's magnetosphere: A global MHD simulation and spacecraft measurements

    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.

  6. A global MHD model of the Jovian magnetosphere

    A new global model of the Jovian magnetosphere is proposed by using the magnetohydrodynamic (MHD) simulation, in which the interaction between a rapidly rotating magnetosphere and the external plasma flow is investigated. In this study, the steady solar wind with/without the interplanetary magnetic filed (IMF) is applied on the inflow boundary. The results clearly indicate that the global structure of the magnetic field lines of force is largely affected by the rotation of Jupiter as well as by the solar wind. It is also found that magnetic reconnection at the tail region generates the plasma outflow and greatly modifies the corotational structure of the plasmas in the magnetosphere. Especially when the relatively high speed solar wind or the solar wind with the northward IMF collides on the magnetosphere, strong magnetic reconnection at the tail induces the tail-dusk flow, which is against the corotation. On the other hand, the evident tail dusk flow does not appear in the case of the low speed solar wind or the case with the southward IMF. Thus, the internal structure of the rotating magnetosphere will be altered even qualitatively in response to the solar wind parameters. (author)

  7. A 3-D Global History from a Glober Identity ina Noncentric and Holistic Perspective

    2014-01-01

    In modern times, global and world history has been understood, writtenand taught from either a national perspective or a regional/cultural/ideological one in the light of centrism such as West/Eurocentrism eversince it was created. The author argues, in a noncentric and holisticperspective, that nationalized global or world histories are nothingbut distorting mirrors full of pride and prejudice. It is argued that bothobjective global history and subjective global history are contingentresultants in a certain special and temporal context. In the early 21stcentury, the author heralds a globalization of global history and "GlobalHistorians of All Countries Unite!" He also suggests that it is necessaryfor a global historian to develop a glober identity beyond her/hisnational identity and that a reasonable and intelligible global historywhich is closer to the objective global history should be a 3-D globalhistory of the glober, by the glober, for the glober.

  8. Local-global alignment for finding 3D similarities in protein structures

    Zemla, Adam T.

    2011-09-20

    A method of finding 3D similarities in protein structures of a first molecule and a second molecule. The method comprises providing preselected information regarding the first molecule and the second molecule. Comparing the first molecule and the second molecule using Longest Continuous Segments (LCS) analysis. Comparing the first molecule and the second molecule using Global Distance Test (GDT) analysis. Comparing the first molecule and the second molecule using Local Global Alignment Scoring function (LGA_S) analysis. Verifying constructed alignment and repeating the steps to find the regions of 3D similarities in protein structures.

  9. Global MHD Simulations of Accretion Disks in Cataclysmic Variables. I. The Importance of Spiral Shocks

    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.

  10. Global regular solutions for the 3D Zakharov-Kuznetsov equation posed on unbounded domains

    Larkin, N. A.

    2015-09-01

    An initial-boundary value problem for the 3D Zakharov-Kuznetsov equation posed on unbounded domains is considered. Existence and uniqueness of a global regular solution as well as exponential decay of the H2-norm for small initial data are proven.

  11. Global and Kinetic MHD Simulation by the Gpic-MHD Code

    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 φ 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.

  12. Characterisation of small-scale heating events in the solar atmosphere from 3D MHD simulations and their potential role in coronal heating

    Haberreiter, M.; Guerreiro, N.; Hansteen, V. H.; Schmutz, W. K.

    2015-12-01

    The physical mechanism that heats the solar corona is one of the still open science questions in solar physics. One of the proposed mechanism for coronal heating are nanoflares. To investigate their role in coronal heating we study the properties of the small-scale heating events in the solar atmosphere using 3D MHD simulations. We present a method to identify and track these heating events in time which allows us to study their life time, energy, and spectral signatures. These spectal signatures will be compared with available spectrosopic observations obtained with IRIS and SUMER. Ultimately, these results will be important for the coordinated scientific exploitation of SPICE and EUI along with other instruments onboard Solar Orbiter to address the coronal heating problem.

  13. Dense 3d Point Cloud Generation from Uav Images from Image Matching and Global Optimazation

    Rhee, S.; Kim, T.

    2016-06-01

    3D spatial information from unmanned aerial vehicles (UAV) images is usually provided in the form of 3D point clouds. For various UAV applications, it is important to generate dense 3D point clouds automatically from over the entire extent of UAV images. In this paper, we aim to apply image matching for generation of local point clouds over a pair or group of images and global optimization to combine local point clouds over the whole region of interest. We tried to apply two types of image matching, an object space-based matching technique and an image space-based matching technique, and to compare the performance of the two techniques. The object space-based matching used here sets a list of candidate height values for a fixed horizontal position in the object space. For each height, its corresponding image point is calculated and similarity is measured by grey-level correlation. The image space-based matching used here is a modified relaxation matching. We devised a global optimization scheme for finding optimal pairs (or groups) to apply image matching, defining local match region in image- or object- space, and merging local point clouds into a global one. For optimal pair selection, tiepoints among images were extracted and stereo coverage network was defined by forming a maximum spanning tree using the tiepoints. From experiments, we confirmed that through image matching and global optimization, 3D point clouds were generated successfully. However, results also revealed some limitations. In case of image-based matching results, we observed some blanks in 3D point clouds. In case of object space-based matching results, we observed more blunders than image-based matching ones and noisy local height variations. We suspect these might be due to inaccurate orientation parameters. The work in this paper is still ongoing. We will further test our approach with more precise orientation parameters.

  14. Signatures of small-scale heating events in EUV spectral lines as modeled from 3D MHD simulations

    Guerreiro, Nuno; Haberreiter, Margit; Hansteen, Viggo; Curdt, Werner; Schmutz, Werner

    2014-05-01

    We aim at understanding the implications of small scale heating events in the solar atmosphere for the variations of the solar spectral irradiance. We present a technique for identification and characterization of these events in 3D simulations of the solar atmosphere. An accurate property determination of these events in time and space will help us to understand how spectral lines, in particular in the EUV, respond to them and which kind of spectral signatures one would expect to find in observations as from SOHO/SUMER and eventually from future space missions, as for example observations by SPICE on board Solar Orbiter.

  15. The 3-D global spatial data model foundation of the spatial data infrastructure

    Burkholder, Earl F

    2008-01-01

    Traditional methods for handling spatial data are encumbered by the assumption of separate origins for horizontal and vertical measurements. Modern measurement systems operate in a 3-D spatial environment. The 3-D Global Spatial Data Model: Foundation of the Spatial Data Infrastructure offers a new model for handling digital spatial data, the global spatial data model or GSDM. The GSDM preserves the integrity of three-dimensional spatial data while also providing additional benefits such as simpler equations, worldwide standardization, and the ability to track spatial data accuracy with greater specificity and convenience. This groundbreaking spatial model incorporates both a functional model and a stochastic model to connect the physical world to the ECEF rectangular system. Combining horizontal and vertical data into a single, three-dimensional database, this authoritative monograph provides a logical development of theoretical concepts and practical tools that can be used to handle spatial data mo...

  16. Global solution for the 3D gravity water waves system above a flat bottom

    Wang, Xuecheng

    2015-01-01

    In this paper, we consider the $3D$ water waves system above a fixed flat bottom (finite depth) with gravity, meanwhile without surface tension (the gravity water waves system). For this system, we prove global existence for suitably small initial data and non-existence of traveling waves below certain level of smallness, which strongly contrasts the behavior of solution of the same system in the $2D$ case.

  17. Effect of the Interplanetary Electric Field on the Magnetopause From Global MHD Simulations

    HUANG Zhaohui; DING Kai; WANG Chi

    2012-01-01

    The north-south component B_z of the Interplanetary Magnetic Field(IMF) and solar wind dynamic pressure P_d are generally treated as the two main factors in the solar wind that determine the geometry of the magnetosphere.By using the 3D global MHD simulations,we investigate the effect of the Interplanetary Electric Field(IEF) on the size and shape of magnetopause quantitatively. Our numerical experiments confirm that the geometry of the magnetopause are mainly determined by P_d and B_z,as expected.However,the dawn-dusk IEFs have great impact on the magnetopause erosion because of the magnetic reconnection,thus affecting the size and shape of the magnetopause.Higher solar wind speed with the same B_z will lead to bigger dawn-dusk IEFs,which means the higher reconnection rate,and then results in more magnetic flux removal from the dayside. Consequently,the dayside magnetopause moves inward and flank magnetopause moves outward.

  18. Pressure driven currents near magnetic islands in 3D MHD equilibria: Effects of pressure variation within flux surfaces and of symmetry

    Reiman, Allan H.

    2016-07-01

    In toroidal, magnetically confined plasmas, the heat and particle transport is strongly anisotropic, with transport along the field lines sufficiently strong relative to cross-field transport that the equilibrium pressure can generally be regarded as constant on the flux surfaces in much of the plasma. The regions near small magnetic islands, and those near the X-lines of larger islands, are exceptions, having a significant variation of the pressure within the flux surfaces. It is shown here that the variation of the equilibrium pressure within the flux surfaces in those regions has significant consequences for the pressure driven currents. It is further shown that the consequences are strongly affected by the symmetry of the magnetic field if the field is invariant under combined reflection in the poloidal and toroidal angles. (This symmetry property is called "stellarator symmetry.") In non-stellarator-symmetric equilibria, the pressure-driven currents have logarithmic singularities at the X-lines. In stellarator-symmetric MHD equilibria, the singular components of the pressure-driven currents vanish. These equilibria are to be contrasted with equilibria having B ṡ∇p =0 , where the singular components of the pressure-driven currents vanish regardless of the symmetry. They are also to be contrasted with 3D MHD equilibrium solutions that are constrained to have simply nested flux surfaces, where the pressure-driven current goes like 1 /x near rational surfaces, where x is the distance from the rational surface, except in the case of quasi-symmetric flux surfaces. For the purpose of calculating the pressure-driven currents near magnetic islands, we work with a closed subset of the MHD equilibrium equations that involves only perpendicular force balance, and is decoupled from parallel force balance. It is not correct to use the parallel component of the conventional MHD force balance equation, B ṡ∇p =0 , near magnetic islands. Small but nonzero values of B

  19. External force back-projective composition and globally deformable optimization for 3-D coronary artery reconstruction

    The clinical value of the 3D reconstruction of a coronary artery is important for the diagnosis and intervention of cardiovascular diseases. This work proposes a method based on a deformable model for reconstructing coronary arteries from two monoplane angiographic images acquired from different angles. First, an external force back-projective composition model is developed to determine the external force, for which the force distributions in different views are back-projected to the 3D space and composited in the same coordinate system based on the perspective projection principle of x-ray imaging. The elasticity and bending forces are composited as an internal force to maintain the smoothness of the deformable curve. Second, the deformable curve evolves rapidly toward the true vascular centerlines in 3D space and angiographic images under the combination of internal and external forces. Third, densely matched correspondence among vessel centerlines is constructed using a curve alignment method. The bundle adjustment method is then utilized for the global optimization of the projection parameters and the 3D structures. The proposed method is validated on phantom data and routine angiographic images with consideration for space and re-projection image errors. Experimental results demonstrate the effectiveness and robustness of the proposed method for the reconstruction of coronary arteries from two monoplane angiographic images. The proposed method can achieve a mean space error of 0.564 mm and a mean re-projection error of 0.349 mm. (paper)

  20. A Global Hypothesis Verification Framework for 3D Object Recognition in Clutter.

    Aldoma, Aitor; Tombari, Federico; Stefano, Luigi Di; Vincze, Markus

    2016-07-01

    Pipelines to recognize 3D objects despite clutter and occlusions usually end up with a final verification stage whereby recognition hypotheses are validated or dismissed based on how well they explain sensor measurements. Unlike previous work, we propose a Global Hypothesis Verification (GHV) approach which regards all hypotheses jointly so as to account for mutual interactions. GHV provides a principled framework to tackle the complexity of our visual world by leveraging on a plurality of recognition paradigms and cues. Accordingly, we present a 3D object recognition pipeline deploying both global and local 3D features as well as shape and color. Thereby, and facilitated by the robustness of the verification process, diverse object hypotheses can be gathered and weak hypotheses need not be suppressed too early to trade sensitivity for specificity. Experiments demonstrate the effectiveness of our proposal, which significantly improves over the state-of-art and attains ideal performance (no false negatives, no false positives) on three out of the six most relevant and challenging benchmark datasets. PMID:26485476

  1. 3-D carotid multi-region MRI segmentation by globally optimal evolution of coupled surfaces.

    Ukwatta, Eranga; Yuan, Jing; Rajchl, Martin; Qiu, Wu; Tessier, David; Fenster, Aaron

    2013-04-01

    In this paper, we propose a novel global optimization based 3-D multi-region segmentation algorithm for T1-weighted black-blood carotid magnetic resonance (MR) images. The proposed algorithm partitions a 3-D carotid MR image into three regions: wall, lumen, and background. The algorithm performs such partitioning by simultaneously evolving two coupled 3-D surfaces of carotid artery adventitia boundary (AB) and lumen-intima boundary (LIB) while preserving their anatomical inter-surface consistency such that the LIB is always located within the AB. In particular, we show that the proposed algorithm results in a fully time implicit scheme that propagates the two linearly ordered surfaces of the AB and LIB to their globally optimal positions during each discrete time frame by convex relaxation. In this regard, we introduce the continuous max-flow model and prove its duality/equivalence to the convex relaxed optimization problem with respect to each evolution step. We then propose a fully parallelized continuous max-flow-based algorithm, which can be readily implemented on a GPU to achieve high computational efficiency. Extensive experiments, with four users using 12 3T MR and 26 1.5T MR images, demonstrate that the proposed algorithm yields high accuracy and low operator variability in computing vessel wall volume. In addition, we show the algorithm outperforms previous methods in terms of high computational efficiency and robustness with fewer user interactions. PMID:23303689

  2. Marine denitrification rates determined from a global 3-D inverse model

    T. DeVries

    2013-04-01

    Full Text Available A major impediment to understanding long-term changes in the marine nitrogen (N cycle is the persistent uncertainty about the rates, distribution, and sensitivity of its largest fluxes in the modern ocean. We use a global ocean circulation model to obtain the first 3-D estimate of marine denitrification rates that is maximally consistent with available observations of nitrate deficits and the nitrogen isotopic ratio of oceanic nitrate. We find a global rate of marine denitrification in suboxic waters and sediments of 120–240 Tg N yr−1, which is lower than many other recent estimates. The difference stems from the ability to represent the 3-D spatial structure of suboxic zones, where denitrification rates of 50–77 Tg N yr−1 result in up to 50% depletion of nitrate. This depletion reduces the effect of local isotopic enrichment on the rest of the ocean, allowing the N isotope ratio of oceanic nitrate to be achieved with a sedimentary denitrification rate about 1.3–2.3 times that of suboxic zones. This balance of N losses between sediments and suboxic zones is shown to obey a simple relationship between isotope fractionation and the degree of nitrate consumption in the core of the suboxic zones. The global denitrification rates derived here suggest that the marine nitrogen budget is likely close to balanced.

  3. Investigation of tetrafluoromethane as a plasma gas in a very high pressure/low current dc batch reactor by means of 3D MHD modelling

    This paper deals with 3D MHD modelling of the behaviour of a tetrafluoromethane (CF4) plasma arc in a batch reactor under peculiar conditions of low current (0.35 A) and very high pressure (50 atm). The first part of the manuscript presents results for a horizontal configuration of the reactor, as is undertaken experimentally. The model has led to the understanding of the instabilities observed experimentally for such unusual operating conditions. The curved shape of the arc and the sliding of the anodic arc root along the electrode have been revealed to be the source of the experimental instabilities. The latter part of the manuscript investigates the effect of two vertical configurations of the reactor; with a cathode at the top and cathode at the bottom to overcome the instabilities. In these reactor configurations, the arc is much more stable and stays centred in the middle of the electrodes. These configurations are more suitable for the stability of the arc discharge, but have to be verified experimentally. (paper)

  4. The global 3-D distribution of tropospheric aerosols as characterized by CALIOP

    D. M. Winker

    2013-03-01

    Full Text Available The CALIOP lidar, carried on the CALIPSO satellite, has been acquiring global atmospheric profiles since June 2006. This dataset now offers the opportunity to characterize the global 3-D distribution of aerosol as well as seasonal and interannual variations, and confront aerosol models with observations in a way that has not been possible before. With that goal in mind, a monthly global gridded dataset of daytime and nighttime aerosol extinction profiles has been constructed, available as a Level 3 aerosol product. Averaged aerosol profiles for cloud-free and all-sky conditions are reported separately. This 6-yr dataset characterizes the global 3-dimensional distribution of tropospheric aerosol. Vertical distributions are seen to vary with season, as both source strengths and transport mechanisms vary. In most regions, clear-sky and all-sky mean aerosol profiles are found to be quite similar, implying a lack of correlation between high semi-transparent cloud and aerosol in the lower troposphere. An initial evaluation of the accuracy of the aerosol extinction profiles is presented. Detection limitations and the representivity of aerosol profiles in the upper troposphere are of particular concern. While results are preliminary, we present evidence that the monthly-mean CALIOP aerosol profiles provide quantitative characterization of elevated aerosol layers in major transport pathways. Aerosol extinction in the free troposphere in clean conditions, where the true aerosol extinction is typically 0.001 km−1 or less, is generally underestimated, however. The work described here forms an initial global 3-D aerosol climatology which we plan to extend and improve over time.

  5. Global gyrokinetic models for energetic particle driven Alfvén instabilities in 3D equilibria

    Spong, Don; Holod, Ihor

    2015-11-01

    The GTC global gyrokinetic PIC model has been adapted to 3D VMEC equilibria and provides a new method for the analysis of Alfvénic instabilities in stellarators, 3D tokamaks, and helical RFP states. The gyrokinetic orderings (k||/k⊥ applicable to a range of energetic particle driven instabilities that have been observed in 3D configurations. Applications of this model to stellarators have indicated that a variety of different Alfvén instabilities can be excited, depending on the toroidal mode number, fast ion average energy and fast ion density profile. Both an LHD discharge where bursting n = 1 Alfvén activity in the TAE gap was observed and a W7-X case have been examined. TAE,/EAE/GAE modes have been found in the simulations, depending on the mode family and fast ion profiles used. The dynamical evolution of the instabilities shows the field period coupling between n and n + Nfp expected for a stellarator. The development of gyrofluid reduced models that can capture relevant physics aspects of the gyrokinetic models will also be discussed. Research sponsored by the U.S. Department of Energy under Contract DE-AC05-00OR22725 with UT-Battelle, LLC and the GSEP SciDAC Center.

  6. Effective solutions to a global 3D visual system in networking environments

    ZHANG; Liqiang; ZHANG; Yan; YANG; Chongjun; LIU; Suhong; R

    2005-01-01

    The development of Web-based global 3D visual system has made progresses. However, there exist few mature and effective ways for managing, transmitting and visualizing massive spatial data. Based on the related work, the paper illustrates the ellipsoidal quadtree technique for rapid access multi-scale and multiple level geographical data, integrates the streaming with level-of-detail rendering method for transmitting the data on the network, and implements large-scale terrain surface simplification using M-band wavelet transforms and multiresolution triangulations. We fulfil a web-based global terrain visual system using COM components on the basis of the above techniques. The system has a good prospect in the military simulation and city plans.

  7. Towards the next generation of global 3D upper mantle Q models

    Gung, Y.; Romanowicz, B.; Capdeville, Y.

    2003-12-01

    Global anelastic tomography can bring important constraints on the thermal structure of the mantle and therefore is dynamics, complementing those provided by elastic tomography. Progress in anelastic tomography has been slow, because of the inherent technical difficulties encountered in discriminating anelastic signal from elastic effects on amplitude data. It has been shown that while the elastic focusing/defocusing effects are not significant at low degrees ( ˜ 8) (e.g. Selby and Woodhouse, 2002; Gung and Romanowicz, 2003), they need to be included to achieve a higher resolution Q model. Ideally, one would use an exact method, such as the Spectral Element Method (SEM) for predicting the focusing effects. SEM is however very heavy computationally. We present a procedure to better constrain the 3D upper mantle Q from 3 component long-period seismic waveforms. In this procedure, the amplitude and phase perturbations due to the 3D elastic structure are corrected for using higher order normal mode asymptotic theory, and applying it to current elastic models. We first evaluate the normal mode asymptotic approach by comparing the corresponding 3D synthetics with those computed using the coupled spectral element/normal mode method (CSEM). 3 normal mode based asymptotic approaches are compared: path average approximation (PAVA), non-linear asymptotic coupling theory (NACT) and NACT+F, an extension of NACT with focusing terms computed using higher order asymptotic theory. Systematic waveform comparison and inversion experiments are implemented. We find that (1) when the anomaly lies on the source-receiver great circle path, the 3 techniques are fairly accurate for fundamental mode surface waves, but NACT and NACT+F provide much better fit for overtone phases and are therefore more powerful in resolving 3D structure in the mid and lower mantle; and (2) the off-great-circle effects, which result in focusing/defocusing and not seen by PAVA or NACT, are well explained by NACT

  8. SALSA3D: Validating a Global 3D P-Velocity Model of the Earth's Crust and Mantle for Improved Event Location

    Begnaud, M. L.; Ballard, S.; Young, C. J.; Hipp, J. R.; Encarnacao, A.; Phillips, W. S.; Chael, E. P.; Rowe, C. A.

    2012-12-01

    We are developing a global 3D P wave velocity model of the Earth's crust and mantle using seismic tomography to assess improvement to seismic event locations obtained using high quality 3D Earth models in lieu of 1D and 2/2.5D models. We present the most recent version of SALSA3D (SAndia LoS Alamos 3D) version 1.9, and demonstrate its ability to reduce mislocations for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth (GT) events. Our model is derived from the latest version of the GT catalog of P/Pn travel-time picks assembled by Los Alamos National Laboratory. For this current version, we employ more robust data quality control measures than previously used, as well as additional global GT data sources. To prevent over-weighting due to ray path redundancy and to reduce the computational burden, we cluster rays into representative rays. The model is represented using the triangular tessellation system described by Ballard et al. (2009), which incorporates variable resolution in both the geographic and radial dimensions. For our starting model, we use a simplified layer crustal model derived from the NNSA Unified model in Eurasia and Crust 2.0 model everywhere else, overlying a uniform ak135 mantle. Sufficient damping is used to reduce velocity adjustments so that ray path changes between iterations are small. We obtain proper model smoothness by using progressive grid refinement, refining the grid only in areas where the data warrant such a refinement. In previous versions, we based this refinement on velocity changes from previous model iterations. For the current version, we utilize the diagonal of the model resolution matrix to control where grid refinement occurs, resulting in more consistent and continuous areas of refinement than before. In addition to the changes in grid refinement, we also employ a more robust convergence criterion between successive grid refinements, allowing a better fit to first broader

  9. Ionospheric energy input as a function of solar wind parameters: global MHD simulation results

    Palmroth, M.; Janhunen, P.; T. I. Pulkkinen; Koskinen, H. E. J.

    2004-01-01

    We examine the global energetics of the solar wind magnetosphere-ionosphere system by using the global MHD simulation code GUMICS-4. We show simulation results for a major magnetospheric storm (6 April 2000) and a moderate substorm (15 August 2001). The ionospheric dissipation is investigated by determining the Joule heating and precipitation powers in the simulation during the two events. The ionospheric dissipation is concentrated largely on the dayside cusp region during ...

  10. Solar wind interaction with Mars' upper atmosphere: Results from 3-D studies using one-way coupling between the Multi-fluid MHD, the M-GITM and the AMPS models

    Dong, C.; Bougher, S. W.; Ma, Y.; Toth, G.; Lee, Y.; Nagy, A. F.; Tenishev, V.; Pawlowski, D. J.; Meng, X.; Combi, M. R.

    2013-12-01

    The study of the solar wind interaction with Mars upper atmosphere/ionosphere has triggered a great of interest in recent years. Among the large number of topics in this research area, the investigation of ion escape fluxes has become increasingly important due to its potential impact on the long-term evolution of Mars atmosphere (e.g., loss of water) over its history. In the present work, we adopt the 3-D Mars cold neutral atmosphere profiles (0~300 km) from the newly developed and validated Mars Global Ionosphere Thermosphere Model (M-GITM) and the 3-D hot oxygen profiles (100km~5RM) from the exosphere Monte Carlo model Adaptive Mesh Particle Simulator (AMPS). We apply these 3-D model outputs fields into the 3-D BATS-R-US Mars multi-fluid MHD model (100km~20RM) that can better simulate the interplay between Mars upper atmosphere and solar wind by considering the dynamics of individual ion species. The multi-fluid model solves separate continuity, momentum and energy equations for each ion species (H+, O+, O2+, CO2+). The M-GITM model together with the AMPS exosphere model take into account the effects of solar cycle and seasonal variations on both cold and hot neutral atmospheres, allowing us to investigate the corresponding effects on the Mars upper atmosphere ion escape by using a one-way coupling approach, i.e., both the M-GITM and AMPS model outputs are used as the inputs for the multi-fluid model and M-GITM is used as input into the AMPS exosphere model. The calculations are carried out for selected cases with different nominal solar wind, solar cycle and crustal field orientation conditions. This work has the potential to provide predictions of ion escape rates for comparison to future data to be returned by the MAVEN primary mission (2014-2016) and thereby improve our understanding of present day escape processes. Acknowledgments: The work presented here was supported by NASA grants NNH10CC04C, NNX09AL26G, NSF grant ATM-0535811.

  11. 3D global estimation and augmented reality visualization of intra-operative X-ray dose.

    Rodas, Nicolas Loy; Padoy, Nicolas

    2014-01-01

    The growing use of image-guided minimally-invasive surgical procedures is confronting clinicians and surgical staff with new radiation exposure risks from X-ray imaging devices. The accurate estimation of intra-operative radiation exposure can increase staff awareness of radiation exposure risks and enable the implementation of well-adapted safety measures. The current surgical practice of wearing a single dosimeter at chest level to measure radiation exposure does not provide a sufficiently accurate estimation of radiation absorption throughout the body. In this paper, we propose an approach that combines data from wireless dosimeters with the simulation of radiation propagation in order to provide a global radiation risk map in the area near the X-ray device. We use a multi-camera RGBD system to obtain a 3D point cloud reconstruction of the room. The positions of the table, C-arm and clinician are then used 1) to simulate the propagation of radiation in a real-world setup and 2) to overlay the resulting 3D risk-map onto the scene in an augmented reality manner. By using real-time wireless dosimeters in our system, we can both calibrate the simulation and validate its accuracy at specific locations in real-time. We demonstrate our system in an operating room equipped with a robotised X-ray imaging device and validate the radiation simulation on several X-ray acquisition setups. PMID:25333145

  12. Ray tracing technique for global 3-D modeling of ionospheric electron density using GNSS measurements

    Alizadeh, Mohamad Mahdi; Schuh, Harald; Schmidt, Michael

    2015-06-01

    For space geodetic techniques, operating in microwave band, ionosphere is a dispersive medium; thus, signals traveling through this medium are in the first approximation, affected proportional to the inverse of the square of their frequencies. This effect allows gaining information about the parameters of the ionosphere in terms of total electron content (TEC) or the electron density (Ne). Making use of this phenomenon, space geodetic techniques have turned into a capable tool for studying the ionosphere in the last decades. Up to now, two-dimensional (2-D) models of Vertical TEC (VTEC) have been widely developed and used by different communities; however, due to the fact that these models provide information about the integral of the whole electron content along the vertical or slant raypath, these maps are not useful when information about the ionosphere at different altitude is required. This paper presents a recent study which aims at developing a global 3-D model of the electron density, using measurements from Global Navigation Satellite Systems and by applying the ray tracing technique to the upper atmosphere. The developed modeling approach represents the horizontal variations of the electron density, with two sets of spherical harmonic expansions of degree and order 15. The height dependency of the electron density is represented by a multilayered Chapman profile function for the bottomside and topside ionosphere, and an appropriate model for the plasmasphere. In addition to the geodetic applications of the developed models, within this study, the 3-D models of electron density can include geophysical parameters like maximum electron density and its corresponding height. High-resolution modeling of these parameters allows an improved geophysical interpretation, which is essential in all studies of the upper atmosphere, space weather, and for the solar-terrestrial environment.

  13. A global 3-D CTM evaluation of black carbon in the Tibetan Plateau

    C. He

    2014-03-01

    Full Text Available We evaluate the black carbon (BC simulations for 2006 over the Tibetan Plateau by a global 3-D chemical transport model using surface observations of BC in surface air and in snow and BC absorption aerosol optical depth (AAOD. Using updated Asian anthropogenic BC emissions (Lu et al., 2011; Zhang et al., 2009 and global biomass burning emissions (Randerson et al., 2012; van der Werf et al., 2010, model results of both surface BC and BC in snow are statistically in good agreement with observations (biases < 15%. Model results capture the seasonal variation of surface BC concentration, but the observed wintertime high values at rural sites in the Indo-Gangetic Plain are absent in the model. Model results are in general agreement with observations (within a factor of two at remote sites. Model simulated BC concentrations in snow are spatiotemporally consistent with observations at most sites. We find that modeled BC AAOD are significantly lower than observations to the northwest of the Plateau and along the southern slopes of the Himalayas during winter and spring, reflecting model deficiencies in emissions, topography and BC mixing state. We find that anthropogenic emissions strongly affect surface BC concentration and AAOD, while the BC aging mainly affects BC in snow over the Plateau.

  14. Extracting 3D Layout From a Single Image Using Global Image Structures

    Z. Lou; T. Gevers; N. Hu

    2015-01-01

    Extracting the pixel-level 3D layout from a single image is important for different applications, such as object localization, image, and video categorization. Traditionally, the 3D layout is derived by solving a pixel-level classification problem. However, the image-level 3D structure can be very b

  15. GLOBAL EXISTENCE AND ASYMPTOTIC BEHAVIOR FOR THE 3D COMPRESSIBLE NON-ISENTROPIC EULER EQUATIONS WITH DAMPING

    张映辉; 吴国春

    2014-01-01

    We investigate the global existence and asymptotic behavior of classical solutions for the 3D compressible non-isentropic damped Euler equations on a periodic domain. The global existence and uniqueness of classical solutions are obtained when the initial data is near an equilibrium. Furthermore, the exponential convergence rates of the pressure and velocity are also proved by delicate energy methods.

  16. Self-organisation in protoplanetary discs. Global, non-stratified Hall-MHD simulations

    Béthune, William; Lesur, Geoffroy; Ferreira, Jonathan

    2016-04-01

    Context. Recent observations have revealed organised structures in protoplanetary discs, such as axisymmetric rings or horseshoe concentrations, evocative of large-scale vortices. These structures are often interpreted as the result of planet-disc interactions. However, these discs are also known to be unstable to the magneto-rotational instability (MRI) which is believed to be one of the dominant angular momentum transport mechanism in these objects. It is therefore natural to ask whether the MRI itself could produce these structures without invoking planets. Aims: The nonlinear evolution of the MRI is strongly affected by the low ionisation fraction in protoplanetary discs. The Hall effect in particular, which is dominant in dense and weakly ionised parts of these objects, has been shown to spontaneously drive self-organising flows in local, shearing box simulations. Here, we investigate the behaviour of global MRI-unstable disc models dominated by the Hall effect and characterise their dynamics. Methods: We validated our implementation of the Hall effect into the PLUTO code with predictions from a spectral method in cylindrical geometry. We then performed 3D unstratified Hall-MHD simulations of Keplerian discs for a broad range of Hall, Ohmic, and ambipolar Elsasser numbers. Results: We confirm the transition from a turbulent to an organised state as the intensity of the Hall effect is increased. We observe the formation of zonal flows, their number depending on the available magnetic flux and on the intensity of the Hall effect. For intermediate Hall intensity, the flow self-organises into long-lived magnetised vortices. Neither the addition of a toroidal field nor Ohmic or ambipolar diffusion change this picture drastically in the range of parameters we have explored. Conclusions: Self-organisation by the Hall effect is a robust phenomenon in global non-stratified simulations. It is able to quench turbulent transport and spontaneously produce axisymmetric

  17. SALSA3D - A Global 3D P-Velocity Model of the Earth's Crust and Mantle for Improved Event Location

    Ballard, S.; Begnaud, M. L.; Young, C. J.; Hipp, J. R.; Chang, M.; Encarnacao, A. V.; Rowe, C. A.; Phillips, W. S.; Steck, L.

    2010-12-01

    To test the hypothesis that high quality 3D Earth models will produce seismic event locations which are more accurate and more precise, we are developing a global 3D P wave velocity model of the Earth’s crust and mantle using seismic tomography. In this paper, we present the most recent version of our model, SALSA3D version 1.5, and demonstrate its ability to reduce mislocations for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth events. Our model is derived from the latest version of the Ground Truth (GT) catalog of P and Pn travel time picks assembled by Los Alamos National Laboratory. To prevent over-weighting due to ray path redundancy and to reduce the computational burden, we cluster rays to produce representative rays. Reduction in the total number of ray paths is ~50%. The model is represented using the triangular tessellation system described by Ballard et al. (2009), which incorporates variable resolution in both the geographic and radial dimensions.. For our starting model, we use a simplified two layer crustal model derived from the Crust 2.0 model over a uniform AK135 mantle. Sufficient damping is used to reduce velocity adjustments so that ray path changes between iterations are small. We obtain proper model smoothness by using progressive grid refinement, refining the grid only around areas with significant velocity changes from the starting model. At each grid refinement level except the last one we limit the number of iterations to prevent convergence thereby preserving aspects of broad features resolved at coarser resolutions. Our approach produces a smooth, multi-resolution model with node density appropriate to both ray coverage and the velocity gradients required by the data. This scheme is computationally expensive, so we use a distributed computing framework based on the Java Parallel Processing Framework, providing us with ~400 processors. Resolution of our model is assessed using a

  18. Towards a Realistic, Data-Driven Thermodynamic MHD Model of the Global Solar Corona

    Downs, Cooper; Roussev, Ilia I.; van der Holst, Bart; Lugaz, Noé; Sokolov, Igor V.; Gombosi, Tamas I.

    2009-01-01

    In this work we describe our implementation of a thermodynamic energy equation into the global corona model of the Space Weather Modeling Framework (SWMF), and its development into the new Lower Corona (LC) model. This work includes the integration of the additional energy transport terms of coronal heating, electron heat conduction, and optically thin radiative cooling into the governing magnetohydrodynamic (MHD) energy equation. We examine two different boundary conditions using this model;...

  19. A hand-held 3D laser scanning with global positioning system of subvoxel precision

    In this paper we propose a hand-held 3D laser scanner composed of an optical head device to extract 3D local surface information and a stereo vision system with subvoxel precision to measure the position and orientation of the 3D optical head. The optical head is manually scanned over the surface object by the operator. The orientation and position of the 3D optical head is determined by a phase-sensitive method using a 2D regular intensity pattern. This phase reference pattern is rigidly fixed to the optical head and allows their 3D location with subvoxel precision in the observation field of the stereo vision system. The 3D resolution achieved by the stereo vision system is about 33 microns at 1.8 m with an observation field of 60cm x 60cm.

  20. A Study of Tropospheric Ozone over China with a 3-D Global CTM Model

    Li Liu

    2007-01-01

    Full Text Available A global 3-D CTM model (OsloCTM2 has been used to study the tropospheric ozone distribution and budget over China. An area covering China and most of East Asia is chosen as the study area. Because of the very nevenly distributed emissions and population in China, the budget study has been done by splitting China into three sub-areas, according to the emission distribution and topography of the country. The model results indicate that in Western China (Area1 dynamic processes are dominating, and the contribution from photochemical ozone production is small. Central and South-East China (Area2 has on average 65% of the photochemical ozone production in China, since more than 80% of the anthropogenic emissions come from this area. Northeast China (Area3 is influenced both by natural and*9nthropogenic emissions. The seasonal variation of ozone budgets was calculated in order to understand how different processes vary with the seasons. The strongest influences of emissions from the continent over the West Pacific region are found in spring, because of the large eastward transport and increased photochemical activities. Most NOx is consumed close to the emission sources; therefore, only 4% of emitted NOx is transported out of China, whereas 70% of the emitted CO is exported. It is calculated that the average net chemical ozone production efficiency by NOx loss is 7.2 in China.

  1. Global small solutions of 2-D incompressible MHD system

    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.

  2. Ladyzhenskaya-Prodi-Serrin type regularity criteria for the 3D incompressible MHD equations in terms of 3 × 3 mixture matrices

    Jia, Xuanji; Zhou, Yong

    2015-09-01

    We prove that a weak solution (u, b) to the MHD equations is smooth on (0, T ] if \\text{M}\\in {{L}α}≤ft(0,T;{{L}γ}≤ft({{{R}}3}\\right)\\right) with 2/α +3/γ =2 , 1≤slant α definition below). As we will explain later, this kind of regularity criteria is more likely to capture the nature of the coupling effects between the fluid velocity and the magnetic field in the evolution of the MHD flows. Moreover, the condition on \\text{M} is scaling invariant, i.e. it is of Ladyzhenskaya-Prodi-Serrin type.

  3. Combined global 2D-local 3D modeling of the industrial Czochralski silicon crystal growth process

    Jung, T.; Seebeck, J.; Friedrich, J.

    2013-04-01

    A global, axisymmetric thermal model of a Czochralski furnace is coupled to an external, local, 3D, time-dependent flow model of the melt via the inclusion of turbulent heat fluxes, extracted from the 3D melt model, into the 2D furnace model. Boundary conditions of the 3D model are updated using results from the 2D model. In the 3D model the boundary layers are resolved by aggressive mesh refinement towards the walls, and the Large Eddy Simulation approach is used to model the turbulent flow in the melt volume on a relatively coarse mesh to minimize calculation times. It is shown that by using this approach it is possible to reproduce fairly good results from Direct Numerical Simulations obtained on much finer meshes, as well as experimental results for interface shape and oxygen concentration in the case of growth of silicon crystals with 210 mm diameter for photovoltaics by the Czochralski method.

  4. Multi-scale self-organisation of edge plasma turbulent transport in 3D global simulations

    Tamain, P.; Ghendrih, Ph; Bufferand, H.; Ciraolo, G.; Colin, C.; Fedorczak, N.; Nace, N.; Schwander, F.; Serre, E.

    2015-05-01

    The 3D global edge turbulence code TOKAM3X is used to study the properties of edge particle turbulent transport in circular limited plasmas, including both closed and open flux surfaces. Turbulence is driven by an incoming particle flux from the core plasma and no scale separation between the equilibrium and the fluctuations is assumed. Simulations show the existence of a complex self-organization of turbulence transport coupling scales ranging from a few Larmor radii up to the machine scale. Particle transport is largely dominated by small scale turbulence with fluctuations forming quasi field-aligned filaments. Radial particle transport is intermittent and associated with the propagation of coherent structures on long distances via avalanches. Long range correlations are also found in the poloidal and toroidal direction. The statistical properties of fluctuations vary with the radial and poloidal directions, with larger fluctuation levels and intermittency found in the outboard scrape-off layer (SOL). Radial turbulent transport is strongly ballooned, with 90% of the flux at the separatrix flowing through the low-field side. One of the main consequences is the existence of quasi-sonic asymmetric parallel flows driving a net rotation of the plasma. Simulations also show the spontaneous onset of an intermittent E × B rotation characterized by a larger shear at the separatrix. Strong correlation is found between the turbulent particle flux and the E × B flow shear in a phenomenology reminiscent of H-mode physics. The poloidal position of the limiter is a key player in the observed dynamics.

  5. Characterizing Vegetation 3D structure Globally using Spaceborne Lidar and Radar.

    Simard, M.; Pinto, N.; Riddick, S.

    2008-12-01

    We characterized global vegetation 3D structure using ICEsat-I/Geoscience Laser Altimeter (GLAS) and improved spatial resolution using ALOS/Phased Array L-band Synthetic Aperture radar (PALSAR) data over 3 sites in the United States. GLAS is a 70m footprint lidar altimeter sampling the ground along-track every 170m with a track separation near the equator around 30km. Forest type classes were initially defined according to the Global Land Cover 2000 map (GLC2000), and 5-degree latitude intervals. This strategy enabled analysis of canopy structure as a function of land cover type and latitude. This produced an irregular grid geographically consistant with GLC2000. To estimate canopy height we removed the ground component from the lidar waveform and computed the centroid of the component due to the forest canopy. Canopy height within a grid cell was produced by computing the weighted mean of the GLAS estimates contained within that cell. The weights were used to reduce the impact of slope on Lidar height estimation errors. Slope is the single most significant source of error when estimating height with a large footprint lidar. It stretches the waveform and causes false estimates of canopy height. The Shuttle Radar Topography Mission (SRTM) elevation data was used to derive slope and weights. Thus, data points located in flat areas were assigned a higher weight than points located in slopes. For each forest type, we modeled the relationship between Lidar-estimated canopy height and five environmental variables: temperature, precipitation, slope, elevation, and anthropogenic disturbance. This ecological model was constructed using the machine learning method Random Forest, due to its flexibility and non-parametric nature. Model accuracy was calculated by subsampling the Lidar data set: using 75% of the data set to produce the map previously described and the remaining 25% for validation. This approach was chosen to characterize individual forest canopy types and their

  6. Local and Global Well-posedness of Strong Solutions to the 3D Primitive Equations with Vertical Eddy Diffusivity

    Cao, Chongsheng; Titi, Edriss S

    2013-01-01

    In this paper, we consider the initial-boundary value problem of the viscous 3D primitive equations for oceanic and atmospheric dynamics with only vertical diffusion in the temperature equation. Local and global well-posedness of strong solutions are established for this system with $H^2$ initial data.

  7. Global Well-posedness of Strong Solutions to the 3D Primitive Equations with Horizontal Eddy Diffusivity

    Cao, Chongsheng; Titi, Edriss S

    2014-01-01

    In this paper, we consider the initial-boundary value problem of the 3D primitive equations for oceanic and atmospheric dynamics with only horizontal diffusion in the temperature equation. Global well-posedness of strong solutions are established with $H^2$ initial data.

  8. Integration of Extended MHD and Kinetic Effects in Global Magnetosphere Models

    Germaschewski, K.; Wang, L.; Maynard, K. R. M.; Raeder, J.; Bhattacharjee, A.

    2015-12-01

    Computational models of Earth's geospace environment are an important tool to investigate the science of the coupled solar-wind -- magnetosphere -- ionosphere system, complementing satellite and ground observations with a global perspective. They are also crucial in understanding and predicting space weather, in particular under extreme conditions. Traditionally, global models have employed the one-fluid MHD approximation, which captures large-scale dynamics quite well. However, in Earth's nearly collisionless plasma environment it breaks down on small scales, where ion and electron dynamics and kinetic effects become important, and greatly change the reconnection dynamics. A number of approaches have recently been taken to advance global modeling, e.g., including multiple ion species, adding Hall physics in a Generalized Ohm's Law, embedding local PIC simulations into a larger fluid domain and also some work on simulating the entire system with hybrid or fully kinetic models, the latter however being to computationally expensive to be run at realistic parameters. We will present an alternate approach, ie., a multi-fluid moment model that is derived rigorously from the Vlasov-Maxwell system. The advantage is that the computational cost remains managable, as we are still solving fluid equations. While the evolution equation for each moment is exact, it depends on the next higher-order moment, so that truncating the hiearchy and closing the system to capture the essential kinetic physics is crucial. We implement 5-moment (density, momentum, scalar pressure) and 10-moment (includes pressure tensor) versions of the model, and use local approximations for the heat flux to close the system. We test these closures by local simulations where we can compare directly to PIC / hybrid codes, and employ them in global simulations using the next-generation OpenGGCM to contrast them to MHD / Hall-MHD results and compare with observations.

  9. A global model simulation for 3-D radiative transfer impact on surface hydrology over Sierra Nevada and Rocky Mountains

    W.-L. Lee

    2014-12-01

    Full Text Available We investigate 3-D mountain effects on solar flux distributions and their impact on surface hydrology over the Western United States, specifically the Rocky Mountains and Sierra Nevada using CCSM4 (CAM4/CLM4 global model with a 0.23° × 0.31° resolution for simulations over 6 years. In 3-D radiative transfer parameterization, we have updated surface topography data from a resolution of 1 km to 90 m to improve parameterization accuracy. In addition, we have also modified the upward-flux deviation [3-D − PP (plane-parallel] adjustment to ensure that energy balance at the surface is conserved in global climate simulations based on 3-D radiation parameterization. We show that deviations of the net surface fluxes are not only affected by 3-D mountains, but also influenced by feedbacks of cloud and snow in association with the long-term simulations. Deviations in sensible heat and surface temperature generally follow the patterns of net surface solar flux. The monthly snow water equivalent (SWE deviations show an increase in lower elevations due to reduced snowmelt, leading to a reduction in cumulative runoff. Over higher elevation areas, negative SWE deviations are found because of increased solar radiation available at the surface. Simulated precipitation increases for lower elevations, while decreases for higher elevations with a minimum in April. Liquid runoff significantly decreases in higher elevations after April due to reduced SWE and precipitation.

  10. 3D Global Coronal Density Structure and Associated Magnetic Field near Solar Maximum

    Kramar, Maxim; Lin, Haosheng

    2016-01-01

    Measurement of the coronal magnetic field is a crucial ingredient in understanding the nature of solar coronal dynamic phenomena at all scales. We employ STEREO/COR1 data obtained near maximum of solar activity in December 2012 (Carrington rotation, CR 2131) to retrieve and analyze the three-dimensional (3D) coronal electron density in the range of heights from $1.5$ to $4\\ \\mathrm{R}_\\odot$ using a tomography method and qualitatively deduce structures of the coronal magnetic field. The 3D electron density analysis is complemented by the 3D STEREO/EUVI emissivity in 195 \\AA \\ band obtained by tomography for the same CR period. We find that the magnetic field configuration during CR 2131 has a tendency to become radially open at heliocentric distances below $\\sim 2.5 \\ \\mathrm{R}_\\odot$. We compared the reconstructed 3D coronal structures over the CR near the solar maximum to the one at deep solar minimum. Results of our 3D density reconstruction will help to constrain solar coronal field models and test the a...

  11. Global MHD modeling of resonant ULF waves: Simulations with and without a plasmasphere

    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.

  12. Radial dependence of solar energetic particles derived from the 15 March 2013 solar energetic particle event and global MHD simulation

    Wu, Chin-Chun; Liou, Kan; Wu, S. T.; Dryer, Murray; Plunkett, Simon

    2016-03-01

    We study an unusual solar energetic particle (SEP) event that was associated with the coronal mass ejection (CME) on March 15, 2013. Enhancements of the SEP fluxes were first detected by the ACE spacecraft at 14:00 UT, ˜7 hours after the onset of the CME (07:00 UT), and the SEP's peak intensities were recorded ˜36 hours after the onset of the CME. Our recent study showed that the CME-driven shock Mach number, based on a global three-dimensional (3-D) magnetohydrodynamic (MHD) simulation, is well correlated with the time-intensity of 10-30 MeV and 30-80 MeV protons. Here we focus on the radial dependence (r-α) of 4He (3.43-41.2 MeV/n) and O (7.30-89.8 MeV/n) energetic particles from ACE/SIS. It is found that the scaling factor (α) ranges between 2 and 4 for most of the energy channels. We also found that the correlation coefficients tend to increase with SEP energies.

  13. Double Dynamo Signatures in a Global MHD Simulation and Mean-field Dynamos

    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.

  14. A global model simulation for 3-D radiative transfer impact on surface hydrology over the Sierra Nevada and Rocky Mountains

    Lee, W.-L.; Gu, Y.; Liou, K. N.; Leung, L. R.; Hsu, H.-H.

    2015-05-01

    We investigate 3-D mountain effects on solar flux distributions and their impact on surface hydrology over the western United States, specifically the Rocky Mountains and the Sierra Nevada, using the global CCSM4 (Community Climate System Model version 4; Community Atmosphere Model/Community Land Model - CAM4/CLM4) with a 0.23° x 0.31° resolution for simulations over 6 years. In a 3-D radiative transfer parameterization, we have updated surface topography data from a resolution of 1 km to 90 m to improve parameterization accuracy. In addition, we have also modified the upward-flux deviation (3-D-PP (plane-parallel)) adjustment to ensure that the energy balance at the surface is conserved in global climate simulations based on 3-D radiation parameterization. We show that deviations in the net surface fluxes are not only affected by 3-D mountains but also influenced by feedbacks of cloud and snow in association with the long-term simulations. Deviations in sensible heat and surface temperature generally follow the patterns of net surface solar flux. The monthly snow water equivalent (SWE) deviations show an increase in lower elevations due to reduced snowmelt, leading to a reduction in cumulative runoff. Over higher-elevation areas, negative SWE deviations are found because of increased solar radiation available at the surface. Simulated precipitation increases for lower elevations, while it decreases for higher elevations, with a minimum in April. Liquid runoff significantly decreases at higher elevations after April due to reduced SWE and precipitation.

  15. Implementation of Inductive Magnetosphere-Ionosphere Coupling and its Effects on Global MHD Magnetospheric Simulations

    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

  16. In-cloud oxalate formation in the global troposphere: a 3-D modeling study

    S. Myriokefalitakis

    2011-01-01

    Full Text Available Organic acids attract increasing attention as contributors to atmospheric acidity, secondary organic aerosol mass and aerosol hygroscopicity. Oxalic acid is globally the most abundant dicarboxylic acid, formed via chemical oxidation of gas-phase precursors in the aqueous phase of aerosols and droplets. Its lifecycle and atmospheric global distribution remain highly uncertain and are the focus of this study. The first global spatial and temporal distribution of oxalate, simulated using a state-of-the-art aqueous phase chemical scheme embedded within the global 3-dimensional chemistry/transport model TM4-ECPL, is here presented. The model accounts for comprehensive gas-phase chemistry and its coupling with major aerosol constituents (including secondary organic aerosol. Model results are consistent with ambient observations of oxalate at rural and remote locations (slope = 0.83 ± 0.06, r2 = 0.67, N = 106 and suggest that aqueous phase chemistry contributes significantly to the global atmospheric burden of secondary organic aerosol. In TM4-ECPL most oxalate is formed in-clouds and less than 10% is produced in aerosol water. About 61% of the oxalate is removed via wet deposition, 35% by in-cloud reaction with hydroxyl radical and 4% by dry deposition. The global oxalate net chemical production is calculated to be about 17–27 Tg yr−1 with almost 91% originating from biogenic hydrocarbons, mainly isoprene. This condensed phase net source of oxalate in conjunction with a global mean turnover time against deposition of about 5 days, maintain oxalate's global tropospheric burden of 0.24–0.39 Tg that is about 13–19% of calculated total organic aerosol burden.

  17. In-cloud oxalate formation in the global troposphere: a 3-D modeling study

    S. Myriokefalitakis

    2011-06-01

    Full Text Available Organic acids attract increasing attention as contributors to atmospheric acidity, secondary organic aerosol mass and aerosol hygroscopicity. Oxalic acid is globally the most abundant dicarboxylic acid, formed via chemical oxidation of gas-phase precursors in the aqueous phase of aerosols and droplets. Its lifecycle and atmospheric global distribution remain highly uncertain and are the focus of this study. The first global spatial and temporal distribution of oxalate, simulated using a state-of-the-art aqueous-phase chemical scheme embedded within the global 3-dimensional chemistry/transport model TM4-ECPL, is here presented. The model accounts for comprehensive gas-phase chemistry and its coupling with major aerosol constituents (including secondary organic aerosol. Model results are consistent with ambient observations of oxalate at rural and remote locations (slope = 1.16 ± 0.14, r2 = 0.36, N = 114 and suggest that aqueous-phase chemistry contributes significantly to the global atmospheric burden of secondary organic aerosol. In TM4-ECPL most oxalate is formed in-cloud and less than 5 % is produced in aerosol water. About 62 % of the oxalate is removed via wet deposition, 30 % by in-cloud reaction with hydroxyl radical, 4 % by in-cloud reaction with nitrate radical and 4 % by dry deposition. The in-cloud global oxalate net chemical production is calculated to be about 21–37 Tg yr−1 with almost 79 % originating from biogenic hydrocarbons, mainly isoprene. This condensed phase net source of oxalate in conjunction with a global mean turnover time against deposition of about 5 days, maintain oxalate's global tropospheric burden of 0.2–0.3 Tg, i.e. 0.05–0.1 Tg-C that is about 5–9 % of model-calculated water soluble organic carbon burden.

  18. Validation of "AW3D" Global Dsm Generated from Alos Prism

    Takaku, Junichi; Tadono, Takeo; Tsutsui, Ken; Ichikawa, Mayumi

    2016-06-01

    Panchromatic Remote-sensing Instrument for Stereo Mapping (PRISM), one of onboard sensors carried by Advanced Land Observing Satellite (ALOS), was designed to generate worldwide topographic data with its optical stereoscopic observation. It has an exclusive ability to perform a triplet stereo observation which views forward, nadir, and backward along the satellite track in 2.5 m ground resolution, and collected its derived images all over the world during the mission life of the satellite from 2006 through 2011. A new project, which generates global elevation datasets with the image archives, was started in 2014. The data is processed in unprecedented 5 m grid spacing utilizing the original triplet stereo images in 2.5 m resolution. As the number of processed data is growing steadily so that the global land areas are almost covered, a trend of global data qualities became apparent. This paper reports on up-to-date results of the validations for the accuracy of data products as well as the status of data coverage in global areas. The accuracies and error characteristics of datasets are analyzed by the comparison with existing global datasets such as Ice, Cloud, and land Elevation Satellite (ICESat) data, as well as ground control points (GCPs) and the reference Digital Elevation Model (DEM) derived from the airborne Light Detection and Ranging (LiDAR).

  19. Dayside Proton Aurora: Comparisons between Global MHD Simulations and Image Observations

    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

  20. Investigating Magnetic Activity in the Galactic Centre by Global MHD Simulation

    Suzuki, Takeru K; Torii, Kazufumi; Machida, Mami; Matsumoto, Ryoji; Kakiuchi, Kensuke

    2016-01-01

    By performing a global magnetohydrodynamical (MHD) simulation for the Milky Way with an axisymmetric gravitational potential, we propose that spatially dependent amplification of magnetic fields possibly explains the observed noncircular motion of the gas in the Galactic centre (GC) region. The radial distribution of the rotation frequency in the bulge region is not monotonic in general. The amplification of the magnetic field is enhanced in regions with stronger differential rotation, because magnetorotational instability and field-line stretching are more effective. The strength of the amplified magnetic field reaches >~ 0.5 mG, and radial flows of the gas are excited by the inhomogeneous transport of angular momentum through turbulent magnetic field that is amplified in a spatially dependent manner. As a result, the simulated position-velocity diagram exhibits a time-dependent asymmetric parallelogram-shape owing to the intermittency of the magnetic turbulence; the present model provides a viable alternati...

  1. Modeling Recent Large Earthquakes Using the 3-D Global Wave Field

    Hjörleifsdóttir, V.; Kanamori, H.; Tromp, J.

    2003-04-01

    We use the spectral-element method (SEM) to accurately compute waveforms at periods of 40 s and longer for three recent large earthquakes using 3D Earth models and finite source models. The M_w~7.6, Jan~26, 2001, Bhuj, India event had a small rupture area and is well modeled at long periods with a point source. We use this event as a calibration event to investigate the effects of 3-D Earth models on the waveforms. The M_w~7.9, Nov~11, 2001, Kunlun, China, event exhibits a large directivity (an asymmetry in the radiation pattern) even at periods longer than 200~s. We used the source time function determined by Kikuchi and Yamanaka (2001) and the overall pattern of slip distribution determined by Lin et al. to guide the wave-form modeling. The large directivity is consistent with a long fault, at least 300 km, and an average rupture speed of 3±0.3~km/s. The directivity at long periods is not sensitive to variations in the rupture speed along strike as long as the average rupture speed is constant. Thus, local variations in rupture speed cannot be ruled out. The rupture speed is a key parameter for estimating the fracture energy of earthquakes. The M_w~8.1, March~25, 1998, event near the Balleny Islands on the Antarctic Plate exhibits large directivity in long period surface waves, similar to the Kunlun event. Many slip models have been obtained from body waves for this earthquake (Kuge et al. (1999), Nettles et al. (1999), Antolik et al. (2000), Henry et al. (2000) and Tsuboi et al. (2000)). We used the slip model from Henry et al. to compute SEM waveforms for this event. The synthetic waveforms show a good fit to the data at periods from 40-200~s, but the amplitude and directivity at longer periods are significantly smaller than observed. Henry et al. suggest that this event comprised two subevents with one triggering the other at a distance of 100 km. To explain the observed directivity however, a significant amount of slip is required between the two subevents

  2. Global Solutions of the Equations of 3D Compressible Magnetohydrodynamics with Zero Resistivity

    Suen, Anthony

    2012-01-01

    We prove the global-in-time existence of H^2 solutions of the equations of compressible magnetohydrodynamics with zero magnetic resistivity in three space dimensions. Initial data are taken to be small in H^2 modulo a constant state and initial densities are positive and essentially bounded. The present work generalizes the results obtained by Kawashima.

  3. Global weak solutions to the inviscid 3D Quasi-Geostrophic equation

    Puel, Marjolaine; Vasseur, Alexis F.

    2014-01-01

    In this article, the authors prove the existence of global weak solutions to the inviscid three-dimensional quasi-geostrophic equation. This equation models the evolution of the temperature on the surface of the earth. It is widely used in geophysics and meteorology.

  4. Propagation of BBFs and Dipolarization Fronts in the Global MHD simulation of February 27, 2009 Substorm

    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.

  5. On the influence of CMEs on the global 3-D coronal electron density

    M. Kramar

    2011-06-01

    Full Text Available In order to analize the influence of a Coronal Mass Ejection (CME on the coronal streamer belt, we made 3-D reconstructions of the electron density in the corona at heliospheric distances from 1.5 to 4 R for periods before and after a CME occured. The reconstructions were performed using a tomography technique. We studied two CME cases: (i a slow CME on 1 June 2008; (ii two fast CMEs on 31 December 2007 and 2 January 2008. For the first case of slow CME, it was found: (i the potential magnetic field configuration in the CME initiation region before the CME does not agree with the coronal density structure while after the CME the agreement between the field and density is much better. This could be manifistation of that that the field was non-potential before the CME and after the CME the field relaxes towards a more potential state. (ii It was shown that the dimming caused by the slow CME is not due to rotation of the corona and a line-of-sight (LOS effect but a streamer blow out effect took place.

  6. Connecting Global Measures of 3D Magnetic Reconnection to Local Kinetic Physics

    Daughton, William Scott [Los Alamos National Laboratory

    2015-07-16

    After giving the motivation for the work, slides present the topic under the following headings: Description of LAPD experiment; Actual simulation setup; Simple kinetic theory of ined-tied tearing; Diagnostics to characterizing 3D reconnection; Example #1 - short-tied system; and Example #2 - long line-tied system. Colorful simulations are shown for quasipotential vs field line exponentiation, field line integrated Ohms Law, and correlation with agyrotopy & energy conversion for example #1; and evolution of current density for largest case, field exponentiation vs quasi-potential, and time evolution of magnetic field lines for example #2. To satisfy line-tied boundary conditions, there is need for superposition of oblique modes--the simple two-mode approximation works surprisingly well. For force-free layers with bg >1, the fastest growing periodic modes are oblique with kxλ ~0.5. This implies a minimum length of Ly > 2πλbg. There are strong correlations between σ → Ξ → A0e (observable with spacecraft). Electron pressure tensor is the dominant non-ideal term.

  7. Ideal and non-ideal MHD regimes of wire array implosion obtained in 3D hybrid simulations and observed during experiments at NTF (Nevada Terawatt Facility)

    Recent 3D hybrid simulation of a plasma current-carrying column revealed two regimes of sausage and kink instability development. In the first regime, with small Hall parameter, development of instabilities leads to appearance of large-scale axial perturbations and eventually to the bending of the plasma column. In the second regime, with five times larger Hall parameter, small-scale perturbations dominated and no bending of the plasma column was observed. Simulation results are compared to recent experimental data, including laser probing, x-ray spectroscopy and time-gated x-ray imaging during wire array implosions at NTF

  8. Higher spin gravity in 3D: Black holes, global charges and thermodynamics

    Pérez, Alfredo, E-mail: aperez@cecs.cl [Centro de Estudios Científicos (CECs), Casilla 1469, Valdivia (Chile); Tempo, David, E-mail: tempo@cecs.cl [Centro de Estudios Científicos (CECs), Casilla 1469, Valdivia (Chile); Troncoso, Ricardo, E-mail: troncoso@cecs.cl [Centro de Estudios Científicos (CECs), Casilla 1469, Valdivia (Chile); Universidad Andrés Bello, Av. República 440, Santiago (Chile)

    2013-10-07

    Global charges and thermodynamic properties of three-dimensional higher spin black holes that have been recently found in the literature are revisited. Since these solutions possess a relaxed asymptotically AdS behavior, following the canonical approach, it is shown that the global charges, and in particular the energy, acquire explicit nontrivial contributions given by nonlinear terms in the deviations with respect to the reference background. It is also found that there are cases for which the first law of thermodynamics can be readily worked out in the canonical ensemble, i.e., without work terms associated to the presence of higher spin fields, and remarkably, the semiclassical higher spin black hole entropy is exactly reproduced from Cardy formula.

  9. Sufficiency Class for Global (in Time) Solutions to the 3D-Navier-Stokes Equations

    Gill, Tepper L.; Zachary, Woodford W

    2007-01-01

    A well-known unsolved problem (in the classical theory of fluid mechanics) is to identify a set of initial velocities, which may depend on the viscosity, the body forces and possibly the boundary of the fluid that will allow global in time solutions to the three-dimensional Navier-Stokes equations. (These equations describe the time evolution of the fluid velocity and pressure of an incompressible viscous homogeneous Newtonian fluid in terms of a given initial velocity and given external body...

  10. Fast and accurate global multiphase arrival tracking: the irregular shortest-path method in a 3-D spherical earth model

    Huang, Guo-Jiao; Bai, Chao-Ying; Greenhalgh, Stewart

    2013-09-01

    The traditional grid/cell-based wavefront expansion algorithms, such as the shortest path algorithm, can only find the first arrivals or multiply reflected (or mode converted) waves transmitted from subsurface interfaces, but cannot calculate the other later reflections/conversions having a minimax time path. In order to overcome the above limitations, we introduce the concept of a stationary minimax time path of Fermat's Principle into the multistage irregular shortest path method. Here we extend it from Cartesian coordinates for a flat earth model to global ray tracing of multiple phases in a 3-D complex spherical earth model. The ray tracing results for 49 different kinds of crustal, mantle and core phases show that the maximum absolute traveltime error is less than 0.12 s and the average absolute traveltime error is within 0.09 s when compared with the AK135 theoretical traveltime tables for a 1-D reference model. Numerical tests in terms of computational accuracy and CPU time consumption indicate that the new scheme is an accurate, efficient and a practical way to perform 3-D multiphase arrival tracking in regional or global traveltime tomography.

  11. Characterisation of global flow and local fluctuations in 3D SPH simulations of protoplanetary discs

    Arena, Serena

    2013-01-01

    A complete and detailed knowledge of the structure of the gaseous component in protoplanetary discs is essential to the study of dust evolution during the early phases of pre-planetesimal formation. The aim of this paper is to determine if three-dimensional accretion discs simulated by the Smoothed Particle Hydrodynamics (SPH) method can reproduce the observational data now available and the expected turbulent nature of protoplanetary discs. The investigation is carried out by setting up a suite of diagnostic tools specifically designed to characterise both the global flow and the fluctuations of the gaseous disc. The main result concerns the role of the artificial viscosity implementation in the SPH method: in addition to the already known ability of SPH artificial viscosity to mimic a physical-like viscosity under specific conditions, we show how the same artificial viscosity prescription behaves like an implicit turbulence model. In fact, we identify a threshold for the parameters in the standard artificia...

  12. Representing Objects using Global 3D Relational Features for Recognition Tasks

    Mustafa, Wail

    2015-01-01

    In robotic systems, visual interpretations of the environment compose an essential element in a variety of applications, especially those involving manipulation of objects. Interpreting the environment is often done in terms of recognition of objects using machine learning approaches. For user...... representations. For representing objects, we derive global descriptors encoding shape using viewpoint-invariant features obtained from multiple sensors observing the scene. Objects are also described using color independently. This allows for combining color and shape when it is required for the task. For more...... to initiate higher-level semantic interpretations of complex scenes. In the object category recognition task, we present a system that is capable of assigning multiple and nested categories for novel objects using a method developed for this purpose. Integrating this method with other multi-label learning...

  13. Model studies with a 3-d global CTM using ECMWF data

    Sundet, Jostein Kandal

    1997-12-31

    It has become evident during the last decades that the large-scale tropospheric distribution of chemical compounds is changing due to man-made emissions. These changes, in turn, lead to changes in the atmospheric oxidation capacity that may significantly affect the environment, such as global warming and increased acid rain. The main objective of this thesis was to use data from a meteorological weather prediction model to drive a global Chemical Tracer (or transport) Model (CTM), the Oslo CTM-2, where the data are generated at the European Centre for Medium-Range Weather Forecasts (ECMWF). It is the kind of data that are taken from the IFS model that is the innovative part of this study. In contrast to other studies, the purpose of the IFS model run is to generate selected data instead of using data from the ECMWF archive. This enabled the extraction of fields useful for atmospheric chemistry studies. The thesis describes the CTM-2 model, presents the data from the IFS model and briefly introduces tropospheric chemistry. Experiments are run to verify the model. The following points for improvements of the model are pointed out: (1) clouds must be included in the photodissociation calculations, (2) better boundary value representing the stratosphere must be obtained for ozone and probably NO and HNO{sub 3} as well, (3) the surface deposition scheme for ozone must be improved, and (4) fluctuations in the chemistry were found and should be eliminated through assessing the fast cycling of OH and HO{sub 2}. Some issues for consideration by the ECMWF are also pointed out. 57 refs., 24 figs., 16 tabs.

  14. Compressing an Ensemble with Statistical Models: An Algorithm for Global 3D Spatio-Temporal Temperature

    Castruccio, Stefano

    2015-04-02

    One of the main challenges when working with modern climate model ensembles is the increasingly larger size of the data produced, and the consequent difficulty in storing large amounts of spatio-temporally resolved information. Many compression algorithms can be used to mitigate this problem, but since they are designed to compress generic scientific data sets, they do not account for the nature of climate model output and they compress only individual simulations. In this work, we propose a different, statistics-based approach that explicitly accounts for the space-time dependence of the data for annual global three-dimensional temperature fields in an initial condition ensemble. The set of estimated parameters is small (compared to the data size) and can be regarded as a summary of the essential structure of the ensemble output; therefore, it can be used to instantaneously reproduce the temperature fields in an ensemble with a substantial saving in storage and time. The statistical model exploits the gridded geometry of the data and parallelization across processors. It is therefore computationally convenient and allows to fit a non-trivial model to a data set of one billion data points with a covariance matrix comprising of 10^18 entries.

  15. Ionospheric energy input as a function of solar wind parameters: global MHD simulation results

    M. Palmroth

    2004-01-01

    Full Text Available We examine the global energetics of the solar wind magnetosphere-ionosphere system by using the global MHD simulation code GUMICS-4. We show simulation results for a major magnetospheric storm (6 April 2000 and a moderate substorm (15 August 2001. The ionospheric dissipation is investigated by determining the Joule heating and precipitation powers in the simulation during the two events. The ionospheric dissipation is concentrated largely on the dayside cusp region during the main phase of the storm period, whereas the nightside oval dominates the ionospheric dissipation during the substorm event. The temporal variations of the precipitation power during the two events are shown to correlate well with the commonly used AE-based proxy of the precipitation power. The temporal variation of the Joule heating power during the substorm event is well-correlated with a commonly used AE-based empirical proxy, whereas during the storm period the simulated Joule heating is different from the empirical proxy. Finally, we derive a power law formula, which gives the total ionospheric dissipation from the solar wind density, velocity and magnetic field z-component and which agrees with the simulation result with more than 80% correlation.

    Key words. Ionosphere (modeling and forecasting – Magnetospheric physics (magnetosphere-ionosphere interactions; storms and substorms

  16. Global MHD Modelling of the ISM - From large towards small scale turbulence

    D'Avillez, M A; Avillez, Miguel A. de; Breitschwerdt, Dieter

    2005-01-01

    Dealing numerically with the turbulent nature and non-linearity of the physical processes involved in the ISM requires the use of sophisticated numerical schemes coupled to HD and MHD mathematical models. SNe are the main drivers of the interstellar turbulence by transferring kinetic energy into the system. This energy is dissipated by shocks (which is more efficient) and by molecular viscosity. We carried out adaptive mesh refinement simulations (with a finest resolution of 0.625 pc) of the turbulent ISM embedded in a magnetic field with mean field components of 2 and 3 $\\mu$G. The time scale of our run was 400 Myr, sufficiently long to avoid memory effects of the initial setup, and to allow for a global dynamical equilibrium to be reached in case of a constant energy input rate. It is found that the longitudinal and transverse turbulent length scales have a time averaged (over a period of 50 Myr) ratio of 0.52-0.6, almost similar to the one expected for isotropic homogeneous turbulence. The mean characteris...

  17. Alpha particle effects on global MHD modes, and alpha particle transport in ignited tokamaks

    The high frequency, low mode number toroidicity-induced Alfven eigenmodes (TAE) are shown to be driven unstable primarily by the circulating α-particles through wave-particle resonances. To destabilize the TAE modes, the inverse Landau damping associated with the α-particle pressure gradient free energy must overcome the velocity space Landau damping due to both the α-particles and the core electrons and ions, as well as Alfven continuum damping. Stability criteria are presented for TFTR, CIT, and ITER tokamaks in terms of the α-particle beta βα, the α-particle pressure gradient parameter (ω*/ωA), where ω* is the α-particle diamagnetic drift frequency, and the α-particle velocity (vα/vA) parameter. Typically the volume averaged α-particle beta threshold is on the order of 10-4. Rough estimates of the TAE mode saturation level give δBr/B ∼ 10-3 for typical D-T tokamak operations. Significant α-particle losses are found when the amplitude of the global MHD modes is large, on the order of (δBr/B) ≥ 10-4. For (δBr/B) = 5 x 10-4, the α-particle loss time is appreciably shorter than the α-particle slowing-down time. 13 refs., 1 fig

  18. Solar wind-magnetosphere energy coupling function fitting: Results from a global MHD simulation

    Wang, C.; Han, J. P.; Li, H.; Peng, Z.; Richardson, J. D.

    2014-08-01

    Quantitatively estimating the energy input from the solar wind into the magnetosphere on a global scale is still an observational challenge. We perform three-dimensional magnetohydrodynamic (MHD) simulations to derive the energy coupling function. Based on 240 numerical test runs, the energy coupling function is given by Ein=3.78×107nsw0.24Vsw1.47BT0.86[sin2.70(θ/2)+0.25]. We study the correlations between the energy coupling function and a wide variety of magnetospheric activity, such as the indices of Dst, Kp, ap, AE, AU, AL, the polar cap index, and the hemispheric auroral power. The results indicate that this energy coupling function gives better correlations than the ɛ function. This result is also applied to a storm event under northward interplanetary magnetic field conditions. About 13% of the solar wind kinetic energy is transferred into the magnetosphere and about 35% of the input energy is dissipated in the ionosphere, consistent with previous studies.

  19. The Influence of Pickup Protons, from Interstellar Neutral Hydrogen, on the Propagation of Interplanetary Shocks from the Halloween 2003 Solar Events to ACE and Ulysses: A 3-D MHD Modeling Study

    Detman, T. R.; Intriligator, D. S.; Dryer, M.; Sun, W.; Deehr, C. S.; Intriligator, J.

    2012-01-01

    We describe our 3-D, time ]dependent, MHD solar wind model that we recently modified to include the physics of pickup protons from interstellar neutral hydrogen. The model has a time-dependent lower boundary condition, at 0.1 AU, that is driven by source surface map files through an empirical interface module. We describe the empirical interface and its parameter tuning to maximize model agreement with background (quiet) solar wind observations at ACE. We then give results of a simulation study of the famous Halloween 2003 series of solar events. We began with shock inputs from the Fearless Forecast real ]time shock arrival prediction study, and then we iteratively adjusted input shock speeds to obtain agreement between observed and simulated shock arrival times at ACE. We then extended the model grid to 5.5 AU and compared those simulation results with Ulysses observations at 5.2 AU. Next we undertook the more difficult tuning of shock speeds and locations to get matching shock arrival times at both ACE and Ulysses. Then we ran this last case again with neutral hydrogen density set to zero, to identify the effect of pickup ions. We show that the speed of interplanetary shocks propagating from the Sun to Ulysses is reduced by the effects of pickup protons. We plan to make further improvements to the model as we continue our benchmarking process to 10 AU, comparing our results with Cassini observations, and eventually on to 100 AU, comparing our results with Voyager 1 and 2 observations.

  20. Intercomparison of modal and sectional aerosol microphysics representations within the same 3-D global chemical transport model

    G. W. Mann

    2012-05-01

    Full Text Available In the most advanced aerosol-climate models it is common to represent the aerosol particle size distribution in terms of several log-normal modes. This approach, motivated by computational efficiency, makes assumptions about the shape of the particle distribution that may not always capture the properties of global aerosol. Here, a global modal aerosol microphysics module (GLOMAP-mode is evaluated and improved by comparing against a sectional version (GLOMAP-bin and observations in the same 3-D global offline chemistry transport model. With both schemes, the model captures the main features of the global particle size distribution, with sub-micron aerosol approximately unimodal in continental regions and bi-modal in marine regions. Initial bin-mode comparisons showed that the current values for two size distribution parameter settings in the modal scheme (mode widths and inter-modal separation sizes resulted in clear biases compared to the sectional scheme. By adjusting these parameters in the modal scheme, much better agreement is achieved against the bin scheme and observations. Annual mean surface-level mass of sulphate, sea-salt, black carbon (BC and organic carbon (OC are within 25% in the two schemes in nearly all regions. Surface level concentrations of condensation nuclei (CN, cloud condensation nuclei (CCN, surface area density and condensation sink also compare within 25% in most regions. However, marine CCN concentrations between 30° N and 30° S are systematically 25–60% higher in the modal model, which we attribute to differences in size-resolved particle growth or cloud-processing. Larger differences also exist in regions or seasons dominated by biomass burning and in free-troposphere and high-latitude regions. Indeed, in the free-troposphere, GLOMAP-mode BC is a factor 2–4 higher than GLOMAP-bin, likely due to differences in size-resolved scavenging. Nevertheless, in most parts of the atmosphere, we conclude that bin

  1. Intercomparison of modal and sectional aerosol microphysics representations within the same 3-D global chemical transport model

    G. W. Mann

    2012-01-01

    Full Text Available A global modal aerosol microphysics module (GLOMAP-mode is evaluated and improved by comparing against a sectional version (GLOMAP-bin and observations in the same 3-D global offline chemistry transport model. With both schemes, the model captures the main features of the global particle size distribution, with sub-micron aerosol approximately unimodal in continental regions and bi-modal in marine regions. Initial bin-mode comparisons showed that various size distribution parameter settings (mode widths and inter-modal separation sizes resulted in clear biases compared to the sectional scheme. By adjusting these parameters in the modal scheme, much better agreement is achieved against the bin scheme and observations. Surface mass of sulphate, sea-salt, black carbon (BC and organic carbon (OC are, on the annual mean, within 25 % in the two schemes in nearly all regions. On the annual mean, surface level concentrations of condensation nuclei (CN, cloud condensation nuclei (CCN, surface area density and condensation sink also compare within 25 % in most regions. However, marine CCN concentrations between 30° N and 30° S are systematically higher in the modal scheme, by 25–60 %, which we attribute to differences in size-resolved particle growth or cloud-processing. Larger differences also exist in regions or seasons dominated by biomass burning and in free-troposphere and high-latitude regions. Indeed, in the free-troposphere, GLOMAP-mode BC is a factor 2–4 higher than GLOMAP-bin, likely due to differences in size-resolved scavenging. Nevertheless, in most parts of the atmosphere, we conclude that bin-mode differences are much less than model-observation differences, although some processes are missing in these runs which may pose a bigger challenge to modal schemes (e.g. boundary layer nucleation, ultra-fine sea-spray. The findings here underline the need for a spectrum of complexity in global models, with size-resolved aerosol properties

  2. 3D MHD Simulations of Spheromak Compression

    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.

  3. Global existence and uniqueness theorem to 2-D incompressible non-resistive MHD system with non-equilibrium background magnetic field

    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.

  4. Three-dimensional Kinetic-MHD Model of the Global Heliosphere with the Heliopause-surface Fitting

    Izmodenov, V. V.; Alexashov, D. B.

    2015-10-01

    This paper provides a detailed description of the latest version of our model of the solar wind (SW) interaction with the local interstellar medium (LISM). This model has already been applied to the analysis of Lyα absorption spectra toward nearby stars and for analyses of Solar and Heliospheric Observatory/SWAN data. Katushkina et al. (this issue) used the model results to analyze IBEX-Lo data. At the same time, the details of this model have not yet been published. This is a three-dimensional (3D) kinetic-magnetohydrodynamical (MHD) model that takes into account SW and interstellar plasmas (including α particles in SW and helium ions in LISM), the solar and interstellar magnetic fields, and interstellar hydrogen atoms. The latitudinal dependence of SW and the actual flow direction of the interstellar gas with respect to the Sun are also taken into account in the model. It was very essential that our numerical code was developed in such a way that any numerical diffusion or reconnection across the heliopause were not allowed in the model. The heliospheric current sheet is a rotational discontinuity in the ideal MHD and can be treated kinematically. In the paper, we focus in particular on the effects of the heliospheric magnetic field and on the heliolatitudinal dependence of SW.

  5. Assessing the performance of community-available global MHD models using key system parameters and empirical relationships

    Gordeev, E.; Sergeev, V.; Honkonen, I.; Kuznetsova, M.; Rastätter, L.; Palmroth, M.; Janhunen, P.; Tóth, G.; Lyon, J.; Wiltberger, M.

    2015-12-01

    Global magnetohydrodynamic (MHD) modeling is a powerful tool in space weather research and predictions. There are several advanced and still developing global MHD (GMHD) models that are publicly available via Community Coordinated Modeling Center's (CCMC) Run on Request system, which allows the users to simulate the magnetospheric response to different solar wind conditions including extraordinary events, like geomagnetic storms. Systematic validation of GMHD models against observations still continues to be a challenge, as well as comparative benchmarking of different models against each other. In this paper we describe and test a new approach in which (i) a set of critical large-scale system parameters is explored/tested, which are produced by (ii) specially designed set of computer runs to simulate realistic statistical distributions of critical solar wind parameters and are compared to (iii) observation-based empirical relationships for these parameters. Being tested in approximately similar conditions (similar inputs, comparable grid resolution, etc.), the four models publicly available at the CCMC predict rather well the absolute values and variations of those key parameters (magnetospheric size, magnetic field, and pressure) which are directly related to the large-scale magnetospheric equilibrium in the outer magnetosphere, for which the MHD is supposed to be a valid approach. At the same time, the models have systematic differences in other parameters, being especially different in predicting the global convection rate, total field-aligned current, and magnetic flux loading into the magnetotail after the north-south interplanetary magnetic field turning. According to validation results, none of the models emerges as an absolute leader. The new approach suggested for the evaluation of the models performance against reality may be used by model users while planning their investigations, as well as by model developers and those interesting to quantitatively

  6. Modelling of Criegee Intermediates using the 3-D global model, STOCHEM-CRI and investigating their global impacts on Secondary Organic Aerosol formation

    Khan, M. Anwar H.; Cooke, Michael; Utembe, Steve; Archibald, Alexander; Derwent, Richard; Jenkin, Mike; Lyons, Kyle; Kent, Adam; Percival, Carl; Shallcross, Dudley E.

    2016-04-01

    Gas phase reactions of ozone with unsaturated compounds form stabilized Criegee intermediates (sCI) which play an important role in controlling the budgets of many tropospheric species including OH, organic acids and secondary organic aerosols (SOA). Recently sCI has been proposed to play a significant role in atmospheric sulfate and nitrate chemistry by forming sulfuric acid (promoter of aerosol formation) and nitrate radical (a powerful oxidizing agent). sCI can also undergo association reactions with water, alcohols, and carboxylic acids to form hydroperoxides and with aldehydes and ketones to form secondary ozonides. The products from these reactions are low volatility compounds which can contribute to the formation of SOA. The importance of plant emitted alkenes (isoprene, monoterpenes, sesquiterpenes) in the production of SOA through sCI formation have already been investigated in laboratory studies. However, the SOA formation from these reactions are absent in current global models. Thus, the formation of SOA has been incorporated in the global model, STOCHEM-CRI, a 3-D global chemistry transport model and the role of CI chemistry in controlling atmospheric composition and climate, and the influence of water vapor has been discussed in the study.

  7. GEM/POPs: a global 3-D dynamic model for semi-volatile persistent organic pollutants – Part 2: Global transports and budgets of PCBs

    L. A. Barrie

    2007-03-01

    Full Text Available Global transports and budgets of three PCBs were investigated with a 3-D dynamic model for semi-volatile persistent organic pollutants – GEM/POPs. Dominant pathways were identified for PCB transports in the atmosphere with a peak transport flux below 8 km and 14 km for gaseous and particulate PCB28, 4 km and 6 km for gaseous and particulate PCB180. The inter-continental transports of PCBs in the Northern Hemisphere (NH are dominated in the zonal direction with their route changes seasonally regulated by the variation of westerly jet. The transport pathways from Europe and North Atlantic to the Arctic contributed the most PCBs over there. Inter-hemispheric transports of PCBs originated from the regions of Europe, Asia and North America in three different flow-paths, accompanying with easterly jet, Asian monsoon winds and trade winds. PCBs from the Southern Hemisphere (SH could export into the NH. According to the PCB emissions of year 2000, Europe, North America and Asia are the three largest sources of the three PCBs, contributing to the global background concentrations in the atmosphere and soil and water. Globally, PCB28 in soil and water has become a comparable source to the anthropogenic emissions while heavier PCBs such as PCB153 and 180 are still transporting into soil and water. It is found that lighter PCBs have more long range transport potentials than their heavier counter-parts in the atmosphere.

  8. Tail current surge: New insights from a global MHD simulation and comparison with satellite observations

    Ohtani, Shin-Ichi; Raeder, Joachim

    2004-01-01

    The present study examines the tailward propagation of substorm-associated variations of the tail current intensity. In the substorm event of 24 November 1996, the Interball and IMP 8 satellites were located in the midnight sector at X = -26 and -36 RE, respectively, and observed an increase and a decrease of the lobe magnetic field strength corresponding to the storage and release of the lobe magnetic energy. Both spacecraft observed BZ to decrease initially and then increase in the course of the decrease in ∣BX∣, a feature that was reported previously as a manifestation of the tailward expansion of the current disruption region. The delay of the signatures between the two satellites confirms that the associated current system moved tailward. Motivated by this fortuitous coordination of the satellite observation, the present study revisits a global MHD simulation previously conducted specifically for this substorm event [, 2001]. The most noticeable feature of the modeled tail dynamics is the repeated occurrence of tail current surges, that is, temporal intensifications of the tail current that propagate tailward. The first tail current surge is accompanied by the stretching of the tail magnetic field, which starts in the inner magnetosphere and extends tailward. The associated tailward flow redistributes the plasma pressure in such a way that the tail current is reduced in its intensity in the near-Earth region, while the pressure gradient increases at the propagation front, which intensifies the local current. The last major tail current surge is caused by the near-Earth reconnection. Inside a plasmoid, the pressure gradient current is intensified on the tailward side of the O-line, and it propagates tailward as the plasmoid grows and is released. For each tail current surge, irrespective of its cause, the intensification of the tail current is followed by the reduction, and its tailward propagation creates the aforementioned phase relationship between BX

  9. Changes in dissolved iron deposition to the oceans driven by human activity: a 3-D global modelling study

    Myriokefalitakis, S.; Daskalakis, N.; Mihalopoulos, N.; Baker, A. R.; Nenes, A.; Kanakidou, M.

    2015-03-01

    The global atmospheric iron (Fe) cycle is parameterized in the global 3-D chemical transport model TM4-ECPL to simulate the proton- and the organic ligand-promoted mineral Fe dissolution as well as the aqueous-phase photochemical reactions between the oxidative states of Fe(III/II). Primary emissions of total (TFe) and dissolved (DFe) Fe associated with dust and combustion processes are also taken into account. TFe emissions are calculated to amount to ~35 Tg Fe yr-1. The model reasonably simulates the available Fe observations, supporting the reliability of the results of this study. Accounting for proton- and organic ligand-promoted Fe-dissolution in present-day TM4-ECPL simulations, the total Fe-dissolution is calculated to be ~0.163 Tg Fe yr-1 that accounts for up to ~50% of the calculated total DFe emissions. The atmospheric burden of DFe is calculated to be ~0.012 Tg Fe. DFe deposition presents strong spatial and temporal variability with an annual deposition flux ~0.489 Tg Fe yr-1 from which about 25% (~0.124 Tg Fe yr-1) are deposited over the ocean. The impact of air-quality on Fe deposition is studied by performing sensitivity simulations using preindustrial (year 1850), present (year 2008) and future (year 2100) emission scenarios. These simulations indicate that an increase (~2 times) in Fe-dissolution may have occurred in the past 150 years due to increasing anthropogenic emissions and thus atmospheric acidity. On the opposite, a decrease (~2 times) of Fe-dissolution is projected for near future, since atmospheric acidity is expected to be lower than present-day due to air-quality regulations of anthropogenic emissions. The organic ligand contribution to Fe dissolution shows inverse relationship to the atmospheric acidity thus its importance has decreased since the preindustrial period but is projected to increase in the future. The calculated changes also show that the atmospheric DFe supply to High-Nutrient-Low-Chlorophyll oceanic areas (HNLC

  10. Changes in dissolved iron deposition to the oceans driven by human activity: a 3-D global modelling study

    S. Myriokefalitakis

    2015-03-01

    Full Text Available The global atmospheric iron (Fe cycle is parameterized in the global 3-D chemical transport model TM4-ECPL to simulate the proton- and the organic ligand-promoted mineral Fe dissolution as well as the aqueous-phase photochemical reactions between the oxidative states of Fe(III/II. Primary emissions of total (TFe and dissolved (DFe Fe associated with dust and combustion processes are also taken into account. TFe emissions are calculated to amount to ~35 Tg Fe yr−1. The model reasonably simulates the available Fe observations, supporting the reliability of the results of this study. Accounting for proton- and organic ligand-promoted Fe-dissolution in present-day TM4-ECPL simulations, the total Fe-dissolution is calculated to be ~0.163 Tg Fe yr−1 that accounts for up to ~50% of the calculated total DFe emissions. The atmospheric burden of DFe is calculated to be ~0.012 Tg Fe. DFe deposition presents strong spatial and temporal variability with an annual deposition flux ~0.489 Tg Fe yr−1 from which about 25% (~0.124 Tg Fe yr−1 are deposited over the ocean. The impact of air-quality on Fe deposition is studied by performing sensitivity simulations using preindustrial (year 1850, present (year 2008 and future (year 2100 emission scenarios. These simulations indicate that an increase (~2 times in Fe-dissolution may have occurred in the past 150 years due to increasing anthropogenic emissions and thus atmospheric acidity. On the opposite, a decrease (~2 times of Fe-dissolution is projected for near future, since atmospheric acidity is expected to be lower than present-day due to air-quality regulations of anthropogenic emissions. The organic ligand contribution to Fe dissolution shows inverse relationship to the atmospheric acidity thus its importance has decreased since the preindustrial period but is projected to increase in the future. The calculated changes also show that the atmospheric DFe supply to High-Nutrient-Low-Chlorophyll oceanic

  11. a Multiple Data Set Joint Inversion Global 3d P-Velocity Model of the Earth's Crust and Mantle for Improved Seismic Event Location

    Ballard, S.; Begnaud, M. L.; Hipp, J. R.; Chael, E. P.; Encarnacao, A.; Maceira, M.; Yang, X.; Young, C. J.; Phillips, W.

    2013-12-01

    SALSA3D is a global 3D P wave velocity model of the Earth's crust and mantle developed specifically to provide seismic event locations that are more accurate and more precise than are locations from 1D and 2.5D models. In this paper, we present the most recent version of our model, for the first time jointly derived from multiple types of data: body wave travel times, surface wave group velocities, and gravity. The latter two are added to provide information in areas with poor body wave coverage, and are down-weighted in areas where body wave coverage is good. To constrain the inversions, we invoked empirical relations among the density, S velocity, and P velocity. We demonstrate the ability of the new SALSA3D model to reduce mislocations and generate statistically robust uncertainty estimates for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth events. We obtain path-dependent travel time prediction uncertainties for our model by computing the full 3D model covariance matrix of our tomographic system and integrating the model slowness variance and covariance along paths of interest. This approach yields very low travel time prediction uncertainties for well-sampled paths through the Earth and higher uncertainties for paths that are poorly represented in the data set used to develop the model. While the calculation of path-dependent prediction uncertainties with this approach is computationally expensive, uncertainties can be pre-computed for a network of stations and stored in 3D lookup tables that can be quickly and efficiently interrogated using GeoTess software.

  12. Interaction of dipolarization fronts within multiple bursty bulk flows in global MHD simulations of a substorm on 27 February 2009

    Ge, Y. S.; Raeder, J.; Angelopoulos, V.; Gilson, M. L.; Runov, A.

    2011-01-01

    We performed a global MHD simulation of a well-studied substorm on 27 February 2009 (Runov et al., 2009) to understand the generation and large-scale evolution of dipolarization fronts within bursty bulk flows (BBFs). Conjugate, well-positioned Time History of Events and Macroscale Interactions During Substorms (THEMIS) observations from space and ground observatories provide significant constraints to the simulation model. The main substorm onset auroral brightening, at 0749 UT, was in the field of view of Fort Smith (FSMI), just poleward of a preexisting auroral arc. Two minutes later, the space probes recorded a sharp dipolarization front moving sunward, passing by THEMIS and traversing ˜10 RE along the magnetotail. Our global MHD model, OpenGGCM, driven by real-time solar wind/interplanetary magnetic field conditions, is able to reproduce the key features of these signatures. We show that the auroral breakup is caused by the strong flow shear and the flow vortices formed by the BBF flows. Rebound oscillations of the intruding BBF (consistent with recent observations by Panov et al. (2010a)) and filamentation of the front into 1 RE size undulations are superimposed on the flow pattern. Further investigation of the interaction of the BBF and the dipolarization fronts (DFs) reveals that an observed bipolar Bz signature ahead of the DF is due to the interaction between two distinct plasmas emanating from multiple X lines: antisunward-moving flux tubes from a reconnection region at ˜13 RE and sunward-moving dipolarization region within a BBF from a midtail reconnection region at ˜23 RE.

  13. Three-dimensional kinetic-MHD model of the global heliosphere with the heliopause-surface fitting

    Izmodenov, Vladislav

    2015-01-01

    This paper provides a detailed description of the latest version of our model of the solar wind (SW) interaction with the local interstellar medium (LISM). This model has already been applied to the analysis of Lyman-alpha absorption spectra toward nearby stars and for analyses of Solar and Heliospheric Observatory/SWAN data. Katushkina et al. (this issue) used the model results to analyze IBEX-Lo data. At the same time, the details of this model have not yet been published. This is a three-dimensional (3D) kinetic-magnetohydrodynamical (MHD) model that takes into account SW and interstellar plasmas (including $\\alpha$ particles in SW and helium ions in LISM), the solar and interstellar magnetic fields, and the interstellar hydrogen atoms. The latitudinal dependence of SW and the actual flow direction of the interstellar gas with respect to the Sun are also taken into account in the model. It was very essential that our numerical code had been developed in such a way that any numerical diffusion or reconnecti...

  14. MHD simulations on an unstructured mesh

    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

  15. Towards OpenGGCM 4.0. Simulating the global magnetosphere with ring current, radiation belts, plasmasphere, and hall MHD

    Complete text of publication follows. The OpenGGCM is a global numerical model of geospace, covering the Earth's magnetosphere, ionosphere, and thermosphere. Although housed at the University of New Hampshire, it is a community model that is available to any researcher for runs on demand at the Community Coordinated Modeling Center (CCMC, http://ccmc.gsfc.nasa.gov). Like any numerical model it is not perfect and thus subject to continuous development work to add more physics, make it more robust, and to improve performance. In this presentation we will outline a number of ongoing improvements to the model. First, the OpenGGCM is being coupled with the RCM (Rice Convection Model) and CRCM (Comprehensive Ring Current Model) models of the inner magnetosphere. This coupling should improve the realism of the model close to Earth, where particle drift physics is important. The coupling is in both directions, where the RC models receive the ionosphere potential, the magnetic topology, and plasma parameters at the boundaries from the OpenGGCM, and in the opposite direction the OpenGGCM receives pressure and density in the inner magnetosphere, as well as field aligned current and e- precipitation. Second, we couple the NASA/GSFC RBM (Radiation Belt Model) to the OpenGGCM. In this case the coupling is only one way, i.e., the RBM receives particle sources and magnetic topology from the OpenGGCM. Third, we replace the current ionosphere-thermosphere module, CTIM, with the new GT-GIP, which also includes a self-consistent plasmasphere model. Finally, we replace the MHD formalism in the magnetosphere with a Hall-MHD formalism, which is primarily of importance for magnetic reconnection. The new upgraded model will eventually become version 4.0 at the CCMC.

  16. GLOBAL EXISTENCE AND CONVERGENCE RATES OF SMOOTH SOLUTIONS FOR THE 3-D COMPRESSIBLE MAGNETOHYDRODYNAMIC EQUATIONS WITHOUT HEAT CONDUCTIVITY

    Zhensheng GAO; Zhong TAN; Guochun WU

    2014-01-01

    In this paper, we are concerned with the global existence and convergence rates of the smooth solutions for the compressible magnetohydrodynamic equations without heat conductivity, which is a hyperbolic-parabolic system. The global solutions are obtained by combining the local existence and a priori estimates if H3-norm of the initial perturbation around a constant states is small enough and its L1-norm is bounded. A priori decay-in-time estimates on the pressure, velocity and magnetic field are used to get the uniform bound of entropy. Moreover, the optimal convergence rates are also obtained.

  17. Three-Dimensional (3-D) Printing: A Cost-Effective Solution for Improving Global Accessibility to Prostheses.

    Silva, Kyle; Rand, Stephanie; Cancel, David; Chen, Yuxi; Kathirithamby, Rani; Stern, Michelle

    2015-12-01

    The lack of access to prostheses is a global problem, partially caused by the high cost associated with the current manufacturing process. Three-dimensional printing is gaining use in the medical field, and one such area is prosthetics. In addition to using cost-effective materials, this technology allows for rapid prototyping, making it an efficient solution for the development of affordable prostheses. If the rehabilitation medicine community embraces this novel technology, we can help alleviate the global disparity of access to prostheses. PMID:26709247

  18. Optimal estimation of the surface fluxes of methyl chloride using a 3-D global chemical transport model

    X. Xiao

    2009-12-01

    Full Text Available Methyl chloride (CH3Cl is a chlorine-containing trace gas in the atmosphere contributing significantly to stratospheric ozone depletion. Large uncertainties in estimates of its source and sink magnitudes and temporal and spatial variations currently exist. GEIA inventories and other bottom-up emission estimates are used to construct a priori maps of the surface fluxes of CH3Cl. The Model of Atmospheric Transport and Chemistry (MATCH, driven by NCEP interannually varying meteorological data, is then used to simulate CH3Cl mole fractions and quantify the time series of sensitivities of the mole fractions at each measurement site to the surface fluxes of various regional and global sources and sinks. We then implement the Kalman filter (with the unit pulse response method to estimate the surface fluxes on regional/global scales with monthly resolution from January 2000 to December 2004. High frequency observations from the AGAGE, SOGE, NIES, and NOAA/ESRL HATS in situ networks and low frequency observations from the NOAA/ESRL HATS flask network are used to constrain the source and sink magnitudes. The inversion results indicate global total emissions around 4100±470 Gg yr−1 with very large emissions of 2200±390 Gg yr−1 from tropical plants, which turn out to be the largest single source in the CH3Cl budget. Relative to their a priori annual estimates, the inversion increases global annual fungal and tropical emissions, and reduces the global oceanic source. The inversion implies greater seasonal and interannual oscillations of the natural sources and sink of CH3Cl compared to the a priori. The inversion also reflects the strong effects of the 2002/2003 globally widespread heat waves and droughts on global emissions from tropical plants, biomass burning and salt marshes, and on the soil sink.

  19. Optimal estimation of the surface fluxes of methyl chloride using a 3-D global chemical transport model

    X. Xiao

    2010-06-01

    Full Text Available Methyl chloride (CH3Cl is a chlorine-containing trace gas in the atmosphere contributing significantly to stratospheric ozone depletion. Large uncertainties in estimates of its source and sink magnitudes and temporal and spatial variations currently exist. GEIA inventories and other bottom-up emission estimates are used to construct a priori maps of the surface fluxes of CH3Cl. The Model of Atmospheric Transport and Chemistry (MATCH, driven by NCEP interannually varying meteorological data, is then used to simulate CH3Cl mole fractions and quantify the time series of sensitivities of the mole fractions at each measurement site to the surface fluxes of various regional and global sources and sinks. We then implement the Kalman filter (with the unit pulse response method to estimate the surface fluxes on regional/global scales with monthly resolution from January 2000 to December 2004. High frequency observations from the AGAGE, SOGE, NIES, and NOAA/ESRL HATS in situ networks and low frequency observations from the NOAA/ESRL HATS flask network are used to constrain the source and sink magnitudes. The inversion results indicate global total emissions around 4100 ± 470 Gg yr−1 with very large emissions of 2200 ± 390 Gg yr−1 from tropical plants, which turn out to be the largest single source in the CH3Cl budget. Relative to their a priori annual estimates, the inversion increases global annual fungal and tropical emissions, and reduces the global oceanic source. The inversion implies greater seasonal and interannual oscillations of the natural sources and sink of CH3Cl compared to the a priori. The inversion also reflects the strong effects of the 2002/2003 globally widespread heat waves and droughts on global emissions from tropical plants, biomass burning and salt marshes, and on the soil sink.

  20. Towards global scale coastal flood hazard in Delta Cities with 30-meter SRTM and 3D_i

    Winsemius, Hessel; Verhoeven, Govert; Van Leeuwen, Elgard; Van der Klis, Hanneke; Van Wesenbeeck, Bregje; Cumiskey, Lydia; Verlaan, Martin; Muis, Sanne; Ward, Philip; Kwadijk, Jaap

    2015-04-01

    Most attempts to globally simulate inundation at the land-coast interface rely on maximum flood level GIS-based flood spreading models. These are generally not mass conservative, do not account for the genesis of tidal and surges in time, and do not include channel geometry and surface roughness. Furthermore, these methods cannot be used to study the impact of hazard reducing intervention measures that increase roughness at the land-coast interface. These measures include breakwaters and coastal ecosystems, such as mangrove forests and shell fish and coral reefs. Recently, new datasets and models are becoming available that allow us to greatly improve simulation of inundation in global deltas in a rapid and computationally feasible way. In this poster we demonstrate the feasibility of modelling all global deltas with strongly urbanised areas explicitly using these datasets and models. This will allow initiatives such as the 100 resilient cities (Rockefeller foundation) and the 'making cities resilient' campaign (UNISDR) to tackle the issue of coastal flood risk efficiently. We propose to use the following materials: A subgrid enabling 1D-2D model code Outputs from a global tidal and storm surge model Open topographical data We demonstrate the feasibility of this approach by modelling the Mississippi delta with: a) a lidar derived topography dataset (www.gis.ms.gov/); and b) the recently released 30 meter elevation dataset from the Shuttle Radar Topography Mission. We use the new 3Di subgrid code to rapidly schematise the vast delta area with a quadtree mesh. We force the model at the boundaries with water level estimates during the Katrina cyclone. We invite scientists working on global scale inundation modelling to visit our poster in order to discuss possibilities and limitations of the proposed methods related to model codes, data quality and calibration.

  1. First simulation results of Titan's atmosphere dynamics with a global 3-D non-hydrostatic circulation model

    I. V. Mingalev

    2006-09-01

    Full Text Available We present the first results of a 3-D General Circulation Model of Titan's atmosphere which differs from traditional models in that the hydrostatic equation is not used and all three components of the neutral gas velocity are obtained from the numerical solution of the Navier-Stokes equation. The current version of our GCM is, however, a simplified version, as it uses a predescribed temperature field in the model region thereby avoiding the complex simulation of radiative transfer based on the energy equation. We present the first simulation results and compare them to the results of existing GCMs and direct wind observations. The wind speeds obtained from our GCM correspond well with data obtained during the Huygens probe descent through Titan's atmosphere. We interpret the most unexpected feature of these data which consist of the presence of a non-monotonicity of the altitude profile of the zonal wind speed between 60 and 75 km.

  2. Global Well-posedness of the 3D Primitive Equations With Partial Vertical Turbulence Mixing Heat Diffusion

    Cao, Chongsheng

    2010-01-01

    The three--dimensional incompressible viscous Boussinesq equations, under the assumption of hydrostatic balance, govern the large scale dynamics of atmospheric and oceanic motion, and are commonly called the primitive equations. To overcome the turbulence mixing a partial vertical diffusion is usually added to the temperature advection (or density stratification) equation. In this paper we prove the global regularity of strong solutions to this model in a three-dimensional infinite horizontal channel, subject to periodic boundary conditions in the horizontal directions, and with no-penetration and stress-free boundary conditions on the solid, top and bottom, boundaries. Specifically, we show that short time strong solutions to the above problem exist globally in time, and that they depend continuously on the initial data.

  3. The GPlates Portal: Cloud-Based Interactive 3D Visualization of Global Geophysical and Geological Data in a Web Browser.

    Müller, R Dietmar; Qin, Xiaodong; Sandwell, David T; Dutkiewicz, Adriana; Williams, Simon E; Flament, Nicolas; Maus, Stefan; Seton, Maria

    2016-01-01

    The pace of scientific discovery is being transformed by the availability of 'big data' and open access, open source software tools. These innovations open up new avenues for how scientists communicate and share data and ideas with each other and with the general public. Here, we describe our efforts to bring to life our studies of the Earth system, both at present day and through deep geological time. The GPlates Portal (portal.gplates.org) is a gateway to a series of virtual globes based on the Cesium Javascript library. The portal allows fast interactive visualization of global geophysical and geological data sets, draped over digital terrain models. The globes use WebGL for hardware-accelerated graphics and are cross-platform and cross-browser compatible with complete camera control. The globes include a visualization of a high-resolution global digital elevation model and the vertical gradient of the global gravity field, highlighting small-scale seafloor fabric such as abyssal hills, fracture zones and seamounts in unprecedented detail. The portal also features globes portraying seafloor geology and a global data set of marine magnetic anomaly identifications. The portal is specifically designed to visualize models of the Earth through geological time. These space-time globes include tectonic reconstructions of the Earth's gravity and magnetic fields, and several models of long-wavelength surface dynamic topography through time, including the interactive plotting of vertical motion histories at selected locations. The globes put the on-the-fly visualization of massive data sets at the fingertips of end-users to stimulate teaching and learning and novel avenues of inquiry. PMID:26960151

  4. Global left ventricular function in cardiac CT. Evaluation of an automated 3D region-growing segmentation algorithm

    Muehlenbruch, Georg; Das, Marco; Hohl, Christian; Wildberger, Joachim E.; Guenther, Rolf W.; Mahnken, Andreas H. [University Hospital (RWTH) Aachen, Department of Diagnostic Radiology, Aachen (Germany); Rinck, Daniel; Flohr, Thomas G. [Siemens Medical Solutions, Forchheim (Germany); Koos, Ralf; Knackstedt, Christian [University Hospital (RWTH) Aachen, Department of Cardiology, Aachen (Germany)

    2006-05-15

    The purpose was to evaluate a new semi-automated 3D region-growing segmentation algorithm for functional analysis of the left ventricle in multislice CT (MSCT) of the heart. Twenty patients underwent contrast-enhanced MSCT of the heart (collimation 16 x 0.75 mm; 120 kV; 550 mAseff). Multiphase image reconstructions with 1-mm axial slices and 8-mm short-axis slices were performed. Left ventricular volume measurements (end-diastolic volume, end-systolic volume, ejection fraction and stroke volume) from manually drawn endocardial contours in the short axis slices were compared to semi-automated region-growing segmentation of the left ventricle from the 1-mm axial slices. The post-processing-time for both methods was recorded. Applying the new region-growing algorithm in 13/20 patients (65%), proper segmentation of the left ventricle was feasible. In these patients, the signal-to-noise ratio was higher than in the remaining patients (3.2{+-}1.0 vs. 2.6{+-}0.6). Volume measurements of both segmentation algorithms showed an excellent correlation (all P{<=}0.0001); the limits of agreement for the ejection fraction were 2.3{+-}8.3 ml. In the patients with proper segmentation the mean post-processing time using the region-growing algorithm was diminished by 44.2%. On the basis of a good contrast-enhanced data set, a left ventricular volume analysis using the new semi-automated region-growing segmentation algorithm is technically feasible, accurate and more time-effective. (orig.)

  5. A global MHD simulation study of the vortices at the magnetosphere boundary under the southward IMF condition

    Park, K.; Ogino, T.; Lee, D.; Walker, R. J.; Kim, K.

    2013-12-01

    One of the significant problems in magnetospheric physics concerns the nature and properties of the processes which occur at the magnetopause boundary; in particular how energy, momentum, and plasma the magnetosphere receives from the solar wind. Basic processes are magnetic reconnection [Dungey, 1961] and viscouslike interaction, such as Kelvin-Helmholtz instability [Dungey 1955, Miura, 1984] and pressure-pulse driven [Sibeck et al. 1989]. In generally, magnetic reconnection occurs efficiently when the IMF is southward and the rate is largest where the magnetosheath magnetic field is antiparallel to the geomagnetic field. [Sonnerup, 1974; Crooker, 1979; Luhmann et al., 1984; Park et al., 2006, 2009]. The Kelvin-Helmholtz instability is driven by the velocity shear at the boundary, which occur frequently when the IMF is northward. Also variation of the magnetic field and the plasma properties is reported to be quasi-periodic with 2-3min [Otto and Fairfield, 2000] and period of vortex train with 3 to 4 minutes by global MHD simulation [Ogino, 2011]. The pressure-pulse is driven by the solar wind. And the observations of the magnetospheric magnetic field response show quasi-periodic with a period of 8 minutes [Sibeck et al., 1989; Kivelson and Chen, 1995]. There have been few studies of the vortices in the magnetospheric boundary under southward IMF condition. However it is not easy to find the generation mechanism and characteristic for vortices in complicated 3-dimensional space. Thus we have performed global MHD simulation for the steady solar wind and southward IMF conditions. From the simulation results, we find that the vortex occurs at R= 11.7Re (IMF Bz = -2 nT) and R= 10.2Re (IMF Bz = -10 nT) in the dayside magnetopause boundary. Also the vortex rotates counterclockwise in duskside magnetopause (clockwise in dawnside) and propagates tailward. Across the vortex, magnetic field and plasma properties clearly show quasi-periodic fluctuations with a period of 8

  6. Realistic primary and new productions in a 3D global biogeochemical model: biological complexity or physical forcing?

    Popova, E. E.; Coward, A. C.

    2003-04-01

    A nitrogen-based, five compartment biological model has been coupled to a one degree OCCAM (Ocean Circulation and Climate Advanced Modelling Project) model with a KPP ("K profile parameterisation") of the vertical mixing. The biological model state variables are Phytoplankton, Zooplankton, Detritus, Nitrate, and Ammonium. A comparison of the solution with global satellite ocean colour shows that the model is capable of a realistic description of the main seasonal and regional patterns of the surface chlorophyll. Agreement is also good for satellite derived estimates of primary production. In situ data available from local study sites (such as BATS, NABE, India, Papa) are used for the detailed comparison of the model output with the observed ecosystem dynamics in different biological provinces. We discuss performance of the physical and biological model in contrasting areas of the World Ocean. In spite of the biological model being a very simple one, we are able to reproduce the major differences between ecosystem dynamics of these areas. We believe that the success of any global biogeochemical model is dependent first of all on the correct representation of the upper mixed layer (UML) dynamics. Without being able to reproduce contrasting UML regimes in different areas of the World Ocean (such as difference between the North Atlantic and Southern Ocean, or North Atlantic and North Pacific), increased complexity biological models are in danger of producing the right results by the wrong reason.

  7. Normal-mode function representation of global 3-D data sets: open-access software for the atmospheric research community

    Žagar, N.; Kasahara, A.; Terasaki, K.; Tribbia, J.; Tanaka, H.

    2015-04-01

    This article presents new software for the analysis of global dynamical fields in (re)analyses, weather forecasts and climate models. A new diagnostic tool, developed within the MODES project, allows one to diagnose properties of balanced and inertio-gravity (IG) circulations across many scales. In particular, the IG spectrum, which has only recently become observable, can be studied simultaneously in the mass and wind fields while considering the whole model depth in contrast to the majority of studies. The paper includes the theory of normal-mode function (NMF) expansion, technical details of the Fortran 90 code, examples of namelists which control the software execution and outputs of the software application on the ERA Interim reanalysis data set. The applied libraries and default compiler are from the open-source domain. A limited understanding of Fortran suffices for the successful implementation of the software. The presented application of the software to the ERA Interim data set reveals several aspects of the large-scale circulation after it has been partitioned into the linearly balanced and IG components. The global energy distribution is dominated by the balanced energy while the IG modes contribute around 10% of the total wave energy. However, on sub-synoptic scales, IG energy dominates and it is associated with the main features of tropical variability on all scales. The presented energy distribution and features of the zonally averaged and equatorial circulation provide a reference for the validation of climate models.

  8. Normal-mode function representation of global 3-D datasets: an open-access software for atmospheric research community

    Žagar, N.; Kasahara, A.; Terasaki, K.; Tribbia, J.; Tanaka, H.

    2014-12-01

    The paper presents new software for the analysis of global dynamical fields in (re)analyses, weather forecasts and climate models. A new diagnostic tool, developed within the MODES project, allows one to diagnose properties of balanced and inertio-gravity (IG) circulation across many scales. In particular, the IG spectrum, which has only recently become observable, can be studied simultaneously in the mass field and wind field and considering the whole model depth in contrary to majority of studies. The paper presentation includes the theory of normal-mode function expansion, technical details of the Fortran 90 code, examples of namelists which control the software execution and outputs of the software application on the reanalysis dataset ERA Interim. The applied libraries and default compiler are from the open-source domain. A limited understanding of Fortran suffices for the successful implementation of the software. The presented application of the software to the ERA Interim dataset show some features of the large-scale circulation after it has been split into the balanced and IG components. The global energy distribution is dominated by the balanced energy with IG modes making less than 10% of the total wave energy. However, on subsynoptic scales IG energy dominates and it is associated with the main features of tropical variability on all scales. The presented energy distribution and features of the zonally-averaged and equatorial circulation provide a reference for the validation of climate models.

  9. Global effects of transmitted shock wave propagation through the Earth's inner magnetosphere: First results from 3-D hybrid kinetic modeling

    Lipatov, A. S.; Sibeck, D. G.

    2016-09-01

    We use a new hybrid kinetic model to simulate the response of ring current, outer radiation belt, and plasmaspheric particle populations to impulsive interplanetary shocks. Since particle distributions attending the interplanetary shock waves and in the ring current and radiation belts are non-Maxwellian, wave-particle interactions play a crucial role in energy transport within the inner magnetosphere. Finite gyroradius effects become important in mass loading the shock waves with the background plasma in the presence of higher energy ring current and radiation belt ions and electrons. Initial results show that shocks cause strong deformations in the global structure of the ring current, radiation belt, and plasmasphere. The ion velocity distribution functions at the shock front, in the ring current, and in the radiation belt help us determine energy transport through the Earth's inner magnetosphere.

  10. The global distribution of tropospheric NO{sub x} estimated by a 3-D chemical tracer model

    Kraus, A.B.; Rohrer, F.; Ehhalt, D.H. [Forschungszentrum Juelich GmbH (Germany). Inst. fuer Atmosphaerische Chemie

    1997-12-31

    The global distribution of NO{sub x} in the troposphere is calculated using a three-dimensional chemical tracer model with a simplified chemistry scheme for the tracers NO{sub x} {identical_to} NO + NO{sub 2} and HNO{sub 3}. At northern mid- and high latitudes, the calculated tropospheric NO{sub x} content is dominated by the surface source fossil fuel combustion. In the tropical free troposphere lightning discharges provide about 80% of the total NO{sub x} throughout the year. The zonally averaged fractional contribution of aircraft emissions strongly depends on season. The NO mixing ratios determined by the model show good overall agreement with corresponding zonal mean values observed during the STRATOZ III aircraft campaign in June. Over Canada, mixing ratios as high as 0.5-1.0 ppbv NO were measured during TROPOZ II, the origin of which is not yet understood. (author) 8 refs.

  11. THE GLOBAL ARTIFICIAL BOUNDARY CONDITIONS FOR NUMERICAL SIMULATIONS OF THE 3D FLOW AROUND A SUBMERGED BODY

    Hou-de Han; Xin Wen

    2003-01-01

    We consider the numerical approximations of the three-dimensional steady potential flow around a body moving in a liquid of finite constant depth at constant speed and distance below a free surface in a channel. One vertical side is introduced as the upstream artificial boundary and two vertical sides are introduced as the downstream artificial boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the original problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. After solving the variational problem by the finite element method, we obtain the numerical approximation of the original problem. The numerical examples show that the artificial boundary conditions given in this paper are very effective.

  12. Free-boundary ideal MHD stability of W7-X divertor equilibria

    Nührenberg, C.

    2016-07-01

    Plasma configurations describing the stellarator experiment Wendelstein 7-X (W7-X) are computationally established taking into account the geometry of the test-divertor unit and the high-heat-flux divertor which will be installed in the vacuum chamber of the device (Gasparotto et al 2014 Fusion Eng. Des. 89 2121). These plasma equilibria are computationally studied for their global ideal magnetohydrodynamic (MHD) stability properties. Results from the ideal MHD stability code cas3d (Nührenberg 1996 Phys. Plasmas 3 2401), stability limits, spatial structures and growth rates are presented for free-boundary perturbations. The work focusses on the exploration of MHD unstable regions of the W7-X configuration space, thereby providing information for future experiments in W7-X aiming at an assessment of the role of ideal MHD in stellarator confinement.

  13. Magnetopause surface waves triggered by a rotating IMF with the global MHD SWMF/BAT-S-RUS model

    Andriyas, T.; Spencer, E. A.

    2010-12-01

    The solar wind driving of magnetopause surface waves is only partly understood. In particular we do not have a picture of the magnetopause surface wave properties and behavior when a magnetic cloud event, which sometimes involves a rotating IMF, impinges on the magnetosphere. Here we investigate the effect of a twisting or rotational IMF under moderate solar wind velocity (about 500 km/s) upon the magnetosphere with the Global MHD BATS-R-US code. Synthetic solar wind data is constructed to simulate the most important features of a magnetic cloud event, but without including shock features. A sinusoidally varying By component accompanied by a cosinusoidally varying Bz component of the IMF is input into the model with magnitude 10-20 nT. The synthetic data is representative of the magnetic cloud event that occurred on October 3-7 2000. We use the results of the simulation to infer the modes, properties, and particularly the phase speed and wavelength of the surface wave structures.

  14. Proceedings of the workshop on nonlinear MHD and extended MHD

    Nonlinear MHD simulations have proven their value in interpreting experimental results over the years. As magnetic fusion experiments reach higher performance regimes, more sophisticated experimental diagnostics coupled with ever expanding computer capabilities have increased both the need for and the feasibility of nonlinear global simulations using models more realistic than regular ideal and resistive MHD. Such extended-MHD nonlinear simulations have already begun to produce useful results. These studies are expected to lead to ever more comprehensive simulation models in the future and to play a vital role in fully understanding fusion plasmas. Topics include the following: (1) current state of nonlinear MHD and extended-MHD simulations; (2) comparisons to experimental data; (3) discussions between experimentalists and theorists; (4) /equations for extended-MHD models, kinetic-based closures; and (5) paths toward more comprehensive simulation models, etc. Selected papers have been indexed separately for inclusion in the Energy Science and Technology Database

  15. Global sensitivity analysis of a 3D street canyon model—Part II: Application and physical insight using sensitivity analysis

    Benson, James; Ziehn, Tilo; Dixon, Nick S.; Tomlin, Alison S.

    In this work global sensitivity studies using Monte Carlo sampling and high dimensional model representations (HDMR) have been carried out on the k- ɛ closure computational fluid dynamic (CFD) model MISKAM, allowing detailed representation of the effects of changing input parameters on the model outputs. The scenario studied is that of a complex street canyon in the city of York, UK. The sensitivity of the turbulence and mean flow fields to the input parameters is detailed both at specific measurement points and in the associated canyon cross-section to aid comparison with field data. This analysis gives insight into how model parameters can influence the predicted outputs. It also shows the relative strength of each parameter in its influence. Four main input parameters are addressed. Three parameters are surface roughness lengths, determining the flow over a surface, and the fourth is the background wind direction. In order to determine the relative importance of each parameter, sensitivity indices are calculated for the canyon cross-section. The sensitivity of the flow structures in and above the canyon to each parameter is found to be very location dependant. In general, at a particular measurement point, it is the closest wall surface that is most influential on the model output. However, due to the complexity of the flow at different wind angles this is not always the case, for example when a re-circulating canyon flow pattern is present. The background wind direction is shown to be an important parameter as it determines the surface features encountered by the flow. The accuracy with which this is specified when modelling a full-scale situation is therefore an important consideration when considering model uncertainty. Overall, the uncertainty due to roughness lengths is small in comparison to the mean outputs, indicating that the model is well defined even with large ranges of input parameter uncertainty.

  16. CYCLIC THERMAL SIGNATURE IN A GLOBAL MHD SIMULATION OF SOLAR CONVECTION

    Global magnetohydrodynamical simulations of the solar convection zone have recently achieved cyclic large-scale axisymmetric magnetic fields undergoing polarity reversals on a decadal time scale. In this Letter, we show that these simulations also display a thermal convective luminosity that varies in-phase with the magnetic cycle, and trace this modulation to deep-seated magnetically mediated changes in convective flow patterns. Within the context of the ongoing debate on the physical origin of the observed 11 yr variations in total solar irradiance, such a signature supports the thesis according to which all, or part, of the variations on decadal time scales and longer could be attributed to a global modulation of the Sun's internal thermal structure by magnetic activity

  17. A 3D global-to-local deformable mesh model based registration and anatomy-constrained segmentation method for image guided prostate radiotherapy

    Zhou Jinghao; Kim, Sung; Jabbour, Salma; Goyal, Sharad; Haffty, Bruce; Chen, Ting; Levinson, Lydia; Metaxas, Dimitris; Yue, Ning J. [Department of Radiation Oncology, UMDNJ-Robert Wood Johnson Medical School, Cancer Institute of New Jersey, New Brunswick, New Jersey 08903 (United States); Department of Bioinformatics, UMDNJ-Robert Wood Johnson Medical School, Cancer Institute of New Jersey, New Brunswick, New Jersey 08903 (United States); Department of Radiation Oncology, UMDNJ-Robert Wood Johnson Medical School, Cancer Institute of New Jersey, New Brunswick, New Jersey 08903 (United States); Department of Computer Science, Rutgers, State University of New Jersey, Piscataway, New Jersey 08854 (United States); Department of Radiation Oncology, UMDNJ-Robert Wood Johnson Medical School, Cancer Institute of New Jersey, New Brunswick, New Jersey 08903 (United States)

    2010-03-15

    Purpose: In the external beam radiation treatment of prostate cancers, successful implementation of adaptive radiotherapy and conformal radiation dose delivery is highly dependent on precise and expeditious segmentation and registration of the prostate volume between the simulation and the treatment images. The purpose of this study is to develop a novel, fast, and accurate segmentation and registration method to increase the computational efficiency to meet the restricted clinical treatment time requirement in image guided radiotherapy. Methods: The method developed in this study used soft tissues to capture the transformation between the 3D planning CT (pCT) images and 3D cone-beam CT (CBCT) treatment images. The method incorporated a global-to-local deformable mesh model based registration framework as well as an automatic anatomy-constrained robust active shape model (ACRASM) based segmentation algorithm in the 3D CBCT images. The global registration was based on the mutual information method, and the local registration was to minimize the Euclidian distance of the corresponding nodal points from the global transformation of deformable mesh models, which implicitly used the information of the segmented target volume. The method was applied on six data sets of prostate cancer patients. Target volumes delineated by the same radiation oncologist on the pCT and CBCT were chosen as the benchmarks and were compared to the segmented and registered results. The distance-based and the volume-based estimators were used to quantitatively evaluate the results of segmentation and registration. Results: The ACRASM segmentation algorithm was compared to the original active shape model (ASM) algorithm by evaluating the values of the distance-based estimators. With respect to the corresponding benchmarks, the mean distance ranged from -0.85 to 0.84 mm for ACRASM and from -1.44 to 1.17 mm for ASM. The mean absolute distance ranged from 1.77 to 3.07 mm for ACRASM and from 2.45 to

  18. A 3D global-to-local deformable mesh model based registration and anatomy-constrained segmentation method for image guided prostate radiotherapy

    Purpose: In the external beam radiation treatment of prostate cancers, successful implementation of adaptive radiotherapy and conformal radiation dose delivery is highly dependent on precise and expeditious segmentation and registration of the prostate volume between the simulation and the treatment images. The purpose of this study is to develop a novel, fast, and accurate segmentation and registration method to increase the computational efficiency to meet the restricted clinical treatment time requirement in image guided radiotherapy. Methods: The method developed in this study used soft tissues to capture the transformation between the 3D planning CT (pCT) images and 3D cone-beam CT (CBCT) treatment images. The method incorporated a global-to-local deformable mesh model based registration framework as well as an automatic anatomy-constrained robust active shape model (ACRASM) based segmentation algorithm in the 3D CBCT images. The global registration was based on the mutual information method, and the local registration was to minimize the Euclidian distance of the corresponding nodal points from the global transformation of deformable mesh models, which implicitly used the information of the segmented target volume. The method was applied on six data sets of prostate cancer patients. Target volumes delineated by the same radiation oncologist on the pCT and CBCT were chosen as the benchmarks and were compared to the segmented and registered results. The distance-based and the volume-based estimators were used to quantitatively evaluate the results of segmentation and registration. Results: The ACRASM segmentation algorithm was compared to the original active shape model (ASM) algorithm by evaluating the values of the distance-based estimators. With respect to the corresponding benchmarks, the mean distance ranged from -0.85 to 0.84 mm for ACRASM and from -1.44 to 1.17 mm for ASM. The mean absolute distance ranged from 1.77 to 3.07 mm for ACRASM and from 2.45 to

  19. Existence of Global Weak Solutions to a Hybrid Vlasov-MHD Model for Magnetized Plasmas

    Cheng, Bin; Tronci, Cesare

    2016-01-01

    We prove the global-in-time existence of large-data finite-energy weak solutions to an incompressible hybrid Vlasov-magnetohydrodynamic model in three space dimensions. The model couples three essential ingredients of magnetized plasmas: a transport equation for the probability density function, which models energetic rarefied particles of one species; the incompressible Navier--Stokes system for the bulk fluid; and a parabolic evolution equation, involving magnetic diffusivity, for the magnetic field. The physical derivation of our model is given. It is also shown that the weak solution, whose existence is established, has nonincreasing total energy, and that it satisfies a number of physically relevant properties, including conservation of the total momentum, conservation of the total mass, and nonnegativity of the probability density function for the energetic particles. The proof is based on a one-level approximation scheme, which is carefully devised to avoid increase of the total energy for the sequence...

  20. Response of Mercury's Magnetosphere to Solar Wind Forcing: Results of Global MHD Simulations with Coupled Planetary Interior

    Jia, Xianzhe; Slavin, James; Poh, Gangkai; Toth, Gabor; Gombosi, Tamas

    2016-04-01

    As the innermost planet, Mercury arguably undergoes the most direct space weathering interactions due to its weak intrinsic magnetic field and its close proximity to the Sun. It has long been suggested that two processes, i.e., erosion of the dayside magnetosphere due to intense magnetopause reconnection and the shielding effect of the induction currents generated at the conducting core, compete against each other in governing the large-scale structure of Mercury's magnetosphere. An outstanding question concerning Mercury's space weather is which of the two processes is more important. To address this question, we have developed a global MHD model in which Mercury's interior is electromagnetically coupled to the surrounding space environment. As demonstrated in Jia et al. (2015), the new modeling capability allows for self-consistently characterizing the dynamical response of the Mercury system to time-varying external conditions. To assess the relative importance of induction and magnetopause reconnection in controlling the magnetospheric configuration, especially under strong solar driving conditions, we have carried out multiple global simulations that adopt a wide range of solar wind dynamic pressure and IMF conditions. We find that, while the magnetopause standoff distance decreases with increasing solar wind pressure, just as expected, its dependence on the solar wind pressure follows closely a power-law relationship with an index of ~ -1/6, rather than a steeper power-law falling-off expected for the case with only induction present. This result suggests that for the range of solar wind conditions examined, the two competing processes, namely induction and reconnection, appear to play equally important roles in determining the global configuration of Mercury's magnetosphere, consistent with the finding obtained by Slavin et al. (2014) based on MESSENGER observations. We also find that the magnetic perturbations produced by the magnetospheric current systems

  1. 3D video

    Lucas, Laurent; Loscos, Céline

    2013-01-01

    While 3D vision has existed for many years, the use of 3D cameras and video-based modeling by the film industry has induced an explosion of interest for 3D acquisition technology, 3D content and 3D displays. As such, 3D video has become one of the new technology trends of this century.The chapters in this book cover a large spectrum of areas connected to 3D video, which are presented both theoretically and technologically, while taking into account both physiological and perceptual aspects. Stepping away from traditional 3D vision, the authors, all currently involved in these areas, provide th

  2. 3D Animation Essentials

    Beane, Andy

    2012-01-01

    The essential fundamentals of 3D animation for aspiring 3D artists 3D is everywhere--video games, movie and television special effects, mobile devices, etc. Many aspiring artists and animators have grown up with 3D and computers, and naturally gravitate to this field as their area of interest. Bringing a blend of studio and classroom experience to offer you thorough coverage of the 3D animation industry, this must-have book shows you what it takes to create compelling and realistic 3D imagery. Serves as the first step to understanding the language of 3D and computer graphics (CG)Covers 3D anim

  3. Low-cost photodynamic therapy devices for global health settings: Characterization of battery-powered LED performance and smartphone imaging in 3D tumor models

    Hempstead, Joshua; Jones, Dustin P.; Ziouche, Abdelali; Cramer, Gwendolyn M.; Rizvi, Imran; Arnason, Stephen; Hasan, Tayyaba; Celli, Jonathan P.

    2015-05-01

    A lack of access to effective cancer therapeutics in resource-limited settings is implicated in global cancer health disparities between developed and developing countries. Photodynamic therapy (PDT) is a light-based treatment modality that has exhibited safety and efficacy in the clinic using wavelengths and irradiances achievable with light-emitting diodes (LEDs) operated on battery power. Here we assess low-cost enabling technology to extend the clinical benefit of PDT to regions with little or no access to electricity or medical infrastructure. We demonstrate the efficacy of a device based on a 635 nm high-output LED powered by three AA disposable alkaline batteries, to achieve strong cytotoxic response in monolayer and 3D cultures of A431 squamous carcinoma cells following photosensitization by administering aminolevulinic acid (ALA) to induce the accumulation of protoporphyrin IX (PpIX). Here we characterize challenges of battery-operated device performance, including battery drain and voltage stability specifically over relevant PDT dose parameters. Further motivated by the well-established capacity of PDT photosensitizers to serve as tumour-selective fluorescence contrast agents, we demonstrate the capability of a consumer smartphone with low-cost add-ons to measure concentration-dependent PpIX fluorescence. This study lays the groundwork for the on-going development of image-guided ALA-PDT treatment technologies for global health applications.

  4. Nonlinear evolution of MHD instabilities

    A 3-D nonlinear MHD computer code was used to study the time evolution of internal instabilities. Velocity vortex cells are observed to persist into the nonlinear evolution. Pressure and density profiles convect around these cells for a weak localized instability, or convect into the wall for a strong instability. (U.S.)

  5. Global 3-D model of oceanic mercury coupled to carbon biogeochemistry and particle dynamics: application to the transport and fate or riverine mercury

    Zhang, Y.; Jacob, D. J.; Dutkiewicz, S.; Amos, H. M.; Long, M. S.; Sunderland, E. M.

    2014-12-01

    Rivers are estimated to deliver 27 Mmol a-1 of mercury (Hg) to ocean margins, which is comparable to the global atmospheric deposition flux of Hg to the ocean. Previous studies presumed that most of this riverine Hg is sequestered by settling to the coastal regions. However, there has been little investigation of the mechanism and efficiency with which this sequestration takes place, and the implications for riverine influence in different ocean regions. Here we develop a global 3-D chemical transport model for Hg in the ocean (MITgcm-Hg) with ecology (DARWIN model). We track offshore export of the discharged Hg from heterogeneous river systems over different ocean regions, and how it is influenced by the interaction of Hg in a variety of geochemical forms with carbon and suspended particles. We constrain our model assumptions with available offshore observations that bear strong riverine signals. Modeling results suggest that some of the riverine Hg is highly refractory, sorbs strongly to particles and does not follow equilibrium partitioning with the dissolved phase. Simulated global Hg evasion from riverine sources is 50 times larger without this refractory particulate pool, which results in a total evasion flux two times larger than our current best estimate. Based on a typology system of global rivers, we calculate that 10% to 60% of the particulate Hg from different rivers settles in ocean margin sediments because of subgrid sedimentation processes. The remaining 7.5 Mmol a-1 (28% of total river discharge) is available for offshore transport, where it undergoes further sedimentation to the shelf (5.3 Mmol a-1) as well as evasion to the atmosphere (0.44 Mmol a-1). Only 1.7 Mmol a-1 (6.4% of the global riverine Hg) reaches the open ocean, although that fraction varies from 2.6% in East Asia because of the blockage of Korean Peninsula to 25% in east North America facilitated by the Gulf Stream. We find large riverine influences over coastal oceans off East Asia

  6. Using Global Plate Velocity Boundary Conditions for Embedded Regional Geodynamic Models: Application to 3-D Modeling of the Early Rifting of the South Atlantic

    Taramón, Jorge M.; Morgan, Jason P.; Pérez-Gussinyé, Marta

    2016-04-01

    The treatment of far-field boundary conditions is one of the most poorly resolved issues for regional modeling of geodynamic processes. In viscous flow, the choice of far-field boundary conditions often strongly shapes the large-scale structure of a geosimulation. The mantle velocity field along the sidewalls and base of a modeling region is typically much more poorly known than the geometry of past global motions of the surface plates as constrained by global plate motion reconstructions. For regional rifting models it has become routine to apply highly simplified 'plate spreading' or 'uniform rifting' boundary conditions to a 3-D model that limits its ability to simulate the geodynamic evolution of a specific rifted margin. One way researchers are exploring the sensitivity of regional models to uncertain boundary conditions is to use a nested modeling approach in which a global model is used to determine a large-scale flow pattern that is imposed as a constraint along the boundaries of the region to be modeled. Here we explore the utility of a different approach that takes advantage of the ability of finite element models to use unstructured meshes than can embed much higher resolution sub-regions. Here we demonstrate the workflow and code tools that we created to generate this unstructured mesh: solver based on springs, guide-mesh and routines to improve the quality, e.g., closeness to a regular tetrahedron, of the tetrahedral elements of the mesh. Note that the same routines are used to generate a new mesh in the remeshing of a distorted Lagrangian mesh. In our initial project to validate this approach, we create a global spherical shell mesh in which a higher resolution sub-region is created around the nascent South Atlantic Rifting Margin. Global Plate motion BCs and plate boundaries are applied for the time of the onset of rifting, continuing through several 10s of Ma of rifting. Thermal, compositional, and melt-related buoyancy forces are only non

  7. MHD model

    The author's goal is to provide a physical understanding of the ideal MHD model which includes: (1) a basic description of the model, (2) a derivation starting from a more fundamental kinetic model, and (3) a discussion of its range of validity. The ideal MHD model is a single-fluid model that describes the effects of magnetic geometry on the macroscopic equilibrium and stability properties of fusion plasmas. The model is derived in a straight forward manner by forming the mass, momentum, and energy moments of the Boltzmann equation. The moment equations reduce to ideal MHD with the introduction of three critical assumptions: high collisionality, small ion gyro radius, and small resistivity. An analysis of the validity conditions shows that the collision-dominated assumption is never satisfied in plasmas of fusion interest. The remaining two conditions are satisfied by a wide margin. A careful examination of the collision-dominated assumption shows that those particular parts of ideal MHD treated inaccurately (i.e., the parallel momentum and energy equations), play little, if any practical role in MHD equilibrium and stability. These equations primarily describe compression and expansion of a plasma whereas most MHD instabilities involve incompressible motions. The model is incorrect only where it does not matter. This realization leads to the introduction of a modified MHD model known as collisionless MHD which makes predictions nearly identical to collision-dominated assumption. It is thus valid for plasmas of fusion interest. The derivation follows from an analysis of single-particle guiding center motion in a collisionless plasma and the subsequent closure of the system by the heuristic assumption that the motions of interest are incompressible

  8. Reduced Extended MHD

    Morrison, P. J.; Abdelhamid, H. M.; Grasso, D.; Hazeltine, R. D.; Lingam, M.; Tassi, E.

    2015-11-01

    Over the years various reduced fluid models have been obtained for modeling plasmas, with the goal of capturing important physics while maintaining computability. Such models have included the physics contained in various generalizations of Ohm's law, including Hall drift and electron inertia. In a recent publication it was shown that full 3D extended MHD is a Hamiltonian system by finding its noncanonical Poisson bracket. Subsequently, this bracket was shown to be derivable from that for Hall MHD by a series of remarkable transformations, which greatly simplifies the proof of the Jacobi identity and allows one to immediately obtain generalizations of the helicity and cross helicity. In this poster we use this structure to obtain exact reduced fluid models with the effects of full two-fluid theory. Results of numerical computations of collisionless reconnection using an exact reduced 4-field model will be presented and analytical comparisons of mode structure of previous reduced models will be made.

  9. A Global Wave-Driven MHD Solar Model with a Unified Treatment of Open and Closed Magnetic Field Topologies

    Oran, R; Landi, E; Jin, M; Sokolov, I V; Gombosi, T I

    2013-01-01

    We describe, analyze and validate the recently developed Alfv\\'en Wave Solar Model (AWSoM), a 3D global model starting from the top of the chromosphere and extending into interplanetary space (up to 1-2 AU). This model solves the extended two temperature magnetohydrodynamics equations coupled to a wave kinetic equation for low frequency Alfv\\'en waves. In this picture, heating and acceleration of the plasma are due to wave dissipation and wave pressure gradients, respectively. The dissipation process is described by a fully developed turbulent cascade of counter-propagating waves. We adopt a unified approach for calculating the wave dissipation in both open and closed magnetic field lines, allowing for a self-consistent treatment of any magnetic topology. Wave dissipation is the only heating mechanism assumed in the model, and no geometric heating functions are invoked. Electron heat conduction and radiative cooling are also included. We demonstrate that the large-scale, steady-state (in the co-rotating frame...

  10. EUROPEANA AND 3D

    D. Pletinckx

    2012-09-01

    Full Text Available The current 3D hype creates a lot of interest in 3D. People go to 3D movies, but are we ready to use 3D in our homes, in our offices, in our communication? Are we ready to deliver real 3D to a general public and use interactive 3D in a meaningful way to enjoy, learn, communicate? The CARARE project is realising this for the moment in the domain of monuments and archaeology, so that real 3D of archaeological sites and European monuments will be available to the general public by 2012. There are several aspects to this endeavour. First of all is the technical aspect of flawlessly delivering 3D content over all platforms and operating systems, without installing software. We have currently a working solution in PDF, but HTML5 will probably be the future. Secondly, there is still little knowledge on how to create 3D learning objects, 3D tourist information or 3D scholarly communication. We are still in a prototype phase when it comes to integrate 3D objects in physical or virtual museums. Nevertheless, Europeana has a tremendous potential as a multi-facetted virtual museum. Finally, 3D has a large potential to act as a hub of information, linking to related 2D imagery, texts, video, sound. We describe how to create such rich, explorable 3D objects that can be used intuitively by the generic Europeana user and what metadata is needed to support the semantic linking.

  11. Solid works 3D

    This book explains modeling of solid works 3D and application of 3D CAD/CAM. The contents of this book are outline of modeling such as CAD and 2D and 3D, solid works composition, method of sketch, writing measurement fixing, selecting projection, choosing condition of restriction, practice of sketch, making parts, reforming parts, modeling 3D, revising 3D modeling, using pattern function, modeling necessaries, assembling, floor plan, 3D modeling method, practice floor plans for industrial engineer data aided manufacturing, processing of CAD/CAM interface.

  12. Study of Electron Acceleration and Multiple Dipolarization Fronts in 3D kinetic models

    Lapenta, Giovanni; Ashour-Abdalla, Maha; Walker, Raymond; El-Alaoui, Mostafa

    2014-05-01

    The THEMIS mission encountered a depolarization front (DF) during a magnetotail crossing in the interval 035600 - 035900 UT on February 15, 2008 [1]. We present the results of an innovative investigative approach: we combine a global MHD model of the full Earth environment with a local PIC simulation. The global MHD view is provided on the UCLA model applied to the conditions for the interval of interest on Feb 15, 2008. At the specific time of 034800UT, a reconnection site first appear at about x=-15RE, y=4RE. We then use this specific MHD state as the initial setup for a fully kinetic PIC simulation, performed with the iPic3D code [2]. We consider a one way coupling where the MHD state is used as initial state and boundary conditions for the kinetic study [3]. In the present case, the time span of the kinetic simulation is short form the perspective of the global MHD simulation and does not require a full coupling where the MHD then process the information received back from the kinetic run [4]. The fields and particles are advanced self-consistently from the MHD state using a completely kinetic treatment. Many features missed by the MHD model emerge. Most notably a fast reconnection pattern develops and an unsteady reconnection process develops. The typical signatures of fast kinetic reconnection (Hall field) are observed and particle acceleration is obtained self consistently in the fields generated by the PIC simulation. The focus of the presentation will be the mechanisms of unsteady reconnection leading to multiple DFs. We observe intense wave activity propagating off the separatrices. We conduct a spectral analysis to isolate the different wave components in the lower hybrid and whistler regime. The unsteady reconnection and multiple DFs are also analysed in their impact on the energy transfer. We track the conversion of magnetic energy to particle energy and Poynting flux. The processes observed in the simulation are then compared with in situ THEMIS data

  13. Open 3D Projects

    Felician ALECU

    2010-01-01

    Full Text Available Many professionals and 3D artists consider Blender as being the best open source solution for 3D computer graphics. The main features are related to modeling, rendering, shading, imaging, compositing, animation, physics and particles and realtime 3D/game creation.

  14. 3d-3d correspondence revisited

    Chung, Hee-Joong; Dimofte, Tudor; Gukov, Sergei; Sułkowski, Piotr

    2016-04-01

    In fivebrane compactifications on 3-manifolds, we point out the importance of all flat connections in the proper definition of the effective 3d {N}=2 theory. The Lagrangians of some theories with the desired properties can be constructed with the help of homological knot invariants that categorify colored Jones polynomials. Higgsing the full 3d theories constructed this way recovers theories found previously by Dimofte-Gaiotto-Gukov. We also consider the cutting and gluing of 3-manifolds along smooth boundaries and the role played by all flat connections in this operation.

  15. IZDELAVA TISKALNIKA 3D

    Brdnik, Lovro

    2015-01-01

    Diplomsko delo analizira trenutno stanje 3D tiskalnikov na trgu. Prikazan je razvoj in principi delovanja 3D tiskalnikov. Predstavljeni so tipi 3D tiskalnikov, njihove prednosti in slabosti. Podrobneje je predstavljena zgradba in delovanje koračnih motorjev. Opravljene so meritve koračnih motorjev. Opisana je programska oprema za rokovanje s 3D tiskalniki in komponente, ki jih potrebujemo za izdelavo. Diploma se oklepa vprašanja, ali je izdelava 3D tiskalnika bolj ekonomična kot pa naložba v ...

  16. 3D and Education

    Meulien Ohlmann, Odile

    2013-02-01

    Today the industry offers a chain of 3D products. Learning to "read" and to "create in 3D" becomes an issue of education of primary importance. 25 years professional experience in France, the United States and Germany, Odile Meulien set up a personal method of initiation to 3D creation that entails the spatial/temporal experience of the holographic visual. She will present some different tools and techniques used for this learning, their advantages and disadvantages, programs and issues of educational policies, constraints and expectations related to the development of new techniques for 3D imaging. Although the creation of display holograms is very much reduced compared to the creation of the 90ies, the holographic concept is spreading in all scientific, social, and artistic activities of our present time. She will also raise many questions: What means 3D? Is it communication? Is it perception? How the seeing and none seeing is interferes? What else has to be taken in consideration to communicate in 3D? How to handle the non visible relations of moving objects with subjects? Does this transform our model of exchange with others? What kind of interaction this has with our everyday life? Then come more practical questions: How to learn creating 3D visualization, to learn 3D grammar, 3D language, 3D thinking? What for? At what level? In which matter? for whom?

  17. Characterisation of the turbulent electromotive force and its magnetically-mediated quenching in a global EULAG-MHD simulation of solar convection

    Simard, C; Dube, C

    2016-01-01

    We perform a mean-field analysis of the EULAG-MHD millenium simulation of global magnetohydrodynamical convection presented in Passos et al. 2014. The turbulent electromotive force operating in the simulation is assumed to be linearly related to the cyclic axisymmetric mean magnetic field and its first spatial derivatives. At every grid point in the simulation's meridional plane, this assumed relationship involves 27 independent tensorial coefficients. Expanding on Racine et al. 2011, we extract these coefficients from the simulation data through a least-squares minimization procedure based on singular value decomposition. The reconstructed alpha-tensor shows good agreement with that obtained by Racine et al. 2011, who did not include derivatives of the mean-field in their fit, as well as with the alpha-tensor extracted by Augustson et al. 2015 from a distinct ASH MHD simulation. The isotropic part of the turbulent magnetic diffusivity tensor beta is positive definite and reaches values of 5.0x10^7 m2s-1 in t...

  18. 3D Printing our future: Now

    Taylor, Andrew; Unver, Ertu

    2015-01-01

    This 3D Printing our Future:Now talk and visual presentation was given to delegates at the IMI 3D Workshop held at 3M Buckley Innovation Centre on 17th March 2015. The event was hosted by 3Mbuckley Innovation Centre for IMI plc a global engineering company, 3M, and leading 3D additive manufacturing technology providers: EOS, Renishaw and HK 3D printing to disseminate and share their experience on the latest 3D additive design and manufacturing technologies available to the engineering an...

  19. Regularity criterion for the 3D Hall-magneto-hydrodynamics

    Dai, Mimi

    2016-07-01

    This paper studies the regularity problem for the 3D incompressible resistive viscous Hall-magneto-hydrodynamic (Hall-MHD) system. The Kolmogorov 41 phenomenological theory of turbulence [14] predicts that there exists a critical wavenumber above which the high frequency part is dominated by the dissipation term in the fluid equation. Inspired by this idea, we apply an approach of splitting the wavenumber combined with an estimate of the energy flux to obtain a new regularity criterion. The regularity condition presented here is weaker than conditions in the existing criteria (Prodi-Serrin type criteria) for the 3D Hall-MHD system.

  20. Magnetic fields in protoplanetary disks: from MHD simulations to ALMA observations

    Bertrang, Gesa H -M; Wolf, Sebastian

    2016-01-01

    Magnetic fields significantly influence the evolution of protoplanetary disks and the formation of planets, following the predictions of numerous magnetohydrodynamic (MHD) simulations. However, these predictions are yet observationally unconstrained. To validate the predictions on the influence of magnetic fields on protoplanetary disks, we apply 3D radiative transfer simulations of the polarized emission of aligned aspherical dust grains that directly link 3D global non-ideal MHD simulations to ALMA observations. Our simulations show that it is feasible to observe the predicted toroidal large-scale magnetic field structures, not only in the ideal observations but also with high-angular resolution ALMA observations. Our results show further that high angular resolution observations by ALMA are able to identify vortices embedded in outer magnetized disk regions.

  1. 一种全球尺度三维大气数据可视化系统%Visualization System of 3D Global Scale Meteorological Data

    梅鸿辉; 陈海东; 肇昕; 刘昊南; 朱标; 陈为

    2016-01-01

    Meteorological data are multisource, multidimensional, large-scaled and multi-scaled. As a result, it is hard to display a complicated 3d scene in meteorology using traditional visualization methods. This paper presents a visualization system, called AVIS for global scale meteorological data. AVIS employs several standard methods for visualizing density fields, vector fields and tensor fields, as well as non-spatial data. AVIS implements spherical volume rendering and hybrid rendering to support the study of internal structures. AVIS also utilizes a cross-platform and parallel framework to support the visualization both in browsers and on other platforms. The framework benefits from a back-end computing cluster and can accelerate the efficiency of computation and rendering. Use cases verify that the presented system can show meteorological data from many aspects and help users analyze multiple types of data comprehensively.%气象数据具有时变、多源、多维度、数据量大和多尺度等特征,常规的气象数据可视化方法难以满足气象预报与气象研究中的需求.描述了一个新的全球尺度三维大气数据可视化系统——AVIS.系统实现了面向各类密度场、向量场、张量场和非空间数据的可视化方法,设计了适用于大气数据的球面体绘制和混合绘制方法.系统还实现了一套跨平台并行可视化及分析构架,支持在浏览器中或其他显示和交互设备上使用,并利用后端计算集群加速数据的计算与绘制过程.案例分析表明,该系统可以全方位地展现气象数据中的信息,帮助用户综合多种气象信息进行分析.

  2. 3D Equilibrium Reconstructions in DIII-D

    Lao, L. L.; Ferraro, N. W.; Strait, E. J.; Turnbull, A. D.; King, J. D.; Hirshman, H. P.; Lazarus, E. A.; Sontag, A. C.; Hanson, J.; Trevisan, G.

    2013-10-01

    Accurate and efficient 3D equilibrium reconstruction is needed in tokamaks for study of 3D magnetic field effects on experimentally reconstructed equilibrium and for analysis of MHD stability experiments with externally imposed magnetic perturbations. A large number of new magnetic probes have been recently installed in DIII-D to improve 3D equilibrium measurements and to facilitate 3D reconstructions. The V3FIT code has been in use in DIII-D to support 3D reconstruction and the new magnetic diagnostic design. V3FIT is based on the 3D equilibrium code VMEC that assumes nested magnetic surfaces. V3FIT uses a pseudo-Newton least-square algorithm to search for the solution vector. In parallel, the EFIT equilibrium reconstruction code is being extended to allow for 3D effects using a perturbation approach based on an expansion of the MHD equations. EFIT uses the cylindrical coordinate system and can include the magnetic island and stochastic effects. Algorithms are being developed to allow EFIT to reconstruct 3D perturbed equilibria directly making use of plasma response to 3D perturbations from the GATO, MARS-F, or M3D-C1 MHD codes. DIII-D 3D reconstruction examples using EFIT and V3FIT and the new 3D magnetic data will be presented. Work supported in part by US DOE under DE-FC02-04ER54698, DE-FG02-95ER54309 and DE-AC05-06OR23100.

  3. The THEMIS Magnetospheric Breach Discovery and an Anomaly in the Global Distribution of Petroglyphs; MHD Instabilities Recorded by Mankind in Antiquity

    Peratt, Anthony; Yao, W. Fay; Bustamante, P.; Tuki, R.

    2009-05-01

    The recent THEMIS spacecraft discovery of two very large holes in the Earth's magnetosphere helps explain an anomaly in the global distribution of petroglyphs on our planet [1]. Previously, we reported a world wide GPS logging of some 4 million of these objects, each a picture of a filamental MHD instability carved in rock [2, 3]. In all cases, the field-of-view of the petroglyphs was true south with an off-horizon inclination between 21 -- 31 degrees. However, in a complete survey of the braided lava tube caves on Easter Island, petroglyphs were also found in long, true-north shafts, 50 m or more in length. This observation had been noted in natural shafts of similar lengths in the Columbia River Basin. 1. W. Li, To be published in the Journal of Geophysical Research. 2. A. L. Peratt et al, Trans. Plasma Sci. 35. 778. 2007. 3. A. L. Peratt and W. F. Yao, Physica Scripta, T130, 2008.

  4. 3D virtuel udstilling

    Tournay, Bruno; Rüdiger, Bjarne

    2006-01-01

    3d digital model af Arkitektskolens gård med virtuel udstilling af afgangsprojekter fra afgangen sommer 2006. 10 s.......3d digital model af Arkitektskolens gård med virtuel udstilling af afgangsprojekter fra afgangen sommer 2006. 10 s....

  5. Underwater 3D filming

    Roberto Rinaldi

    2014-12-01

    Full Text Available After an experimental phase of many years, 3D filming is now effective and successful. Improvements are still possible, but the film industry achieved memorable success on 3D movie’s box offices due to the overall quality of its products. Special environments such as space (“Gravity” and the underwater realm look perfect to be reproduced in 3D. “Filming in space” was possible in “Gravity” using special effects and computer graphic. The underwater realm is still difficult to be handled. Underwater filming in 3D was not that easy and effective as filming in 2D, since not long ago. After almost 3 years of research, a French, Austrian and Italian team realized a perfect tool to film underwater, in 3D, without any constrains. This allows filmmakers to bring the audience deep inside an environment where they most probably will never have the chance to be.

  6. FISH: A 3D parallel MHD code for astrophysical applications

    Kaeppeli, R; Scheidegger, S; Pen, U -L; Liebendörfer, M

    2009-01-01

    FISH is a fast and simple ideal magneto-hydrodynamics code that scales to ~10 000 processes for a Cartesian computational domain of ~1000^3 cells. The simplicity of FISH has been achieved by the rigorous application of the operator splitting technique, while second order accuracy is maintained by the symmetric ordering of the operators. Between directional sweeps, the three-dimensional data is rotated in memory so that the sweep is always performed in a cache-efficient way along the direction of contiguous memory. Hence, the code only requires a one-dimensional description of the conservation equations to be solved. This approach also enable an elegant novel parallelisation of the code that is based on persistent communications with MPI for cubic domain decomposition on machines with distributed memory. This scheme is then combined with an additional OpenMP parallelisation of different sweeps that can take advantage of clusters of shared memory. We document the detailed implementation of a second order TVD ad...

  7. 3D MHD simulations of subsurface convection in OB stars

    Cantiello, Matteo; Brandenburg, Axel; Del Sordo, Fabio; Käpylä, Petri; Langer, Norbert

    2010-01-01

    During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These p...

  8. Blender 3D cookbook

    Valenza, Enrico

    2015-01-01

    This book is aimed at the professionals that already have good 3D CGI experience with commercial packages and have now decided to try the open source Blender and want to experiment with something more complex than the average tutorials on the web. However, it's also aimed at the intermediate Blender users who simply want to go some steps further.It's taken for granted that you already know how to move inside the Blender interface, that you already have 3D modeling knowledge, and also that of basic 3D modeling and rendering concepts, for example, edge-loops, n-gons, or samples. In any case, it'

  9. Enhanced MHD transport in astrophysical accretion flows: turbulence, winds and jets

    Dobbie, Peter B; Bicknell, Geoffrey V; Salmeron, Raquel

    2009-01-01

    Astrophysical accretion is arguably the most prevalent physical process in the Universe; it occurs during the birth and death of individual stars and plays a pivotal role in the evolution of entire galaxies. Accretion onto a black hole, in particular, is also the most efficient mechanism known in nature, converting up to 40% of accreting rest mass energy into spectacular forms such as high-energy (X-ray and gamma-ray) emission and relativistic jets. Whilst magnetic fields are thought to be ultimately responsible for these phenomena, our understanding of the microphysics of MHD turbulence in accretion flows as well as large-scale MHD outflows remains far from complete. We present a new theoretical model for astrophysical disk accretion which considers enhanced vertical transport of momentum and energy by MHD winds and jets, as well as transport resulting from MHD turbulence. We also describe new global, 3D simulations that we are currently developing to investigate the extent to which non-ideal MHD effects may...

  10. 3D Digital Modelling

    Hundebøl, Jesper

    wave of new building information modelling tools demands further investigation, not least because of industry representatives' somewhat coarse parlance: Now the word is spreading -3D digital modelling is nothing less than a revolution, a shift of paradigm, a new alphabet... Research qeustions. Based...... on empirical probes (interviews, observations, written inscriptions) within the Danish construction industry this paper explores the organizational and managerial dynamics of 3D Digital Modelling. The paper intends to - Illustrate how the network of (non-)human actors engaged in the promotion (and arrest) of 3......D Modelling (in Denmark) stabilizes - Examine how 3D Modelling manifests itself in the early design phases of a construction project with a view to discuss the effects hereof for i.a. the management of the building process. Structure. The paper introduces a few, basic methodological concepts...

  11. 3D Coronal Density Reconstruction and Retrieving the Magnetic Field Structure during Solar Minimum

    Kramar, M; Mikić, Z; Davila, J

    2014-01-01

    Measurement of the coronal magnetic field is a crucial ingredient in understanding the nature of solar coronal phenomena at all scales. We employed STEREO/COR1 data obtained during a deep minimum of solar activity in February 2008 (Carrington rotation CR 2066) to retrieve and analyze the three-dimensional (3D) coronal electron density in the range of heights from 1.5 to 4 Rsun using a tomography method. With this, we qualitatively deduced structures of the coronal magnetic field. The 3D electron density analysis is complemented by the 3D STEREO/EUVI emissivity in the 195 A band obtained by tomography for the same CR. A global 3D MHD model of the solar corona was used to relate the reconstructed 3D density and emissivity to open/closed magnetic field structures. We show that the density maximum locations can serve as an indicator of current sheet position, while the locations of the density gradient maximum can be a reliable indicator of coronal hole boundaries. We find that the magnetic field configuration du...

  12. Professional Papervision3D

    Lively, Michael

    2010-01-01

    Professional Papervision3D describes how Papervision3D works and how real world applications are built, with a clear look at essential topics such as building websites and games, creating virtual tours, and Adobe's Flash 10. Readers learn important techniques through hands-on applications, and build on those skills as the book progresses. The companion website contains all code examples, video step-by-step explanations, and a collada repository.

  13. Some Recent Results on Compressible MHD Equations%关于可压缩MHD方程组的若干研究进展

    张剑文

    2011-01-01

    In this paper, we recall some recent mathematical results on the compressible magnehydrodynamics ( MHD) ,including the global exisntece and uniqeness of classical solution and the vanishing shear viscosity limit to the 1-D MHD equations, the global existence of solution and the incomressible limit to the 3-D MHD equations,and the blow-up criteria of strong solution for the 3-D MHD equations with vacuum states.%主要介绍近年来关于可压缩磁流体力学(MHD)方程组的若干研究进展,主要包括:一维可压缩MHD方程组古典解的存在唯一性和剪切粘性极限,三维可压缩MHD方程组的整体解存在性和不可压极限,以及三维可压缩MHD方程组整体强解的爆破准则.

  14. The M3D-C1 Approach to Simulating 3D 2-fluid Magnetohydrodynamics in Magnetic Fusion Experiments

    A new approach for solving the 3D MHD equations in a strongly magnetized toroidal plasma is presented which uses high-order 2D finite elements with C1 continuity. The vector fields use a physics-based decomposition. An efficient implicit time advance separates the velocity and field advance. ITAPS (SCOREC) adaptivity software and TOPS solvers are used

  15. General Concept of 3D SLAM

    Zhang, Peter; Millos, Evangelous; Gu, Jason

    2009-01-01

    This chapter established an approach to solve the full 3D SLAM problem, applied to an underwater environment. First, a general approach to the 3D SLAM problem was presented, which included the models in 3D case, data association and estimation algorithm. For an underwater mobile robot, a new measurement system was designed for large area's globally-consistent SLAM: buoys for long-range estimation, and camera for short-range estimation and map building. Globally-consistent results could be obt...

  16. 3D Spectroscopic Instrumentation

    Bershady, Matthew A

    2009-01-01

    In this Chapter we review the challenges of, and opportunities for, 3D spectroscopy, and how these have lead to new and different approaches to sampling astronomical information. We describe and categorize existing instruments on 4m and 10m telescopes. Our primary focus is on grating-dispersed spectrographs. We discuss how to optimize dispersive elements, such as VPH gratings, to achieve adequate spectral resolution, high throughput, and efficient data packing to maximize spatial sampling for 3D spectroscopy. We review and compare the various coupling methods that make these spectrographs ``3D,'' including fibers, lenslets, slicers, and filtered multi-slits. We also describe Fabry-Perot and spatial-heterodyne interferometers, pointing out their advantages as field-widened systems relative to conventional, grating-dispersed spectrographs. We explore the parameter space all these instruments sample, highlighting regimes open for exploitation. Present instruments provide a foil for future development. We give an...

  17. 3D Projection Installations

    Halskov, Kim; Johansen, Stine Liv; Bach Mikkelsen, Michelle

    2014-01-01

    Three-dimensional projection installations are particular kinds of augmented spaces in which a digital 3-D model is projected onto a physical three-dimensional object, thereby fusing the digital content and the physical object. Based on interaction design research and media studies, this article...... contributes to the understanding of the distinctive characteristics of such a new medium, and identifies three strategies for designing 3-D projection installations: establishing space; interplay between the digital and the physical; and transformation of materiality. The principal empirical case, From...... Fingerplan to Loop City, is a 3-D projection installation presenting the history and future of city planning for the Copenhagen area in Denmark. The installation was presented as part of the 12th Architecture Biennale in Venice in 2010....

  18. Magnetohydrodynamic (MHD) power generation

    The concept of MHD power generation, principles of operation of the MHD generator, its design, types, MHD generator cycles, technological problems to be overcome, the current state of the art in USA and USSR are described. Progress of India's experimental 5 Mw water-gas fired open cycle MHD power generator project is reported in brief. (M.G.B.)

  19. Herramientas SIG 3D

    Francisco R. Feito Higueruela

    2010-04-01

    Full Text Available Applications of Geographical Information Systems on several Archeology fields have been increasing during the last years. Recent avances in these technologies make possible to work with more realistic 3D models. In this paper we introduce a new paradigm for this system, the GIS Thetrahedron, in which we define the fundamental elements of GIS, in order to provide a better understanding of their capabilities. At the same time the basic 3D characteristics of some comercial and open source software are described, as well as the application to some samples on archeological researchs

  20. Bootstrapping 3D fermions

    Iliesiu, Luca; Kos, Filip; Poland, David; Pufu, Silviu S.; Simmons-Duffin, David; Yacoby, Ran

    2016-03-01

    We study the conformal bootstrap for a 4-point function of fermions in 3D. We first introduce an embedding formalism for 3D spinors and compute the conformal blocks appearing in fermion 4-point functions. Using these results, we find general bounds on the dimensions of operators appearing in the ψ × ψ OPE, and also on the central charge C T . We observe features in our bounds that coincide with scaling dimensions in the GrossNeveu models at large N . We also speculate that other features could coincide with a fermionic CFT containing no relevant scalar operators.

  1. TOWARDS: 3D INTERNET

    Ms. Swapnali R. Ghadge

    2013-01-01

    In today’s ever-shifting media landscape, it can be a complex task to find effective ways to reach your desired audience. As traditional media such as television continue to lose audience share, one venue in particular stands out for its ability to attract highly motivated audiences and for its tremendous growth potential the 3D Internet. The concept of '3D Internet' has recently come into the spotlight in the R&D arena, catching the attention of many people, and leading to a lot o...

  2. 3D Dental Scanner

    Kotek, L.

    2015-01-01

    This paper is about 3D scan of plaster dental casts. The main aim of the work is a hardware and software proposition of 3D scan system for scanning of dental casts. There were used camera, projector and rotate table for this scanning system. Surface triangulation was used, taking benefits of projections of structured light on object, which is being scanned. The rotate table is controlled by PC. The camera, projector and rotate table are synchronized by PC. Controlling of stepper motor is prov...

  3. Interaktiv 3D design

    Villaume, René Domine; Ørstrup, Finn Rude

    2002-01-01

    Projektet undersøger potentialet for interaktiv 3D design via Internettet. Arkitekt Jørn Utzons projekt til Espansiva blev udviklet som et byggesystem med det mål, at kunne skabe mangfoldige planmuligheder og mangfoldige facade- og rumudformninger. Systemets bygningskomponenter er digitaliseret som...... 3D elementer og gjort tilgængelige. Via Internettet er det nu muligt at sammenstille og afprøve en uendelig  række bygningstyper som  systemet blev tænkt og udviklet til....

  4. 3D Harmonic Echocardiography:

    M.M. Voormolen

    2007-01-01

    textabstractThree dimensional (3D) echocardiography has recently developed from an experimental technique in the ’90 towards an imaging modality for the daily clinical practice. This dissertation describes the considerations, implementation, validation and clinical application of a unique

  5. Tangible 3D Modelling

    Hejlesen, Aske K.; Ovesen, Nis

    2012-01-01

    This paper presents an experimental approach to teaching 3D modelling techniques in an Industrial Design programme. The approach includes the use of tangible free form models as tools for improving the overall learning. The paper is based on lecturer and student experiences obtained through...

  6. Shaping 3-D boxes

    Stenholt, Rasmus; Madsen, Claus B.

    2011-01-01

    Enabling users to shape 3-D boxes in immersive virtual environments is a non-trivial problem. In this paper, a new family of techniques for creating rectangular boxes of arbitrary position, orientation, and size is presented and evaluated. These new techniques are based solely on position data...

  7. 3D Equilibrium Effects Due to RMP Application on DIII-D

    Lazerson, S.; Lazarus, E.; Hudson, S.; Pablant, N.; Gates, D.

    2012-06-20

    The mitigation and suppression of edge localized modes (ELMs) through application of resonant magnetic perturbations (RMPs) in Tokamak plasmas is a well documented phenomenon. Vacuum calculations suggest the formation of edge islands and stochastic regions when RMPs are applied to the axisymmetric equilibria. Self-consistent calculations of the plasma equilibrium with the VMEC and SPEC codes have been performed for an up-down symmetric shot in DIII-D. In these codes, a self-consistent calculation of the plasma response due to the RMP coils is calculated. The VMEC code globally enforces the constraints of ideal MHD; consequently, a continuously nested family of flux surfaces is enforced throughout the plasma domain. This approach necessarily precludes the observation of islands or field-line chaos. The SPEC code relaxes the constraints of ideal MHD locally, and allows for islands and field line chaos at or near the rational surfaces. Equilibria with finite pressure gradients are approximated by a set of discrete "ideal-interfaces" at the most irrational flux surfaces and where the strongest pressure gradients are observed. Both the VMEC and SPEC calculations are initialized from EFIT reconstructions of the plasma that are consistent with the experimental pressure and current profiles. A 3D reconstruction using the STELLOPT code, which fits VMEC equilibria to experimental measurements, has also been performed. Comparisons between the equilibria generated by the 3D codes and between STELLOPT and EFIT are presented.

  8. 3D Equilibrium Effects Due to RMP Application on DIII-D

    S. Lazerson, E. Lazarus, S. Hudson, N. Pablant and D. Gates

    2012-06-20

    The mitigation and suppression of edge localized modes (ELMs) through application of resonant magnetic perturbations (RMPs) in Tokamak plasmas is a well documented phenomenon [1]. Vacuum calculations suggest the formation of edge islands and stochastic regions when RMPs are applied to the axisymmetric equilibria. Self-consistent calculations of the plasma equilibrium with the VMEC [2] and SPEC [3] codes have been performed for an up-down symmetric shot (142603) in DIII-D. In these codes, a self-consistent calculation of the plasma response due to the RMP coils is calculated. The VMEC code globally enforces the constraints of ideal MHD; consequently, a continuously nested family of flux surfaces is enforced throughout the plasma domain. This approach necessarily precludes the observation of islands or field-line chaos. The SPEC code relaxes the constraints of ideal MHD locally, and allows for islands and field line chaos at or near the rational surfaces. Equilibria with finite pressure gradients are approximated by a set of discrete "ideal-interfaces" at the most irrational flux surfaces and where the strongest pressure gradients are observed. Both the VMEC and SPEC calculations are initialized from EFIT reconstructions of the plasma that are consistent with the experimental pressure and current profiles. A 3D reconstruction using the STELLOPT code, which fits VMEC equilibria to experimental measurements, has also been performed. Comparisons between the equilibria generated by the 3D codes and between STELLOPT and EFIT are presented.

  9. Modelling of the concentration-time relationship based on global diffusion-charge transfer parameters in a flow-by reactor with a 3D electrode

    Nava, J.L. [Universidad Autonoma Metropolitana-Iztapalapa, Departamento de Quimica, Av. San Rafael Atlixco 186, A.P. 55-534, C.P. 09340, Mexico D.F. (Mexico); Sosa, E. [Instituto Mexicano del Petroleo, Programa de Investigacion en Ingenieria Molecular, Eje Central 152, C.P. 07730, Mexico D.F. (Mexico); Carreno, G. [Universidad de Guanajuato, Facultad de Ingenieria en Geomatica e Hidraulica, Av. Juarez 77, C.P. 36000, Guanajuato, Gto. (Mexico); Ponce-de-Leon, C. [Electrochemical Engineering Group, School of Engineering Sciences, University of Southampton, Highfield, Southampton SO17 1BJ (United Kingdom)]. E-mail: capla@soton.ac.uk; Oropeza, M.T. [Centro de Graduados e Investigacion del Instituto Tecnologico de Tijuana, Blvd. Industrial, s/n, C.P. 22500, Tijuana B.C. (Mexico)

    2006-05-25

    A concentration versus time relationship model based on the isothermal diffusion-charge transfer mechanism was developed for a flow-by reactor with a three-dimensional (3D) reticulated vitreous carbon (RVC) electrode. The relationship was based on the effectiveness factor ({eta}) which lead to the simulation of the concentration decay at different electrode polarisation conditions, i.e. -0.1, -0.3 and -0.59 V versus SCE; the charge transfer process was used for the former and mix and a mass transport control was used for the latter. Charge transfer and mass transport parameters were estimated from experimental data using Electrochemical Impedance Spectroscopy (EIS) and Linear Voltammetry (LV) techniques, respectively.

  10. Simulations of drift waves in 3D magnetic configurations

    Drift waves are commonly held responsible for anomalous transport in tokamak configurations and in particular for the anomalously high heat loss. The next generation of stellarators on the other hand are hoped to be characterized by a much smaller neo-classical transport and by particle confinement close to that of tokamaks. There is nevertheless a strong interest in the stellarator community to study the properties of drift waves in 3D magnetic configurations. To serve this interest we have developed the first global gyrokinetic code, EUTERPE, aimed at the investigation of linear drift wave stability in general toroidal geometry. The physical model assumes electrostatic waves and adiabatic electrons. EUTERPE is a particle-in-cell (PIC) code in which the gyrokinetic Poisson equation is discretized with the finite element method defined in the PEST -1 system of magnetic coordinates. The magnetic geometry is provided by the magnetohydrodynamic (MHD) equilibrium code VMEC. The complete 3D model has been successfully validated in toroidal axisymmetric and straight helical geometries and has permitted the first simulation of unstable global ITG driven modes in non-axisymmetric toroidal configurations. As a first application, two configurations have been studied, the Quasi-Axially symmetric Stellarator with three fields periods (QAS3) currently one system under consideration at the Princeton Plasma Physics Laboratory and the Helically Symmetric experiment (HSX) which has recently started operation at the University of Wisconsin. QAS3 is characterized by a tokamak-Iike field in the outer part of the torus. In this structure the drift waves are mainly affected by the magnetic shear and barely by the shape of the plasma. Also, the results are very close to those obtained for a tokamak. On the other hand, results for the HSX configuration, which is characterized by a dominant helical magnetic field, show a clear 3D effect, namely a strong toroidal variation of the drift wave