Numerical modelling of gravel unconstrained flow experiments with the DAN3D and RASH3D codes
NASA Astrophysics Data System (ADS)
Sauthier, Claire; Pirulli, Marina; Pisani, Gabriele; Scavia, Claudio; Labiouse, Vincent
2015-12-01
Landslide continuum dynamic models have improved considerably in the last years, but a consensus on the best method of calibrating the input resistance parameter values for predictive analyses has not yet emerged. In the present paper, numerical simulations of a series of laboratory experiments performed at the Laboratory for Rock Mechanics of the EPF Lausanne were undertaken with the RASH3D and DAN3D numerical codes. They aimed at analysing the possibility to use calibrated ranges of parameters (1) in a code different from that they were obtained from and (2) to simulate potential-events made of a material with the same characteristics as back-analysed past-events, but involving a different volume and propagation path. For this purpose, one of the four benchmark laboratory tests was used as past-event to calibrate the dynamic basal friction angle assuming a Coulomb-type behaviour of the sliding mass, and this back-analysed value was then used to simulate the three other experiments, assumed as potential-events. The computational findings show good correspondence with experimental results in terms of characteristics of the final deposits (i.e., runout, length and width). Furthermore, the obtained best fit values of the dynamic basal friction angle for the two codes turn out to be close to each other and within the range of values measured with pseudo-dynamic tilting tests.
Code and Solution Verification of 3D Numerical Modeling of Flow in the Gust Erosion Chamber
NASA Astrophysics Data System (ADS)
Yuen, A.; Bombardelli, F. A.
2014-12-01
Erosion microcosms are devices commonly used to investigate the erosion and transport characteristics of sediments at the bed of rivers, lakes, or estuaries. In order to understand the results these devices provide, the bed shear stress and flow field need to be accurately described. In this research, the UMCES Gust Erosion Microcosm System (U-GEMS) is numerically modeled using Finite Volume Method. The primary aims are to simulate the bed shear stress distribution at the surface of the sediment core/bottom of the microcosm, and to validate the U-GEMS produces uniform bed shear stress at the bottom of the microcosm. The mathematical model equations are solved by on a Cartesian non-uniform grid. Multiple numerical runs were developed with different input conditions and configurations. Prior to developing the U-GEMS model, the General Moving Objects (GMO) model and different momentum algorithms in the code were verified. Code verification of these solvers was done via simulating the flow inside the top wall driven square cavity on different mesh sizes to obtain order of convergence. The GMO model was used to simulate the top wall in the top wall driven square cavity as well as the rotating disk in the U-GEMS. Components simulated with the GMO model were rigid bodies that could have any type of motion. In addition cross-verification was conducted as results were compared with numerical results by Ghia et al. (1982), and good agreement was found. Next, CFD results were validated by simulating the flow within the conventional microcosm system without suction and injection. Good agreement was found when the experimental results by Khalili et al. (2008) were compared. After the ability of the CFD solver was proved through the above code verification steps. The model was utilized to simulate the U-GEMS. The solution was verified via classic mesh convergence study on four consecutive mesh sizes, in addition to that Grid Convergence Index (GCI) was calculated and based on
ERIC Educational Resources Information Center
Sack, Jacqueline J.
2013-01-01
This article explicates the development of top-view numeric coding of 3-D cube structures within a design research project focused on 3-D visualization skills for elementary grades children. It describes children's conceptual development of 3-D cube structures using concrete models, conventional 2-D pictures and abstract top-view numeric…
Numerical simulation of jet aerodynamics using the three-dimensional Navier-Stokes code PAB3D
NASA Technical Reports Server (NTRS)
Pao, S. Paul; Abdol-Hamid, Khaled S.
1996-01-01
This report presents a unified method for subsonic and supersonic jet analysis using the three-dimensional Navier-Stokes code PAB3D. The Navier-Stokes code was used to obtain solutions for axisymmetric jets with on-design operating conditions at Mach numbers ranging from 0.6 to 3.0, supersonic jets containing weak shocks and Mach disks, and supersonic jets with nonaxisymmetric nozzle exit geometries. This report discusses computational methods, code implementation, computed results, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions. The Navier-Stokes method using the standard Jones-Launder two-equation kappa-epsilon turbulence model can accurately predict jet flow, and such predictions are made without any modification to the published constants for the turbulence model.
Numerical Simulation of Two-grid Ion Optics Using a 3D Code
NASA Technical Reports Server (NTRS)
Anderson, John R.; Katz, Ira; Goebel, Dan
2004-01-01
A three-dimensional ion optics code has been developed under NASA's Project Prometheus to model two grid ion optics systems. The code computes the flow of positive ions from the discharge chamber through the ion optics and into the beam downstream of the thruster. The rate at which beam ions interact with background neutral gas to form charge exchange ions is also computed. Charge exchange ion trajectories are computed to determine where they strike the ion optics grid surfaces and to determine the extent of sputter erosion they cause. The code has been used to compute predictions of the erosion pattern and wear rate on the NSTAR ion optics system; the code predicts the shape of the eroded pattern but overestimates the initial wear rate by about 50%. An example of use of the code to estimate the NEXIS thruster accelerator grid life is also presented.
Numerical model of water flow and solute accumulation in vertisols using HYDRUS 2D/3D code
NASA Astrophysics Data System (ADS)
Weiss, Tomáš; Dahan, Ofer; Turkeltub, Tuvia
2015-04-01
Keywords: dessication-crack-induced-salinization, preferential flow, conceptual model, numerical model, vadose zone, vertisols, soil water retention function, HYDRUS 2D/3D Vertisols cover a hydrologically very significant area of semi-arid regions often through which water infiltrates to groundwater aquifers. Understanding of water flow and solute accumulation is thus very relevant to agricultural activity and water resources management. Previous works suggest a conceptual model of dessication-crack-induced-salinization where salinization of sediment in the deep section of the vadose zone (up to 4 m) is induced by subsurface evaporation due to convective air flow in the dessication cracks. It suggests that the salinization is induced by the hydraulic gradient between the dry sediment in the vicinity of cracks (low potential) and the relatively wet sediment further from the main cracks (high potential). This paper presents a modified previously suggested conceptual model and a numerical model. The model uses a simple uniform flow approach but unconventionally prescribes the boundary conditions and the hydraulic parameters of soil. The numerical model is bound to one location close to a dairy farm waste lagoon, but the application of the suggested conceptual model could be possibly extended to all semi-arid regions with vertisols. Simulations were conducted using several modeling approaches with an ultimate goal of fitting the simulation results to the controlling variables measured in the field: temporal variation in water content across thick layer of unsaturated clay sediment (>10 m), sediment salinity and salinity the water draining down the vadose zone to the water table. The development of the model was engineered in several steps; all computed as forward solutions by try-and-error approach. The model suggests very deep instant infiltration of fresh water up to 12 m, which is also supported by the field data. The paper suggests prescribing a special atmospheric
NASA Astrophysics Data System (ADS)
Yoshida, Hiroyuki; Misawa, Takeharu; Takase, Kazuyuki
Two-fluid model can simulate two-phase flow by computational cost less than detailed two-phase flow simulation method such as interface tracking method or particle interaction method. Therefore, two-fluid model is useful for thermal hydraulic analysis in large-scale domain such as a rod bundle. Japan Atomic Energy Agency (JAEA) develops three dimensional two-fluid model analysis code ACE-3D that adopts boundary fitted coordinate system in order to simulate complex shape flow channel. In this paper, boiling two-phase flow analysis in a tight-lattice rod bundle was performed by ACE-3D code. The parallel computation using 126 CPUs was applied to this analysis. In the results, the void fraction, which distributes in outermost region of rod bundle, is lower than that in center region of rod bundle. The tendency of void fraction distribution agreed with the measurement results by neutron radiography qualitatively. To evaluate effects of two-phase flow model used in ACE-3D code, numerical simulation of boiling two-phase in tight-lattice rod bundle with no lift force model was also performed. From the comparison of calculated results, it was concluded that the effects of lift force model were not so large for overall void fraction distribution of tight-lattice rod bundle. However, the lift force model is important for local void fraction distribution of fuel bundles.
NASA Astrophysics Data System (ADS)
Bates, Jason; Schmitt, Andrew; Zalesak, Steve
2015-11-01
The ablative Rayleigh-Taylor (RT) instability is a key factor in the performance of directly-drive inertial-confinement-fusion (ICF) targets. Although this subject has been studied for quite some time, the accurate simulation of the ablative RT instability has proven to be a challenging task for many radiation hydrodynamics codes, particularly when it comes to capturing the ablatively-stabilized region of the linear dispersion spectrum and modeling ab initio perturbations. In this poster, we present results from recent two-dimensional numerical simulations of the ablative RT instability that were performed using the Eulerian code FastRad3D at the U.S. Naval Research Laboratory. We consider both planar and spherical geometries, low and moderate-Z target materials, different laser wavelengths and where possible, compare our findings with experiment data, linearized theory and/or results from other radiation hydrodynamics codes. Overall, we find that FastRad3D is capable of simulating the ablative RT instability quite accurately, although some uncertainties/discrepancies persist. We discuss these issues, as well as some of the numerical challenges associated with modeling this class of problems. Work supported by U.S. DOE/NNSA.
TACO3D. 3-D Finite Element Heat Transfer Code
Mason, W.E.
1992-03-04
TACO3D is a three-dimensional, finite-element program for heat transfer analysis. An extension of the two-dimensional TACO program, it can perform linear and nonlinear analyses and can be used to solve either transient or steady-state problems. The program accepts time-dependent or temperature-dependent material properties, and materials may be isotropic or orthotropic. A variety of time-dependent and temperature-dependent boundary conditions and loadings are available including temperature, flux, convection, and radiation boundary conditions and internal heat generation. Additional specialized features treat enclosure radiation, bulk nodes, and master/slave internal surface conditions (e.g., contact resistance). Data input via a free-field format is provided. A user subprogram feature allows for any type of functional representation of any independent variable. A profile (bandwidth) minimization option is available. The code is limited to implicit time integration for transient solutions. TACO3D has no general mesh generation capability. Rows of evenly-spaced nodes and rows of sequential elements may be generated, but the program relies on separate mesh generators for complex zoning. TACO3D does not have the ability to calculate view factors internally. Graphical representation of data in the form of time history and spatial plots is provided through links to the POSTACO and GRAPE postprocessor codes.
NASA Astrophysics Data System (ADS)
Hsu, Kung-Chuan; Brun, Todd
Transversal circuits are important components of fault-tolerant quantum computation. Several classes of quantum error-correcting codes are known to have transversal implementations of any logical Clifford operation. However, to achieve universal quantum computation, it would be helpful to have high-performance error-correcting codes that have a transversal implementation of some logical non-Clifford operation. The 3-D color codes are a class of topological codes that permit transversal implementation of the logical π / 8 -gate. The decoding problem of a 3-D color code can be understood as a graph-matching problem on a three-dimensional lattice. Whether this class of codes will be useful in terms of performance is still an open question. We investigate the decoding problem of 3-D color codes and analyze the performance of some possible decoders.
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.
FARGO3D: Hydrodynamics/magnetohydrodynamics code
NASA Astrophysics Data System (ADS)
Benítez Llambay, Pablo; Masset, Frédéric
2015-09-01
A successor of FARGO (ascl:1102.017), FARGO3D is a versatile HD/MHD code that runs on clusters of CPUs or GPUs, with special emphasis on protoplanetary disks. FARGO3D offers Cartesian, cylindrical or spherical geometry; 1-, 2- or 3-dimensional calculations; and orbital advection (aka FARGO) for HD and MHD calculations. As in FARGO, a simple Runge-Kutta N-body solver may be used to describe the orbital evolution of embedded point-like objects. There is no need to know CUDA; users can develop new functions in C and have them translated to CUDA automatically to run on GPUs.
Multirate 3-D subband coding of video.
Taubman, D; Zakhor, A
1994-01-01
We propose a full color video compression strategy, based on 3-D subband coding with camera pan compensation, to generate a single embedded bit stream supporting multiple decoder display formats and a wide, finely gradated range of bit rates. An experimental implementation of our algorithm produces a single bit stream, from which suitable subsets are extracted to be compatible with many decoder frame sizes and frame rates and to satisfy transmission bandwidth constraints ranging from several tens of kilobits per second to several megabits per second. Reconstructed video quality from any of these bit stream subsets is often found to exceed that obtained from an MPEG-1 implementation, operated with equivalent bit rate constraints, in both perceptual quality and mean squared error. In addition, when restricted to 2-D, the algorithm produces some of the best results available in still image compression. PMID:18291953
3D Multigroup Sn Neutron Transport Code
Energy Science and Technology Software Center (ESTSC)
2001-02-14
ATTILA is a 3D multigroup transport code with arbitrary order ansotropic scatter. The transport equation is solved in first order form using a tri-linear discontinuous spatial differencing on an arbitrary tetrahedral mesh. The overall solution technique is source iteration with DSA acceleration of the scattering source. Anisotropic boundary and internal sources may be entered in the form of spherical harmonics moments. Alpha and k eigenvalue problems are allowed, as well as fixed source problems. Forwardmore » and adjoint solutions are available. Reflective, vacumn, and source boundary conditions are available. ATTILA can perform charged particle transport calculations using slowing down (CSD) terms. ATTILA can also be used to peform infra-red steady-state calculations for radiative transfer purposes.« less
3D Multigroup Sn Neutron Transport Code
McGee, John; Wareing, Todd; Pautz, Shawn
2001-02-14
ATTILA is a 3D multigroup transport code with arbitrary order ansotropic scatter. The transport equation is solved in first order form using a tri-linear discontinuous spatial differencing on an arbitrary tetrahedral mesh. The overall solution technique is source iteration with DSA acceleration of the scattering source. Anisotropic boundary and internal sources may be entered in the form of spherical harmonics moments. Alpha and k eigenvalue problems are allowed, as well as fixed source problems. Forward and adjoint solutions are available. Reflective, vacumn, and source boundary conditions are available. ATTILA can perform charged particle transport calculations using slowing down (CSD) terms. ATTILA can also be used to peform infra-red steady-state calculations for radiative transfer purposes.
Recent update of the RPLUS2D/3D codes
NASA Technical Reports Server (NTRS)
Tsai, Y.-L. Peter
1991-01-01
The development of the RPLUS2D/3D codes is summarized. These codes utilize LU algorithms to solve chemical non-equilibrium flows in a body-fitted coordinate system. The motivation behind the development of these codes is the need to numerically predict chemical non-equilibrium flows for the National AeroSpace Plane Program. Recent improvements include vectorization method, blocking algorithms for geometric flexibility, out-of-core storage for large-size problems, and an LU-SW/UP combination for CPU-time efficiency and solution quality.
3D Elastic Seismic Wave Propagation Code
Energy Science and Technology Software Center (ESTSC)
1998-09-23
E3D is capable of simulating seismic wave propagation in a 3D heterogeneous earth. Seismic waves are initiated by earthquake, explosive, and/or other sources. These waves propagate through a 3D geologic model, and are simulated as synthetic seismograms or other graphical output.
3-D Numerical Field Calculations of CESR's Upgraded Superconducting Magnets
NASA Astrophysics Data System (ADS)
Greenwald, Zipi; Greenwald, Shlomo
1997-05-01
A 3-D numerical code( Z. Greenwald, ``BST.c 3-D Magnetic Field Calculation Numerical Code'', Cornell University Note 96-09) was used to calculate the spatial magnetic fields generated by a current carrying wire. In particular, the code calculates the fields of wire loops wrapped on a pipe similar to superconductive magnet structures. The arrangement and dimensions of the loops can be easily modified to create dipoles, quadrupoles, skew magnets etc., and combinations of the above. In this paper we show the calculated 3-D fields of ironless superconducting quadrupole dipole combination designed for CESR phase III upgrade (which will be manufactured by TESLA). Since the magnet poles are made of loops, the fields at the edges are not only distorted but have a component, B_z, in the z direction as well. This Bz field can cause X-Y coupling of the beam. In order to calculate the coupling, the particle trajectories through the whole magnet were computed. The code is also used to calculate local fields errors due to possible manufacturing imperfections. An example of a rotational error of one pole, and an example of an error in the winding width are shown.
Numerical study on 3D composite morphing actuators
NASA Astrophysics Data System (ADS)
Oishi, Kazuma; Saito, Makoto; Anandan, Nishita; Kadooka, Kevin; Taya, Minoru
2015-04-01
There are a number of actuators using the deformation of electroactive polymer (EAP), where fewer papers seem to have focused on the performance of 3D morphing actuators based on the analytical approach, due mainly to their complexity. The present paper introduces a numerical analysis approach on the large scale deformation and motion of a 3D half dome shaped actuator composed of thin soft membrane (passive material) and EAP strip actuators (EAP active coupon with electrodes on both surfaces), where the locations of the active EAP strips is a key parameter. Simulia/Abaqus Static and Implicit analysis code, whose main feature is the high precision contact analysis capability among structures, are used focusing on the whole process of the membrane to touch and wrap around the object. The unidirectional properties of the EAP coupon actuator are used as input data set for the material properties for the simulation and the verification of our numerical model, where the verification is made as compared to the existing 2D solution. The numerical results can demonstrate the whole deformation process of the membrane to wrap around not only smooth shaped objects like a sphere or an egg, but also irregularly shaped objects. A parametric study reveals the proper placement of the EAP coupon actuators, with the modification of the dome shape to induce the relevant large scale deformation. The numerical simulation for the 3D soft actuators shown in this paper could be applied to a wider range of soft 3D morphing actuators.
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
NASA Astrophysics Data System (ADS)
Audigane, Pascal; Chiaberge, Christophe; Mathurin, Frédéric; Lions, Julie; Picot-Colbeaux, Géraldine
2011-04-01
This paper is addressed to the TOUGH2 user community. It presents a new tool for handling simulations run with the TOUGH2 code with specific application to CO 2 geological storage. This tool is composed of separate FORTRAN subroutines (or modules) that can be run independently, using input and output files in ASCII format for TOUGH2. These modules have been developed specifically for modeling of carbon dioxide geological storage and their use with TOUGH2 and the Equation of State module ECO2N, dedicated to CO 2-water-salt mixture systems, with TOUGHREACT, which is an adaptation of TOUGH2 with ECO2N and geochemical fluid-rock interactions, and with TOUGH2 and the EOS7C module dedicated to CO 2-CH 4 gas mixture is described. The objective is to save time for the pre-processing, execution and visualization of complex geometry for geological system representation. The workflow is rapid and user-friendly and future implementation to other TOUGH2 EOS modules for other contexts (e.g. nuclear waste disposal, geothermal production) is straightforward. Three examples are shown for validation: (i) leakage of CO 2 up through an abandoned well; (ii) 3D reactive transport modeling of CO 2 in a sandy aquifer formation in the Sleipner gas Field, (North Sea, Norway); and (iii) an estimation of enhanced gas recovery technology using CO 2 as the injected and stored gas to produce methane in the K12B Gas Field (North Sea, Denmark).
Streamlining of the RELAP5-3D Code
Mesina, George L; Hykes, Joshua; Guillen, Donna Post
2007-11-01
RELAP5-3D is widely used by the nuclear community to simulate general thermal hydraulic systems and has proven to be so versatile that the spectrum of transient two-phase problems that can be analyzed has increased substantially over time. To accommodate the many new types of problems that are analyzed by RELAP5-3D, both the physics and numerical methods of the code have been continuously improved. In the area of computational methods and mathematical techniques, many upgrades and improvements have been made decrease code run time and increase solution accuracy. These include vectorization, parallelization, use of improved equation solvers for thermal hydraulics and neutron kinetics, and incorporation of improved library utilities. In the area of applied nuclear engineering, expanded capabilities include boron and level tracking models, radiation/conduction enclosure model, feedwater heater and compressor components, fluids and corresponding correlations for modeling Generation IV reactor designs, and coupling to computational fluid dynamics solvers. Ongoing and proposed future developments include improvements to the two-phase pump model, conversion to FORTRAN 90, and coupling to more computer programs. This paper summarizes the general improvements made to RELAP5-3D, with an emphasis on streamlining the code infrastructure for improved maintenance and development. With all these past, present and planned developments, it is necessary to modify the code infrastructure to incorporate modifications in a consistent and maintainable manner. Modifying a complex code such as RELAP5-3D to incorporate new models, upgrade numerics, and optimize existing code becomes more difficult as the code grows larger. The difficulty of this as well as the chance of introducing errors is significantly reduced when the code is structured. To streamline the code into a structured program, a commercial restructuring tool, FOR_STRUCT, was applied to the RELAP5-3D source files. The
The CONV-3D code for DNS CFD calculation
NASA Astrophysics Data System (ADS)
Chudanov, Vladimir; ALCF ThermHydraX Team
2014-03-01
The CONV-3D code for DNS CFD calculation of thermal and hydrodynamics on Fast Reactor with use of supercomputers is developed. This code is highly effective in a scalability at the high performance computers such as ``Chebyshev'', ``Lomonosov'' (Moscow State University, Russia), Blue Gene/Q(ALCF MIRA, ANL). The scalability is reached up to 106 processors. The code was validated on a series of the well known tests in a wide range of Rayleigh (106-1016) and Reynolds (103-105. Such code was validated on the blind tests OECD/NEA of the turbulent intermixing in horizontal subchannels of the fuel assembly at normal pressure and temperature (Matis-H), of the flows in T-junction and the report IBRAE/ANL was published. The good coincidence of numerical predictions with experimental data was reached, that specifies applicability of the developed approach for a prediction of thermal and hydrodynamics in a boundary layer at small Prandtl that is characteristic of the liquid metal reactors. Project Name: ThermHydraX. Project Title: U.S.-Russia Collaboration on Cross-Verification and Validation in Thermal Hydraulics.
Interface requirements to couple thermal-hydraulic codes to 3D neutronic codes
Langenbuch, S.; Austregesilo, H.; Velkov, K.
1997-07-01
The present situation of thermalhydraulics codes and 3D neutronics codes is briefly described and general considerations for coupling of these codes are discussed. Two different basic approaches of coupling are identified and their relative advantages and disadvantages are discussed. The implementation of the coupling for 3D neutronics codes in the system ATHLET is presented. Meanwhile, this interface is used for coupling three different 3D neutronics codes.
Numerical analysis of 3-D potential flow in centrifugal turbomachines
NASA Astrophysics Data System (ADS)
Daiguji, H.
1983-09-01
A numerical method is developed for analysing a three-dimensional steady incompressible potential flow through an impeller in centrifugal turbomachines. The method is the same as the previous method which was developed for the axial flow turbomachines, except for some treatments in the downstream region. In order to clarify the validity and limitation of the method, a comparison with the existing experimental data and numerical results is made for radial flow compressor impellers. The calculated blade surface pressure distributions almost coincide with the quasi-3-D calculation by Krimerman and Adler (1978), but are different partly from the quasi-3-D calculation using one meridional flow analysis. It is suggested from this comparison that the flow through an impeller with high efficiency near the design point can be predicted by this fully 3-D numerical method.
RELAP5-3D code validation for RBMK phenomena
Fisher, J.E.
1999-09-01
The RELAP5-3D thermal-hydraulic code was assessed against Japanese Safety Experiment Loop (SEL) and Heat Transfer Loop (HTL) tests. These tests were chosen because the phenomena present are applicable to analyses of Russian RBMK reactor designs. The assessment cases included parallel channel flow fluctuation tests at reduced and normal water levels, a channel inlet pipe rupture test, and a high power, density wave oscillation test. The results showed that RELAP5-3D has the capability to adequately represent these RBMK-related phenomena.
RELAP5-3D Code Validation for RBMK Phenomena
Fisher, James Ebberly
1999-09-01
The RELAP5-3D thermal-hydraulic code was assessed against Japanese Safety Experiment Loop (SEL) and Heat Transfer Loop (HTL) tests. These tests were chosen because the phenomena present are applicable to analyses of Russian RBMK reactor designs. The assessment cases included parallel channel flow fluctuation tests at reduced and normal water levels, a channel inlet pipe rupture test, and a high power, density wave oscillation test. The results showed that RELAP5-3D has the capability to adequately represent these RBMK-related phenomena.
3D numerical model for NGC 6888 Nebula
NASA Astrophysics Data System (ADS)
Reyes-Iturbide, J.; Velázquez, P. F.; Rosado, M.
We present 3D numerical simulations of the NGC6888 nebula considering the proper motion and the evolution of the star, from the red supergiant (RSG) to the Wolf-Rayet (WR) phase. Our simulations reproduce the limb-brightened morphology observed in [OIII] and X-ray emission maps. The synthetic maps computed by the numerical simulations show filamentary and clumpy structures produced by instabilities triggered in the interaction between the WR wind and the RSG shell.
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.
Beam Optics Analysis - An Advanced 3D Trajectory Code
Ives, R. Lawrence; Bui, Thuc; Vogler, William; Neilson, Jeff; Read, Mike; Shephard, Mark; Bauer, Andrew; Datta, Dibyendu; Beal, Mark
2006-01-03
Calabazas Creek Research, Inc. has completed initial development of an advanced, 3D program for modeling electron trajectories in electromagnetic fields. The code is being used to design complex guns and collectors. Beam Optics Analysis (BOA) is a fully relativistic, charged particle code using adaptive, finite element meshing. Geometrical input is imported from CAD programs generating ACIS-formatted files. Parametric data is inputted using an intuitive, graphical user interface (GUI), which also provides control of convergence, accuracy, and post processing. The program includes a magnetic field solver, and magnetic information can be imported from Maxwell 2D/3D and other programs. The program supports thermionic emission and injected beams. Secondary electron emission is also supported, including multiple generations. Work on field emission is in progress as well as implementation of computer optimization of both the geometry and operating parameters. The principle features of the program and its capabilities are presented.
Beam Optics Analysis — An Advanced 3D Trajectory Code
NASA Astrophysics Data System (ADS)
Ives, R. Lawrence; Bui, Thuc; Vogler, William; Neilson, Jeff; Read, Mike; Shephard, Mark; Bauer, Andrew; Datta, Dibyendu; Beal, Mark
2006-01-01
Calabazas Creek Research, Inc. has completed initial development of an advanced, 3D program for modeling electron trajectories in electromagnetic fields. The code is being used to design complex guns and collectors. Beam Optics Analysis (BOA) is a fully relativistic, charged particle code using adaptive, finite element meshing. Geometrical input is imported from CAD programs generating ACIS-formatted files. Parametric data is inputted using an intuitive, graphical user interface (GUI), which also provides control of convergence, accuracy, and post processing. The program includes a magnetic field solver, and magnetic information can be imported from Maxwell 2D/3D and other programs. The program supports thermionic emission and injected beams. Secondary electron emission is also supported, including multiple generations. Work on field emission is in progress as well as implementation of computer optimization of both the geometry and operating parameters. The principle features of the program and its capabilities are presented.
The numerical measure of symmetry for 3D stick creatures.
Jaśkowski, Wojciech; Komosinski, Maciej
2008-01-01
This work introduces a numerical, continuous measure of symmetry for 3D stick creatures and solid 3D objects. Background information about the property of symmetry is provided, and motivations for developing a symmetry measure are described. Three approaches are mentioned, and two of them are presented in detail using formal mathematical language. The best approach is used to sort a set of creatures according to their symmetry. Experiments with a mixed set of 84 individuals originating from both human design and evolution are performed to examine symmetry within these two sources, and to determine if human designers and evolutionary processes prefer symmetry or asymmetry. PMID:18573069
Towards a 3D Space Radiation Transport Code
NASA Technical Reports Server (NTRS)
Wilson, J. W.; Tripathl, R. K.; Cicomptta, F. A.; Heinbockel, J. H.; Tweed, J.
2002-01-01
High-speed computational procedures for space radiation shielding have relied on asymptotic expansions in terms of the off-axis scatter and replacement of the general geometry problem by a collection of flat plates. This type of solution was derived for application to human rated systems in which the radius of the shielded volume is large compared to the off-axis diffusion limiting leakage at lateral boundaries. Over the decades these computational codes are relatively complete and lateral diffusion effects are now being added. The analysis for developing a practical full 3D space shielding code is presented.
CALTRANS: A parallel, deterministic, 3D neutronics code
Carson, L.; Ferguson, J.; Rogers, J.
1994-04-01
Our efforts to parallelize the deterministic solution of the neutron transport equation has culminated in a new neutronics code CALTRANS, which has full 3D capability. In this article, we describe the layout and algorithms of CALTRANS and present performance measurements of the code on a variety of platforms. Explicit implementation of the parallel algorithms of CALTRANS using both the function calls of the Parallel Virtual Machine software package (PVM 3.2) and the Meiko CS-2 tagged message passing library (based on the Intel NX/2 interface) are provided in appendices.
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.
Axisymmetric Implementation for 3D-Based DSMC Codes
NASA Technical Reports Server (NTRS)
Stewart, Benedicte; Lumpkin, F. E.; LeBeau, G. J.
2011-01-01
The primary objective in developing NASA s DSMC Analysis Code (DAC) was to provide a high fidelity modeling tool for 3D rarefied flows such as vacuum plume impingement and hypersonic re-entry flows [1]. The initial implementation has been expanded over time to offer other capabilities including a novel axisymmetric implementation. Because of the inherently 3D nature of DAC, this axisymmetric implementation uses a 3D Cartesian domain and 3D surfaces. Molecules are moved in all three dimensions but their movements are limited by physical walls to a small wedge centered on the plane of symmetry (Figure 1). Unfortunately, far from the axis of symmetry, the cell size in the direction perpendicular to the plane of symmetry (the Z-direction) may become large compared to the flow mean free path. This frequently results in inaccuracies in these regions of the domain. A new axisymmetric implementation is presented which aims to solve this issue by using Bird s approach for the molecular movement while preserving the 3D nature of the DAC software [2]. First, the computational domain is similar to that previously used such that a wedge must still be used to define the inflow surface and solid walls within the domain. As before molecules are created inside the inflow wedge triangles but they are now rotated back to the symmetry plane. During the move step, molecules are moved in 3D but instead of interacting with the wedge walls, the molecules are rotated back to the plane of symmetry at the end of the move step. This new implementation was tested for multiple flows over axisymmetric shapes, including a sphere, a cone, a double cone and a hollow cylinder. Comparisons to previous DSMC solutions and experiments, when available, are made.
3D numerical modeling of India-Asia-like collision
NASA Astrophysics Data System (ADS)
-Erika Püsök, Adina; Kaus, Boris; Popov, Anton
2013-04-01
above a strong mantle lithosphere - the jelly sandwich model (Burov and Watts, 2006). 3D models are thus needed to investigate these hypotheses. However, fully 3D models of the dynamics of continent collision zones have only been developed very recently, and presently most research groups have relied on certain explicit assumptions for their codes. Here, we employ the parallel 3D code LaMEM (Lithosphere and Mantle Evolution Model), with a finite difference staggered grid solver, which is capable of simulating lithospheric deformation while simultaneously taking mantle flow and a free surface into account. We here report on first lithospheric and upper-mantle scale simulations in which the Indian lithosphere is indented into Asia. Acknowledgements. Funding was provided by the European Research Council under the European Community's Seventh Framework Program (FP7/2007-2013) / ERC Grant agreement #258830. Numerical computations have been performed on JUQUEEN of the Jülich high-performance computing center. • Beaumont, C., Jamieson, R.A., Nguyen, M.H., Medvedev, S.E., 2004. Crustal channel flows: 1. Numerical models with applications to the tectonics of the Himalayan-Tibetan orogeny. J. Geophys. Res. 109, B06406. • Burov, E. & Watts, W.S., 2006. The long-term strength of continental lithosphere: "jelly sandwich" or "crème brûlée"?. GSA Today, 16, doi: 10.1130/1052-5173(2006)1016<1134:TLTSOC>1132.1130.CO;1132. • England P., Houseman, G., 1986. Finite strain calculations of continental deformation. 2. Comparison with the India-Asia collision zone. J. Geophys. Res.- Solid Earth and Planets 91 (B3), 3664-3676. • Jackson, J., 2002. Strength of the continental lithosphere: time to abandon the jelly sandwich?. GSA Today, September, 4-10. • Lechmann, S.M., May, D.A., Kaus, B.J.P., Schmalholz, S.M., 2011. Comparing thin-sheet models with 3D multilayer models for continental collision. Geophy. Int. J. doi: 10.1111/j.1365-246X.2011.05164.x • Royden, L.H., Burchfiel, B
3D Finite Element Trajectory Code with Adaptive Meshing
NASA Astrophysics Data System (ADS)
Ives, Lawrence; Bui, Thuc; Vogler, William; Bauer, Andy; Shephard, Mark; Beal, Mark; Tran, Hien
2004-11-01
Beam Optics Analysis, a new, 3D charged particle program is available and in use for the design of complex, 3D electron guns and charged particle devices. The code reads files directly from most CAD and solid modeling programs, includes an intuitive Graphical User Interface (GUI), and a robust mesh generator that is fully automatic. Complex problems can be set up, and analysis initiated in minutes. The program includes a user-friendly post processor for displaying field and trajectory data using 3D plots and images. The electrostatic solver is based on the standard nodal finite element method. The magnetostatic field solver is based on the vector finite element method and is also called during the trajectory simulation process to solve for self magnetic fields. The user imports the geometry from essentially any commercial CAD program and uses the GUI to assign parameters (voltages, currents, dielectric constant) and designate emitters (including work function, emitter temperature, and number of trajectories). The the mesh is generated automatically and analysis is performed, including mesh adaptation to improve accuracy and optimize computational resources. This presentation will provide information on the basic structure of the code, its operation, and it's capabilities.
Depth-controlled 3D TV image coding
NASA Astrophysics Data System (ADS)
Chiari, Armando; Ciciani, Bruno; Romero, Milton; Rossi, Ricardo
1998-04-01
Conventional 3D-TV codecs processing one down-compatible (either left, or right) channel may optionally include the extraction of the disparity field associated with the stereo-pairs to support the coding of the complementary channel. A two-fold improvement over such approaches is proposed in this paper by exploiting 3D features retained in the stereo-pairs to reduce the redundancies in both channels, and according to their visual sensitiveness. Through an a-priori disparity field analysis, our coding scheme separates a region of interest from the foreground/background in the volume space reproduced in order to code them selectively based on their visual relevance. Such a region of interest is here identified as the one which is focused by the shooting device. By suitably scaling the DCT coefficient n such a way that precision is reduced for the image blocks lying on less relevant areas, our approach aims at reducing the signal energy in the background/foreground patterns, while retaining finer details on the more relevant image portions. From an implementation point of view, it is worth noticing that the system proposed keeps its surplus processing power on the encoder side only. Simulation results show such improvements as a better image quality for a given transmission bit rate, or a graceful quality degradation of the reconstructed images with decreasing data-rates.
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.
Code portability and data management considerations in the SAS3D LMFBR accident-analysis code
Dunn, F.E.
1981-01-01
The SAS3D code was produced from a predecessor in order to reduce or eliminate interrelated problems in the areas of code portability, the large size of the code, inflexibility in the use of memory and the size of cases that can be run, code maintenance, and running speed. Many conventional solutions, such as variable dimensioning, disk storage, virtual memory, and existing code-maintenance utilities were not feasible or did not help in this case. A new data management scheme was developed, coding standards and procedures were adopted, special machine-dependent routines were written, and a portable source code processing code was written. The resulting code is quite portable, quite flexible in the use of memory and the size of cases that can be run, much easier to maintain, and faster running. SAS3D is still a large, long running code that only runs well if sufficient main memory is available.
FARGO3D: A New GPU-oriented MHD Code
NASA Astrophysics Data System (ADS)
Benítez-Llambay, Pablo; Masset, Frédéric S.
2016-03-01
We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on the physics of protoplanetary disks and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite-difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve the divergence-free property of the magnetic field to machine accuracy, coupled to a method of characteristics for the evaluation of electromotive forces and Lorentz forces. Orbital advection is implemented, and an N-body solver is included to simulate planets or stars interacting with the gas. We present our implementation in detail and present a number of widely known tests for comparison purposes. One strength of FARGO3D is that it can run on either graphical processing units (GPUs) or central processing units (CPUs), achieving large speed-up with respect to CPU cores. We describe our implementation choices, which allow a user with no prior knowledge of GPU programming to develop new routines for CPUs, and have them translated automatically for GPUs.
Implementation of a kappa-epsilon turbulence model to RPLUS3D code
NASA Technical Reports Server (NTRS)
Chitsomboon, Tawit
1992-01-01
The RPLUS3D code has been developed at the NASA Lewis Research Center to support the National Aerospace Plane (NASP) project. The code has the ability to solve three dimensional flowfields with finite rate combustion of hydrogen and air. The combustion process of the hydrogen-air system are simulated by an 18 reaction path, 8 species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though a fully implicit numerical algorithm is used. A k-epsilon turbulence model has recently been incorporated into the code. Initial validations have been conducted for a flow over a flat plate. Results of the validation studies are shown. Some difficulties in implementing the k-epsilon equations to the code are also discussed.
MOM3D/EM-ANIMATE - MOM3D WITH ANIMATION CODE
NASA Technical Reports Server (NTRS)
Shaeffer, J. F.
1994-01-01
compare surface-current distribution due to various initial excitation directions or electric field orientations. The program can accept up to 50 planes of field data consisting of a grid of 100 by 100 field points. These planes of data are user selectable and can be viewed individually or concurrently. With these preset limits, the program requires 55 megabytes of core memory to run. These limits can be changed in the header files to accommodate the available core memory of an individual workstation. An estimate of memory required can be made as follows: approximate memory in bytes equals (number of nodes times number of surfaces times 14 variables times bytes per word, typically 4 bytes per floating point) plus (number of field planes times number of nodes per plane times 21 variables times bytes per word). This gives the approximate memory size required to store the field and surface-current data. The total memory size is approximately 400,000 bytes plus the data memory size. The animation calculations are performed in real time at any user set time step. For Silicon Graphics Workstations that have multiple processors, this program has been optimized to perform these calculations on multiple processors to increase animation rates. The optimized program uses the SGI PFA (Power FORTRAN Accelerator) library. On single processor machines, the parallelization directives are seen as comments to the program and will have no effect on compilation or execution. MOM3D and EM-ANIMATE are written in FORTRAN 77 for interactive or batch execution on SGI series computers running IRIX 3.0 or later. The RAM requirements for these programs vary with the size of the problem being solved. A minimum of 30Mb of RAM is required for execution of EM-ANIMATE; however, the code may be modified to accommodate the available memory of an individual workstation. For EM-ANIMATE, twenty-four bit, double-buffered color capability is suggested, but not required. Sample executables and sample input and
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.
Improvements to the RELAP5-3D Nearly-Implicit Numerical Scheme
Richard A. Riemke; Walter L. Weaver; RIchard R. Schultz
2005-05-01
The RELAP5-3D computer program has been improved with regard to its nearly-implicit numerical scheme for twophase flow and single-phase flow. Changes were made to the nearly-implicit numerical scheme finite difference momentum equations as follows: (1) added the velocity flip-flop mass/energy error mitigation logic, (2) added the modified Henry-Fauske choking model, (3) used the new time void fraction in the horizontal stratification force terms and gravity head, and (4) used an implicit form of the artificial viscosity. The code modifications allow the nearly-implicit numerical scheme to be more implicit and lead to enhanced numerical stability.
RHALE: A 3-D MMALE code for unstructured grids
Peery, J.S.; Budge, K.G.; Wong, M.K.W.; Trucano, T.G.
1993-08-01
This paper describes RHALE, a multi-material arbitrary Lagrangian-Eulerian (MMALE) shock physics code. RHALE is the successor to CTH, Sandia`s 3-D Eulerian shock physics code, and will be capable of solving problems that CTH cannot adequately address. We discuss the Lagrangian solid mechanics capabilities of RHALE, which include arbitrary mesh connectivity, superior artificial viscosity, and improved material models. We discuss the MMALE algorithms that have been extended for arbitrary grids in both two- and three-dimensions. The MMALE addition to RHALE provides the accuracy of a Lagrangian code while allowing a calculation to proceed under very large material distortions. Coupling an arbitrary quadrilateral or hexahedral grid to the MMALE solution facilitates modeling of complex shapes with a greatly reduced number of computational cells. RHALE allows regions of a problem to be modeled with Lagrangian, Eulerian or ALE meshes. In addition, regions can switch from Lagrangian to ALE to Eulerian based on user input or mesh distortion. For ALE meshes, new node locations are determined with a variety of element based equipotential schemes. Element quantities are advected with donor, van Leer, or Super-B algorithms. Nodal quantities are advected with the second order SHALE or HIS algorithms. Material interfaces are determined with a modified Young`s high resolution interface tracker or the SLIC algorithm. RHALE has been used to model many problems of interest to the mechanics, hypervelocity impact, and shock physics communities. Results of a sampling of these problems are presented in this paper.
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
3-D TECATE/BREW: Thermal, stress, and birefringent ray-tracing codes for solid-state laser design
Gelinas, R.J.; Doss, S.K.; Nelson, R.G.
1994-07-20
This report describes the physics, code formulations, and numerics that are used in the TECATE (totally Eulerian code for anisotropic thermo-elasticity) and BREW (birefringent ray-tracing of electromagnetic waves) codes for laser design. These codes resolve thermal, stress, and birefringent optical effects in 3-D stationary solid-state systems. This suite of three constituent codes is a package referred to as LASRPAK.
Collision of continental corner from 3-D numerical modeling
NASA Astrophysics Data System (ADS)
Li, Zhong-Hai; Xu, Zhiqin; Gerya, Taras; Burg, Jean-Pierre
2013-10-01
Continental collision has been extensively investigated with 2-D numerical models assuming infinitely wide plates or insignificant along-strike deformation in the third dimension. However, the corners of natural collision zones normally have structural characteristics that differ from linear parts of mountain belt. We conducted 3-D high-resolution numerical simulations to study the dynamics of a continental corner (lateral continental/oceanic transition zone) during subduction/collision. The results demonstrate different modes between the oceanic subduction side (continuous subduction and retreating trench) and the continental collision side (slab break-off and topography uplift). Slab break-off occurs at a depth (⩽100 km to ˜300 km) that depends on the convergence velocity. The numerical models produce lateral extrusion of the overriding crust from the collisional side to the subduction side, which is also a phenomenon recognized around natural collision of continental corners, for instance around the western corner of the Arabia-Asia collision zone and around the eastern corner of the India-Asia collision zone. Modeling results also indicate that extrusion tectonics may be driven both from above by the topography and gravitational potentials and from below by the trench retreat and asthenospheric mantle return flow, which supports the link between deep mantle dynamics and shallower crustal deformation.
3-D numerical modeling of plume-induced subduction initiation
NASA Astrophysics Data System (ADS)
Baes, Marzieh; Gerya, taras; Sobolev, Stephan
2016-04-01
Investigation of mechanisms involved in formation of a new subduction zone can help us to better understand plate tectonics. Despite numerous previous studies, it is still unclear how and where an old oceanic plate starts to subduct beneath the other plate. One of the proposed scenarios for nucleation of subduction is plume-induced subduction initiation, which was investigated in detail, using 2-D models, by Ueda et al. (2008). Recently. Gerya et al. (2015), using 3D numerical models, proposed that plume-lithosphere interaction in the Archean led to the subduction initiation and onset of plate tectonic. In this study, we aim to pursue work of Ueda et al. (2008) by incorporation of 3-D thermo-mechanical models to investigate conditions leading to oceanic subduction initiation as a result of thermal-chemical mantle plume-lithosphere interaction in the modern earth. Results of our experiments show four different deformation regimes in response to plume-lithosphere interaction, that are a) self-sustaining subduction initiation where subduction becomes self-sustained, b) freezing subduction initiation where subduction stops at shallow depths, c) slab break-off where subducting circular slab breaks off soon after formation and d) plume underplating where plume does not pass through the lithosphere but spreads beneath it (failed subduction initiation). These different regimes depend on several parameters such as plume's size, composition and temperature, lithospheric brittle/plastic strength, age of the oceanic lithosphere and presence/absence of lithospheric heterogeneities. Results show that subduction initiates and becomes self-sustained when lithosphere is older than 10 Myr and non-dimensional ratio of the plume buoyancy force and lithospheric strength above the plume is higher than 2.
NASA Astrophysics Data System (ADS)
Froger, Etienne
1993-05-01
A description of the electromagnetic behavior of a satellite subjected to an electric discharge is given using a specially developed numerical code. One of the particularities of vacuum discharges, obtained by irradiation of polymers, is the intense emission of electrons into the spacecraft environment. Electromagnetic radiation, associated with the trajectories of the particles around the spacecraft, is considered as the main source of the interference observed. In the absence of accurate orbital data and realistic ground tests, the assessment of these effects requires numerical simulation of the interaction between this electron source and the spacecraft. This is done by the GEODE particle code which is applied to characteristic configurations in order to estimate the spacecraft response to a discharge, which is simulated from a vacuum discharge model designed in laboratory. The spacecraft response to a current injection is simulated by the ALICE numerical three dimensional code. The comparison between discharge and injection effects, from the results given by the two codes, illustrates the representativity of electromagnetic susceptibility tests and the main parameters for their definition.
3D numerical analysis of crack propagation of heterogeneous notched rock under uniaxial tension
NASA Astrophysics Data System (ADS)
Wang, S. Y.; Sloan, S. W.; Sheng, D. C.; Tang, C. A.
2016-05-01
Macroscopic notches play an important role in evaluating the fracture process zone (FPZ) and the strengths of a heterogeneous rock mass. Crack initiation, propagation and coalescence for unnotched, single-notched and double-notched rock specimens are numerically simulated in a 3-D numerical model (RFPA3D). A feature of the code RFPA3D is that it can numerically simulate the evolution of cracks in three-dimensional space, as well as the heterogeneity of the rock mass. For the unnotched case, special attention is given to the complete stress-strain curve and the corresponding AE events for the failure process of rock specimen. By comparing with published experimental results, the simulation results from RFPA3D are found to be satisfactory. For the single-notched case, the effect of the length and the depth of the single notch and the thickness of the specimen on the failure mode and peak stress are evaluated. The 3D FPZ is very different from that in two dimensions. For the double-notched case, the effects of the separation distance and overlap distance of the double notches, as well as influence of the homogeneity index (m) are also investigated. As the overlap distance increases, the direction of the principal tensile stress at each notch-end changes from a perpendicular direction (tensile stress field) to a nearly parallel direction (compressive stress field), which affects the evolution of the cracks from the two notches.
3D Structured Grid Generation Codes for Turbomachinery
NASA Technical Reports Server (NTRS)
Loellbach, James; Tsung, Fu-Lin
1999-01-01
This report describes the research tasks during the past year. The research was mainly in the area of computational grid generation in support of CFD analyses of turbomachinery components. In addition to the grid generation work, a numerical simulation was obtained for the flow through a centrifugal gas compressor using an unstructured Navier-Stokes solver. Other tasks involved many different turbomachinery component analyses. These analyses were performed for NASA projects or for industrial applications. The work includes both centrifugal and axial machines, single and multiple blade rows, and steady and unsteady analyses. Over the past five years, a set of structured grid generation codes were developed that allow grids to be obtained fairly quickly for the large majority of configurations we encounter. These codes do not comprise a generalized grid generation package; they are noninteractive codes specifically designed for turbomachinery blade row geometries. But because of this limited scope, the codes are small, fast, and portable, and they can be run in the batch mode on small workstations. During the past year, these programs were used to generate computational grids were modified for a wide variety of configurations. In particular, the codes or wrote supplementary codes to improve our grid generation capabilities for multiple blade row configurations. This involves generating separate grids for each blade row, and then making them match and overlap by a few grid points at their common interface so that fluid properties are communicated across the interface. Unsteady rotor/stator analyses were performed for an axial turbine, a centrifugal compressor, and a centrifugal pump. Steady-state single-blade-row analyses were made for a study of blade sweep in transonic compressors. There was also cooperation on the application of an unstructured Navier-Stokes solver for turbomachinery flow simulations. In particular, the unstructured solver was used to analyze the
3D visualization of numeric planetary data using JMARS
NASA Astrophysics Data System (ADS)
Dickenshied, S.; Christensen, P. R.; Anwar, S.; Carter, S.; Hagee, W.; Noss, D.
2013-12-01
JMARS (Java Mission-planning and Analysis for Remote Sensing) is a free geospatial application developed by the Mars Space Flight Facility at Arizona State University. Originally written as a mission planning tool for the THEMIS instrument on board the MARS Odyssey Spacecraft, it was released as an analysis tool to the general public in 2003. Since then it has expanded to be used for mission planning and scientific data analysis by additional NASA missions to Mars, the Moon, and Vesta, and it has come to be used by scientists, researchers and students of all ages from more than 40 countries around the world. The public version of JMARS now also includes remote sensing data for Mercury, Venus, Earth, the Moon, Mars, and a number of the moons of Jupiter and Saturn. Additional datasets for asteroids and other smaller bodies are being added as they becomes available and time permits. In addition to visualizing multiple datasets in context with one another, significant effort has been put into on-the-fly projection of georegistered data over surface topography. This functionality allows a user to easily create and modify 3D visualizations of any regional scene where elevation data is available in JMARS. This can be accomplished through the use of global topographic maps or regional numeric data such as HiRISE or HRSC DTMs. Users can also upload their own regional or global topographic dataset and use it as an elevation source for 3D rendering of their scene. The 3D Layer in JMARS allows the user to exaggerate the z-scale of any elevation source to emphasize the vertical variance throughout a scene. In addition, the user can rotate, tilt, and zoom the scene to any desired angle and then illuminate it with an artificial light source. This scene can be easily overlain with additional JMARS datasets such as maps, images, shapefiles, contour lines, or scale bars, and the scene can be easily saved as a graphic image for use in presentations or publications.
Current status of the WHAMS-3D code
Belytschko, T.; Kennedy, J.M.
1987-03-01
The program WHAMS-3D is an explicit time integration program which can be used for frames, shells, plates and continua in three dimensions. Both material nonlinearities due to elasto-plastic behavior and geometric nonlinearities due to large displacements can be treated. The program has been developed to serve as a test-bed for research into methods for nonlinear structural dynamics, but it can also be used for production calculations. The program is quite compact, so it can be coupled with other codes. The program employs a finite element format, so that is possesses considerable versatility in modeling complex shapes and boundary conditions. The element library consists of the following: quadrilateral and triangular plate-shell elements, a beam element, a spring element and a hexahedral continuum element. In addition, a rigid linkage is included which permits the efficient modeling of very stiff portions of a structure, such as the bottom ring of a core barrel. In a rigid linkage, the motion of a master node defines the motion of all slave nodes linked to the master node. This option is also useful for eccentrically connected elements where the modlines of the connected elements do not coincide, as for example, in stiffeners. Time integration is performed by the central difference method. The mass matrix is diagonal (lumped), so no equations need be solved. Different time steps can be used in different parts of the mesh.
3D flow past transonic turbine cascade SE 1050 — Experiment and numerical simulations
NASA Astrophysics Data System (ADS)
Šimurda, D.; Fürst, J.; Luxa, M.
2013-08-01
This paper is concerned with experimental and numerical research on 3D flow past prismatic turbine cascade SE1050 (known in QNET network as open test case SE1050). The primary goal was to assess the influence of the inlet velocity profile on the flow structures in the interblade channel and on the flow field parameters at the cascade exit and to compare these findings to results of numerical simulations. Investigations of 3D flow past the cascade with non-uniform inlet velocity profile were carried out both experimentally and numerically at subsonic ( M 2is = 0.8) and at transonic ( M 2is = 1.2) regime at design angle of incidence. Experimental data was obtained using a traversing device with a five-hole conical probe. Numerically, the 3D flow was simulated by open source code OpenFOAM and in-house code. Analyses of experimental data and CFD simulations have revealed the development of distinctive vortex structures resulting from non-uniform inlet velocity profile. Origin of these structures results in increased loss of kinetic energy and spanwise shift of kinetic energy loss coefficient distribution. Differences found between the subsonic and the transonic case confirm earlier findings available in the literature. Results of CFD and experiments agree reasonably well.
Numerical model of sonic boom in 3D kinematic turbulence
NASA Astrophysics Data System (ADS)
Coulouvrat, François; Luquet, David; Marchiano, Régis
2015-10-01
Sonic boom is one of the key issues to be considered in the development of future supersonic or hypersonic civil aircraft concepts. The classical sonic boom, typical for Concorde with an N-wave shape and a ground amplitude of the order of 100 Pa, prevents overland flight. Future concepts target carefully shaped sonic booms with low amplitude weak shocks. However, sonic boom when perceived at the ground level is influenced not only by the aircraft characteristics, but also by atmospheric propagation. In particular, the effect of atmospheric turbulence on sonic boom propagation near the ground is not well characterized. Flight tests performed as early as the 1960s demonstrated that classical sonic booms are sensitive to atmospheric turbulence. However, this sensitivity remains only partially understood. This is related to the fact that i) turbulence is a random process that requires a statistical approach, ii) standard methods used to predict sonic booms, mainly geometrical acoustics based on ray tracing, are inadequate within the turbulent planetary boundary layer. Moreover, the ray theory fails to predict the acoustical field in many areas of interest, such as caustics or shadow zones. These zones are of major interest for sonic boom acceptability (highest levels, lateral extent of zone of impact). These limitations outline the need for a numerical approach that is sufficiently efficient to perform a large number of realizations for a statistical approach, but that goes beyond the limitations of ray theory. With this in view, a 3D one-way numerical method solving a nonlinear scalar wave equation established for heterogeneous, moving and absorbing atmosphere, is used to assess the effects of a 3D kinematic turbulence on sonic boom in various configurations. First, a plane N-wave is propagated in the free field through random realizations of kinematic fluctuations. Then the case of a more realistic Atmospheric Boundary Layer (ABL) is investigated, with a mean
ATHENA 3D: A finite element code for ultrasonic wave propagation
NASA Astrophysics Data System (ADS)
Rose, C.; Rupin, F.; Fouquet, T.; Chassignole, B.
2014-04-01
The understanding of wave propagation phenomena requires use of robust numerical models. 3D finite element (FE) models are generally prohibitively time consuming. However, advances in computing processor speed and memory allow them to be more and more competitive. In this context, EDF R&D developed the 3D version of the well-validated FE code ATHENA2D. The code is dedicated to the simulation of wave propagation in all kinds of elastic media and in particular, heterogeneous and anisotropic materials like welds. It is based on solving elastodynamic equations in the calculation zone expressed in terms of stress and particle velocities. The particularity of the code relies on the fact that the discretization of the calculation domain uses a Cartesian regular 3D mesh while the defect of complex geometry can be described using a separate (2D) mesh using the fictitious domains method. This allows combining the rapidity of regular meshes computation with the capability of modelling arbitrary shaped defects. Furthermore, the calculation domain is discretized with a quasi-explicit time evolution scheme. Thereby only local linear systems of small size have to be solved. The final step to reduce the computation time relies on the fact that ATHENA3D has been parallelized and adapted to the use of HPC resources. In this paper, the validation of the 3D FE model is discussed. A cross-validation of ATHENA 3D and CIVA is proposed for several inspection configurations. The performances in terms of calculation time are also presented in the cases of both local computer and computation cluster use.
A 3D numerical model for Kepler's supernova remnant
NASA Astrophysics Data System (ADS)
Toledo-Roy, J. C.; Esquivel, A.; Velázquez, P. F.; Reynoso, E. M.
2014-07-01
We present new 3D numerical simulations for Kepler's supernova remnant. In this work we revisit the possibility that the asymmetric shape of the remnant in X-rays is the product of a Type Ia supernova explosion which occurs inside the wind bubble previously created by an AGB companion star. Due to the large peculiar velocity of the system, the interaction of the strong AGB wind with the interstellar medium results in a bow shock structure. In this new model we propose that the AGB wind is anisotropic, with properties such as mass-loss rate and density having a latitude dependence, and that the orientation of the polar axis of the AGB star is not aligned with the direction of motion. The ejecta from the Type Ia supernova explosion is modelled using a power-law density profile, and we let the remnant evolve for 400 yr. We computed synthetic X-ray maps from the numerical results. We find that the estimated size and peculiar X-ray morphology of Kepler's supernova remnant are well reproduced by considering an AGB mass-loss rate of 10-5 M⊙ yr-1, a wind terminal velocity of 10 km s-1, an ambient medium density of 10-3 cm-3 and an explosion energy of 7 × 1050 erg. The obtained total X-ray luminosity of the remnant in this model reaches 6 × 1050 erg, which is within a factor of 2 of the observed value, and the time evolution of the luminosity shows a rate of decrease in recent decades of ˜2.4 per cent yr-1 that is consistent with the observations.
Spacecraft charging analysis with the implicit particle-in-cell code iPic3D
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.
Ababou, R.
1996-12-31
Subsurface flow processes are inherently three-dimensional and heterogeneous over many scales. Taking this into account, for instance assuming random heterogeneity in 3-D space, puts heavy constraints on numerical models. An efficient numerical code has been developed for solving the porous media flow equations, appropriately generalized to account for 3-D, random-like heterogeneity. The code is based on implicit finite differences (or finite volumes), and uses specialized versions of pre-conditioned iterative solvers that take advantage of sparseness. With Diagonally Scaled Conjugate Gradients, in particular, large systems on the order of several million equations, with randomly variable coefficients, have been solved efficiently on Cray-2 and Cray-Y/MP8 machines, in serial mode as well as parallel mode (autotasking). The present work addresses, first, the numerical aspects and computational issues associated with detailed 3-D flow simulations, and secondly, presents a specific application related to the conductivity homogenization problem (identifying a macroscale conduction law, and an equivalent or effective conductivity). Analytical expressions of effective conductivities are compared with empirical values obtained from several large scale simulations conducted for single realizations of random porous media.
Convective instability in sedimentation: 3-D numerical study
NASA Astrophysics Data System (ADS)
Yu, Xiao; Hsu, Tian-Jian; Balachandar, S.
2014-11-01
To provide a probable explanation on the field observed rapid sedimentation process near river mouths, we investigate the convective sedimentation in stably stratified saltwater using 3-D numerical simulations. Guided by the linear stability analysis, this study focuses on the nonlinear interactions of several mechanisms, which lead to various sediment finger patterns, and the effective settling velocity for sediment ranging from clay (single-particle settling velocity V0 = 0.0036 and 0.0144 mm/s, or particle diameter d = 2 and 4 μm) to silt (V0 = 0.36 mm/s, or d = 20 μm). For very fine sediment with V0 = 0.0036 mm/s, the convective instability is dominated by double diffusion, characterized by millimeter-scale fingers. Gravitational settling slightly increases the growth rate; however, it has notable effect on the downward development of vertical mixing shortly after the sediment interface migrates below the salt interface. For sediment with V0 = 0.0144 mm/s, Rayleigh-Taylor instabilities become dominant before double-diffusive modes grow sufficiently large. Centimeter-scale and highly asymmetric sediment fingers are obtained due to nonlinear interactions between different modes. For sediment with V0 = 0.36 mm/s, Rayleigh-Taylor mechanism dominates and the resulting centimeter-scale sediment fingers show a plume-like structure. The flow pattern is similar to that without ambient salt stratification. Rapid sedimentation with effective settling velocity on the order of 1 cm/s is likely driven by convective sedimentation for sediment with V0 greater than 0.1 mm/s at concentration greater than 10-20 g/L.
3D neutronic codes coupled with thermal-hydraulic system codes for PWR, and BWR and VVER reactors
Langenbuch, S.; Velkov, K.; Lizorkin, M.
1997-07-01
This paper describes the objectives of code development for coupling 3D neutronics codes with thermal-hydraulic system codes. The present status of coupling ATHLET with three 3D neutronics codes for VVER- and LWR-reactors is presented. After describing the basic features of the 3D neutronic codes BIPR-8 from Kurchatov-Institute, DYN3D from Research Center Rossendorf and QUABOX/CUBBOX from GRS, first applications of coupled codes for different transient and accident scenarios are presented. The need of further investigations is discussed.
International "Intercomparison of 3-Dimensional (3D) Radiation Codes" (13RC)
NASA Technical Reports Server (NTRS)
Cahalan, Robert F.; Einaudi, Franco (Technical Monitor)
2000-01-01
An international "Intercomparison of 3-dimensional (3D) Radiation Codes" 13RC) has been initiated. It is endorsed by the GEWEX Radiation Panel, and funded jointly by the United States Department of Energy ARM program, and by the National Aeronautics and Space Administration Radiation Sciences program. It is a 3-phase effort that has as its goals to: (1) understand the errors and limits of 3D methods; (2) provide 'baseline' cases for future 3D code development; (3) promote sharing of 3D tools; (4) derive guidelines for 3D tool selection; and (5) improve atmospheric science education in 3D radiation.
A Magnetic Diagnostic Code for 3D Fusion Equilibria
Samuel Aaron Lazerson
2012-07-27
A synthetic magnetic diagnostics code for fusion equilibria is presented. This code calculates the response of various magnetic diagnostics to the equilibria produced by the VMEC and PIES codes. This allows for treatment of equilibria with both good nested flux surfaces and those with stochastic regions. DIAGNO v2.0 builds upon previous codes through the implementation of a virtual casing principle. The codes is validated against a vacuum shot on the Large Helical Device where the vertical field was ramped. As an exercise of the code, the diagnostic response for various equilibria are calculated on the Large Helical Device (LHD).
A Magnetic Diagnostic Code for 3D Fusion Equilibria
Samuel A. Lazerson, S. Sakakibara and Y. Suzuki
2013-03-12
A synthetic magnetic diagnostics code for fusion equilibria is presented. This code calculates the response of various magnetic diagnostics to the equilibria produced by the VMEC and PIES codes. This allows for treatment of equilibria with both good nested flux surfaces and those with stochastic regions. DIAGNO v2.0 builds upon previous codes through the implementation of a virtual casing principle. The code is validated against a vacuum shot on the Large Helical Device (LHD) where the vertical field was ramped. As an exercise of the code, the diagnostic response for various equilibria are calculated on the LHD.
Numerical Results of 3-D Modeling of Moon Accumulation
NASA Astrophysics Data System (ADS)
Khachay, Yurie; Anfilogov, Vsevolod; Antipin, Alexandr
2014-05-01
For the last time for the model of the Moon usually had been used the model of mega impact in which the forming of the Earth and its sputnik had been the consequence of the Earth's collision with the body of Mercurial mass. But all dynamical models of the Earth's accumulation and the estimations after the Pb-Pb system, lead to the conclusion that the duration of the planet accumulation was about 1 milliard years. But isotopic results after the W-Hf system testify about a very early (5-10) million years, dividing of the geochemical reservoirs of the core and mantle. In [1,2] it is shown, that the account of energy dissipating by the decay of short living radioactive elements and first of all Al26,it is sufficient for heating even small bodies with dimensions about (50-100) km up to the iron melting temperature and can be realized a principal new differentiation mechanism. The inner parts of the melted preplanets can join and they are mainly of iron content, but the cold silicate fragments return to the supply zone and additionally change the content of Moon forming to silicates. Only after the increasing of the gravitational radius of the Earth, the growing area of the future Earth's core can save also the silicate envelope fragments [3]. For understanding the further system Earth-Moon evolution it is significant to trace the origin and evolution of heterogeneities, which occur on its accumulation stage.In that paper we are modeling the changing of temperature,pressure,velocity of matter flowing in a block of 3d spherical body with a growing radius. The boundary problem is solved by the finite-difference method for the system of equations, which include equations which describe the process of accumulation, the Safronov equation, the equation of impulse balance, equation Navier-Stocks, equation for above litho static pressure and heat conductivity in velocity-pressure variables using the Businesque approach.The numerical algorithm of the problem solution in velocity
A 3D multi-block structured version of the KIVA 2 code
NASA Astrophysics Data System (ADS)
Habachi, C.; Torres, A.
A numerical procedure is developed in the KIVA 2 code for calculating flows in complex geometries. Those geometries consist of an arbitrary number of 3D secondary domains which are connected with any angle to a main region. In this procedure, the governing equations are discretized on a system of partial overlapping structured grids. Calculations are performed in the different meshes of the computation domain which are linked by a fully conservative algorithm. By this numerical technique, calculations in those geometries are possible with a reasonable number of inactive cells involved by a structured code like KIVA 2. This algorithm was validated on an 1D analytical case and a 2D experimental case. It was then used for modeling an industrial problem, a two stroke engine with ports and moving boundaries.
Quantitative analysis of accuracy of seismic wave-propagation codes in 3D random scattering media
NASA Astrophysics Data System (ADS)
Galis, Martin; Imperatori, Walter; Mai, P. Martin
2013-04-01
Several recent verification studies (e.g. Day et al., 2001; Bielak et al., 2010, Chaljub et al., 2010) have demonstrated the importance of assessing the accuracy of available numerical tools at low frequency in presence of large-scale features (basins, topography, etc.). The fast progress in high-performance computing, including efficient optimization of numerical codes on petascale supercomputers, has permitted the simulation of 3D seismic wave propagation at frequencies of engineering interest (up to 10Hz) in highly heterogeneous media (e.g. Hartzell et al, 2010; Imperatori and Mai, 2013). However, high frequency numerical simulations involving random scattering media, characterized by small-scale heterogeneities, are much more challenging for most numerical methods, and their verification may therefore be even more crucial than in the low-frequency case. Our goal is to quantitatively compare the accuracy and the behavior of three different numerical codes for seismic wave propagation in 3D random scattering media at high frequency. We deploy a point source with omega-squared spectrum, and focus on the near-source region, being of great interest in strong motion seismology. We use two codes based on finite-difference method (FD1 and FD2) and one code based on support-operator method (SO). Both FD1 and FD2 are 4-th order staggered-grid finite-difference codes (for FD1 see Olsen et al., 2009; for FD2 see Moczo et al., 2007). The FD1 and FD2 codes are characterized by slightly different medium representations, since FD1 uses point values of material parameters in each FD-cell, while FD2 uses the effective material parameters at each grid-point (Moczo et al., 2002). SO is 2-nd order support-operator method (Ely et al., 2008). We considered models with random velocity perturbations described by van Karman correlation function with different correlation lengths and different standard deviations. Our results show significant variability in both phase and amplitude as
An easy implementation of displacement calculations in 3D discrete dislocation dynamics codes
NASA Astrophysics Data System (ADS)
Fivel, Marc; Depres, Christophe
2014-10-01
Barnett's coordinate-free expression of the displacement field of a triangular loop in an isotropic media is revisited in a view to be implemented in 3D discrete dislocation dynamics codes. A general meshing procedure solving the problems of non-planar loops is presented. The method is user-friendly and can be used in numerical simulations since it gives the contribution of each dislocation segment to the global displacement field without defining the connectivity of closed loops. Easy to implement in parallel calculations, this method is successfully applied to large-scale simulations.
GPU-accelerated 3D neutron diffusion code based on finite difference method
Xu, Q.; Yu, G.; Wang, K.
2012-07-01
Finite difference method, as a traditional numerical solution to neutron diffusion equation, although considered simpler and more precise than the coarse mesh nodal methods, has a bottle neck to be widely applied caused by the huge memory and unendurable computation time it requires. In recent years, the concept of General-Purpose computation on GPUs has provided us with a powerful computational engine for scientific research. In this study, a GPU-Accelerated multi-group 3D neutron diffusion code based on finite difference method was developed. First, a clean-sheet neutron diffusion code (3DFD-CPU) was written in C++ on the CPU architecture, and later ported to GPUs under NVIDIA's CUDA platform (3DFD-GPU). The IAEA 3D PWR benchmark problem was calculated in the numerical test, where three different codes, including the original CPU-based sequential code, the HYPRE (High Performance Pre-conditioners)-based diffusion code and CITATION, were used as counterpoints to test the efficiency and accuracy of the GPU-based program. The results demonstrate both high efficiency and adequate accuracy of the GPU implementation for neutron diffusion equation. A speedup factor of about 46 times was obtained, using NVIDIA's Geforce GTX470 GPU card against a 2.50 GHz Intel Quad Q9300 CPU processor. Compared with the HYPRE-based code performing in parallel on an 8-core tower server, the speedup of about 2 still could be observed. More encouragingly, without any mathematical acceleration technology, the GPU implementation ran about 5 times faster than CITATION which was speeded up by using the SOR method and Chebyshev extrapolation technique. (authors)
Parameterized code SHARM-3D for radiative transfer over inhomogeneous surfaces
NASA Astrophysics Data System (ADS)
Lyapustin, Alexei; Wang, Yujie
2005-12-01
The code SHARM-3D, developed for fast and accurate simulations of the monochromatic radiance at the top of the atmosphere over spatially variable surfaces with Lambertian or anisotropic reflectance, is described. The atmosphere is assumed to be laterally uniform across the image and to consist of two layers with aerosols contained in the bottom layer. The SHARM-3D code performs simultaneous calculations for all specified incidence-view geometries and multiple wavelengths in one run. The numerical efficiency of the current version of code is close to its potential limit and is achieved by means of two innovations. The first is the development of a comprehensive precomputed lookup table of the three-dimensional atmospheric optical transfer function for various atmospheric conditions. The second is the use of a linear kernel model of the land surface bidirectional reflectance factor (BRF) in our algorithm that has led to a fully parameterized solution in terms of the surface BRF parameters. The code is also able to model inland lakes and rivers. The water pixels are described with the Nakajima-Tanaka BRF model of wind-roughened water surface with a Lambertian offset, which is designed to model approximately the reflectance of suspended matter and of a shallow lake or river bottom.
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.
Numerical Results of Earth's Core Accumulation 3-D Modelling
NASA Astrophysics Data System (ADS)
Khachay, Yurie; Anfilogov, Vsevolod
2013-04-01
For a long time as a most convenient had been the model of mega impact in which the early forming of the Earth's core and mantle had been the consequence of formed protoplanet collision with the body of Mercurial mass. But all dynamical models of the Earth's accumulation and the estimations after the Pb-Pb system, lead to the conclusion that the duration of the planet accumulation was about 1 milliard years. But isotopic results after the W-Hf system testify about a very early (5-10) million years, dividing of the geochemical reservoirs of the core and mantle. In [1,3] it is shown, that the account of energy dissipating by the decay of short living radioactive elements and first of all Al,it is sufficient for heating even small bodies with dimensions about (50-100) km up to the iron melting temperature and can be realized a principal new differentiation mechanism. The inner parts of the melted preplanets can join and they are mainly of iron content, but the cold silicate fragments return to the supply zone. Only after the increasing of the gravitational radius, the growing area of the future core can save also the silicate envelope fragments. All existing dynamical accumulation models are constructed by using a spherical-symmetrical model. Hence for understanding the further planet evolution it is significant to trace the origin and evolution of heterogeneities, which occur on the planet accumulation stage. In that paper we are modeling distributions of temperature, pressure, velocity of matter flowing in a block of 3D- spherical body with a growing radius. The boundary problem is solved by the finite-difference method for the system of equations, which include equations which describe the process of accumulation, the Safronov equation, the equation of impulse balance, equation Navier-Stocks, equation for above litho static pressure and heat conductivity in velocity-pressure variables using the Businesque approach. The numerical algorithm of the problem solution in
Video coding and transmission standards for 3D television — a survey
NASA Astrophysics Data System (ADS)
Buchowicz, A.
2013-03-01
The emerging 3D television systems require effective techniques for transmission and storage of data representing a 3-D scene. The 3-D scene representations based on multiple video sequences or multiple views plus depth maps are especially important since they can be processed with existing video technologies. The review of the video coding and transmission techniques is presented in this paper.
Numerical modeling of Tibetan Plateau formation: Thin-sheet versus fully 3D models
NASA Astrophysics Data System (ADS)
Lechmann, S. M.; Schmalholz, S. M.; Kaus, B. J. P.
2009-04-01
Knowledge about the tectonic evolution of the Tibetan Plateau is still incomplete and many open questions remain concerning the deformation style of the crustal thickening, causing the abnormally high elevation of the Tibetan Plateau. Different models have been suggested explaining the crustal thickening by (1) homogeneous, continuous deformation using thin-sheet models, (2) discrete movement along thrusts developing crustal wedges and (3) lateral crustal flow due to pressure gradients resulting from topography. Most existing models are not fully three-dimensional (3D) models (e.g. thin-sheet models) and assume a certain deformation style a priori, which makes it difficult to judge the applicability of such constrained models to the formation of the Tibetan Plateau. We present a comparison of deformation styles during continent indentation resulting from a fully 3D numerical model and a thin-sheet model. The rheology for both models is power-law. The 3D model consists of four layers representing a simplified lithosphere: strong upper crust, weak lower crust, strong upper mantle and weak lower mantle. From the effective viscosity distribution of the 3D model a vertically averaged effective viscosity is calculated and used for the thin-sheet model to make direct comparisons between the two models. Simulating indentation is achieved by assigning free slip at one lateral side of the model, and fixing two other sides. The boundary at which indentation is taking place, exhibits a tripartite velocity profile: Next to the free slip side a section with constant horizontal velocity is applied. The velocity then gradually decreases towards zero, applying a cosine-function. The last section of the indenting boundary next to the fixed side is also fixed. The 3D model additionally exhibits a free surface and a bottom boundary allowing free slip. The 3D code employs the finite element method with a mixed velocity-pressure formulation to simulate incompressible flow. A Lagrangian
Verification and Validation of the k-kL Turbulence Model in FUN3D and CFL3D Codes
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.; Carlson, Jan-Renee; Rumsey, Christopher L.
2015-01-01
The implementation of the k-kL turbulence model using multiple computational uid dy- namics (CFD) codes is reported herein. The k-kL model is a two-equation turbulence model based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's equation is well suited for term-by-term mod- eling and displays useful features compared to other two-equation models. An important di erence is that this formulation leads to the inclusion of higher-order velocity derivatives in the source terms of the scale equations. This can enhance the ability of the Reynolds- averaged Navier-Stokes (RANS) solvers to simulate unsteady ows. The present report documents the formulation of the model as implemented in the CFD codes Fun3D and CFL3D. Methodology, veri cation and validation examples are shown. Attached and sepa- rated ow cases are documented and compared with experimental data. The results show generally very good comparisons with canonical and experimental data, as well as matching results code-to-code. The results from this formulation are similar or better than results using the SST turbulence model.
Newly-Developed 3D GRMHD Code and its Application to Jet Formation
NASA Technical Reports Server (NTRS)
Mizuno, Y.; Nishikawa, K.-I.; Koide, S.; Hardee, P.; Fishman, G. J.
2006-01-01
We have developed a new three-dimensional general relativistic magnetohydrodynamic code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-interpolated constrained transport scheme is used to maintain a divergence-free magnetic field. We have performed various 1-dimensional test problems in both special and general relativity by using several reconstruction methods and found that the new 3D GRMHD code shows substantial improvements over our previous model. The . preliminary results show the jet formations from a geometrically thin accretion disk near a non-rotating and a rotating black hole. We will discuss the jet properties depended on the rotation of a black hole and the magnetic field strength.
3D-radiative transfer in terrestrial atmosphere: An efficient parallel numerical procedure
NASA Astrophysics Data System (ADS)
Bass, L. P.; Germogenova, T. A.; Nikolaeva, O. V.; Kokhanovsky, A. A.; Kuznetsov, V. S.
2003-04-01
Light propagation and scattering in terrestrial atmosphere is usually studied in the framework of the 1D radiative transfer theory [1]. However, in reality particles (e.g., ice crystals, solid and liquid aerosols, cloud droplets) are randomly distributed in 3D space. In particular, their concentrations vary both in vertical and horizontal directions. Therefore, 3D effects influence modern cloud and aerosol retrieval procedures, which are currently based on the 1D radiative transfer theory. It should be pointed out that the standard radiative transfer equation allows to study these more complex situations as well [2]. In recent year the parallel version of the 2D and 3D RADUGA code has been developed. This version is successfully used in gammas and neutrons transport problems [3]. Applications of this code to radiative transfer in atmosphere problems are contained in [4]. Possibilities of code RADUGA are presented in [5]. The RADUGA code system is an universal solver of radiative transfer problems for complicated models, including 2D and 3D aerosol and cloud fields with arbitrary scattering anisotropy, light absorption, inhomogeneous underlying surface and topography. Both delta type and distributed light sources can be accounted for in the framework of the algorithm developed. The accurate numerical procedure is based on the new discrete ordinate SWDD scheme [6]. The algorithm is specifically designed for parallel supercomputers. The version RADUGA 5.1(P) can run on MBC1000M [7] (768 processors with 10 Gb of hard disc memory for each processor). The peak productivity is equal 1 Tfl. Corresponding scalar version RADUGA 5.1 is working on PC. As a first example of application of the algorithm developed, we have studied the shadowing effects of clouds on neighboring cloudless atmosphere, depending on the cloud optical thickness, surface albedo, and illumination conditions. This is of importance for modern satellite aerosol retrieval algorithms development. [1] Sobolev
MOM3D method of moments code theory manual
NASA Astrophysics Data System (ADS)
Shaeffer, John F.
1992-03-01
MOM3D is a FORTRAN algorithm that solves Maxwell's equations as expressed via the electric field integral equation for the electromagnetic response of open or closed three dimensional surfaces modeled with triangle patches. Two joined triangles (couples) form the vector current unknowns for the surface. Boundary conditions are for perfectly conducting or resistive surfaces. The impedance matrix represents the fundamental electromagnetic interaction of the body with itself. A variety of electromagnetic analysis options are possible once the impedance matrix is computed including backscatter radar cross section (RCS), bistatic RCS, antenna pattern prediction for user specified body voltage excitation ports, RCS image projection showing RCS scattering center locations, surface currents excited on the body as induced by specified plane wave excitation, and near field computation for the electric field on or near the body.
MOM3D method of moments code theory manual
NASA Technical Reports Server (NTRS)
Shaeffer, John F.
1992-01-01
MOM3D is a FORTRAN algorithm that solves Maxwell's equations as expressed via the electric field integral equation for the electromagnetic response of open or closed three dimensional surfaces modeled with triangle patches. Two joined triangles (couples) form the vector current unknowns for the surface. Boundary conditions are for perfectly conducting or resistive surfaces. The impedance matrix represents the fundamental electromagnetic interaction of the body with itself. A variety of electromagnetic analysis options are possible once the impedance matrix is computed including backscatter radar cross section (RCS), bistatic RCS, antenna pattern prediction for user specified body voltage excitation ports, RCS image projection showing RCS scattering center locations, surface currents excited on the body as induced by specified plane wave excitation, and near field computation for the electric field on or near the body.
3D unstructured-mesh radiation transport codes
Morel, J.
1997-12-31
Three unstructured-mesh radiation transport codes are currently being developed at Los Alamos National Laboratory. The first code is ATTILA, which uses an unstructured tetrahedral mesh in conjunction with standard Sn (discrete-ordinates) angular discretization, standard multigroup energy discretization, and linear-discontinuous spatial differencing. ATTILA solves the standard first-order form of the transport equation using source iteration in conjunction with diffusion-synthetic acceleration of the within-group source iterations. DANTE is designed to run primarily on workstations. The second code is DANTE, which uses a hybrid finite-element mesh consisting of arbitrary combinations of hexahedra, wedges, pyramids, and tetrahedra. DANTE solves several second-order self-adjoint forms of the transport equation including the even-parity equation, the odd-parity equation, and a new equation called the self-adjoint angular flux equation. DANTE also offers three angular discretization options: $S{_}n$ (discrete-ordinates), $P{_}n$ (spherical harmonics), and $SP{_}n$ (simplified spherical harmonics). DANTE is designed to run primarily on massively parallel message-passing machines, such as the ASCI-Blue machines at LANL and LLNL. The third code is PERICLES, which uses the same hybrid finite-element mesh as DANTE, but solves the standard first-order form of the transport equation rather than a second-order self-adjoint form. DANTE uses a standard $S{_}n$ discretization in angle in conjunction with trilinear-discontinuous spatial differencing, and diffusion-synthetic acceleration of the within-group source iterations. PERICLES was initially designed to run on workstations, but a version for massively parallel message-passing machines will be built. The three codes will be described in detail and computational results will be presented.
A numerical solution of 3D inviscid rotational flow in turbines and ducts
NASA Astrophysics Data System (ADS)
Oktay, Erdal; Akmandor, Sinan; Üçer, Ahmet
1998-04-01
The numerical solutions of inviscid rotational (Euler) flows were obtained using an explicit hexahedral unstructured cell vertex finite volume method. A second-order-accurate, one-step Lax-Wendroff scheme was used to solve the unsteady governing equations discretized in conservative form. The transonic circular bump, in which the location and the strength of the captured shock are well predicted, was used as the first test case. The nozzle guide vanes of the VKI low-speed turbine facility were used to validate the Euler code in highly 3D environment. Despite the high turning and the secondary flows which develop, close agreements have been obtained with experimental and numerical results associated with these test cases.
Filice, Luigino; Gagliardi, Francesco; Umbrello, Domenico; Shivpuri, Rajiv
2007-05-17
Metallic foams represent one of the most exciting materials introduced in the manufacturing scenario in the last years. In the study here addressed, the experimental and numerical investigations on the forging process of a simple foam billet shaped into complex sculptured parts were carried out. In particular, the deformation behavior of metallic foams and the development of density gradients were investigated through a series of experimental forging tests in order to produce a selected portion of a hip prosthesis. The human bone replacement was chosen as case study due to its industrial demand and for its particular 3D complex shape. A finite element code (Deform 3D) was utilized for modeling the foam behavior during the forging process and an accurate material rheology description was used based on a porous material model which includes the measured local density. Once the effectiveness of the utilized Finite Element model was verified through the comparison with the experimental evidences, a numerical study of the influence of the foam density was investigated. The obtained numerical results shown as the initial billet density plays an important role on the prediction of the final shape, the optimization of the flash as well as the estimation of the punch load.
Spatial parallelism of a 3D finite difference, velocity-stress elastic wave propagation code
Minkoff, S.E.
1999-12-01
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately, finite difference simulations for 3D elastic wave propagation are expensive. The authors model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MPI library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speedup. Because I/O is handled largely outside of the time-step loop (the most expensive part of the simulation) the authors have opted for straight-forward broadcast and reduce operations to handle I/O. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ghost cells. When this communication is balanced against computation by allocating subdomains of reasonable size, they observe excellent scaled speedup. Allocating subdomains of size 25 x 25 x 25 on each node, they achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
Spatial Parallelism of a 3D Finite Difference, Velocity-Stress Elastic Wave Propagation Code
MINKOFF,SUSAN E.
1999-12-09
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately. finite difference simulations for 3D elastic wave propagation are expensive. We model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MP1 library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speed up. Because i/o is handled largely outside of the time-step loop (the most expensive part of the simulation) we have opted for straight-forward broadcast and reduce operations to handle i/o. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ''ghost cells''. When this communication is balanced against computation by allocating subdomains of reasonable size, we observe excellent scaled speed up. Allocating subdomains of size 25 x 25 x 25 on each node, we achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
Implementation of wall boundary conditions for transpiration in F3D thin-layer Navier-Stokes code
NASA Technical Reports Server (NTRS)
Kandula, M.; Martin, F. W., Jr.
1991-01-01
Numerical boundary conditions for mass injection/suction at the wall are incorporated in the thin-layer Navier-Stokes code, F3D. The accuracy of the boundary conditions and the code is assessed by a detailed comparison of the predictions of velocity distributions and skin-friction coefficients with exact similarity solutions for laminar flow over a flat plate with variable blowing/suction, and measurements for turbulent flow past a flat plate with uniform blowing. In laminar flow, F3D predictions for friction coefficient compare well with exact similarity solution with and without suction, but produces large errors at moderate-to-large values of blowing. A slight Mach number dependence of skin-friction coefficient due to blowing in turbulent flow is computed by F3D code. Predicted surface pressures for turbulent flow past an airfoil with mass injection are in qualitative agreement with measurements for a flat plate.
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
3D visualization for the MARS14 Code
Rzepecki, Jaroslaw P.; Kostin, Mikhail A; Mokhov, Nikolai V.
2003-01-23
A new three-dimensional visualization engine has been developed for the MARS14 code system. It is based on the OPENINVENTOR graphics library and integrated with the MARS built-in two-dimensional Graphical-User Interface, MARS-GUI-SLICE. The integrated package allows thorough checking of complex geometry systems and their fragments, materials, magnetic fields, particle tracks along with a visualization of calculated 2-D histograms. The algorithms and their optimization are described for two geometry classes along with examples in accelerator and detector applications.
Numerical MHD codes for modeling astrophysical flows
NASA Astrophysics Data System (ADS)
Koldoba, A. V.; Ustyugova, G. V.; Lii, P. S.; Comins, M. L.; Dyda, S.; Romanova, M. M.; Lovelace, R. V. E.
2016-05-01
We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar coordinates, and fully 3D cylindrical coordinates. Viscosity and diffusivity are implemented in the code to control the accretion rate in the disk and the rate of penetration of the disk matter through the magnetic field lines. The code has been utilized for the numerical investigations of a number of different astrophysical problems, several examples of which are shown.
Code verification for unsteady 3-D fluid-solid interaction problems
NASA Astrophysics Data System (ADS)
Yu, Kintak Raymond; Étienne, Stéphane; Hay, Alexander; Pelletier, Dominique
2015-12-01
This paper describes a procedure to synthesize Manufactured Solutions for Code Verification of an important class of Fluid-Structure Interaction (FSI) problems whose behaviors can be modeled as rigid body vibrations in incompressible fluids. We refer this class of FSI problems as Fluid-Solid Interaction problems, which can be found in many practical engineering applications. The methodology can be utilized to develop Manufactured Solutions for both 2-D and 3-D cases. We demonstrate the procedure with our numerical code. We present details of the formulation and methodology. We also provide the reasonings behind our proposed approach. Results from grid and time step refinement studies confirm the verification of our solver and demonstrate the versatility of the simple synthesis procedure. In addition, the results also demonstrate that the modified decoupled approach to verify flow problems with high-order time-stepping schemes can be employed equally well to verify code for multi-physics problems (here, those of the Fluid-Solid Interaction) when the numerical discretization is based on the Method of Lines.
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.
Extending ALE3D, an Arbitrarily Connected hexahedral 3D Code, to Very Large Problem Size (U)
Nichols, A L
2010-12-15
As the number of compute units increases on the ASC computers, the prospect of running previously unimaginably large problems is becoming a reality. In an arbitrarily connected 3D finite element code, like ALE3D, one must provide a unique identification number for every node, element, face, and edge. This is required for a number of reasons, including defining the global connectivity array required for domain decomposition, identifying appropriate communication patterns after domain decomposition, and determining the appropriate load locations for implicit solvers, for example. In most codes, the unique identification number is defined as a 32-bit integer. Thus the maximum value available is 231, or roughly 2.1 billion. For a 3D geometry consisting of arbitrarily connected hexahedral elements, there are approximately 3 faces for every element, and 3 edges for every node. Since the nodes and faces need id numbers, using 32-bit integers puts a hard limit on the number of elements in a problem at roughly 700 million. The first solution to this problem would be to replace 32-bit signed integers with 32-bit unsigned integers. This would increase the maximum size of a problem by a factor of 2. This provides some head room, but almost certainly not one that will last long. Another solution would be to replace all 32-bit int declarations with 64-bit long long declarations. (long is either a 32-bit or a 64-bit integer, depending on the OS). The problem with this approach is that there are only a few arrays that actually need to extended size, and thus this would increase the size of the problem unnecessarily. In a future computing environment where CPUs are abundant but memory relatively scarce, this is probably the wrong approach. Based on these considerations, we have chosen to replace only the global identifiers with the appropriate 64-bit integer. The problem with this approach is finding all the places where data that is specified as a 32-bit integer needs to be
3D Direct Simulation Monte Carlo Code Which Solves for Geometrics
Energy Science and Technology Software Center (ESTSC)
1998-01-13
Pegasus is a 3D Direct Simulation Monte Carlo Code which solves for geometries which can be represented by bodies of revolution. Included are all the surface chemistry enhancements in the 2D code Icarus as well as a real vacuum pump model. The code includes multiple species transport.
PEGASUS. 3D Direct Simulation Monte Carlo Code Which Solves for Geometrics
Bartel, T.J.
1998-12-01
Pegasus is a 3D Direct Simulation Monte Carlo Code which solves for geometries which can be represented by bodies of revolution. Included are all the surface chemistry enhancements in the 2D code Icarus as well as a real vacuum pump model. The code includes multiple species transport.
Wall-touching kink mode calculations with the M3D code
Breslau, J. A. Bhattacharjee, A.
2015-06-15
This paper seeks to address a controversy regarding the applicability of the 3D nonlinear extended MHD code M3D [W. Park et al., Phys. Plasmas 6, 1796 (1999)] and similar codes to calculations of the electromagnetic interaction of a disrupting tokamak plasma with the surrounding vessel structures. M3D is applied to a simple test problem involving an external kink mode in an ideal cylindrical plasma, used also by the Disruption Simulation Code (DSC) as a model case for illustrating the nature of transient vessel currents during a major disruption. While comparison of the results with those of the DSC is complicated by effects arising from the higher dimensionality and complexity of M3D, we verify that M3D is capable of reproducing both the correct saturation behavior of the free boundary kink and the “Hiro” currents arising when the kink interacts with a conducting tile surface interior to the ideal wall.
Wall-touching kink mode calculations with the M3D code
NASA Astrophysics Data System (ADS)
Breslau, J. A.; Bhattacharjee, A.
2015-06-01
This paper seeks to address a controversy regarding the applicability of the 3D nonlinear extended MHD code M3D [W. Park et al., Phys. Plasmas 6, 1796 (1999)] and similar codes to calculations of the electromagnetic interaction of a disrupting tokamak plasma with the surrounding vessel structures. M3D is applied to a simple test problem involving an external kink mode in an ideal cylindrical plasma, used also by the Disruption Simulation Code (DSC) as a model case for illustrating the nature of transient vessel currents during a major disruption. While comparison of the results with those of the DSC is complicated by effects arising from the higher dimensionality and complexity of M3D, we verify that M3D is capable of reproducing both the correct saturation behavior of the free boundary kink and the "Hiro" currents arising when the kink interacts with a conducting tile surface interior to the ideal wall.
EM modeling for GPIR using 3D FDTD modeling codes
Nelson, S.D.
1994-10-01
An analysis of the one-, two-, and three-dimensional electrical characteristics of structural cement and concrete is presented. This work connects experimental efforts in characterizing cement and concrete in the frequency and time domains with the Finite Difference Time Domain (FDTD) modeling efforts of these substances. These efforts include Electromagnetic (EM) modeling of simple lossless homogeneous materials with aggregate and targets and the modeling dispersive and lossy materials with aggregate and complex target geometries for Ground Penetrating Imaging Radar (GPIR). Two- and three-dimensional FDTD codes (developed at LLNL) where used for the modeling efforts. Purpose of the experimental and modeling efforts is to gain knowledge about the electrical properties of concrete typically used in the construction industry for bridges and other load bearing structures. The goal is to optimize the performance of a high-sample-rate impulse radar and data acquisition system and to design an antenna system to match the characteristics of this material. Results show agreement to within 2 dB of the amplitudes of the experimental and modeled data while the frequency peaks correlate to within 10% the differences being due to the unknown exact nature of the aggregate placement.
Numerical simulation of 3-D Benard convection with gravitational modulation
NASA Technical Reports Server (NTRS)
Biringen, S.; Peltier, L. J.
1990-01-01
In this numerical study, randomly and sinusoidally modulated gravitational fields imposed on three-dimensional Rayleigh-Benard convection are investigated in an effort to understand the effects of vibration (G-Jitter) on fluid systems. The time-dependent, Navier-Stokes equations and the energy equation with Boussinesq approximations are solved by a semi-implicit, pseudospectral procedure. An analysis of energy balances indicates that with increasing modulation amplitude, transition from synchronous to relaxation oscillation goes through the subharmonic response. Random modulations are found to be less stabilizing than sinusoidal and are shown to impose three-dimensionality on the flow for some parameter ranges both at terrestrial and zero base gravity conditions.
3D numerical simulations of vesicle and inextensible capsule dynamics
NASA Astrophysics Data System (ADS)
Farutin, Alexander; Biben, Thierry; Misbah, Chaouqi
2014-10-01
Vesicles are locally-inextensible fluid membranes, capsules are endowed with in-plane shear elasticity mimicking the cytoskeleton of red blood cells (RBCs), but are extensible, while RBCs are inextensible. We use boundary integral (BI) methods based on the Green function techniques to model and solve numerically their dynamics. We regularize the single layer integral by subtraction of exact identities for the terms involving the normal and the tangential components of the force. The stability and precision of BI calculation is enhanced by taking advantage of additional quadrature nodes located in vertices of an auxiliary mesh, constructed by a standard refinement procedure from the main mesh. We extend the partition of unity technique to boundary integral calculation on triangular meshes. The proposed algorithm offers the same treatment of near-singular integration regardless whether the source and the target points belong to the same surface or not. Bending forces are calculated by using expressions derived from differential geometry. Membrane incompressibility is handled by using two penalization parameters per suspended entity: one for deviation of the global area from prescribed value and another for the sum of squares of local strains defined on each vertex. Extensible or inextensible capsules, a model of RBC, are studied by storing the position in the reference configuration for each vertex. The elastic force is then calculated by direct variation of the elastic energy. Various nonequilibrium physical examples on vesicles and capsules will be presented and the convergence and precision tests highlighted. Overall, a good convergence is observed with numerical error inversely proportional to the number of vertices used for surface discretization, the highest order of convergence allowed by piece-wise linear interpolation of the surface.
3-D numerical evaluation of density effects on tracer tests.
Beinhorn, M; Dietrich, P; Kolditz, O
2005-12-01
In this paper we present numerical simulations carried out to assess the importance of density-dependent flow on tracer plume development. The scenario considered in the study is characterized by a short-term tracer injection phase into a fully penetrating well and a natural hydraulic gradient. The scenario is thought to be typical for tracer tests conducted in the field. Using a reference case as a starting point, different model parameters were changed in order to determine their importance to density effects. The study is based on a three-dimensional model domain. Results were interpreted using concentration contours and a first moment analysis. Tracer injections of 0.036 kg per meter of saturated aquifer thickness do not cause significant density effects assuming hydraulic gradients of at least 0.1%. Higher tracer input masses, as used for geoelectrical investigations, may lead to buoyancy-induced flow in the early phase of a tracer test which in turn impacts further plume development. This also holds true for shallow aquifers. Results of simulations with different tracer injection rates and durations imply that the tracer input scenario has a negligible effect on density flow. Employing model cases with different realizations of a log conductivity random field, it could be shown that small variations of hydraulic conductivity in the vicinity of the tracer injection well have a major control on the local tracer distribution but do not mask effects of buoyancy-induced flow. PMID:16183165
Numerical homogenization for seismic wave propagation in 3D geological media
NASA Astrophysics Data System (ADS)
Cupillard, P.; Capdeville, Y.; Botella, A.
2014-12-01
Despite the important increase of the computational power in the last decades, simulating the seismic wave propagation through realistic geological models is still a challenge. By realistic models we here mean 3D media in which a broad variety (in terms of amplitude and extent) of heterogeneities lies, including discontinuities with complex geometry such as faulted and folded horizons, intrusive geological contacts and fault systems. To perform accurate numerical simulations, these discontinuities require complicated meshes which usually contain extremely small elements, yielding large, sometimes prohibitive, computation costs. Fortunately, the recent development of the non-periodic homogenization technique now enables to overcome this problem by computing smooth equivalent models for which a coarse mesh is sufficient to get an accurate wavefield. In this work, we present an efficient implementation of the technique which now allows for the homogenization of large 3D geological models. This implementation relies on a tetrahedral finite-element solution of the elasto-static equation behind the homogenization problem. Because this equation is time-independent, solving it is numerically cheaper than solving the wave equation, but it nevertheless requires some care because of the large size of the stiffness matrix arising from the fine mesh of realistic geological structures. A domain decomposition is therefore adopted. In our strategy, the obtained sub-domains overlap but they are independent so the solution within each of them can be computed either in series or in parallel. In addition, well-balanced loads, efficient search algorithms and multithreading are implemented to speed up the computation. The resulting code enables the homogenization of 3D elastic media in a time that is neglectable with respect to the simulation time of the wave propagation within. This is illustrated through a sub-surface model of the Furfooz karstic region, Belgium.
Effect of Frictions on the Ballistic Performance of a 3D Warp Interlock Fabric: Numerical Analysis
NASA Astrophysics Data System (ADS)
Ha-Minh, Cuong; Boussu, François; Kanit, Toufik; Crépin, David; Imad, Abdellatif
2012-06-01
3D interlock woven fabrics are promising materials to replace the 2D structures in the field of ballistic protection. The structural complexity of this material caused many difficulties in numerical modeling. This paper presents a new tool that permits to generate a geometry model of any woven fabric, then, mesh this model in shell or solid elements, and apply the mechanical properties of yarns to them. The tool shows many advantages over existing software. It is very handy in use with an organization of the functions in menu and using a graphic interface. It can describe correctly the geometry of all textile woven fabrics. With this tool, the orientation of the local axes of finite elements following the yarn direction facilitates defining the yarn mechanical properties in a numerical model. This tool can be largely applied because it is compatible with popular finite element codes such as Abaqus, Ansys, Radioss etc. Thanks to this tool, a finite element model was carried out to describe a ballistic impact on a 3D warp interlock Kevlar KM2® fabric. This work focuses on studying the effect of friction onto the ballistic impact behavior of this textile interlock structure. Results showed that the friction among yarns affects considerably on the impact behavior of this fabric. The effect of the friction between projectile and yarn is less important. The friction plays an important role in keeping the fabric structural stability during the impact event. This phenomenon explained why the projectile is easier to penetrate this 3D warp interlock fabric in the no-friction case. This result also indicates that the ballistic performance of the interlock woven fabrics can be improved by using fibers with great friction coefficients.
GATOR: A 3-D time-dependent simulation code for helix TWTs
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.
Standards-based approaches to 3D and multiview video coding
NASA Astrophysics Data System (ADS)
Sullivan, Gary J.
2009-08-01
The extension of video applications to enable 3D perception, which typically is considered to include a stereo viewing experience, is emerging as a mass market phenomenon, as is evident from the recent prevalence of 3D major cinema title releases. For high quality 3D video to become a commonplace user experience beyond limited cinema distribution, adoption of an interoperable coded 3D digital video format will be needed. Stereo-view video can also be studied as a special case of the more general technologies of multiview and "free-viewpoint" video systems. The history of standardization work on this topic is actually richer than people may typically realize. The ISO/IEC Moving Picture Experts Group (MPEG), in particular, has been developing interoperability standards to specify various such coding schemes since the advent of digital video as we know it. More recently, the ITU-T Visual Coding Experts Group (VCEG) has been involved as well in the Joint Video Team (JVT) work on development of 3D features for H.264/14496-10 Advanced Video Coding, including Multiview Video Coding (MVC) extensions. This paper surveys the prior, ongoing, and anticipated future standardization efforts on this subject to provide an overview and historical perspective on feasible approaches to 3D and multiview video coding.
Terascale direct numerical simulations of turbulent combustion using S3D.
Sankaran, Ramanan; Mellor-Crummy, J.; DeVries, M.; Yoo, Chun Sang; Ma, K. L.; Podhorski, N.; Liao, W. K.; Klasky, S.; de Supinski, B.; Choudhary, A.; Hawkes, Evatt R.; Chen, Jacqueline H.; Shende, Sameer
2008-08-01
Computational science is paramount to the understanding of underlying processes in internal combustion engines of the future that will utilize non-petroleum-based alternative fuels, including carbon-neutral biofuels, and burn in new combustion regimes that will attain high efficiency while minimizing emissions of particulates and nitrogen oxides. Next-generation engines will likely operate at higher pressures, with greater amounts of dilution and utilize alternative fuels that exhibit a wide range of chemical and physical properties. Therefore, there is a significant role for high-fidelity simulations, direct numerical simulations (DNS), specifically designed to capture key turbulence-chemistry interactions in these relatively uncharted combustion regimes, and in particular, that can discriminate the effects of differences in fuel properties. In DNS, all of the relevant turbulence and flame scales are resolved numerically using high-order accurate numerical algorithms. As a consequence terascale DNS are computationally intensive, require massive amounts of computing power and generate tens of terabytes of data. Recent results from terascale DNS of turbulent flames are presented here, illustrating its role in elucidating flame stabilization mechanisms in a lifted turbulent hydrogen/air jet flame in a hot air co-flow, and the flame structure of a fuel-lean turbulent premixed jet flame. Computing at this scale requires close collaborations between computer and combustion scientists to provide optimized scaleable algorithms and software for terascale simulations, efficient collective parallel I/O, tools for volume visualization of multiscale, multivariate data and automating the combustion workflow. The enabling computer science, applied to combustion science, is also required in many other terascale physics and engineering simulations. In particular, performance monitoring is used to identify the performance of key kernels in the DNS code, S3D and especially memory
Terascale direct numerical simulations of turbulent combustion using S3D
NASA Astrophysics Data System (ADS)
Chen, J. H.; Choudhary, A.; de Supinski, B.; DeVries, M.; Hawkes, E. R.; Klasky, S.; Liao, W. K.; Ma, K. L.; Mellor-Crummey, J.; Podhorszki, N.; Sankaran, R.; Shende, S.; Yoo, C. S.
2009-01-01
Computational science is paramount to the understanding of underlying processes in internal combustion engines of the future that will utilize non-petroleum-based alternative fuels, including carbon-neutral biofuels, and burn in new combustion regimes that will attain high efficiency while minimizing emissions of particulates and nitrogen oxides. Next-generation engines will likely operate at higher pressures, with greater amounts of dilution and utilize alternative fuels that exhibit a wide range of chemical and physical properties. Therefore, there is a significant role for high-fidelity simulations, direct numerical simulations (DNS), specifically designed to capture key turbulence-chemistry interactions in these relatively uncharted combustion regimes, and in particular, that can discriminate the effects of differences in fuel properties. In DNS, all of the relevant turbulence and flame scales are resolved numerically using high-order accurate numerical algorithms. As a consequence terascale DNS are computationally intensive, require massive amounts of computing power and generate tens of terabytes of data. Recent results from terascale DNS of turbulent flames are presented here, illustrating its role in elucidating flame stabilization mechanisms in a lifted turbulent hydrogen/air jet flame in a hot air coflow, and the flame structure of a fuel-lean turbulent premixed jet flame. Computing at this scale requires close collaborations between computer and combustion scientists to provide optimized scaleable algorithms and software for terascale simulations, efficient collective parallel I/O, tools for volume visualization of multiscale, multivariate data and automating the combustion workflow. The enabling computer science, applied to combustion science, is also required in many other terascale physics and engineering simulations. In particular, performance monitoring is used to identify the performance of key kernels in the DNS code, S3D and especially memory
3D numerical investigation on landslide generated tsunamis around a conical island
NASA Astrophysics Data System (ADS)
Montagna, Francesca; Bellotti, Giorgio
2010-05-01
This paper presents numerical computations of tsunamis generated by subaerial and submerged landslides falling along the flank of a conical island. The study is inspired by the tsunamis that on 30th December 2002 attacked the coast of the volcanic island of Stromboli (South Tyrrhenian sea, Italy). In particular this paper analyzes the important feature of the lateral spreading of landside generated tsunamis and the associated flooding hazard. The numerical model used in this study is the full three dimensional commercial code FLOW-3D. The model has already been successfully used (Choi et al., 2007; 2008; Chopakatla et al, 2008) to study the interaction of waves and structures. In the simulations carried out in this work a particular feature of the code has been employed: the GMO (General Moving Object) algorithm. It allows to reproduce the interaction between moving objects, as a landslide, and the water. FLOW-3D has been firstly validated using available 3D experiments reproducing tsunamis generated by landslides at the flank of a conical island. The experiments have been carried out in the LIC laboratory of the Polytechnic of Bari, Italy (Di Risio et al., 2009). Numerical and experimental time series of run-up and sea level recorded at gauges located at the flanks of the island and offshore have been successfully compared. This analysis shows that the model can accurately represent the generation, the propagation and the inundation of landslide generated tsunamis and suggests the use of the numerical model as a tool for preparing inundation maps. At the conference we will present the validation of the model and parametric analyses aimed to investigate how wave properties depend on the landslide kinematic and on further parameters such as the landslide volume and shape, as well as the radius of the island. The expected final results of the research are precomputed inundation maps that depend on the characteristics of the landslide and of the island. Finally we
NASA Astrophysics Data System (ADS)
Skála, J.; Baruffa, F.; Büchner, J.; Rampp, M.
2015-08-01
Context. The numerical simulation of turbulence and flows in almost ideal astrophysical plasmas with large Reynolds numbers motivates the implementation of magnetohydrodynamical (MHD) computer codes with low resistivity. They need to be computationally efficient and scale well with large numbers of CPU cores, allow obtaining a high grid resolution over large simulation domains, and be easily and modularly extensible, for instance, to new initial and boundary conditions. Aims: Our aims are the implementation, optimization, and verification of a computationally efficient, highly scalable, and easily extensible low-dissipative MHD simulation code for the numerical investigation of the dynamics of astrophysical plasmas with large Reynolds numbers in three dimensions (3D). Methods: The new GOEMHD3 code discretizes the ideal part of the MHD equations using a fast and efficient leap-frog scheme that is second-order accurate in space and time and whose initial and boundary conditions can easily be modified. For the investigation of diffusive and dissipative processes the corresponding terms are discretized by a DuFort-Frankel scheme. To always fulfill the Courant-Friedrichs-Lewy stability criterion, the time step of the code is adapted dynamically. Numerically induced local oscillations are suppressed by explicit, externally controlled diffusion terms. Non-equidistant grids are implemented, which enhance the spatial resolution, where needed. GOEMHD3 is parallelized based on the hybrid MPI-OpenMP programing paradigm, adopting a standard two-dimensional domain-decomposition approach. Results: The ideal part of the equation solver is verified by performing numerical tests of the evolution of the well-understood Kelvin-Helmholtz instability and of Orszag-Tang vortices. The accuracy of solving the (resistive) induction equation is tested by simulating the decay of a cylindrical current column. Furthermore, we show that the computational performance of the code scales very
3-D field computation: The near-triumph of commerical codes
Turner, L.R.
1995-07-01
In recent years, more and more of those who design and analyze magnets and other devices are using commercial codes rather than developing their own. This paper considers the commercial codes and the features available with them. Other recent trends with 3-D field computation include parallel computation and visualization methods such as virtual reality systems.
NASA Astrophysics Data System (ADS)
Filice, Luigino; Gagliardi, Francesco; Shivpuri, Rajiv; Umbrello, Domenico
2007-05-01
Metallic foams represent one of the most exciting materials introduced in the manufacturing scenario in the last years. In the study here addressed, the experimental and numerical investigations on the forging process of a simple foam billet shaped into complex sculptured parts were carried out. In particular, the deformation behavior of metallic foams and the development of density gradients were investigated through a series of experimental forging tests in order to produce a selected portion of a hip prosthesis. The human bone replacement was chosen as case study due to its industrial demand and for its particular 3D complex shape. A finite element code (Deform 3D®) was utilized for modeling the foam behavior during the forging process and an accurate material rheology description was used based on a porous material model which includes the measured local density. Once the effectiveness of the utilized Finite Element model was verified through the comparison with the experimental evidences, a numerical study of the influence of the foam density was investigated. The obtained numerical results shown as the initial billet density plays an important role on the prediction of the final shape, the optimization of the flash as well as the estimation of the punch load.
NASA Astrophysics Data System (ADS)
Ceccuzzi, Silvio; Maggiora, Riccardo; Milanesio, Daniele; Mirizzi, Francesco; Panaccione, Luigi
2011-12-01
The present work compares and experimentally validates the results coming out from the following three Lower Hybrid (LH) coupling codes: Brambilla code (M. Brambilla), GRILL3D-U (Mikhail Irzak, A. F. Ioffe Physico-Technical Institute, Russia) and TOPLHA (Politecnico di Torino, Italy). The conventional grill antenna, operating in FTU in different scenarios, is used as benchmark. The validation with experimental data is carried out with respect to the average reflection coefficients at the input of a row of the grill, considering two different phasings between adjacent waveguides: -90 ° and -75 °. A comparison between calculated power spectra is also presented. Good agreement can be observed for all the simulated plasma profiles and waveguide phasings between experimental data and codes, in particular for the most recent numerical tools, namely GRILL3D-U and TOPLHA.
Development of the PARVMEC Code for Rapid Analysis of 3D MHD Equilibrium
NASA Astrophysics Data System (ADS)
Seal, Sudip; Hirshman, Steven; Cianciosa, Mark; Wingen, Andreas; Unterberg, Ezekiel; Wilcox, Robert; ORNL Collaboration
2015-11-01
The VMEC three-dimensional (3D) MHD equilibrium has been used extensively for designing stellarator experiments and analyzing experimental data in such strongly 3D systems. Recent applications of VMEC include 2D systems such as tokamaks (in particular, the D3D experiment), where application of very small (delB/B ~ 10-3) 3D resonant magnetic field perturbations render the underlying assumption of axisymmetry invalid. In order to facilitate the rapid analysis of such equilibria (for example, for reconstruction purposes), we have undertaken the task of parallelizing the VMEC code (PARVMEC) to produce a scalable and temporally rapidly convergent equilibrium code for use on parallel distributed memory platforms. The parallelization task naturally splits into three distinct parts 1) radial surfaces in the fixed-boundary part of the calculation; 2) two 2D angular meshes needed to compute the Green's function integrals over the plasma boundary for the free-boundary part of the code; and 3) block tridiagonal matrix needed to compute the full (3D) pre-conditioner near the final equilibrium state. Preliminary results show that scalability is achieved for tasks 1 and 3, with task 2 still nearing completion. The impact of this work on the rapid reconstruction of D3D plasmas using PARVMEC in the V3FIT code will be discussed. Work supported by U.S. DOE under Contract DE-AC05-00OR22725 with UT-Battelle, LLC.
3D video coding: an overview of present and upcoming standards
NASA Astrophysics Data System (ADS)
Merkle, Philipp; Müller, Karsten; Wiegand, Thomas
2010-07-01
An overview of existing and upcoming 3D video coding standards is given. Various different 3D video formats are available, each with individual pros and cons. The 3D video formats can be separated into two classes: video-only formats (such as stereo and multiview video) and depth-enhanced formats (such as video plus depth and multiview video plus depth). Since all these formats exist of at least two video sequences and possibly additional depth data, efficient compression is essential for the success of 3D video applications and technologies. For the video-only formats the H.264 family of coding standards already provides efficient and widely established compression algorithms: H.264/AVC simulcast, H.264/AVC stereo SEI message, and H.264/MVC. For the depth-enhanced formats standardized coding algorithms are currently being developed. New and specially adapted coding approaches are necessary, as the depth or disparity information included in these formats has significantly different characteristics than video and is not displayed directly, but used for rendering. Motivated by evolving market needs, MPEG has started an activity to develop a generic 3D video standard within the 3DVC ad-hoc group. Key features of the standard are efficient and flexible compression of depth-enhanced 3D video representations and decoupling of content creation and display requirements.
Optimising GPR modelling: A practical, multi-threaded approach to 3D FDTD numerical modelling
NASA Astrophysics Data System (ADS)
Millington, T. M.; Cassidy, N. J.
2010-09-01
The demand for advanced interpretational tools has lead to the development of highly sophisticated, computationally demanding, 3D GPR processing and modelling techniques. Many of these methods solve very large problems with stepwise methods that utilise numerically similar functions within iterative computational loops. Problems of this nature are readily parallelised by splitting the computational domain into smaller, independent chunks for direct use on cluster-style, multi-processor supercomputers. Unfortunately, the implications of running such facilities, as well as time investment needed to develop the parallel codes, means that for most researchers, the use of these advanced methods is too impractical. In this paper, we propose an alternative method of parallelisation which exploits the capabilities of the modern multi-core processors (upon which today's desktop PCs are built) by multi-threading the calculation of a problem's individual sub-solutions. To illustrate the approach, we have applied it to an advanced, 3D, finite-difference time-domain (FDTD) GPR modelling tool in which the calculation of the individual vector field components is multi-threaded. To be of practical use, the FDTD scheme must be able to deliver accurate results with short execution times and we, therefore, show that the performance benefits of our approach can deliver runtimes less than half those of the more conventional, serial programming techniques. We evaluate implementations of the technique using different programming languages (e.g., Matlab, Java, C++), which will facilitate the construction of a flexible modelling tool for use in future GPR research. The implementations are compared on a variety of typical hardware platforms, having between one and eight processing cores available, and also a modern Graphical Processing Unit (GPU)-based computer. Our results show that a multi-threaded xyz modelling approach is easy to implement and delivers excellent results when implemented
INS3D: An incompressible Navier-Stokes code in generalized three-dimensional coordinates
NASA Technical Reports Server (NTRS)
Rogers, S. E.; Kwak, D.; Chang, J. L. C.
1987-01-01
The operation of the INS3D code, which computes steady-state solutions to the incompressible Navier-Stokes equations, is described. The flow solver utilizes a pseudocompressibility approach combined with an approximate factorization scheme. This manual describes key operating features to orient new users. This includes the organization of the code, description of the input parameters, description of each subroutine, and sample problems. Details for more extended operations, including possible code modifications, are given in the appendix.
Quantum self-correction in the 3D cubic code model.
Bravyi, Sergey; Haah, Jeongwan
2013-11-15
A big open question in the quantum information theory concerns the feasibility of a self-correcting quantum memory. A quantum state recorded in such memory can be stored reliably for a macroscopic time without need for active error correction, if the memory is in contact with a cold enough thermal bath. Here we report analytic and numerical evidence for self-correcting behavior in the quantum spin lattice model known as the 3D cubic code. We prove that its memory time is at least L(cβ), where L is the lattice size, β is the inverse temperature of the bath, and c>0 is a constant coefficient. However, this bound applies only if the lattice size L does not exceed a critical value which grows exponentially with β. In that sense, the model can be called a partially self-correcting memory. We also report a Monte Carlo simulation indicating that our analytic bounds on the memory time are tight up to constant coefficients. To model the readout step we introduce a new decoding algorithm, which can be implemented efficiently for any topological stabilizer code. A longer version of this work can be found in Bravyi and Haah, arXiv:1112.3252. PMID:24289671
Quantum Self-Correction in the 3D Cubic Code Model
NASA Astrophysics Data System (ADS)
Bravyi, Sergey; Haah, Jeongwan
2013-11-01
A big open question in the quantum information theory concerns the feasibility of a self-correcting quantum memory. A quantum state recorded in such memory can be stored reliably for a macroscopic time without need for active error correction, if the memory is in contact with a cold enough thermal bath. Here we report analytic and numerical evidence for self-correcting behavior in the quantum spin lattice model known as the 3D cubic code. We prove that its memory time is at least Lcβ, where L is the lattice size, β is the inverse temperature of the bath, and c>0 is a constant coefficient. However, this bound applies only if the lattice size L does not exceed a critical value which grows exponentially with β. In that sense, the model can be called a partially self-correcting memory. We also report a Monte Carlo simulation indicating that our analytic bounds on the memory time are tight up to constant coefficients. To model the readout step we introduce a new decoding algorithm, which can be implemented efficiently for any topological stabilizer code. A longer version of this work can be found in Bravyi and Haah, arXiv:1112.3252.
Status and future of the 3D MAFIA group of codes
NASA Astrophysics Data System (ADS)
Ebeling, F.; Klatt, R.; Krawzcyk, F.; Lawinsky, E.; Weiland, T.; Wipf, S. G.; Steffen, B.; Barts, T.; Browman, J.; Cooper, R. K.; Rodenz, G.
1988-12-01
The group of fully three dimensional computer codes for solving Maxwell's equations for a wide range of applications, MAFIA, is already well established. Extensive comparisons with measurements have demonstrated the accuracy of the computations. A large numer of components have been designed for accelerators, such as kicker magnets, non cyclindrical cavities, ferrite loaded cavities, vacuum chambers with slots and transitions, etc. The latest additions to the system include a new static solver that can calculate 3D magneto- and electrostatic fields, and a self consistent version of the 2D-BCI that solves the field equations and the equations of motion in parallel. Work on new eddy current modules has started, which will allow treatment of laminated and/or solid iron cores excited by low frequency currents. Based on our experience with the present releases 1 and 2, we have started a complete revision of the whole user interface and data structure, which will make the codes even more user-friendly and flexible.
Depth-based coding of MVD data for 3D video extension of H.264/AVC
NASA Astrophysics Data System (ADS)
Rusanovskyy, Dmytro; Hannuksela, Miska M.; Su, Wenyi
2013-06-01
This paper describes a novel approach of using depth information for advanced coding of associated video data in Multiview Video plus Depth (MVD)-based 3D video systems. As a possible implementation of this conception, we describe two coding tools that have been developed for H.264/AVC based 3D Video Codec as response to Moving Picture Experts Group (MPEG) Call for Proposals (CfP). These tools are Depth-based Motion Vector Prediction (DMVP) and Backward View Synthesis Prediction (BVSP). Simulation results conducted under JCT-3V/MPEG 3DV Common Test Conditions show, that proposed in this paper tools reduce bit rate of coded video data by 15% of average delta bit rate reduction, which results in 13% of bit rate savings on total for the MVD data over the state-of-the-art MVC+D coding. Moreover, presented in this paper conception of depth-based coding of video has been further developed by MPEG 3DV and JCT-3V and this work resulted in even higher compression efficiency, bringing about 20% of delta bit rate reduction on total for coded MVD data over the reference MVC+D coding. Considering significant gains, proposed in this paper coding approach can be beneficial for development of new 3D video coding standards. [Figure not available: see fulltext.
Reactor Dosimetry Applications Using RAPTOR-M3G:. a New Parallel 3-D Radiation Transport Code
NASA Astrophysics Data System (ADS)
Longoni, Gianluca; Anderson, Stanwood L.
2009-08-01
The numerical solution of the Linearized Boltzmann Equation (LBE) via the Discrete Ordinates method (SN) requires extensive computational resources for large 3-D neutron and gamma transport applications due to the concurrent discretization of the angular, spatial, and energy domains. This paper will discuss the development RAPTOR-M3G (RApid Parallel Transport Of Radiation - Multiple 3D Geometries), a new 3-D parallel radiation transport code, and its application to the calculation of ex-vessel neutron dosimetry responses in the cavity of a commercial 2-loop Pressurized Water Reactor (PWR). RAPTOR-M3G is based domain decomposition algorithms, where the spatial and angular domains are allocated and processed on multi-processor computer architectures. As compared to traditional single-processor applications, this approach reduces the computational load as well as the memory requirement per processor, yielding an efficient solution methodology for large 3-D problems. Measured neutron dosimetry responses in the reactor cavity air gap will be compared to the RAPTOR-M3G predictions. This paper is organized as follows: Section 1 discusses the RAPTOR-M3G methodology; Section 2 describes the 2-loop PWR model and the numerical results obtained. Section 3 addresses the parallel performance of the code, and Section 4 concludes this paper with final remarks and future work.
Insights from 3D numerical simulations on the dynamics of the India-Asia collision zone
NASA Astrophysics Data System (ADS)
Pusok, A. E.; Kaus, B.; Popov, A.
2013-12-01
The dynamics of the India-Asia collision zone remains one of the most remarkable topics of the current research interest: the transition from subduction to collision and uplift, followed by the rise of the abnormally thick Tibetan plateau, and the deformation at its Eastern and Western syntaxes, are processes still not fully understood. Models that have addressed this topic include wholescale underthrusting of Indian lithospheric mantle under Tibet, distributed homogeneous shortening or the thin-sheet model, slip-line field model for lateral extrusion or lower crustal flow models for the exhumation of the Himalayan units and lateral spreading of the Tibetan plateau. Of these, the thin-sheet model has successfully illustrated some of the basic physics of continental collision and has the advantage of a 3D model being reduced to 2D, but one of its major shortcomings is that it cannot simultaneously represent channel flow and gravitational collapse of the mantle lithosphere, since these mechanisms require the lithosphere to interact with the underlying mantle, or to have a vertically non-homogeneous rheology. As a consequence, 3D models are emerging as powerful tools to understand the dynamics of coupled systems. However, because of yet recent developments and various complexities, the current 3D models simulating the dynamics of continent collision zones have relied on certain explicit assumptions, such as replacing part of the asthenosphere with various types of boundary conditions that mimic the effect of mantle flow, in order to focus on the lithospheric/crustal deformation. Here, we employ the parallel 3D code LaMEM (Lithosphere and Mantle Evolution Model), with a finite difference staggered grid solver, which is capable of simulating lithospheric deformation while simultaneously taking mantle flow and a free surface into account. We present qualitative results on lithospheric and upper-mantle scale simulations in which the Indian lithosphere is subducted and
NASA Astrophysics Data System (ADS)
Crosta, G.; Imposimato, S.; Roddeman, D.; Frattini, P.
2012-04-01
Fast moving landslides can be originated along slopes in mountainous terrains with natural and artificial lakes, or fjords at the slope foot. This landslides can reach extremely high speed and the impact with the immobile reservoir water can be influenced by the local topography and the landslide mass profile. The impact can generate large impulse waves and landslide tsunami. Initiation, propagation and runup are the three phases that need to be considered. The landslide evolution and the consequent wave can be controlled by the initial mass position (subaerial, partially or completely submerged), the landslide speed, the type of material, the subaerial and subaqueous slope geometry, the landslide depth and length at the impact, and the water depth. Extreme events have been caused by subaerial landslides: the 1963 Vajont rockslide (Italy), the 1958 Lituya Bay event (Alaska), the Tafjord and the Loen multiple events event (Norway), also from volcanic collapses (Hawaii and Canary islands). Various researchers completed a systematic experimental work on 2D and 3D wave generation and propagation (Kamphuis and Bowering, 1970; Huber, 1980; Müller, 1995; Huber and Hager, 1997; Fritz, 2002; Zweifel, 2004; Panizzo et al., 2005; Heller, 2007; Heller and Kinnear, 2010; Sælevik et al., 2009), using both rigid blocks and deformable granular" masses. Model data and results have been used to calibrate and validate numerical modelling tools (Harbitz, 1992; Jiang and LeBlond, 1993; Grilli et al., 2002; Grilli and Watts, 2005; Lynett and Liu, 2005; Tinti et al., 2006; Abadie et al., 2010) generally considering simplified rheologies (e.g. viscous rheologies) for subaerial subaqueous spreading. We use a FEM code (Roddeman, 2011; Crosta et al., 2006, 2009, 2010, 2011) adopting an Eulerian-Lagrangian approach to give accurate results for large deformations. We model both 2D and fully 3D events considering different settings. The material is considered as a fully deformable elasto
User's manual for PELE3D: a computer code for three-dimensional incompressible fluid dynamics
McMaster, W H
1982-05-07
The PELE3D code is a three-dimensional semi-implicit Eulerian hydrodynamics computer program for the solution of incompressible fluid flow coupled to a structure. The fluid and coupling algorithms have been adapted from the previously developed two-dimensional code PELE-IC. The PELE3D code is written in both plane and cylindrical coordinates. The coupling algorithm is general enough to handle a variety of structural shapes. The free surface algorithm is able to accommodate a top surface and several independent bubbles. The code is in a developmental status since all the intended options have not been fully implemented and tested. Development of this code ended in 1980 upon termination of the contract with the Nuclear Regulatory Commission.
Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers
NASA Technical Reports Server (NTRS)
Liu, Chaoqun; Liu, Zhining
1993-01-01
A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.
An ALE Based FE Formulation for the 3D Numerical Simulation of Fineblanking Processes
NASA Astrophysics Data System (ADS)
Manopulo, Niko; Tong, Longchang; Hora, Pavel
2010-06-01
Fineblanking is a manufacturing process which allows the mass production of blanked products with superior surface quality. The 3D numerical simulation of this particularly precise process is however challenging. This is because quality-critical tool features such as the die clearance and the shape of the cutting edges have dimensions up to two orders of magnitude smaller than the average part dimensions. If conventional Updated Lagrange codes are used, a very high FE mesh resolution becomes a must in order to accurately represent the surface evolution along the edge, which in turn makes the computation unfeasible. The methodology presented in this paper makes use of the Arbitrary Lagrangian Eulerian FE Formulation in order to keep control over the mesh region in contact with the tools. This way an optimal FE mesh can be guaranteed throughout the computation. This not only reduces the computational cost considerably, but also avoids mesh distortion along the cutting edge, allowing an accurate representation of the tool features. This approach will be used in conjunction to the stress limit criterion delineated in order to predict material failure in fine blanked products. Numerical results will be validated against the experiments carried out with a specially designed fineblanking tool in use at our institute.
Numerical simulation of unsteady flow characteristics for cavitation around a 3-D hydrofoil
NASA Astrophysics Data System (ADS)
Ahn, S. H.; Xiao, Y. X.; Wang, Z. W.
2015-01-01
At present it is possible to predict more accurately by various numerical methods established for cavitation simulation around a hydrofoil. However, for the solution of the complex unsteady cavity flow, it is still marginal. In this paper, numerical method is adopted to simulate cavitation around 3-D NACA0015 hydrofoil with homogeneous two-phase flow calculation using commercial code CFX-solver with two turbulence models, the standard RNG k-epsilon turbulence model and the modified RNG k-epsilon turbulence model respectively. First, pressure coefficient for non-cavitating flow, time averaged values of unsteady cavity flow around a hydrofoil are verified to simulate more closely to an actual cavity flow. And then frequency analysis is performed with Fast Fourier Transform. The results show that the calculation results with modified RNG k-epsilon turbulence model agree with experimental results in terms of mean cavity length and pressure drop, but the unsteady flow characteristics of oscillating cavitation still deviate slightly in terms of unsteady cavity flow.
Numerical study of the 3-D effect on FEL performance and its application to the APS LEUTL FEL
Chae, Y.C.
1998-09-01
A Low-Energy Undulator Test Line (LEUTL) is under construction at the Advanced Photon Source (APS). In LEUTL periodic focusing is provided by external quadrupoles. This results in an elliptical beam with its betatron oscillation envelope varying along the undulators. The free-electron laser (FEL) interaction with such a beam will exhibit truly 3-D effects. Thus the investigation of 3-D effects is important in optimizing the FEL performance. The programs GINGER and TDA3D, coupled with theoretically known facts, have been used for this purpose. Both programs are fully 3-D in moving the particle, but model the interaction between particles and axially symmetric electromagnetic waves. Even though TDA3D can include a few azimuthal modes in the interaction, it is still not a fully 3-D FEL code. However, they show that these 2-D programs can still be used for an elliptical beam whose aspect ratio is within certain limits. The author presents numerical results of FEL performance for the circular beam, the elliptical beam, and finally for the beam in the realistic LEUTL lattice.
NASA Astrophysics Data System (ADS)
Miensopust, Marion P.; Queralt, Pilar; Jones, Alan G.; 3D MT modellers
2013-06-01
Over the last half decade the need for, and importance of, three-dimensional (3-D) modelling of magnetotelluric (MT) data have increased dramatically and various 3-D forward and inversion codes are in use and some have become commonly available. Comparison of forward responses and inversion results is an important step for code testing and validation prior to `production' use. The various codes use different mathematical approximations to the problem (finite differences, finite elements or integral equations), various orientations of the coordinate system, different sign conventions for the time dependence and various inversion strategies. Additionally, the obtained results are dependent on data analysis, selection and correction as well as on the chosen mesh, inversion parameters and regularization adopted, and therefore, a careful and knowledge-based use of the codes is essential. In 2008 and 2011, during two workshops at the Dublin Institute for Advanced Studies over 40 people from academia (scientists and students) and industry from around the world met to discuss 3-D MT inversion. These workshops brought together a mix of code writers as well as code users to assess the current status of 3-D modelling, to compare the results of different codes, and to discuss and think about future improvements and new aims in 3-D modelling. To test the numerical forward solutions, two 3-D models were designed to compare the responses obtained by different codes and/or users. Furthermore, inversion results of these two data sets and two additional data sets obtained from unknown models (secret models) were also compared. In this manuscript the test models and data sets are described (supplementary files are available) and comparisons of the results are shown. Details regarding the used data, forward and inversion parameters as well as computational power are summarized for each case, and the main discussion points of the workshops are reviewed. In general, the responses
Kou, Wen; Li, Decong; Qiao, Jiyan; Chen, Li; Ding, Yansheng; Sjögren, Göran
2011-02-01
In this study, a numerical approach to the fracture behavior in a three-unit zirconia-based fixed partial denture (FPD) framework was made under mechanical loading using a newly developed three-dimensional (3D) numerical modeling code. All the materials studied were treated heterogeneously and Weibull distribution law was applied to describe the heterogeneity. The Mohr-Coulomb failure criterion with tensile strength cut-off was utilized to judge whether the material was in an elastic or failed state. For validation, the fracture pattern obtained from the numerical modeling was compared with a laboratory test; they largely correlated with each other. Similar fracture initiation sites were detected both in the numerical simulation and in an earlier fractographic analysis. The numerical simulation applied in this study clearly described the stress distribution and fracture process of zirconia-based FPD frameworks, information that could not be gained from the laboratory tests alone. Thus, the newly developed 3D numerical modeling code seems to be an efficient tool for prediction of the fracture process in ceramic FPD frameworks. PMID:21210519
Description of FEL3D: A three dimensional simulation code for TOK and FEL
Dutt, S.; Friedman, A.; Gover, A.
1988-10-20
FEL3D is a three dimensional simulation code, written for the purpose of calculating the parameters of coherent radiation emitted by electrons in an undulator. The program was written predominantly for simulating the coherent super-radiant harmonic frequency emission of electrons which are being bunched by an external laser beam while propagating in an undulator magnet. This super-radiant emission is to be studied in the TOK (transverse optical klystron) experiment, which is under construction in the NSLS department at Brookhaven National Laboratory. The program can also calculate the stimulated emission radiometric properties of a free electron laser (FEL) taking into account three dimensional effects. While this application is presently limited to the small gain operation regime of FEL's, extension to the high gain regime is expected to be relatively easy. The code is based on a semi-analytical concept. Instead of a full numerical solution of the Maxwell-Lorentz equations, the trajectories of the electron in the wiggler field are calculated analytically, and the radiation fields are expanded in terms of free space eigen-modes. This approach permits efficient computation, with a computation time of about 0.1 sec/electron on the BNL IBM 3090. The code reflects the important three dimensional features of the electron beam, the modulating laser beam, and the emitted radiation field. The statistical approach is based on averaging over the electron initial conditions according to a given distribution function in phase space, rather than via Monte-Carlo simulation. The present version of the program is written for uniform periodic wiggler field, but extension to nonuniform fields is straightforward. 4 figs., 5 tabs.
Impact of packet losses in scalable 3D holoscopic video coding
NASA Astrophysics Data System (ADS)
Conti, Caroline; Nunes, Paulo; Ducla Soares, Luís.
2014-05-01
Holoscopic imaging became a prospective glassless 3D technology to provide more natural 3D viewing experiences to the end user. Additionally, holoscopic systems also allow new post-production degrees of freedom, such as controlling the plane of focus or the viewing angle presented to the user. However, to successfully introduce this technology into the consumer market, a display scalable coding approach is essential to achieve backward compatibility with legacy 2D and 3D displays. Moreover, to effectively transmit 3D holoscopic content over error-prone networks, e.g., wireless networks or the Internet, error resilience techniques are required to mitigate the impact of data impairments in the user quality perception. Therefore, it is essential to deeply understand the impact of packet losses in terms of decoding video quality for the specific case of 3D holoscopic content, notably when a scalable approach is used. In this context, this paper studies the impact of packet losses when using a three-layer display scalable 3D holoscopic video coding architecture previously proposed, where each layer represents a different level of display scalability (i.e., L0 - 2D, L1 - stereo or multiview, and L2 - full 3D holoscopic). For this, a simple error concealment algorithm is used, which makes use of inter-layer redundancy between multiview and 3D holoscopic content and the inherent correlation of the 3D holoscopic content to estimate lost data. Furthermore, a study of the influence of 2D views generation parameters used in lower layers on the performance of the used error concealment algorithm is also presented.
NASA Astrophysics Data System (ADS)
Fernandez, N.; Kaus, B. J. P.
2012-04-01
Many fold-and-thrust belts are dominated by crustal scale folding that exhibits fairly regular fold spacing. For example, the Fars region in the Zagros Mountains shows a fold spacing with a normal distribution around a dominant wavelength of 14 Km ± 3 Km, yet having a wide variability of aspect ratios (length to wavelength ratios; Yamato et al., 2011). To which extend this is consistent with a crustal-scale folding instability or how the regional spacing of folding can be used to constrain regional rheological parameters are not fully resolved questions. To get insights into these problems we have investigated the dominant wavelength selection and evolution in a true multilayer system (Schmid and Podlachikov, 2006) with three different viscosities: lower salt layer (ηs), and overlying weak layers (ηw) and competent layers (ηc). This has been done by means of two tools: a semi-analytical solution and numerical models. The 2D semi-analytical approach was applied to derive mechanical phase diagrams that can be used to distinguish different folding modes using two viscosity ratios (R1= ηc/ ηs and R2= ηc/ ηw). To test the validity of the phase diagrams beyond the initial stages of folding for which the analytical approach is valid, we performed several 3D high-resolution forward numerical runs using a finite element code (LaMEM). Additionally, irregular bottom topography was implemented in the numerical runs in order to account for variable salt thickness distribution and consequently study its effect on the wavelength selection. A straight but gradual salt thickness variation, sudden thickness variations due to a basement step or an arc shaped salt basin among other cases could be investigated. It was observed that the bottom topography exerts an impact on the velocity field of the different folding modes and as a result, its influence can be observed on the resulting topography. However, not all the folding modes exhibit an initial wavelength that is dependent
Multitasking the INS3D-LU code on the Cray Y-MP
NASA Technical Reports Server (NTRS)
Fatoohi, Rod; Yoon, Seokkwan
1991-01-01
This paper presents the results of multitasking the INS3D-LU code on eight processors. The code is a full Navier-Stokes solver for incompressible fluid in three dimensional generalized coordinates using a lower-upper symmetric-Gauss-Seidel implicit scheme. This code has been fully vectorized on oblique planes of sweep and parallelized using autotasking with some directives and minor modifications. The timing results for five grid sizes are presented and analyzed. The code has achieved a processing rate of over one Gflops.
RELAP5-3D Code for Supercritical-Pressure Light-Water-Cooled Reactors
Riemke, Richard Allan; Davis, Cliff Bybee; Schultz, Richard Raphael
2003-04-01
The RELAP5-3D computer program has been improved for analysis of supercritical-pressure, light-water-cooled reactors. Several code modifications were implemented to correct code execution failures. Changes were made to the steam table generation, steam table interpolation, metastable states, interfacial heat transfer coefficients, and transport properties (viscosity and thermal conductivity). The code modifications now allow the code to run slow transients above the critical pressure as well as blowdown transients (modified Edwards pipe and modified existing pressurized water reactor model) that pass near the critical point.
Description of a parallel, 3D, finite element, hydrodynamics-diffusion code
Milovich, J L; Prasad, M K; Shestakov, A I
1999-04-11
We describe a parallel, 3D, unstructured grid finite element, hydrodynamic diffusion code for inertial confinement fusion (ICF) applications and the ancillary software used to run it. The code system is divided into two entities, a controller and a stand-alone physics code. The code system may reside on different computers; the controller on the user's workstation and the physics code on a supercomputer. The physics code is composed of separate hydrodynamic, equation-of-state, laser energy deposition, heat conduction, and radiation transport packages and is parallelized for distributed memory architectures. For parallelization, a SPMD model is adopted; the domain is decomposed into a disjoint collection of subdomains, one per processing element (PE). The PEs communicate using MPI. The code is used to simulate the hydrodynamic implosion of a spherical bubble.
Validation of Numerical Codes to Compute Tsunami Runup And Inundation
NASA Astrophysics Data System (ADS)
Velioğlu, Deniz; Cevdet Yalçıner, Ahmet; Kian, Rozita; Zaytsev, Andrey
2015-04-01
FLOW 3D and NAMI DANCE are two numerical codes which can be applied to analysis of flow and motion of long waves. Flow 3D simulates linear and nonlinear propagating surface waves as well as irregular waves including long waves. NAMI DANCE uses finite difference computational method to solve nonlinear shallow water equations (NSWE) in long wave problems, specifically tsunamis. Both codes can be applied to tsunami simulations and visualization of long waves. Both codes are capable of solving flooding problems. However, FLOW 3D is designed mainly to solve flooding problem from land and NAMI DANCE is designed to solve flooding problem from the sea. These numerical codes are applied to some benchmark problems for validation and verification. One useful benchmark problem is the runup of solitary waves which is investigated analytically and experimentally by Synolakis (1987). Since 1970s, solitary waves have commonly been used to model tsunamis especially in experimental and numerical studies. In this respect, a benchmark problem on runup of solitary waves is a relevant choice to assess the capability and validity of the numerical codes on amplification of tsunamis. In this study both codes have been tested, compared and validated by applying to the analytical benchmark problem of solitary wave runup on a sloping beach. Comparison of the results showed that both codes are in good agreement with the analytical and experimental results and thus can be proposed to be used in inundation of long waves and tsunami hazard analysis.
Planet-Disk Interaction on the GPU: The FARGO3D code
NASA Astrophysics Data System (ADS)
Masset, F. S.; Benítez-Llambay, P.
2015-10-01
We present the new code FARGO3D. It is a finite difference code that solves the equations of hydrodynamics or magnetohydrodynamics on a Cartesian, cylindrical or spherical mesh. It features orbital advection, conserves mass and (angular) momentum to machine accuracy. Special emphasis is put on the description of planet disk tidal interactions. It is parallelized with MPI, and it can run indistinctly on CPUs or GPUs, without the need to program in a GPU oriented language.
Numerical Simulations of High-Frequency Respiratory Flows in 2D and 3D Lung Bifurcation Models
NASA Astrophysics Data System (ADS)
Chen, Zixi; Parameswaran, Shamini; Hu, Yingying; He, Zhaoming; Raj, Rishi; Parameswaran, Siva
2014-07-01
To better understand the human pulmonary system and optimize the high-frequency oscillatory ventilation (HFOV) design, numerical simulations were conducted under normal breathing frequency and HFOV condition using a CFD code Ansys Fluent and its user-defined C programs. 2D and 3D double bifurcating lung models were created, and the geometry corresponds to fifth to seventh generations of airways with the dimensions based on the Weibel's pulmonary model. Computations were carried out for different Reynolds numbers (Re = 400 and 1000) and Womersley numbers (α = 4 and 16) to study the air flow fields, gas transportation, and wall shear stresses in the lung airways. Flow structure was compared with experimental results. Both 2D and 3D numerical models successfully reproduced many results observed in the experiment. The oxygen concentration distribution in the lung model was investigated to analyze the influence of flow oscillation on gas transport inside the lung model.
3D EM imaging from a single borehole; a numerical feasibility study
Alumbaugh, D.L.; Wilt, M.J.
1998-07-01
Often the drilling of an oil well is followed by a logging process to characterize the region immediately surrounding the well bore. The electromagnetic (EM) induction tool, which provides the formation resistivity, is among the most frequently run logs. A preliminary study has been conducted to analyze the feasibility of three dimensional (3D) electromagnetic (EM) imaging from a single borehole. The logging tool consists of a vertical magnetic dipole source and multiple 3 component magnetic field receivers offset at different distances from the source. Synthetic data calculated with a 3D finite difference code demonstrate that the phase of the horizontal magnetic fields provides the critical information on the three dimensionality of the medium. A 3D inversion algorithm is then employed to demonstrate the plausibility of 3D inversion using 3 component magnetic field data. Finally, problems associated with introducing biased noise into the horizontal components of the field through misalignment of the logging tool is discussed.
Turbomachinery Heat Transfer and Loss Modeling for 3D Navier-Stokes Codes
NASA Technical Reports Server (NTRS)
DeWitt, Kenneth; Ameri, Ali
2005-01-01
This report's contents focus on making use of NASA Glenn on-site computational facilities,to develop, validate, and apply models for use in advanced 3D Navier-Stokes Computational Fluid Dynamics (CFD) codes to enhance the capability to compute heat transfer and losses in turbomachiney.
3D Neutron Transport PWR Full-core Calculation with RMC code
NASA Astrophysics Data System (ADS)
Qiu, Yishu; She, Ding; Fan, Xiao; Wang, Kan; Li, Zeguang; Liang, Jingang; Leroyer, Hadrien
2014-06-01
Nowadays, there are more and more interests in the use of Monte Carlo codes to calculate the detailed power density distributions in full-core reactors. With the Inspur TS1000 HPC Server of Tsinghua University, several calculations have been done based on the EDF 3D Neutron Transport PWR Full-core benchmark through large-scale parallelism. To investigate and compare the results of the deterministic method and Monte Carlo method, EDF R&D and Department of Engineering Physics of Tsinghua University are having a collaboration to make code to code verification. So in this paper, two codes are used. One is the code COCAGNE developed by the EDF R&D, a deterministic core code, and the other is the Monte Carlo code RMC developed by Department of Engineering Physics in Tsinghua University. First, the full-core model is described and a 26-group calculation was performed by these two codes using the same 26-group cross-section library provided by EDF R&D. Then the parallel and tally performance of RMC is discussed. RMC employs a novel algorithm which can cut down most of the communications. It can be seen clearly that the speedup ratio almost linearly increases with the nodes. Furthermore the cell-mapping method applied by RMC consumes little time to tally even millions of cells. The results of the codes COCAGNE and RMC are compared in three ways. The results of these two codes agree well with each other. It can be concluded that both COCAGNE and RMC are able to provide 3D-transport solutions associated with detailed power density distributions calculation in PWR full-core reactors. Finally, to investigate how many histories are needed to obtain a given standard deviation for a full 3D solution, the non-symmetrized condensed 2-group fluxes of RMC are discussed.
LINFLUX-AE: A Turbomachinery Aeroelastic Code Based on a 3-D Linearized Euler Solver
NASA Technical Reports Server (NTRS)
Reddy, T. S. R.; Bakhle, M. A.; Trudell, J. J.; Mehmed, O.; Stefko, G. L.
2004-01-01
This report describes the development and validation of LINFLUX-AE, a turbomachinery aeroelastic code based on the linearized unsteady 3-D Euler solver, LINFLUX. A helical fan with flat plate geometry is selected as the test case for numerical validation. The steady solution required by LINFLUX is obtained from the nonlinear Euler/Navier Stokes solver TURBO-AE. The report briefly describes the salient features of LINFLUX and the details of the aeroelastic extension. The aeroelastic formulation is based on a modal approach. An eigenvalue formulation is used for flutter analysis. The unsteady aerodynamic forces required for flutter are obtained by running LINFLUX for each mode, interblade phase angle and frequency of interest. The unsteady aerodynamic forces for forced response analysis are obtained from LINFLUX for the prescribed excitation, interblade phase angle, and frequency. The forced response amplitude is calculated from the modal summation of the generalized displacements. The unsteady pressures, work done per cycle, eigenvalues and forced response amplitudes obtained from LINFLUX are compared with those obtained from LINSUB, TURBO-AE, ASTROP2, and ANSYS.
NASA Astrophysics Data System (ADS)
McFall, B. C.; Fritz, H. M.; Horrillo, J. J.; Mohammed, F.
2014-12-01
Landslide generated tsunamis such as Lituya Bay, Alaska 1958 account for some of highest recorded tsunami runup heights. Source and runup scenarios based on real world events are physically modeled using generalized Froude similarity in the three dimensional NEES tsunami wave basin at Oregon State University. A novel pneumatic landslide tsunami generator (LTG) was deployed to simulate landslides with varying geometry and kinematics. The bathymetric and topographic scenarios tested with the LTG are the basin-wide propagation and runup, fjord, curved headland fjord and a conical island setting representing a landslide off an island or a volcano flank collapse. The LTG consists of a sliding box filled with 1,350 kg of landslide material which is accelerated by pneumatic pistons down slope. Two different landslide materials are used to study the granulometry effects: naturally rounded river gravel and cobble mixtures. Water surface elevations are recorded by an array of resistance wave gauges. The landslide deformation is measured from above and underwater camera recordings. The landslide deposit is measured on the basin floor with a multiple transducer acoustic array (MTA). Landslide surface reconstruction and kinematics are determined with a stereo particle image velocimetry (PIV) system. Wave runup is recorded with resistance wave gauges along the slope and verified with video image processing. The measured landslide and wave parameters are compared between the planar hill slope used in various scenarios and the convex hill slope of the conical island. The energy conversion rates from the landslide motion to the wave train is quantified for the planar and convex hill slopes. The wave runup data on the opposing headland is analyzed and evaluated with wave theories. The measured landslide and tsunami data serve to validate and advance three-dimensional numerical landslide tsunami prediction models. Two 3D Navier-Stokes models were tested, the commercial code FLOW-3D
Users manual for the NASA Lewis three-dimensional ice accretion code (LEWICE 3D)
NASA Technical Reports Server (NTRS)
Bidwell, Colin S.; Potapczuk, Mark G.
1993-01-01
A description of the methodology, the algorithms, and the input and output data along with an example case for the NASA Lewis 3D ice accretion code (LEWICE3D) has been produced. The manual has been designed to help the user understand the capabilities, the methodologies, and the use of the code. The LEWICE3D code is a conglomeration of several codes for the purpose of calculating ice shapes on three-dimensional external surfaces. A three-dimensional external flow panel code is incorporated which has the capability of calculating flow about arbitrary 3D lifting and nonlifting bodies with external flow. A fourth order Runge-Kutta integration scheme is used to calculate arbitrary streamlines. An Adams type predictor-corrector trajectory integration scheme has been included to calculate arbitrary trajectories. Schemes for calculating tangent trajectories, collection efficiencies, and concentration factors for arbitrary regions of interest for single droplets or droplet distributions have been incorporated. A LEWICE 2D based heat transfer algorithm can be used to calculate ice accretions along surface streamlines. A geometry modification scheme is incorporated which calculates the new geometry based on the ice accretions generated at each section of interest. The three-dimensional ice accretion calculation is based on the LEWICE 2D calculation. Both codes calculate the flow, pressure distribution, and collection efficiency distribution along surface streamlines. For both codes the heat transfer calculation is divided into two regions, one above the stagnation point and one below the stagnation point, and solved for each region assuming a flat plate with pressure distribution. Water is assumed to follow the surface streamlines, hence starting at the stagnation zone any water that is not frozen out at a control volume is assumed to run back into the next control volume. After the amount of frozen water at each control volume has been calculated the geometry is modified by
Computational time analysis of the numerical solution of 3D electrostatic Poisson's equation
NASA Astrophysics Data System (ADS)
Kamboh, Shakeel Ahmed; Labadin, Jane; Rigit, Andrew Ragai Henri; Ling, Tech Chaw; Amur, Khuda Bux; Chaudhary, Muhammad Tayyab
2015-05-01
3D Poisson's equation is solved numerically to simulate the electric potential in a prototype design of electrohydrodynamic (EHD) ion-drag micropump. Finite difference method (FDM) is employed to discretize the governing equation. The system of linear equations resulting from FDM is solved iteratively by using the sequential Jacobi (SJ) and sequential Gauss-Seidel (SGS) methods, simulation results are also compared to examine the difference between the results. The main objective was to analyze the computational time required by both the methods with respect to different grid sizes and parallelize the Jacobi method to reduce the computational time. In common, the SGS method is faster than the SJ method but the data parallelism of Jacobi method may produce good speedup over SGS method. In this study, the feasibility of using parallel Jacobi (PJ) method is attempted in relation to SGS method. MATLAB Parallel/Distributed computing environment is used and a parallel code for SJ method is implemented. It was found that for small grid size the SGS method remains dominant over SJ method and PJ method while for large grid size both the sequential methods may take nearly too much processing time to converge. Yet, the PJ method reduces computational time to some extent for large grid sizes.
Numerical Calculations of 3-D High-Lift Flows and Comparison with Experiment
NASA Technical Reports Server (NTRS)
Compton, William B, III
2015-01-01
Solutions were obtained with the Navier-Stokes CFD code TLNS3D to predict the flow about the NASA Trapezoidal Wing, a high-lift wing composed of three elements: the main-wing element, a deployed leading-edge slat, and a deployed trailing-edge flap. Turbulence was modeled by the Spalart-Allmaras one-equation turbulence model. One case with massive separation was repeated using Menter's two-equation SST (Menter's Shear Stress Transport) k-omega turbulence model in an attempt to improve the agreement with experiment. The investigation was conducted at a free stream Mach number of 0.2, and at angles of attack ranging from 10.004 degrees to 34.858 degrees. The Reynolds number based on the mean aerodynamic chord of the wing was 4.3 x 10 (sup 6). Compared to experiment, the numerical procedure predicted the surface pressures very well at angles of attack in the linear range of the lift. However, computed maximum lift was 5% low. Drag was mainly under predicted. The procedure correctly predicted several well-known trends and features of high-lift flows, such as off-body separation. The two turbulence models yielded significantly different solutions for the repeated case.
A new 3-D integral code for computation of accelerator magnets
Turner, L.R.; Kettunen, L.
1991-01-01
For computing accelerator magnets, integral codes have several advantages over finite element codes; far-field boundaries are treated automatically, and computed field in the bore region satisfy Maxwell's equations exactly. A new integral code employing edge elements rather than nodal elements has overcome the difficulties associated with earlier integral codes. By the use of field integrals (potential differences) as solution variables, the number of unknowns is reduced to one less than the number of nodes. Two examples, a hollow iron sphere and the dipole magnet of Advanced Photon Source injector synchrotron, show the capability of the code. The CPU time requirements are comparable to those of three-dimensional (3-D) finite-element codes. Experiments show that in practice it can realize much of the potential CPU time saving that parallel processing makes possible. 8 refs., 4 figs., 1 tab.
NASA Astrophysics Data System (ADS)
Mandumpala Devassy, B.; Edelbauer, W.; Greif, D.
2015-12-01
Cavitation and its effect on spray formation and its dispersion play a crucial role in proper engine combustion and controlled emission. This study focuses on these effects in a typical common rail 6-hole diesel injector accounting for 3D needle movement and flow compressibility effects. Coupled numerical simulations using 1D and 3D CFD codes are used for this investigation. Previous studies in this direction have already presented a detailed structure of the adopted methodology. Compared to the previous analysis, the present study investigates the effect of 3D needle movement and cavitation on the spray formation for pilot and main injection events for a typical diesel engine operating point. The present setup performs a 3D compressible multiphase simulation coupled with a standalone 1D high pressure flow simulation. The simulation proceeds by the mutual communication between 1D and 3D solvers. In this work a typical common rail injector with a mini-sac nozzle is studied. The lateral and radial movement of the needle and its effect on the cavitation generation and the subsequent spray penetration are analyzed. The result indicates the effect of compressibility of the liquid on damping the needle forces, and also the difference in the spray penetration levels due to the asymmetrical flow field. Therefore, this work intends to provide an efficient and user-friendly engineering tool for simulating a complete fuel injector including spray propagation.
RELAP5-3D Code Includes Athena Features and Models
Richard A. Riemke; Cliff B. Davis; Richard R. Schultz
2006-07-01
Version 2.3 of the RELAP5-3D computer program includes all features and models previously available only in the ATHENA version of the code. These include the addition of new working fluids (i.e., ammonia, blood, carbon dioxide, glycerol, helium, hydrogen, lead-bismuth, lithium, lithium-lead, nitrogen, potassium, sodium, and sodium-potassium) and a magnetohydrodynamic model that expands the capability of the code to model many more thermal-hydraulic systems. In addition to the new working fluids along with the standard working fluid water, one or more noncondensable gases (e.g., air, argon, carbon dioxide, carbon monoxide, helium, hydrogen, krypton, nitrogen, oxygen, sf6, xenon) can be specified as part of the vapor/gas phase of the working fluid. These noncondensable gases were in previous versions of RELAP5- 3D. Recently four molten salts have been added as working fluids to RELAP5-3D Version 2.4, which has had limited release. These molten salts will be in RELAP5-3D Version 2.5, which will have a general release like RELAP5-3D Version 2.3. Applications that use these new features and models are discussed in this paper.
RELAP5-3D Code Includes ATHENA Features and Models
Riemke, Richard A.; Davis, Cliff B.; Schultz, Richard R.
2006-07-01
Version 2.3 of the RELAP5-3D computer program includes all features and models previously available only in the ATHENA version of the code. These include the addition of new working fluids (i.e., ammonia, blood, carbon dioxide, glycerol, helium, hydrogen, lead-bismuth, lithium, lithium-lead, nitrogen, potassium, sodium, and sodium-potassium) and a magnetohydrodynamic model that expands the capability of the code to model many more thermal-hydraulic systems. In addition to the new working fluids along with the standard working fluid water, one or more noncondensable gases (e.g., air, argon, carbon dioxide, carbon monoxide, helium, hydrogen, krypton, nitrogen, oxygen, SF{sub 6}, xenon) can be specified as part of the vapor/gas phase of the working fluid. These noncondensable gases were in previous versions of RELAP5-3D. Recently four molten salts have been added as working fluids to RELAP5-3D Version 2.4, which has had limited release. These molten salts will be in RELAP5-3D Version 2.5, which will have a general release like RELAP5-3D Version 2.3. Applications that use these new features and models are discussed in this paper. (authors)
Edge Transport Modeling using the 3D EMC3-Eirene code on Tokamaks and Stellarators
NASA Astrophysics Data System (ADS)
Lore, J. D.; Ahn, J. W.; Briesemeister, A.; Ferraro, N.; Labombard, B.; McLean, A.; Reinke, M.; Shafer, M.; Terry, J.
2015-11-01
The fluid plasma edge transport code EMC3-Eirene has been applied to aid data interpretation and understanding the results of experiments with 3D effects on several tokamaks. These include applied and intrinsic 3D magnetic fields, 3D plasma facing components, and toroidally and poloidally localized heat and particle sources. On Alcator C-Mod, a series of experiments explored the impact of toroidally and poloidally localized impurity gas injection on core confinement and asymmetries in the divertor fluxes, with the differences between the asymmetry in L-mode and H-mode qualitatively reproduced in the simulations due to changes in the impurity ionization in the private flux region. Modeling of NSTX experiments on the effect of 3D fields on detachment matched the trend of a higher density at which the detachment occurs when 3D fields are applied. On DIII-D, different magnetic field models were used in the simulation and compared against the 2D Thomson scattering diagnostic. In simulating each device different aspects of the code model are tested pointing to areas where the model must be further developed. The application to stellarator experiments will also be discussed. Work supported by U.S. DOE: DE-AC05-00OR22725, DE AC02-09CH11466, DE-FC02-99ER54512, and DE-FC02-04ER54698.
Inheritance of pre-existing weakness in continental breakup: 3D numerical modeling
NASA Astrophysics Data System (ADS)
Liao, Jie; Gerya, Taras
2013-04-01
breakup order of crust and mantle (Huismans and Beaumont, 2011). However, the inheritance of pre-existing lithospheric weakness in the evolution of continental rifts and oceanic ridge is not well studied. We use 3D numerical modeling to study this problem, by changing the weak zone position and geometry, and the rheological structure of the model. In our study, we find that: 1).3D continental breakup and seafloor spreading patterns are controlled by (a) crust-mantle rheological coupling and (b) geometry and position of the pre-existing weak zones. 2).Three spreading patterns are obtained: (a) straight ridges, (b) curved ridges and (c) overlapping ridges. 3).When crust and mantle are decoupled, abandoned rift structures often form.
The Monte Carlo SRNA-VOX code for 3D proton dose distribution in voxelized geometry using CT data
NASA Astrophysics Data System (ADS)
Ilic, Radovan D.; Spasic-Jokic, Vesna; Belicev, Petar; Dragovic, Milos
2005-03-01
This paper describes the application of the SRNA Monte Carlo package for proton transport simulations in complex geometry and different material compositions. The SRNA package was developed for 3D dose distribution calculation in proton therapy and dosimetry and it was based on the theory of multiple scattering. The decay of proton induced compound nuclei was simulated by the Russian MSDM model and our own using ICRU 63 data. The developed package consists of two codes: the SRNA-2KG, which simulates proton transport in combinatorial geometry and the SRNA-VOX, which uses the voxelized geometry using the CT data and conversion of the Hounsfield's data to tissue elemental composition. Transition probabilities for both codes are prepared by the SRNADAT code. The simulation of the proton beam characterization by multi-layer Faraday cup, spatial distribution of positron emitters obtained by the SRNA-2KG code and intercomparison of computational codes in radiation dosimetry, indicate immediate application of the Monte Carlo techniques in clinical practice. In this paper, we briefly present the physical model implemented in the SRNA package, the ISTAR proton dose planning software, as well as the results of the numerical experiments with proton beams to obtain 3D dose distribution in the eye and breast tumour.
The Monte Carlo SRNA-VOX code for 3D proton dose distribution in voxelized geometry using CT data.
Ilić, Radovan D; Spasić-Jokić, Vesna; Belicev, Petar; Dragović, Milos
2005-03-01
This paper describes the application of the SRNA Monte Carlo package for proton transport simulations in complex geometry and different material compositions. The SRNA package was developed for 3D dose distribution calculation in proton therapy and dosimetry and it was based on the theory of multiple scattering. The decay of proton induced compound nuclei was simulated by the Russian MSDM model and our own using ICRU 63 data. The developed package consists of two codes: the SRNA-2KG, which simulates proton transport in combinatorial geometry and the SRNA-VOX, which uses the voxelized geometry using the CT data and conversion of the Hounsfield's data to tissue elemental composition. Transition probabilities for both codes are prepared by the SRNADAT code. The simulation of the proton beam characterization by multi-layer Faraday cup, spatial distribution of positron emitters obtained by the SRNA-2KG code and intercomparison of computational codes in radiation dosimetry, indicate immediate application of the Monte Carlo techniques in clinical practice. In this paper, we briefly present the physical model implemented in the SRNA package, the ISTAR proton dose planning software, as well as the results of the numerical experiments with proton beams to obtain 3D dose distribution in the eye and breast tumour. PMID:15798273
Development of Unsteady Aerodynamic and Aeroelastic Reduced-Order Models Using the FUN3D Code
NASA Technical Reports Server (NTRS)
Silva, Walter A.; Vatsa, Veer N.; Biedron, Robert T.
2009-01-01
Recent significant improvements to the development of CFD-based unsteady aerodynamic reduced-order models (ROMs) are implemented into the FUN3D unstructured flow solver. These improvements include the simultaneous excitation of the structural modes of the CFD-based unsteady aerodynamic system via a single CFD solution, minimization of the error between the full CFD and the ROM unsteady aero- dynamic solution, and computation of a root locus plot of the aeroelastic ROM. Results are presented for a viscous version of the two-dimensional Benchmark Active Controls Technology (BACT) model and an inviscid version of the AGARD 445.6 aeroelastic wing using the FUN3D code.
Ui, Atsushi; Miyaji, Takamasa
2004-10-15
The best-estimate coupled three-dimensional (3-D) core and thermal-hydraulic code system TRAC-BF1/COS3D has been developed. COS3D, based on a modified one-group neutronic model, is a 3-D core simulator used for licensing analyses and core management of commercial boiling water reactor (BWR) plants in Japan. TRAC-BF1 is a plant simulator based on a two-fluid model. TRAC-BF1/COS3D is a coupled system of both codes, which are connected using a parallel computing tool. This code system was applied to the OECD/NRC BWR Turbine Trip Benchmark. Since the two-group cross-section tables are provided by the benchmark team, COS3D was modified to apply to this specification. Three best-estimate scenarios and four hypothetical scenarios were calculated using this code system. In the best-estimate scenario, the predicted core power with TRAC-BF1/COS3D is slightly underestimated compared with the measured data. The reason seems to be a slight difference in the core boundary conditions, that is, pressure changes and the core inlet flow distribution, because the peak in this analysis is sensitive to them. However, the results of this benchmark analysis show that TRAC-BF1/COS3D gives good precision for the prediction of the actual BWR transient behavior on the whole. Furthermore, the results with the modified one-group model and the two-group model were compared to verify the application of the modified one-group model to this benchmark. This comparison shows that the results of the modified one-group model are appropriate and sufficiently precise.
A 3D-CFD code for accurate prediction of fluid flows and fluid forces in seals
NASA Astrophysics Data System (ADS)
Athavale, M. M.; Przekwas, A. J.; Hendricks, R. C.
1994-01-01
Current and future turbomachinery requires advanced seal configurations to control leakage, inhibit mixing of incompatible fluids and to control the rotodynamic response. In recognition of a deficiency in the existing predictive methodology for seals, a seven year effort was established in 1990 by NASA's Office of Aeronautics Exploration and Technology, under the Earth-to-Orbit Propulsion program, to develop validated Computational Fluid Dynamics (CFD) concepts, codes and analyses for seals. The effort will provide NASA and the U.S. Aerospace Industry with advanced CFD scientific codes and industrial codes for analyzing and designing turbomachinery seals. An advanced 3D CFD cylindrical seal code has been developed, incorporating state-of-the-art computational methodology for flow analysis in straight, tapered and stepped seals. Relevant computational features of the code include: stationary/rotating coordinates, cylindrical and general Body Fitted Coordinates (BFC) systems, high order differencing schemes, colocated variable arrangement, advanced turbulence models, incompressible/compressible flows, and moving grids. This paper presents the current status of code development, code demonstration for predicting rotordynamic coefficients, numerical parametric study of entrance loss coefficients for generic annular seals, and plans for code extensions to labyrinth, damping, and other seal configurations.
Li, Yong Gang; Yang, Yang; Short, Michael P.; Ding, Ze Jun; Zeng, Zhi; Li, Ju
2015-01-01
SRIM-like codes have limitations in describing general 3D geometries, for modeling radiation displacements and damage in nanostructured materials. A universal, computationally efficient and massively parallel 3D Monte Carlo code, IM3D, has been developed with excellent parallel scaling performance. IM3D is based on fast indexing of scattering integrals and the SRIM stopping power database, and allows the user a choice of Constructive Solid Geometry (CSG) or Finite Element Triangle Mesh (FETM) method for constructing 3D shapes and microstructures. For 2D films and multilayers, IM3D perfectly reproduces SRIM results, and can be ∼102 times faster in serial execution and > 104 times faster using parallel computation. For 3D problems, it provides a fast approach for analyzing the spatial distributions of primary displacements and defect generation under ion irradiation. Herein we also provide a detailed discussion of our open-source collision cascade physics engine, revealing the true meaning and limitations of the “Quick Kinchin-Pease” and “Full Cascades” options. The issues of femtosecond to picosecond timescales in defining displacement versus damage, the limitation of the displacements per atom (DPA) unit in quantifying radiation damage (such as inadequacy in quantifying degree of chemical mixing), are discussed. PMID:26658477
NASA Astrophysics Data System (ADS)
Li, Yong Gang; Yang, Yang; Short, Michael P.; Ding, Ze Jun; Zeng, Zhi; Li, Ju
2015-12-01
SRIM-like codes have limitations in describing general 3D geometries, for modeling radiation displacements and damage in nanostructured materials. A universal, computationally efficient and massively parallel 3D Monte Carlo code, IM3D, has been developed with excellent parallel scaling performance. IM3D is based on fast indexing of scattering integrals and the SRIM stopping power database, and allows the user a choice of Constructive Solid Geometry (CSG) or Finite Element Triangle Mesh (FETM) method for constructing 3D shapes and microstructures. For 2D films and multilayers, IM3D perfectly reproduces SRIM results, and can be ∼102 times faster in serial execution and > 104 times faster using parallel computation. For 3D problems, it provides a fast approach for analyzing the spatial distributions of primary displacements and defect generation under ion irradiation. Herein we also provide a detailed discussion of our open-source collision cascade physics engine, revealing the true meaning and limitations of the “Quick Kinchin-Pease” and “Full Cascades” options. The issues of femtosecond to picosecond timescales in defining displacement versus damage, the limitation of the displacements per atom (DPA) unit in quantifying radiation damage (such as inadequacy in quantifying degree of chemical mixing), are discussed.
Equation-of-State Test Suite for the DYNA3D Code
Benjamin, Russell D.
2015-11-05
This document describes the creation and implementation of a test suite for the Equationof- State models in the DYNA3D code. A customized input deck has been created for each model, as well as a script that extracts the relevant data from the high-speed edit file created by DYNA3D. Each equation-of-state model is broken apart and individual elements of the model are tested, as well as testing the entire model. The input deck for each model is described and the results of the tests are discussed. The intent of this work is to add this test suite to the validation suite presently used for DYNA3D.
Vay, J.-L.; Furman, M.A.; Azevedo, A.W.; Cohen, R.H.; Friedman, A.; Grote, D.P.; Stoltz, P.H.
2004-04-19
We have integrated the electron-cloud code POSINST [1] with WARP [2]--a 3-D parallel Particle-In-Cell accelerator code developed for Heavy Ion Inertial Fusion--so that the two can interoperate. Both codes are run in the same process, communicate through a Python interpreter (already used in WARP), and share certain key arrays (so far, particle positions and velocities). Currently, POSINST provides primary and secondary sources of electrons, beam bunch kicks, a particle mover, and diagnostics. WARP provides the field solvers and diagnostics. Secondary emission routines are provided by the Tech-X package CMEE.
Object-adaptive depth compensated inter prediction for depth video coding in 3D video system
NASA Astrophysics Data System (ADS)
Kang, Min-Koo; Lee, Jaejoon; Lim, Ilsoon; Ho, Yo-Sung
2011-01-01
Nowadays, the 3D video system using the MVD (multi-view video plus depth) data format is being actively studied. The system has many advantages with respect to virtual view synthesis such as an auto-stereoscopic functionality, but compression of huge input data remains a problem. Therefore, efficient 3D data compression is extremely important in the system, and problems of low temporal consistency and viewpoint correlation should be resolved for efficient depth video coding. In this paper, we propose an object-adaptive depth compensated inter prediction method to resolve the problems where object-adaptive mean-depth difference between a current block, to be coded, and a reference block are compensated during inter prediction. In addition, unique properties of depth video are exploited to reduce side information required for signaling decoder to conduct the same process. To evaluate the coding performance, we have implemented the proposed method into MVC (multiview video coding) reference software, JMVC 8.2. Experimental results have demonstrated that our proposed method is especially efficient for depth videos estimated by DERS (depth estimation reference software) discussed in the MPEG 3DV coding group. The coding gain was up to 11.69% bit-saving, and it was even increased when we evaluated it on synthesized views of virtual viewpoints.
Numerical classification of coding sequences
NASA Technical Reports Server (NTRS)
Collins, D. W.; Liu, C. C.; Jukes, T. H.
1992-01-01
DNA sequences coding for protein may be represented by counts of nucleotides or codons. A complete reading frame may be abbreviated by its base count, e.g. A76C158G121T74, or with the corresponding codon table, e.g. (AAA)0(AAC)1(AAG)9 ... (TTT)0. We propose that these numerical designations be used to augment current methods of sequence annotation. Because base counts and codon tables do not require revision as knowledge of function evolves, they are well-suited to act as cross-references, for example to identify redundant GenBank entries. These descriptors may be compared, in place of DNA sequences, to extract homologous genes from large databases. This approach permits rapid searching with good selectivity.
NASA Astrophysics Data System (ADS)
Velioǧlu, Deniz; Cevdet Yalçıner, Ahmet; Zaytsev, Andrey
2016-04-01
Tsunamis are huge waves with long wave periods and wave lengths that can cause great devastation and loss of life when they strike a coast. The interest in experimental and numerical modeling of tsunami propagation and inundation increased considerably after the 2011 Great East Japan earthquake. In this study, two numerical codes, FLOW 3D and NAMI DANCE, that analyze tsunami propagation and inundation patterns are considered. Flow 3D simulates linear and nonlinear propagating surface waves as well as long waves by solving three-dimensional Navier-Stokes (3D-NS) equations. NAMI DANCE uses finite difference computational method to solve 2D depth-averaged linear and nonlinear forms of shallow water equations (NSWE) in long wave problems, specifically tsunamis. In order to validate these two codes and analyze the differences between 3D-NS and 2D depth-averaged NSWE equations, two benchmark problems are applied. One benchmark problem investigates the runup of long waves over a complex 3D beach. The experimental setup is a 1:400 scale model of Monai Valley located on the west coast of Okushiri Island, Japan. Other benchmark problem is discussed in 2015 National Tsunami Hazard Mitigation Program (NTHMP) Annual meeting in Portland, USA. It is a field dataset, recording the Japan 2011 tsunami in Hilo Harbor, Hawaii. The computed water surface elevation and velocity data are compared with the measured data. The comparisons showed that both codes are in fairly good agreement with each other and benchmark data. The differences between 3D-NS and 2D depth-averaged NSWE equations are highlighted. All results are presented with discussions and comparisons. Acknowledgements: Partial support by Japan-Turkey Joint Research Project by JICA on earthquakes and tsunamis in Marmara Region (JICA SATREPS - MarDiM Project), 603839 ASTARTE Project of EU, UDAP-C-12-14 project of AFAD Turkey, 108Y227, 113M556 and 213M534 projects of TUBITAK Turkey, RAPSODI (CONCERT_Dis-021) of CONCERT
Numerical 3D models support two distinct hydrothermal circulation systems at fast spreading ridges
NASA Astrophysics Data System (ADS)
Hasenclever, Jörg; Theissen-Krah, Sonja; Rüpke, Lars
2013-04-01
We present 3D numerical calculations of hydrothermal fluid flow at fast spreading ridges. The setup of the 3D models is based our previous 2D studies, in which we have coupled numerical models for crustal accretion and hydrothermal fluid flow. One result of these calculations is a crustal permeability field that leads to a thermal structure in the crust that matches seismic tomography data of the East Pacific Rise (EPR). The 1000°C isotherm obtained from the 2D results is now used as the lower boundary of the 3D model domain, while the upper boundary is a smoothed bathymetry of the EPR. The same permeability field as in the 2D models is used, with the highest permeability at the ridge axis and a decrease with both depth and distance to the ridge. Permeability is also reduced linearly between 600 and 1000°C. Using a newly developed parallel finite element code written in Matlab that solves for thermal evolution, fluid pressure and Darcy flow, we simulate the flow patterns of hydrothermal circulation in a segment of 5000m along-axis, 10000m across-axis and up to 5000m depth. We observe two distinct hydrothermal circulation systems: An on-axis system forming a series of vents with a spacing ranging from 100 to 500m that is recharged by nearby (100-200m) downflows on both sides of the ridge axis. Simultaneously a second system with much broader extensions both laterally and vertically exists off-axis. It is recharged by fluids intruding between 1500m to 5000m off-axis and sampling both upper and lower crust. These fluids are channeled in the deepest and hottest regions with high permeability and migrate up-slope following the 600°C isotherm until reaching the edge of the melt lens. Depending on the width of the melt lens these off-axis fluids either merge with the on-axis hydrothermal system or form separate vents. We observe separate off-axis vent fields if the magma lens half-width exceeds 1000m and confluence of both systems for half-widths smaller than 500m. For
Wall touching kink mode calculations with the M3D code
NASA Astrophysics Data System (ADS)
Breslau, J. A.
2014-10-01
In recent years there have been a number of results published concerning the transient vessel currents and forces occurring during a tokamak VDE, as predicted by simulations with the nonlinear MHD code M3D. The nature of the simulations is such that these currents and forces occur at the boundary of the computational domain, making the proper choice of boundary conditions critical to the reliability of the results. The M3D boundary condition includes the prescription that the normal component of the velocity vanish at the wall. It has been argued that this prescription invalidates the calculations because it would seem to rule out the possibility of advection of plasma surface currents into the wall. This claim has been tested by applying M3D to an idealized case - a kink-unstable plasma column - in order to abstract the essential physics from the complications involved in the attempt to model real devices. While comparison of the results is complicated by effects arising from the higher dimensionality and complexity of M3D, we have verified that M3D is capable of reproducing both the correct saturation behavior of the free boundary kink and the ``Hiro'' currents arising when the kink interacts with a conducting tile surface interior to the ideal wall.
Comparison of a 3-D GPU-Assisted Maxwell Code and Ray Tracing for Reflectometry on ITER
NASA Astrophysics Data System (ADS)
Gady, Sarah; Kubota, Shigeyuki; Johnson, Irena
2015-11-01
Electromagnetic wave propagation and scattering in magnetized plasmas are important diagnostics for high temperature plasmas. 1-D and 2-D full-wave codes are standard tools for measurements of the electron density profile and fluctuations; however, ray tracing results have shown that beam propagation in tokamak plasmas is inherently a 3-D problem. The GPU-Assisted Maxwell Code utilizes the FDTD (Finite-Difference Time-Domain) method for solving the Maxwell equations with the cold plasma approximation in a 3-D geometry. Parallel processing with GPGPU (General-Purpose computing on Graphics Processing Units) is used to accelerate the computation. Previously, we reported on initial comparisons of the code results to 1-D numerical and analytical solutions, where the size of the computational grid was limited by the on-board memory of the GPU. In the current study, this limitation is overcome by using domain decomposition and an additional GPU. As a practical application, this code is used to study the current design of the ITER Low Field Side Reflectometer (LSFR) for the Equatorial Port Plug 11 (EPP11). A detailed examination of Gaussian beam propagation in the ITER edge plasma will be presented, as well as comparisons with ray tracing. This work was made possible by funding from the Department of Energy for the Summer Undergraduate Laboratory Internship (SULI) program. This work is supported by the US DOE Contract No.DE-AC02-09CH11466 and DE-FG02-99-ER54527.
Comprehensive study of numerical anisotropy and dispersion in 3-D TLM meshes
NASA Astrophysics Data System (ADS)
Berini, Pierre; Wu, Ke
1995-05-01
This paper presents a comprehensive analysis of the numerical anisotropy and dispersion of 3-D TLM meshes constructed using several generalized symmetrical condensed TLM nodes. The dispersion analysis is performed in isotropic lossless, isotropic lossy and anisotropic lossless media and yields a comparison of the simulation accuracy for the different TLM nodes. The effect of mesh grading on the numerical dispersion is also determined. The results compare meshes constructed with Johns' symmetrical condensed node (SCN), two hybrid symmetrical condensed nodes (HSCN) and two frequency domain symmetrical condensed nodes (FDSCN). It has been found that under certain circumstances, the time domain nodes may introduce numerical anisotropy when modelling isotropic media.
Preliminary results of 3D dose calculations with MCNP-4B code from a SPECT image.
Rodríguez Gual, M; Lima, F F; Sospedra Alfonso, R; González González, J; Calderón Marín, C
2004-01-01
Interface software was developed to generate the input file to run Monte Carlo MCNP-4B code from medical image in Interfile format version 3.3. The software was tested using a spherical phantom of tomography slides with known cumulated activity distribution in Interfile format generated with IMAGAMMA medical image processing system. The 3D dose calculation obtained with Monte Carlo MCNP-4B code was compared with the voxel S factor method. The results show a relative error between both methods less than 1 %. PMID:15625058
Simulations of implosions with a 3D, parallel, unstructured-grid, radiation-hydrodynamics code
Kaiser, T B; Milovich, J L; Prasad, M K; Rathkopf, J; Shestakov, A I
1998-12-28
An unstructured-grid, radiation-hydrodynamics code is used to simulate implosions. Although most of the problems are spherically symmetric, they are run on 3D, unstructured grids in order to test the code's ability to maintain spherical symmetry of the converging waves. Three problems, of increasing complexity, are presented. In the first, a cold, spherical, ideal gas bubble is imploded by an enclosing high pressure source. For the second, we add non-linear heat conduction and drive the implosion with twelve laser beams centered on the vertices of an icosahedron. In the third problem, a NIF capsule is driven with a Planckian radiation source.
Dynamic coupling between fluid flow and vein growth in fractures: a 3D numerical model
NASA Astrophysics Data System (ADS)
Schwarz, J.-O.; Enzmann, F.
2012-04-01
Fluid flow is one of the main mass transport mechanisms in the Earth's crust and abundant mineral vein networks are important indicators for fluid flow and fluid rock interaction. These systems are dynamic and part of the so called RTM processes (reaction-transport-mechanics). Understanding of mineral vein systems requires coupling of these processes. Here we present a conceptional model for dynamic vein growth of syntaxial, posttectonic veins generated by advective fluid flow and show first results of a numerical model for this scenario. Vein generation requires three processes to occur: (i) fracture generation by mechanical stress e.g. hydro-fracturing, (ii) flow of a supersaturated fluid on that fracture and (iii) crystallization of phase(s) on or in the fracture. 3D synthetic fractures are generated with the SynFrac code (Ogilvie, et al. 2006). Subsequently solutions of the Navier-Stokes equation for this fracture are computed by a computational fluid dynamics code called GeoDict (Wiegmann 2007). Transport (advective and diffusive) of chemical species to growth sites in the fracture and vein growth are computed by a self-written MATLAB script. The numerical model discretizes the wall rock and fracture geometry by volumetric pixels (voxels). Based on this representation, the model computes the three basic functions for vein generation: (a) nucleation, (b) fluid flow with transport of chemical species and (c) growth. The following conditions were chosen for these three modules. Nucleation is heterogeneous and occurs instantaneously at the wall rock/fracture interface. Advective and diffusive flow of a supersaturated fluid and related transport of chemical species occurs according to the computed fluid flow field by GeoDict. Concentration of chemical species at the inflow is constant, representing external fluid buffering. Changes/decrease in the concentration of chemical species occurs only due to vein growth. Growth of nuclei is limited either by transport of
The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications
NASA Technical Reports Server (NTRS)
Bravo, Ramiro H.; Chen, Ching-Jen
1992-01-01
In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.
3D numerical simulation analysis of passive drag near free surface in swimming
NASA Astrophysics Data System (ADS)
Zhan, Jie-min; Li, Tian-zeng; Chen, Xue-bin; Li, Yok-sheung; Wai, Wing-hong Onyx
2015-04-01
The aim of this work is to build a 3D numerical model to study the characteristics of passive drag on competitive swimmers taking into account the impact of the free surface. This model solves the 3D incompressible Navier-Stokes equations using RNG k- ɛ turbulence closure. The volume of fluid (VOF) method is used to locate the free surface. The 3D virtual model is created by Computer Aided Industrial Design (CAID) software, Rhinoceros. Firstly, a specific posture of swimming is studied. The simulation results are in good agreement with the data from mannequin towing experiments. The effects of a swimmer's arms and legs positions on swimming performance are then studied. Finally, it is demonstrated that the present method is capable of simulating gliding near the free surface.
Impact of 3D root uptake on solute transport: a numerical study
NASA Astrophysics Data System (ADS)
Schröder, N.; Javaux, M.; Vanderborght, J.; Steffen, B.; Vereecken, H.
2011-12-01
Plant transpiration is an important component of the hydrological cycle. Through root water uptake, plants do not only affect the 3D soil water flow velocity distribution, but also solute movement in soil. This numerical study aims at investigating how solute fate is impacted by root uptake using the 3D biophysical model R-SWMS (Javaux et al., 2008). This model solves the Richards equation in 3D in the soil and the flow equation within the plant root xylem vessels. Furthermore, for solute transport simulations, the 3D particle tracker PARTRACE (Bechtold et al., 2011) was used. . We generated 3D virtual steady-state breakthrough curves (BTC) experiments in soils with transpiring plants. The averaged BTCs were then fitted with a 1D numerical flow model under steady-state conditions to obtain apparent CDE parameters. Two types of root architecture, a fibrous and a taprooted structure, were compared in virtual 3D experiments. The solute uptake type or the transpiration rate were also modified and we analyzed how these parameters affected apparent disperisivity and velocity profiles. Our simulation results show, that both, apparent velocity and dispersivity length are affected by water and solute root uptake. In addition, under high exclusion processes (slight or no active uptake), solute accumulates around roots and generates a long tailing to the breakthrough curves, which cannot be reproduced by 1D models that simulate root water uptake with solute exclusion. This observation may have an important impact on how to model pollutant mass transfer to groundwater at larger scales. Javaux, M., T. Schröder, J. Vanderborght, and H. Vereecken. 2008. Use of a three-dimensional detailed modeling approach for predicting root water uptake. Vadose Zone J. 7:1079-1088.doi: 10.2136/vzj2007.0115. Bechtold, M., S. Haber-Pohlmeier, J. Vanderborght, A. Pohlmeier, P.A. Ferre, and H. Vereecken. 2011. Near-surface solute redistribution during evaporation. Submitted to Geophys. Res. Lett
The Transient 3-D Transport Coupled Code TORT-TD/ATTICA3D for High-Fidelity Pebble-Bed HTGR Analyses
NASA Astrophysics Data System (ADS)
Seubert, Armin; Sureda, Antonio; Lapins, Janis; Bader, Johannes; Laurien, Eckart
2012-01-01
This article describes the 3D discrete ordinates-based coupled code system TORT-TD/ATTICA3D that aims at steady state and transient analyses of pebble-bed high-temperature gas cooled reactors. In view of increasing computing power, the application of time-dependent neutron transport methods becomes feasible for best estimate evaluations of safety margins. The calculation capabilities of TORT-TD/ATTICA3D are presented along with the coupling approach, with focus on the time-dependent neutron transport features of TORT-TD. Results obtained for the OECD/NEA/NSC PBMR-400 benchmark demonstrate the transient capabilities of TORT-TD/ATTICA3D.
Numerical Optimization Strategy for Determining 3D Flow Fields in Microfluidics
NASA Astrophysics Data System (ADS)
Eden, Alex; Sigurdson, Marin; Mezic, Igor; Meinhart, Carl
2015-11-01
We present a hybrid experimental-numerical method for generating 3D flow fields from 2D PIV experimental data. An optimization algorithm is applied to a theory-based simulation of an alternating current electrothermal (ACET) micromixer in conjunction with 2D PIV data to generate an improved representation of 3D steady state flow conditions. These results can be used to investigate mixing phenomena. Experimental conditions were simulated using COMSOL Multiphysics to solve the temperature and velocity fields, as well as the quasi-static electric fields. The governing equations were based on a theoretical model for ac electrothermal flows. A Nelder-Mead optimization algorithm was used to achieve a better fit by minimizing the error between 2D PIV experimental velocity data and numerical simulation results at the measurement plane. By applying this hybrid method, the normalized RMS velocity error between the simulation and experimental results was reduced by more than an order of magnitude. The optimization algorithm altered 3D fluid circulation patterns considerably, providing a more accurate representation of the 3D experimental flow field. This method can be generalized to a wide variety of flow problems. This research was supported by the Institute for Collaborative Biotechnologies through grant W911NF-09-0001 from the U.S. Army Research Office.
User Guide for the R5EXEC Coupling Interface in the RELAP5-3D Code
Forsmann, J. Hope; Weaver, Walter L.
2015-04-01
This report describes the R5EXEC coupling interface in the RELAP5-3D computer code from the users perspective. The information in the report is intended for users who want to couple RELAP5-3D to other thermal-hydraulic, neutron kinetics, or control system simulation codes.
Development of a GPU-Accelerated 3-D Full-Wave Code for Reflectometry Simulations
NASA Astrophysics Data System (ADS)
Reuther, K. S.; Kubota, S.; Feibush, E.; Johnson, I.
2013-10-01
1-D and 2-D full-wave codes used as synthetic diagnostics in microwave reflectometry are standard tools for understanding electron density fluctuations in fusion plasmas. The accuracy of the code depends on how well the wave properties along the ignored dimensions can be pre-specified or neglected. In a toroidal magnetic geometry, such assumptions are never strictly correct and ray tracing has shown that beam propagation is inherently a 3-D problem. Previously, we reported on the application of GPGPU's (General-Purpose computing on Graphics Processing Units) to a 2-D FDTD (Finite-Difference Time-Domain) code ported to utilize the parallel processing capabilities of the NVIDIA C870 and C1060. Here, we report on the development of a FDTD code for 3-D problems. Initial tests will use NVIDIA's M2070 GPU and concentrate on the launching and propagation of Gaussian beams in free space. If available, results using a plasma target will also be presented. Performance will be compared with previous generations of GPGPU cards as well as with NVIDIA's newest K20C GPU. Finally, the possibility of utilizing multiple GPGPU cards in a cluster environment or in a single node will also be discussed. Supported by U.S. DoE Grants DE-FG02-99-ER54527 and DE-AC02-09CH11466 and the DoE National Undergraduate Fusion Fellowship.
A 3-D Vortex Code for Parachute Flow Predictions: VIPAR Version 1.0
STRICKLAND, JAMES H.; HOMICZ, GREGORY F.; PORTER, VICKI L.; GOSSLER, ALBERT A.
2002-07-01
This report describes a 3-D fluid mechanics code for predicting flow past bluff bodies whose surfaces can be assumed to be made up of shell elements that are simply connected. Version 1.0 of the VIPAR code (Vortex Inflation PARachute code) is described herein. This version contains several first order algorithms that we are in the process of replacing with higher order ones. These enhancements will appear in the next version of VIPAR. The present code contains a motion generator that can be used to produce a large class of rigid body motions. The present code has also been fully coupled to a structural dynamics code in which the geometry undergoes large time dependent deformations. Initial surface geometry is generated from triangular shell elements using a code such as Patran and is written into an ExodusII database file for subsequent input into VIPAR. Surface and wake variable information is output into two ExodusII files that can be post processed and viewed using software such as EnSight{trademark}.
2D/1D approximations to the 3D neutron transport equation. II: Numerical comparisons
Kelley, B. W.; Collins, B.; Larsen, E. W.
2013-07-01
In a companion paper [1], (i) several new '2D/1D equations' are introduced as accurate approximations to the 3D Boltzmann transport equation, (ii) the simplest of these approximate equations is systematically discretized, and (iii) a theoretically stable iteration scheme is developed to solve the discrete equations. In this paper, numerical results are presented that confirm the theoretical predictions made in [1]. (authors)
PRONTO3D users` instructions: A transient dynamic code for nonlinear structural analysis
Attaway, S.W.; Mello, F.J.; Heinstein, M.W.; Swegle, J.W.; Ratner, J.A.; Zadoks, R.I.
1998-06-01
This report provides an updated set of users` instructions for PRONTO3D. PRONTO3D is a three-dimensional, transient, solid dynamics code for analyzing large deformations of highly nonlinear materials subjected to extremely high strain rates. This Lagrangian finite element program uses an explicit time integration operator to integrate the equations of motion. Eight-node, uniform strain, hexahedral elements and four-node, quadrilateral, uniform strain shells are used in the finite element formulation. An adaptive time step control algorithm is used to improve stability and performance in plasticity problems. Hourglass distortions can be eliminated without disturbing the finite element solution using either the Flanagan-Belytschko hourglass control scheme or an assumed strain hourglass control scheme. All constitutive models in PRONTO3D are cast in an unrotated configuration defined using the rotation determined from the polar decomposition of the deformation gradient. A robust contact algorithm allows for the impact and interaction of deforming contact surfaces of quite general geometry. The Smooth Particle Hydrodynamics method has been embedded into PRONTO3D using the contact algorithm to couple it with the finite element method.
A 3d particle simulation code for heavy ion fusion accelerator studies
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.
The future of 3D and video coding in mobile and the internet
NASA Astrophysics Data System (ADS)
Bivolarski, Lazar
2013-09-01
The current Internet success has already changed our social and economic world and is still continuing to revolutionize the information exchange. The exponential increase of amount and types of data that is currently exchanged on the Internet represents significant challenge for the design of future architectures and solutions. This paper reviews the current status and trends in the design of solutions and research activities in the future Internet from point of view of managing the growth of bandwidth requirements and complexity of the multimedia that is being created and shared. Outlines the challenges that are present before the video coding and approaches to the design of standardized media formats and protocols while considering the expected convergence of multimedia formats and exchange interfaces. The rapid growth of connected mobile devices adds to the current and the future challenges in combination with the expected, in near future, arrival of multitude of connected devices. The new Internet technologies connecting the Internet of Things with wireless visual sensor networks and 3D virtual worlds requires conceptually new approaches of media content handling from acquisition to presentation in the 3D Media Internet. Accounting for the entire transmission system properties and enabling adaptation in real-time to context and content throughout the media proceeding path will be paramount in enabling the new media architectures as well as the new applications and services. The common video coding formats will need to be conceptually redesigned to allow for the implementation of the necessary 3D Media Internet features.
3-D Numerical Modeling of MHD Flows in Variable Magnetic Field
NASA Astrophysics Data System (ADS)
Abdullina, K. I.; Bogovalov, S. V.
3-D numerical simulation of the liquid metal flow affected by the electromagnetic field in the magnetohydrodynamic (MHD) devices is performed. Software package ANSYS has been used for the numerical calculations. The non-stationary problem has been solved taking into account the influence of the metal flow on the electromagnetic field and nonlinear magnetic permeability of the ferromagnetic cores. Simplified calculations with constant magnetic permeability of the ferromagnetic cores have been performed as well. Comparison of these calculations shows that the simulation of the MHD pump can be performed in the linear approximation. The pump performance curve has been derived in this approximation.
FURN3D: A computer code for radiative heat transfer in pulverized coal furnaces
Ahluwalia, R.K.; Im, K.H.
1992-08-01
A computer code FURN3D has been developed for assessing the impact of burning different coals on heat absorption pattern in pulverized coal furnaces. The code is unique in its ability to conduct detailed spectral calculations of radiation transport in furnaces fully accounting for the size distributions of char, soot and ash particles, ash content, and ash composition. The code uses a hybrid technique of solving the three-dimensional radiation transport equation for absorbing, emitting and anisotropically scattering media. The technique achieves an optimal mix of computational speed and accuracy by combining the discrete ordinate method (S[sub 4]), modified differential approximation (MDA) and P, approximation in different range of optical thicknesses. The code uses spectroscopic data for estimating the absorption coefficients of participating gases C0[sub 2], H[sub 2]0 and CO. It invokes Mie theory for determining the extinction and scattering coefficients of combustion particulates. The optical constants of char, soot and ash are obtained from dispersion relations derived from reflectivity, transmissivity and extinction measurements. A control-volume formulation is adopted for determining the temperature field inside the furnace. A simple char burnout model is employed for estimating heat release and evolution of particle size distribution. The code is written in Fortran 77, has modular form, and is machine-independent. The computer memory required by the code depends upon the number of grid points specified and whether the transport calculations are performed on spectral or gray basis.
FURN3D: A computer code for radiative heat transfer in pulverized coal furnaces
Ahluwalia, R.K.; Im, K.H.
1992-08-01
A computer code FURN3D has been developed for assessing the impact of burning different coals on heat absorption pattern in pulverized coal furnaces. The code is unique in its ability to conduct detailed spectral calculations of radiation transport in furnaces fully accounting for the size distributions of char, soot and ash particles, ash content, and ash composition. The code uses a hybrid technique of solving the three-dimensional radiation transport equation for absorbing, emitting and anisotropically scattering media. The technique achieves an optimal mix of computational speed and accuracy by combining the discrete ordinate method (S{sub 4}), modified differential approximation (MDA) and P, approximation in different range of optical thicknesses. The code uses spectroscopic data for estimating the absorption coefficients of participating gases C0{sub 2}, H{sub 2}0 and CO. It invokes Mie theory for determining the extinction and scattering coefficients of combustion particulates. The optical constants of char, soot and ash are obtained from dispersion relations derived from reflectivity, transmissivity and extinction measurements. A control-volume formulation is adopted for determining the temperature field inside the furnace. A simple char burnout model is employed for estimating heat release and evolution of particle size distribution. The code is written in Fortran 77, has modular form, and is machine-independent. The computer memory required by the code depends upon the number of grid points specified and whether the transport calculations are performed on spectral or gray basis.
Mechanical Modelling of Pultrusion Process: 2D and 3D Numerical Approaches
NASA Astrophysics Data System (ADS)
Baran, Ismet; Hattel, Jesper H.; Akkerman, Remko; Tutum, Cem C.
2015-02-01
The process induced variations such as residual stresses and distortions are a critical issue in pultrusion, since they affect the structural behavior as well as the mechanical properties and geometrical precision of the final product. In order to capture and investigate these variations, a mechanical analysis should be performed. In the present work, the two dimensional (2D) quasi-static plane strain mechanical model for the pultrusion of a thick square profile developed by the authors is further improved using generalized plane strain elements. In addition to that, a more advanced 3D thermo-chemical-mechanical analysis is carried out using 3D quadratic elements which is a novel application for the numerical modelling of the pultrusion process. It is found that the 2D mechanical models give relatively reasonable and accurate stress and displacement evolutions in the transverse direction as compared to the 3D model. Moreover, the generalized plane strain model predicts the longitudinal process induced stresses more similar to the ones calculated in the 3D model as compared with the plane strain model.
Validation of CATHARE 3D Code Against UPTF TRAM C3 Transients
NASA Astrophysics Data System (ADS)
Glantz, Tony; Freitas, Roberto
Within the nuclear reactor safety analysis, one of the events that could potentially lead to a re-criticality accident in case of a Small Break Loss of Coolant Accident (SBLOCA) in a Pressurized Water Reactor (PWR) is a boron dilution scenario followed by a coolant mixing transient. Some UPTF experiments can be interpreted as generic boron dilution experiments. In fact, the UPTF experiments were originally designed to conduct separate effects studies focused on multi-dimensional thermal hydraulic phenomena. However, in the case of experimental program TRAM, some studies are realized on the boron mixing: tests C3. Some of these tests have been used for the validation and assessment of the 3D module of CATHARE code. Results are very satisfying; CATHARE 3D code is able to reproduce correctly the main features of the UPTF TRAM C3 tests, the temperature mixing in the cold leg, the formation of a strong stratification in the upper downcomer, the perfect mixing temperature in the lower downcomer and the strong stratification in the lower plenum. These results are also compared with the CFX5 and TRIO-U codes results on these tests.
3D and 4D Simulations of the Dynamics of the Radiation Belts using VERB code
NASA Astrophysics Data System (ADS)
Shprits, Yuri; Kellerman, Adam; Drozdov, Alexander; Orlova, Ksenia
2015-04-01
Modeling and understanding of ring current and higher energy radiation belts has been a grand challenge since the beginning of the space age. In this study we show long term simulations with a 3D VERB code of modeling the radiation belts with boundary conditions derived from observations around geosynchronous orbit. We also present 4D VERB simulations that include convective transport, radial diffusion, pitch angle scattering and local acceleration. We show that while lower energy radial transport is dominated by the convection and higher energy transport is dominated by the diffusive radial transport. We also show there exists an intermediate range of energies for electrons for which both processes work simultaneously.
A 3D numerical study of antimicrobial persistence in heterogeneous multi-species biofilms.
Zhao, Jia; Shen, Ya; Haapasalo, Markus; Wang, Zhejun; Wang, Qi
2016-03-01
We develop a 3D hydrodynamic model to investigate the mechanism of antimicrobial persistence in a multi-species oral biofilm and its recovery after being treated by bisbiguanide chlorhexidine gluconate (CHX). In addition to the hydrodynamic transport in the spatially heterogeneous biofilm, the model also includes mechanisms of solvent-biomass interaction, bacterial phenotype conversion, and bacteria-drug interaction. A numerical solver for the model is developed using a second order numerical scheme in 3D space and time and implemented on GPUs for high-performance computing. The model is calibrated against a set of experimental data obtained using confocal laser scan microscopy (CLSM) on multi-species oral biofilms, where a quantitative agreement is reached. Our numerical results reveal that quorum sensing molecules and growth factors in this model are instrumental in biofilm formation and recovery after the antimicrobial treatment. In particular, we show that (i) young biofilms are more susceptible to the antimicrobial treatment than the mature ones, (ii) this phenomenon is strongly correlated with volume fractions of the persister and EPS in the biofilm being treated. This suggests that antimicrobial treatment should be best administered to biofilms earlier before they mature to produce a thick protective EPS layer. In addition, the numerical study also indicates that an antimicrobial effect can be achieved should a proper mechanism be devised to minimize the conversion of susceptible bacteria to persisters during and even after the treatment. PMID:26739374
Dong, Xiaoqing; Fang, Yiliang; Wang, Kejing; Zhu, Lijuan; Wang, Ke; Huang, Tao
2016-01-01
With the development of new technologies in transcriptome and epigenetics, RNAs have been identified to play more and more important roles in life processes. Consequently, various methods have been proposed to assess the biological functions of RNAs and thus classify them functionally, among which comparative study of RNA structures is perhaps the most important one. To measure the structural similarity of RNAs and classify them, we propose a novel three dimensional (3D) graphical representation of RNA secondary structure, in which an RNA secondary structure is first transformed into a characteristic sequence based on chemical property of nucleic acids; a dynamic 3D graph is then constructed for the characteristic sequence; and lastly a numerical characterization of the 3D graph is used to represent the RNA secondary structure. We tested our algorithm on three datasets: (1) Dataset I consisting of nine RNA secondary structures of viruses, (2) Dataset II consisting of complex RNA secondary structures including pseudo-knots, and (3) Dataset III consisting of 18 non-coding RNA families. We also compare our method with other nine existing methods using Dataset II and III. The results demonstrate that our method is better than other methods in similarity measurement and classification of RNA secondary structures. PMID:27213271
Zhang, Yi; Huang, Haiyun; Dong, Xiaoqing; Fang, Yiliang; Wang, Kejing; Zhu, Lijuan; Wang, Ke; Huang, Tao; Yang, Jialiang
2016-01-01
With the development of new technologies in transcriptome and epigenetics, RNAs have been identified to play more and more important roles in life processes. Consequently, various methods have been proposed to assess the biological functions of RNAs and thus classify them functionally, among which comparative study of RNA structures is perhaps the most important one. To measure the structural similarity of RNAs and classify them, we propose a novel three dimensional (3D) graphical representation of RNA secondary structure, in which an RNA secondary structure is first transformed into a characteristic sequence based on chemical property of nucleic acids; a dynamic 3D graph is then constructed for the characteristic sequence; and lastly a numerical characterization of the 3D graph is used to represent the RNA secondary structure. We tested our algorithm on three datasets: (1) Dataset I consisting of nine RNA secondary structures of viruses, (2) Dataset II consisting of complex RNA secondary structures including pseudo-knots, and (3) Dataset III consisting of 18 non-coding RNA families. We also compare our method with other nine existing methods using Dataset II and III. The results demonstrate that our method is better than other methods in similarity measurement and classification of RNA secondary structures. PMID:27213271
Numerical simulation of internal and external inviscid and viscous 3-D flow fields
NASA Astrophysics Data System (ADS)
Leicher, Stefan
1986-11-01
A numerical method for solving the 3-D Euler equations in geometrical complex domains was developed. The approach divides the computational space into multiple blocks whose structure follows the natural lines of the conficuration. A systematic, multi-block grid generation scheme is used to produce the grid. The flow solutions are obtained by solving the Euler equations by a finite volume discretization and a Runge-Kutta time stepping scheme. The main advantage of this method is the applicability to complex geometries, for example complete aircraft configurations including wing, fuselage, canard and tail. The coupling with a 3-D boundary layer method allows to account for viscous effects. Another application for the method was the simulation of flows in the presence of a propeller.
Numerical solution of 3-D magnetotelluric using vector finite element method
NASA Astrophysics Data System (ADS)
Prihantoro, Rudy; Sutarno, Doddy; Nurhasan
2015-09-01
Magnetotelluric (MT) is a passive electromagnetic (EM) method which measure natural variations of electric and magnetic vector fields at the Earth surface to map subsurface electrical conductivity/resistivity structure. In this study, we obtained numerical solution of three-dimensional (3-D) MT using vector finite element method by solving second order Maxwell differential equation describing diffusion of plane wave through the conductive earth. Rather than the nodes of the element, the edges of the element is used as a vector basis to overcome the occurrence of nonphysical solutions that usually faced by scalar (node based) finite element method. Electric vector fields formulation was used and the resulting system of equation was solved using direct solution method to obtain the electric vector field distribution throughout the earth resistivity model structure. The resulting MT response functions was verified with 1-D layered Earth and 3-D2 COMMEMI outcropping structure. Good agreement is achieved for both structure models.
The 3D modeling of high numerical aperture imaging in thin films
NASA Technical Reports Server (NTRS)
Flagello, D. G.; Milster, Tom
1992-01-01
A modelling technique is described which is used to explore three dimensional (3D) image irradiance distributions formed by high numerical aperture (NA is greater than 0.5) lenses in homogeneous, linear films. This work uses a 3D modelling approach that is based on a plane-wave decomposition in the exit pupil. Each plane wave component is weighted by factors due to polarization, aberration, and input amplitude and phase terms. This is combined with a modified thin-film matrix technique to derive the total field amplitude at each point in a film by a coherent vector sum over all plane waves. Then the total irradiance is calculated. The model is used to show how asymmetries present in the polarized image change with the influence of a thin film through varying degrees of focus.
FDFD: A 3D Finite-Difference Frequency-Domain Code for Electromagnetic Induction Tomography
NASA Astrophysics Data System (ADS)
Champagne, Nathan J.; Berryman, James G.; Buettner, H. Michael
2001-07-01
A new 3D code for electromagnetic induction tomography with intended applications to environmental imaging problems has been developed. The approach consists of calculating the fields within a volume using an implicit finite-difference frequency-domain formulation. The volume is terminated by an anisotropic perfectly matched layer region that simulates an infinite domain by absorbing outgoing waves. Extensive validation of this code has been done using analytical and semianalytical results from other codes, and some of those results are presented in this paper. The new code is written in Fortran 90 and is designed to be easily parallelized. Finally, an adjoint field method of data inversion, developed in parallel for solving the fully nonlinear inverse problem for electrical conductivity imaging (e.g., for mapping underground conducting plumes), uses this code to provide solvers for both forward and adjoint fields. Results obtained from this inversion method for high-contrast media are encouraging and provide a significant improvement over those obtained from linearized inversion methods.
Automated design of coupled RF cavities using 2-D and 3-D codes
Smith, Peter; Christiansen, D. W.; Greninger, P. T.; Spalek, G.
2001-01-01
Coupled RF cavities in the Accelerator Production of Tritium Project have been designed using a procedure in which a 2-D code (CCT) searches for a design that meets frequency and coupling requirements, while a 3-D code (HFSS) is used to obtain empirical factors used by CCT to characterize the coupling slot between cavities. Using assumed values of the empirical factors, CCT runs the Superfish code iteratively to solve for a trial cavity design that has a specified frequency and coupling. The frequency shifts and the coupling constant k of the slot are modeled in CCT using a perturbation theory, the results of which are adjusted using the empirical factors. Given a trial design, HFSS is run using periodic boundary conditions to obtain a mode spectrum. The mode spectrum is processed using the DISPER code to obtain values of the coupling and the frequencies with slots. These results are used to calculate a new set of empirical factors, which are fed back into CCT for another design iteration. Cold models have been fabricated and tested to validate the codes, and results will be presented.
3-D localization of gamma ray sources with coded apertures for medical applications
NASA Astrophysics Data System (ADS)
Kaissas, I.; Papadimitropoulos, C.; Karafasoulis, K.; Potiriadis, C.; Lambropoulos, C. P.
2015-09-01
Several small gamma cameras for radioguided surgery using CdTe or CdZnTe have parallel or pinhole collimators. Coded aperture imaging is a well-known method for gamma ray source directional identification, applied in astrophysics mainly. The increase in efficiency due to the substitution of the collimators by the coded masks renders the method attractive for gamma probes used in radioguided surgery. We have constructed and operationally verified a setup consisting of two CdTe gamma cameras with Modified Uniform Redundant Array (MURA) coded aperture masks of rank 7 and 19 and a video camera. The 3-D position of point-like radioactive sources is estimated via triangulation using decoded images acquired by the gamma cameras. We have also developed code for both fast and detailed simulations and we have verified the agreement between experimental results and simulations. In this paper we present a simulation study for the spatial localization of two point sources using coded aperture masks with rank 7 and 19.
NASA Astrophysics Data System (ADS)
Krebs, Isabel; Jardin, Stephen C.; Igochine, Valentin; Guenter, Sibylle; Hoelzl, Matthias; ASDEX Upgrade Team
2014-10-01
We study sawtooth reconnection in ASDEX Upgrade tokamak plasmas by means of 3D non-linear two-fluid MHD simulations in toroidal geometry using the high-order finite element code M3D-C1. Parameters and equilibrium of the simulations are based on typical sawtoothing ASDEX Upgrade discharges. The simulation results are compared to features of the experimental observations such as the sawtooth crash time and frequency, the evolution of the safety factor profile and the 3D evolution of the temperature. 2D ECE imaging measurements during sawtooth crashes in ASDEX Upgrade indicate that the heat is transported out of the core through a narrow poloidally localized region. We investigate if incomplete sawtooth reconnection can be seen in the simulations which is suggested by soft X-ray tomography measurements in ASDEX Upgrade showing that an (m = 1, n = 1) perturbation is typically observed to survive the sawtooth crash and approximately maintain its radial position.
Assessment of 3D Codes for Predicting Liner Attenuation in Flow Ducts
NASA Technical Reports Server (NTRS)
Watson, W. R.; Nark, D. M.; Jones, M. G.
2008-01-01
This paper presents comparisons of seven propagation codes for predicting liner attenuation in ducts with flow. The selected codes span the spectrum of methods available (finite element, parabolic approximation, and pseudo-time domain) and are collectively representative of the state-of-art in the liner industry. These codes are included because they have two-dimensional and three-dimensional versions and can be exported to NASA's Columbia Supercomputer. The basic assumptions, governing differential equations, boundary conditions, and numerical methods underlying each code are briefly reviewed and an assessment is performed based on two predefined metrics. The two metrics used in the assessment are the accuracy of the predicted attenuation and the amount of wall clock time to predict the attenuation. The assessment is performed over a range of frequencies, mean flow rates, and grazing flow liner impedances commonly used in the liner industry. The primary conclusions of the study are (1) predicted attenuations are in good agreement for rigid wall ducts, (2) the majority of codes compare well to each other and to approximate results from mode theory for soft wall ducts, (3) most codes compare well to measured data on a statistical basis, (4) only the finite element codes with cubic Hermite polynomials capture extremely large attenuations, and (5) wall clock time increases by an order of magnitude or more are observed for a three-dimensional code relative to the corresponding two-dimensional version of the same code.
Implementation of the 3D edge plasma code EMC3-EIRENE on NSTX
Lore, J. D.; Canik, J. M.; Feng, Y.; Ahn, J. -W.; Maingi, R.; Soukhanovskii, V.
2012-05-09
The 3D edge transport code EMC3-EIRENE has been applied for the first time to the NSTX spherical tokamak. A new disconnected double null grid has been developed to allow the simulation of plasma where the radial separation of the inner and outer separatrix is less than characteristic widths (e.g. heat flux width) at the midplane. Modelling results are presented for both an axisymmetric case and a case where 3D magnetic field is applied in an n = 3 configuration. In the vacuum approximation, the perturbed field consists of a wide region of destroyed flux surfaces and helical lobes which are a mixture of long and short connection length field lines formed by the separatrix manifolds. This structure is reflected in coupled 3D plasma fluid (EMC3) and kinetic neutral particle (EIRENE) simulations. The helical lobes extending inside of the unperturbed separatrix are filled in by hot plasma from the core. The intersection of the lobes with the divertor results in a striated flux footprint pattern on the target plates. As a result, profiles of divertor heat and particle fluxes are compared with experimental data, and possible sources of discrepancy are discussed.
Implementation of the 3D edge plasma code EMC3-EIRENE on NSTX
Lore, J. D.; Canik, J. M.; Feng, Y.; Ahn, J. -W.; Maingi, R.; Soukhanovskii, V.
2012-05-09
The 3D edge transport code EMC3-EIRENE has been applied for the first time to the NSTX spherical tokamak. A new disconnected double null grid has been developed to allow the simulation of plasma where the radial separation of the inner and outer separatrix is less than characteristic widths (e.g. heat flux width) at the midplane. Modelling results are presented for both an axisymmetric case and a case where 3D magnetic field is applied in an n = 3 configuration. In the vacuum approximation, the perturbed field consists of a wide region of destroyed flux surfaces and helical lobes which aremore » a mixture of long and short connection length field lines formed by the separatrix manifolds. This structure is reflected in coupled 3D plasma fluid (EMC3) and kinetic neutral particle (EIRENE) simulations. The helical lobes extending inside of the unperturbed separatrix are filled in by hot plasma from the core. The intersection of the lobes with the divertor results in a striated flux footprint pattern on the target plates. As a result, profiles of divertor heat and particle fluxes are compared with experimental data, and possible sources of discrepancy are discussed.« less
Numerical and measured data from the 3D salt canopy physical modeling project
Bradley, C.; House, L.; Fehler, M.; Pearson, J.; TenCate, J.; Wiley, R.
1997-11-01
The evolution of salt structures in the Gulf of Mexico have been shown to provide a mechanism for the trapping of significant hydrocarbon reserves. Most of these structures have complex geometries relative to the surrounding sedimentary layers. This aspect in addition to high velocities within the salt tend to scatter and defocus seismic energy and make imaging of subsalt lithology extremely difficult. An ongoing program the SEG/EAEG modeling project (Aminzadeh et al. 1994a: Aminzadeh et al. 1994b: Aminzadeh et al. 1995), and a follow-up project funded as part of the Advanced Computational Technology Initiative (ACTI) (House et al. 1996) have sought to investigate problems with imaging beneath complex salt structures using numerical modeling and more recently, construction of a physical model patterned after the numerical subsalt model (Wiley and McKnight. 1996). To date, no direct comparison of the numerical and physical aspects of these models has been attempted. We present the results of forward modeling a numerical realization of the 3D salt canopy physical model with the French Petroleum Institute (IFP) acoustic finite difference algorithm used in the numerical subsalt tests. We compare the results from the physical salt canopy model, the acoustic modeling of the physical/numerical model and the original numerical SEG/EAEG Salt Model. We will be testing the sensitivity of migration to the presence of converted shear waves and acquisition geometry.
2D and 3D numerical models on compositionally buoyant diapirs in the mantle wedge
NASA Astrophysics Data System (ADS)
Hasenclever, Jörg; Morgan, Jason Phipps; Hort, Matthias; Rüpke, Lars H.
2011-11-01
We present 2D and 3D numerical model calculations that focus on the physics of compositionally buoyant diapirs rising within a mantle wedge corner flow. Compositional buoyancy is assumed to arise from slab dehydration during which water-rich volatiles enter the mantle wedge and form a wet, less dense boundary layer on top of the slab. Slab dehydration is prescribed to occur in the 80-180 km deep slab interval, and the water transport is treated as a diffusion-like process. In this study, the mantle's rheology is modeled as being isoviscous for the benefit of easier-to-interpret feedbacks between water migration and buoyant viscous flow of the mantle. We use a simple subduction geometry that does not change during the numerical calculation. In a large set of 2D calculations we have identified that five different flow regimes can form, in which the position, number, and formation time of the diapirs vary as a function of four parameters: subduction angle, subduction rate, water diffusivity (mobility), and mantle viscosity. Using the same numerical method and numerical resolution we also conducted a suite of 3D calculations for 16 selected parameter combinations. Comparing the 2D and 3D results for the same model parameters reveals that the 2D models can only give limited insights into the inherently 3D problem of mantle wedge diapirism. While often correctly predicting the position and onset time of the first diapir(s), the 2D models fail to capture the dynamics of diapir ascent as well as the formation of secondary diapirs that result from boundary layer perturbations caused by previous diapirs. Of greatest importance for physically correct results is the numerical resolution in the region where diapirs nucleate, which must be high enough to accurately capture the growth of the thin wet boundary layer on top of the slab and, subsequently, the formation, morphology, and ascent of diapirs. Here 2D models can be very useful to quantify the required resolution, which we
3D thermo-chemical-mechanical simulation of power ramps with ALCYONE fuel code
NASA Astrophysics Data System (ADS)
Baurens, B.; Sercombe, J.; Riglet-Martial, C.; Desgranges, L.; Trotignon, L.; Maugis, P.
2014-09-01
This paper presents the coupling of the fuel performance code ALCYONE with the thermochemical code ANGE and its application to Iodine-Stress Corrosion Cracking (I-SCC). The coupling is illustrated by a 3D simulation of a power ramp. The release of chemically active gases (CsI(g), Tex(1
3-D model-based frame interpolation for distributed video coding of static scenes.
Maitre, Matthieu; Guillemot, Christine; Morin, Luce
2007-05-01
This paper addresses the problem of side information extraction for distributed coding of videos captured by a camera moving in a 3-D static environment. Examples of targeted applications are augmented reality, remote-controlled robots operating in hazardous environments, or remote exploration by drones. It explores the benefits of the structure-from-motion paradigm for distributed coding of this type of video content. Two interpolation methods constrained by the scene geometry, based either on block matching along epipolar lines or on 3-D mesh fitting, are first developed. These techniques are based on a robust algorithm for sub-pel matching of feature points, which leads to semi-dense correspondences between key frames. However, their rate-distortion (RD) performances are limited by misalignments between the side information and the actual Wyner-Ziv (WZ) frames due to the assumption of linear motion between key frames. To cope with this problem, two feature point tracking techniques are introduced, which recover the camera parameters of the WZ frames. A first technique, in which the frames remain encoded separately, performs tracking at the decoder and leads to significant RD performance gains. A second technique further improves the RD performances by allowing a limited tracking at the encoder. As an additional benefit, statistics on tracks allow the encoder to adapt the key frame frequency to the video motion content. PMID:17491456
Radiation Coupling with the FUN3D Unstructured-Grid CFD Code
NASA Technical Reports Server (NTRS)
Wood, William A.
2012-01-01
The HARA radiation code is fully-coupled to the FUN3D unstructured-grid CFD code for the purpose of simulating high-energy hypersonic flows. The radiation energy source terms and surface heat transfer, under the tangent slab approximation, are included within the fluid dynamic ow solver. The Fire II flight test, at the Mach-31 1643-second trajectory point, is used as a demonstration case. Comparisons are made with an existing structured-grid capability, the LAURA/HARA coupling. The radiative surface heat transfer rates from the present approach match the benchmark values within 6%. Although radiation coupling is the focus of the present work, convective surface heat transfer rates are also reported, and are seen to vary depending upon the choice of mesh connectivity and FUN3D ux reconstruction algorithm. On a tetrahedral-element mesh the convective heating matches the benchmark at the stagnation point, but under-predicts by 15% on the Fire II shoulder. Conversely, on a mixed-element mesh the convective heating over-predicts at the stagnation point by 20%, but matches the benchmark away from the stagnation region.
NASA Astrophysics Data System (ADS)
Inogamov, Nail A.; Zhakhovsky, Vasily V.
2016-02-01
There are many important applications in which the ultrashort diffraction-limited and therefore tightly focused laser pulses irradiates metal films mounted on dielectric substrate. Here we present the detailed picture of laser peeling and 3D structure formation of the thin (relative to a depth of a heat affected zone in the bulk targets) gold films on glass substrate. The underlying physics of such diffraction-limited laser peeling was not well understood previously. Our approach is based on a physical model which takes into consideration the new calculations of the two-temperature (2T) equation of state (2T EoS) and the two-temperature transport coefficients together with the coupling parameter between electron and ion subsystems. The usage of the 2T EoS and the kinetic coefficients is required because absorption of an ultrashort pulse with duration of 10-1000 fs excites electron subsystem of metal and transfers substance into the 2T state with hot electrons (typical electron temperatures 1-3 eV) and much colder ions. It is shown that formation of submicrometer-sized 3D structures is a result of the electron-ion energy transfer, melting, and delamination of film from substrate under combined action of electron and ion pressures, capillary deceleration of the delaminated liquid metal or semiconductor, and ultrafast freezing of molten material. We found that the freezing is going in non-equilibrium regime with strongly overcooled liquid phase. In this case the Stefan approximation is non-applicable because the solidification front speed is limited by the diffusion rate of atoms in the molten material. To solve the problem we have developed the 2T Lagrangian code including all this reach physics in. We also used the high-performance combined Monte- Carlo and molecular dynamics code for simulation of surface 3D nanostructuring at later times after completion of electron-ion relaxation.
NASA Astrophysics Data System (ADS)
Wang, Lijun; Huang, Xiaolong; Jia, Shenli; Deng, Jie; Qian, Zhonghao; Shi, Zongqian; Schellenkens, H.; Godechot, X.
2015-06-01
A time-dependent 3D numerical model considering anode evaporation is developed for the high current vacuum arc (VA) under a realistic spatial magnetic field. The simulation work contains steady state 3D numerical simulation of high current VA considering anode evaporation at nine discrete moments of first half wave of 50 Hz AC current, transient numerical simulation of anode activity, and realistic spatial magnetic field calculation of commercial cup-shaped electrodes. In the simulation, contact opening and arc diffusion processes are also considered. Due to the effect of electrode slots, the simulation results of magnetic field and temperature of anode plate exhibit six leaves shape (SLS). During 6-8 ms, the strong evaporation of anode surface seriously influence the parameter distributions of VA. Ions emitted from anode penetrate into arc column and the axial velocity distribution on the anode side exhibits SLS. The ions emitted from anode surface have the same temperature with anode surface, which cool the arc plasma and lead to a relative low temperature area formed. The seriously evaporation of anode leads to the accumulation of ions near the anode, and then the current density is more uniform.
Development and preliminary verification of the 3D core neutronic code: COCO
Lu, H.; Mo, K.; Li, W.; Bai, N.; Li, J.
2012-07-01
As the recent blooming economic growth and following environmental concerns (China)) is proactively pushing forward nuclear power development and encouraging the tapping of clean energy. Under this situation, CGNPC, as one of the largest energy enterprises in China, is planning to develop its own nuclear related technology in order to support more and more nuclear plants either under construction or being operation. This paper introduces the recent progress in software development for CGNPC. The focus is placed on the physical models and preliminary verification results during the recent development of the 3D Core Neutronic Code: COCO. In the COCO code, the non-linear Green's function method is employed to calculate the neutron flux. In order to use the discontinuity factor, the Neumann (second kind) boundary condition is utilized in the Green's function nodal method. Additionally, the COCO code also includes the necessary physical models, e.g. single-channel thermal-hydraulic module, burnup module, pin power reconstruction module and cross-section interpolation module. The preliminary verification result shows that the COCO code is sufficient for reactor core design and analysis for pressurized water reactor (PWR). (authors)
Numerical investigation of band gaps in 3D printed cantilever-in-mass metamaterials
NASA Astrophysics Data System (ADS)
Qureshi, Awais; Li, Bing; Tan, K. T.
2016-06-01
In this research, the negative effective mass behavior of elastic/mechanical metamaterials is exhibited by a cantilever-in-mass structure as a proposed design for creating frequency stopping band gaps, based on local resonance of the internal structure. The mass-in-mass unit cell model is transformed into a cantilever-in-mass model using the Bernoulli-Euler beam theory. An analytical model of the cantilever-in-mass structure is derived and the effects of geometrical dimensions and material parameters to create frequency band gaps are examined. A two-dimensional finite element model is created to validate the analytical results, and excellent agreement is achieved. The analytical model establishes an easily tunable metamaterial design to realize wave attenuation based on locally resonant frequency. To demonstrate feasibility for 3D printing, the analytical model is employed to design and fabricate 3D printable mechanical metamaterial. A three-dimensional numerical experiment is performed using COMSOL Multiphysics to validate the wave attenuation performance. Results show that the cantilever-in-mass metamaterial is capable of mitigating stress waves at the desired resonance frequency. Our study successfully presents the use of one constituent material to create a 3D printed cantilever-in-mass metamaterial with negative effective mass density for stress wave mitigation purposes.
Numerical investigation of band gaps in 3D printed cantilever-in-mass metamaterials.
Qureshi, Awais; Li, Bing; Tan, K T
2016-01-01
In this research, the negative effective mass behavior of elastic/mechanical metamaterials is exhibited by a cantilever-in-mass structure as a proposed design for creating frequency stopping band gaps, based on local resonance of the internal structure. The mass-in-mass unit cell model is transformed into a cantilever-in-mass model using the Bernoulli-Euler beam theory. An analytical model of the cantilever-in-mass structure is derived and the effects of geometrical dimensions and material parameters to create frequency band gaps are examined. A two-dimensional finite element model is created to validate the analytical results, and excellent agreement is achieved. The analytical model establishes an easily tunable metamaterial design to realize wave attenuation based on locally resonant frequency. To demonstrate feasibility for 3D printing, the analytical model is employed to design and fabricate 3D printable mechanical metamaterial. A three-dimensional numerical experiment is performed using COMSOL Multiphysics to validate the wave attenuation performance. Results show that the cantilever-in-mass metamaterial is capable of mitigating stress waves at the desired resonance frequency. Our study successfully presents the use of one constituent material to create a 3D printed cantilever-in-mass metamaterial with negative effective mass density for stress wave mitigation purposes. PMID:27329828
NuSol - Numerical solver for the 3D stationary nuclear Schrödinger equation
NASA Astrophysics Data System (ADS)
Graen, Timo; Grubmüller, Helmut
2016-01-01
The classification of short hydrogen bonds depends on several factors including the shape and energy spacing between the nuclear eigenstates of the hydrogen. Here, we describe the NuSol program in which three classes of algorithms were implemented to solve the 1D, 2D and 3D time independent nuclear Schrödinger equation. The Schrödinger equation was solved using the finite differences based Numerov's method which was extended to higher dimensions, the more accurate pseudo-spectral Chebyshev collocation method and the sinc discrete variable representation by Colbert and Miller. NuSol can be applied to solve the Schrödinger equation for arbitrary analytical or numerical potentials with focus on nuclei bound by the potential of their molecular environment. We validated the methods against literature values for the 2D Henon-Heiles potential, the 3D linearly coupled sextic oscillators and applied them to study hydrogen bonding in the malonaldehyde derivate 4-cyano-2,2,6,6-tetramethyl-3,5-heptanedione. With NuSol, the extent of nuclear delocalization in a given molecular potential can directly be calculated without relying on linear reaction coordinates in 3D molecular space.
Numerical investigation of band gaps in 3D printed cantilever-in-mass metamaterials
Qureshi, Awais; Li, Bing; Tan, K. T.
2016-01-01
In this research, the negative effective mass behavior of elastic/mechanical metamaterials is exhibited by a cantilever-in-mass structure as a proposed design for creating frequency stopping band gaps, based on local resonance of the internal structure. The mass-in-mass unit cell model is transformed into a cantilever-in-mass model using the Bernoulli-Euler beam theory. An analytical model of the cantilever-in-mass structure is derived and the effects of geometrical dimensions and material parameters to create frequency band gaps are examined. A two-dimensional finite element model is created to validate the analytical results, and excellent agreement is achieved. The analytical model establishes an easily tunable metamaterial design to realize wave attenuation based on locally resonant frequency. To demonstrate feasibility for 3D printing, the analytical model is employed to design and fabricate 3D printable mechanical metamaterial. A three-dimensional numerical experiment is performed using COMSOL Multiphysics to validate the wave attenuation performance. Results show that the cantilever-in-mass metamaterial is capable of mitigating stress waves at the desired resonance frequency. Our study successfully presents the use of one constituent material to create a 3D printed cantilever-in-mass metamaterial with negative effective mass density for stress wave mitigation purposes. PMID:27329828
Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model
Baudron, Anne-Marie; Riahi, Mohamed Kamel; Salomon, Julien
2014-12-15
In this paper we present a time-parallel algorithm for the 3D neutrons calculation of a transient model in a nuclear reactor core. The neutrons calculation consists in numerically solving the time dependent diffusion approximation equation, which is a simplified transport equation. The numerical resolution is done with finite elements method based on a tetrahedral meshing of the computational domain, representing the reactor core, and time discretization is achieved using a θ-scheme. The transient model presents moving control rods during the time of the reaction. Therefore, cross-sections (piecewise constants) are taken into account by interpolations with respect to the velocity of the control rods. The parallelism across the time is achieved by an adequate use of the parareal in time algorithm to the handled problem. This parallel method is a predictor corrector scheme that iteratively combines the use of two kinds of numerical propagators, one coarse and one fine. Our method is made efficient by means of a coarse solver defined with large time step and fixed position control rods model, while the fine propagator is assumed to be a high order numerical approximation of the full model. The parallel implementation of our method provides a good scalability of the algorithm. Numerical results show the efficiency of the parareal method on large light water reactor transient model corresponding to the Langenbuch–Maurer–Werner benchmark.
Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model
NASA Astrophysics Data System (ADS)
Baudron, Anne-Marie; Lautard, Jean-Jacques; Maday, Yvon; Riahi, Mohamed Kamel; Salomon, Julien
2014-12-01
In this paper we present a time-parallel algorithm for the 3D neutrons calculation of a transient model in a nuclear reactor core. The neutrons calculation consists in numerically solving the time dependent diffusion approximation equation, which is a simplified transport equation. The numerical resolution is done with finite elements method based on a tetrahedral meshing of the computational domain, representing the reactor core, and time discretization is achieved using a θ-scheme. The transient model presents moving control rods during the time of the reaction. Therefore, cross-sections (piecewise constants) are taken into account by interpolations with respect to the velocity of the control rods. The parallelism across the time is achieved by an adequate use of the parareal in time algorithm to the handled problem. This parallel method is a predictor corrector scheme that iteratively combines the use of two kinds of numerical propagators, one coarse and one fine. Our method is made efficient by means of a coarse solver defined with large time step and fixed position control rods model, while the fine propagator is assumed to be a high order numerical approximation of the full model. The parallel implementation of our method provides a good scalability of the algorithm. Numerical results show the efficiency of the parareal method on large light water reactor transient model corresponding to the Langenbuch-Maurer-Werner benchmark.
Numerical Simulation of 3-D Supersonic Viscous Flow in an Experimental MHD Channel
NASA Technical Reports Server (NTRS)
Kato, Hiromasa; Tannehill, John C.; Gupta, Sumeet; Mehta, Unmeel B.
2004-01-01
The 3-D supersonic viscous flow in an experimental MHD channel has been numerically simulated. The experimental MHD channel is currently in operation at NASA Ames Research Center. The channel contains a nozzle section, a center section, and an accelerator section where magnetic and electric fields can be imposed on the flow. In recent tests, velocity increases of up to 40% have been achieved in the accelerator section. The flow in the channel is numerically computed using a new 3-D parabolized Navier-Stokes (PNS) algorithm that has been developed to efficiently compute MHD flows in the low magnetic Reynolds number regime. The MHD effects are modeled by introducing source terms into the PNS equations which can then be solved in a very e5uent manner. To account for upstream (elliptic) effects, the flowfield can be computed using multiple streamwise sweeps with an iterated PNS algorithm. The new algorithm has been used to compute two test cases that match the experimental conditions. In both cases, magnetic and electric fields are applied to the flow. The computed results are in good agreement with the available experimental data.
3-D-numerical approach to simulate an avalanche impact into a reservoir
NASA Astrophysics Data System (ADS)
Gabl, R.; Seibl, J.; Gems, B.; Aufleger, M.
2015-06-01
The impact of an avalanche into a reservoir induces an impulse wave, which poses a threat to population and infrastructure. For a good approximation of the generated wave height and length as well as the resulting outflow volume over structures and dams, formulas, which base on different simplifying assumptions, can be used. Further project-specific investigations by means of a scale model test or numerical simulations are advisable for complex reservoirs as well as the inclusion of hydraulic structures such as spillways. The paper presents a new approach for a 3-D-numerical simulation of an avalanche impact into a reservoir. In this model concept the energy and mass of the avalanche are represented by accelerated water on the real hill slope. Instead of snow, only water and air are used to simulate the moving avalanche with the software FLOW-3D. A significant advantage of this assumption is the self-adaptation of the model avalanche onto the terrain. In order to reach good comparability of the results with existing research at the ETH Zürich, a simplified reservoir geometry is investigated. Thus, a reference case has been analysed including a variation of three geometry parameters (still water depth in the reservoir, freeboard of the dam and reservoir width).
NASA Astrophysics Data System (ADS)
Beretta, S.; Moia, F.; Guandalini, R.; Cappelletti, F.
2012-04-01
The research activities carried out by the Environment and Sustainable Development Department of RSE S.p.A. aim to evaluate the feasibility of CO2 geological sequestration in Italy, with particular reference to the storage into saline aquifers. The identification and geological characterization of the Italian potential storage sites, together with the study of the temporal and spatial evolution of the CO2 plume within the caprock-reservoir system, are performed using different modelling tools available in the Integrated Analysis Modelling System (SIAM) entirely powered in RSE. The numerical modelling approach is the only one that allows to investigate the behaviour of the injected CO2 regarding the fluid dynamic, geochemical and geomechanical aspects and effects due to its spread, in order to verify the safety of the process. The SIAM tools allow: - Selection of potential Italian storage sites through geological and geophysical data collected in the GIS-CO2 web database; - Characterization of caprock and aquifer parameters, seismic risk and environmental link for the selected site; - Creation of the 3D simulation model for the selected domain, using the modeller METHODRdS powered by RSE and the mesh generator GMSH; - Simulation of the injection and the displacement of CO2: multiphase fluid 3D dynamics is based on the modified version of TOUGH2 model; - Evaluation of geochemical reaction effects; - Evaluation of geomechanic effects, using the coupled 3D CANT-SD finite elements code; - Detailed local analysis through the use of open source auxiliary tools, such as SHEMAT and FEHM. - 3D graphic analysis of the results. These numerical tools have been successfully used for simulating the injection and the spread of CO2 into several real Italian reservoirs and have allowed to achieve accurate results in terms of effective storage capacity and safety analysis. The 3D geological models represent the high geological complexity of the Italian subsoil, where reservoirs are
NASA Astrophysics Data System (ADS)
Conti, Caroline; Nunes, Paulo; Ducla Soares, Luís.
2013-09-01
Holoscopic imaging, also known as integral imaging, has been recently attracting the attention of the research community, as a promising glassless 3D technology due to its ability to create a more realistic depth illusion than the current stereoscopic or multiview solutions. However, in order to gradually introduce this technology into the consumer market and to efficiently deliver 3D holoscopic content to end-users, backward compatibility with legacy displays is essential. Consequently, to enable 3D holoscopic content to be delivered and presented on legacy displays, a display scalable 3D holoscopic coding approach is required. Hence, this paper presents a display scalable architecture for 3D holoscopic video coding with a three-layer approach, where each layer represents a different level of display scalability: Layer 0 - a single 2D view; Layer 1 - 3D stereo or multiview; and Layer 2 - the full 3D holoscopic content. In this context, a prediction method is proposed, which combines inter-layer prediction, aiming to exploit the existing redundancy between the multiview and the 3D holoscopic layers, with self-similarity compensated prediction (previously proposed by the authors for non-scalable 3D holoscopic video coding), aiming to exploit the spatial redundancy inherent to the 3D holoscopic enhancement layer. Experimental results show that the proposed combined prediction can improve significantly the rate-distortion performance of scalable 3D holoscopic video coding with respect to the authors' previously proposed solutions, where only inter-layer or only self-similarity prediction is used.
Li, Shengtai; Li, Hui
2012-06-14
the position of the planet, we adopt the corotating frame that allows the planet moving only in radial direction if only one planet is present. This code has been extensively tested on a number of problems. For the earthmass planet with constant aspect ratio h = 0.05, the torque calculated using our code matches quite well with the the 3D linear theory results by Tanaka et al. (2002). The code is fully parallelized via message-passing interface (MPI) and has very high parallel efficiency. Several numerical examples for both fixed planet and moving planet are provided to demonstrate the efficacy of the numerical method and code.
Simulation of a Synthetic Jet in Quiescent Air Using TLNS3D Flow Code
NASA Technical Reports Server (NTRS)
Vatsa, Veer N.; Turkel, Eli
2007-01-01
Although the actuator geometry is highly three-dimensional, the outer flowfield is nominally two-dimensional because of the high aspect ratio of the rectangular slot. For the present study, this configuration is modeled as a two-dimensional problem. A multi-block structured grid available at the CFDVAL2004 website is used as a baseline grid. The periodic motion of the diaphragm is simulated by specifying a sinusoidal velocity at the diaphragm surface with a frequency of 450 Hz, corresponding to the experimental setup. The amplitude is chosen so that the maximum Mach number at the jet exit is approximately 0.1, to replicate the experimental conditions. At the solid walls zero slip, zero injection, adiabatic temperature and zero pressure gradient conditions are imposed. In the external region, symmetry conditions are imposed on the side (vertical) boundaries and far-field conditions are imposed on the top boundary. A nominal free-stream Mach number of 0.001 is imposed in the free stream to simulate incompressible flow conditions in the TLNS3D code, which solves compressible flow equations. The code was run in unsteady (URANS) mode until the periodicity was established. The time-mean quantities were obtained by running the code for at least another 15 periods and averaging the flow quantities over these periods. The phase-locked average of flow quantities were assumed to be coincident with their values during the last full time period.
A 3D Parallel Beam Dynamics Code for Modeling High Brightness Beams in Photoinjectors
Qiang, Ji; Lidia, S.; Ryne, R.D.; Limborg, C.; /SLAC
2006-02-13
In this paper we report on IMPACT-T, a 3D beam dynamics code for modeling high brightness beams in photoinjectors and rf linacs. IMPACT-T is one of the few codes used in the photoinjector community that has a parallel implementation, making it very useful for high statistics simulations of beam halos and beam diagnostics. It has a comprehensive set of beamline elements, and furthermore allows arbitrary overlap of their fields. It is unique in its use of space-charge solvers based on an integrated Green function to efficiently and accurately treat beams with large aspect ratio, and a shifted Green function to efficiently treat image charge effects of a cathode. It is also unique in its inclusion of energy binning in the space-charge calculation to model beams with large energy spread. Together, all these features make IMPACT-T a powerful and versatile tool for modeling beams in photoinjectors and other systems. In this paper we describe the code features and present results of IMPACT-T simulations of the LCLS photoinjectors. We also include a comparison of IMPACT-T and PARMELA results.
A 3d Parallel Beam Dynamics Code for Modeling High BrightnessBeams in Photoinjectors
Qiang, J.; Lidia, S.; Ryne, R.; Limborg, C.
2005-05-16
In this paper we report on IMPACT-T, a 3D beam dynamics code for modeling high brightness beams in photoinjectors and rf linacs. IMPACT-T is one of the few codes used in the photoinjector community that has a parallel implementation, making it very useful for high statistics simulations of beam halos and beam diagnostics. It has a comprehensive set of beamline elements, and furthermore allows arbitrary overlap of their fields. It is unique in its use of space-charge solvers based on an integrated Green function to efficiently and accurately treat beams with large aspect ratio, and a shifted Green function to efficiently treat image charge effects of a cathode. It is also unique in its inclusion of energy binning in the space-charge calculation to model beams with large energy spread. Together, all these features make IMPACT-T a powerful and versatile tool for modeling beams in photoinjectors and other systems. In this paper we describe the code features and present results of IMPACT-T simulations of the LCLS photoinjectors. We also include a comparison of IMPACT-T and PARMELA results.
Implementation of a 3d numerical model of a folded multilayer carbonate aquifer
NASA Astrophysics Data System (ADS)
Di Salvo, Cristina; Guyennon, Nicolas; Romano, Emanuele; Bruna Petrangeli, Anna; Preziosi, Elisabetta
2016-04-01
The main objective of this research is to present a case study of the numerical model implementation of a complex carbonate, structurally folded aquifer, with a finite difference, porous equivalent model. The case study aquifer (which extends over 235 km2 in the Apennine chain, Central Italy) provides a long term average of 3.5 m3/s of good quality groundwater to the surface river network, sustaining the minimum vital flow, and it is planned to be exploited in the next years for public water supply. In the downstream part of the river in the study area, a "Site of Community Importance" include the Nera River for its valuable aquatic fauna. However, the possible negative effects of the foreseen exploitation on groundwater dependent ecosystems are a great concern and model grounded scenarios are needed. This multilayer aquifer was conceptualized as five hydrostratigraphic units: three main aquifers (the uppermost unconfined, the central and the deepest partly confined), are separated by two locally discontinuous aquitards. The Nera river cuts through the two upper aquifers and acts as the main natural sink for groundwater. An equivalent porous medium approach was chosen. The complex tectonic structure of the aquifer requires several steps in defining the conceptual model; the presence of strongly dipping layers with very heterogeneous hydraulic conductivity, results in different thicknesses of saturated portions. Aquifers can have both unconfined or confined zones; drying and rewetting must be allowed when considering recharge/discharge cycles. All these characteristics can be included in the conceptual and numerical model; however, being the number of flow and head target scarce, the over-parametrization of the model must be avoided. Following the principle of parsimony, three steady state numerical models were developed, starting from a simple model, and then adding complexity: 2D (single layer), QUASI -3D (with leackage term simulating flow through aquitards) and
Liner Optimization Studies Using the Ducted Fan Noise Prediction Code TBIEM3D
NASA Technical Reports Server (NTRS)
Dunn, M. H.; Farassat, F.
1998-01-01
In this paper we demonstrate the usefulness of the ducted fan noise prediction code TBIEM3D as a liner optimization design tool. Boundary conditions on the interior duct wall allow for hard walls or a locally reacting liner with axially segmented, circumferentially uniform impedance. Two liner optimization studies are considered in which farfield noise attenuation due to the presence of a liner is maximized by adjusting the liner impedance. In the first example, the dependence of optimal liner impedance on frequency and liner length is examined. Results show that both the optimal impedance and attenuation levels are significantly influenced by liner length and frequency. In the second example, TBIEM3D is used to compare radiated sound pressure levels between optimal and non-optimal liner cases at conditions designed to simulate take-off. It is shown that significant noise reduction is achieved for most of the sound field by selecting the optimal or near optimal liner impedance. Our results also indicate that there is relatively large region of the impedance plane over which optimal or near optimal liner behavior is attainable. This is an important conclusion for the designer since there are variations in liner characteristics due to manufacturing imprecisions.
Gray coded trapezoidal fringes for 3-D surface-shape measurement
NASA Astrophysics Data System (ADS)
Pérez, Oscar G.; Flores, Jorge L.; García-Torales, G.; Muñoz-G, J. A.; Soto, Horacio; Balderas, Sandra E.
2014-09-01
We propose a two-step trapezoidal-pattern phase-shifting method for 3-D surface-shape measurements. Shape measurements by trapezoidal phase-shifting methods require high-quality trapezoidal patterns. Furthermore, most of the video projectors are nonlinear, making it difficult to generate high quality phase without nonlinearity calibration and correction. To overcome the limitations, we propose a method for synthesizing trapezoidal intensity fringes as a way to solve the problems caused by projector/camera gamma nonlinearity. The fringe generation technique consists of projecting and acquiring a temporal sequence of strictly binary color patterns (Gray code), whose (adequately weighted) average leads to trapezoidal fringe patterns with the required number of bits, which allows a reliable three-dimensional profile reconstruction using phase-shifting methods. Validation experiments are presented.
Numerical 3D Hydrodynamics Study of Gravitational Instabilities in a Circumbinary Disk
NASA Astrophysics Data System (ADS)
Desai, Karna Mahadev; Steiman-Cameron, Thomas Y.; Michael, Scott; Cai, Kai; Durisen, Richard H.
2016-01-01
We present a 3D hydrodynamical study of gravitational instabilities (GIs) in a circumbinary protoplanetary disk around a Solar mass star and a brown dwarf companion (0.02 M⊙). GIs can play an important, and at times dominant, role in driving the structural evolution of protoplanetary disks. The reported simulations were performed employing CHYMERA, a radiative 3D hydrodynamics code developed by the Indiana University Hydrodynamics Group. The simulations include disk self-gravity and radiative cooling governed by realistic dust opacities. We examine the role of GIs in modulating the thermodynamic state of the disks, and determine the strengths of GI-induced density waves, non-axisymmetric density structures, radial mass transport, and gravitational torques. The principal goal of this study is to determine how the presence of the companion affects the nature and strength of GIs. Results are compared with a parallel simulation of a protoplanetary disk without the presence of the brown dwarf binary companion. We detect no fragmentation in either disk. A persistent vortex forms in the inner region of both disks. The vortex seems to be stabilized by the presence of the binary companion.
A 3D-PNS computer code for the calculation of supersonic combusting flows
NASA Technical Reports Server (NTRS)
Chitsomboon, Tawit; Northam, G. Burton
1988-01-01
A computer code has been developed based on the three-dimensional parabolized Navier-Stokes (PNS) equations which govern the supersonic combusting flow of the hydrogen-air system. The finite difference algorithm employed was a hybrid of the Schiff-Steger algorithm and the Vigneron, et al., algorithm which is fully implicit and fully coupled. The combustion of hydrogen and air was modeled by the finite-rate two-step combustion model of Rogers-Chinitz. A new dependent variable vector was introduced to simplify the numerical algorithm. Robustness of the algorithm was considerably enhanced by introducing an adjustable parameter. The computer code was used to solve a premixed shock-induced combustion problem and the results were compared with those of a full Navier-Stokes code. Reasonably good agreement was obtained at a fraction of the cost of the full Navier-Stokes procedure.
PORTA: A Massively Parallel Code for 3D Non-LTE Polarized Radiative Transfer
NASA Astrophysics Data System (ADS)
Štěpán, J.
2014-10-01
The interpretation of the Stokes profiles of the solar (stellar) spectral line radiation requires solving a non-LTE radiative transfer problem that can be very complex, especially when the main interest lies in modeling the linear polarization signals produced by scattering processes and their modification by the Hanle effect. One of the main difficulties is due to the fact that the plasma of a stellar atmosphere can be highly inhomogeneous and dynamic, which implies the need to solve the non-equilibrium problem of generation and transfer of polarized radiation in realistic three-dimensional stellar atmospheric models. Here we present PORTA, a computer program we have developed for solving, in three-dimensional (3D) models of stellar atmospheres, the problem of the generation and transfer of spectral line polarization taking into account anisotropic radiation pumping and the Hanle and Zeeman effects in multilevel atoms. The numerical method of solution is based on a highly convergent iterative algorithm, whose convergence rate is insensitive to the grid size, and on an accurate short-characteristics formal solver of the Stokes-vector transfer equation which uses monotonic Bezier interpolation. In addition to the iterative method and the 3D formal solver, another important feature of PORTA is a novel parallelization strategy suitable for taking advantage of massively parallel computers. Linear scaling of the solution with the number of processors allows to reduce the solution time by several orders of magnitude. We present useful benchmarks and a few illustrations of applications using a 3D model of the solar chromosphere resulting from MHD simulations. Finally, we present our conclusions with a view to future research. For more details see Štěpán & Trujillo Bueno (2013).
NASA Astrophysics Data System (ADS)
Bi, J.; Zhou, X. P.; Qian, Q. H.
2016-05-01
General particle dynamics (GPD), which is a novel meshless numerical method, is proposed to simulate the initiation, propagation and coalescence of 3D pre-existing penetrating and embedded flaws under biaxial compression. The failure process for rock-like materials subjected to biaxial compressive loads is investigated using the numerical code GPD3D. Moreover, internal crack evolution processes are successfully simulated using GPD3D. With increasing lateral stress, the secondary cracks keep growing in the samples, while the growth of the wing cracks is restrained. The samples are mainly split into fragments in a shear failure mode under biaxial compression, which is different from the splitting failure of the samples subjected to uniaxial compression. For specimens with macroscopic pre-existing flaws, the simulated types of cracks, the simulated coalescence types and the simulated failure modes are in good agreement with the experimental results.
NASA Astrophysics Data System (ADS)
Gillespie, K. M.; Speirs, D. C.; Ronald, K.; McConville, S. L.; Phelps, A. D. R.; Bingham, R.; Cross, A. W.; Robertson, C. W.; Whyte, C. G.; He, W.; Vorgul, I.; Cairns, R. A.; Kellett, B. J.
2008-12-01
Auroral Kilometric Radiation (AKR), occurs naturally in the polar regions of the Earth's magnetosphere where electrons are accelerated by electric fields into the increasing planetary magnetic dipole. Here conservation of the magnetic moment converts axial to rotational momentum forming a horseshoe distribution in velocity phase space. This distribution is unstable to cyclotron emission with radiation emitted in the X-mode. In a scaled laboratory reproduction of this process, a 75-85 keV electron beam of 5-40 A was magnetically compressed by a system of solenoids and emissions were observed for cyclotron frequencies of 4.42 GHz and 11.7 GHz resonating with near cut-off TE0,1 and TE0,3 modes, respectively. Here we compare these measurements with numerical predictions from the 3D PiC code KARAT. The 3D simulations accurately predicted the radiation modes and frequencies produced by the experiment. The predicted conversion efficiency between electron kinetic and wave field energy of around 1% is close to the experimental measurements and broadly consistent with quasi-linear theoretical analysis and geophysical observations.
3D numerical simulation of the evolutionary process of aeolian downsized crescent-shaped dunes
NASA Astrophysics Data System (ADS)
Zhou, Xiaosi; Zhang, Yang; Wang, Yuan; Li, Min
2016-06-01
A dune constitutive model was coupled with a large eddy simulation (LES) with the Smagorinsky subgrid-scale (SGS) model to accurately describe the evolutionary process of dunes from the macroscopic perspective of morphological dynamics. A 3D numerical simulation of the evolution of aeolian downsized crescent-shaped dunes was then performed. The evolution of the 3D structure of Gaussian-shaped dunes was simulated under the influence of gravity modulation, which was the same with the vertical oscillation of the sand bed to adjust the threshold of sand grain liftoff in wind tunnel experiments under the same wind speed. The influence of gravity modulation intensity on the characteristic scale parameter of the dune was discussed. Results indicated that the crescent shape of the dune was reproduced with the action of gravity during regulation of the saturation of wind-sand flow at specific times. The crescent shape was not dynamically maintained as time passed, and the dunes dwindled until they reached final decomposition because of wind erosion. The height of the dunes decreased over time, and the height-time curve converged as the intensity of modulation increased linearly. The results qualitatively agreed with those obtained from wind tunnel experiments.
3D numerical simulation of the transport of chemical signature compounds from buried landmines
NASA Astrophysics Data System (ADS)
Irrazabal, Maik; Borrero, Ernesto; Briano, Julio G.; Castro, Miguel; Hernandez, Samuel P.
2005-06-01
The transport of the chemical signature compounds from buried landmines in a three-dimensional (3D) array has been numerically modeled using the finite-volume technique. Compounds such as trinitrotoluene, dinitrotoluene, and their degradation products, are semi volatile and somewhat soluble in water. Furthermore, they can strongly adsorb to the soil and undergo chemical and biological degradation. Consequently, the spatial and temporal concentration distributions of such chemicals depend on the mobility of the water and gaseous phases, their molecular and mechanical diffusion, adsorption characteristics, soil water content, compaction, and environmental factors. A 3D framework is required since two-dimensional (2D) symmetry may easily fade due to terrain topography: non-flat surfaces, soil heterogeneity, or underground fractures. The spatial and temporal distribution of the chemical-signature-compounds, in an inclined grid has been obtained. The fact that the chemicals may migrate horizontally, giving higher surface concentrations at positions not directly on top of the objects, emphasizes the need for understanding the transport mechanism when a chemical detector is used. Deformation in the concentration contours after rainfall is observed in the inclined surface and is attributed to both: the advective flux, and to the water flux at the surface caused by the slope. The analysis of the displacements in the position of the maximum concentrations at the surface, respect to the actual location of the mine, in an inclined system, is presented.
Ruh, Dominic; Tränkle, Benjamin; Rohrbach, Alexander
2011-10-24
Multi-dimensional, correlated particle tracking is a key technology to reveal dynamic processes in living and synthetic soft matter systems. In this paper we present a new method for tracking micron-sized beads in parallel and in all three dimensions - faster and more precise than existing techniques. Using an acousto-optic deflector and two quadrant-photo-diodes, we can track numerous optically trapped beads at up to tens of kHz with a precision of a few nanometers by back-focal plane interferometry. By time-multiplexing the laser focus, we can calibrate individually all traps and all tracking signals in a few seconds and in 3D. We show 3D histograms and calibration constants for nine beads in a quadratic arrangement, although trapping and tracking is easily possible for more beads also in arbitrary 2D arrangements. As an application, we investigate the hydrodynamic coupling and diffusion anomalies of spheres trapped in a 3 × 3 arrangement. PMID:22109012
Energy Science and Technology Software Center (ESTSC)
2013-06-24
Version 07 TART2012 is a coupled neutron-photon Monte Carlo transport code designed to use three-dimensional (3-D) combinatorial geometry. Neutron and/or photon sources as well as neutron induced photon production can be tracked. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART2012 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared tomore » other similar codes. Use of the entire system can save you a great deal of time and energy. TART2012 extends the general utility of the code to even more areas of application than available in previous releases by concentrating on improving the physics, particularly with regard to improved treatment of neutron fission, resonance self-shielding, molecular binding, and extending input options used by the code. Several utilities are included for creating input files and displaying TART results and data. TART2012 uses the latest ENDF/B-VI, Release 8, data. New for TART2012 is the use of continuous energy neutron cross sections, in addition to its traditional multigroup cross sections. For neutron interaction, the data are derived using ENDF-ENDL2005 and include both continuous energy cross sections and 700 group neutron data derived using a combination of ENDF/B-VI, Release 8, and ENDL data. The 700 group structure extends from 10-5 eV up to 1 GeV. Presently nuclear data are only available up to 20 MeV, so that only 616 of the groups are currently used. For photon interaction, 701 point photon data were derived using the Livermore EPDL97 file. The new 701 point structure extends from 100 eV up to 1 GeV, and is currently used over this entire energy range. TART2012 completely supersedes all older versions of TART, and it is strongly recommended that one use only the most recent version of TART2012 and its data files. Check authors homepage for related information: http
Validation of Heat Transfer and Film Cooling Capabilities of the 3-D RANS Code TURBO
NASA Technical Reports Server (NTRS)
Shyam, Vikram; Ameri, Ali; Chen, Jen-Ping
2010-01-01
The capabilities of the 3-D unsteady RANS code TURBO have been extended to include heat transfer and film cooling applications. The results of simulations performed with the modified code are compared to experiment and to theory, where applicable. Wilcox s k-turbulence model has been implemented to close the RANS equations. Two simulations are conducted: (1) flow over a flat plate and (2) flow over an adiabatic flat plate cooled by one hole inclined at 35 to the free stream. For (1) agreement with theory is found to be excellent for heat transfer, represented by local Nusselt number, and quite good for momentum, as represented by the local skin friction coefficient. This report compares the local skin friction coefficients and Nusselt numbers on a flat plate obtained using Wilcox's k-model with the theory of Blasius. The study looks at laminar and turbulent flows over an adiabatic flat plate and over an isothermal flat plate for two different wall temperatures. It is shown that TURBO is able to accurately predict heat transfer on a flat plate. For (2) TURBO shows good qualitative agreement with film cooling experiments performed on a flat plate with one cooling hole. Quantitatively, film effectiveness is under predicted downstream of the hole.
Studies of coupled cavity LINAC (CCL) accelerating structures with 3-D codes
Spalek, G.
2000-08-01
The cw CCL being designed for the Accelerator Production of Tritium (APT) project accelerates protons from 96 MeV to 211 MeV. It consists of 99 segments each containing up to seven accelerating cavities. Segments are coupled by intersegment coupling cavities and grouped into supermodules. The design method needs to address not only basic cavity sizing for a given coupling and pi/2 mode frequency, but also the effects of high power densities on the cavity frequency, mechanical stresses, and the structure's stop band during operation. On the APT project, 3-D RF (Ansoft Corp.'s HFSS) and coupled RF/structural (Ansys Inc.'s ANSYS) codes are being used. to develop tools to address the above issues and guide cooling channel design. The code's predictions are being checked against available low power Aluminum models. Stop band behavior under power will be checked once the tools are extended to CCDTL structures that have been tested at high power. A summary of calculations made to date and agreement with measured results will be presented.
ORPHEE research reactor: 3D core depletion calculation using Monte-Carlo code TRIPOLI-4®
NASA Astrophysics Data System (ADS)
Damian, F.; Brun, E.
2014-06-01
ORPHEE is a research reactor located at CEA Saclay. It aims at producing neutron beams for experiments. This is a pool-type reactor (heavy water), and the core is cooled by light water. Its thermal power is 14 MW. ORPHEE core is 90 cm height and has a cross section of 27x27 cm2. It is loaded with eight fuel assemblies characterized by a various number of fuel plates. The fuel plate is composed of aluminium and High Enriched Uranium (HEU). It is a once through core with a fuel cycle length of approximately 100 Equivalent Full Power Days (EFPD) and with a maximum burnup of 40%. Various analyses under progress at CEA concern the determination of the core neutronic parameters during irradiation. Taking into consideration the geometrical complexity of the core and the quasi absence of thermal feedback for nominal operation, the 3D core depletion calculations are performed using the Monte-Carlo code TRIPOLI-4® [1,2,3]. A preliminary validation of the depletion calculation was performed on a 2D core configuration by comparison with the deterministic transport code APOLLO2 [4]. The analysis showed the reliability of TRIPOLI-4® to calculate a complex core configuration using a large number of depleting regions with a high level of confidence.
Slab detachment in laterally varying subduction zones: 3-D numerical modeling
NASA Astrophysics Data System (ADS)
Duretz, T.; Gerya, T. V.; Spakman, W.
2014-03-01
Understanding the three-dimensional (3-D) dynamics of subduction-collision systems is a longstanding challenge in geodynamics. We investigate the impact of slab detachment in collision systems that are subjected to along-trench variations. High-resolution thermomechanical numerical models, encompassing experimentally derived flow laws and a pseudo free surface, are employed to unravel lithospheric and topographic evolutions. First, we consider coeval subduction of adjacent continental and oceanic lithospheres (SCO). This configuration yields to two-stage slab detachment during collision, topographic buildup and extrusion, variable along-trench convergence rates, and associated trench deformation. The second setting considers a convergent margin, which is laterally limited by a transform boundary (STB). Such collisional system is affected by a single slab detachment, little trench deformation, and moderately confined upper plate topography. The effect of initial thermal slab age on SCO and STB models are explored. Similarities with natural analogs along the Arabia-Eurasia collision are discussed.
Temperature distributions in the laser-heated diamond anvil cell from 3-D numerical modeling
Rainey, E. S. G.; Kavner, A.; Hernlund, J. W.
2013-11-28
We present TempDAC, a 3-D numerical model for calculating the steady-state temperature distribution for continuous wave laser-heated experiments in the diamond anvil cell. TempDAC solves the steady heat conduction equation in three dimensions over the sample chamber, gasket, and diamond anvils and includes material-, temperature-, and direction-dependent thermal conductivity, while allowing for flexible sample geometries, laser beam intensity profile, and laser absorption properties. The model has been validated against an axisymmetric analytic solution for the temperature distribution within a laser-heated sample. Example calculations illustrate the importance of considering heat flow in three dimensions for the laser-heated diamond anvil cell. In particular, we show that a “flat top” input laser beam profile does not lead to a more uniform temperature distribution or flatter temperature gradients than a wide Gaussian laser beam.
Numerical 3D study of FZ growth: dependence on growth parameters and melt instability
NASA Astrophysics Data System (ADS)
Ratnieks, G.; Muižnieks, A.; Mühlbauer, A.; Raming, G.
2001-08-01
Three-dimensional modelling of the floating zone (needle-eye) crystal growth process is carried out to analyse numerically the stability of the melt flow and the influence of the crystal rotation rate and inductor slit width on the 3D flow field and on the grown crystal resistivity. The unsteadiness of the melt is simulated and it is found that for the considered growth parameters a steady-state flow can be a reasonable approximation to the unsteady melt motion. The parametric studies have shown that increasing the rotation rate essentially changes the flow pattern and weakens the rotational striations, while the inductor slit width has a more local influence on these characteristics.
FERM3D: A finite element R-matrix electron molecule scattering code
NASA Astrophysics Data System (ADS)
Tonzani, Stefano
2007-01-01
FERM3D is a three-dimensional finite element program, for the elastic scattering of a low energy electron from a general polyatomic molecule, which is converted to a potential scattering problem. The code is based on tricubic polynomials in spherical coordinates. The electron-molecule interaction is treated as a sum of three terms: electrostatic, exchange, and polarization. The electrostatic term can be extracted directly from ab initio codes ( GAUSSIAN 98 in the work described here), while the exchange term is approximated using a local density functional. A local polarization potential based on density functional theory [C. Lee, W. Yang, R.G. Parr, Phys. Rev. B 37 (1988) 785] describes the long range attraction to the molecular target induced by the scattering electron. Photoionization calculations are also possible and illustrated in the present work. The generality and simplicity of the approach is important in extending electron-scattering calculations to more complex targets than it is possible with other methods. Program summaryTitle of program:FERM3D Catalogue identifier:ADYL_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYL_v1_0 Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Computer for which the program is designed and others on which it has been tested:Intel Xeon, AMD Opteron 64 bit, Compaq Alpha Operating systems or monitors under which the program has been tested:HP Tru64 Unix v5.1, Red Hat Linux Enterprise 3 Programming language used:Fortran 90 Memory required to execute with typical data:900 MB (neutral CO 2), 2.3 GB (ionic CO 2), 1.4 GB (benzene) No. of bits in a word:32 No. of processors used:1 Has the code been vectorized?:No No. of lines in distributed program, including test data, etc.:58 383 No. of bytes in distributed program, including test data, etc.:561 653 Distribution format:tar.gzip file CPC Program library subprograms used:ADDA, ACDP Nature of physical problem:Scattering of an
A study of the earth radiation budget using a 3D Monte-Carlo radiative transer code
NASA Astrophysics Data System (ADS)
Okata, M.; Nakajima, T.; Sato, Y.; Inoue, T.; Donovan, D. P.
2013-12-01
The purpose of this study is to evaluate the earth's radiation budget when data are available from satellite-borne active sensors, i.e. cloud profiling radar (CPR) and lidar, and a multi-spectral imager (MSI) in the project of the Earth Explorer/EarthCARE mission. For this purpose, we first developed forward and backward 3D Monte Carlo radiative transfer codes that can treat a broadband solar flux calculation including thermal infrared emission calculation by k-distribution parameters of Sekiguchi and Nakajima (2008). In order to construct the 3D cloud field, we tried the following three methods: 1) stochastic cloud generated by randomized optical thickness each layer distribution and regularly-distributed tilted clouds, 2) numerical simulations by a non-hydrostatic model with bin cloud microphysics model and 3) Minimum cloud Information Deviation Profiling Method (MIDPM) as explained later. As for the method-2 (numerical modeling method), we employed numerical simulation results of Californian summer stratus clouds simulated by a non-hydrostatic atmospheric model with a bin-type cloud microphysics model based on the JMA NHM model (Iguchi et al., 2008; Sato et al., 2009, 2012) with horizontal (vertical) grid spacing of 100m (20m) and 300m (20m) in a domain of 30km (x), 30km (y), 1.5km (z) and with a horizontally periodic lateral boundary condition. Two different cell systems were simulated depending on the cloud condensation nuclei (CCN) concentration. In the case of horizontal resolution of 100m, regionally averaged cloud optical thickness,
Tidal dynamics of the Terminos Lagoon, Mexico: observations and 3D numerical modelling
NASA Astrophysics Data System (ADS)
Contreras Ruiz Esparza, Adolfo; Douillet, Pascal; Zavala-Hidalgo, Jorge
2014-09-01
The tidal circulation patterns in the Terminos Lagoon were studied based on the analysis of 1 year of measurements and numerical simulations using a baroclinic 3D hydrodynamic model, the MARS3D. A gauging network was installed consisting of six self-recording pressure-temperature sensors, a tide gauge station and two current profilers, with pressure and temperature sensors moored in the main lagoon inlets. Model simulations were validated against current and sea level observations and were used to analyse the circulation patterns caused by the tidal forcing. The numerical model was forced with eight harmonic components, four diurnal ( K 1, O 1, P 1, Q 1) and four semi-diurnal ( M 2, S 2, N 2, K 2), extracted from the TPX0.7 database. The tidal patterns in the study area vary from mixed, mainly diurnal in the two main inlets of the lagoon, to diurnal in its interior. The tidal residual circulation inside the lagoon is dominated by a cyclonic gyre. The results indicate a net flux from the southwest Ciudad del Carmen inlet (CdC) towards the northeast Puerto Real inlet (PtR) along the southern side of the lagoon and the opposite in the northern side. The results indicate two areas of strong currents in the vicinity of the inlets and weak currents inside the lagoon. The area of strong currents in the vicinity of the CdC inlet is larger than that observed in the PtR inlet. Nevertheless, the current analysis indicates that the highest current speeds, which can reach a magnitude of 1.9 m s-1, occurred in PtR. A further analysis of the tide distortion in the inlets revealed that both passages are ebb dominated.
NASA Astrophysics Data System (ADS)
Zhang, Yujia; Yilmaz, Alper
2016-06-01
Surface reconstruction using coded structured light is considered one of the most reliable techniques for high-quality 3D scanning. With a calibrated projector-camera stereo system, a light pattern is projected onto the scene and imaged by the camera. Correspondences between projected and recovered patterns are computed in the decoding process, which is used to generate 3D point cloud of the surface. However, the indirect illumination effects on the surface, such as subsurface scattering and interreflections, will raise the difficulties in reconstruction. In this paper, we apply maximum min-SW gray code to reduce the indirect illumination effects of the specular surface. We also analysis the errors when comparing the maximum min-SW gray code and the conventional gray code, which justifies that the maximum min-SW gray code has significant superiority to reduce the indirect illumination effects. To achieve sub-pixel accuracy, we project high frequency sinusoidal patterns onto the scene simultaneously. But for specular surface, the high frequency patterns are susceptible to decoding errors. Incorrect decoding of high frequency patterns will result in a loss of depth resolution. Our method to resolve this problem is combining the low frequency maximum min-SW gray code and the high frequency phase shifting code, which achieves dense 3D reconstruction for specular surface. Our contributions include: (i) A complete setup of the structured light based 3D scanning system; (ii) A novel combination technique of the maximum min-SW gray code and phase shifting code. First, phase shifting decoding with sub-pixel accuracy. Then, the maximum min-SW gray code is used to resolve the ambiguity resolution. According to the experimental results and data analysis, our structured light based 3D scanning system enables high quality dense reconstruction of scenes with a small number of images. Qualitative and quantitative comparisons are performed to extract the advantages of our new
3D-radiation hydro simulations of disk-planet interactions. I. Numerical algorithm and test cases
NASA Astrophysics Data System (ADS)
Klahr, H.; Kley, W.
2006-01-01
We study the evolution of an embedded protoplanet in a circumstellar disk using the 3D-Radiation Hydro code TRAMP, and treat the thermodynamics of the gas properly in three dimensions. The primary interest of this work lies in the demonstration and testing of the numerical method. We show how far numerical parameters can influence the simulations of gap opening. We study a standard reference model under various numerical approximations. Then we compare the commonly used locally isothermal approximation to the radiation hydro simulation using an equation for the internal energy. Models with different treatments of the mass accretion process are compared. Often mass accumulates in the Roche lobe of the planet creating a hydrostatic atmosphere around the planet. The gravitational torques induced by the spiral pattern of the disk onto the planet are not strongly affected in the average magnitude, but the short time scale fluctuations are stronger in the radiation hydro models. An interesting result of this work lies in the analysis of the temperature structure around the planet. The most striking effect of treating the thermodynamics properly is the formation of a hot pressure-supported bubble around the planet with a pressure scale height of H/R ≈ 0.5 rather than a thin Keplerian circumplanetary accretion disk.
PFLOW: A 3-D Numerical Modeling Tool for Calculating Fluid-Pressure Diffusion from Coulomb Strain
NASA Astrophysics Data System (ADS)
Wolf, L. W.; Lee, M.; Meir, A.; Dyer, G.; Ma, K.; Chan, C.
2009-12-01
A new 3D time-dependent pore-pressure diffusion model PFLOW is developed to investigate the response of pore fluids to the crustal deformation generated by strong earthquakes in heterogeneous geologic media. Given crustal strain generated by changes in Coulomb stress, this MATLAB-based code uses Skempton's coefficient to calculate resulting changes fluid pressure. Pore-pressure diffusion can be tracked over time in a user-defined model space with user-prescribed Neumann or Dirchilet boundary conditions and with spatially variable values of permeability. PFLOW employs linear or quadratic finite elements for spatial discretization and first order or second order, explicit or implicit finite difference discretization in time. PFLOW is easily interfaced with output from deformation modeling programs such as Coulomb (Toda et al., 2007) or 3D-DEF (Gomberg and Ellis, 1994). The code is useful for investigating to first-order the evolution of pore pressure changes induced by changes in Coulomb stress and their possible relation to water-level changes in wells or changes in stream discharge. It can also be used for student research and classroom instruction. As an example application, we calculate the coseismic pore pressure changes and diffusion induced by volumetric strain associated with the 1999 Chi-Chi earthquake (Mw = 7.6) in Taiwan. The Chi-Chi earthquake provides an unique opportunity to investigate the spatial and time-dependent poroelastic response of near-field rocks and sediments because there exist extensive observational data of water-level changes and crustal deformation. The integrated model allows us to explore whether changes in Coulomb stress can adequately explain hydrologic anomalies observed in areas such as Taiwan’s western foothills and the Choshui River alluvial plain. To calculate coseismic strain, we use the carefully calibrated finite fault-rupture model of Ma et al. (2005) and the deformation modeling code Coulomb 3.1 (Toda et al., 2007