Science.gov

Sample records for 3d particle simulations

  1. 3-D Particle Simulation of Current Sheet Instabilities

    NASA Astrophysics Data System (ADS)

    Wang, Zhenyu; Lin, Yu; Wang, Xueyi; Tummel, Kurt; Chen, Liu

    2015-11-01

    The electrostatic (ES) and electromagnetic (EM) instabilities of a Harris current sheet are investigated using a 3-D linearized (δf) gyrokinetic (GK) electron and fully kinetic (FK) ion (GeFi) particle simulation code. The equilibrium magnetic field consists of an asymptotic anti-parallel Bx 0 and a guide field BG. The ES simulations show the excitation of lower-hybrid drift instability (LHDI) at the current sheet edge. The growth rate of the 3-D LHDI is scanned through the (kx ,ky) space. The most unstable modes are found to be at k∥ = 0 for smaller ky. As ky increases, the growth rate shows two peaks at k∥ ≠ 0 , consistent with analytical GK theory. The eigenmode structure and growth rate of LHDI obtained from the GeFi simulation agree well with those obtained from the FK PIC simulation. Decreasing BG, the asymptotic βe 0, or background density can destabilize the LHDI. In the EM simulation, tearing mode instability is dominant in the cases with ky kx , there exist two unstable modes: a kink-like (LHDI) mode at the current sheet edge and a sausage-like mode at the sheet center. The results are compared with the GK eigenmode theory and the FK simulation.

  2. Development of a 3D particle treecode for plasma simulations

    NASA Astrophysics Data System (ADS)

    Ong, Benjamin; Christlieb, Andrew; Krasny, Robert

    2008-11-01

    In this work we present a fully 3-D Boundary Integral Treecode (BIT). We apply the method to several classic problems such as sheath formation and 3D simulations of a Penning trap. In addition, we investigate the ability of the solver to naturally capture Coloumb scattering. A key point in the investigation is to understand the effect of different types of regularizations, and how to appropriately incorporate the regularization in the BIT framework. This work builds on substantial efforts in 1- and 2-D. [1] R. Krasny and K. Lindsay, A particle method and adaptive treecode for vortex sheet motion in 3-D flow, JCP, Vol. 172, No. 2, 879-907 [2] K. Matyash, R. Schneider, R. Sydora, and F. Taccogna, Application of a Grid-Free Kinetic Model to the Collisionless Sheath, Contrib. Plasma Phys, Vol. 48, No. 1-3, 116-120 (2008) [3] K. Cartwright and A. Christlieb, Boundary Integral Corrected Particle in Cell, SIAM Journal on Sci. Comput., submitted [4] A. Christlieb, R. Krasny, B. Ong and J. Qiu, A Step Towards Addressing Temporal Multi-scale Problems in Plasma Physics, in prep.

  3. 3-D Spreadsheet Simulation of a Modern Particle Detector

    ERIC Educational Resources Information Center

    Scott, Alan J.

    2004-01-01

    A spreadsheet simulation of a modern particle detector has been developed and can be readily used as an instructional tool in the physics classroom. The spreadsheet creates a three-dimensional model that can be rotated and helical trajectories can be highlighted. An associated student worksheet is also presented.

  4. Particle Acceleration in the Low Corona Over Broad Longitudes: Coupling MHD and 3D Particle Simulations

    NASA Astrophysics Data System (ADS)

    Gorby, M.; Schwadron, N.; Torok, T.; Downs, C.; Lionello, R.; Linker, J.; Titov, V. S.; Mikic, Z.; Riley, P.; Desai, M. I.; Dayeh, M. A.

    2014-12-01

    Recent work on the coupling between the Energetic Particle Radiation Environment Module (EPREM, a 3D energetic particle model) and Magnetohydrodynamics Around a Sphere (MAS, an MHD code developed at Predictive Science, Inc.) has demonstrated the efficacy of compression regions around fast coronal mass ejections (CMEs) for particle acceleration low in the corona (˜ 3 - 6 solar radii). These couplings show rapid particle acceleration over a broad longitudinal extent (˜ 80 degrees) resulting from the pile-up of magnetic flux in the compression regions and their subsequent expansion. The challenge for forming large SEP events in such compression-acceleration scenarios is to have enhanced scattering within the acceleration region while also allowing for efficient escape of accelerated particles downstream (away from the Sun) from the compression region. We present here the most recent simulation results including energetic particle and CME plasma profiles, the subsequent flux and dosages at 1AU, and an analysis of the compressional regions as efficient accelerators.

  5. Effect of Model Scale and Particle Size Distribution on PFC3D Simulation Results

    NASA Astrophysics Data System (ADS)

    Ding, Xiaobin; Zhang, Lianyang; Zhu, Hehua; Zhang, Qi

    2014-11-01

    This paper investigates the effect of model scale and particle size distribution on the simulated macroscopic mechanical properties, unconfined compressive strength (UCS), Young's modulus and Poisson's ratio, using the three-dimensional particle flow code (PFC3D). Four different maximum to minimum particle size ( d max/ d min) ratios, all having a continuous uniform size distribution, were considered and seven model (specimen) diameter to median particle size ratios ( L/ d) were studied for each d max/ d min ratio. The results indicate that the coefficients of variation (COVs) of the simulated macroscopic mechanical properties using PFC3D decrease significantly as L/ d increases. The results also indicate that the simulated mechanical properties using PFC3D show much lower COVs than those in PFC2D at all model scales. The average simulated UCS and Young's modulus using the default PFC3D procedure keep increasing with larger L/ d, although the rate of increase decreases with larger L/ d. This is mainly caused by the decrease of model porosity with larger L/ d associated with the default PFC3D method and the better balanced contact force chains at larger L/ d. After the effect of model porosity is eliminated, the results on the net model scale effect indicate that the average simulated UCS still increases with larger L/ d but the rate is much smaller, the average simulated Young's modulus decreases with larger L/ d instead, and the average simulated Poisson's ratio versus L/ d relationship remains about the same. Particle size distribution also affects the simulated macroscopic mechanical properties, larger d max/ d min leading to greater average simulated UCS and Young's modulus and smaller average simulated Poisson's ratio, and the changing rates become smaller at larger d max/ d min. This study shows that it is important to properly consider the effect of model scale and particle size distribution in PFC3D simulations.

  6. M3D-K simulations of sawteeth and energetic particle transport in tokamak plasmas

    NASA Astrophysics Data System (ADS)

    Shen, Wei; Fu, G. Y.; Sheng, Zheng-Mao; Breslau, J. A.; Wang, Feng

    2014-09-01

    Nonlinear simulations of sawteeth and related energetic particle transport are carried out using the kinetic/magnetohydrodynamic (MHD) hybrid code M3D-K. MHD simulations show repeated sawtooth cycles for a model tokamak equilibrium. Furthermore, test particle simulations are carried out to study the energetic particle transport due to a sawtooth crash. The results show that energetic particles are redistributed radially in the plasma core, depending on pitch angle and energy. For trapped particles, the redistribution occurs for particle energy below a critical value in agreement with existing theories. For co-passing particles, the redistribution is strong with little dependence on particle energy. In contrast, the redistribution level of counter-passing particles decreases with increasing particle energy.

  7. M3D-K Simulations of Sawteeth and Energetic Particle Transport in Tokamak Plasmas

    NASA Astrophysics Data System (ADS)

    Shen, Wei; Fu, Guoyong; Sheng, Zhengmao; Breslau, Joshua; Wang, Feng

    2013-10-01

    Nonlinear simulations of Sawteeth and energetic particle transport are carried out using the kinetic/MHD hybrid code M3D-K. MHD simulations show repeated sawtooth cycles due to a resistive (1,1) internal kink mode for a model tokamak equilibrium. Furthermore, test particle simulations are carried out to study the energetic particle transport due to a sawtooth crash. The results show that energetic particles are redistributed radially in plasma core depending on pitch angle and energy. For trapped particles, the redistribution occurs for particle energy below a critical value in agreement with previous theory. For co-passing particles, the redistribution is strong with little dependence on particle energy. In contrast, the redistribution level of counter-passing particles decreases as particle energy becomes large.

  8. M3D-K simulations of sawteeth and energetic particle transport in tokamak plasmas

    SciTech Connect

    Shen, Wei; Sheng, Zheng-Mao; Fu, G. Y.; Breslau, J. A.; Wang, Feng

    2014-09-15

    Nonlinear simulations of sawteeth and related energetic particle transport are carried out using the kinetic/magnetohydrodynamic (MHD) hybrid code M3D-K. MHD simulations show repeated sawtooth cycles for a model tokamak equilibrium. Furthermore, test particle simulations are carried out to study the energetic particle transport due to a sawtooth crash. The results show that energetic particles are redistributed radially in the plasma core, depending on pitch angle and energy. For trapped particles, the redistribution occurs for particle energy below a critical value in agreement with existing theories. For co-passing particles, the redistribution is strong with little dependence on particle energy. In contrast, the redistribution level of counter-passing particles decreases with increasing particle energy.

  9. Simulation of surface tension in 2D and 3D with smoothed particle hydrodynamics method

    NASA Astrophysics Data System (ADS)

    Zhang, Mingyu

    2010-09-01

    The methods for simulating surface tension with smoothed particle hydrodynamics (SPH) method in two dimensions and three dimensions are developed. In 2D surface tension model, the SPH particle on the boundary in 2D is detected dynamically according to the algorithm developed by Dilts [G.A. Dilts, Moving least-squares particle hydrodynamics II: conservation and boundaries, International Journal for Numerical Methods in Engineering 48 (2000) 1503-1524]. The boundary curve in 2D is reconstructed locally with Lagrangian interpolation polynomial. In 3D surface tension model, the SPH particle on the boundary in 3D is detected dynamically according to the algorithm developed by Haque and Dilts [A. Haque, G.A. Dilts, Three-dimensional boundary detection for particle methods, Journal of Computational Physics 226 (2007) 1710-1730]. The boundary surface in 3D is reconstructed locally with moving least squares (MLS) method. By transforming the coordinate system, it is guaranteed that the interface function is one-valued in the local coordinate system. The normal vector and curvature of the boundary surface are calculated according to the reconstructed boundary surface and then surface tension force can be calculated. Surface tension force acts only on the boundary particle. Density correction is applied to the boundary particle in order to remove the boundary inconsistency. The surface tension models in 2D and 3D have been applied to benchmark tests for surface tension. The ability of the current method applying to the simulation of surface tension in 2D and 3D is proved.

  10. Test-particle Orbit Simulations in Fields from a Realistic 3D MHD Simulation

    NASA Astrophysics Data System (ADS)

    Decker, R. B.; Opher, M.; Hill, M. E.

    2007-05-01

    Models designed to explore the global structure of the heliosphere have become increasing sophisticated. Incentives to increase and to further explore the predictive capabilities of such models include the entry of the Voyager spacecraft into the foreshock region of the termination shock (TS), Voyager 1 in mid-2002 and Voyager 2 in late 2004, and the crossing of the TS and passage into the heliosheath (HSH) of Voyager 1 in 2004 day 351. Using the electric and magnetic fields generated by a MHD model of a 3D, asymmetric heliosphere [Opher et al., Ap. J. L., 640, 2006], we have developed full-particle and adiabatic-orbit codes to simulate the motion of test particles in the solar wind, TS, and HSH environments. The full-particle orbits are necessary to investigate energetic ion (e.g., anomalous and galactic cosmic ray) motion at the TS and within the heliospheric current sheet that is included in the MHD model. Adiabatic orbits are used to study particle motion in the much larger volume of the HSH where the non-homogeneous model fields produce complex guiding center motions, including mirroring in local field compressions. We will present results from these orbit computations, which are intended to provide an initial, albeit simplified, look at the propagation of high-energy charged particles, in the scatter-free limit, in the best model of the TS/HSH field configurations currently available. We will also display drift paths of high-energy ions in the HSH fields using the guiding center drift equations that are applicable in the limit of diffusive propagation.

  11. GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method

    SciTech Connect

    Gong Chunye; Liu Jie; Chi Lihua; Huang Haowei; Fang Jingyue; Gong Zhenghu

    2011-07-01

    Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates (S{sub n}) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.

  12. GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method

    NASA Astrophysics Data System (ADS)

    Gong, Chunye; Liu, Jie; Chi, Lihua; Huang, Haowei; Fang, Jingyue; Gong, Zhenghu

    2011-07-01

    Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates ( Sn) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.

  13. 3-D Particl-in-Cell Simulations of Transport Driven Currents

    NASA Astrophysics Data System (ADS)

    Tsung, F. S.; Dawson, J. M.

    1997-11-01

    In the advanced tokamak regime, transport phenomena can account for a signficant fraction of the toroidal current, possibly over that driven directly by the ohmic heating electric fields. Although bootstrap theory accounts for contributions of the collisional modification of banana orbits on the toroidal currents, the corresponding transport theory does not accurately predict the transport of particles and heat in present-day tokamak experiments. Furthermore, in our previous simulations in 21/2-D, currents were spontaneously generated in both the cylindrical and the toroidal geometries, contrary to neoclassical predictions. In these calculations, it was believed that the driving mechanism is the preferential loss of particles whose initial velocity is opposite to that of the plasma current. Because the preferential loss mechanism assumes the conservation of toroidal angular momentum, we have extended these simulations to three dimensions to study the effects of toroidal assymetries. A parallel, 3-D electromagnetic PIC code running on the IBM SP, with a localized field-solver has been developed to investigate the effects of perturbations parallel to the field lines, and direct comparisons has been made between the 21/2-D and 3-D simulations, and we have found good agreements between the 2 1/2-D calculations and the 3-D results. We will present these results at the meeting.

  14. A 3d particle simulation code for heavy ion fusion accelerator studies

    SciTech Connect

    Friedman, A.; Bangerter, R.O.; Callahan, D.A.; Grote, D.P.; Langdon, A.B. ); Haber, I. )

    1990-06-08

    We describe WARP, a new particle-in-cell code being developed and optimized for ion beam studies in true geometry. We seek to model transport around bends, axial compression with strong focusing, multiple beamlet interaction, and other inherently 3d processes that affect emittance growth. Constraints imposed by memory and running time are severe. Thus, we employ only two 3d field arrays ({rho} and {phi}), and difference {phi} directly on each particle to get E, rather than interpolating E from three meshes; use of a single 3d array is feasible. A new method for PIC simulation of bent beams follows the beam particles in a family of rotated laboratory frames, thus straightening'' the bends. We are also incorporating an envelope calculation, an (r, z) model, and 1d (axial) model within WARP. The BASIS development and run-time system is used, providing a powerful interactive environment in which the user has access to all variables in the code database. 10 refs., 3 figs.

  15. Earth's Magnetosphere 3D Simulation by Coupling Particle-In-Cell and Magnetohydrodynamics Models: Parametric Study

    NASA Astrophysics Data System (ADS)

    Baraka, S. M.; Ben-Jaffel, L. B.

    2014-12-01

    We use particle-in-cell PIC 3D Electromagnetic, relativistic global code to address large-scale problems in magnetosphere electrodynamics. Terrestrial bow shock is simulated as an example. 3D Magnetohydrodynamics model ,MHD GUMICS in CCMC project, have been used in parallel with PIC under same scaled Solar wind (SW) and IMF conditions. We report new results from the coupling between the two models. Further investigations are required for confirmations of these results. In both codes the Earth's bow shock position is found at ~14.8 RE along the Sun-Earth line, and ~29 RE on the dusk side which is consistent with past in situ observation. Both simulations reproduce the theoretical jump conditions at the shock. However, PIC code density and temperature distributions are inflated and slightly shifted sunward when compared to MHD results. Reflected ions upstream of the bow shock may cause this sunward shift for density and temperature. Distribution of reflected ions and electrons are shown in the foreshock region, within the transition of the shock and in the downstream. The current version of PIC code can be run under modest computing facilities and resources. Additionally, existing MHD simulations should be useful to calibrate scaled properties of plasma resulting from PIC simulations for comparison with observations. Similarities and drawbacks of the results obtained by the two models are listed. The ultimate goal of using these different models in a complimentary manner rather than competitive is to better understand the macrostructure of the magnetosphere

  16. 3D Printing Meets Computational Astrophysics: Deciphering the Structure of Eta Carinae’s Colliding Winds Using 3D Prints of Smoothed Particle Hydrodynamics Simulations

    NASA Astrophysics Data System (ADS)

    Madura, Thomas; Gull, Theodore R.; Clementel, Nicola; Paardekooper, Jan-Pieter; Kruip, Chael; Corcoran, Michael F.; Hamaguchi, Kenji; Teodoro, Mairan

    2015-01-01

    We present the first 3D prints of output from a supercomputer simulation of a complex astrophysical system, the colliding stellar winds in the massive (>120 MSun), highly eccentric (e ~ 0.9) binary Eta Carinae. Using a consumer-grade 3D printer (Makerbot Replicator 2X), we successfully printed 3D smoothed particle hydrodynamics simulations of Eta Carinae's inner (r ~110 AU) wind-wind collision interface at multiple orbital phases. These 3D prints reveal important, previously unknown 'finger-like' structures at orbital phases shortly after periastron (φ ~1.045) that protrude radially outward from the spiral wind-wind collision region. We speculate that these fingers are related to instabilities (e.g. Rayleigh-Taylor) that arise at the interface between the radiatively-cooled layer of dense post-shock primary-star wind and the hot, adiabatic post-shock companion-star wind. The success of our work and easy identification of previously unknown physical features highlight the important role 3D printing can play in the visualization and understanding of complex 3D time-dependent numerical simulations of astrophysical phenomena.

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

    NASA Astrophysics Data System (ADS)

    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.

  18. 3D Particle-In-Cell (PIC) simulations of plasma sheath formation above lunar craters

    NASA Astrophysics Data System (ADS)

    Likhanskii, A.; Poppe, A. R.; Piquette, M.; Amyx, K.; Messmer, P.; Horanyi, M.

    2010-12-01

    Comprehensive investigation of plasma sheath formation and consequent dust levitation on lunar surface is important for interpretation of results of future lunar missions (such as LADEE and ARTEMIS). Until recently, most of such studies were done in experimental laboratories at reduced scales. Due to the complexity and nonlinearity of the problem, only simplified theories, describing this effect, were developed. However, recent progress in high-performance kinetic plasma simulations allowed tackling the problem of plasma sheath formation numerically. In this poster we will present the simulation results of plasma sheath formation above the lunar craters in presence of solar wind and photoelectron emission. These results were obtained using 3D Particle-In-Cell (PIC) code VORPAL. In the simulations we considered plasma sheath formation for normal, 45 and 90 degree incidence solar wind. Sample distribution of electric field in plasma sheath is shown in Figure 1. In the second part of the poster, we will present results of simulations of the LASP (Laboratory for Atmospheric and Space Physics at University of Colorado) experiments on study of plasma sheath formation above hemispherical isolated dimple. Figure 1. Electric field distribution in the plasma sheath above the lunar crater

  19. 3D and r,z particle simulations of heavy ion fusion beams

    NASA Astrophysics Data System (ADS)

    Friedman, A.; Grote, D. P.; Callahan, D. A.; Langdon, A. B.; Haber, I.

    1992-08-01

    The space-charge-dominated beams in a heavy ion beam driven inertial fusion (HIF) accelerator must be focused onto small (few mm) spots at the fusion target, and so preservation of a small emittance is crucial. The nonlinear beam self-fields can lead to emittance growth; thus, a self-consistent field description is necessary. We have developed a multi-dimensional time-dependent discrete particle simulation code, WARP, and are using it to study the behavior of HIF beams. The code's 3d package combines features of an accelerator code and a particle-in-cell (PIC) plasma simulation. Novel techniques allow it to follow beams through many accelerator elements over long distances and around bends. We have used the code to understand the emittance growth observed in the MBE4 experiment at Lawrence Berkeley Laboratory (LBL) under conditions of aggressive drift-compression. We are currently applying it to LBL's planned ILSE experiments, and (most recently) to an ESQ injector option being evaluated for ILSE. The code's r, z package is being used to study the axial confinement afforded by the shaped ends of the accelerating pulses, and to study longitudinal instability induced by induction module impedance.

  20. 3D Surgical Simulation

    PubMed Central

    Cevidanes, Lucia; Tucker, Scott; Styner, Martin; Kim, Hyungmin; Chapuis, Jonas; Reyes, Mauricio; Proffit, William; Turvey, Timothy; Jaskolka, Michael

    2009-01-01

    This paper discusses the development of methods for computer-aided jaw surgery. Computer-aided jaw surgery allows us to incorporate the high level of precision necessary for transferring virtual plans into the operating room. We also present a complete computer-aided surgery (CAS) system developed in close collaboration with surgeons. Surgery planning and simulation include construction of 3D surface models from Cone-beam CT (CBCT), dynamic cephalometry, semi-automatic mirroring, interactive cutting of bone and bony segment repositioning. A virtual setup can be used to manufacture positioning splints for intra-operative guidance. The system provides further intra-operative assistance with the help of a computer display showing jaw positions and 3D positioning guides updated in real-time during the surgical procedure. The CAS system aids in dealing with complex cases with benefits for the patient, with surgical practice, and for orthodontic finishing. Advanced software tools for diagnosis and treatment planning allow preparation of detailed operative plans, osteotomy repositioning, bone reconstructions, surgical resident training and assessing the difficulties of the surgical procedures prior to the surgery. CAS has the potential to make the elaboration of the surgical plan a more flexible process, increase the level of detail and accuracy of the plan, yield higher operative precision and control, and enhance documentation of cases. Supported by NIDCR DE017727, and DE018962 PMID:20816308

  1. 3D Unsteady Multiphase Simulation of Uranium Tetrafluoride Particle Fluorination in Fluidized Bed Pilot

    NASA Astrophysics Data System (ADS)

    Konan, N. A.; Neau, H.; Simonin, O.; Dupoizat, M.; Le Geaziou, T.

    This paper investigates Eulerian simulation approach of uranium hexafluoride production in fluidized bed pilot. Mass transfer is modeled by using the shrinking particle model. The model successfully predicts expected amount of uranium hexafluoride. As heat transfers with wall are neglected, temperature increases within the reactor and there is no gradient in wall-normal direction of reactor. External diffusion model of reactive gas around particle is developed. Effect of fluorine diffusion within nitrogen is found to be negligible under the simulation conditions. Moreover, inter-particle radiative heat transfer between particles of dense phase in the bed is investigated in the frame of Rosseland approximation.

  2. 3D hybrid simulations with gyrokinetic particle ions and fluid electrons

    SciTech Connect

    Belova, E.V.; Park, W.; Fu, G.Y.; Strauss, H.R.; Sugiyama, L.E.

    1998-12-31

    The previous hybrid MHD/particle model (MH3D-K code) represented energetic ions as gyrokinetic (or drift-kinetic) particles coupled to MHD equations using the pressure or current coupling scheme. A small energetic to bulk ion density ratio was assumed, n{sub h}/n{sub b} {much_lt} 1, allowing the neglect of the energetic ion perpendicular inertia in the momentum equation and the use of MHD Ohm`s law E = {minus}v{sub b} {times} B. A generalization of this model in which all ions are treated as gyrokinetic/drift-kinetic particles and fluid description is used for the electron dynamics is considered in this paper.

  3. Simulation of the 3D viscoelastic free surface flow by a parallel corrected particle scheme

    NASA Astrophysics Data System (ADS)

    Jin-Lian, Ren; Tao, Jiang

    2016-02-01

    In this work, the behavior of the three-dimensional (3D) jet coiling based on the viscoelastic Oldroyd-B model is investigated by a corrected particle scheme, which is named the smoothed particle hydrodynamics with corrected symmetric kernel gradient and shifting particle technique (SPH_CS_SP) method. The accuracy and stability of SPH_CS_SP method is first tested by solving Poiseuille flow and Taylor-Green flow. Then the capacity for the SPH_CS_SP method to solve the viscoelastic fluid is verified by the polymer flow through a periodic array of cylinders. Moreover, the convergence of the SPH_CS_SP method is also investigated. Finally, the proposed method is further applied to the 3D viscoelastic jet coiling problem, and the influences of macroscopic parameters on the jet coiling are discussed. The numerical results show that the SPH_CS_SP method has higher accuracy and better stability than the traditional SPH method and other corrected SPH method, and can improve the tensile instability. Project supported by the Natural Science Foundation of Jiangsu Province, China (Grant Nos. BK20130436 and BK20150436) and the Natural Science Foundation of the Higher Education Institutions of Jiangsu Province, China (Grant No. 15KJB110025).

  4. Enabling Lorentz boosted frame particle-in-cell simulations of laser wakefield acceleration in quasi-3D geometry

    NASA Astrophysics Data System (ADS)

    Yu, Peicheng; Xu, Xinlu; Davidson, Asher; Tableman, Adam; Dalichaouch, Thamine; Li, Fei; Meyers, Michael D.; An, Weiming; Tsung, Frank S.; Decyk, Viktor K.; Fiuza, Frederico; Vieira, Jorge; Fonseca, Ricardo A.; Lu, Wei; Silva, Luis O.; Mori, Warren B.

    2016-07-01

    When modeling laser wakefield acceleration (LWFA) using the particle-in-cell (PIC) algorithm in a Lorentz boosted frame, the plasma is drifting relativistically at βb c towards the laser, which can lead to a computational speedup of ∼ γb2 = (1 - βb2)-1. Meanwhile, when LWFA is modeled in the quasi-3D geometry in which the electromagnetic fields and current are decomposed into a limited number of azimuthal harmonics, speedups are achieved by modeling three dimensional (3D) problems with the computational loads on the order of two dimensional r - z simulations. Here, we describe a method to combine the speedups from the Lorentz boosted frame and quasi-3D algorithms. The key to the combination is the use of a hybrid Yee-FFT solver in the quasi-3D geometry that significantly mitigates the Numerical Cerenkov Instability (NCI) which inevitably arises in a Lorentz boosted frame due to the unphysical coupling of Langmuir modes and EM modes of the relativistically drifting plasma in these simulations. In addition, based on the space-time distribution of the LWFA data in the lab and boosted frame, we propose to use a moving window to follow the drifting plasma, instead of following the laser driver as is done in the LWFA lab frame simulations, in order to further reduce the computational loads. We describe the details of how the NCI is mitigated for the quasi-3D geometry, the setups for simulations which combine the Lorentz boosted frame, quasi-3D geometry, and the use of a moving window, and compare the results from these simulations against their corresponding lab frame cases. Good agreement is obtained among these sample simulations, particularly when there is no self-trapping, which demonstrates it is possible to combine the Lorentz boosted frame and the quasi-3D algorithms when modeling LWFA. We also discuss the preliminary speedups achieved in these sample simulations.

  5. Numerical simulation of inhaled aerosol particle deposition within 3D realistic human upper respiratory tract

    NASA Astrophysics Data System (ADS)

    Lin, J.; Fan, J. R.; Zheng, Y. Q.; Hu, G. L.; Pan, D.

    2010-03-01

    Computational fluid dynamics (CFD) simulations of airflow and particle deposition in the upper respiratory tract (URT) were conducted in this paper. Based on the CT (Computerized Tomography) scanned images of a 19-years-old healthy boy, a realistic geometric model of URT from oral cavity to the upper six-generation bronchial is rebuilt. To investigate airflow and particle deposition in the obtained realistic human upper respiratory tract, RNG k-ɛ turbulence model was used to describe the primary flow and particle deposition under three breathing intensity such as 15 L/min, 30 L/min and 60 L/min. The particle is tracked and analyzed in the Lagrangian frame. The velocity fields of airflow under different airflow rates were computed and discussed. In order to study the characteristics of particles movement and the effect of particles diameter on the deposition pattern, eleven kinds of sphere particles with different diameters are selected as research object. The diameters of selected particles as follows: 0.1 μm, 0.5 μm, 1 μm, 2.5 μm, 3 μm, 3.5 μm, 4 μm, 4.5 μm, 5 μm, 6.5 μm and 8 μm. The variation of inhalable particles deposition in realistic human upper respiratory tract with respiratory intensity and particle size was researched and compared. Furthermore, the more real inhalable particles with Rosin-Rammler mass distribution are used to study the effect of particles size. The deposition rate of particles with the different diameter scope in the different part of upper respiratory tract was summarized. The geometrical model based images technology promises to provide more real results of airflow field and particle deposition in the URT.

  6. Some Progress in Large-Eddy Simulation using the 3-D Vortex Particle Method

    NASA Technical Reports Server (NTRS)

    Winckelmans, G. S.

    1995-01-01

    This two-month visit at CTR was devoted to investigating possibilities in LES modeling in the context of the 3-D vortex particle method (=vortex element method, VEM) for unbounded flows. A dedicated code was developed for that purpose. Although O(N(sup 2)) and thus slow, it offers the advantage that it can easily be modified to try out many ideas on problems involving up to N approx. 10(exp 4) particles. Energy spectrums (which require O(N(sup 2)) operations per wavenumber) are also computed. Progress was realized in the following areas: particle redistribution schemes, relaxation schemes to maintain the solenoidal condition on the particle vorticity field, simple LES models and their VEM extension, possible new avenues in LES. Model problems that involve strong interaction between vortex tubes were computed, together with diagnostics: total vorticity, linear and angular impulse, energy and energy spectrum, enstrophy. More work is needed, however, especially regarding relaxation schemes and further validation and development of LES models for VEM. Finally, what works well will eventually have to be incorporated into the fast parallel tree code.

  7. MHD-Epic: Embedded Particle-in-Cell Simulations of Reconnection in Global 3D Extended MHD Simulations

    NASA Astrophysics Data System (ADS)

    Daldorff, L. K. S.; Toth, G.; Borovikov, D.; Gombosi, T. I.; Lapenta, G.

    2014-12-01

    With the new modeling capability in the Space Weather Modeling Framework (SWMF) of embedding an implicit Particle-in-Cell (PIC) model iPIC3D into the BATS-R-US magnetohydrodynamics model (Daldorff et al. 2014, JCP, 268, 236) we are ready to locally handle the full physics of the reconnection and its implications on the full system where globally, away from the reconnection region, a magnetohydrodynamic description is satisfactory. As magnetic reconnection is one of the main drivers in magnetospheric and heliospheric plasma dynamics, the self-consistent description of the electron dynamics in the coupled MHD-EPIC model is well suited for investigating the nature of these systems. We will compare the new embedded MHD-EPIC model with pure MHD and Hall MHD simulations of the Earth's magnetosphere.

  8. Soot particle size modelling in 3D simulations of diesel engine combustion

    NASA Astrophysics Data System (ADS)

    Fraioli, V.; Beatrice, C.; Lazzaro, M.

    2011-12-01

    The present work is focused on multi-dimensional simulations of combustion in diesel engines. The primary objective was to test, in a diesel engine framework, a soot particle size model to represent the carbon particle formation and calculate the corresponding size distribution function. Simulations are performed by means of a parallel version of the KIVA3V numerical code, modified to adopt detailed kinetics reaction mechanisms. A skeletal reaction scheme for n-heptane autoignition has been extended, to include PAH kinetics and carbonaceous particle formation and consumption rates: the full reaction set is made up of 82 gas species and 50 species accounting for the particles, thus the complete reaction scheme comprises 132 species and 2206 reaction steps. Four different engine operative conditions, varying engine speed and load, are taken into account and experimentally tested on a single cylinder diesel engine fuelling pure n-heptane. Computed particle size distribution functions are compared with corresponding measurements at the exhaust, performed by a differential mobility spectrometer. A satisfying agreement between computed and measured combustion profiles is obtained in all the conditions. A reasonable aerosol evolution can be obtained, yet in all the cases the model exhibits the tendency to overestimate the number of particles within the range 5-160 nm. Moreover calculations predict a nucleation mode not detected by the available instrument. According to the simulations, the total number and size of the nascent particles would not depend on the operative conditions, while the features of the larger aggregates distinctly vary with the engine functioning.

  9. 3-D RPIC simulations of relativistic jets: Particle acceleration, magnetic field generation, and emission

    NASA Technical Reports Server (NTRS)

    Nishikawa, K.-I.

    2006-01-01

    Nonthermal radiation observed from astrophysical systems containing (relativistic) jets and shocks, e.g., supernova remnants, active galactic nuclei (AGNs), gamma-ray bursts (GRBs), and Galactic microquasar systems usually have power-law emission spectra. Fermi acceleration is the mechanism usually assumed for the acceleration of particles in astrophysical environments. Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that acceleration occurs within the downstream jet, rather than by the scattering of particles back and forth across the shock as in Fermi acceleration. Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) created in the .shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. We will review recent PIC simulations which show particle acceleration in jets.

  10. Fully 3D Particle-in-Cell Simulation of Double Post-Hole Convolute on PTS Facility

    NASA Astrophysics Data System (ADS)

    Zhao, Hailong; Dong, Ye; Zhou, Haijing; Zou, Wenkang; Institute of Fluid Physics Collaboration; Institute of Applied Physics; Computational Mathematics Collaboration

    2015-11-01

    In order to get better understand of energy transforming and converging process during High Energy Density Physics (HEDP) experiments, fully 3D particle-in-cell (PIC) simulation code NEPTUNE3D was used to provide numerical approach towards parameters which could hardly be acquired through diagnostics. Cubic region (34cm × 34cm × 18cm) including the double post-hole convolute (DPHC) on the primary test stand (PTS) facility was chosen to perform a series of fully 3D PIC simulations, calculating ability of codes were tested and preliminary simulation results about DPHC on PTS facility were discussed. Taking advantages of 3D simulation codes and large-scale parallel computation, massive data (~ 250GB) could be acquired in less than 5 hours and clear process of current transforming and electron emission in DPHC were demonstrated with the help of visualization tools. Cold-chamber tests were performed during which only cathode electron emission was considered without temperature rise or ion emission, current loss efficiency was estimated to be 0.46% ~ 0.48% by comparisons between output magnetic field profiles with or without electron emission. Project supported by the National Natural Science Foundation of China (Grant No. 11205145, 11305015, 11475155).

  11. Particle Acceleration in 3D Magnetic Reconnection

    NASA Astrophysics Data System (ADS)

    Dahlin, J.; Drake, J. F.; Swisdak, M.

    2015-12-01

    Magnetic reconnection is an important driver of energetic particles in phenomena such as magnetospheric storms and solar flares. Using kinetic particle-in-cell (PIC) simulations, we show that the stochastic magnetic field structure which develops during 3D reconnection plays a vital role in particle acceleration and transport. In a 2D system, electrons are trapped in magnetic islands which limits their energy gain. In a 3D system, however, the stochastic magnetic field enables the energetic electrons to access volume-filling acceleration regions and therefore gain energy much more efficiently than in the 2D system. We also examine the relative roles of two important acceleration drivers: parallel electric fields and a Fermi mechanism associated with reflection of charged particles from contracting field lines. We find that parallel electric fields are most important for accelerating low energy particles, whereas Fermi reflection dominates energetic particle production. We also find that proton energization is reduced in the 3D system.

  12. 3-D RPIC Simulations of Relativistic Jets: Particle Acceleration, Magnetic Field Generation, and Emission

    NASA Technical Reports Server (NTRS)

    Nishikawa, K.-I.; Mizuno, Y.; Hardee, P.; Hededal, C. B.; Fishman, G. J.

    2006-01-01

    Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets into ambient plasmas show that acceleration occurs in relativistic shocks. The Weibel instability created in shocks is responsible for particle acceleration, and generation and amplification of highly inhomogeneous, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection in relativistic jets. The "jitter" radiation from deflected electrons has different properties than the synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understand the complex time evolution and spectral structure in relativistic jets and gamma-ray bursts. We will present recent PIC simulations which show particle acceleration and magnetic field generation. We will also calculate associated self-consistent emission from relativistic shocks.

  13. 3D particle simulation of beams using the WARP code: Transport around bends

    SciTech Connect

    Friedman, A.; Grote, D.P.; Callahan, D.A.; Langdon, A.B. ); Haber, I. )

    1990-11-30

    WARP is a discrete-particle simulation program which was developed for studies of space charge dominated ion beams. It combines features of an accelerator code and a particle-in-cell plasma simulation. The code architecture, and techniques employed to enhance efficiency, are briefly described. Current applications are reviewed. In this paper we emphasize the physics of transport of three-dimensional beams around bends. We present a simple bent-beam PIC algorithm. Using this model, we have followed a long, thin beam around a bend in a simple racetrack system (assuming straight-pipe self-fields). Results on beam dynamics are presented; no transverse emittance growth (at mid-pulse) is observed. 11 refs., 5 figs.

  14. 3D Boltzmann Simulation of the Io's Plasma Environment with Adaptive Mesh and Particle Refinement

    NASA Astrophysics Data System (ADS)

    Lipatov, A. S.; Combi, M. R.

    2002-12-01

    The global dynamics of the ionized and neutral components in the environment of Io plays an important role in the interaction of Jupiter's corotating magnetospheric plasma with Io [Combi et al., 2002; 1998; Kabin et al., 2001]. The stationary simulation of this problem was done in the MHD [Combi et al., 1998; Linker et al, 1998; Kabin et al., 2001] and the electrodynamic [Saur et al., 1999] approaches. In this report, we develop a method of kinetic ion-neutral simulation, which is based on a multiscale adaptive mesh, particle and algorithm refinement. This method employs the fluid description for electrons whereas for ions the drift-kinetic and particle approaches are used. This method takes into account charge-exchange and photoionization processes. The first results of such simulation of the dynamics of ions in the Io's environment are discussed in this report. ~ M R Combi et al., J. Geophys. Res., 103, 9071, 1998. M R Combi, T I Gombosi, K Kabin, Atmospheres in the Solar System: Comparative\\ Aeronomy. Geophys. Monograph Series, 130, 151, 2002. K Kabin et al., Planetary and Space Sci., 49, 337, 2001. J A Linker et al., J. Geophys. Res., 103(E9), 19867, 1998. J Saur et al., J. Geophys. Res., 104, 25105, 1999.

  15. Particle sedimentation in curved tubes: A 3D simulation and optimization for treatment of vestibular vertigo

    NASA Astrophysics Data System (ADS)

    White, Brian; Squires, Todd M.; Hain, Timothy C.; Stone, Howard A.

    2003-11-01

    Benign paroxysmal positional vertigo (BPPV) is a mechanical disorder of the vestibular system where micron-size crystals abnormally drift into the semicircular canals of the inner ear that sense angular motion of the head. Sedimentation of these crystals causes sensation of motion after true head motion has stopped: vertigo results. The usual clinical treatment is through a series of head maneuvers designed to move the particles into a less sensitive region of the canal system. We present a three-dimensional model to simulate treatment of BPPV by determining the complete hydrodynamic motion of the particles through the course of a therapeutic maneuver while using a realistic representation of the actual geometry. Analyses of clinical maneuvers show the parameter range for which they are effective, and indicate inefficiencies in current practice. In addition, an optimization process determines the most effective head maneuver, which significantly differs from those currently in practice.

  16. Grid dependent noise and entropy growth in anisotropic 3d particle-in-cell simulation of high intensity beams

    NASA Astrophysics Data System (ADS)

    Hofmann, I.; Boine-Frankenheim, O.

    2014-12-01

    The numerical noise inherent to particle-in-cell (PIC) simulation of 3d anisotropic high intensity bunched beams in periodic focusing is compared with the analytical model by Struckmeier [Part. Accel. 45, 229 (1994)]. The latter assumes that entropy growth can be related to Markov type stochastic processes due to temperature anisotropy and the artificial "collisions" caused by using macro-particles and calculating the space charge effect. The PIC simulations are carried out with the tracewin code widely used for high intensity beam simulation. The resulting noise can lead to growth of the six-dimensional rms emittance. The logarithm of the latter is shown to qualify as rms-based entropy. We confirm the dependence of this growth on the bunch temperature anisotropy as predicted by Struckmeier. However, we also find a grid and focusing dependent component of noise not predicted by Struckmeier. Although commonalities exist with well-established models for collision effects in PIC-simulation of extended plasmas, a distinctive feature is the presence of a periodic focusing potential, wherein the beam one-component plasma extends only over relatively few Debye lengths. Our findings are applied in particular to noise in high current linac beam simulation, where they help for optimization of the balance between the number of simulation particles and the grid resolution.

  17. 3D Electromagnetic Particle-in-Cell simulations of the solar wind interaction with lunar magnetic anomalies

    NASA Astrophysics Data System (ADS)

    Deca, J.; Lapenta, G.; Divin, A. V.; Lembege, B.; Markidis, S.

    2013-12-01

    Unlike the Earth and Mercury, our Moon has no global magnetic field and is therefore not shielded from the impinging solar wind by a magnetosphere. However, lunar magnetic field measurements made by the Apollo missions provided direct evidence that the Moon has regions of small-scale crustal magnetic fields, ranging up to a few 100km in scale size with surface magnetic field strengths up to hundreds of nanoTeslas. More recently, the Lunar Prospector spacecraft has provided high-resolution observations allowing to construct magnetic field maps of the entire Moon, confirming the earlier results from Apollo, but also showing that the lunar plasma environment is much richer than earlier believed. Typically the small-scale magnetic fields are non-dipolar and rather tiny compared to the lunar radius and mainly clustered on the far side of the moon. Using iPic3D we present the first 3D fully kinetic and electromagnetic Particle-in-Cell simulations of the solar wind interaction with lunar magnetic anomalies. We study the behaviour of a dipole model with variable surface magnetic field strength under changing solar wind conditions and confirm that lunar crustal magnetic fields may indeed be strong enough to stand off the solar wind and form a mini-magnetosphere, as suggested by MHD and hybrid simulations and spacecraft observations. 3D-PIC simulations reveal to be very helpful to analyze the diversion/braking of the particle flux and the characteristics of the resulting particles accumulation. The particle flux to the surface is significantly reduced at the magnetic anomaly, surrounded by a region of enhanced density due to the magnetic mirror effect. Second, the ability of iPic3D to resolve all plasma components (heavy ions, protons and electrons) allows to discuss in detail the electron physics leading to the highly non-adiabatic interactions expected as well as the implications for solar wind shielding of the lunar surface, depending on the scale size (solar wind protons

  18. Guiding-centre and full-Lorentz orbit solving in 3D magnetic coordinates for fast particle simulations

    NASA Astrophysics Data System (ADS)

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

    2014-10-01

    Designed to accurately solve the motion of energetic particles in the presence of 3D magnetic fields, the VENUS-LEVIS code leans on a non-canonical general coordinate Lagrangian formulation of the equations of motion. It switches between full-orbit particle following and guiding-centre tracing by verifying the perpendicular variation of magnetic vector field, not only including gradients and curvature terms but also the shearing of field-lines. The criteria is particularly relevant for the study of fast ion redistribution in the kinked core of hybrid plasmas, where the compression of flux-surfaces against the axisymmetric outer mantle creates strongly varying magnetic field-lines and large parallel currents. Slowing-down simulations of NBI fast ions show that co-passing particles helically align in the opposite side of the plasma deformation whereas counter-passing particles are barely affected by the kinked structure. Results are compared with experimental neutron camera traces and FIDA measurements during long-lived modes (LLM).

  19. Particle entry through sash in the magnetopause with a dawndard IMF as simulated by a 3-D EM particle code

    NASA Astrophysics Data System (ADS)

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

    2003-12-01

    We report a new progress in the long-term effort to represent the global interaction of the solar wind with the Earth's magnetosphere using a three-dimensional electromagnetic particle code with the improved resolutions using the HPF Tristan code. After a quasi-steady state is established with an unmagnetized solar wind we gradually switch on a northward interplanetary magnetic field (IMF), which causes a magnetic reconnection at the nightside cusps and the magnetosphere to be depolarized. In the case that the northward IMF is switched gradually to dawnward, there is no signature of reconnection in the near-Earth magnetotail as in the case with the southward turning. On the contrary analysis of magnetic fields in the magnetopause confirms a signature of magnetic reconnection at both the dawnside and duskside. And the plasma sheet in the near-Earth magnetotail clearly thins as in the case of southward turning. Arrival of dawnward IMF to the magnetopause creates a reconnection groove which cause particle entry into the deep region of the magnetosphere via field lines that go near the magnetopause. This deep connection is more fully recognized tailward of Earth. The flank weak-field fan joins onto the plasma sheet and the current sheet to form a geometrical feature called the cross-tail S that structurally integrates the magnetopause and the tail interior. This structure contributes to direct plasma entry between the magnetosheath to the inner magnetosphere and plasma sheet, in which the entry process heats the magnetosheath plasma to plasma sheet temperatures. These phenomena have been found by Cluster observations. Further investigation with Cluster observations will provide new insights for unsolved problems such as hot flow anomalies (HFAs), substorms, and storm-substorm relationship. 3-D movies with sash structure will be presented at the meeting.

  20. 3-D Full-kinetic Simulations of the Solar Wind Interaction with Lunar Magnetic Anomalies: Particle Behaviour

    NASA Astrophysics Data System (ADS)

    Deca, J.; Divin, A. V.; Wang, X.; Lembege, B.; Markidis, S.; Lapenta, G.; Horanyi, M.

    2015-12-01

    We present three-dimensional full-kinetic electromagnetic simulations of the solar wind interaction with lunar crustal magnetic anomalies (LMAs). Using the implicit particle-in-cell code iPic3D, we confirm that LMAs may indeed be strong enough to stand off the solar wind from directly impacting the lunar surface forming a mini-magnetosphere, as suggested by spacecraft observations and theory. In contrast to earlier MHD and hybrid simulations, the full-kinetic nature of iPic3D allows to self-consistently investigate space charge effects, and in particular the electron dynamics dominating the near-surface lunar plasma environment. We describe the general mechanism of the interaction of both a horizontal and vertical dipole model embedded just below the lunar surface focussing on the ion and electron kinetic behaviour of the system. It is shown that the configurations are largely dominated by electron motion, because the LMA scale size is small with respect to the gyro-radius of the solar wind ions. The formation of mini-magnetospheres is an electrostatic effect. Additionally, we discuss typical particle trajectories as well as complete particle distribution functions covering thermal and suprathermal energies, within the interaction region and on viable spacecraft altitudes. Our work opens new frontiers of research toward a deeper understanding of LMAs and is ideally suited to be compared with field or particle observations from spacecraft such as Kaguya (SELENE), Lunar Prospector or ARTEMIS. The ability to evaluate the implications for future lunar exploration as well as lunar science in general hinges on a better understanding of LMAs.This research has received funding from the European Commission's FP7 Program with the grant agreement EHEROES (project 284461, www.eheroes.eu). The simulations were conducted on the computational resources provided by the PRACE Tier-0 project 2013091928 (SuperMUC). This research was supported by the Swedish National Space Board

  1. 3-D Relativistic MHD Simulations

    NASA Astrophysics Data System (ADS)

    Nishikawa, K.-I.; Frank, J.; Koide, S.; Sakai, J.-I.; Christodoulou, D. M.; Sol, H.; Mutel, R. L.

    1998-12-01

    We present 3-D numerical simulations of moderately hot, supersonic jets propagating initially along or obliquely to the field lines of a denser magnetized background medium with Lorentz factors of W = 4.56 and evolving in a four-dimensional spacetime. The new results are understood as follows: Relativistic simulations have consistently shown that these jets are effectively heavy and so they do not suffer substantial momentum losses and are not decelerated as efficiently as their nonrelativistic counterparts. In addition, the ambient magnetic field, however strong, can be pushed aside with relative ease by the beam, provided that the degrees of freedom associated with all three spatial dimensions are followed self-consistently in the simulations. This effect is analogous to pushing Japanese ``noren'' or vertical Venetian blinds out of the way while the slats are allowed to bend in 3-D space rather than as a 2-D slab structure.

  2. 3-D Relativistic MHD Simulations

    NASA Astrophysics Data System (ADS)

    Nishikaw, K.-I.; Frank, J.; Christodoulou, D. M.; Koide, S.; Sakai, J.-I.; Sol, H.; Mutel, R. L.

    1998-12-01

    We present 3-D numerical simulations of moderately hot, supersonic jets propagating initially along or obliquely to the field lines of a denser magnetized background medium with Lorentz factors of W=4.56 and evolving in a four-dimensional spacetime. The new results are understood as follows: Relativistic simulations have consistently shown that these jets are effectively heavy and so they do not suffer substantial momentum losses and are not decelerated as efficiently as their nonrelativistic counterparts. In addition, the ambient magnetic field, however strong, can be pushed aside with relative ease by the beam, provided that the degrees of freedom associated with all three spatial dimensions are followed self-consistently in the simulations. This effect is analogous to pushing Japanese ``noren'' or vertical Venetian blinds out of the way while the slats are allowed to bend in 3-D space rather than as a 2-D slab structure. We also simulate jets with the more realistic initial conditions for injecting jets for helical mangetic field, perturbed density, velocity, and internal energy, which are supposed to be caused in the process of jet generation. Three possible explanations for the observed variability are (i) tidal disruption of a star falling into the black hole, (ii) instabilities in the relativistic accretion disk, and (iii) jet-related PRocesses. New results will be reported at the meeting.

  3. Prediction of the effect of an extreme solar event on the Martian environment using a 3-D, self consistent hybrid model supplemented by test particle simulations.

    NASA Astrophysics Data System (ADS)

    McKenna-Lawlor, Susan; Kallio, Esa; Dyadechkin, Sergey; Jarvinen, Riku; Janhunen, Pekka

    2010-05-01

    Energetic particle data recorded by the SLED and LET instruments aboard the Phobos-2 spacecraft while in circular orbit about Mars showed the presence from 6 > 26 March, 1989, in association with an extreme solar event, of intense ambient particle radiation (> 30 MeV) punctuated by traveling interplanetary shocks. The response of the Martian environment to the March 1989 solar disturbances is modeled using a 3-D, self-consistent, hybrid model (HYB-Mars) supplemented by test particle simulations. In HYB-Mars ions are modeled as particles while electrons form a massless, charge neutralizing, fluid. The magnetic and electric fields present during the March 1989 activity are each derived from HYB- Mars while the high energy ion populations are analyzed using test particle simulations. Finally, the predictions of the model are validated by comparing the simulated properties of the disturbed Martian environment with the in situ observations.

  4. 3D Finite Element Analysis of Particle-Reinforced Aluminum

    NASA Technical Reports Server (NTRS)

    Shen, H.; Lissenden, C. J.

    2002-01-01

    Deformation in particle-reinforced aluminum has been simulated using three distinct types of finite element model: a three-dimensional repeating unit cell, a three-dimensional multi-particle model, and two-dimensional multi-particle models. The repeating unit cell model represents a fictitious periodic cubic array of particles. The 3D multi-particle (3D-MP) model represents randomly placed and oriented particles. The 2D generalized plane strain multi-particle models were obtained from planar sections through the 3D-MP model. These models were used to study the tensile macroscopic stress-strain response and the associated stress and strain distributions in an elastoplastic matrix. The results indicate that the 2D model having a particle area fraction equal to the particle representative volume fraction of the 3D models predicted the same macroscopic stress-strain response as the 3D models. However, there are fluctuations in the particle area fraction in a representative volume element. As expected, predictions from 2D models having different particle area fractions do not agree with predictions from 3D models. More importantly, it was found that the microscopic stress and strain distributions from the 2D models do not agree with those from the 3D-MP model. Specifically, the plastic strain distribution predicted by the 2D model is banded along lines inclined at 45 deg from the loading axis while the 3D model prediction is not. Additionally, the triaxial stress and maximum principal stress distributions predicted by 2D and 3D models do not agree. Thus, it appears necessary to use a multi-particle 3D model to accurately predict material responses that depend on local effects, such as strain-to-failure, fracture toughness, and fatigue life.

  5. 3D Magnetron simulation with CST STUDIO SUITE

    SciTech Connect

    Balk, Monika C.

    2011-07-01

    The modeling of magnetrons compared to other tubes is more difficult since it requires 3D modeling rather than a 2D investigation. This is not only due to the geometry which can include complicated details to be modeled in 3D but also due to the interaction process itself. The electric field, magnetic field and particle movement span a 3D space. In this paper 3D simulations of a strapped magnetron with CSTSTUDIO SUITE{sup TM} are presented. (author)

  6. Self-consistent Nonlinear Analysis and 3D Particle-In-Cell Simulation of a W-band Gyro-TWT

    NASA Astrophysics Data System (ADS)

    Tang, Yong; Luo, Yong; Xu, Yong; Yan, Ran

    2014-10-01

    The self-consistent nonlinear analysis and CST 3D particle-in-cell (PIC) simulation of a W-band gyrotron traveling wave tube (gyro-TWT) are presented in this paper. Both the simulation results of the two codes are excellent agreement with each other. The gyro-TWT loaded with periodic lossy dielectric in the circuit for suppressing potential spurious oscillations. It is driven by a 70kV, 10A gyrating electron beam with velocity ratio of 1.0. PIC simulation results are: the maximum peak output power of 198kW, statured gain of 62.3dB and efficiency of 28.3% at 92.5GHz. Only the operating mode TE 01 is observed in the CST 3D simulation and the potential competing backward wave oscillations are effectively suppressed. The CST simulation also predicts that the device works stably under the condition of the beam current lower than 14A and B 0 /B g lower than 1.05. The simulated bandwidth with peak power greater than 100kW is 6.8GHz without axial velocity spread, and 4.1GHz with 6% axial velocity spread.

  7. Self consistent particles dynamics in/out of the cusp region by using back tracking technics; a global 3D PIC simulation approach

    NASA Astrophysics Data System (ADS)

    Esmaeili, A.; Cai, D.; Lembege, B.; Nishikawa, K.

    2013-12-01

    Large scale three dimensionbal PIC (particle in cell) simulations are presently used in order to analyze the global solar wind-terrestrial magnetosphere intreraction within a full self-consistent approach, and where both electrons and ions are treated as an assembly of individual particles. This 3D kinetic approach allows us to analyze in particular the dynamics and the fine structures of the cusp region when including self consistently not only its whole neighborhood (in the terrestrial magnetosphere) but also the impact of the solar wind and the interplanetary field (IMF) features. Herein, we focuss our attention on the cusp region and in particular on the acceleration and precipitation of particles (both ions and electrons) within the cusp. In present simulations, the IMF is chosen northward, (i.e. where the X -reconnection region is just above the cusp, in the meridian plane). Back-trackings of self-consistent particles are analyzed in details in order to determine (i) which particles (just above the cusp) are precipitated deeply into the cusp, (ii) which populations are injected from the cusp into the nearby tail, (iii) where the particles suffer the largest energisation along their self-consistent trajectories, (iv) where these populations accumulate, and (v) where the most energetic particles are originally coming from. This approach allows to make a traking of particles within the scenario "solar wind-magnetosheath- cusp -nearbytail"; moreover it strongly differs from the standard test particles technics and allows to provide informations not accessible when using full MHD approach. Keywords: Tracing Particles, Particle In Cell (PIC) simulation, double cusp, test particles method, IMF, Solar wind, Magnetosphere

  8. Non-isothermal 3D SDPD Simulations

    NASA Astrophysics Data System (ADS)

    Yang, Jun; Potami, Raffaele; Gatsonis, Nikolaos

    2012-11-01

    The study of fluids at micro and nanoscale requires new modeling and computational approaches. Smooth Particle Dissipative Dynamics (SDPD) is a mesh-free method that provides a bridge between the continuum equations of hydrodynamics embedded in the Smooth Particle Hydrodynamics approach and the molecular nature embedded in the DPD approach. SDPD is thermodynamically consistent, does not rely on arbitrary coefficients for its thermostat, involves realistic transport coefficients, and includes fluctuation terms. SDPD is implemented in our work for arbitrary 3D geometries with a methodology to model solid wall boundary conditions. We present simulations for isothermal flows for verification of our approach. The entropy equation is implemented with a velocity-entropy Verlet integration algorithm Flows with heat transfer are simulated for verification of the SDPD. We present also the self-diffusion coefficient derived from SDPD simulations for gases and liquids. Results show the scale dependence of self-diffusion coefficient on SDPD particle size. Computational Mathematics Program of the Air Force Office of Scientific Research under grant/contract number FA9550-06-1-0236.

  9. 3D Simulation: Microgravity Environments and Applications

    NASA Technical Reports Server (NTRS)

    Hunter, Steve L.; Dischinger, Charles; Estes, Samantha; Parker, Nelson C. (Technical Monitor)

    2001-01-01

    Most, if not all, 3-D and Virtual Reality (VR) software programs are designed for one-G gravity applications. Space environments simulations require gravity effects of one one-thousandth to one one-million of that of the Earth's surface (10(exp -3) - 10(exp -6) G), thus one must be able to generate simulations that replicate those microgravity effects upon simulated astronauts. Unfortunately, the software programs utilized by the National Aeronautical and Space Administration does not have the ability to readily neutralize the one-G gravity effect. This pre-programmed situation causes the engineer or analysis difficulty during micro-gravity simulations. Therefore, microgravity simulations require special techniques or additional code in order to apply the power of 3D graphic simulation to space related applications. This paper discusses the problem and possible solutions to allow microgravity 3-D/VR simulations to be completed successfully without program code modifications.

  10. 3D View of Mars Particle

    NASA Technical Reports Server (NTRS)

    2008-01-01

    This is a 3D representation of the pits seen in the first Atomic Force Microscope, or AFM, images sent back from NASA's Phoenix Mars Lander. Red represents the highest point and purple represents the lowest point.

    The particle in the upper left corner shown at the highest magnification ever seen from another world is a rounded particle about one micrometer, or one millionth of a meter, across. It is a particle of the dust that cloaks Mars. Such dust particles color the Martian sky pink, feed storms that regularly envelop the planet and produce Mars' distinctive red soil.

    The particle was part of a sample informally called 'Sorceress' delivered to the AFM on the 38th Martian day, or sol, of the mission (July 2, 2008). The AFM is part of Phoenix's microscopic station called MECA, or the Microscopy, Electrochemistry, and Conductivity Analyzer.

    The AFM was developed by a Swiss-led consortium, with Imperial College London producing the silicon substrate that holds sampled particles.

    The Phoenix Mission is led by the University of Arizona, Tucson, on behalf of NASA. Project management of the mission is by NASA's Jet Propulsion Laboratory, Pasadena, Calif. Spacecraft development is by Lockheed Martin Space Systems, Denver.

  11. A 3D immersed finite element method with non-homogeneous interface flux jump for applications in particle-in-cell simulations of plasma-lunar surface interactions

    NASA Astrophysics Data System (ADS)

    Han, Daoru; Wang, Pu; He, Xiaoming; Lin, Tao; Wang, Joseph

    2016-09-01

    Motivated by the need to handle complex boundary conditions efficiently and accurately in particle-in-cell (PIC) simulations, this paper presents a three-dimensional (3D) linear immersed finite element (IFE) method with non-homogeneous flux jump conditions for solving electrostatic field involving complex boundary conditions using structured meshes independent of the interface. This method treats an object boundary as part of the simulation domain and solves the electric field at the boundary as an interface problem. In order to resolve charging on a dielectric surface, a new 3D linear IFE basis function is designed for each interface element to capture the electric field jump on the interface. Numerical experiments are provided to demonstrate the optimal convergence rates in L2 and H1 norms of the IFE solution. This new IFE method is integrated into a PIC method for simulations involving charging of a complex dielectric surface in a plasma. A numerical study of plasma-surface interactions at the lunar terminator is presented to demonstrate the applicability of the new method.

  12. Fractal and fractional calculus to model hydrological processes with application to particle-based 2D and 3D landslide simulation

    NASA Astrophysics Data System (ADS)

    Martelloni, Gianluca; Bagnoli, Franco; Di Cintio, Pierfrancesco

    2015-04-01

    We integrate existing soil infiltration modeling with particle based methods in order to simulate two and three-dimensional setups of triggered landslides. Commonly, the infiltration models are based on continuum schemes (e.g. Eulerian approach) by means of which it is possible to define the field of the pore pressure within a soil. By contrast, the particle based methods follow a Lagrangian scheme that allows one to identify the particle trajectories and their dynamical properties. In this work, in order to simulate the triggering mechanism, we apply the classical, fractal and fractional Richards equations and the Mohr-Coulomb failure criterion, adapted to the molecular dynamics technique. In our scheme the (local) positive pore pressure is simply implemented as a perturbation of the rest state of each grain. Therefore, the pore pressure function can be interpreted as a time-space dependent scalar field acting on each particle. To initialize the system we generate, using a molecular dynamics based algorithm, a mechanically stable disk (2D) or sphere (3D) packing simulating the consolidated soil. In this way, we can built the micro and macro pore structure related to different infiltration time scales. The inter-particle interactions are modeled with a Lennard-Jones like potential. The particle positions are updated in time, after and during a rainfall, with standard molecular dynamics. We analyze the sensitivity of the model with respect to the variation of some parameters such as hydraulic conductivity, cohesion, slope and friction angle, soil depth and fractional order of the generalized infiltration model. In addition, we consider both regular and random particle configurations. The results of our simulations are found to be in agreement with real landslides. In particular, the mean velocity patterns of the simulated landslides appear extremely similar to the observed ones. Moreover, it is possible to apply the method of the inverse surface displacement

  13. M3D project for simulation studies of plasmas

    SciTech Connect

    Park, W.; Belova, E.V.; Fu, G.Y.; Strauss, H.R.; Sugiyama, L.E.

    1998-12-31

    The M3D (Multi-level 3D) project carries out simulation studies of plasmas of various regimes using multi-levels of physics, geometry, and mesh schemes in one code package. This paper and papers by Strauss, Sugiyama, and Belova in this workshop describe the project, and present examples of current applications. The currently available physics models of the M3D project are MHD, two-fluids, gyrokinetic hot particle/MHD hybrid, and gyrokinetic particle ion/two-fluid hybrid models. The code can be run with both structured and unstructured meshes.

  14. Development of the 3D Parallel Particle-In-Cell Code IMPACT to Simulate the Ion Beam Transport System of VENUS (Abstract)

    SciTech Connect

    Qiang, J.; Leitner, D.; Todd, D.S.; Ryne, R.D.

    2005-03-15

    The superconducting ECR ion source VENUS serves as the prototype injector ion source for the Rare Isotope Accelerator (RIA) driver linac. The RIA driver linac requires a great variety of high charge state ion beams with up to an order of magnitude higher intensity than currently achievable with conventional ECR ion sources. In order to design the beam line optics of the low energy beam line for the RIA front end for the wide parameter range required for the RIA driver accelerator, reliable simulations of the ion beam extraction from the ECR ion source through the ion mass analyzing system are essential. The RIA low energy beam transport line must be able to transport intense beams (up to 10 mA) of light and heavy ions at 30 keV.For this purpose, LBNL is developing the parallel 3D particle-in-cell code IMPACT to simulate the ion beam transport from the ECR extraction aperture through the analyzing section of the low energy transport system. IMPACT, a parallel, particle-in-cell code, is currently used to model the superconducting RF linac section of RIA and is being modified in order to simulate DC beams from the ECR ion source extraction. By using the high performance of parallel supercomputing we will be able to account consistently for the changing space charge in the extraction region and the analyzing section. A progress report and early results in the modeling of the VENUS source will be presented.

  15. Development of the 3D Parallel Particle-In-Cell Code IMPACT to Simulate the Ion Beam Transport System of VENUS (Abstract)

    NASA Astrophysics Data System (ADS)

    Qiang, J.; Leitner, D.; Todd, D. S.; Ryne, R. D.

    2005-03-01

    The superconducting ECR ion source VENUS serves as the prototype injector ion source for the Rare Isotope Accelerator (RIA) driver linac. The RIA driver linac requires a great variety of high charge state ion beams with up to an order of magnitude higher intensity than currently achievable with conventional ECR ion sources. In order to design the beam line optics of the low energy beam line for the RIA front end for the wide parameter range required for the RIA driver accelerator, reliable simulations of the ion beam extraction from the ECR ion source through the ion mass analyzing system are essential. The RIA low energy beam transport line must be able to transport intense beams (up to 10 mA) of light and heavy ions at 30 keV. For this purpose, LBNL is developing the parallel 3D particle-in-cell code IMPACT to simulate the ion beam transport from the ECR extraction aperture through the analyzing section of the low energy transport system. IMPACT, a parallel, particle-in-cell code, is currently used to model the superconducting RF linac section of RIA and is being modified in order to simulate DC beams from the ECR ion source extraction. By using the high performance of parallel supercomputing we will be able to account consistently for the changing space charge in the extraction region and the analyzing section. A progress report and early results in the modeling of the VENUS source will be presented.

  16. Comparison of quasi-3D and full-3D laser wakefield PIC simulations using azimuthal mode decomposition

    NASA Astrophysics Data System (ADS)

    Dalichaouch, Thamine; Yu, Peicheng; Davidson, Asher; Mori, Warren; Vieira, Jorge; Fonseca, Ricardo

    2015-11-01

    Laser wakefield acceleration (LWFA) has attracted a lot of interest as a possible compact particle accelerator. However, 3D simulations of plasma-based accelerators are computationally intensive, sometimes taking millions of core hours on today's computers. A quasi-3D particle-In-cell (PIC) approach has been developed to take advantage of azimuthal symmetry in LWFA (and PWFA) simulations by using a particle-in-cell description in r-z and a Fourier description in φ. Quasi-3D simulations of LWFA are computationally more efficient and faster than Full-3D simulations because only first few azimuthal harmonics are needed to capture the physics of the problem. We have developed a cylindrical mode decomposition diagnostic for 3D Cartesian geometry simulations to analyze the agreement between full-3D and quasi-3D PIC simulations of laser and beam-plasma interactions. The diagnostic interpolates field data from Full-3D PIC simulations onto an irregular cylindrical grid (r , φ , z). A Fourier decomposition is then performed on the interpolated 3D simulation data along the azimuthal direction. This diagnostic has the added advantage of separating out the wakefields from the laser field. Preliminary results for this diagnostic of LWFA and PWFA simulations with symmetric and nearly symmetric spot sizes as well as of laser-plasma interactions using lasers with orbital angular momentum (higher order Laguerre-Gaussian modes) will be presented.

  17. Loop heating by D.C. electric current and electromagnetic wave emissions simulated by 3-D EM particle zone

    NASA Technical Reports Server (NTRS)

    Sakai, J. I.; Zhao, J.; Nishikawa, K.-I.

    1994-01-01

    We have shown that a current-carrying plasma loop can be heated by magnetic pinch driven by the pressure imbalance between inside and outside the loop, using a 3-dimensional electromagnetic (EM) particle code. Both electrons and ions in the loop can be heated in the direction perpendicular to the ambient magnetic field, therefore the perpendicular temperature can be increased about 10 times compared with the parallel temperature. This temperature anisotropy produced by the magnetic pinch heating can induce a plasma instability, by which high-frequency electromagnetic waves can be excited. The plasma current which is enhanced by the magnetic pinch can also excite a kinetic kink instability, which can heat ions perpendicular to the magnetic field. The heating mechanism of ions as well as the electromagnetic emission could be important for an understanding of the coronal loop heating and the electromagnetic wave emissions from active coronal regions.

  18. Large-eddy simulations of 3D Taylor-Green vortex: comparison of Smoothed Particle Hydrodynamics, Lattice Boltzmann and Finite Volume methods

    NASA Astrophysics Data System (ADS)

    Kajzer, A.; Pozorski, J.; Szewc, K.

    2014-08-01

    In the paper we present Large-eddy simulation (LES) results of 3D Taylor- Green vortex obtained by the three different computational approaches: Smoothed Particle Hydrodynamics (SPH), Lattice Boltzmann Method (LBM) and Finite Volume Method (FVM). The Smagorinsky model was chosen as a subgrid-scale closure in LES for all considered methods and a selection of spatial resolutions have been investigated. The SPH and LBM computations have been carried out with the use of the in-house codes executed on GPU and compared, for validation purposes, with the FVM results obtained using the open-source CFD software OpenFOAM. A comparative study in terms of one-point statistics and turbulent energy spectra shows a good agreement of LES results for all methods. An analysis of the GPU code efficiency and implementation difficulties has been made. It is shown that both SPH and LBM may offer a significant advantage over mesh-based CFD methods.

  19. INCORPORATING DYNAMIC 3D SIMULATION INTO PRA

    SciTech Connect

    Steven R Prescott; Curtis Smith

    2011-07-01

    Through continued advancement in computational resources, development that was previously done by trial and error production is now performed through computer simulation. These virtual physical representations have the potential to provide accurate and valid modeling results and are being used in many different technical fields. Risk assessment now has the opportunity to use 3D simulation to improve analysis results and insights, especially for external event analysis. By using simulations, the modeler only has to determine the likelihood of an event without having to also predict the results of that event. The 3D simulation automatically determines not only the outcome of the event, but when those failures occur. How can we effectively incorporate 3D simulation into traditional PRA? Most PRA plant modeling is made up of components with different failure modes, probabilities, and rates. Typically, these components are grouped into various systems and then are modeled together (in different combinations) as a “system” with logic structures to form fault trees. Applicable fault trees are combined through scenarios, typically represented by event tree models. Though this method gives us failure results for a given model, it has limitations when it comes to time-based dependencies or dependencies that are coupled to physical processes which may themselves be space- or time-dependent. Since, failures from a 3D simulation are naturally time related, they should be used in that manner. In our simulation approach, traditional static models are converted into an equivalent state diagram representation with start states, probabilistic driven movements between states and terminal states. As the state model is run repeatedly, it converges to the same results as the PRA model in cases where time-related factors are not important. In cases where timing considerations are important (e.g., when events are dependent upon each other), then the simulation approach will typically

  20. HIGH-TEMPERATURE PROCESSING OF SOLIDS THROUGH SOLAR NEBULAR BOW SHOCKS: 3D RADIATION HYDRODYNAMICS SIMULATIONS WITH PARTICLES

    SciTech Connect

    Boley, A. C.; Morris, M. A.; Desch, S. J.

    2013-10-20

    A fundamental, unsolved problem in solar system formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks have been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals and planetary embryos on eccentric orbits can produce bow shocks as they move supersonically through the disk gas, and are one possible source of chondrule-melting shocks. We investigate chondrule formation in bow shocks around planetoids through three-dimensional radiation hydrodynamics simulations. A new radiation transport algorithm that combines elements of flux-limited diffusion and Monte Carlo methods is used to capture the complexity of radiative transport around bow shocks. An equation of state that includes the rotational, vibrational, and dissociation modes of H{sub 2} is also used. Solids are followed directly in the simulations and their thermal histories are recorded. Adiabatic expansion creates rapid cooling of the gas, and tail shocks behind the embryo can cause secondary heating events. Radiative transport is efficient, and bow shocks around planetoids can have luminosities ∼few× 10{sup –8} L{sub ☉}. While barred and radial chondrule textures could be produced in the radiative shocks explored here, porphyritic chondrules may only be possible in the adiabatic limit. We present a series of predicted cooling curves that merit investigation in laboratory experiments to determine whether the solids produced by bow shocks are represented in the meteoritic record by chondrules or other solids.

  1. 3D laser traking of a particle in 3DFM

    NASA Astrophysics Data System (ADS)

    Desai, Kalpit; Welch, Gregory; Bishop, Gary; Taylor, Russell; Superfine, Richard

    2003-11-01

    The principal goal of 3D tracking in our home-built 3D Magnetic Force Microscope is to monitor movement of the particle with respect to laser beam waist and keep the particle at the center of laser beam. The sensory element is a Quadrant Photo Diode (QPD) which captures scattering of light caused by particle motion with bandwidth up to 40 KHz. XYZ translation stage is the driver element which moves particle back in the center of the laser with accuracy of couple of nanometers and with bandwidth up to 300 Hz. Since our particles vary in size, composition and shape, instead of using a priori model we use standard system identification techniques to have optimal approximation to the relationship between particle motion and QPD response. We have developed position feedback control system software that is capable of 3-dimensional tracking of beads that are attached to cilia on living cells which are beating at up to 15Hz. We have also modeled the control system of instrument to simulate performance of 3D particle tracking for different experimental conditions. Given operational level of nanometers, noise poses a great challenge for the tracking system. We propose to use stochastic control theory approaches to increase robustness of tracking.

  2. 3D MHD Simulations of Tokamak Disruptions

    NASA Astrophysics Data System (ADS)

    Woodruff, Simon; Stuber, James

    2014-10-01

    Two disruption scenarios are modeled numerically by use of the CORSICA 2D equilibrium and NIMROD 3D MHD codes. The work follows the simulations of pressure-driven modes in DIII-D and VDEs in ITER. The aim of the work is to provide starting points for simulation of tokamak disruption mitigation techniques currently in the CDR phase for ITER. Pressure-driven instability growth rates previously observed in simulations of DIIID are verified; Halo and Hiro currents produced during vertical displacements are observed in simulations of ITER with implementation of resistive walls in NIMROD. We discuss plans to exercise new code capabilities and validation.

  3. Crashworthiness simulations with DYNA3D

    SciTech Connect

    Schauer, D.A.; Hoover, C.G.; Kay, G.J.; Lee, A.S.; De Groot, A.J.

    1996-04-01

    Current progress in parallel algorithm research and applications in vehicle crash simulation is described for the explicit, finite element algorithms in DYNA3D. Problem partitioning methods and parallel algorithms for contact at material interfaces are the two challenging algorithm research problems that are addressed. Two prototype parallel contact algorithms have been developed for treating the cases of local and arbitrary contact. Demonstration problems for local contact are crashworthiness simulations with 222 locally defined contact surfaces and a vehicle/barrier collision modeled with arbitrary contact. A simulation of crash tests conducted for a vehicle impacting a U-channel small sign post embedded in soil has been run on both the serial and parallel versions of DYNA3D. A significant reduction in computational time has been observed when running these problems on the parallel version. However, to achieve maximum efficiency, complex problems must be appropriately partitioned, especially when contact dominates the computation.

  4. DRACO development for 3D simulations

    NASA Astrophysics Data System (ADS)

    Fatenejad, Milad; Moses, Gregory

    2006-10-01

    The DRACO (r-z) lagrangian radiation-hydrodynamics laser fusion simulation code is being extended to model 3D hydrodynamics in (x-y-z) coordinates with hexahedral cells on a structured grid. The equation of motion is solved with a lagrangian update with optional rezoning. The fluid equations are solved using an explicit scheme based on (Schulz, 1964) while the SALE-3D algorithm (Amsden, 1981) is used as a template for computing cell volumes and other quantities. A second order rezoner has been added which uses linear interpolation of the underlying continuous functions to preserve accuracy (Van Leer, 1976). Artificial restoring force terms and smoothing algorithms are used to avoid grid distortion in high aspect ratio cells. These include alternate node couplers along with a rotational restoring force based on the Tensor Code (Maenchen, 1964). Electron and ion thermal conduction is modeled using an extension of Kershaw's method (Kershaw, 1981) to 3D geometry. Test problem simulations will be presented to demonstrate the applicability of this new version of DRACO to the study of fluid instabilities in three dimensions.

  5. Analytical augmentation of 3D simulation environments

    NASA Astrophysics Data System (ADS)

    Loughran, Julia J.; Stahl, Marchelle M.

    1998-05-01

    This paper describes an approach for augmenting three- dimensional (3D) virtual environments (VEs) with analytic information and multimedia annotations to enhance training and education applications. Analytic or symbolic information in VEs is presented as bar charts, text, graphical overlays, or with the use of color. Analytic results can be computed and displayed in the VE at run-time or, more likely, while replaying a simulation. These annotations would typically include computations of pre-defined Measures of Performance (MOPs) or Measures of Effectiveness (MOEs) associated with the training or educational goals of the simulation. Multimedia annotations are inserted into the VE by the user and may include: a drawing or whiteboarding capability, enabling participants to insert written text and/or graphics into the two-dimensional (2D) or 3D world; audio comments, and/or video recordings. These annotations can clarify a point, capture teacher feedback, or elaborate on the student's perspective or understanding of the experience. The annotations are captured in the VE either synchronously or asynchronously from the users (students and instructors), during simulation execution or afterward during a replay. When replaying or reviewing the simulation, the embedded annotations can be reviewed by a single user or by multiple users through the use of collaboration technologies. By augmenting 3D virtual environments with analytic and multimedia annotations, the education and training experience may be enhanced. The annotations can offer more effective feedback, enhance understanding, and increase participation. They may also support distance learning by promoting student/teacher interaction without co-location.

  6. 3-D simulations of multiple beam klystrons

    SciTech Connect

    Smithe, David N.; Bettenhausen, Mike; Ludeking, Larry; Caryotakis, G.; Sprehn, Daryl; Scheitrum, Glenn

    1999-05-07

    The MAGIC3D simulation code is being used to assess the multi-dimensional physics issues relating to the design and operation of multiple beam klystrons. Investigations, to date, include a detailed study of the mode structure of the cavities in the 19-beam hexagonally packed geometry and a study of the velocity spread caused by the cavity mode's field profile. Some attempts to minimize this effect are investigated. Additional simulations have provided quantification of the beam loading Q in a dual input cavity, and optimization of a dual output cavity. An important goal of the simulations is an accurate picture of beam transport along the length of the MBK. We have quantified the magnitude and spatial variation of the beam-line space charge interactions within a cavity gap. Present simulations have demonstrated the transport of the beam through three cavities (the present limits of our simulation size) without difficulty; additional length simulations are expected. We have also examined unbalanced beam-line scenarios, e.g., one beam-line suppressed, and find little disturbance to the transport in individual cavity tests, with results for multiple cavity transport expected.

  7. 3D Hydrodynamic Simulations of Relativistic Jets

    NASA Astrophysics Data System (ADS)

    Hughes, P. A.; Miller, M. A.; Duncan, G. C.; Swift, C. M.

    1998-12-01

    We present the results of validation runs and the first extragalactic jet simulations performed with a 3D relativistic numerical hydrodynamic code employing a solver of the RHLLE type and using adaptive mesh refinement (AMR; Duncan & Hughes, 1994, Ap. J., 436, L119). Test problems include the shock tube, blast wave and spherical shock reflection (implosion). Trials with the code show that as a consequence of AMR it is viable to perform exploratory runs on workstation class machines (with no more than 128Mb of memory) prior to production runs. In the former case we achieve a resolution not much less than that normally regarded as the minimum needed to capture the essential physics of a problem, which means that such runs can provide valuable guidance allowing the optimum use of supercomputer resources. We present initial results from a program to explore the 3D stability properties of flows previously studied using a 2D axisymmetric code, and our first attempt to explore the structure and morphology of a relativistic jet encountering an ambient density gradient that mimics an ambient inhomogeneity or cloud.

  8. Simulating coronal condensation dynamics in 3D

    NASA Astrophysics Data System (ADS)

    Moschou, S. P.; Keppens, R.; Xia, C.; Fang, X.

    2015-12-01

    We present numerical simulations in 3D settings where coronal rain phenomena take place in a magnetic configuration of a quadrupolar arcade system. Our simulation is a magnetohydrodynamic simulation including anisotropic thermal conduction, optically thin radiative losses, and parametrised heating as main thermodynamical features to construct a realistic arcade configuration from chromospheric to coronal heights. The plasma evaporation from chromospheric and transition region heights eventually causes localised runaway condensation events and we witness the formation of plasma blobs due to thermal instability, that evolve dynamically in the heated arcade part and move gradually downwards due to interchange type dynamics. Unlike earlier 2.5D simulations, in this case there is no large scale prominence formation observed, but a continuous coronal rain develops which shows clear indications of Rayleigh-Taylor or interchange instability, that causes the denser plasma located above the transition region to fall down, as the system moves towards a more stable state. Linear stability analysis is used in the non-linear regime for gaining insight and giving a prediction of the system's evolution. After the plasma blobs descend through interchange, they follow the magnetic field topology more closely in the lower coronal regions, where they are guided by the magnetic dips.

  9. Elastically deformable 3D organs for haptic surgical simulation.

    PubMed

    Webster, Roger; Haluck, Randy; Ravenscroft, Rob; Mohler, Betty; Crouthamel, Eric; Frack, Tyson; Terlecki, Steve; Sheaffer, Jeremy

    2002-01-01

    This paper describes a technique for incorporating real-time elastically deformable 3D organs in haptic surgical simulators. Our system is a physically based particle model utilizing a mass-springs-damper connectivity with an implicit predictor to speed up calculations during each time step. The solution involves repeated application of Newton's 2ndd Law of motion: F = ma using an implicit solver for numerically solving the differential equations. PMID:15458154

  10. Flare Particle Escape in 3D Solar Eruptive Events

    NASA Astrophysics Data System (ADS)

    Antiochos, Spiro K.; Masson, Sophie; DeVore, C. R.

    2015-04-01

    Among the most important, but least understood forms of space weather are the so-called Impulsive Solar Energetic Particle (SEP) events, which can be especially hazardous to deep-space astronauts. These energetic particles are generally believed to be produced by the flare reconnection that is the primary driver of solar eruptive events (SEE). A key point is that in the standard model of SEEs, the particles should remain trapped in the coronal flare loops and in the ejected plasmoid, the CME. However, flare-accelerated particles frequently reach the Earth long before the CME does. In previous 2.5D calculations we showed how the external reconnection that is an essential element of the breakout model for CME initiation could lead to the escape of flare-accelerated particles. The problem, however, is that in 2.5D this reconnection also tends to destroy the plasmoid, which disagrees with the observation that SEP events are often associated with well-defined plasmoids at 1 AU known as “magnetic clouds”. Consequently, we have extended our model to a fully 3D topology that includes a multi-polar coronal field suitable for a breakout SEE near a coronal hole region. We performed high-resolution 3D MHD numerical simulations with the Adaptively Refined MHD Solver (ARMS). Our results demonstrate that the model allows for the effective escape of energetic particles from deep within an ejecting well-defined plasmoid. We show how the complex interactions between the flare and breakout reconnection reproduce all the main observational features of SEEs and SEPs. We discuss the implications of our calculations for the upcoming Solar Orbiter and Solar Probe Plus missions, which will measure SEEs and SEPs near the Sun, thereby, mitigating propagation effects.This research was supported, in part, by the NASA SR&T and TR&T Programs.

  11. Resolution improvement by 3D particle averaging in localization microscopy

    PubMed Central

    Broeken, Jordi; Johnson, Hannah; Lidke, Diane S.; Liu, Sheng; Nieuwenhuizen, Robert P.J.; Stallinga, Sjoerd; Lidke, Keith A.; Rieger, Bernd

    2015-01-01

    Inspired by recent developments in localization microscopy that applied averaging of identical particles in 2D for increasing the resolution even further, we discuss considerations for alignment (registration) methods for particles in general and for 3D in particular. We detail that traditional techniques for particle registration from cryo electron microscopy based on cross-correlation are not suitable, as the underlying image formation process is fundamentally different. We argue that only localizations, i.e. a set of coordinates with associated uncertainties, are recorded and not a continuous intensity distribution. We present a method that owes to this fact and that is inspired by the field of statistical pattern recognition. In particular we suggest to use an adapted version of the Bhattacharyya distance as a merit function for registration. We evaluate the method in simulations and demonstrate it on three-dimensional super-resolution data of Alexa 647 labelled to the Nup133 protein in the nuclear pore complex of Hela cells. From the simulations we find suggestions that for successful registration the localization uncertainty must be smaller than the distance between labeling sites on a particle. These suggestions are supported by theoretical considerations concerning the attainable resolution in localization microscopy and its scaling behavior as a function of labeling density and localization precision. PMID:25866640

  12. 3D flare particle model for ShipIR/NTCS

    NASA Astrophysics Data System (ADS)

    Ramaswamy, Srinivasan; Vaitekunas, David A.

    2016-05-01

    A key component in any soft-kill response to an incoming guided missile is the flare /chaff decoy used to distract or seduce the seeker homing system away from the naval platform. This paper describes a new 3D flare particle model in the naval threat countermeasure simulator (NTCS) of the NATO-standard ship signature model (ShipIR), which provides independent control over the size and radial distribution of its signature. The 3D particles of each flare sub-munition are modelled stochastically and rendered using OpenGL z-buffering, 2D projection, and alpha-blending to produce a unique and time varying signature. A sensitivity analysis on each input parameter provides the data and methods needed to synthesize a model from an IR measurement of a decoy. The new model also eliminated artifacts and deficiencies in our previous model which prevented reliable tracks from the adaptive track gate algorithm already presented by Ramaswamy and Vaitekunas (2015). A sequence of scenarios are used to test and demonstrate the new flare model during a missile engagement.

  13. 3D Simulations of the Beehive Proplyd

    NASA Astrophysics Data System (ADS)

    Feitosa, J. A.; Vasconcelos, M. J.; Cerqueira, A. H.

    2014-10-01

    Some star formation regions, like the Orion nebula, have stars of different masses, from massive stars, responsible for strong ionizing winds and HII regions, to low-mass stars, which spend a long time in the protostellar phase, and are frequently associated with protostellar disks and jets. Massive O or B stars emit a great deal of UV radiation, able to dissociate the hydrogen molecule (FUV radiation, energies between 6-13 eV), to ionize the atomic hydrogen (EUV radiation, energies greater than 13.6 eV) and heat the gas. Around these stars, a large and hot (10^{4}K) region is formed, known as HII region. T-Tauri stars inside HII regions produce a type of young stellar object, a proplyd, described with accuracy in O'Dell et al. (1993). Proplyds exhibit a cometary shape from which we can distinguish a central low-mass star with an accretion disk, an ionization front, a photodissociation region and, sometimes, an external bow shock and a protostellar jet. Its morphological characteristics depends on the distance between the low-mass star and the source of the ionizing radiation. The Beehive, a giant proplyd in Orion Nebula, has attracted attention due to its exotic system of rings coaxial to the HH540 jet's axis. Bally et al. (2005) suggested that the rings are perturbations due to the crossing of the ionization front by the jet. In this work, we test this hypothesis making 3D hydrodynamic numerical simulations over an adaptive grid, using the Yguazú-A code (Raga et al., 2000), properly adapted for the Beehive conditions. Our results show that the jet causes a perturbation in the ionization front of the proplyd, but is necessary to adjust carefully some parameters of the jet like its velocity and ejection frequency in order to have the results matching the observations.

  14. Optimizing Stellarators for Energetic Particle Confinement using BEAMS3D

    NASA Astrophysics Data System (ADS)

    Bolgert, Peter; Drevlak, Michael; Lazerson, Sam; Gates, David; White, Roscoe

    2015-11-01

    Energetic particle (EP) loss has been called the ``Achilles heel of stellarators,'' (Helander, Rep. Prog. Phys. 77 087001 (2014)) and there is a great need for magnetic configurations with improved EP confinement. In this study we utilize a newly developed capability of the stellarator optimization code STELLOPT: the ability to optimize EP confinement via an interface with guiding center code BEAMS3D (McMillan et al., Plasma Phys. Control. Fusion 56, 095019 (2014)). Using this new tool, optimizations of the W7-X experiment and ARIES-CS reactor are performed where the EP loss fraction is one of many target functions to be minimized. In W7-X, we simulate the experimental NBI system using realistic beam geometry and beam deposition physics. The goal is to find configurations with improved neutral beam deposition and energetic particle confinement. These calculations are compared to previous studies of W7-X NBI deposition. In ARIES-CS, we launch 3.5 MeV alpha particles from a near-axis flux surface using a uniform grid in toroidal and poloidal angle. As these particles are born from D-T reactions, we consider an isotropic distribution in velocity space. This research is supported by DoE Contract Number DE-AC02-09CH11466.

  15. Optofluidic fabrication for 3D-shaped particles.

    PubMed

    Paulsen, Kevin S; Di Carlo, Dino; Chung, Aram J

    2015-01-01

    Complex three-dimensional (3D)-shaped particles could play unique roles in biotechnology, structural mechanics and self-assembly. Current methods of fabricating 3D-shaped particles such as 3D printing, injection moulding or photolithography are limited because of low-resolution, low-throughput or complicated/expensive procedures. Here, we present a novel method called optofluidic fabrication for the generation of complex 3D-shaped polymer particles based on two coupled processes: inertial flow shaping and ultraviolet (UV) light polymerization. Pillars within fluidic platforms are used to deterministically deform photosensitive precursor fluid streams. The channels are then illuminated with patterned UV light to polymerize the photosensitive fluid, creating particles with multi-scale 3D geometries. The fundamental advantages of optofluidic fabrication include high-resolution, multi-scalability, dynamic tunability, simple operation and great potential for bulk fabrication with full automation. Through different combinations of pillar configurations, flow rates and UV light patterns, an infinite set of 3D-shaped particles is available, and a variety are demonstrated. PMID:25904062

  16. Optofluidic fabrication for 3D-shaped particles

    NASA Astrophysics Data System (ADS)

    Paulsen, Kevin S.; di Carlo, Dino; Chung, Aram J.

    2015-04-01

    Complex three-dimensional (3D)-shaped particles could play unique roles in biotechnology, structural mechanics and self-assembly. Current methods of fabricating 3D-shaped particles such as 3D printing, injection moulding or photolithography are limited because of low-resolution, low-throughput or complicated/expensive procedures. Here, we present a novel method called optofluidic fabrication for the generation of complex 3D-shaped polymer particles based on two coupled processes: inertial flow shaping and ultraviolet (UV) light polymerization. Pillars within fluidic platforms are used to deterministically deform photosensitive precursor fluid streams. The channels are then illuminated with patterned UV light to polymerize the photosensitive fluid, creating particles with multi-scale 3D geometries. The fundamental advantages of optofluidic fabrication include high-resolution, multi-scalability, dynamic tunability, simple operation and great potential for bulk fabrication with full automation. Through different combinations of pillar configurations, flow rates and UV light patterns, an infinite set of 3D-shaped particles is available, and a variety are demonstrated.

  17. 3D Numerical simulations of oblique subduction

    NASA Astrophysics Data System (ADS)

    Malatesta, C.; Gerya, T.; Scambelluri, M.; Crispini, L.; Federico, L.; Capponi, G.

    2012-04-01

    In the past 2D numerical studies (e.g. Gerya et al., 2002; Gorczyk et al., 2007; Malatesta et al., 2012) provided evidence that during intraoceanic subduction a serpentinite channel forms above the downgoing plate. This channel forms as a result of hydration of the mantle wedge by uprising slab-fluids. Rocks buried at high depths are finally exhumed within this buoyant low-viscosity medium. Convergence rate in these 2D models was described by a trench-normal component of velocity. Several present and past subduction zones worldwide are however driven by oblique convergence between the plates, where trench-normal motion of the subducting slab is coupled with trench-parallel displacement of the plates. Can the exhumation mechanism and the exhumation rates of high-pressure rocks be affected by the shear component of subduction? And how uprise of these rocks can vary along the plate margin? We tried to address these questions performing 3D numerical models that simulate an intraoceanic oblique subduction. The models are based on thermo-mechanical equations that are solved with finite differences method and marker-in-cell techniques combined with multigrid approach (Gerya, 2010). In most of the models a narrow oceanic basin (500 km-wide) surrounded by continental margins is depicted. The basin is floored by either layered or heterogeneous oceanic lithosphere with gabbro as discrete bodies in serpentinized peridotite and a basaltic layer on the top. A weak zone in the mantle is prescribed to control the location of subduction initiation and therefore the plate margins geometry. Finally, addition of a third dimension in the simulations allowed us to test the role of different plate margin geometries on oblique subduction dynamics. In particular in each model we modified the dip angle of the weak zone and its "lateral" geometry (e.g. continuous, segmented). We consider "continuous" weak zones either parallel or increasingly moving away from the continental margins

  18. Simulation of Current Generation in a 3-D Plasma Model

    NASA Astrophysics Data System (ADS)

    Tsung, F. S.; Dawson, J. M.

    1996-11-01

    In the advanced tokamak regime, transport phenomena can account for a signficant fraction of the toroidal current, possibly over that driven directly by the ohmic heating electric fields. Bootstrap theory accounts for contributions of the collisional modification of banana orbits on the toroidal currents. In our previous simulations in 21/2-D, currents were spontaneously generated in both the cylindrical and the toroidal geometries, contrary to neoclassical predictions. In these calculations, it was believed that the driving mechanism is the preferential loss of particles whose initial velocity is opposite to that of the plasma current. We are extending these simulations to three dimensions. A parallel 3-D electromagnetic PIC code running on the IBM SP2, with a localized field-solver has been developed to investigate the effects of perturbations parallel to the field lines, and direct comparisons has been made between the 21/2-D and 3-D simulations and we have found good agreements between the 2 1/2-D calculations and the 3-D results. We will present our new results at the meeting. Research partially supported by NSF and DOE.

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

    SciTech Connect

    Park, W.; Chang, Z.; Fredrickson, E.; Fu, G.Y.

    1996-12-31

    The M3D (Multi-level 3D) tokamak simulation project aims at the simulation of tokamak plasmas using a multi-level tokamak code package. Several current applications using MHD and Extended-MHD models are presented; high-{beta} disruption studies in reversed shear plasmas using the MHD level MH3D code, {omega}{sub *i} stabilization and nonlinear island saturation of TAE mode using the hybrid particle/MHD level MH3D-K code, and unstructured mesh MH3D{sup ++} code studies. In particular, three internal mode disruption mechanisms are identified from simulation results which agree which agree well with experimental data.

  20. Light Attenuation Method for 3D data acquisition (LAM3D) of bottom particle deposits

    NASA Astrophysics Data System (ADS)

    Er, Jenn Wei; Law, Adrian W. K.; Adams, E. Eric; Yang, Yang

    2015-11-01

    We have developed a novel experimental technique, Light Attenuation Method for 3D data acquisition (LAM3D), to acquire three-dimensional spatial characteristics and temporal development of bottom particle deposits. The new technique performs data acquisition with higher spatial and temporal resolution than existing approaches with laser and ultrasonic 3D profilers, and is therefore ideal for laboratory investigations with fast varying changes in the sediment bed, such as the developing deposition profile from sediment clouds commonly formed during dredging or land reclamation projects and the dynamic evolution in movable bed processes in rivers. The principle of the technique is based on the analysis of the light attenuation due to multiple light scattering through the particle deposits layer compared to the clear water column. With appropriate calibration, the particles size and distribution thickness can be quantified by the transmitted light spectrum. In the presentation, we will first show our measurement setup with a light panel for calibrated illumination and a system of DSLR cameras for the light capturing. Subsequently, we shall present the experimental results of fast evolving deposition profile of a barge-disposed sediment cloud upon its bottom impact on the sea bed.

  1. 3D simulations of an electrostatic quadrupole injector

    SciTech Connect

    Grote, D.P. |; Friedman, A.; Yu, S.

    1993-02-01

    Analysis of the dynamics of a space charge dominated beam in a lattice of electrostatic focusing structures requires a full three-dimensional conic that includes self-consistent space charge fields and the fields from the complex conductor shapes. The existing WARP3d code, a particle simulation code which has been developed for heavy-ion fusion (HIF) applications contains machinery for handling particles in three-dimensional fields. A successive overrelaxation field solver with subgrid-scale placement of boundaries for rounded surface and four-fold symmetry has been added to the code. The electrostatic quadrupole (ESQ) injector for the ILSE accelerator facility being planned at Lawrence Berkeley Laboratory is shown as an application. The issue of concern is possible emittance degradation because the focusing voltages are a significant fraction of the particles` energy and because there are significant nonlinear fields arising from the shapes of the quadrupole structures.

  2. Numerical simulation of 3D breaking waves

    NASA Astrophysics Data System (ADS)

    Fraunie, Philippe; Golay, Frederic

    2015-04-01

    Numerical methods dealing with two phase flows basically can be classified in two ways : the "interface tracking" methods when the two phases are resolved separately including boundary conditions fixed at the interface and the "interface capturing" methods when a single flow is considered with variable density. Physical and numerical properties of the two approaches are discussed, based on some numerical experiments performed concerning 3D breaking waves. Acknowledgements : This research was supported by the Modtercom program of Region PACA.

  3. Ultrafine particle emissions from desktop 3D printers

    NASA Astrophysics Data System (ADS)

    Stephens, Brent; Azimi, Parham; El Orch, Zeineb; Ramos, Tiffanie

    2013-11-01

    The development of low-cost desktop versions of three-dimensional (3D) printers has made these devices widely accessible for rapid prototyping and small-scale manufacturing in home and office settings. Many desktop 3D printers rely on heated thermoplastic extrusion and deposition, which is a process that has been shown to have significant aerosol emissions in industrial environments. However, we are not aware of any data on particle emissions from commercially available desktop 3D printers. Therefore, we report on measurements of size-resolved and total ultrafine particle (UFP) concentrations resulting from the operation of two types of commercially available desktop 3D printers inside a commercial office space. We also estimate size-resolved (11.5 nm-116 nm) and total UFP (<100 nm) emission rates and compare them to emission rates from other desktop devices and indoor activities known to emit fine and ultrafine particles. Estimates of emission rates of total UFPs were large, ranging from ˜2.0 × 1010 # min-1 for a 3D printer utilizing a polylactic acid (PLA) feedstock to ˜1.9 × 1011 # min-1 for the same type of 3D printer utilizing a higher temperature acrylonitrile butadiene styrene (ABS) thermoplastic feedstock. Because most of these devices are currently sold as standalone devices without any exhaust ventilation or filtration accessories, results herein suggest caution should be used when operating in inadequately ventilated or unfiltered indoor environments. Additionally, these results suggest that more controlled experiments should be conducted to more fundamentally evaluate particle emissions from a wider arrange of desktop 3D printers.

  4. Light shaping along 3D curves and particle manipulation

    NASA Astrophysics Data System (ADS)

    Rodrigo, José A.; Alieva, Tatiana

    2015-03-01

    We present a non-iterative holographic technique for efficient and versatile laser beam shaping along arbitrary 3D curves. Light beams with intensity shaped for several 3D curves: Tilted ring, Viviani's curve, Archimedean spiral, and trefoil-knotted curve have been experimentally generated and applied for optical trapping of micrometer-sized dielectric particles. The high intensity gradients and independent phase control prescribed along the curve make this kind of laser trap attractive for multiple particle manipulation and allow for forward and backward motion to the light source. Indeed, different configurations of tractor beam traps are experimentally demonstrated. This technique can also be applied for laser micro-machining.

  5. Iterative Reconstruction of Volumetric Particle Distribution for 3D Velocimetry

    NASA Astrophysics Data System (ADS)

    Wieneke, Bernhard; Neal, Douglas

    2011-11-01

    A number of different volumetric flow measurement techniques exist for following the motion of illuminated particles. For experiments that have lower seeding densities, 3D-PTV uses recorded images from typically 3-4 cameras and then tracks the individual particles in space and time. This technique is effective in flows that have lower seeding densities. For flows that have a higher seeding density, tomographic PIV uses a tomographic reconstruction algorithm (e.g. MART) to reconstruct voxel intensities of the recorded volume followed by the cross-correlation of subvolumes to provide the instantaneous 3D vector fields on a regular grid. A new hybrid algorithm is presented which iteratively reconstructs the 3D-particle distribution directly using particles with certain imaging properties instead of voxels as base functions. It is shown with synthetic data that this method is capable of reconstructing densely seeded flows up to 0.05 particles per pixel (ppp) with the same or higher accuracy than 3D-PTV and tomographic PIV. Finally, this new method is validated using experimental data on a turbulent jet.

  6. 3-D numerical simulations of volcanic ash transport and deposition

    NASA Astrophysics Data System (ADS)

    Suzuki, Y. J.; Koyaguchi, T.

    2012-12-01

    During an explosive volcanic eruption, volcanic gas and pyroclasts are ejected from the volcanic vent. The pyroclasts are carried up within a convective plume, advected by the surrounding wind field, and sediment on the ground depending on their terminal velocity. The fine ash are expected to have atmospheric residence, whereas the coarser particles form fall deposits. Accurate modeling of particle transport and deposition is of critical importance from the viewpoint of disaster prevention. Previously, some particle-tracking models (e.g., PUFF) and advection-diffusion models (e.g., TEPHRA2 and FALL3D) tried to forecast particle concentration in the atmosphere and particle loading at ground level. However, these models assumed source conditions (the grain-size distribution, plume height, and mass release location) based on the simple 1-D model of convective plume. In this study, we aim to develop a new 3-D model which reproduces both of the dynamics of convective plume and the ash transport. The model is designed to describe the injection of eruption cloud and marker particles from a circular vent above a flat surface into the stratified atmosphere. Because the advection is the predominant mechanism of particle transport near the volcano, the diffusive process is not taken into account in this model. The distribution of wind velocity is given as an initial condition. The model of the eruption cloud dynamics is based on the 3-D time-dependent model of Suzuki et al. (2005). We apply a pseudo-gas model to calculate the eruption cloud dynamics: the effect of particle separation on the cloud dynamics is not considered. In order to reproduce the drastic change of eruption cloud density, we change the effective gas constant and heat capacity of the mixture in the equation of state for ideal gases with the mixing ratio between the ejected material and entrained air. In order to calculate the location and movement of ash particles, the present model employs Lagrangian marker

  7. Vertical Flow Lithography for Fabrication of 3D Anisotropic Particles.

    PubMed

    Habasaki, Shohei; Lee, Won Chul; Yoshida, Shotaro; Takeuchi, Shoji

    2015-12-22

    A microfluidics-based method for the 3D fabrication of anisotropic particles is reported. The method uses a vertical microchannel where tunable light patterns solidify photocurable resins for stacking multiple layers of the resins, thus enabling an application of stereolithography concepts to conventional flow lithography. Multilayered, tapered, and angular compartmental microparticles are demonstrated. PMID:26551590

  8. 3D visualization of port simulation.

    SciTech Connect

    Horsthemke, W. H.; Macal, C. M.; Nevins, M. R.

    1999-06-14

    Affordable and realistic three dimensional visualization technology can be applied to large scale constructive simulations such as the port simulation model, PORTSIM. These visualization tools enhance the experienced planner's ability to form mental models of how seaport operations will unfold when the simulation model is implemented and executed. They also offer unique opportunities to train new planners not only in the use of the simulation model but on the layout and design of seaports. Simulation visualization capabilities are enhanced by borrowing from work on interface design, camera control, and data presentation. Using selective fidelity, the designers of these visualization systems can reduce their time and efforts by concentrating on those features which yield the most value for their simulation. Offering the user various observational tools allows the freedom to simply watch or engage in the simulation without getting lost. Identifying the underlying infrastructure or cargo items with labels can provide useful information at the risk of some visual clutter. The PortVis visualization expands the PORTSIM user base which can benefit from the results provided by this capability, especially in strategic planning, mission rehearsal, and training. Strategic planners will immediately reap the benefits of seeing the impact of increased throughput visually without keeping track of statistical data. Mission rehearsal and training users will have an effective training tool to supplement their operational training exercises which are limited in number because of their high costs. Having another effective training modality in this visualization system allows more training to take place and more personnel to gain an understanding of seaport operations. This simulation and visualization training can be accomplished at lower cost than would be possible for the operational training exercises alone. The application of PORTSIM and PortVis will lead to more efficient

  9. 3D implicit PIC simulations of solar wind - moon interactions

    NASA Astrophysics Data System (ADS)

    Deca, J.; Markidis, S.; Divin, A.; Lapenta, G.; Vapirev, A.

    2012-04-01

    We present three-dimensional Particle-in-Cell simulations of an unmagnetized insulating Moon-sized body immersed in the solar wind. The simulations are performed using the implicit electromagnetic Particle-in-Cell code iPIC3D [Markidis, 2009]. Multiscale kinetic physics is resolved for all plasma components (heavy ions, protons and electrons) in the code, recently updated with a set of open boundary conditions designed for solar wind - body interaction studies. Particles are injected at the inflow side of the computational domain and absorbed at all others. A bow shock is not formed upstream of the body, but the obstacle generates faint dispersive waves propagating parallel to the magnetic field lines, in agreement with numerical simulations done in MHD approach. Polarization electric field is generated in the wake. In addition, plasma flows filling the wake tend to excite streaming instabilities, which lead to bipolar signatures in the parallel electric field. Our future work includes updating the physical model to include photoionization and re-emission at the object's surface.

  10. Simulation of 3D infrared scenes using random fields model

    NASA Astrophysics Data System (ADS)

    Shao, Xiaopeng; Zhang, Jianqi

    2001-09-01

    Analysis and simulation of smart munitions requires imagery for the munition's sensor to view. The traditional infrared background simulations are always limited in the plane scene studies. A new method is described to synthesize the images in 3D view and with various terrains texture. We develop the random fields model and temperature fields to simulate 3D infrared scenes. Generalized long-correlation (GLC) model, one of random field models, will generate both the 3D terrains skeleton data and the terrains texture in this work. To build the terrain mesh with the random fields, digital elevation models (DEM) are introduced in the paper. And texture mapping technology will perform the task of pasting the texture in the concavo-convex surfaces of the 3D scene. The simulation using random fields model is a very available method to produce 3D infrared scene with great randomicity and reality.

  11. 3-D PARTICLE TRANSPORT WITHIN THE HUMAN UPPER RESPIRATORY TRACT

    EPA Science Inventory

    In this study trajectories of inhaled particulate matter (PM) were simulated within a three-dimensional (3-D) computer model of the human upper respiratory tract (URT). The airways were described by computer-reconstructed images of a silicone rubber cast of the human head, throat...

  12. 3D MHD Simulations of Spheromak Compression

    NASA Astrophysics Data System (ADS)

    Stuber, James E.; Woodruff, Simon; O'Bryan, John; Romero-Talamas, Carlos A.; Darpa Spheromak Team

    2015-11-01

    The adiabatic compression of compact tori could lead to a compact and hence low cost fusion energy system. The critical scientific issues in spheromak compression relate both to confinement properties and to the stability of the configuration undergoing compression. We present results from the NIMROD code modified with the addition of magnetic field coils that allow us to examine the role of rotation on the stability and confinement of the spheromak (extending prior work for the FRC). We present results from a scan in initial rotation, from 0 to 100km/s. We show that strong rotational shear (10km/s over 1cm) occurs. We compare the simulation results with analytic scaling relations for adiabatic compression. Work performed under DARPA grant N66001-14-1-4044.

  13. Rockfall assessment of natural hazards by 3D-simulation

    NASA Astrophysics Data System (ADS)

    Tobler, D.; Graf, K.; Krummenacher, B.

    2009-04-01

    Keywords: 3D rockfall modelling, protection measures, case study, ROFMOD 4.1 model In mountain areas many residential areas as well as important lifelines are exposed to potential rockfall events. The assessment of the hazard potential in a complex morphology with changing surface parameters and inhomogeneous forestation, call for an in-depth analysis of the physical processes involved. Therefore the different problems concerning the assessment of rockfall hazards have to be solved adequately. Moreover, due to the actual discussion about cost efficiency the planning of protective measures calls for (more) detailed information about the intensity and probability of expected rockfall incidents along a typical trajectory or for a given area. Today simulation models like GEOTEST's - ROFMOD 4.1 (2D and 3D) are important tools for hazard assessment. This simulation model produces profiles and maps with intensity and range and allows to calculate decisive parameters for the dimensioning optimized positioning of protection measures. Within the impact analysis the relevant scenarios according to their intensities are determined using a 3D rockfall model. It calculates the trajectories of rock-bodies with realistic forms based on physical laws from defined starting points within the mapped rock area as well as steep slopes with loose blocks. The movement is of any chosen particle with defined axis length ratio, rounding and specific weight is strongly influenced by the spatial resolution of the digital terrain model (DTM as well as the interaction with the underground and possible tree-impacts). Therefore, the following parameters are implemented in the rockfall model: geometry of blocks, underground plasticity (value of attenuation), surface roughness as well as forest quality (amount of trees per area (ha) and trunk diameter). A very important aspect is the geometry of the blocks. For realistic results only a simulation with real block shapes will be sucessful. For these

  14. 3D radiative transfer in colliding wind binaries: Application of the SimpleX algorithm to 3D SPH simulations

    NASA Astrophysics Data System (ADS)

    Madura, Thomas; Clementel, Nicola; Kruip, Chael; Icke, Vincent; Gull, Theodore

    2014-09-01

    We present the first results of full 3D radiative transfer simulations of the colliding stellar winds in a massive binary system. We accomplish this by applying the SIMPLEX algorithm for 3D radiative transfer on an unstructured Delaunay grid to recent 3D smoothed particle hydrodynamics (SPH) simulations of the colliding winds in the binary system η Carinae. We use SIMPLEX to obtain detailed ionization fractions of hydrogen and helium, in 3D, at the resolution of the original SPH simulations. We show how the SIMPLEX simulations can be used to generate synthetic spectral data cubes for comparison to data obtained with the Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph as part of a multi-cycle program to map changes in η Car's extended interacting wind structures across one binary cycle. Comparison of the HST observations to the SIMPLEX models can help lead to more accurate constraints on the orbital, stellar, and wind parameters of the η Car system, such as the primary's mass-loss rate and the companion's temperature and luminosity. While we initially focus specifically on the η Car binary, the numerical methods employed can be applied to numerous other colliding wind (WR140, WR137, WR19) and dusty 'pinwheel' (WR104, WR98a) binary systems. One of the biggest remaining mysteries is how dust can form and survive in such systems that contain a hot, luminous O star. Coupled with 3D hydrodynamical simulations, SIMPLEX simulations have the potential to help determine the regions where dust can form and survive in these unique objects.

  15. Design of 3D simulation engine for oilfield safety training

    NASA Astrophysics Data System (ADS)

    Li, Hua-Ming; Kang, Bao-Sheng

    2015-03-01

    Aiming at the demand for rapid custom development of 3D simulation system for oilfield safety training, this paper designs and implements a 3D simulation engine based on script-driven method, multi-layer structure, pre-defined entity objects and high-level tools such as scene editor, script editor, program loader. A scripting language been defined to control the system's progress, events and operating results. Training teacher can use this engine to edit 3D virtual scenes, set the properties of entity objects, define the logic script of task, and produce a 3D simulation training system without any skills of programming. Through expanding entity class, this engine can be quickly applied to other virtual training areas.

  16. An Evaluative Review of Simulated Dynamic Smart 3d Objects

    NASA Astrophysics Data System (ADS)

    Romeijn, H.; Sheth, F.; Pettit, C. J.

    2012-07-01

    Three-dimensional (3D) modelling of plants can be an asset for creating agricultural based visualisation products. The continuum of 3D plants models ranges from static to dynamic objects, also known as smart 3D objects. There is an increasing requirement for smarter simulated 3D objects that are attributed mathematically and/or from biological inputs. A systematic approach to plant simulation offers significant advantages to applications in agricultural research, particularly in simulating plant behaviour and the influences of external environmental factors. This approach of 3D plant object visualisation is primarily evident from the visualisation of plants using photographed billboarded images, to more advanced procedural models that come closer to simulating realistic virtual plants. However, few programs model physical reactions of plants to external factors and even fewer are able to grow plants based on mathematical and/or biological parameters. In this paper, we undertake an evaluation of plant-based object simulation programs currently available, with a focus upon the components and techniques involved in producing these objects. Through an analytical review process we consider the strengths and weaknesses of several program packages, the features and use of these programs and the possible opportunities in deploying these for creating smart 3D plant-based objects to support agricultural research and natural resource management. In creating smart 3D objects the model needs to be informed by both plant physiology and phenology. Expert knowledge will frame the parameters and procedures that will attribute the object and allow the simulation of dynamic virtual plants. Ultimately, biologically smart 3D virtual plants that react to changes within an environment could be an effective medium to visually represent landscapes and communicate land management scenarios and practices to planners and decision-makers.

  17. 3D Ultrasonic Wave Simulations for Structural Health Monitoring

    NASA Technical Reports Server (NTRS)

    Campbell, Leckey Cara A/; Miler, Corey A.; Hinders, Mark K.

    2011-01-01

    Structural health monitoring (SHM) for the detection of damage in aerospace materials is an important area of research at NASA. Ultrasonic guided Lamb waves are a promising SHM damage detection technique since the waves can propagate long distances. For complicated flaw geometries experimental signals can be difficult to interpret. High performance computing can now handle full 3-dimensional (3D) simulations of elastic wave propagation in materials. We have developed and implemented parallel 3D elastodynamic finite integration technique (3D EFIT) code to investigate ultrasound scattering from flaws in materials. EFIT results have been compared to experimental data and the simulations provide unique insight into details of the wave behavior. This type of insight is useful for developing optimized experimental SHM techniques. 3D EFIT can also be expanded to model wave propagation and scattering in anisotropic composite materials.

  18. Digital In-Line Holography System for 3D-3C Particle Tracking Velocimetry

    NASA Astrophysics Data System (ADS)

    Malek, Mokrane; Lebrun, Denis; Allano, Daniel

    Digital in-line holography is a suitable method for measuring three dimensional (3D) velocity fields. Such a system records directly on a charge-coupled device (CCD) camera a couple of diffraction patterns produced by small particles illuminated by a modulated laser diode. The numerical reconstruction is based on the wavelet transformation method. A 3D particle field is reconstructed by computing the wavelet components for different scale parameters. The scale parameter is directly related to the axial distance between a given particle and the CCD camera. The particle images are identified and localized by analyzing the maximum of the wavelet transform modulus (WTMM) and the equivalent diameter of the particle image (Deq). Afterwards, a 3D point-matching (PM) algorithm is applied to the pair of sets containing the 3D particle locations. In the PM algorithm, the displacement of the particles is modeled by an affine transformation. This affine transformation is based on the use of the dual number quaternions. Afterwards, the velocity-field extraction is performed. This system is tested with simulated particle field displacements and the feasibility is checked with an experimental displacement.

  19. Computationally efficient storage of 3D particle intensity and position data for use in 3D PIV and 3D PTV

    NASA Astrophysics Data System (ADS)

    Atkinson, C.; Buchmann, N. A.; Soria, J.

    2013-11-01

    Three-dimensional (3D) volumetric velocity measurement techniques, such as tomographic or holographic particle image velocimetry (PIV), rely upon the computationally intensive formation, storage and localized interrogation of multiple 3D particle intensity fields. Calculation of a single velocity field typically requires the extraction of particle intensities into tens of thousands of 3D sub-volumes or discrete particle clusters, the processing of which can significantly affect the performance of 3D cross-correlation based PIV and 3D particle tracking velocimetry (PTV). In this paper, a series of popular and customized volumetric data formats are presented and investigated using synthetic particle volumes and experimental data arising from tomographic PIV measurements of a turbulent boundary layer. Results show that the use of a sub-grid ordered non-zero intensity format with a sub-grid size of 16 × 16 × 16 points provides the best performance for cross-correlation based PIV analysis, while a particle clustered non-zero intensity format provides the best format for PTV applications. In practical tomographic PIV measurements the sub-grid ordered non-zero intensity format offered a 29% improvement in reconstruction times, while providing a 93% reduction in volume data requirements and a 28% overall improvement in cross-correlation based velocity analysis and validation times.

  20. First 3-D simulations of meteor plasma dynamics and turbulence

    NASA Astrophysics Data System (ADS)

    Oppenheim, Meers M.; Dimant, Yakov S.

    2015-02-01

    Millions of small but detectable meteors hit the Earth's atmosphere every second, creating trails of hot plasma that turbulently diffuse into the background atmosphere. For over 60 years, radars have detected meteor plasmas and used these signals to infer characteristics of the meteoroid population and upper atmosphere, but, despite the importance of meteor radar measurements, the complex processes by which these plasmas evolve have never been thoroughly explained or modeled. In this paper, we present the first fully 3-D simulations of meteor evolution, showing meteor plasmas developing instabilities, becoming turbulent, and inhomogeneously diffusing into the background ionosphere. These instabilities explain the characteristics and strength of many radar observations, in particular the high-resolution nonspecular echoes made by large radars. The simulations reveal how meteors create strong electric fields that dig out deep plasma channels along the Earth's magnetic fields. They also allow researchers to explore the impacts of the intense winds and wind shears, commonly found at these altitudes, on meteor plasma evolution. This study will allow the development of more sophisticated models of meteor radar signals, enabling the extraction of detailed information about the properties of meteoroid particles and the atmosphere.

  1. Simulations of Coulomb systems with slab geometry using an efficient 3D Ewald summation method

    NASA Astrophysics Data System (ADS)

    dos Santos, Alexandre P.; Girotto, Matheus; Levin, Yan

    2016-04-01

    We present a new approach to efficiently simulate electrolytes confined between infinite charged walls using a 3d Ewald summation method. The optimal performance is achieved by separating the electrostatic potential produced by the charged walls from the electrostatic potential of electrolyte. The electric field produced by the 3d periodic images of the walls is constant inside the simulation cell, with the field produced by the transverse images of the charged plates canceling out. The non-neutral confined electrolyte in an external potential can be simulated using 3d Ewald summation with a suitable renormalization of the electrostatic energy, to remove a divergence, and a correction that accounts for the conditional convergence of the resulting lattice sum. The new algorithm is at least an order of magnitude more rapid than the usual simulation methods for the slab geometry and can be further sped up by adopting a particle-particle particle-mesh approach.

  2. 3d visualization of atomistic simulations on every desktop

    NASA Astrophysics Data System (ADS)

    Peled, Dan; Silverman, Amihai; Adler, Joan

    2013-08-01

    Once upon a time, after making simulations, one had to go to a visualization center with fancy SGI machines to run a GL visualization and make a movie. More recently, OpenGL and its mesa clone have let us create 3D on simple desktops (or laptops), whether or not a Z-buffer card is present. Today, 3D a la Avatar is a commodity technique, presented in cinemas and sold for home TV. However, only a few special research centers have systems large enough for entire classes to view 3D, or special immersive facilities like visualization CAVEs or walls, and not everyone finds 3D immersion easy to view. For maximum physics with minimum effort a 3D system must come to each researcher and student. So how do we create 3D visualization cheaply on every desktop for atomistic simulations? After several months of attempts to select commodity equipment for a whole room system, we selected an approach that goes back a long time, even predating GL. The old concept of anaglyphic stereo relies on two images, slightly displaced, and viewed through colored glasses, or two squares of cellophane from a regular screen/projector or poster. We have added this capability to our AViz atomistic visualization code in its new, 6.1 version, which is RedHat, CentOS and Ubuntu compatible. Examples using data from our own research and that of other groups will be given.

  3. 3-D laser radar simulation for autonomous spacecraft landing

    NASA Technical Reports Server (NTRS)

    Reiley, Michael F.; Carmer, Dwayne C.; Pont, W. F.

    1991-01-01

    A sophisticated 3D laser radar sensor simulation, developed and applied to the task of autonomous hazard detection and avoidance, is presented. This simulation includes a backward ray trace to sensor subpixels, incoherent subpixel integration, range dependent noise, sensor point spread function effects, digitization noise, and AM-CW modulation. Specific sensor parameters, spacecraft lander trajectory, and terrain type have been selected to generate simulated sensor data.

  4. Current loop coalescence studied by 3-D electromagnetic particle code

    NASA Technical Reports Server (NTRS)

    Nishikawa, Ken-Ichi; Sakai, Jun-Ichi; Koide, Shinji; Buneman, O.; Neubert, T.

    1993-01-01

    Solar flare plasma data from the Yohkoh satellite is analyzed. The interactions of current loops were observed in the active regions on the Sun. This observation pointed out the importance of the idea that the solar flare is generated by the coalescence of current loops. The three dimensional electromagnetic particle simulations are to help in understanding the global interaction between two current loops including the evolution of the twist of loops due to instabilities. Associated rapid dynamics of current loop coalescence such as reconnection, shock waves and associated kinetic processes such as energy transfer, acceleration of particles, and electromagnetic emissions are to be studied by the code to complement analytical theories and magnetohydrodynamic simulations of the current loop coalescence. The simulation results show the strong interactions between two current loops, beam and whistler instabilities, and associated parallel and perpendicular particle heating.

  5. 3D Fiber Orientation Simulation for Plastic Injection Molding

    NASA Astrophysics Data System (ADS)

    Lin, Baojiu; Jin, Xiaoshi; Zheng, Rong; Costa, Franco S.; Fan, Zhiliang

    2004-06-01

    Glass fiber reinforced polymer is widely used in the products made using injection molding processing. The distribution of fiber orientation inside plastic parts has direct effects on quality of molded parts. Using computer simulation to predict fiber orientation distribution is one of most efficient ways to assist engineers to do warpage analysis and to find a good design solution to produce high quality plastic parts. Fiber orientation simulation software based on 2-1/2D (midplane /Dual domain mesh) techniques has been used in industry for a decade. However, the 2-1/2D technique is based on the planar Hele-Shaw approximation and it is not suitable when the geometry has complex three-dimensional features which cannot be well approximated by 2D shells. Recently, a full 3D simulation software for fiber orientation has been developed and integrated into Moldflow Plastics Insight 3D simulation software. The theory for this new 3D fiber orientation calculation module is described in this paper. Several examples are also presented to show the benefit in using 3D fiber orientation simulation.

  6. 3D mapping and simulation of Geneva Lake environmental data

    NASA Astrophysics Data System (ADS)

    Villard, Roch; Maignan, Michel; Kanevski, Mikhail; Rapin, Francois; Klein, Audrey

    2010-05-01

    The Geneva Lake is the biggest alpine and subalpine lake in central Europe. The depth of this lake is 309 meters and its total volume of water is 89 billions m3. It takes, on average, around twelve years so that waters of the lake are completely brewed. Furthermore the Geneva lake waters are rich in dissolved substances as carbonate, sulfate. The quantity of particles in suspension in the lake, which mainly arrived from the Rhône, is nowadays around height million of tones. The International Commission for the Leman Lake (CIPEL) works about the improvement of the quality of this lake since 1962. In the present study three dimensional environmental data (temperature, oxygen and nitrate) which cover the period from 1954 to 2008, for a total of 27'500 cases are investigated. We are interested to study the evolution of the temperature of the lake because there is an impact on the reproduction of fishes and also because the winter brewing of the water makes the re-oxygenation of deep-water. In order that biological balance is maintained in a lake, there must be enough oxygen in the water. Moreover, we work on nitrate distribution and evolution because contributions in fertilizers cause eutrophication of lake. The data are very numerous when we consider the time series, some of them with more than 300 occurrences, but there are between 2 and 15 data available for spatial cartography. The basic methodology used for the analysis, mapping and simulations of 3D patterns of environmental data is based on geostatistical predictions (family of kriging models) and conditional stochastic simulations. Spatial and temporal variability, 3D monitoring networks changing over time, make this study challenging. An important problem is also to make interpolation/simulations over a long period of time, like ten years. One way used to overcome this problem, consists in using a weighted average of ten variograms during this period. 3D mapping was carried out using environment data for

  7. 3D two-fluid simulations of turbulence in LAPD

    NASA Astrophysics Data System (ADS)

    Fisher, Dustin M.

    The Large Plasma Device (LAPD) is modeled using a modified version of the 3D Global Braginskii Solver code (GBS) for a nominal Helium plasma. The unbiased low-flow regime is explored in simulations where there is an intrinsic E x B rotation of the plasma. In the simulations this rotation is caused primarily by sheath effects with the Reynolds stress and J x B torque due to a cross-field Pederson conductivity having little effect. Explicit biasing simulations are also explored for the first time where the intrinsic rotation of the plasma is modified through boundary conditions that mimic the biasable limiter used in LAPD. Comparisons to experimental measurements in the unbiased case show strong qualitative agreement with the data, particularly the radial dependence of the density fluctuations, cross-correlation lengths, radial flux dependence outside of the cathode edge, and camera imagery. Kelvin Helmholtz (KH) turbulence at relatively large scales is the dominant driver of cross-field transport in these simulations with smaller-scale drift waves and sheath modes playing a secondary role. Plasma holes and blobs arising from KH vortices are consistent with the scale sizes and overall appearance of those in LAPD camera images. The addition of ion-neutral collisions in the unbiased simulations at previously theorized values reduces the radial particle flux due to a modest stabilizing contribution of the collisions on the KH-modes driving the turbulent transport. In the biased runs the ion-neutral collisions have a much smaller effect due to the modification of the potential from sheath terms. In biasing the plasma to increase the intrinsic rotation, simulations show the emergence of a nonlinearly saturated coherent mode of order m = 6. In addition, the plasma inside of the cathode edge becomes quiescent due to the strong influence of the wall bias in setting up the equilibrium plasma potential. Biasing in the direction opposite to the intrinsic flow reduces the

  8. Parallelization of Program to Optimize Simulated Trajectories (POST3D)

    NASA Technical Reports Server (NTRS)

    Hammond, Dana P.; Korte, John J. (Technical Monitor)

    2001-01-01

    This paper describes the parallelization of the Program to Optimize Simulated Trajectories (POST3D). POST3D uses a gradient-based optimization algorithm that reaches an optimum design point by moving from one design point to the next. The gradient calculations required to complete the optimization process, dominate the computational time and have been parallelized using a Single Program Multiple Data (SPMD) on a distributed memory NUMA (non-uniform memory access) architecture. The Origin2000 was used for the tests presented.

  9. 3D structural fluctuation of IgG1 antibody revealed by individual particle electron tomography

    SciTech Connect

    Zhang, Xing; Zhang, Lei; Tong, Huimin; Peng, Bo; Rames, Matthew J.; Zhang, Shengli; Ren, Gang

    2015-05-05

    Commonly used methods for determining protein structure, including X-ray crystallography and single-particle reconstruction, often provide a single and unique three-dimensional (3D) structure. However, in these methods, the protein dynamics and flexibility/fluctuation remain mostly unknown. Here, we utilized advances in electron tomography (ET) to study the antibody flexibility and fluctuation through structural determination of individual antibody particles rather than averaging multiple antibody particles together. Through individual-particle electron tomography (IPET) 3D reconstruction from negatively-stained ET images, we obtained 120 ab-initio 3D density maps at an intermediate resolution (~1–3 nm) from 120 individual IgG1 antibody particles. Using these maps as a constraint, we derived 120 conformations of the antibody via structural flexible docking of the crystal structure to these maps by targeted molecular dynamics simulations. Statistical analysis of the various conformations disclosed the antibody 3D conformational flexibility through the distribution of its domain distances and orientations. This blueprint approach, if extended to other flexible proteins, may serve as a useful methodology towards understanding protein dynamics and functions.

  10. 3D structural fluctuation of IgG1 antibody revealed by individual particle electron tomography

    DOE PAGESBeta

    Zhang, Xing; Zhang, Lei; Tong, Huimin; Peng, Bo; Rames, Matthew J.; Zhang, Shengli; Ren, Gang

    2015-05-05

    Commonly used methods for determining protein structure, including X-ray crystallography and single-particle reconstruction, often provide a single and unique three-dimensional (3D) structure. However, in these methods, the protein dynamics and flexibility/fluctuation remain mostly unknown. Here, we utilized advances in electron tomography (ET) to study the antibody flexibility and fluctuation through structural determination of individual antibody particles rather than averaging multiple antibody particles together. Through individual-particle electron tomography (IPET) 3D reconstruction from negatively-stained ET images, we obtained 120 ab-initio 3D density maps at an intermediate resolution (~1–3 nm) from 120 individual IgG1 antibody particles. Using these maps as a constraint, wemore » derived 120 conformations of the antibody via structural flexible docking of the crystal structure to these maps by targeted molecular dynamics simulations. Statistical analysis of the various conformations disclosed the antibody 3D conformational flexibility through the distribution of its domain distances and orientations. This blueprint approach, if extended to other flexible proteins, may serve as a useful methodology towards understanding protein dynamics and functions.« less

  11. 3D Structural Fluctuation of IgG1 Antibody Revealed by Individual Particle Electron Tomography

    PubMed Central

    Zhang, Xing; Zhang, Lei; Tong, Huimin; Peng, Bo; Rames, Matthew J.; Zhang, Shengli; Ren, Gang

    2015-01-01

    Commonly used methods for determining protein structure, including X-ray crystallography and single-particle reconstruction, often provide a single and unique three-dimensional (3D) structure. However, in these methods, the protein dynamics and flexibility/fluctuation remain mostly unknown. Here, we utilized advances in electron tomography (ET) to study the antibody flexibility and fluctuation through structural determination of individual antibody particles rather than averaging multiple antibody particles together. Through individual-particle electron tomography (IPET) 3D reconstruction from negatively-stained ET images, we obtained 120 ab-initio 3D density maps at an intermediate resolution (~1–3 nm) from 120 individual IgG1 antibody particles. Using these maps as a constraint, we derived 120 conformations of the antibody via structural flexible docking of the crystal structure to these maps by targeted molecular dynamics simulations. Statistical analysis of the various conformations disclosed the antibody 3D conformational flexibility through the distribution of its domain distances and orientations. This blueprint approach, if extended to other flexible proteins, may serve as a useful methodology towards understanding protein dynamics and functions. PMID:25940394

  12. ORBXYZ: a 3D single-particle orbit code for following charged-particle trajectories in equilibrium magnetic fields

    SciTech Connect

    Anderson, D.V.; Cohen, R.H.; Ferguson, J.R.; Johnston, B.M.; Sharp, C.B.; Willmann, P.A.

    1981-06-30

    The single particle orbit code, TIBRO, has been modified extensively to improve the interpolation methods used and to allow use of vector potential fields in the simulation of charged particle orbits on a 3D domain. A 3D cubic B-spline algorithm is used to generate spline coefficients used in the interpolation. Smooth and accurate field representations are obtained. When vector potential fields are used, the 3D cubic spline interpolation formula analytically generates the magnetic field used to push the particles. This field has del.BETA = 0 to computer roundoff. When magnetic induction is used the interpolation allows del.BETA does not equal 0, which can lead to significant nonphysical results. Presently the code assumes quadrupole symmetry, but this is not an essential feature of the code and could be easily removed for other applications. Many details pertaining to this code are given on microfiche accompanying this report.

  13. Continuous dielectrophoretic particle separation using a microfluidic device with 3D electrodes and vaulted obstacles.

    PubMed

    Jia, Yankai; Ren, Yukun; Jiang, Hongyuan

    2015-08-01

    This paper reports a microfluidic separation device combining 3D electrodes and vaulted obstacles to continuously separate particles experiencing strong positive dielectrophoresis (DEP) from particles experiencing weak positive DEP, or from particles experiencing negative DEP. The separation is achieved by first focusing the particle mixture into a narrow stream by a hydrodynamic focusing flow, and then deviating them into different outlets by AC DEP. The vaulted obstacles facilitate the separation by both increasing the non-uniformity of the electric field, and influencing the particles to move in regions strongly affected by DEP. The 3D electrodes give rise to a spatially non-uniform electric field and extend DEP effect to the channel height. Numerical simulations are performed to investigate the effects of the obstacles on electric field distribution and particle trajectories so as to optimize the obstacle height and compare with the experimental results. The performance of the device is assessed by separating 25 μm gold-coated particles from 10 μm particles in different flow rates by positive DEP and negative DEP, and also separating 25 μm gold-coated particles from yeast cells using only positive DEP. The experimental observation shows a reasonable agreement with numerical simulation results. PMID:25962351

  14. VISRAD, 3-D Target Design and Radiation Simulation Code

    NASA Astrophysics Data System (ADS)

    Li, Yingjie; Macfarlane, Joseph; Golovkin, Igor

    2015-11-01

    The 3-D view factor code VISRAD is widely used in designing HEDP experiments at major laser and pulsed-power facilities, including NIF, OMEGA, OMEGA-EP, ORION, LMJ, Z, and PLX. It simulates target designs by generating a 3-D grid of surface elements, utilizing a variety of 3-D primitives and surface removal algorithms, and can be used to compute the radiation flux throughout the surface element grid by computing element-to-element view factors and solving power balance equations. Target set-up and beam pointing are facilitated by allowing users to specify positions and angular orientations using a variety of coordinates systems (e.g., that of any laser beam, target component, or diagnostic port). Analytic modeling for laser beam spatial profiles for OMEGA DPPs and NIF CPPs is used to compute laser intensity profiles throughout the grid of surface elements. We will discuss recent improvements to the software package and plans for future developments.

  15. 3D finite element simulations of high velocity projectile impact

    NASA Astrophysics Data System (ADS)

    Ožbolt, Joško; İrhan, Barış; Ruta, Daniela

    2015-09-01

    An explicit three-dimensional (3D) finite element (FE) code is developed for the simulation of high velocity impact and fragmentation events. The rate sensitive microplane material model, which accounts for large deformations and rate effects, is used as a constitutive law. In the code large deformation frictional contact is treated by forward incremental Lagrange multiplier method. To handle highly distorted and damaged elements the approach based on the element deletion is employed. The code is then used in 3D FE simulations of high velocity projectile impact. The results of the numerical simulations are evaluated and compared with experimental results. It is shown that it realistically predicts failure mode and exit velocities for different geometries of plain concrete slab. Moreover, the importance of some relevant parameters, such as contact friction, rate sensitivity, bulk viscosity and deletion criteria are addressed.

  16. 3-D simulation of nanopore structure for DNA sequencing.

    PubMed

    Park, Jun-Mo; Pak, Y Eugene; Chun, Honggu; Lee, Jong-Ho

    2012-07-01

    In this paper, we propose a method for simulating nanopore structure by using conventional 3-D simulation tool to mimic the I-V behavior of the nanopore structure. In the simulation, we use lightly doped silicon for ionic solution where some parameters like electron affinity and dielectric constant are fitted to consider the ionic solution. By using this method, we can simulate the I-V behavior of nanopore structure depending on the location and the size of the sphere shaped silicon oxide which is considered to be an indicator of a DNA base. In addition, we simulate an Ionic Field Effect Transistor (IFET) which has basically the nanopore structure, and show that the simulated curves follow sufficiently the I-V behavior of the measurement data. Therefore, we think it is reasonable to apply parameter modeling mentioned above to simulate nanopore structure. The key idea is to modify electron affinity of silicon which is used to mimic the KCl solution to avoid band bending and depletion inside the nanopore. We could efficiently utilize conventional 3-D simulation tool to simulate the I-V behavior of nanopore structures. PMID:22966538

  17. Using the symmetries in 3D Smoothed Particle Hydrodynamics

    NASA Astrophysics Data System (ADS)

    Cruz-Pérez, Juan P.; González, José A.

    2011-10-01

    We apply the cartoon SPH technique to solve numerically the Newtonian Euler's equations in scenarios with spherical symmetry with the standard Smoothed Particle Hydrodynamics method, using the equations written in cartesian coordinates. This implementation allow us to increase the resolution of the simulations in order to obtain accurate results. We test it using the Shock tube and the isothermal collapse of a perfect fluid in spherical symmetry.

  18. Fast 3D stereo flood simulations in urban areas

    NASA Astrophysics Data System (ADS)

    Hoes, O.; de Haan, G.; Stelling, G.; van Leeuwen, E.; van Dam, A.; Pleumeekers, O.; Schuurmans, W.

    2012-04-01

    Flood propagation models are essential to study floods as it is problematic to collect data during actual floods. Moreover, models are needed to explore the consequences of additional scenarios above the actual flood itself. The results of these model studies are generally graphs with water levels over time for certain locations or maps with the flood extent in an area for different return periods. A main difficulty with these depictions of flood information is that they do not reflect the seriousness of flood impacts well in terms of life-like human experience. Typically, one needs a (near) flood before measures are implemented. Apparently, a graph or map is not the proper material to convince politicians and policy makers, even if they live in the threatened area. The recent introduction of commercially available 3D stereo projectors and high resolution elevation data make it possible to build life-like visualizations of simulations. In our research we explored using 3D stereo, the recently collected elevation data of the Netherlands (20 laser points per m2!) in combination with aerial photographs, and a new fast 2D flood propagation calculation scheme. This scheme (under construction) is able to simulate floods using such high amounts of data points. The model simulates flood propagation on an irregular grid; at locations with large elevation differences (e.g. in urban areas) and fast flowing water, smaller cells are used compared to flat surfaces where the water is not or hardly flowing. The result of our combination is a very detailed flood simulation model that can be used to simulate floods within a fraction of the current calculation time. The opportunities of models and their results increase enormously with fast calculations and visualizations combined. For example, the model allows on the spot exploration of measures during a flood, with the 3D visualization ensuring that flood impacts become clear for decision makers. We will show the preliminary

  19. Automated 3D trajectory measuring of large numbers of moving particles.

    PubMed

    Wu, Hai Shan; Zhao, Qi; Zou, Danping; Chen, Yan Qiu

    2011-04-11

    Complex dynamics of natural particle systems, such as insect swarms, bird flocks, fish schools, has attracted great attention of scientists for years. Measuring 3D trajectory of each individual in a group is vital for quantitative study of their dynamic properties, yet such empirical data is rare mainly due to the challenges of maintaining the identities of large numbers of individuals with similar visual features and frequent occlusions. We here present an automatic and efficient algorithm to track 3D motion trajectories of large numbers of moving particles using two video cameras. Our method solves this problem by formulating it as three linear assignment problems (LAP). For each video sequence, the first LAP obtains 2D tracks of moving targets and is able to maintain target identities in the presence of occlusions; the second one matches the visually similar targets across two views via a novel technique named maximum epipolar co-motion length (MECL), which is not only able to effectively reduce matching ambiguity but also further diminish the influence of frequent occlusions; the last one links 3D track segments into complete trajectories via computing a globally optimal assignment based on temporal and kinematic cues. Experiment results on simulated particle swarms with various particle densities validated the accuracy and robustness of the proposed method. As real-world case, our method successfully acquired 3D flight paths of fruit fly (Drosophila melanogaster) group comprising hundreds of freely flying individuals. PMID:21503074

  20. Novel Discrete Element Method for 3D non-spherical granular particles.

    NASA Astrophysics Data System (ADS)

    Seelen, Luuk; Padding, Johan; Kuipers, Hans

    2015-11-01

    Granular materials are common in many industries and nature. The different properties from solid behavior to fluid like behavior are well known but less well understood. The main aim of our work is to develop a discrete element method (DEM) to simulate non-spherical granular particles. The non-spherical shape of particles is important, as it controls the behavior of the granular materials in many situations, such as static systems of packed particles. In such systems the packing fraction is determined by the particle shape. We developed a novel 3D discrete element method that simulates the particle-particle interactions for a wide variety of shapes. The model can simulate quadratic shapes such as spheres, ellipsoids, cylinders. More importantly, any convex polyhedron can be used as a granular particle shape. These polyhedrons are very well suited to represent non-rounded sand particles. The main difficulty of any non-spherical DEM is the determination of particle-particle overlap. Our model uses two iterative geometric algorithms to determine the overlap. The algorithms are robust and can also determine multiple contact points which can occur for these shapes. With this method we are able to study different applications such as the discharging of a hopper or silo. Another application the creation of a random close packing, to determine the solid volume fraction as a function of the particle shape.

  1. Comparative visual analysis of 3D urban wind simulations

    NASA Astrophysics Data System (ADS)

    Röber, Niklas; Salim, Mohamed; Grawe, David; Leitl, Bernd; Böttinger, Michael; Schlünzen, Heinke

    2016-04-01

    Climate simulations are conducted in large quantity for a variety of different applications. Many of these simulations focus on global developments and study the Earth's climate system using a coupled atmosphere ocean model. Other simulations are performed on much smaller regional scales, to study very small fine grained climatic effects. These microscale climate simulations pose similar, yet also different, challenges for the visualization and the analysis of the simulation data. Modern interactive visualization and data analysis techniques are very powerful tools to assist the researcher in answering and communicating complex research questions. This presentation discusses comparative visualization for several different wind simulations, which were created using the microscale climate model MITRAS. The simulations differ in wind direction and speed, but are all centered on the same simulation domain: An area of Hamburg-Wilhelmsburg that hosted the IGA/IBA exhibition in 2013. The experiments contain a scenario case to analyze the effects of single buildings, as well as examine the impact of the Coriolis force within the simulation. The scenario case is additionally compared with real measurements from a wind tunnel experiment to ascertain the accuracy of the simulation and the model itself. We also compare different approaches for tree modeling and evaluate the stability of the model. In this presentation, we describe not only our workflow to efficiently and effectively visualize microscale climate simulation data using common 3D visualization and data analysis techniques, but also discuss how to compare variations of a simulation and how to highlight the subtle differences in between them. For the visualizations we use a range of different 3D tools that feature techniques for statistical data analysis, data selection, as well as linking and brushing.

  2. DREAM3D simulations of inner-belt dynamics

    SciTech Connect

    Cunningham, Gregory Scott

    2015-05-26

    A 1973 paper by Lyons and Thorne explains the two-belt structure for electrons in the inner magnetosphere as a balance between inward radial diffusion and loss to the atmosphere, where the loss to the atmosphere is enabled by pitch-angle scattering from Coulomb and wave-particle interactions. In the 1973 paper, equilibrium solutions to a decoupled set of 1D radial diffusion equations, one for each value of the first invariant of motion, μ, were computed to produce the equilibrium two-belt structure. Each 1D radial diffusion equation incorporated an L-and μ-dependent `lifetime' due to the Coulomb and wave-particle interactions. This decoupling of the problem is appropriate under the assumption that radial diffusion is slow in comparison to pitch-angle scattering. However, for some values of μ and L the lifetime associated with pitch-angle scattering is comparable to the timescale associated with radial diffusion, suggesting that the true equilibrium solutions might reflect `coupled modes' involving pitch-angle scattering and radial diffusion and thus requiring a 3D diffusion model. In the work we show here, we have computed the equilibrium solutions using our 3D diffusion model, DREAM3D, that allows for such coupling. We find that the 3D equilibrium solutions are quite similar to the solutions shown in the 1973 paper when we use the same physical models for radial diffusion and pitch-angle scattering from hiss. However, we show that the equilibrium solutions are quite sensitive to various aspects of the physics model employed in the 1973 paper that can be improved, suggesting that additional work needs to be done to understand the two-belt structure.

  3. 3D Smoothed Particle Hydrodynamics Models of Betelgeuse's Bow Shock

    NASA Astrophysics Data System (ADS)

    Mohamed, S.; Mackey, J.; Langer, N.

    2013-05-01

    Betelgeuse, the bright red supergiant (RSG) in Orion, is a runaway star. Its supersonic motion through the interstellar medium has resulted in the formation of a bow shock, a cometary structure pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. We show that the bow shock morphology depends substantially on the growth timescale for Rayleigh-Taylor versus Kelvin-Helmholtz instabilities. We discuss our models in light of the recent Herschel, GALEX and VLA observations. If the mass in the bow shock shell is low (~few × 10-3 M⊙), as seems to be implied by the AKARI and Herschel observations, then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular, smooth bow shock shell is consistent with this conclusion. We further discuss the implications of our results, in particular, the possibility that Betelgeuse may have only recently entered the RSG phase.

  4. 3D simulation of the Cluster-Cluster Aggregation model

    NASA Astrophysics Data System (ADS)

    Li, Chao; Xiong, Hailing

    2014-12-01

    We write a program to implement the Cluster-Cluster Aggregation (CCA) model with java programming language. By using the simulation program, the fractal aggregation growth process can be displayed dynamically in the form of a three-dimensional (3D) figure. Meanwhile, the related kinetics data of aggregation simulation can be also recorded dynamically. Compared to the traditional programs, the program has better real-time performance and is more helpful to observe the fractal growth process, which contributes to the scientific study in fractal aggregation. Besides, because of adopting java programming language, the program has very good cross-platform performance.

  5. 3D Simulation of External Flooding Events for the RISMC Pathway

    SciTech Connect

    Prescott, Steven; Mandelli, Diego; Sampath, Ramprasad; Smith, Curtis; Lin, Linyu

    2015-09-01

    Incorporating 3D simulations as part of the Risk-Informed Safety Margins Characterization (RISMIC) Toolkit allows analysts to obtain a more complete picture of complex system behavior for events including external plant hazards. External events such as flooding have become more important recently – however these can be analyzed with existing and validated simulated physics toolkits. In this report, we describe these approaches specific to flooding-based analysis using an approach called Smoothed Particle Hydrodynamics. The theory, validation, and example applications of the 3D flooding simulation are described. Integrating these 3D simulation methods into computational risk analysis provides a spatial/visual aspect to the design, improves the realism of results, and can prove visual understanding to validate the analysis of flooding.

  6. Virtual environment display for a 3D audio room simulation

    NASA Technical Reports Server (NTRS)

    Chapin, William L.; Foster, Scott H.

    1992-01-01

    The development of a virtual environment simulation system integrating a 3D acoustic audio model with an immersive 3D visual scene is discussed. The system complements the acoustic model and is specified to: allow the listener to freely move about the space, a room of manipulable size, shape, and audio character, while interactively relocating the sound sources; reinforce the listener's feeling of telepresence in the acoustical environment with visual and proprioceptive sensations; enhance the audio with the graphic and interactive components, rather than overwhelm or reduce it; and serve as a research testbed and technology transfer demonstration. The hardware/software design of two demonstration systems, one installed and one portable, are discussed through the development of four iterative configurations.

  7. 3D-CANVENT: An interactive mine ventilation simulator

    SciTech Connect

    Hardcastle, S.G.

    1995-12-31

    3D-CANVENT is a software package that integrates advanced computer aided design (ACAD) true 3D graphics with a mine ventilation simulator. The package runs as a Windows{trademark} application to access its printer drivers environment and does not need third party CAD software. It is composed of two primary modules: DMVENT and MINEDESIGNER. DMVENT is a traditional Fortran coded Hardy-Cross iterative ventilation network solver written in 1980 with thermodynamic capabilities. This module is relatively unchanged with the traditional data input options for branch type, specified or calculated resistances, fixed flows, and fixed or variable pressure fans. MINEDESIGNER is the graphics engine that optimizes the ventilation design process. It performs the front-end transformation of input data entered in the graphical interface into the correct format for the solver. At the back-end it reconverts the historically standard tabular data output from the solver into an easily viewed graphical format. ACAD features of MINEDESIGNER are used to generate a 3D wire-frame node and branch network of the mine`s ventilation system. The network can be displayed in up to 4 views orientated to XYZ planes or a 3D view. AU the views have zoom, pan, slice and rotate options. The graphical interface efficiently permits data entry and editing via a mouse with pick-and-point item selection. Branches can be found or added with {open_quotes}search{close_quotes} and {open_quotes}join{close_quotes} options. Visual interpretation is enhanced by the 16 colour options for branches and numerous graphical attributes. Network locations are readily identified by alpha-numeric names for branches, junctions and fans, and also the logical numbering of junctions. The program is also readily expandable for pollutant simulation and control/monitoring applications.

  8. Acoustic simulation in realistic 3D virtual scenes

    NASA Astrophysics Data System (ADS)

    Gozard, Patrick; Le Goff, Alain; Naz, Pierre; Cathala, Thierry; Latger, Jean

    2003-09-01

    The simulation workshop CHORALE developed in collaboration with OKTAL SE company for the French MoD is used by government services and industrial companies for weapon system validation and qualification trials in the infrared domain. The main operational reference for CHORALE is the assessment of the infrared guidance system of the Storm Shadow missile French version, called Scalp. The use of CHORALE workshop is now extended to the acoustic domain. The main objective is the simulation of the detection of moving vehicles in realistic 3D virtual scenes. This article briefly describes the acoustic model in CHORALE. The 3D scene is described by a set of polygons. Each polygon is characterized by its acoustic resistivity or its complex impedance. Sound sources are associated with moving vehicles and are characterized by their spectra and directivities. A microphone sensor is defined by its position, its frequency band and its sensitivity. The purpose of the acoustic simulation is to calculate the incoming acoustic pressure on microphone sensors. CHORALE is based on a generic ray tracing kernel. This kernel possesses original capabilities: computation time is nearly independent on the scene complexity, especially the number of polygons, databases are enhanced with precise physical data, special mechanisms of antialiasing have been developed that enable to manage very accurate details. The ray tracer takes into account the wave geometrical divergence and the atmospheric transmission. The sound wave refraction is simulated and rays cast in the 3D scene are curved according to air temperature gradient. Finally, sound diffraction by edges (hill, wall,...) is also taken into account.

  9. 3-D Imaging and Simulation for Nephron Sparing Surgical Training.

    PubMed

    Ahmadi, Hamed; Liu, Jen-Jane

    2016-08-01

    Minimally invasive partial nephrectomy (MIPN) is now considered the procedure of choice for small renal masses largely based on functional advantages over traditional open surgery. Lack of haptic feedback, the need for spatial understanding of tumor borders, and advanced operative techniques to minimize ischemia time or achieve zero-ischemia PN are among factors that make MIPN a technically demanding operation with a steep learning curve for inexperienced surgeons. Surgical simulation has emerged as a useful training adjunct in residency programs to facilitate the acquisition of these complex operative skills in the setting of restricted work hours and limited operating room time and autonomy. However, the majority of available surgical simulators focus on basic surgical skills, and procedure-specific simulation is needed for optimal surgical training. Advances in 3-dimensional (3-D) imaging have also enhanced the surgeon's ability to localize tumors intraoperatively. This article focuses on recent procedure-specific simulation models for laparoscopic and robotic-assisted PN and advanced 3-D imaging techniques as part of pre- and some cases, intraoperative surgical planning. PMID:27314271

  10. Using Delft3D to Simulate Current Energy Conversion

    NASA Astrophysics Data System (ADS)

    James, S. C.; Chartrand, C.; Roberts, J.

    2015-12-01

    As public concern with renewable energy increases, current energy conversion (CEC) technology is being developed to optimize energy output and minimize environmental impact. CEC turbines generate energy from tidal and current systems and create wakes that interact with turbines located downstream of a device. The placement of devices can greatly influence power generation and structural reliability. CECs can also alter the ecosystem process surrounding the turbines, such as flow regimes, sediment dynamics, and water quality. Software is needed to investigate specific CEC sites to simulate power generation and hydrodynamic responses of a flow through a CEC turbine array. This work validates Delft3D against several flume experiments by simulating the power generation and hydrodynamic response of flow through a turbine or actuator disc(s). Model parameters are then calibrated against these data sets to reproduce momentum removal and wake recovery data with 3-D flow simulations. Simulated wake profiles and turbulence intensities compare favorably to the experimental data and demonstrate the utility and accuracy of a fast-running tool for future siting and analysis of CEC arrays in complex domains.

  11. A Bayesian approach for suppression of limited angular sampling artifacts in single particle 3D reconstruction.

    PubMed

    Moriya, Toshio; Acar, Erman; Cheng, R Holland; Ruotsalainen, Ulla

    2015-09-01

    In the single particle reconstruction, the initial 3D structure often suffers from the limited angular sampling artifact. Selecting 2D class averages of particle images generally improves the accuracy and efficiency of the reference-free 3D angle estimation, but causes an insufficient angular sampling to fill the information of the target object in the 3D frequency space. Similarly, the initial 3D structure by the random-conical tilt reconstruction has the well-known "missing cone" artifact. Here, we attempted to solve the limited angular sampling problem by sequentially applying maximum a posteriori estimate with expectation maximization algorithm (sMAP-EM). Using both simulated and experimental cryo-electron microscope images, the sMAP-EM was compared to the direct Fourier method on the basis of reconstruction error and resolution. To establish selection criteria of the final regularization weight for the sMAP-EM, the effects of noise level and sampling sparseness on the reconstructions were examined with evenly distributed sampling simulations. The frequency information filled in the missing cone of the conical tilt sampling simulations was assessed by developing new quantitative measurements. All the results of visual and numerical evaluations showed the sMAP-EM performed better than the direct Fourier method, regardless of the sampling method, noise level, and sampling sparseness. Furthermore, the frequency domain analysis demonstrated that the sMAP-EM can fill the meaningful information in the unmeasured angular space without detailed a priori knowledge of the objects. The current research demonstrated that the sMAP-EM has a high potential to facilitate the determination of 3D protein structures at near atomic-resolution. PMID:26193484

  12. 3D simulation of plant and living tissue superficial lesions

    NASA Astrophysics Data System (ADS)

    Bratchenko, Ivan A.; Sindyaeva, Alexandra R.; Zakharov, Valery P.

    2008-06-01

    The analytic schemes of calculated absorbed and scattered radiation spatial distribution in multilayer plant and living tissues and diagnostic of their physical state are presented. The correct realization of these tasks was obtained with 3D Monte Carlo simulation of optical radiation propagation through multiple scattering medium in TracePro environment. Analysis of simulation data was made by differential backscattering method, which allows to investigate general backscattered radiation dependences on optical and geometrical parameters of living tissue. It was shown that obtained results formed the basis for developing an algorithm of optical superficial inhomogeneous registration and spatial localization. Such diagnosis can be executed in tissues of any arbitrary surface structure. Designed scheme is intended to utilize in contactless macro diagnostics device. The same approach was used for simulation of optical spectra of healthy and diseased virtual leaves for plant tissue pathological changes revealing.

  13. 3-D simulations of solar and stellar convection and magnetoconvection

    NASA Astrophysics Data System (ADS)

    Nordlund, Å.; Stein, R. F.

    1990-05-01

    We present the key components of a 3-D code designed for simulating the hydrodynamics and magnetohydrodynamics of stellar atmospheres and envelopes. Some particular properties of the code are: (1) the ability to handle strong stratification (extensive simulations with bottom/top pressure ratios of 3×104 have been performed, and simulations with pressure ratios of 5×106 are being initiated); (2) a detailed treatment of the radiating surface; (3) a functional form of the subgrid-scale diffusion designed to minimize the influence on resolved motions; (4) boundary conditions open to flows. The top boundary allows the transmission of short period waves, while the bottom boundary condition was designed to enforce a displacement node for radial pressure modes.

  14. Unstructured 3D grid toolbox for modeling and simulation

    SciTech Connect

    George, D.

    1997-11-01

    Computable 3D grids that accurately represent complex multimaterial geometries are essential for both static and time-dependent modeling and simulation. LaGriT, the grid toolbox developed at Los Alamos provides a sophisticated set of initial grid generation, grid maintenance and grid optimization tools. We present example grids that demonstrate the flexibility of the grid generator. Additionally, we present the results of an electrostatic calculation and a grain growth problem that illustrate the grid optimization features and the utility of the grid server architecture.

  15. Numerical simulation of vortex breakdown via 3-D Euler equations

    NASA Astrophysics Data System (ADS)

    Le, T. H.; Mege, P.; Morchoisne, Y.

    1990-06-01

    The long term goal is the modeling of vortex breakdown that occurs in some aerodynamic configurations at high angle of attack, (i.e., fighters with highly swept delta wings or missiles). A numerical simulation was made based on solving the 3-D Euler equations for an usteady incompressible flow. Preliminary results were obtained using a pressure-velocity formulation with periodic boundary conditions, the Euler equations being discretized by 2nd order finite difference schemes. The continuation to this work by implementing more realistic boundary conditions and 4th order finite difference discretization schemes are presented.

  16. Finite-element 3D simulation tools for high-current relativistic electron beams

    NASA Astrophysics Data System (ADS)

    Humphries, Stanley; Ekdahl, Carl

    2002-08-01

    The DARHT second-axis injector is a challenge for computer simulations. Electrons are subject to strong beam-generated forces. The fields are fully three-dimensional and accurate calculations at surfaces are critical. We describe methods applied in OmniTrak, a 3D finite-element code suite that can address DARHT and the full range of charged-particle devices. The system handles mesh generation, electrostatics, magnetostatics and self-consistent particle orbits. The MetaMesh program generates meshes of conformal hexahedrons to fit any user geometry. The code has the unique ability to create structured conformal meshes with cubic logic. Organized meshes offer advantages in speed and memory utilization in the orbit and field solutions. OmniTrak is a versatile charged-particle code that handles 3D electric and magnetic field solutions on independent meshes. The program can update both 3D field solutions from the calculated beam space-charge and current-density. We shall describe numerical methods for orbit tracking on a hexahedron mesh. Topics include: 1) identification of elements along the particle trajectory, 2) fast searches and adaptive field calculations, 3) interpolation methods to terminate orbits on material surfaces, 4) automatic particle generation on multiple emission surfaces to model space-charge-limited emission and field emission, 5) flexible Child law algorithms, 6) implementation of the dual potential model for 3D magnetostatics, and 7) assignment of charge and current from model particle orbits for self-consistent fields.

  17. 3D measurement of the position of gold particles via evanescent digital holographic particle tracking velocimetry

    NASA Astrophysics Data System (ADS)

    Satake, Shin-ichi; Unno, Noriyuki; Nakata, Shuichiro; Taniguchi, Jun

    2016-08-01

    A new technique based on digital holography and evanescent waves was developed for 3D measurements of the position of gold nanoparticles in water. In this technique, an intensity profile is taken from a holographic image of a gold particle. To detect the position of the gold particle with high accuracy, its holographic image is recorded on a nanosized step made of MEXFLON, which has a refractive index close to that of water, and the position of the particle is reconstructed by means of digital holography. The height of the nanosized step was measured by using a profilometer and the digitally reconstructed height of the glass substrate had good agreement with the measured value. Furthermore, this method can be used to accurately track the 3D position of a gold particle in water.

  18. Advanced system for 3D dental anatomy reconstruction and 3D tooth movement simulation during orthodontic treatment

    NASA Astrophysics Data System (ADS)

    Monserrat, Carlos; Alcaniz-Raya, Mariano L.; Juan, M. Carmen; Grau Colomer, Vincente; Albalat, Salvador E.

    1997-05-01

    This paper describes a new method for 3D orthodontics treatment simulation developed for an orthodontics planning system (MAGALLANES). We develop an original system for 3D capturing and reconstruction of dental anatomy that avoid use of dental casts in orthodontic treatments. Two original techniques are presented, one direct in which data are acquired directly form patient's mouth by mean of low cost 3D digitizers, and one mixed in which data are obtained by 3D digitizing of hydrocollids molds. FOr this purpose we have designed and manufactured an optimized optical measuring system based on laser structured light. We apply these 3D dental models to simulate 3D movement of teeth, including rotations, during orthodontic treatment. The proposed algorithms enable to quantify the effect of orthodontic appliance on tooth movement. The developed techniques has been integrated in a system named MAGALLANES. This original system present several tools for 3D simulation and planning of orthodontic treatments. The prototype system has been tested in several orthodontic clinic with very good results.

  19. 3D FEM Simulation of Flank Wear in Turning

    NASA Astrophysics Data System (ADS)

    Attanasio, Aldo; Ceretti, Elisabetta; Giardini, Claudio

    2011-05-01

    This work deals with tool wear simulation. Studying the influence of tool wear on tool life, tool substitution policy and influence on final part quality, surface integrity, cutting forces and power consumption it is important to reduce the global process costs. Adhesion, abrasion, erosion, diffusion, corrosion and fracture are some of the phenomena responsible of the tool wear depending on the selected cutting parameters: cutting velocity, feed rate, depth of cut, …. In some cases these wear mechanisms are described by analytical models as a function of process variables (temperature, pressure and sliding velocity along the cutting surface). These analytical models are suitable to be implemented in FEM codes and they can be utilized to simulate the tool wear. In the present paper a commercial 3D FEM software has been customized to simulate the tool wear during turning operations when cutting AISI 1045 carbon steel with uncoated tungsten carbide tip. The FEM software was improved by means of a suitable subroutine able to modify the tool geometry on the basis of the estimated tool wear as the simulation goes on. Since for the considered couple of tool-workpiece material the main phenomena generating wear are the abrasive and the diffusive ones, the tool wear model implemented into the subroutine was obtained as combination between the Usui's and the Takeyama and Murata's models. A comparison between experimental and simulated flank tool wear curves is reported demonstrating that it is possible to simulate the tool wear development.

  20. Simulation of 3D Global Wave Propagation Through Geodynamic Models

    NASA Astrophysics Data System (ADS)

    Schuberth, B.; Piazzoni, A.; Bunge, H.; Igel, H.; Steinle-Neumann, G.

    2005-12-01

    This project aims at a better understanding of the forward problem of global 3D wave propagation. We use the spectral element program "SPECFEM3D" (Komatitsch and Tromp, 2002a,b) with varying input models of seismic velocities derived from mantle convection simulations (Bunge et al., 2002). The purpose of this approach is to obtain seismic velocity models independently from seismological studies. In this way one can test the effects of varying parameters of the mantle convection models on the seismic wave field. In order to obtain the seismic velocities from the temperature field of the geodynamical simulations we follow a mineral physics approach. Assuming a certain mantle composition (e.g. pyrolite with CMASF composition) we compute the stable phases for each depth (i.e. pressure) and temperature by system Gibbs free energy minimization. Elastic moduli and density are calculated from the equations of state of the stable mineral phases. For this we use a mineral physics database derived from calorimetric experiments (enthalphy and entropy of formation, heat capacity) and EOS parameters.

  1. GBS: Global 3D simulation of tokamak edge region

    NASA Astrophysics Data System (ADS)

    Zhu, Ben; Fisher, Dustin; Rogers, Barrett; Ricci, Paolo

    2012-10-01

    A 3D two-fluid global code, namely Global Braginskii Solver (GBS), is being developed to explore the physics of turbulent transport, confinement, self-consistent profile formation, pedestal scaling and related phenomena in the edge region of tokamaks. Aimed at solving drift-reduced Braginskii equations [1] in complex magnetic geometry, the GBS is used for turbulence simulation in SOL region. In the recent upgrade, the simulation domain is expanded into close flux region with twist-shift boundary conditions. Hence, the new GBS code is able to explore global transport physics in an annular full-torus domain from the top of the pedestal into the far SOL. We are in the process of identifying and analyzing the linear and nonlinear instabilities in the system using the new GBS code. Preliminary results will be presented and compared with other codes if possible.[4pt] [1] A. Zeiler, J. F. Drake and B. Rogers, Phys. Plasmas 4, 2134 (1997)

  2. Coniferous Canopy BRF Simulation Based on 3-D Realistic Scene

    NASA Technical Reports Server (NTRS)

    Wang, Xin-yun; Guo, Zhi-feng; Qin, Wen-han; Sun, Guo-qing

    2011-01-01

    It is difficulties for the computer simulation method to study radiation regime at large-scale. Simplified coniferous model was investigate d in the present study. It makes the computer simulation methods such as L-systems and radiosity-graphics combined method (RGM) more powerf ul in remote sensing of heterogeneous coniferous forests over a large -scale region. L-systems is applied to render 3-D coniferous forest scenarios: and RGM model was used to calculate BRF (bidirectional refle ctance factor) in visible and near-infrared regions. Results in this study show that in most cases both agreed well. Meanwhiie at a tree and forest level. the results are also good.

  3. Code System to Simulate 3D Tracer Dispersion in Atmosphere.

    Energy Science and Technology Software Center (ESTSC)

    2002-01-25

    Version 00 SHREDI is a shielding code system which executes removal-diffusion computations for bi-dimensional shields in r-z or x-y geometries. It may also deal with monodimensional problems (infinitely high cylinders or slabs). MESYST can simulate 3D tracer dispersion in the atmosphere. Three programs are part of this system: CRE_TOPO prepares the terrain data for MESYST. NOABL calculates three-dimensional free divergence windfields over complex terrain. PAS computes tracer concentrations and depositions on a given domain. Themore » purpose of this work is to develop a reliable simulation tool for pollutant atmospheric dispersion, which gives a realistic approach and allows one to compute the pollutant concentrations over complex terrains with good accuracy. The factional brownian model, which furnishes more accurate concentration values, is introduced to calculate pollutant atmospheric dispersion. The model was validated on SIESTA international experiments.« less

  4. Coniferous canopy BRF simulation based on 3-D realistic scene.

    PubMed

    Wang, Xin-Yun; Guo, Zhi-Feng; Qin, Wen-Han; Sun, Guo-Qing

    2011-09-01

    It is difficulties for the computer simulation method to study radiation regime at large-scale. Simplified coniferous model was investigated in the present study. It makes the computer simulation methods such as L-systems and radiosity-graphics combined method (RGM) more powerful in remote sensing of heterogeneous coniferous forests over a large-scale region. L-systems is applied to render 3-D coniferous forest scenarios, and RGM model was used to calculate BRF (bidirectional reflectance factor) in visible and near-infrared regions. Results in this study show that in most cases both agreed well. Meanwhile at a tree and forest level, the results are also good. PMID:22097856

  5. A Parallelized 3D Particle-In-Cell Method With Magnetostatic Field Solver And Its Applications

    NASA Astrophysics Data System (ADS)

    Hsu, Kuo-Hsien; Chen, Yen-Sen; Wu, Men-Zan Bill; Wu, Jong-Shinn

    2008-10-01

    A parallelized 3D self-consistent electrostatic particle-in-cell finite element (PIC-FEM) code using an unstructured tetrahedral mesh was developed. For simulating some applications with external permanent magnet set, the distribution of the magnetostatic field usually also need to be considered and determined accurately. In this paper, we will firstly present the development of a 3D magnetostatic field solver with an unstructured mesh for the flexibility of modeling objects with complex geometry. The vector Poisson equation for magnetostatic field is formulated using the Galerkin nodal finite element method and the resulting matrix is solved by parallel conjugate gradient method. A parallel adaptive mesh refinement module is coupled to this solver for better resolution. Completed solver is then verified by simulating a permanent magnet array with results comparable to previous experimental observations and simulations. By taking the advantage of the same unstructured grid format of this solver, the developed PIC-FEM code could directly and easily read the magnetostatic field for particle simulation. In the upcoming conference, magnetron is simulated and presented for demonstrating the capability of this code.

  6. Physics and 3D in Flash Simulations: Open Source Reality

    NASA Astrophysics Data System (ADS)

    Harold, J. B.; Dusenbery, P.

    2009-12-01

    Over the last decade our ability to deliver simulations over the web has steadily advanced. The improvements in speed of the Adobe Flash engine, and the development of open source tools to expand it, allow us to deliver increasingly sophisticated simulation based games through the browser, with no additional downloads required. In this paper we will present activities we are developing as part of two asteroids education projects: Finding NEO (funded through NSF and NASA SMD), and Asteroids! (funded through NSF). The first activity is Rubble!, an asteroids deflection game built on the open source Box2D physics engine. This game challenges players to push asteroids in to safe orbits before they crash in to the Earth. The Box2D engine allows us to go well beyond simple 2-body orbital calculations and incorporate “rubble piles”. These objects, which are representative of many asteroids, are composed of 50 or more individual rocks which gravitationally bind and separate in realistic ways. Even bombs can be modeled with sufficient physical accuracy to convince players of the hazards of trying to “blow up” incoming asteroids. The ability to easily build games based on underlying physical models allows us to address physical misconceptions in a natural way: by having the player operate in a world that directly collides with those misconceptions. Rubble! provides a particularly compelling example of this due to the variety of well documented misconceptions regarding gravity. The second activity is a Light Curve challenge, which uses the open source PaperVision3D tools to analyze 3D asteroid models. The goal of this activity is to introduce the player to the concept of “light curves”, measurements of asteroid brightness over time which are used to calculate the asteroid’s period. These measurements can even be inverted to generate three dimensional models of asteroids that are otherwise too small and distant to directly image. Through the use of the Paper

  7. Spacecraft charging analysis with the implicit particle-in-cell code iPic3D

    SciTech Connect

    Deca, J.; Lapenta, G.; Marchand, R.; Markidis, S.

    2013-10-15

    We present the first results on the analysis of spacecraft charging with the implicit particle-in-cell code iPic3D, designed for running on massively parallel supercomputers. The numerical algorithm is presented, highlighting the implementation of the electrostatic solver and the immersed boundary algorithm; the latter which creates the possibility to handle complex spacecraft geometries. As a first step in the verification process, a comparison is made between the floating potential obtained with iPic3D and with Orbital Motion Limited theory for a spherical particle in a uniform stationary plasma. Second, the numerical model is verified for a CubeSat benchmark by comparing simulation results with those of PTetra for space environment conditions with increasing levels of complexity. In particular, we consider spacecraft charging from plasma particle collection, photoelectron and secondary electron emission. The influence of a background magnetic field on the floating potential profile near the spacecraft is also considered. Although the numerical approaches in iPic3D and PTetra are rather different, good agreement is found between the two models, raising the level of confidence in both codes to predict and evaluate the complex plasma environment around spacecraft.

  8. Optimizing prostate needle biopsy through 3D simulation

    NASA Astrophysics Data System (ADS)

    Zeng, Jianchao; Kaplan, Charles; Xuan, Jian Hua; Sesterhenn, Isabell A.; Lynch, John H.; Freedman, Matthew T.; Mun, Seong K.

    1998-06-01

    Prostate needle biopsy is used for the detection of prostate cancer. The protocol of needle biopsy that is currently routinely used in the clinical environment is the systematic sextant technique, which defines six symmetric locations on the prostate surface for needle insertion. However, this protocol has been developed based on the long-term observation and experience of urologists. Little quantitative or scientific evidence supports the use of this biopsy technique. In this research, we aim at developing a statistically optimized new prostate needle biopsy protocol to improve the quality of diagnosis of prostate cancer. This new protocol will be developed by using a three-dimensional (3-D) computer- based probability map of prostate cancer. For this purpose, we have developed a computer-based 3-D visualization and simulation system with prostate models constructed from the digitized prostate specimens, in which the process of prostate needle biopsy can be simulated automatically by the computer. In this paper, we first develop an interactive biopsy simulation mode in the system, and evaluate the performance of the automatic biopsy simulation with the sextant biopsy protocol by comparing the results by the urologist using the interactive simulation mode with respect to 53 prostate models. This is required to confirm that the automatic simulation is accurate and reliable enough for the simulation with respect to a large number of prostate models. Then we compare the performance of the existing protocols using the automatic biopsy simulation system with respect to 107 prostate models, which will statistically identify if one protocol is better than another. Since the estimation of tumor volume is extremely important in determining the significance of a tumor and in deciding appropriate treatment methods, we further investigate correlation between the tumor volume and the positive core volume with 89 prostate models. This is done in order to develop a method to

  9. Coherent Vortex Simulations of 3D isotropic turbulence

    NASA Astrophysics Data System (ADS)

    Goldstein, Daniel E.; Vasilyev, Oleg V.; Kevlahan, Nicholas K.-R.

    2006-11-01

    This is the first of three talks on the wavelet filter based dynamically adaptive eddy capturing computational methodology that unifies variable fidelity simulation approaches such as wavelet-based DNS, Coherent Vortex Simulation (CVS), and Stochastic Coherent Adaptive Large Eddy Simulation. The commonality of these approaches is their ability to identify and ``track" on an adaptive mesh energetic coherent vortical structures. In CVS the velocity field is decomposed into two orthogonal parts: a coherent, inhomogeneous, non-Gaussian component and an incoherent, homogeneous, Gaussian component. This separation of coherent and incoherent components is achieved by wavelet thresholding which can be viewed as a non-linear filter that depends on each flow realization. The essence of the CVS approach is to solve for the coherent non-Gaussian component of a turbulent flow field. It has been shown previously that second generation bi-orthogonal wavelet threshold filtering is able to decompose a turbulent velocity field such that the total resulting SGS dissipation is approximately zero. This physically allows a CVS simulation to recover low order statistics with no SGS model. In this work CVS simulations of decaying incompressible 3D isotropic turbulence are compared to DNS results. -6pt

  10. Simulated Photoevaporative Mass Loss from Hot Jupiters in 3D

    NASA Astrophysics Data System (ADS)

    Tripathi, Anjali; Kratter, Kaitlin M.; Murray-Clay, Ruth A.; Krumholz, Mark R.

    2015-08-01

    Ionizing stellar photons heat the upper regions of planetary atmospheres, driving atmospheric mass loss. Gas escaping from several hot, hydrogen-rich planets has been detected using UV and X-ray transmission spectroscopy. Because these planets are tidally locked, and thus asymmetrically irradiated, escaping gas is unlikely to be spherically symmetric. In this paper, we focus on the effects of asymmetric heating on local outflow structure. We use the Athena code for hydrodynamics to produce 3D simulations of hot Jupiter mass loss that jointly model wind launching and stellar heating via photoionization. Our fiducial planet is an inflated, hot Jupiter with radius {R}{{p}}=2.14{R}{Jup} and mass {M}{{p}}=0.53{M}{Jup}. We irradiate the initially neutral, atomic hydrogen atmosphere with 13.6 eV photons and compute the outflow’s ionization structure. There are clear asymmetries in the atmospheric outflow, including a neutral shadow on the planet’s nightside. Given an incident ionizing UV flux comparable to that of the Sun, we find a steady-state mass loss rate of ˜ 2× {10}10 g s-1. The total mass loss rate and the outflow substructure along the substellar ray show good agreement with earlier 1D models, for two different fluxes. Our 3D data cube can be used to generate the outflow’s extinction spectrum during transit. As a proof of concept, we find absorption of stellar Lyα at Doppler-shifted velocities of up to ±50 km s-1. Our work provides a starting point for further 3D models that can be used to predict observable signatures of hot Jupiter mass loss.

  11. A 3D diamond detector for particle tracking

    NASA Astrophysics Data System (ADS)

    Artuso, M.; Bachmair, F.; Bäni, L.; Bartosik, M.; Beacham, J.; Bellini, V.; Belyaev, V.; Bentele, B.; Berdermann, E.; Bergonzo, P.; Bes, A.; Brom, J.-M.; Bruzzi, M.; Cerv, M.; Chau, C.; Chiodini, G.; Chren, D.; Cindro, V.; Claus, G.; Collot, J.; Costa, S.; Cumalat, J.; Dabrowski, A.; D`Alessandro, R.; de Boer, W.; Dehning, B.; Dobos, D.; Dünser, M.; Eremin, V.; Eusebi, R.; Forcolin, G.; Forneris, J.; Frais-Kölbl, H.; Gan, K. K.; Gastal, M.; Goffe, M.; Goldstein, J.; Golubev, A.; Gonella, L.; Gorišek, A.; Graber, L.; Grigoriev, E.; Grosse-Knetter, J.; Gui, B.; Guthoff, M.; Haughton, I.; Hidas, D.; Hits, D.; Hoeferkamp, M.; Hofmann, T.; Hosslet, J.; Hostachy, J.-Y.; Hügging, F.; Jansen, H.; Janssen, J.; Kagan, H.; Kanxheri, K.; Kasieczka, G.; Kass, R.; Kassel, F.; Kis, M.; Kramberger, G.; Kuleshov, S.; Lacoste, A.; Lagomarsino, S.; Lo Giudice, A.; Maazouzi, C.; Mandic, I.; Mathieu, C.; McFadden, N.; McGoldrick, G.; Menichelli, M.; Mikuž, M.; Morozzi, A.; Moss, J.; Mountain, R.; Murphy, S.; Oh, A.; Olivero, P.; Parrini, G.; Passeri, D.; Pauluzzi, M.; Pernegger, H.; Perrino, R.; Picollo, F.; Pomorski, M.; Potenza, R.; Quadt, A.; Re, A.; Riley, G.; Roe, S.; Sapinski, M.; Scaringella, M.; Schnetzer, S.; Schreiner, T.; Sciortino, S.; Scorzoni, A.; Seidel, S.; Servoli, L.; Sfyrla, A.; Shimchuk, G.; Smith, D. S.; Sopko, B.; Sopko, V.; Spagnolo, S.; Spanier, S.; Stenson, K.; Stone, R.; Sutera, C.; Taylor, A.; Traeger, M.; Tromson, D.; Trischuk, W.; Tuve, C.; Uplegger, L.; Velthuis, J.; Venturi, N.; Vittone, E.; Wagner, S.; Wallny, R.; Wang, J. C.; Weilhammer, P.; Weingarten, J.; Weiss, C.; Wengler, T.; Wermes, N.; Yamouni, M.; Zavrtanik, M.

    2016-07-01

    In the present study, results towards the development of a 3D diamond sensor are presented. Conductive channels are produced inside the sensor bulk using a femtosecond laser. This electrode geometry allows full charge collection even for low quality diamond sensors. Results from testbeam show that charge is collected by these electrodes. In order to understand the channel growth parameters, with the goal of producing low resistivity channels, the conductive channels produced with a different laser setup are evaluated by Raman spectroscopy.

  12. Study, simulation and design of a 3D clinostat

    NASA Astrophysics Data System (ADS)

    Pavone, Valentina; Guarnieri, Vincenzo; Lobascio, Cesare; Soma, Aurelio; Bosso, Nicola; Lamantea, Matteo Maria

    High cost and limited number of physically executable experiments in space have introduced the need for ground simulation systems that enable preparing experiments to be carried out on board, identifying phenomena associated with the altered gravity conditions, and taking advantage of these conditions, as in Biotechnology. Among systems developed to simulate microgravity, especially for life sciences experiments, different types of clinostats were realized. This work deals with mechanical design of a three-dimensional clinostat and simulation of the dynamic behavior of the system by varying the operating parameters. The design and simulation phase was preceded by a careful analysis of the state of art and by the review of the most recent results, in particular from the major investigators of Life Sciences in Space. The mechanical design is quite innovative by adoption of a structure entirely in aluminum, which allows robustness while reducing the overall weight. The transmission system of motion has been optimized by means of brushless DC micro motors, light and compact, which helped to reduce weight, dimensions, power consumption and increase the reliability and durability of the system. The study of the dynamic behavior using SIMPACK, a multibody simulation software, led to results in line with those found in the most important and recent scientific publications. This model was also appropriately configured to represent any desired operating condition, and for eventual system scalability. It would be interesting to generate simulated hypogravity - e.g.: 0.38-g (Mars) or 0.17-g (Moon). This would allow to investigate how terrestrial life forms can grow in other planetary habitats, or to determine the gravity threshold response of different organisms. At the moment, such a system can only be achieved by centrifuges in real microgravity. We are confident that simulation and associated tests with our 3D clinostat can help adjusting the parameters allowing variable g

  13. 3D Discrete Element Model with 1 Million Particles: an Example of Hydro-fracturing

    NASA Astrophysics Data System (ADS)

    Liu, C.; Pollard, D. D.

    2013-12-01

    The Discrete Element Method (DEM) permits large relative motion and breakage of elements, and does not require re-meshing, for example as would the Finite Element Method. DEM has a wide range of applications in the fields of solid-earth geophysics, geomechanics, mining engineering, and structural geology. However, due to the computational cost, particle numbers of discrete element models are generally less than a few tens of thousands, which limits the applications. A new 3D DEM system 'MatDEM' can complete dynamic simulations of one million particles. The conversion formulas between particle parameters and model mechanical properties were derived, and the conversion of energy in DEM can be simulated. In a recent paper (Liu et al., 2013, JGR), the analytical solutions of elastic properties and failure modes of a 2D close-packed discrete element model were proposed. Based on these theoretical results, it is easy to create materials using DEM, which have similar mechanical properties to rock. Given the mechanical properties and state of stress, geologists and engineers can investigate the characteristics of rock deformation and failure under different conditions. MatDEM provides an alternative way to study the micro-macro relationships of rock and soil, and the evolution of geologic structures. As an example, MatDEM was used to investigate the generation and development of fluid driven fractures around a micro pore. The simulation result of fractures of an anisotropic 3D model, which includes 1 million particles, is demonstrated. Via parallel computing technology, MatDEM may handle tens of millions of particles in near future. Left: Fluid pressure is applied in the pore to generate fractures. Right: Simulation results (black segments represent fractures).

  14. MPSalsa 3D Simulations of Chemically Reacting Flows

    DOE Data Explorer

    Many important scientific and engineering applications require a detailed analysis of complex systems with coupled fluid flow, thermal energy transfer, mass transfer and nonequilibrium chemical reactions. Currently, computer simulations of these complex reacting flow problems are limited to idealized systems in one or two spatial dimensions when coupled with a detailed, fundamental chemistry model. The goal of our research is to develop, analyze and implement advanced MP numerical algorithms that will allow high resolution 3D simulations with an equal emphasis on fluid flow and chemical kinetics modeling. In our research, we focus on the development of new, fully coupled, implicit solution strategies that are based on robust MP iterative solution methods (copied from http://www.cs.sandia.gov/CRF/MPSalsa/). These simulations are needed for scientific and technical areas such as: combustion research for transportation, atmospheric chemistry modeling for pollution studies, chemically reacting flow models for analysis and control of manufacturing processes, surface catalytic reactors for methane to methanol conversion and chemical vapor deposition (CVD) process modeling for production of advanced semiconductor materials (http://www.cs.sandia.gov/CRF/MPSalsa/).

    This project website provides six QuickTime videos of these simulations, along with a small image gallery and slideshow animations. A list of related publications and conference presentations is also made available.

  15. 3D Simulation Modeling of the Tooth Wear Process

    PubMed Central

    Dai, Ning; Hu, Jian; Liu, Hao

    2015-01-01

    Severe tooth wear is the most common non-caries dental disease, and it can seriously affect oral health. Studying the tooth wear process is time-consuming and difficult, and technological tools are frequently lacking. This paper presents a novel method of digital simulation modeling that represents a new way to study tooth wear. First, a feature extraction algorithm is used to obtain anatomical feature points of the tooth without attrition. Second, after the alignment of non-attrition areas, the initial homogeneous surface is generated by means of the RBF (Radial Basic Function) implicit surface and then deformed to the final homogeneous by the contraction and bounding algorithm. Finally, the method of bilinear interpolation based on Laplacian coordinates between tooth with attrition and without attrition is used to inversely reconstruct the sequence of changes of the 3D tooth morphology during gradual tooth wear process. This method can also be used to generate a process simulation of nonlinear tooth wear by means of fitting an attrition curve to the statistical data of attrition index in a certain region. The effectiveness and efficiency of the attrition simulation algorithm are verified through experimental simulation. PMID:26241942

  16. Three-dimensional parallel UNIPIC-3D code for simulations of high-power microwave devices

    NASA Astrophysics Data System (ADS)

    Wang, Jianguo; Chen, Zaigao; Wang, Yue; Zhang, Dianhui; Liu, Chunliang; Li, Yongdong; Wang, Hongguang; Qiao, Hailiang; Fu, Meiyan; Yuan, Yuan

    2010-07-01

    This paper introduces a self-developed, three-dimensional parallel fully electromagnetic particle simulation code UNIPIC-3D. In this code, the electromagnetic fields are updated using the second-order, finite-difference time-domain method, and the particles are moved using the relativistic Newton-Lorentz force equation. The electromagnetic field and particles are coupled through the current term in Maxwell's equations. Two numerical examples are used to verify the algorithms adopted in this code, numerical results agree well with theoretical ones. This code can be used to simulate the high-power microwave (HPM) devices, such as the relativistic backward wave oscillator, coaxial vircator, and magnetically insulated line oscillator, etc. UNIPIC-3D is written in the object-oriented C++ language and can be run on a variety of platforms including WINDOWS, LINUX, and UNIX. Users can use the graphical user's interface to create the complex geometric structures of the simulated HPM devices, which can be automatically meshed by UNIPIC-3D code. This code has a powerful postprocessor which can display the electric field, magnetic field, current, voltage, power, spectrum, momentum of particles, etc. For the sake of comparison, the results computed by using the two-and-a-half-dimensional UNIPIC code are also provided for the same parameters of HPM devices, the numerical results computed from these two codes agree well with each other.

  17. Colossal Tooling Design: 3D Simulation for Ergonomic Analysis

    NASA Technical Reports Server (NTRS)

    Hunter, Steve L.; Dischinger, Charles; Thomas, Robert E.; Babai, Majid

    2003-01-01

    The application of high-level 3D simulation software to the design phase of colossal mandrel tooling for composite aerospace fuel tanks was accomplished to discover and resolve safety and human engineering problems. The analyses were conducted to determine safety, ergonomic and human engineering aspects of the disassembly process of the fuel tank composite shell mandrel. Three-dimensional graphics high-level software, incorporating various ergonomic analysis algorithms, was utilized to determine if the process was within safety and health boundaries for the workers carrying out these tasks. In addition, the graphical software was extremely helpful in the identification of material handling equipment and devices for the mandrel tooling assembly/disassembly process.

  18. 3D Simulations of Supernovae into the Young Remnant Phase

    NASA Astrophysics Data System (ADS)

    Ellinger, Carola I.; Young, P. A.; Fryer, C.; Rockefeller, G.; Park, S.

    2013-01-01

    The explosion of massive stars as core-collapse supernovae is an inherently three dimensional phenomenon. Observations of the young, ejecta dominated remnants of those explosions unambiguously demonstrate that asymmetry on large and small scales is the rule, rather than the exception. Numerical models of supernova remnants connect the observed remnants to models of the exploding stellar system and thus facilitate both improved interpretations of the observations as well as improve our understanding of the explosion mechanism. We present first 3D simulations of core- collapse supernovae evolving into supernova remnants calculated with SNSPH. The calculations were started from 1D collapsed models of 2 progenitor stars of different types, and follow the explosion from revival of the shock wave to shock break out in 3D. Two different interstellar media, a cold neutral medium and a dense molecular cloud, as well as a red supergiant stellar wind profile, were added to the explosion calculations shortly before shock breakout, so that the blast wave stays in the simulation. With this setup we can follow the dispersal of the nucleosynthesis products from the explosion into the Sedov stage of the supernova remnant evolution starting from realistic initial conditions for the supernova ejecta. We will present a first investigation in the mixing between stellar and interstellar matter as the supernova evolves into the young supernova remnant phase, and contrasts differences that are observed between the scenarios that are investigated. One of the goals is to distinguish between features that arose in instabilities during the explosion from those that were created in the interaction with the surrounding medium.

  19. The Vajont disaster: a 3D numerical simulation for the slide and the waves

    NASA Astrophysics Data System (ADS)

    Rubino, Angelo; Androsov, Alexey; Vacondio, Renato; Zanchettin, Davide; Voltzinger, Naum

    2016-04-01

    A very high resolution O(5 m), 3D hydrostatic nonlinear numerical model was used to simulate the dynamics of both the slide and the surface waves produced during the Vajont disaster (north Italy, 1963), one of the major landslide-induced tsunamis ever documented. Different simulated wave phenomena like, e.g., maximum run-up on the opposite shore, maximum height, and water velocity were analyzed and compared with data available in literature, including the results of a fully 3D simulation obtained with a Smoothed Particle Hydrodynamic code. The difference between measured and simulated after-slide bathymetries was calculated and used in an attempt to quantify the relative magnitude and extension of rigid and fluid motion components during the event.

  20. PAB3D Simulations for the CAWAPI F-16XL

    NASA Technical Reports Server (NTRS)

    Elmiligui, Alaa; Abdol-Hamid, K. S.; Massey, Steven J.

    2007-01-01

    Numerical simulations of the flow around F-16XL are performed as a contribution to the Cranked Arrow Wing Aerodynamic Project International (CAWAPI) using the PAB3D CFD code. Two turbulence models are used in the calculations: a standard k-! model, and the Shih-Zhu-Lumley (SZL) algebraic stress model. Seven flight conditions are simulated for the flow around the F-16XL where the free stream Mach number varies from 0.242 to 0.97. The range of angles of attack varies from 0deg to 20deg. Computational results, surface static pressure, boundary layer velocity profiles, and skin friction are presented and compared with flight data. Numerical results are generally in good agreement with flight data, considering that only one grid resolution is utilized for the different flight conditions simulated in this study. The ASM results are closer to the flight data than the k-! model results. The ASM predicted a stronger primary vortex, however, the origin of the vortex and footprint is approximately the same as in the k-! predictions.

  1. 3D Dynamic Earthquake Fracture Simulation (Test Case)

    NASA Astrophysics Data System (ADS)

    Korkusuz Öztürk, Yasemin; Meral Özel, Nurcan; Ando, Ryosuke

    2016-04-01

    A 3D dynamic earthquake fracture simulation is being developed for the fault structures which are non-planar to understand heterogeneous stress states in the Marmara Sea. Locating in a seismic gap, a large earthquake is expected in the center of the Sea of Marmara. Concerning the fact that more than 14 million inhabitants of İstanbul, located very closely to the Marmara Sea, the importance of the analysis of the Central Marmara Sea is extremely high. A few 3D dynamic earthquake fracture studies have been already done in the Sea of Marmara for pure right lateral strike-slip stress regimes (Oglesby and Mai, 2012; Aochi and Ulrich, 2015). In this study, a 3D dynamic earthquake fracture model with heterogeneous stress patches from the TPV5, a SCEC code validation case, is adapted. In this test model, the fault and the ground surfaces are gridded by a scalene triangulation technique using GMSH program. For a grid size changing between 0.616 km and 1.050 km the number of elements for the fault surface is 1984 and for the ground surface is 1216. When these results are compared with Kaneko's results for TPV5 from SPECFEM3D, reliable findings could be observed for the first 6.5 seconds (stations on the fault) although a stability problem is encountered after this time threshold. To solve this problem grid sizes are made smaller, so the number of elements increase 7986 for the fault surface and 4867 for the ground surface. On the other hand, computational problems arise in that case, since the computation time is directly proportional to the number of total elements and the required memory also increases with the square of that. Therefore, it is expected that this method can be adapted for less coarse grid cases, regarding the main difficulty coming from the necessity of an effective supercomputer and run time limitations. The main objective of this research is to obtain 3D dynamic earthquake rupture scenarios, concerning not only planar and non-planar faults but also

  2. Freestanding 3D supramolecular particle bridges: fabrication and mechanical behavior.

    PubMed

    Ling, Xing Yi; Phang, In Yee; Schönherr, Holger; Reinhoudt, David N; Vancso, G Julius; Huskens, Jurriaan

    2009-06-01

    Freestanding particle bridges with controlled composition and macroscopic robustness are demonstrated by the use of supramolecular nanoparticle assembly. Self-assembly of nanoparticles, templating, and supramolecular glue infiltration are combined to form stable and ordered three-dimensional polystyrene particle composites on a polydimethylsiloxane stamp. Freestanding hybrid polystyrene nanoparticle bridges are obtained by transfer printing of the hybrid structures onto topographically patterned substrates via host-guest interactions. The mechanical robustness and rigidity of the particle bridges can be controlled by manipulating the layer-by-layer cycles of supramolecular glues of gold nanoparticles and dendrimers. Atomic force microscopy-based microbending results, in particular the location and force-dependent deflection behavior, confirm that the particle bridge fulfills the classical supported-beam characteristics. As estimated from classical beam theory, the bending moduli of the particle bridges vary between 0.8 and 1.1 GPa, depending on the degree of filling by the supramolecular glues. Failure analysis on the particle structure indicates linear elastic behavior and a plastic deformation upon failure. PMID:19373830

  3. Simulation of underwater explosion benchmark experiments with ALE3D

    SciTech Connect

    Couch, R.; Faux, D.

    1997-05-19

    Some code improvements have been made during the course of this study. One immediately obvious need was for more flexibility in the constitutive representation for materials in shell elements. To remedy this situation, a model with a tabular representation of stress versus strain and rate dependent effects was implemented. This was required in order to obtain reasonable results in the IED cylinder simulation. Another deficiency was in the ability to extract and plot variables associated with shell elements. The pipe whip analysis required the development of a scheme to tally and plot time dependent shell quantities such as stresses and strains. This capability had previously existed only for solid elements. Work was initiated to provide the same range of plotting capability for structural elements that exist with the DYNA3D/TAURUS tools. One of the characteristics of these problems is the disparity in zoning required in the vicinity of the charge and bubble compared to that needed in the far field. This disparity can cause the equipotential relaxation logic to provide a less than optimal solution. Various approaches were utilized to bias the relaxation to obtain more optimal meshing during relaxation. Extensions of these techniques have been developed to provide more powerful options, but more work still needs to be done. The results presented here are representative of what can be produced with an ALE code structured like ALE3D. They are not necessarily the best results that could have been obtained. More experience in assessing sensitivities to meshing and boundary conditions would be very useful. A number of code deficiencies discovered in the course of this work have been corrected and are available for any future investigations.

  4. 3-D Relativistic MHD Simulations of Extragalactic Jets

    NASA Astrophysics Data System (ADS)

    Nishikawa, K.-I.; Koide, S.; Sakai, J.-I.; Frank, J.; Christodoulou, D. M.; Sol, H.; Mutel, R. L.

    1997-12-01

    We present the numerical simulations of relativistic jets propagating initially oblique to the field lines of a magnetized ambient medium. Our simulations incorporate relativistic MHD in a four-dimensional spacetime and clearly show that (a) relatively weak, oblique fields (at 1/16 of the equipartition value) have only a negligible influence on the propagating jet and they are passively pushed away by the relativistically moving head; (b) oblique fields in equipartition with the ambient plasma provide more resistance and cause bending at the jet head, but the magnitude of this deflection and the associated backflow are small compared to those identified by previous studies with a 2-D slab model. The new results are understood as follows: Relativistic simulations have consistently shown that these jets are effectively heavy and so they do not suffer substantial momentum losses and are not decelerated as efficiently as their nonrelativistic counterparts. In addition, the ambient magnetic field, however strong, can be pushed aside with relative ease by the beam, provided that the degrees of freedom associated with all three spatial dimensions are followed self-consistently during the simulations. The effect is analogous to pushing Japanese ``noren'' or vertical Venetian blinds out of the way while the slats are allowed to bend in 3-D space rather than as a 2-D slab structure. Applied to relativistic extragalactic jets from blazars, the new results are encouraging since superluminal outflows exhibit bending near their sources and their environments are profoundly magnetized---but observations do not provide support for irregular kinematics such as large-scale vortical motions and pronounced reverse flows near the points of origin.

  5. Virtual environment display for a 3D audio room simulation

    NASA Astrophysics Data System (ADS)

    Chapin, William L.; Foster, Scott

    1992-06-01

    Recent developments in virtual 3D audio and synthetic aural environments have produced a complex acoustical room simulation. The acoustical simulation models a room with walls, ceiling, and floor of selected sound reflecting/absorbing characteristics and unlimited independent localizable sound sources. This non-visual acoustic simulation, implemented with 4 audio ConvolvotronsTM by Crystal River Engineering and coupled to the listener with a Poihemus IsotrakTM, tracking the listener's head position and orientation, and stereo headphones returning binaural sound, is quite compelling to most listeners with eyes closed. This immersive effect should be reinforced when properly integrated into a full, multi-sensory virtual environment presentation. This paper discusses the design of an interactive, visual virtual environment, complementing the acoustic model and specified to: 1) allow the listener to freely move about the space, a room of manipulable size, shape, and audio character, while interactively relocating the sound sources; 2) reinforce the listener's feeling of telepresence into the acoustical environment with visual and proprioceptive sensations; 3) enhance the audio with the graphic and interactive components, rather than overwhelm or reduce it; and 4) serve as a research testbed and technology transfer demonstration. The hardware/software design of two demonstration systems, one installed and one portable, are discussed through the development of four iterative configurations. The installed system implements a head-coupled, wide-angle, stereo-optic tracker/viewer and multi-computer simulation control. The portable demonstration system implements a head-mounted wide-angle, stereo-optic display, separate head and pointer electro-magnetic position trackers, a heterogeneous parallel graphics processing system, and object oriented C++ program code.

  6. Online reconstruction of 3D magnetic particle imaging data.

    PubMed

    Knopp, T; Hofmann, M

    2016-06-01

    Magnetic particle imaging is a quantitative functional imaging technique that allows imaging of the spatial distribution of super-paramagnetic iron oxide particles at high temporal resolution. The raw data acquisition can be performed at frame rates of more than 40 volumes s(-1). However, to date image reconstruction is performed in an offline step and thus no direct feedback is available during the experiment. Considering potential interventional applications such direct feedback would be mandatory. In this work, an online reconstruction framework is implemented that allows direct visualization of the particle distribution on the screen of the acquisition computer with a latency of about 2 s. The reconstruction process is adaptive and performs block-averaging in order to optimize the signal quality for a given amount of reconstruction time. PMID:27182668

  7. Online reconstruction of 3D magnetic particle imaging data

    NASA Astrophysics Data System (ADS)

    Knopp, T.; Hofmann, M.

    2016-06-01

    Magnetic particle imaging is a quantitative functional imaging technique that allows imaging of the spatial distribution of super-paramagnetic iron oxide particles at high temporal resolution. The raw data acquisition can be performed at frame rates of more than 40 volumes s‑1. However, to date image reconstruction is performed in an offline step and thus no direct feedback is available during the experiment. Considering potential interventional applications such direct feedback would be mandatory. In this work, an online reconstruction framework is implemented that allows direct visualization of the particle distribution on the screen of the acquisition computer with a latency of about 2 s. The reconstruction process is adaptive and performs block-averaging in order to optimize the signal quality for a given amount of reconstruction time.

  8. 3D thermal modeling of TRISO fuel coupled with neutronic simulation

    SciTech Connect

    Hu, Jianwei; Uddin, Rizwan

    2010-01-01

    The Very High Temperature Gas Reactor (VHTR) is widely considered as one of the top candidates identified in the Next Generation Nuclear Power-plant (NGNP) Technology Roadmap under the U.S . Depanment of Energy's Generation IV program. TRlSO particle is a common element among different VHTR designs and its performance is critical to the safety and reliability of the whole reactor. A TRISO particle experiences complex thermo-mechanical changes during reactor operation in high temperature and high burnup conditions. TRISO fuel performance analysis requires evaluation of these changes on micro scale. Since most of these changes are temperature dependent, 3D thermal modeling of TRISO fuel is a crucial step of the whole analysis package. In this paper, a 3D numerical thermal model was developed to calculate temperature distribution inside TRISO and pebble under different scenarios. 3D simulation is required because pebbles or TRISOs are always subjected to asymmetric thermal conditions since they are randomly packed together. The numerical model was developed using finite difference method and it was benchmarked against ID analytical results and also results reported from literature. Monte-Carlo models were set up to calculate radial power density profile. Complex convective boundary condition was applied on the pebble outer surface. Three reactors were simulated using this model to calculate temperature distribution under different power levels. Two asymmetric boundary conditions were applied to the pebble to test the 3D capabilities. A gas bubble was hypothesized inside the TRISO kernel and 3D simulation was also carried out under this scenario. Intuition-coherent results were obtained and reported in this paper.

  9. Particle Acceleration in the Low Corona Over Broad Longitudes: Coupling Between 3D Magnetohydrodynamic and Energetic Particle Models

    NASA Astrophysics Data System (ADS)

    Gorby, M.; Schwadron, N.; Lee, M. A.; Booth, A. C.; Spence, H.; Torok, T.; Downs, C.; Lionello, R.; Linker, J.; Titov, V. S.; Mikic, Z.; Riley, P.; Desai, M. I.; Dayeh, M. A.; Kozarev, K. A.

    2013-12-01

    Recent work on the coupling between 3D energetic particle models (e.g., the Energetic Particle Radiation Environment Model, EPREM) and magnetohydrodynamic (MHD) models of Coronal Mass Ejections (CMEs, e.g., the PSI MAS model) has demonstrated the efficacy of compression regions in front of fast CMEs for particle acceleration from remarkably low in the corona (3-6 solar radii). Typically particle acceleration becomes rapid beyond 3Rs and often in regions where compression regions have not yet formed active shocks. The challenge for forming large SEP events in such compression-acceleration scenarios is to have enhanced scattering within the acceleration region while also allowing for efficient escape of accelerated particles downstream (away from the Sun) from the compression region. Simulations show rapid particle acceleration in the range of 3-8 Rs over a broad longitudinal region (80°) resulting from the pile-up of magnetic flux in the compression and the subsequent expansion of these fields. These results have important implications for multi-instrument observations that will allow Solar Probe Plus and Solar Orbiter to test the developing paradigm of SEP acceleration and transport from coronal compressions. We present here recent coupled simulations for SEP acceleration and transport, including energetic particle and CME plasma profiles. The broadness of the longitudinal profile from such events may be a key observational test of compression acceleration in the low corona.

  10. 3-D simulation of posterior fossa reduction in Chiari I.

    PubMed

    Fernandes, Yvens Barbosa; Perestrelo, Pedro Fábio Mendonça; Noritomi, Pedro Yoshito; Mathias, Roger Neves; Silva, Jorge Vicente Lopes da; Joaquim, Andrei Fernandes

    2016-05-01

    We proposed a 3D model to evaluate the role of platybasia and clivus length in the development of Chiari I (CI). Using a computer aided design software, two DICOM files of a normal CT scan and MR were used to simulate different clivus lengths (CL) and also different basal angles (BA). The final posterior fossa volume (PFV) was obtained for each variation and the percentage of the volumetric change was acquired with the same method. The initial normal values of CL and BA were 35.65 mm and 112.66º respectively, with a total PFV of 209 ml. Ranging the CL from 34.65 to 29.65 - 24.65 - 19.65, there was a PFV decrease of 0.47% - 1.12% - 1.69%, respectively. Ranging the BA from 122.66º to 127.66º - 142.66º, the PFV decreased 0.69% - 3.23%, respectively. Our model highlights the importance of the basal angle and clivus length to the development of CI. PMID:27191237

  11. 3D finite element simulation of TIG weld pool

    NASA Astrophysics Data System (ADS)

    Kong, X.; Asserin, O.; Gounand, S.; Gilles, P.; Bergheau, J. M.; Medale, M.

    2012-07-01

    The aim of this paper is to propose a three-dimensional weld pool model for the moving gas tungsten arc welding (GTAW) process, in order to understand the main factors that limit the weld quality and improve the productivity, especially with respect to the welding speed. Simulation is a very powerful tool to help in understanding the physical phenomena in the weld process. A 3D finite element model of heat and fluid flow in weld pool considering free surface of the pool and traveling speed has been developed for the GTAW process. Cast3M software is used to compute all the governing equations. The free surface of the weld pool is calculated by minimizing the total surface energy. The combined effects of surface tension gradient, buoyancy force, arc pressure, arc drag force to drive the fluid flow is included in our model. The deformation of the weld pool surface and the welding speed affect fluid flow, heat flow and thus temperature gradients and molten pool dimensions. Welding trials study is presented to compare our numerical results with macrograph of the molten pool.

  12. Simulations of soluble surfactants in 3D multiphase flow

    NASA Astrophysics Data System (ADS)

    Muradoglu, Metin; Tryggvason, Gretar

    2014-10-01

    A finite-difference/front-tracking method is developed for simulations of soluble surfactants in 3D multiphase flows. The interfacial and bulk surfactant concentration evolution equations are solved fully coupled with the incompressible Navier-Stokes equations. A non-linear equation of state is used to relate interfacial surface tension to surfactant concentration at the interface. Simple test cases are designed to validate different parts of the numerical algorithm and the computational results are found to be in a good agreement with the analytical solutions. The numerical algorithm is parallelized using a domain-decomposition method. It is then applied to study the effects of soluble surfactants on the motion of buoyancy-driven bubbles in a straight square channel in nearly undeformable (spherical) and deformable (ellipsoidal) regimes. Finally the method is used to examine the effects of soluble surfactants on the lateral migration of bubbles in a pressure-driven channel flow. It is found that surfactant-induced Marangoni stresses counteract the shear-induced lift force and can reverse the lateral bubble migration completely, i.e., the contaminated bubble drifts away from the channel wall and stabilizes at the center of the channel when the surfactant-induced Marangoni stresses are sufficiently large.

  13. Two new methods for simulating photolithography development in 3D

    SciTech Connect

    Helmsen, J.; Colella, P.; Dorr, M.; Puckett, E.G.

    1997-01-30

    Two methods are presented for simulating the development of photolithographic profiles during the resist dissolution phase. These algorithms are the volume-of-fluid algorithm, and the steady level-set algorithm. They are compared with the ray-trace, cell, and level-set techniques employed in SAMPLE-3D. The volume-of-fluid algorithm employs an Euclidean Grid with volume fractions. At each time step, the surface is reconstructed by computing an approximation of the tangent plane of the surface in each cell that contains a value between 0 and 1. The geometry constructed in this manner is used to determine flow velocity vectors and the flux across each edge. The material is then advanced by a split advection scheme. The steady Level Set algorithm is an extension of the Iterative Level Set algorithm. The steady Level Set algorithm combines Fast Level Set concepts and a technique for finding zero residual solutions to the ( ) function. The etch time for each cell is calculated in a time ordered manner. Use of heap sorting data structures allows the algorithm to execute extremely quickly. Comparisons of the methods have been performed and results shown.

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

    NASA Astrophysics Data System (ADS)

    Schneiter, E. M.; Velázquez, P. F.; Reynoso, E. M.; Esquivel, A.; De Colle, F.

    2015-05-01

    We use three-dimensional magnetohydrodynamic simulations to model the supernova remnant SN 1006. From our numerical results, we have carried out a polarization study, obtaining synthetic maps of the polarized intensity, the Stokes parameter Q, and the polar-referenced angle, which can be compared with observational results. Synthetic maps were computed considering two possible particle acceleration mechanisms: quasi-parallel and quasi-perpendicular. The comparison of synthetic maps of the Stokes parameter Q maps with observations proves to be a valuable tool to discern unambiguously which mechanism is taking place in the remnant of SN 1006, giving strong support to the quasi-parallel model.

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

    NASA Astrophysics Data System (ADS)

    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⊥ << 1, ω/Ωci << 1, ρEP/L << 1) are 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.

  16. 3D Simulations of methane convective storms on Titan's atmosphere

    NASA Astrophysics Data System (ADS)

    Hueso, R.; Sánchez-Lavega, A.

    2005-08-01

    The arrival of the Cassini/Huygens mission to Titan has opened an unprecedented opportunity to study the atmosphere of this satellite. Under the pressure-temperature conditions on Titan, methane, a large atmospheric component amounting perhaps to a 3-5% of the atmosphere, is close to its triple point, potentially playing a similar role as water on Earth. The Huygens probe has shown a terrain shaped by erosion of probably liquid origin, suggestive of past rain. On the other hand, Voyager IRIS spectroscopic observations of Titan imply a saturated atmosphere of methane (amounting perhaps to 150 covered by methane clouds, if we think on Earth meteorology. However, observations from Earth and Cassini have shown that clouds are localized, transient and fast evolving, in particular in the South Pole (currently in its summer season). This might imply a lack of widespread presence on Titan of nuclei where methane could initiate condensation and particle growth with subsequent precipitation. We investigate different scenarios of moist convective storms on Titan using a complete 3D atmospheric model that incorporates a full microphysics treatment required to study cloud formation processes under a saturated atmosphere with low concentration of condensation nuclei. We study local convective development under a variety of atmospheric conditions: sub-saturation, super-saturation, abundances of condensation nuclei fall, condensation nuclei lifted from the ground or gently falling from the stratosphere. We show that under the appropriate circumstances, precipitation rates comparable to typical tropical storms on Earth can be found. Acknowledgements: This work has been funded by Spanish MCYT PNAYA2003-03216, fondos FEDER and Grupos UPV 15946/2004. R. Hueso acknowledges a post-doc fellowship from Gobierno Vasco.

  17. Parallel implementation of 3D FFT with volumetric decomposition schemes for efficient molecular dynamics simulations

    NASA Astrophysics Data System (ADS)

    Jung, Jaewoon; Kobayashi, Chigusa; Imamura, Toshiyuki; Sugita, Yuji

    2016-03-01

    Three-dimensional Fast Fourier Transform (3D FFT) plays an important role in a wide variety of computer simulations and data analyses, including molecular dynamics (MD) simulations. In this study, we develop hybrid (MPI+OpenMP) parallelization schemes of 3D FFT based on two new volumetric decompositions, mainly for the particle mesh Ewald (PME) calculation in MD simulations. In one scheme, (1d_Alltoall), five all-to-all communications in one dimension are carried out, and in the other, (2d_Alltoall), one two-dimensional all-to-all communication is combined with two all-to-all communications in one dimension. 2d_Alltoall is similar to the conventional volumetric decomposition scheme. We performed benchmark tests of 3D FFT for the systems with different grid sizes using a large number of processors on the K computer in RIKEN AICS. The two schemes show comparable performances, and are better than existing 3D FFTs. The performances of 1d_Alltoall and 2d_Alltoall depend on the supercomputer network system and number of processors in each dimension. There is enough leeway for users to optimize performance for their conditions. In the PME method, short-range real-space interactions as well as long-range reciprocal-space interactions are calculated. Our volumetric decomposition schemes are particularly useful when used in conjunction with the recently developed midpoint cell method for short-range interactions, due to the same decompositions of real and reciprocal spaces. The 1d_Alltoall scheme of 3D FFT takes 4.7 ms to simulate one MD cycle for a virus system containing more than 1 million atoms using 32,768 cores on the K computer.

  18. Metos3D: a marine ecosystem toolkit for optimization and simulation in 3-D - Simulation Package v0.2

    NASA Astrophysics Data System (ADS)

    Piwonski, J.; Slawig, T.

    2015-06-01

    A general programming interface for parameter identification for marine ecosystem models is introduced. A comprehensive solver software for periodic steady-states is implemented that includes a fixed point iteration (spin-up) and a Newton solver. The software is based on the Portable, Extensible Toolkit for Scientific Computation (PETSc) library and uses transport matrices for efficient off-line simulation in 3-D. In addition to the usage of PETSc's parallel data structures and PETSc's Newton solver, an own load balancing algorithm is implemented. A simple verification is carried out using a well investigated biogeochemical model for phosphate (PO4) and dissolved organic phosphorous (DOP) with 7 parameters. The model is coupled via the interface to transport matrices that correspond to a longitudinal and latitudinal resolution of 2.8125° and 15 vertical layers. Initial tests show that both solvers and the load balancing algorithm work correctly. Further experiments demonstrate the robustness of the Newton solver with respect to parameter variations. Moreover, the numerical tests reveal that, with optimal control settings, the Newton solver converges at least 6 times faster towards a solution than the spin-up. However, additional twin experiments reveal differences between both solvers regarding a derivative-based black-box optimization. Whereas an optimization run with spin-up-based model evaluations is capable to identify model parameters of a reference solution, Newton-based model evaluations result in an inaccurate gradient approximation.

  19. Towards implementing plate tectonics in 3D mantle convection simulations

    NASA Astrophysics Data System (ADS)

    Bollada, Peter; Davies, Huw

    2010-05-01

    One of the great challenges in numerical mantle convection simulations is to achieve models that naturally develop plate tectonic like behaviour at the surface. In this work we are looking to achieve such models by investigating the set of models where a single consistent rheology is used for the whole model. We have started by investigating a viscoelastic rheology, related to the Oldroyd-B model from the field of polymers. The goal will be to have the parameter that controls the relaxation between elastic and viscous behaviour to depend upon temperature, pressure and strain-rate. With an appropriate choice of this dependence we have, on the near surface, high viscous/elastic regions interfaced with lower, pure viscous, regions of high strain-rate; while it also becomes more viscous at depth in the interior. In this way we hope to obtain plate like behaviour at the surface which naturally progresses to viscous convective behaviour in the interior. We have started to implement this model in the established mantle 3D finite element spherical mantle convection code TERRA (Baumgardner, 1984). Some parts of the model have been implemented as a force (to be combined with the gravitational body force) on the right hand side. The work has required us to develop and code in TERRA: (i) methods to overcome the continuity problem of the stress field stemming from the fact that the velocity field is represented by linear finite elements; (ii) new operators to handle stress and its gradients; (iii) methods to analyse plate-like behaviour at the surface (iv) the necessary functional dependence of viscosity and elastic relaxation time on temperature, strain-rate and pressure We will present the background to the work, its implementation and results.

  20. The ESyS_Particle: A New 3-D Discrete Element Model with Single Particle Rotation

    NASA Astrophysics Data System (ADS)

    Wang, Yucang; Mora, Peter

    In this paper, the Discrete Element Model (DEM) is reviewed, and the ESyS_Particle, our new version of DEM, is introduced. We particularly highlight some of the major physical concerns about DEMs and major differences between our model and most current DEMs. In the new model, single particle rotation is introduced and represented by a unit quaternion. For each 3-D particle, six degrees of freedom are employed: three for translational motion, and three for orientation. Six kinds of relative motions are permitted between two neighboring particles, and six interactions are transferred, i.e., radial, two shearing forces, twisting and two bending torques. The relative rotation between two particles is decomposed into two sequence-independent rotations such that all interactions due to the relative motions between interactive rigid bodies can be uniquely determined. This algorithm can give more accurate results because physical principles are obeyed. A theoretical analysis about how to choose the model parameters is presented. Several numerical tests have been carried out, the results indicate that most laboratory tests can be well reproduced using our model.

  1. 3D Hybrid Simulation of the Titan's Plasma Environment

    NASA Astrophysics Data System (ADS)

    Lipatov, A. S.; Sittler, E. C.; Hartle, R. E.

    2007-05-01

    Titan plays an important role as a simulation laboratory for multiscale kinetic plasma processes which are key processes in space and laboratory plasmas. A development of multiscale combined numerical methods allows us to use more realistic plasma models at Titan. In this report, we describe a Particle-Ion--Fluid-Ion--Fluid- Electron method of kinetic ion-neutral simulation code (see, e.g. [Lipatov, 2002]). This method takes into account charge-exchange and photoionization processes. The model of atmosphere of Titan was based on a paper by Sittler, Hartle, Vinas et al., [2005]. The background ions H+, O+ and pickup ions H2+, CH4+ and N2+ are described in a kinetic approximation, where the electrons are approximated as a fluid. In this report we study the coupling between background ions and pickup ions on the multiple space scales determined by the ion gyroradiis. The first results of such a simulation of the dynamics of ions near Titan are discussed in this report and compared with recent measurements made by the Cassini Plasma Spectrometer (CAPS, [Hartle, Sittler et al., 2006]). E C Sittler Jr., R E Hartle, A F Vinas, R E Johnson, H T Smith and I Mueller-Wodarg, J. Geophys. Res., 110, A09302, 2005. R. E. Hartle, E. C. Sittler, F. M. Neubauer, R. E. Johnson, et al., Planet. Space Sci., 54, 1211, 2006. A S Lipatov, The Hybrid Multiscale Simulation Technology. An Introduction with Application to Astrophysical and Laboratory Plasmas, Springer-Verlag, Berlin, Heidelberg, New York, 2002, p.p. 1-403.

  2. A Software System for Filling Complex Holes in 3D Meshes by Flexible Interacting Particles

    NASA Astrophysics Data System (ADS)

    Yamazaki, Daisuke; Savchenko, Vladimir

    3D meshes generated by acquisition devices such as laser range scanners often contain holes due to occlusion, etc. In practice, these holes are extremely geometrically and topologically complex. We propose a heuristic hole filling technique using particle systems to fill complex holes with arbitrary topology in 3D meshes. Our approach includes the following steps: hole identification, base surface creation, particle distribution, triangulation, and mesh refinement. We demonstrate the functionality of the proposed surface retouching system on synthetic and real data.

  3. Application of 3D hydrodynamic and particle tracking models for better environmental management of finfish culture

    NASA Astrophysics Data System (ADS)

    Moreno Navas, Juan; Telfer, Trevor C.; Ross, Lindsay G.

    2011-04-01

    Hydrographic conditions, and particularly current speeds, have a strong influence on the management of fish cage culture. These hydrodynamic conditions can be used to predict particle movement within the water column and the results used to optimise environmental conditions for effective site selection, setting of environmental quality standards, waste dispersion, and potential disease transfer. To this end, a 3D hydrodynamic model, MOHID, has been coupled to a particle tracking model to study the effects of mean current speed, quiescent water periods and bulk water circulation in Mulroy Bay, Co. Donegal Ireland, an Irish fjard (shallow fjordic system) important to the aquaculture industry. A Lagangrian method simulated the instantaneous release of "particles" emulating discharge from finfish cages to show the behaviour of waste in terms of water circulation and water exchange. The 3D spatial models were used to identify areas of mixed and stratified water using a version of the Simpson-Hunter criteria, and to use this in conjunction with models of current flow for appropriate site selection for salmon aquaculture. The modelled outcomes for stratification were in good agreement with the direct measurements of water column stratification based on observed density profiles. Calculations of the Simpson-Hunter tidal parameter indicated that most of Mulroy Bay was potentially stratified with a well mixed region over the shallow channels where the water is faster flowing. The fjard was characterised by areas of both very low and high mean current speeds, with some areas having long periods of quiescent water. The residual current and the particle tracking animations created through the models revealed an anticlockwise eddy that may influence waste dispersion and potential for disease transfer, among salmon cages and which ensures that the retention time of waste substances from cages is extended. The hydrodynamic model results were incorporated into the ArcView TM GIS

  4. 3D positional tracking of ellipsoidal particles in a microtube flow using holographic microscopy

    NASA Astrophysics Data System (ADS)

    Byeon, Hyeok Jun; Seo, Kyung Won; Lee, Sang Joon

    2014-11-01

    Understanding of micro-scale flow phenomena is getting large attention under advances in micro-scale measurement technologies. Especially, the dynamics of particles suspended in a fluid is essential in both scientific and industrial fields. Moreover, most particles handled in research and industrial fields have non-spherical shapes rather than a simple spherical shape. Under various flow conditions, these non-spherical particles exhibit unique dynamic behaviors. To analyze these dynamic behaviors in a fluid flow, 3D positional information of the particles should be measured accurately. In this study, digital holographic microscopy (DHM) is employed to measure the 3D positional information of non-spherical particles, which are fabricated by stretching spherical polystyrene particles. 3D motions of those particles are obtained by interpreting the holograms captured from particles. Ellipsoidal particles with known size and shape are observed to verify the performance of the DHM technique. In addition, 3D positions of particles in a microtube flow are traced. This DHM technique exhibits promising potential in the analysis of dynamic behaviors of non-spherical particles suspended in micro-scale fluid flows.

  5. Simulation of 3D Seismic Wave Propagation with Volcano Topography

    NASA Astrophysics Data System (ADS)

    Ripperger, J.; Igel, H.; Wassermann, J.

    2001-12-01

    We investigate the possibilities of using three-dimensional finite difference (FD) methods for numerical simulation of the seismic wave field at active volcanoes. We put special emphasis on the implementation of the boundary conditions for free surface topography. We compare two different approaches to solve the free surface boundary conditions. The algorithms are implemented on parallel hardware and have been tested for correctness and stability. We apply them to smooth artificial topographies and to the real topography of Mount Merapi, Indonesia. We conclude, that grid stretching type methods (e.g. Hestholm & Ruud, 1994) are not well suited for realistic volcano topography as they tend to become unstable for large topographic gradients. The representation of topography through staircase shaped grids (Ohminato & Chouet, 1997) results in stable calculations, while demanding very fine gridding. The simulations show the effects of a three-dimensional surface topography on elastic wave propagation. Ground motion at the surface is severely affected by topography. If neglected, this may jeopardize attempts to determine source location by analyzing particle motion. Numerical studies like this can help to understand wave propagation phenomena observed on field recordings in volcano seismology. Future studies will aim at separating the wave effects of internal scattering, topography and sources (tremors, tectonic events, pyroclastic flows).

  6. Measurements of the solid-body rotation of anisotropic particles in 3D turbulence

    NASA Astrophysics Data System (ADS)

    Marcus, Guy G.; Parsa, Shima; Kramel, Stefan; Ni, Rui; Voth, Greg A.

    2014-10-01

    We introduce a new method to measure Lagrangian vorticity and the rotational dynamics of anisotropic particles in a turbulent fluid flow. We use 3D printing technology to fabricate crosses (two perpendicular rods) and jacks (three mutually perpendicular rods). Time-resolved measurements of their orientation and solid-body rotation rate are obtained from four video images of their motion in a turbulent flow between oscillating grids with {{R}λ } = 91. The advected particles have a largest dimension of 6 times the Kolmogorov length, making them a good approximation to anisotropic tracer particles. Crosses rotate like disks and jacks rotate like spheres, so these measurements, combined with previous measurements of tracer rods, allow experimental study of axisymmetric ellipsoids across the full range of aspect ratios. The measured mean square tumbling rate, < {{\\dot{p}}i}{{\\dot{p}}i}> , confirms previous direct numerical simulations that indicate that disks tumble much more rapidly than rods. Measurements of the alignment of a unit vector defining the orientation of crosses with the direction of their solid-body rotation rate vector provide the first direct observation of the alignment of anisotropic particles by the velocity gradients in a turbulent flow.

  7. A Framework for 3D Model-Based Visual Tracking Using a GPU-Accelerated Particle Filter.

    PubMed

    Brown, J A; Capson, D W

    2012-01-01

    A novel framework for acceleration of particle filtering approaches to 3D model-based, markerless visual tracking in monocular video is described. Specifically, we present a methodology for partitioning and mapping the computationally expensive weight-update stage of a particle filter to a graphics processing unit (GPU) to achieve particle- and pixel-level parallelism. Nvidia CUDA and Direct3D are employed to harness the massively parallel computational power of modern GPUs for simulation (3D model rendering) and evaluation (segmentation, feature extraction, and weight calculation) of hundreds of particles at high speeds. The proposed framework addresses the computational intensity that is intrinsic to all particle filter approaches, including those that have been modified to minimize the number of particles required for a particular task. Performance and tracking quality results for rigid object and articulated hand tracking experiments demonstrate markerless, model-based visual tracking on consumer-grade graphics hardware with pixel-level accuracy up to 95 percent at 60+ frames per second. The framework accelerates particle evaluation up to 49 times over a comparable CPU-only implementation, providing an increased particle count while maintaining real-time frame rates. PMID:21301027

  8. Correlated 3D Nanoscale Mapping and Simulation of Coupled Plasmonic Nanoparticles

    PubMed Central

    2015-01-01

    Electron tomography in combination with electron energy-loss spectroscopy (EELS) experiments and simulations was used to unravel the interplay between structure and plasmonic properties of a silver nanocuboid dimer. The precise 3D geometry of the particles fabricated by means of electron beam lithography was reconstructed through electron tomography, and the full three-dimensional information was used as an input for simulations of energy-loss spectra and plasmon resonance maps. Excellent agreement between experiment and theory was found throughout, bringing the comparison between EELS imaging and simulations to a quantitative and correlative level. In addition, interface mode patterns, normally masked by the projection nature of a transmission microscopy investigation, could be unambiguously identified through tomographic reconstruction. This work overcomes the need for geometrical assumptions or symmetry restrictions of the sample in simulations and paves the way for detailed investigations of realistic and complex plasmonic nanostructures. PMID:26495933

  9. 3-D turbulent particle dispersion submodel development. Quarterly progress report No. 1, 5 April--5 July 1991

    SciTech Connect

    Smith, P.J.

    1991-12-31

    The lack of a mathematical description of the interactions of fluid turbulence with other physics-chemical processes is a major obstacle in modeling many industrial program. Turbulent two-phase flow is a phenomenon that is of significant practical importance to coal combustion as well as other disciplines. The interactions of fluid turbulence with the particulate phase has yet to be accurately and efficiently modeled for these industrial applications. On 15 May 1991 work was initiated to cover four major tasks toward the development of a computational submodel for turbulent particle dispersion that would be applicable to coal combustion simulations. Those four tasks are: 1. A critical evaluation of the 2-D Lagrangian particle dispersion submodel, 2. Development of a 3-D submodel for turbulent particle dispersion, 3. Evaluation of the 3-D submodel for turbulent particle dispersion, 4.Exploration of extensions of the Lagrangian dispersion theory to other applications including chemistry-turbulence interactions.

  10. 3D Simulations for a Micron-Scale, Dielectric-Based Acceleration Experiment

    SciTech Connect

    Yoder, R. B.; Travish, G.; Xu Jin; Rosenzweig, J. B.

    2009-01-22

    An experimental program to demonstrate a dielectric, slab-symmetric accelerator structure has been underway for the past two years. These resonant devices are driven by a side-coupled 800-nm laser and can be configured to maintain the field profile necessary for synchronous acceleration and focusing of relativistic or nonrelativistic particles. We present 3D simulations of various versions of the structure geometry, including a metal-walled structure relevant to ongoing cold tests on resonant properties, and an all-dielectric structure to be constructed for a proof-of-principle acceleration experiment.

  11. 3D reconstruction and particle acceleration properties of Coronal Shock Waves During Powerful Solar Particle Events

    NASA Astrophysics Data System (ADS)

    Plotnikov, Illya; Vourlidas, Angelos; Tylka, Allan J.; Pinto, Rui; Rouillard, Alexis; Tirole, Margot

    2016-07-01

    Identifying the physical mechanisms that produce the most energetic particles is a long-standing observational and theoretical challenge in astrophysics. Strong pressure waves have been proposed as efficient accelerators both in the solar and astrophysical contexts via various mechanisms such as diffusive-shock/shock-drift acceleration and betatron effects. In diffusive-shock acceleration, the efficacy of the process relies on shock waves being super-critical or moving several times faster than the characteristic speed of the medium they propagate through (a high Alfven Mach number) and on the orientation of the magnetic field upstream of the shock front. High-cadence, multipoint imaging using the NASA STEREO, SOHO and SDO spacecrafts now permits the 3-D reconstruction of pressure waves formed during the eruption of coronal mass ejections. Using these unprecedented capabilities, some recent studies have provided new insights on the timing and longitudinal extent of solar energetic particles, including the first derivations of the time-dependent 3-dimensional distribution of the expansion speed and Mach numbers of coronal shock waves. We will review these recent developments by focusing on particle events that occurred between 2011 and 2015. These new techniques also provide the opportunity to investigate the enigmatic long-duration gamma ray events.

  12. Multiscale Processes Energizing Plasmas during Reconnection: 3D Simulations in preparation for the MMS mission

    NASA Astrophysics Data System (ADS)

    Lapenta, Giovanni; Goldman, Martin; Newman, David

    2015-11-01

    Magnetic reconnection is a mechanism to convert magnetic energy to particle energy in the form of heat and directed flows. We study here where reconnection and particle energization are found in full 3D models of a reconnecting plasma sheet. Three regions emerge as the primary loci of energy conversion: the separatrices, the dipolarization fronts and the electron diffusion region near x-points. We consider two scenarios: one where the exhaust from multiple x-lines forms a plasmoid (a flux rope in 3D) and one where the exhaust encounters pristine unreconnected plasma and forms a pile-up front. A key process intrinsically 3D, not present in 2D, is the development of an instability in the outflow leading to the formation of secondary reconnection sites that further enhance energy conversion. The MMS mission of NASA was launched on March 12 of this year with the stated goal of finding these regions. We will soon know if we are right in predicting these additional regions of dissipation in the reconnection outflow. Work supported by the NASA MMS Program and by the BOF funds of the KU Leuven. Simulations conducted on NASA Computing facilities, PRACE Tier-0 computing and at NERSC (DOE Office of Science User Contract No. DE-AC02-05CH11231).

  13. Using 3-D Numerical Weather Data in Piloted Simulations

    NASA Technical Reports Server (NTRS)

    Daniels, Taumi S.

    2016-01-01

    This report describes the process of acquiring and using 3-D numerical model weather data sets in NASA Langley's Research Flight Deck (RFD). A set of software tools implement the process and can be used for other purposes as well. Given time and location information of a weather phenomenon of interest, the user can download associated numerical weather model data. These data are created by the National Oceanic and Atmospheric Administration (NOAA) High Resolution Rapid Refresh (HRRR) model, and are then processed using a set of Mathworks' Matlab(TradeMark) scripts to create the usable 3-D weather data sets. Each data set includes radar re ectivity, water vapor, component winds, temperature, supercooled liquid water, turbulence, pressure, altitude, land elevation, relative humidity, and water phases. An open-source data processing program, wgrib2, is available from NOAA online, and is used along with Matlab scripts. These scripts are described with sucient detail to make future modi cations. These software tools have been used to generate 3-D weather data for various RFD experiments.

  14. Method and simulation to study 3D crosstalk perception

    NASA Astrophysics Data System (ADS)

    Khaustova, Dar'ya; Blondé, Laurent; Huynh-Thu, Quan; Vienne, Cyril; Doyen, Didier

    2012-03-01

    To various degrees, all modern 3DTV displays suffer from crosstalk, which can lead to a decrease of both visual quality and visual comfort, and also affect perception of depth. In the absence of a perfect 3D display technology, crosstalk has to be taken into account when studying perception of 3D stereoscopic content. In order to improve 3D presentation systems and understand how to efficiently eliminate crosstalk, it is necessary to understand its impact on human perception. In this paper, we present a practical method to study the perception of crosstalk. The approach consists of four steps: (1) physical measurements of a 3DTV, (2) building of a crosstalk surface based on those measurements and representing specifically the behavior of that 3TV, (3) manipulation of the crosstalk function and application on reference images to produce test images degraded by crosstalk in various ways, and (4) psychophysical tests. Our approach allows both a realistic representation of the behavior of a 3DTV and the easy manipulation of its resulting crosstalk in order to conduct psycho-visual experiments. Our approach can be used in all studies requiring the understanding of how crosstalk affects perception of stereoscopic content and how it can be corrected efficiently.

  15. 3D quantification of brain microvessels exposed to heavy particle radiation

    NASA Astrophysics Data System (ADS)

    Hintermüller, C.; Coats, J. S.; Obenaus, A.; Nelson, G.; Krucker, T.; Stampanoni, M.

    2009-09-01

    Space radiation with high energy particles and cosmic rays presents a significant hazard to spaceflight crews. Recent reviews of the health risk to astronauts from ionizing radiation concluded to establish a level of risk which may indicate the possible performance decrements and decreased latency of late dysfunction syndromes (LDS) of the brain. A hierarchical imaging approach developed at ETH Zürich and PSI, which relies on synchrotron based X-ray Tomographic Microscopy (SRXTM), was used to visualize and analyze 3D vascular structures down to the capillary level in their precise anatomical context. Various morphological parameters, such as overall vessel volume, vessel thickness and spacing, are extracted to characterize the vascular structure within a region of interest. For a first quantification of the effect of high energy particles on the vasculature we scanned a set of 6 animals, all of same age. The animals were irradiated with 1 Gy, 2 Gy and 4 Gy of 600MeV 56Fe heavy particles simulating the space radiation environment. We found that with increasing dose the diameter of vessels and the overall vessel volume are decreased whereas the vessel spacing is increased. As these parameters reflect blood flow in three-dimensional space they can be used as indicators for the degree of vascular efficiency which can have an impact on the function and development of lung tissue or tumors.

  16. Regional airflow and particle distribution in the lung with a 3D-1D coupled subject-specific boundary condition

    NASA Astrophysics Data System (ADS)

    Choi, Jiwoong; Yin, Youbing; Hoffman, Eric; Tawhai, Merryn; Lin, Ching-Long

    2010-11-01

    Correct prediction of regional distribution of inhaled aerosol particles is vital to improve pulmonary medicine. Physiologically consistent regional ventilations of airflow and aerosol particles are simulated with a 3D-1D coupled subject-specific boundary condition (BC). In 3D CT-resolved 7-generation airways, large eddy simulations are performed to capture detailed airflow characteristics and Lagrangian particle simulations are carried to track the particle transport and deposition. Results are compared with two traditional outlet BCs: uniform velocity and uniform pressure. Proposed BC is eligible for physiologically consistent airflow distribution in the lung, while the others are not. The regional ventilation and deposition of particles reflect the regional ventilation of airflow. In this study, two traditional BCs yield up to 98% (334%) over-prediction in lobar particle ventilation (deposition) fraction. Upper to lower particle ventilation ratios of both left and right lungs read ˜0.4 with the proposed BC, while those for the other two BCs vary with the error up to 73%.

  17. Characterizing heterogeneity among virus particles by stochastic 3D signal reconstruction

    NASA Astrophysics Data System (ADS)

    Xu, Nan; Gong, Yunye; Wang, Qiu; Zheng, Yili; Doerschuk, Peter C.

    2015-09-01

    In single-particle cryo electron microscopy, many electron microscope images each of a single instance of a biological particle such as a virus or a ribosome are measured and the 3-D electron scattering intensity of the particle is reconstructed by computation. Because each instance of the particle is imaged separately, it should be possible to characterize the heterogeneity of the different instances of the particle as well as a nominal reconstruction of the particle. In this paper, such an algorithm is described and demonstrated on the bacteriophage Hong Kong 97. The algorithm is a statistical maximum likelihood estimator computed by an expectation maximization algorithm implemented in Matlab software.

  18. 3-D Simulations of Plasma Wakefield Acceleration with Non-Idealized Plasmas and Beams

    SciTech Connect

    Deng, S.; Katsouleas, T.; Lee, S.; Muggli, P.; Mori, W.B.; Hemker, R.; Ren, C.; Huang, C.; Dodd, E.; Blue, B.E.; Clayton, C.E.; Joshi, C.; Wang, S.; Decker, F.J.; Hogan, M.J.; Iverson, R.H.; O'Connell, C.; Raimondi, P.; Walz, D.; /SLAC

    2005-09-27

    3-D Particle-in-cell OSIRIS simulations of the current E-162 Plasma Wakefield Accelerator Experiment are presented in which a number of non-ideal conditions are modeled simultaneously. These include tilts on the beam in both planes, asymmetric beam emittance, beam energy spread and plasma inhomogeneities both longitudinally and transverse to the beam axis. The relative importance of the non-ideal conditions is discussed and a worst case estimate of the effect of these on energy gain is obtained. The simulation output is then propagated through the downstream optics, drift spaces and apertures leading to the experimental diagnostics to provide insight into the differences between actual beam conditions and what is measured. The work represents a milestone in the level of detail of simulation comparisons to plasma experiments.

  19. Methods for Measuring the Orientation and Rotation Rate of 3D-printed Particles in Turbulence.

    PubMed

    Cole, Brendan C; Marcus, Guy G; Parsa, Shima; Kramel, Stefan; Ni, Rui; Voth, Greg A

    2016-01-01

    Experimental methods are presented for measuring the rotational and translational motion of anisotropic particles in turbulent fluid flows. 3D printing technology is used to fabricate particles with slender arms connected at a common center. Shapes explored are crosses (two perpendicular rods), jacks (three perpendicular rods), triads (three rods in triangular planar symmetry), and tetrads (four arms in tetrahedral symmetry). Methods for producing on the order of 10,000 fluorescently dyed particles are described. Time-resolved measurements of their orientation and solid-body rotation rate are obtained from four synchronized videos of their motion in a turbulent flow between oscillating grids with Rλ = 91. In this relatively low-Reynolds number flow, the advected particles are small enough that they approximate ellipsoidal tracer particles. We present results of time-resolved 3D trajectories of position and orientation of the particles as well as measurements of their rotation rates. PMID:27404898

  20. Development of a 3D Digital Particle Image Thermometry and Velocimetry (3DDPITV) System

    NASA Astrophysics Data System (ADS)

    Schmitt, David; Rixon, Greg; Dabiri, Dana

    2006-11-01

    A novel 3D Digital Particle Image Thermometry and Velocimetry (3DDPITV) system has been designed and fabricated. By combining 3D Digital Particle Image Velocimetry (3DDPIV) and Digital Particle Image Thermometry (DPIT) into one system, this technique provides simultaneous temperature and velocity data in a volume of ˜1x1x0.5 in^3 using temperature sensitive liquid crystal particles as flow sensors. Two high-intensity xenon flashlamps were used as illumination sources. The imaging system consists of six CCD cameras, three allocated for measuring velocity, based on particle motion, and three for measuring temperature, based on particle color. The cameras were optically aligned using a precision grid and high-resolution translation stages. Temperature calibration was then performed using a precision thermometer and a temperature-controlled bath. Results from proof-of-concept experiments will be presented and discussed.

  1. Global 3-D Hybrid Simulations of Mars and the Effect of Crustal Magnetic Fields

    NASA Astrophysics Data System (ADS)

    Brecht, S. H.; Ledvina, S. A.; Luhmann, J. G.

    2001-12-01

    Mars is not protected from the solar wind by a strong intrinsic magnetic field. Further, Mars is a very small planet. These two aspects of Mars mean that the solar wind interaction with the planet is both direct and kinetic in behavior. The large gyroradius of the incoming solar wind and the large gyroradius of the pick up ions make the solar wind interaction with Mars very unique. Over the years the authors using a 3-D kinetic hybrid particle code have studied this planet. Mars Global Surveyor has produced many exciting discoveries. Probably the most significant and unexpected was the presence of the crustal magnetic fields on the surface of the planet. The presence of these strong crustal fields suggests that the solar wind interaction with Mars will be altered. In fact, it suggests that the loss rate of pick up ions and the shape of the Martian magnetosphere will be substantially changed. Results will be presented from our latest 3-D hybrid particle code simulations of Mars where models for the crustal magnetic field have been inserted into the simulations. Particular emphasis will be placed on changes to the magnetic field topology and the change in the rate of pick up ions.

  2. 3-D SIMULATIONS OF AIRWAYS WITHIN HUMAN LUNGS

    EPA Science Inventory

    Information regarding the deposition patterns of inhaled particles has important application to the fields of toxicology and medicine. The former concerns the risk assessment of inhaled air pollutants (inhalation toxicology); the latter concerns the targeted delivery of inhaled ...

  3. Enhancing Pre-Service Teachers' Awareness to Pupils' Test-Anxiety with 3D Immersive Simulation

    ERIC Educational Resources Information Center

    Passig, David; Moshe, Ronit

    2008-01-01

    This study investigated whether participating in a 3D immersive virtual reality world simulating the experience of test-anxiety would affect preservice teachers' awareness to the phenomenon. Ninety subjects participated in this study, and were divided into three groups. The experimental group experienced a 3D immersive simulation which made…

  4. Particle-based optical pressure sensors for 3D pressure mapping.

    PubMed

    Banerjee, Niladri; Xie, Yan; Chalaseni, Sandeep; Mastrangelo, Carlos H

    2015-10-01

    This paper presents particle-based optical pressure sensors for in-flow pressure sensing, especially for microfluidic environments. Three generations of pressure sensitive particles have been developed- flat planar particles, particles with integrated retroreflectors and spherical microballoon particles. The first two versions suffer from pressure measurement dependence on particles orientation in 3D space and angle of interrogation. The third generation of microspherical particles with spherical symmetry solves these problems making particle-based manometry in microfluidic environment a viable and efficient methodology. Static and dynamic pressure measurements have been performed in liquid medium for long periods of time in a pressure range of atmospheric to 40 psi. Spherical particles with radius of 12 μm and balloon-wall thickness of 0.5 μm are effective for more than 5 h in this pressure range with an error of less than 5%. PMID:26342493

  5. Full-Scale 3D Simulation of a sputtering magnetron

    NASA Astrophysics Data System (ADS)

    Walton, C. C.; Wilks, S. C.; Ayyaswamy, V.; Verboncoeur, J. P.; Parks, P. B.; Wu, W.; Zhou, C. D.; Stoltz, P. H.

    2010-11-01

    PIC simulations have been used to study ion energy distributions in magnetron plasmas, and coupled with other simulations to relate plasma processes to properties of sputtered films. The plasma is weakly ionized and exchanges heat with the background gas by scattering and charge-exchange reactions. Resulting heating of neutral background gas up to ˜1200K, leading to ˜5X rarefaction and increased plasma impedance, was studied with coupled PIC and Direct Simulation Monte Carlo (DSMC) simulations. Effects of scaling the PIC simulations from 0.1X to 1X physical size, and modifying the plasma potential by a dc substrate bias, will be presented. Comparison to experimental I-V relations and importance for roughness and density of sputtered films will be discussed.

  6. 3D kinetic simulations of the global interaction between the solar wind and the magnetosphere

    NASA Astrophysics Data System (ADS)

    Amaya, Jorge; Maneva, Yana; Deca, Jan; Lapenta, Giovanni

    2015-04-01

    We performed three dimensional simulations of the interaction between the solar wind and the magnetosphere, using the self-consistent fully kinetic code iPic3D. The main objective of our simulations is to link the global interaction phenomena to the local turbulence and reconnection processes in the magnetosphere. Other numerical approaches have been used before to study this problem, including MHD, hybrid and Vlasov codes. However, only particle-in-cell codes offer the possibility to study the kinetic effects of the diffusion regions of the Earth environment that drive the energy transfer from the solar wind to the magnetosphere. Previous attempts to perform such kinds of simulations were limited to unphysical thermal velocities of the ion and electron species, small simulation boxes or cell sizes that do not capture the local kinetic effects at the magnetopause. Using the implicit moment Particle-in-Cell approach we performed simulations that can capture these small scale effects and, at the same time, allow to study large scale phenomena such as the bow shock and the development of the magnetotail. We expect that these results will be used to maximize the impact of future space missions, such as THOR, MMS and BepiColombo, by improving our understanding of the planetary environment, from the conditions observed in the solar wind to the turbulence and reconnection processes downstream of the bow shock.

  7. 3D CFDTD PIC Simulation Study on Low-Frequency Oscillations in a Gyrotron

    NASA Astrophysics Data System (ADS)

    Lin, M. C.; Smithe, D. N.

    2011-10-01

    Low-frequency oscillations (LFOs) have been observed in a high average power gyrotron and the trapped electron population contributing to the oscillation has been measured. As high average power gyrotrons are the most promising millimeter wave source for thermonuclear fusion research, it is important to get a better understanding of this parasitic phenomenon to avoid any deterioration of the electron beam quality thus reducing the gyrotron efficiency. However, understanding of the LFOs remains incomplete and a full picture of this parasitic phenomenon has not been seen yet. In this work, we use a 3D conformal finite-difference time-domain (CFDTD) particle-in-cell (PIC) method to accurately and efficiently study the LFOs in a magnetron injection gun (MIG) of a high average power gyrotron. Employing a highly parallelized computation, the model can be simulated in time domain more realistically. LFOs have been obtained in a 3D time domain simulation for the first time. From our preliminary simulation studies, it is found that not only magnetic compression profile but initial velocity or velocity ratio play an important role in the operation of a MIG electron gun. In addition, the secondary emission effects on the LFOs are also studied. Detailed results will be presented. Work supported by the U.S. Department of Energy under Grant No. DE-SC0004436.

  8. Simple 3-D stimulus for motion parallax and its simulation.

    PubMed

    Ono, Hiroshi; Chornenkyy, Yevgen; D'Amour, Sarah

    2013-01-01

    Simulation of a given stimulus situation should produce the same perception as the original. Rogers et al (2009 Perception 38 907-911) simulated Wheeler's (1982, PhD thesis, Rutgers University, NJ) motion parallax stimulus and obtained quite different perceptions. Wheeler's observers were unable to reliably report the correct direction of depth, whereas Rogers's were. With three experiments we explored the possible reasons for the discrepancy. Our results suggest that Rogers was able to see depth from the simulation partly due to his experience seeing depth with random dot surfaces. PMID:23964382

  9. Tensor3D: A computer graphics program to simulate 3D real-time deformation and visualization of geometric bodies

    NASA Astrophysics Data System (ADS)

    Pallozzi Lavorante, Luca; Dirk Ebert, Hans

    2008-07-01

    Tensor3D is a geometric modeling program with the capacity to simulate and visualize in real-time the deformation, specified through a tensor matrix and applied to triangulated models representing geological bodies. 3D visualization allows the study of deformational processes that are traditionally conducted in 2D, such as simple and pure shears. Besides geometric objects that are immediately available in the program window, the program can read other models from disk, thus being able to import objects created with different open-source or proprietary programs. A strain ellipsoid and a bounding box are simultaneously shown and instantly deformed with the main object. The principal axes of strain are visualized as well to provide graphical information about the orientation of the tensor's normal components. The deformed models can also be saved, retrieved later and deformed again, in order to study different steps of progressive strain, or to make this data available to other programs. The shape of stress ellipsoids and the corresponding Mohr circles defined by any stress tensor can also be represented. The application was written using the Visualization ToolKit, a powerful scientific visualization library in the public domain. This development choice, allied to the use of the Tcl/Tk programming language, which is independent on the host computational platform, makes the program a useful tool for the study of geometric deformations directly in three dimensions in teaching as well as research activities.

  10. 3-D Hybrid Simulation of Interaction Between Solar Wind Discontinuity and Magnetosphere

    NASA Astrophysics Data System (ADS)

    Pang, Y.; Lin, Y.; Deng, X.; Wang, X.

    2008-12-01

    Previous simulations and observations indicate that interaction of interplanetary directional tangential discontinuities (TD), across which only the direction of magnetic field changes,with the bow shock may initiate magnetic reconnection as the TD is transmitted into the magnetosheath. We employ a 3-D electromagnetic, hybrid simulation to study such interaction between the TD and the bow shock-magnetosheath-magnetopause system, while the hybrid model treats the ions kinetically via particle-in-cell methods and the electrons as a massless fluid. We present results corresponding to different orientations of the initial IMF, TDs with various thicknesses ranging from 10 to 30 ion skin depths, and polarizations of magnetic field and senses of field rotation across the TD. Our results indicate that the reconnection rate and structure can be influenced by the width and the structure of the TD. The kinetic structure and evolution of FTEs produced by the magnetosheath reconnection, as they propagate to the magnetopause, will be studied.

  11. GATOR: A 3-D time-dependent simulation code for helix TWTs

    SciTech Connect

    Zaidman, E.G.; Freund, H.P.

    1996-12-31

    A 3D nonlinear analysis of helix TWTs is presented. The analysis and simulation code is based upon a spectral decomposition using the vacuum sheath helix modes. The field equations are integrated on a grid and advanced in time using a MacCormack predictor-corrector scheme, and the electron orbit equations are integrated using a fourth order Runge-Kutta algorithm. Charge is accumulated on the grid and the field is interpolated to the particle location by a linear map. The effect of dielectric liners on the vacuum sheath helix dispersion is included in the analysis. Several numerical cases are considered. Simulation of the injection of a DC beam and a signal at a single frequency is compared with a linear field theory of the helix TWT interaction, and good agreement is found.

  12. 3D Hot Test Simulations of a 220 GHz Folded Waveguide Traveling Wave Tube Using a CFDTD PIC Method

    NASA Astrophysics Data System (ADS)

    Lin, Ming-Chieh; Song, Heather

    2015-11-01

    Millimeter or sub-THz wave sources centered at 220 GHz is of interest due to the potential for its commercial and military applications including high resolution radar, remote sensing, and high-data-rate communications. It has been demonstrated via 3D cold test finite element method (FEM) simulations that a folded waveguide traveling wave tube (FWTWT) can be designed and optimized at this frequency range with a small signal gain of 18 dB over a comparatively broad (-3 dB) bandwidth of ~ 10%. On the other hand, 3D hot test simulations of a V-band ladder TWT have been successfully demonstrated using a conformal finite-difference time-domain (CFDTD) particle-in-cell (PIC) method for center frequency of 50 GHz. In the present work, the 220 GHz FWTWT designs have been reviewed and studied. 3D Cold test simulations using both the CFDTD and FEM methods have been carried out and compared with each other as basis for 3D hot test CFDTD PIC simulations. The preliminary simulation result shows that the gain-bandwidth features at 220 GHz are achievable while carefully avoiding beam interceptions. Our study shows that the interaction characteristics are very sensitive to the operating beam parameters. Detail simulation results and discussions will be presented.

  13. 3-D dynamic rupture simulations by a finite volume method

    NASA Astrophysics Data System (ADS)

    Benjemaa, M.; Glinsky-Olivier, N.; Cruz-Atienza, V. M.; Virieux, J.

    2009-07-01

    Dynamic rupture of a 3-D spontaneous crack of arbitrary shape is investigated using a finite volume (FV) approach. The full domain is decomposed in tetrahedra whereas the surface, on which the rupture takes place, is discretized with triangles that are faces of tetrahedra. First of all, the elastodynamic equations are described into a pseudo-conservative form for an easy application of the FV discretization. Explicit boundary conditions are given using criteria based on the conservation of discrete energy through the crack surface. Using a stress-threshold criterion, these conditions specify fluxes through those triangles that have suffered rupture. On these broken surfaces, stress follows a linear slip-weakening law, although other friction laws can be implemented. For The Problem Version 3 of the dynamic-rupture code verification exercise conducted by the SCEC/USGS, numerical solutions on a planar fault exhibit a very high convergence rate and are in good agreement with the reference one provided by a finite difference (FD) technique. For a non-planar fault of parabolic shape, numerical solutions agree satisfactorily well with those obtained with a semi-analytical boundary integral method in terms of shear stress amplitudes, stopping phases arrival times and stress overshoots. Differences between solutions are attributed to the low-order interpolation of the FV approach, whose results are particularly sensitive to the mesh regularity (structured/unstructured). We expect this method, which is well adapted for multiprocessor parallel computing, to be competitive with others for solving large scale dynamic ruptures scenarios of seismic sources in the near future.

  14. SIERRA - A 3-D device simulator for reliability modeling

    NASA Astrophysics Data System (ADS)

    Chern, Jue-Hsien; Arledge, Lawrence A., Jr.; Yang, Ping; Maeda, John T.

    1989-05-01

    SIERRA is a three-dimensional general-purpose semiconductor-device simulation program which serves as a foundation for investigating integrated-circuit (IC) device and reliability issues. This program solves the Poisson and continuity equations in silicon under dc, transient, and small-signal conditions. Executing on a vector/parallel minisupercomputer, SIERRA utilizes a matrix solver which uses an incomplete LU (ILU) preconditioned conjugate gradient square (CGS, BCG) method. The ILU-CGS method provides a good compromise between memory size and convergence rate. The authors have observed a 5x to 7x speedup over standard direct methods in simulations of transient problems containing highly coupled Poisson and continuity equations such as those found in reliability-oriented simulations. The application of SIERRA to parasitic CMOS latchup and dynamic random-access memory single-event-upset studies is described.

  15. Vectors in Use in a 3D Juggling Game Simulation

    ERIC Educational Resources Information Center

    Kynigos, Chronis; Latsi, Maria

    2006-01-01

    The new representations enabled by the educational computer game the "Juggler" can place vectors in a central role both for controlling and measuring the behaviours of objects in a virtual environment simulating motion in three-dimensional spaces. The mathematical meanings constructed by 13 year-old students in relation to vectors as objects, as a…

  16. Hierarchical Approach to 'Atomistic' 3-D MOSFET Simulation

    NASA Technical Reports Server (NTRS)

    Asenov, Asen; Brown, Andrew R.; Davies, John H.; Saini, Subhash

    1999-01-01

    We present a hierarchical approach to the 'atomistic' simulation of aggressively scaled sub-0.1 micron MOSFET's. These devices are so small that their characteristics depend on the precise location of dopant atoms within them, not just on their average density. A full-scale three-dimensional drift-diffusion atomistic simulation approach is first described and used to verify more economical, but restricted, options. To reduce processor time and memory requirements at high drain voltage, we have developed a self-consistent option based on a solution of the current continuity equation restricted to a thin slab of the channel. This is coupled to the solution of the Poisson equation in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison to the full self-consistent solution. At low drain voltage, a single solution of the nonlinear Poisson equation is sufficient to extract the current with satisfactory accuracy. In this case, the current is calculated by solving the current continuity equation in a drift approximation only, also in a thin slab containing the MOSFET channel. The regions of applicability for the different components of this hierarchical approach are illustrated in example simulations covering the random dopant-induced threshold voltage fluctuations, threshold voltage lowering, threshold voltage asymmetry, and drain current fluctuations.

  17. Advanced 3D Photocathode Modeling and Simulations Final Report

    SciTech Connect

    Dimitre A Dimitrov; David L Bruhwiler

    2005-06-06

    High brightness electron beams required by the proposed Next Linear Collider demand strong advances in photocathode electron gun performance. Significant improvement in the production of such beams with rf photocathode electron guns is hampered by the lack high-fidelity simulations. The critical missing piece in existing gun codes is a physics-based, detailed treatment of the very complex and highly nonlinear photoemission process.

  18. 3D Numerical simulations of the C-start of a Bluegill Sunfish

    NASA Astrophysics Data System (ADS)

    Narayanan, Venkat R. T.; Borazjani, Iman; Sotiropoulos, Fotis; Tytell, Eric D.; Lauder, George V.

    2009-11-01

    Obtaining the 3D flow field, forces, and power produced during the fast start maneuvers of fish is essential for studying this behavior from the hydrodynamics perspective. During a typical fast start, which is typically referred to as the C-start, the fish initially bends its body in a C shape manner and then with a fast stroke bends out of the C shape. We carry out high-resolution, 3D simulations of a bluegill sunfish performing a C-start maneuver. The body geometry and motion during the C-start are obtained from the experimental. We used high-speed video and particle image velocimetry to quantify body motion and flows produced during the C-start. We carry out simulations both with the entire motion prescribed and by prescribing only the deformation of the body but predicting the motion of the fish center of mass via a fluid-structure interaction approach. The computed results are compared with experimental observations and analyzed to further elucidate dynamics and three-dimensional structure of the C-start flowfields.

  19. Large-scale 3D simulations of ICF and HEDP targets

    NASA Astrophysics Data System (ADS)

    Marinak, Michael M.

    2000-10-01

    The radiation hydrodynamics code HYDRA continues to be developed and applied to 3D simulations of a variety of targets for both inertial confinement fusion (ICF) and high energy density physics. Several packages have been added enabling this code to perform ICF target simulations with similar accuracy as two-dimensional codes of long-time historical use. These include a laser ray trace and deposition package, a heavy ion deposition package, implicit Monte Carlo photonics, and non-LTE opacities, derived from XSN or the linearized response matrix approach.(R. More, T. Kato, Phys. Rev. Lett. 81, 814 (1998), S. Libby, F. Graziani, R. More, T. Kato, Proceedings of the 13th International Conference on Laser Interactions and Related Plasma Phenomena, (AIP, New York, 1997).) LTE opacities can also be calculated for arbitrary mixtures online by combining tabular values generated by different opacity codes. Thermonuclear burn, charged particle transport, neutron energy deposition, electron-ion coupling and conduction, and multigroup radiation diffusion packages are also installed. HYDRA can employ ALE hydrodynamics; a number of grid motion algorithms are available. Multi-material flows are resolved using material interface reconstruction. Results from large-scale simulations run on up to 1680 processors, using a combination of massively parallel processing and symmetric multiprocessing, will be described. A large solid angle simulation of Rayleigh-Taylor instability growth in a NIF ignition capsule has resolved simultaneously the full spectrum of the most dangerous modes that grow from surface roughness. Simulations of a NIF hohlraum illuminated with the initial 96 beam configuration have also been performed. The effect of the hohlraum’s 3D intrinsic drive asymmetry on the capsule implosion will be considered. We will also discuss results from a Nova experiment in which a copper sphere is crushed by a planar shock. Several interacting hydrodynamic instabilities, including

  20. 3D imaging of particle tracks in Solid State Nuclear Track Detectors

    NASA Astrophysics Data System (ADS)

    Wertheim, D.; Gillmore, G.; Brown, L.; Petford, N.

    2009-04-01

    Inhalation of radon gas (222Rn) and associated ionizing decay products is known to cause lung cancer in human. In the U.K., it has been suggested that 3 to 5 % of total lung cancer deaths can be linked to elevated radon concentrations in the home and/or workplace. Radon monitoring in buildings is therefore routinely undertaken in areas of known risk. Indeed, some organisations such as the Radon Council in the UK and the Environmental Protection Agency in the USA, advocate a ‘to test is best' policy. Radon gas occurs naturally, emanating from the decay of 238U in rock and soils. Its concentration can be measured using CR?39 plastic detectors which conventionally are assessed by 2D image analysis of the surface; however there can be some variation in outcomes / readings even in closely spaced detectors. A number of radon measurement methods are currently in use (for examples, activated carbon and electrets) but the most widely used are CR?39 solid state nuclear track?etch detectors (SSNTDs). In this technique, heavily ionizing alpha particles leave tracks in the form of radiation damage (via interaction between alpha particles and the atoms making up the CR?39 polymer). 3D imaging of the tracks has the potential to provide information relating to angle and energy of alpha particles but this could be time consuming. Here we describe a new method for rapid high resolution 3D imaging of SSNTDs. A ‘LEXT' OLS3100 confocal laser scanning microscope was used in confocal mode to successfully obtain 3D image data on four CR?39 plastic detectors. 3D visualisation and image analysis enabled characterisation of track features. This method may provide a means of rapid and detailed 3D analysis of SSNTDs. Keywords: Radon; SSNTDs; confocal laser scanning microscope; 3D imaging; LEXT

  1. Photorealistic 3D omni-directional stereo simulator

    NASA Astrophysics Data System (ADS)

    Reiners, Dirk; Cruz-Neira, Carolina; Neumann, Carsten

    2015-03-01

    While a lot of areas in VR have made significant advances, visual rendering in VR is often not quite keeping up with the state of the art. There are many reasons for this, but one way to alleviate some of the issues is by using ray tracing instead of rasterization for image generation. Contrary to popular belief, ray tracing is a realistic, competitive technology nowadays. This paper looks at the pros and cons of using ray tracing and demonstrates the feasibility of employing it using the example of a helicopter flight simulator image generator.

  2. Linking continuum mechanics and 3D discrete dislocation simulations

    SciTech Connect

    El-Azab, A. A.; Fivel, M.

    1998-10-18

    A technique is developed for linking the methods of discrete dislocation dynamics simulation and finite element to treat elasto-plasticity problems. The overall formulation views the plastically deforming crystal as an elastic crystal with continuously changing dislocation microstructure which is tracked by the numerical dynamics simulation. The FEM code needed in this regard is based on linear elasticity only. This formulation presented here is focused on a continuous updating of the outer shape of the crystal, for possible regeneration of the FEM mesh, and adjustment of the surface geometry, in particular the surface normal. The method is expected to be potentially applicable to the nano- indentation experiments, where the zone around the indenter-crystal contact undergoes significant permanent deformation, the rigorous determination of which is very important to the calculation of the indentation print area and in turn, the surface hardness. Furthermore, the technique is expected to account for the plastic history of the surface displacement under the indenter. Other potential applications are mentioned in the text.

  3. Holographic particle velocimetry - A 3D measurement technique for vortex interactions, coherent structures and turbulence

    NASA Astrophysics Data System (ADS)

    Meng, Hui; Hussain, Fazle

    1991-10-01

    To understand the topology and dynamics of coherent structures (CS), the interactions of CS with fine-scale turbulence, and the effects of CS on entrainment, mixing and combustion, experimental tools are needed that can measure velocity (preferably vorticity) vector fields in both 3D space and time. While traditional measurement techniques are not able to serve this purpose, holographic particle velocimetry (HPV) appears to be promising. In a demonstration experiment, the instantaneous 3D velocity vector fields in some simple vortical flows have been obtained using the HPV technique. In this preliminary report, the principles of the HPV technique are illustrated and the key issues in its implementation are discussed.

  4. Hollow Cone Electron Imaging for Single Particle 3D Reconstruction of Proteins

    NASA Astrophysics Data System (ADS)

    Tsai, Chun-Ying; Chang, Yuan-Chih; Lobato, Ivan; van Dyck, Dirk; Chen, Fu-Rong

    2016-06-01

    The main bottlenecks for high-resolution biological imaging in electron microscopy are radiation sensitivity and low contrast. The phase contrast at low spatial frequencies can be enhanced by using a large defocus but this strongly reduces the resolution. Recently, phase plates have been developed to enhance the contrast at small defocus but electrical charging remains a problem. Single particle cryo-electron microscopy is mostly used to minimize the radiation damage and to enhance the resolution of the 3D reconstructions but it requires averaging images of a massive number of individual particles. Here we present a new route to achieve the same goals by hollow cone dark field imaging using thermal diffuse scattered electrons giving about a 4 times contrast increase as compared to bright field imaging. We demonstrate the 3D reconstruction of a stained GroEL particle can yield about 13.5 Å resolution but using a strongly reduced number of images.

  5. Hollow Cone Electron Imaging for Single Particle 3D Reconstruction of Proteins.

    PubMed

    Tsai, Chun-Ying; Chang, Yuan-Chih; Lobato, Ivan; Van Dyck, Dirk; Chen, Fu-Rong

    2016-01-01

    The main bottlenecks for high-resolution biological imaging in electron microscopy are radiation sensitivity and low contrast. The phase contrast at low spatial frequencies can be enhanced by using a large defocus but this strongly reduces the resolution. Recently, phase plates have been developed to enhance the contrast at small defocus but electrical charging remains a problem. Single particle cryo-electron microscopy is mostly used to minimize the radiation damage and to enhance the resolution of the 3D reconstructions but it requires averaging images of a massive number of individual particles. Here we present a new route to achieve the same goals by hollow cone dark field imaging using thermal diffuse scattered electrons giving about a 4 times contrast increase as compared to bright field imaging. We demonstrate the 3D reconstruction of a stained GroEL particle can yield about 13.5 Å resolution but using a strongly reduced number of images. PMID:27292544

  6. Hollow Cone Electron Imaging for Single Particle 3D Reconstruction of Proteins

    PubMed Central

    Tsai, Chun-Ying; Chang, Yuan-Chih; Lobato, Ivan; Van Dyck, Dirk; Chen, Fu-Rong

    2016-01-01

    The main bottlenecks for high-resolution biological imaging in electron microscopy are radiation sensitivity and low contrast. The phase contrast at low spatial frequencies can be enhanced by using a large defocus but this strongly reduces the resolution. Recently, phase plates have been developed to enhance the contrast at small defocus but electrical charging remains a problem. Single particle cryo-electron microscopy is mostly used to minimize the radiation damage and to enhance the resolution of the 3D reconstructions but it requires averaging images of a massive number of individual particles. Here we present a new route to achieve the same goals by hollow cone dark field imaging using thermal diffuse scattered electrons giving about a 4 times contrast increase as compared to bright field imaging. We demonstrate the 3D reconstruction of a stained GroEL particle can yield about 13.5 Å resolution but using a strongly reduced number of images. PMID:27292544

  7. 3-D reconstruction of pre-characterized lithium and tungsten dust particle trajectories in NSTX

    NASA Astrophysics Data System (ADS)

    Nichols, J.; Roquemore, A. L.; Davis, W.; Mansfield, D. K.; Skinner, C. H.; Feibush, E.; Boeglin, W.; Patel, R.; Abolafia, D.; Hartzfeld, K.; Maqueda, R.

    2011-08-01

    Calibrated amounts of 40 μm lithium dust and 10 μm tungsten powder have been dropped from above into the SOL of the National Spherical Torus Experiment (NSTX) to benchmark modeling of dust dynamics and transport. By combining the output from two visible-range fast cameras, 3-D trajectories are reliably obtained and have resulted in the generation of several hundred individual particle tracks. Particles are observed to undergo a variety of accelerations both parallel and perpendicular to the magnetic field, as well as abrupt large-angle changes in direction. All tracks obtained to date display particle motion that is constrained to within a few centimeters of the last closed flux surface. The 3-D trajectories are presented and compared to the location of the last closed flux surface as determined by EFIT.

  8. Multisensor fusion for 3D target tracking using track-before-detect particle filter

    NASA Astrophysics Data System (ADS)

    Moshtagh, Nima; Romberg, Paul M.; Chan, Moses W.

    2015-05-01

    This work presents a novel fusion mechanism for estimating the three-dimensional trajectory of a moving target using images collected by multiple imaging sensors. The proposed projective particle filter avoids the explicit target detection prior to fusion. In projective particle filter, particles that represent the posterior density (of target state in a high-dimensional space) are projected onto the lower-dimensional observation space. Measurements are generated directly in the observation space (image plane) and a marginal (sensor) likelihood is computed. The particles states and their weights are updated using the joint likelihood computed from all the sensors. The 3D state estimate of target (system track) is then generated from the states of the particles. This approach is similar to track-before-detect particle filters that are known to perform well in tracking dim and stealthy targets in image collections. Our approach extends the track-before-detect approach to 3D tracking using the projective particle filter. The performance of this measurement-level fusion method is compared with that of a track-level fusion algorithm using the projective particle filter. In the track-level fusion algorithm, the 2D sensor tracks are generated separately and transmitted to a fusion center, where they are treated as measurements to the state estimator. The 2D sensor tracks are then fused to reconstruct the system track. A realistic synthetic scenario with a boosting target was generated, and used to study the performance of the fusion mechanisms.

  9. Dynamic 3D simulations of earthquakes on en echelon faults

    USGS Publications Warehouse

    Harris, R.A.; Day, S.M.

    1999-01-01

    One of the mysteries of earthquake mechanics is why earthquakes stop. This process determines the difference between small and devastating ruptures. One possibility is that fault geometry controls earthquake size. We test this hypothesis using a numerical algorithm that simulates spontaneous rupture propagation in a three-dimensional medium and apply our knowledge to two California fault zones. We find that the size difference between the 1934 and 1966 Parkfield, California, earthquakes may be the product of a stepover at the southern end of the 1934 earthquake and show how the 1992 Landers, California, earthquake followed physically reasonable expectations when it jumped across en echelon faults to become a large event. If there are no linking structures, such as transfer faults, then strike-slip earthquakes are unlikely to propagate through stepovers >5 km wide. Copyright 1999 by the American Geophysical Union.

  10. 3-D Coldtest Simulations of W-Band Devices

    NASA Astrophysics Data System (ADS)

    Petillo, John J.

    1998-04-01

    Recent advances in the simulation and modeling of vacuum electronic devices has opened up the door for the much-anticipated ``first-pass design success" of these devices by computational means. For W-band, this is an exciting prospect since many current designs begin with scaling from lower-frequency devices. The next step is trial and error by fabrication and experimental testing, which demands high precision machining practices as well as reproducible assembly procedures. This can be a very time-consuming, expensive, and error-prone undertaking. Aside from the unit-to-unit variances in dimensions between prototypes, the repeatability of material properties (e.g., dielectric constant and resistivity) in this frequency band is difficult to measure, predict, and ensure. Computers have no such repeatability problems, and run-to-run differences can be predefined, and subsequently analyzed. The computational limitations include setup difficulty, computer time, computer power, and fidelity of the calculation. For W-band devices with 1% or lower bandwidth, fidelity of the numerical solution is critical in the search for ``first-pass design success" capability. The codes have to not only be up to the task, but the appropriate application of these codes becomes an essential ingredient to predicting the accuracy of the results. These issues as well as cross validation between computational models will be discussed. Several recent simulation results of W-band devices will be presented, along with comparisons with data. These will include the NRL/Litton/CPI gyroklystron input coupler and several CPI Coupled Cavity Tube designs.

  11. 3D Direct Simulation Monte Carlo Modeling of the Spacecraft Environment of Rosetta

    NASA Astrophysics Data System (ADS)

    Bieler, A. M.; Tenishev, V.; Fougere, N.; Gombosi, T. I.; Hansen, K. C.; Combi, M. R.; Huang, Z.; Jia, X.; Toth, G.; Altwegg, K.; Wurz, P.; Jäckel, A.; Le Roy, L.; Gasc, S.; Calmonte, U.; Rubin, M.; Tzou, C. Y.; Hässig, M.; Fuselier, S.; De Keyser, J.; Berthelier, J. J.; Mall, U. A.; Rème, H.; Fiethe, B.; Balsiger, H.

    2014-12-01

    The European Space Agency's Rosetta mission is the first to escort a comet over an extended time as the comet makes its way through the inner solar system. The ROSINA instrument suite consisting of a double focusing mass spectrometer, a time of flight mass spectrometer and a pressure sensor, will provide temporally and spatially resolved data on the comet's volatile inventory. The effect of spacecraft outgassing is well known and has been measured with the ROSINA instruments onboard Rosetta throughout the cruise phase. The flux of released neutral gas originating from the spacecraft cannot be distinguished from the cometary signal by the mass spectrometers and varies significantly with solar illumination conditions. For accurate interpretation of the instrument data, a good understanding of spacecraft outgassing is necessary. In this talk we present results simulating the spacecraft environment with the Adaptive Mesh Particle Simulator (AMPS) code. AMPS is a direct simulation monte carlo code that includes multiple species in a 3D adaptive mesh to describe a full scale model of the spacecraft environment. We use the triangulated surface model of the spacecraft to implement realistic outgassing rates for different areas on the surface and take shadowing effects in consideration. The resulting particle fluxes are compared to the measurements of the ROSINA experiment and implications for ROSINA measurements and data analysis are discussed. Spacecraft outgassing has implications for future space missions to rarefied atmospheres as it imposes a limit on the detection of various species.

  12. The Virtual Radiopharmacy Laboratory: A 3-D Simulation for Distance Learning

    ERIC Educational Resources Information Center

    Alexiou, Antonios; Bouras, Christos; Giannaka, Eri; Kapoulas, Vaggelis; Nani, Maria; Tsiatsos, Thrasivoulos

    2004-01-01

    This article presents Virtual Radiopharmacy Laboratory (VR LAB), a virtual laboratory accessible through the Internet. VR LAB is designed and implemented in the framework of the VirRAD European project. This laboratory represents a 3D simulation of a radio-pharmacy laboratory, where learners, represented by 3D avatars, can experiment on…

  13. Measuring the orientation and rotation rate of 3D printed particles in turbulent flow

    NASA Astrophysics Data System (ADS)

    Voth, Greg; Marcus, Guy G.; Parsa, Shima; Kramel, Stefan; Ni, Rui; Cole, Brendan

    2014-11-01

    The orientation distribution and rotations of anisotropic particles plays a key role in many applications ranging from icy clouds to papermaking and drag reduction in pipe flow. Experimental access to time resolved orientations of anisotropic particles has not been easy to achieve. We have found that 3D printing technology can be used to fabricate a wide range of particle shapes with smallest dimension down to 300 μm. So far we have studied rods, crosses, jacks, tetrads, and helical shapes. We extract the particle orientations from stereoscopic video images using a method of least squares optimization in Euler angle space. We find that in turbulence the orientation and rotation rate of many particles can be understood using a simple picture of alignment of both the vorticity and a long axis of the particle with the Lagrangian stretching direction of the flow. This research is supported by NSF Grant DMR-1208990.

  14. Measuring the orientation and rotation rate of 3D printed particles in turbulent flow

    NASA Astrophysics Data System (ADS)

    Voth, Greg; Kramel, Stefan; Cole, Brendan

    2015-03-01

    The orientation distribution and rotations of anisotropic particles plays a key role in many applications ranging from icy clouds to papermaking and drag reduction in pipe flow. Experimental access to time resolved orientations of anisotropic particles has not been easy to achieve. We have found that 3D printing technology can be used to fabricate a wide range of particle shapes with smallest dimension down to 300 ?m. So far we have studied rods, crosses, jacks, tetrads, and helical shapes. We extract the particle orientations from stereoscopic video images using a method of least squares optimization in Euler angle space. We find that in turbulence the orientation and rotation rate of many particles can be understood using a simple picture of alignment of both the vorticity and a long axis of the particle with the Lagrangian stretching direction of the flow.

  15. 3D Plasma Clusters: Analysis of dynamical evolution and individual particle interaction

    SciTech Connect

    Antonova, T.; Thomas, H. M.; Morfill, G. E.; Annaratone, B. M.

    2008-09-07

    3D plasma clusters (up to 100 particles) have been built inside small (32 mm{sup 3}) plasma volume in gravity. It has been estimated that the external confinement has a negligible influence on the processes inside the clusters. At such conditions the analysis of dynamical evolution and individual particle interactions have shown that the binary interaction among particles in addition to the repelling Coulomb force exhibits also an attractive part. The tendency of the systems to approach the state with minimum energy by rearranging particles inside has been detected. The measured 63 particles' cluster vibrations are in close agreement with vibrations of a drop with surface tension. This indicates that even a 63 particle cluster already exhibits properties normally associated with the cooperative regime.

  16. Development of a 3D to 1D Particle Transport Model to Predict Deposition in the Lungs

    NASA Astrophysics Data System (ADS)

    Oakes, Jessica M.; Grandmont, Celine; Shadden, Shawn C.; Vignon-Clementel, Irene E.

    2014-11-01

    Aerosolized particles are commonly used for therapeutic drug delivery as they can be delivered to the body systemically or be used to treat lung diseases. Recent advances in computational resources have allowed for sophisticated pulmonary simulations, however it is currently impossible to solve for airflow and particle transport for all length and time scales of the lung. Instead, multi-scale methods must be used. In our recent work, where computational methods were employed to solve for airflow and particle transport in the rat airways (Oakes et al. (2014), Annals of Biomedical Engineering 42, 899), the number of particles to exit downstream of the 3D domain was determined. In this current work, the time-dependent Lagrangian description of particles was used to numerically solve a 1D convection-diffusion model (trumpet model, Taulbee and Yu (1975), Journal of Applied Physiology, 38, 77) parameterized specifically for the lung. The expansion of the airway dimensions was determined based on data collected from our aerosol exposure experiments (Oakes et al. (2014), Journal of Applied Physiology, 116, 1561). This 3D-1D framework enables us to predict the fate of particles in the whole lung. This work was supported by the Whitaker Foundation at the IIE, a INRIA Associated Team Postdoc Grant, and a UC Presidential Fellowship.

  17. Recent advances in 3D computed tomography techniques for simulation and navigation in hepatobiliary pancreatic surgery.

    PubMed

    Uchida, Masafumi

    2014-04-01

    A few years ago it could take several hours to complete a 3D image using a 3D workstation. Thanks to advances in computer science, obtaining results of interest now requires only a few minutes. Many recent 3D workstations or multimedia computers are equipped with onboard 3D virtual patient modeling software, which enables patient-specific preoperative assessment and virtual planning, navigation, and tool positioning. Although medical 3D imaging can now be conducted using various modalities, including computed tomography (CT), magnetic resonance imaging (MRI), positron emission tomography (PET), and ultrasonography (US) among others, the highest quality images are obtained using CT data, and CT images are now the most commonly used source of data for 3D simulation and navigation image. If the 2D source image is bad, no amount of 3D image manipulation in software will provide a quality 3D image. In this exhibition, the recent advances in CT imaging technique and 3D visualization of the hepatobiliary and pancreatic abnormalities are featured, including scan and image reconstruction technique, contrast-enhanced techniques, new application of advanced CT scan techniques, and new virtual reality simulation and navigation imaging. PMID:24464989

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

    NASA Astrophysics Data System (ADS)

    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.

  19. Transfer-printing and host-guest properties of 3D supramolecular particle structures.

    PubMed

    Ling, Xing Yi; Phang, In Yee; Reinhoudt, David N; Vancso, G Julius; Huskens, Jurriaan

    2009-04-01

    Mechanically robust and crystalline supramolecular particle structures have been constructed by decoupling nanoparticle assembly and supramolecular glue infiltration into a sequential process. First, beta-cyclodextrin (CD)-functionalized polystyrene particles (d approximately 500 nm) were assembled on a CD-functionalized surface via convective assembly to form highly ordered, but mechanically unstable, particle crystals. Subsequently, the crystals were infiltrated by a solution of adamantyl-functionalized dendrimers, functioning as a supramolecular glue to bind neighboring particles together and to couple the entire particle crystal to the CD surface, both in a noncovalent manner. The supramolecular particle crystals are highly robust, as witnessed by their ability to withstand agitation by ultrasonication. When assembled on a poly(dimethylsiloxane) (PDMS) stamp, the dendrimer-infiltrated particle crystals could be transfer-printed onto a CD-functionalized target surface. By variation of the geometry and size of the PDMS stamps, single particle lines, interconnected particle rings, and V-shaped particle assemblies were obtained. The particle structures served as 3D receptors for the binding of (multiple) complementary guest molecules, indicating that the supramolecular host functionalities of the particle crystals were retained throughout the fabrication process. PMID:20356024

  20. Ultra-high-speed 3D astigmatic particle tracking velocimetry: application to particle-laden supersonic impinging jets

    NASA Astrophysics Data System (ADS)

    Buchmann, N. A.; Cierpka, C.; Kähler, C. J.; Soria, J.

    2014-11-01

    The paper demonstrates ultra-high-speed three-component, three-dimensional (3C3D) velocity measurements of micron-sized particles suspended in a supersonic impinging jet flow. Understanding the dynamics of individual particles in such flows is important for the design of particle impactors for drug delivery or cold gas dynamic spray processing. The underexpanded jet flow is produced via a converging nozzle, and micron-sized particles ( d p = 110 μm) are introduced into the gas flow. The supersonic jet impinges onto a flat surface, and the particle impact velocity and particle impact angle are studied for a range of flow conditions and impingement distances. The imaging system consists of an ultra-high-speed digital camera (Shimadzu HPV-1) capable of recording rates of up to 1 Mfps. Astigmatism particle tracking velocimetry (APTV) is used to measure the 3D particle position (Cierpka et al., Meas Sci Technol 21(045401):13, 2010) by coding the particle depth location in the 2D images by adding a cylindrical lens to the high-speed imaging system. Based on the reconstructed 3D particle positions, the particle trajectories are obtained via a higher-order tracking scheme that takes advantage of the high temporal resolution to increase robustness and accuracy of the measurement. It is shown that the particle velocity and impingement angle are affected by the gas flow in a manner depending on the nozzle pressure ratio and stand-off distance where higher pressure ratios and stand-off distances lead to higher impact velocities and larger impact angles.

  1. 2D and 3D Mass Transfer Simulations in β Lyrae System

    NASA Astrophysics Data System (ADS)

    Nazarenko, V. V.; Glazunova, L. V.; Karetnikov, V. G.

    2001-12-01

    2D and 3D mass transfer simulations of the mass transfer in β Lyrae binary system. We have received that from a point L3 40 per cent of mass transfer from L1-point is lost.The structure of a gas envelope, around system is calculated.3-D mass transfer simulations has shown presence the spiral shock in the disk around primary star's and a jet-like structures (a mass flow in vertical direction) over a stream.

  2. 3D simulations of multipacting in the 56 MHz SRF cavity

    SciTech Connect

    Wu Q.; Belomestnykh, S.; Ge, L.; Ko, K.; Li, Z.; Ng, C.; Xiao, L.

    2012-05-20

    The 56 MHz SRF Quarter-Wave Resonator (QWR) is designed for RHIC as a storage cavity to improve the collider performance. 2D multipacting simulation has been done for the cavity alone. Ripples were added to the outer body of the cavity for multipacting suppression based on the simulation findings. During operation, there will be four higher order mode (HOM) couplers. All of these components will be exposed to high RF fields. In this paper we compare 2D and 3D codes simulation results for multipacting in the cavity. We also report 3D simulation results for multipacting simulation at the couplers.

  3. Characterization of an SRF gun: a 3D full wave simulation

    SciTech Connect

    Wang, E.; Ben-Zvi, I.; Wang, J.

    2011-03-28

    We characterized a BNL 1.3GHz half-cell SRF gun is tested for GaAs photocathode. The gun already was simulated several years ago via two-dimensional (2D) numerical codes (i.e., Superfish and Parmela) with and without the beam. In this paper, we discuss our investigation of its characteristics using a three dimensional (3D) full-wave code (CST STUDIO SUITE{trademark}).The input/pickup couplers are sited symmetrically on the same side of the gun at an angle of 180{sup o}. In particular, the inner conductor of the pickup coupler is considerably shorter than that of the input coupler. We evaluated the cross-talk between the beam (trajectory) and the signal on the input coupler compared our findings with published results based on analytical models. The CST STUDIO SUITE{trademark} also was used to predict the field within the cavity; particularly, a combination of transient/eigenmode solvers was employed to accurately construct the RF field for the particles, which also includes the effects of the couplers. Finally, we explored the beam's dynamics with a particle in cell (PIC) simulation, validated the results and compare them with 2D code result.

  4. 3D real-time visualization of blood flow in cerebral aneurysms by light field particle image velocimetry

    NASA Astrophysics Data System (ADS)

    Carlsohn, Matthias F.; Kemmling, André; Petersen, Arne; Wietzke, Lennart

    2016-04-01

    Cerebral aneurysms require endovascular treatment to eliminate potentially lethal hemorrhagic rupture by hemostasis of blood flow within the aneurysm. Devices (e.g. coils and flow diverters) promote homeostasis, however, measurement of blood flow within an aneurysm or cerebral vessel before and after device placement on a microscopic level has not been possible so far. This would allow better individualized treatment planning and improve manufacture design of devices. For experimental analysis, direct measurement of real-time microscopic cerebrovascular flow in micro-structures may be an alternative to computed flow simulations. An application of microscopic aneurysm flow measurement on a regular basis to empirically assess a high number of different anatomic shapes and the corresponding effect of different devices would require a fast and reliable method at low cost with high throughout assessment. Transparent three dimensional 3D models of brain vessels and aneurysms may be used for microscopic flow measurements by particle image velocimetry (PIV), however, up to now the size of structures has set the limits for conventional 3D-imaging camera set-ups. On line flow assessment requires additional computational power to cope with the processing large amounts of data generated by sequences of multi-view stereo images, e.g. generated by a light field camera capturing the 3D information by plenoptic imaging of complex flow processes. Recently, a fast and low cost workflow for producing patient specific three dimensional models of cerebral arteries has been established by stereo-lithographic (SLA) 3D printing. These 3D arterial models are transparent an exhibit a replication precision within a submillimeter range required for accurate flow measurements under physiological conditions. We therefore test the feasibility of microscopic flow measurements by PIV analysis using a plenoptic camera system capturing light field image sequences. Averaging across a sequence of

  5. 3D numerical simulations of dispersion of volcanic ash using a Lagrangian model

    NASA Astrophysics Data System (ADS)

    Suzuki, Yujiro; Koyaguchi, Takehiro

    2014-05-01

    Dispersion of volcanic ash largely depends on the atmospheric wind speed and eruption intensity. In general, when the atmospheric wind is weak and/or eruption intensity is strong (i.e., magma discharge rate is small), the volcanic plume is characterized by the formation of umbrella cloud and the particles (i.e., volcanic ashes) are transported by the gravity current of umbrella cloud. On the other hand, if the wind is strong and/or eruption intensity is weak, the volcanic plume tends to be distorted by wind and the particles are drifted mainly by the wind. Because these effects of gravity current and wind also change depending on the particle size, it is difficult to quantitatively predict the distributions of particles suspended in the atmosphere and those deposited on the ground. In this study, we are developing a 3-D numerical model which directly simulates the ash transport and deposition. The model is designed to simulate the injection of a mixture of solid pyroclasts and volcanic gas from a circular vent above a flat surface in a stratified atmosphere, using a combination of a pseudo-gas model for fluid motion and a Lagrangian model for particle motion. During fluid dynamics calculations, we ignore the separation of solid pyroclasts from the eruption cloud, treating an eruption cloud as a single gas with a density calculated using a mixing ratio between ejected material and entrained air (Suzuki et al., 2005, JGR). In order to calculate the location and movement of ash particles, we employ Lagrangian marker particles of various sizes and densities. The marker particles are ejected from the vent with the same velocity of the eruption cloud every 2 sec. The particles are accelerated or decelerated by the drag force on the spheres and fall to the ground with their terminal velocities. We carried out a simulation of a small-scale eruption in the strong wind fields with the magma discharge rate of 2.5 x 106 kg/s. The rising plume is largely distorted by wind and

  6. Functional Stereology for 3D Particle Size Distributions from 2D Observations: a Practical Approach

    NASA Astrophysics Data System (ADS)

    Proussevitch, A. A.; Sahagian, D. L.; Jutzeler, M.

    2010-12-01

    Functional stereology applies known deconvolution techniques to obtain 3D size distributions from 2D cross-section data based on an assumption that both 2D and 3D statistics have known distribution functions with unknown parameters. A new stereological approach enables us to solve this problem by utilizing function minimization to find parameters of the distribution functions. There is no limit to continuous distribution function types that could be used, such as Gaussian, Logistic, Weibull, Gamma, and others. As compared to previously known 2D to 3D conversion methods (e.g. Sahagian and Proussevitch, 1998), functional stereology has much greater practical application to non-spherical particles/objects because it is free of uncontrollable error propagation for all particles shapes. The new practical method of functional stereology has been implemented in Stereonet software adapted for both a) direct logarithmic scales of particle/voids volumes, and b) Phi units of linear dimensions (-log2 of size). Applications of the method include distribution of voids/bubbles in all types of volcanic rocks, pore networks in sedimentary rocks, mineral and grain sizes, volcanic clasts, breccia, and texture features of a wide range of rock formations. Such applications demonstrate utility of this functional stereology approach.

  7. The Effects of 3D Computer Simulation on Biology Students' Achievement and Memory Retention

    ERIC Educational Resources Information Center

    Elangovan, Tavasuria; Ismail, Zurida

    2014-01-01

    A quasi experimental study was conducted for six weeks to determine the effectiveness of two different 3D computer simulation based teaching methods, that is, realistic simulation and non-realistic simulation on Form Four Biology students' achievement and memory retention in Perak, Malaysia. A sample of 136 Form Four Biology students in Perak,…

  8. Fast and efficient particle reconstruction on a 3D grid using sparsity

    NASA Astrophysics Data System (ADS)

    Cornic, P.; Champagnat, F.; Cheminet, A.; Leclaire, B.; Le Besnerais, G.

    2015-03-01

    We propose an approach for efficient localization and intensity reconstruction of particles on a 3D grid based on sparsity principles. The computational complexity of the method is limited by using the particle volume reconstruction paradigm (Champagnat et al. in Meas Sci Technol 25, 2014) and a reduction in the problem dimension. Tests on synthetic and experimental data show that the proposed method leads to more efficient detections and to reconstructions of higher quality than classical tomoPIV approaches on a large range of seeding densities, up to ppp ≈ 0.12.

  9. 3D multiple-point statistics simulation using 2D training images

    NASA Astrophysics Data System (ADS)

    Comunian, A.; Renard, P.; Straubhaar, J.

    2012-03-01

    One of the main issues in the application of multiple-point statistics (MPS) to the simulation of three-dimensional (3D) blocks is the lack of a suitable 3D training image. In this work, we compare three methods of overcoming this issue using information coming from bidimensional (2D) training images. One approach is based on the aggregation of probabilities. The other approaches are novel. One relies on merging the lists obtained using the impala algorithm from diverse 2D training images, creating a list of compatible data events that is then used for the MPS simulation. The other (s2Dcd) is based on sequential simulations of 2D slices constrained by the conditioning data computed at the previous simulation steps. These three methods are tested on the reproduction of two 3D images that are used as references, and on a real case study where two training images of sedimentary structures are considered. The tests show that it is possible to obtain 3D MPS simulations with at least two 2D training images. The simulations obtained, in particular those obtained with the s2Dcd method, are close to the references, according to a number of comparison criteria. The CPU time required to simulate with the method s2Dcd is from two to four orders of magnitude smaller than the one required by a MPS simulation performed using a 3D training image, while the results obtained are comparable. This computational efficiency and the possibility of using MPS for 3D simulation without the need for a 3D training image facilitates the inclusion of MPS in Monte Carlo, uncertainty evaluation, and stochastic inverse problems frameworks.

  10. Local characterization of hindered Brownian motion by using digital video microscopy and 3D particle tracking.

    PubMed

    Dettmer, Simon L; Keyser, Ulrich F; Pagliara, Stefano

    2014-02-01

    In this article we present methods for measuring hindered Brownian motion in the confinement of complex 3D geometries using digital video microscopy. Here we discuss essential features of automated 3D particle tracking as well as diffusion data analysis. By introducing local mean squared displacement-vs-time curves, we are able to simultaneously measure the spatial dependence of diffusion coefficients, tracking accuracies and drift velocities. Such local measurements allow a more detailed and appropriate description of strongly heterogeneous systems as opposed to global measurements. Finite size effects of the tracking region on measuring mean squared displacements are also discussed. The use of these methods was crucial for the measurement of the diffusive behavior of spherical polystyrene particles (505 nm diameter) in a microfluidic chip. The particles explored an array of parallel channels with different cross sections as well as the bulk reservoirs. For this experiment we present the measurement of local tracking accuracies in all three axial directions as well as the diffusivity parallel to the channel axis while we observed no significant flow but purely Brownian motion. Finally, the presented algorithm is suitable also for tracking of fluorescently labeled particles and particles driven by an external force, e.g., electrokinetic or dielectrophoretic forces. PMID:24593372