Numerical simulation of 3D backward facing step flows at various Reynolds numbers

Directory of Open Access Journals (Sweden)

Louda Petr

2015-01-01

Full Text Available The work deals with the numerical simulation of 3D turbulent flow over backward facing step in a narrow channel. The mathematical model is based on the RANS equations with an explicit algebraic Reynolds stress model (EARSM. The numerical method uses implicit finite volume upwind discretization. While the eddy viscosity models fail in predicting complex 3D flows, the EARSM model is shown to provide results which agree well with experimental PIV data. The reference experimental data provide the 3D flow field. The simulations are compared with experiment for 3 values of Reynolds number.

3-D numerical modelling of flow around a groin

DEFF Research Database (Denmark)

Miller, R.; Roulund, A.; Sumer, B. Mutlu

2003-01-01

A 3-D flow code, EllipSys3D, has been implemented to simulate the 3-D flow around a groin in steady current. The k  turbulence model has been used for closure. Two kinds of groins are considered: (1) A vertical-wall groin, and (2) A groin with a side slope. Steady-flow simulations were conducted...

Gravitational radiation and 3D numerical relativity

International Nuclear Information System (INIS)

Nakamura, T.

1986-01-01

Study of Numerical Relativity in Kyoto is reviewed. Main topics discussed are 2D rotating collapse, phase cancellation effects and perturbation calculation of the gravitational radiation from a particle falling into a black hole. New numerical results on 3D time evolution of pure gravitational waves are also presented

Numerical modeling of time domain 3-D problems in accelerator physics

International Nuclear Information System (INIS)

Harfoush, F.A.; Jurgens, T.G.

1990-06-01

Time domain analysis is relevant in particle accelerators to study the electromagnetic field interaction of a moving source particle on a lagging test particle as the particles pass an accelerating cavity or some other structure. These fields are called wake fields. The travelling beam inside a beam pipe may undergo more complicated interactions with its environment due to the presence of other irregularities like wires, thin slots, joints and other types of obstacles. Analytical solutions of such problems is impossible and one has to resort to a numerical method. In this paper we present results of our first attempt to model these problems in 3-D using our finite difference time domain (FDTD) code. 10 refs., 9 figs

3D numerical modeling of longwall mining with top-coal caving

Energy Technology Data Exchange (ETDEWEB)

Yasitli, N.E.; Unver, B. [Hacettepe University, Ankara (Turkey). Dept. of Mining Engineering

2005-02-01

Top-coal caving is the key factor affecting the efficiency of production at thick-coal seams. During production of top coal by caving behind the face not only a significant amount of coal is lost in the goaf but the coal drawn by means of caving is diluted considerably with surrounding rock. Therefore, it is not possible to carry out an efficient production operation unless caving of top coal behind the face is optimized. In this paper, results of 3D modeling of the top-coal-caving mechanism by using the finite difference code FLAC3{sup D} at the M3 longwall panel of the Omerler Underground Mine located at Tuncbilek (Turkey) are presented. According to the modeling results, maximum vertical abutment stresses were formed at a distance of 7m in front of the face. An analysis of the conditions of top coal has revealed that a 1.5 m thick layer of coal just above the shield supports is well fractured. However, a 3.5 m thick layer of coal above the fractured part is either not fractured or is fractured in the form of large blocks leading to obstruction of windows of shields during coal drawing. It is concluded that, in order to decrease dilution and increase extraction ratio and efficiency of operation, top coal should be as uniformly fractured as possible. Hence, an efficient and continuous coal flowing behind the face can be maintained. A special pre-fracture blasting strategy just sufficient enough to form cracks in the top coal is suggested by means of comparing with the results of numerical modeling.

2D and 3D numerical modeling of seismic waves from explosion sources

International Nuclear Information System (INIS)

McLaughlin, K.L.; Stevens, J.L.; Barker, T.G.; Shkoller, B.; Day, S.M.

1993-01-01

Over the last decade, nonlinear and linear 2D axisymmetric finite difference codes have been used in conjunction with far-field seismic Green's functions to simulate seismic waves from a variety of sources. In this paper we briefly review some of the results and conclusions that have resulted from numerical simulations and explosion modeling in support of treaty verification research at S-CUBED in the last decade. We then describe in more detail the results from two recent projects. Our goal is to provide a flavor for the kinds of problems that can be examined with numerical methods for modeling excitation of seismic waves from explosions. Two classes of problems have been addressed; nonlinear and linear near-source interactions. In both classes of problems displacements and tractions are saved on a closed surface in the linear region and the representation theorem is used to propagate the seismic waves to the far-field

Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model

Energy Technology Data Exchange (ETDEWEB)

Baudron, Anne-Marie, E-mail: anne-marie.baudron@cea.fr [Laboratoire de Recherche Conventionné MANON, CEA/DEN/DANS/DM2S and UPMC-CNRS/LJLL (France); CEA-DRN/DMT/SERMA, CEN-Saclay, 91191 Gif sur Yvette Cedex (France); Lautard, Jean-Jacques, E-mail: jean-jacques.lautard@cea.fr [Laboratoire de Recherche Conventionné MANON, CEA/DEN/DANS/DM2S and UPMC-CNRS/LJLL (France); CEA-DRN/DMT/SERMA, CEN-Saclay, 91191 Gif sur Yvette Cedex (France); Maday, Yvon, E-mail: maday@ann.jussieu.fr [Sorbonne Universités, UPMC Univ Paris 06, UMR 7598, Laboratoire Jacques-Louis Lions and Institut Universitaire de France, F-75005, Paris (France); Laboratoire de Recherche Conventionné MANON, CEA/DEN/DANS/DM2S and UPMC-CNRS/LJLL (France); Brown Univ, Division of Applied Maths, Providence, RI (United States); Riahi, Mohamed Kamel, E-mail: riahi@cmap.polytechnique.fr [Laboratoire de Recherche Conventionné MANON, CEA/DEN/DANS/DM2S and UPMC-CNRS/LJLL (France); CMAP, Inria-Saclay and X-Ecole Polytechnique, Route de Saclay, 91128 Palaiseau Cedex (France); Salomon, Julien, E-mail: salomon@ceremade.dauphine.fr [CEREMADE, Univ Paris-Dauphine, Pl. du Mal. de Lattre de Tassigny, F-75016, Paris (France)

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.

An innovative 3-D numerical modelling procedure for simulating repository-scale excavations in rock - SAFETI

Energy Technology Data Exchange (ETDEWEB)

Young, R. P.; Collins, D.; Hazzard, J.; Heath, A. [Department of Earth Sciences, Liverpool University, 4 Brownlow street, UK-0 L69 3GP Liverpool (United Kingdom); Pettitt, W.; Baker, C. [Applied Seismology Consultants LTD, 10 Belmont, Shropshire, UK-S41 ITE Shrewsbury (United Kingdom); Billaux, D.; Cundall, P.; Potyondy, D.; Dedecker, F. [Itasca Consultants S.A., Centre Scientifique A. Moiroux, 64, chemin des Mouilles, F69130 Ecully (France); Svemar, C. [Svensk Karnbranslemantering AB, SKB, Aspo Hard Rock Laboratory, PL 300, S-57295 Figeholm (Sweden); Lebon, P. [ANDRA, Parc de la Croix Blanche, 7, rue Jean Monnet, F-92298 Chatenay-Malabry (France)

2004-07-01

This paper presents current results from work performed within the European Commission project SAFETI. The main objective of SAFETI is to develop and test an innovative 3D numerical modelling procedure that will enable the 3-D simulation of nuclear waste repositories in rock. The modelling code is called AC/DC (Adaptive Continuum/ Dis-Continuum) and is partially based on Itasca Consulting Group's Particle Flow Code (PFC). Results are presented from the laboratory validation study where algorithms and procedures have been developed and tested to allow accurate 'Models for Rock' to be produced. Preliminary results are also presented on the use of AC/DC with parallel processors and adaptive logic. During the final year of the project a detailed model of the Prototype Repository Experiment at SKB's Hard Rock Laboratory will be produced using up to 128 processors on the parallel super computing facility at Liverpool University. (authors)

Exploring tectonomagmatic controls on mid-ocean ridge faulting and morphology with 3-D numerical models

Science.gov (United States)

Howell, S. M.; Ito, G.; Behn, M. D.; Olive, J. A. L.; Kaus, B.; Popov, A.; Mittelstaedt, E. L.; Morrow, T. A.

2016-12-01

Previous two-dimensional (2-D) modeling studies of abyssal-hill scale fault generation and evolution at mid-ocean ridges have predicted that M, the ratio of magmatic to total extension, strongly influences the total slip, spacing, and rotation of large faults, as well as the morphology of the ridge axis. Scaling relations derived from these 2-D models broadly explain the globally observed decrease in abyssal hill spacing with increasing ridge spreading rate, as well as the formation of large-offset faults close to the ends of slow-spreading ridge segments. However, these scaling relations do not explain some higher resolution observations of segment-scale variability in fault spacing along the Chile Ridge and the Mid-Atlantic Ridge, where fault spacing shows no obvious correlation with M. This discrepancy between observations and 2-D model predictions illuminates the need for three-dimensional (3-D) numerical models that incorporate the effects of along-axis variations in lithospheric structure and magmatic accretion. To this end, we use the geodynamic modeling software LaMEM to simulate 3-D tectono-magmatic interactions in a visco-elasto-plastic lithosphere under extension. We model a single ridge segment subjected to an along-axis gradient in the rate of magma injection, which is simulated by imposing a mass source in a plane of model finite volumes beneath the ridge axis. Outputs of interest include characteristic fault offset, spacing, and along-axis gradients in seafloor morphology. We also examine the effects of along-axis variations in lithospheric thickness and off-axis thickening rate. The main objectives of this study are to quantify the relative importance of the amount of magmatic extension and the local lithospheric structure at a given along-axis location, versus the importance of along-axis communication of lithospheric stresses on the 3-D fault evolution and morphology of intermediate-spreading-rate ridges.

Gis-Based Smart Cartography Using 3d Modeling

Science.gov (United States)

Malinverni, E. S.; Tassetti, A. N.

2013-08-01

3D City Models have evolved to be important tools for urban decision processes and information systems, especially in planning, simulation, analysis, documentation and heritage management. On the other hand existing and in use numerical cartography is often not suitable to be used in GIS because not geometrically and topologically correctly structured. The research aim is to 3D structure and organize a numeric cartography for GIS and turn it into CityGML standardized features. The work is framed around a first phase of methodological analysis aimed to underline which existing standard (like ISO and OGC rules) can be used to improve the quality requirement of a cartographic structure. Subsequently, from this technical specifics, it has been investigated the translation in formal contents, using an owner interchange software (SketchUp), to support some guide lines implementations to generate a GIS3D structured in GML3. It has been therefore predisposed a test three-dimensional numerical cartography (scale 1:500, generated from range data captured by 3D laser scanner), tested on its quality according to the previous standard and edited when and where necessary. Cad files and shapefiles are converted into a final 3D model (Google SketchUp model) and then exported into a 3D city model (CityGML LoD1/LoD2). The GIS3D structure has been managed in a GIS environment to run further spatial analysis and energy performance estimate, not achievable in a 2D environment. In particular geometrical building parameters (footprint, volume etc.) are computed and building envelop thermal characteristics are derived from. Lastly, a simulation is carried out to deal with asbestos and home renovating charges and show how the built 3D city model can support municipal managers with risk diagnosis of the present situation and development of strategies for a sustainable redevelop.

Numerical difficulties to obtain 3-d critical exponents from platonic solids

International Nuclear Information System (INIS)

Alcaraz, F.C.; Herrmann, H.J.

1985-01-01

The possibility to extract critical exponents of 3-d systems exploring the mass gap amplitudes of platonic solids is tested. For the Ising model the proposed method does not work for numerical reasons. (Author) [pt

Brayton Cycle Numerical Modeling using the RELAP5-3D code, version 4.3.4

Energy Technology Data Exchange (ETDEWEB)

Longhini, Eduardo P.; Lobo, Paulo D.C.; Guimarães, Lamartine N.F.; Filho, Francisco A.B.; Ribeiro, Guilherme B., E-mail: edu_longhini@yahoo.com.br [Instituto de Estudos Avançados (IEAv), São José dos Campos, SP (Brazil). Divisão de Energia Nuclear

2017-07-01

This work contributes to enable and develop technologies to mount fast micro reactors, to generate heat and electric energy, for the purpose to warm and to supply electrically spacecraft equipment and, also, the production of nuclear space propulsion effect. So, for this purpose, the Brayton Cycle demonstrates to be an optimum approach for space nuclear power. The Brayton thermal cycle gas has as characteristic to be a closed cycle, with two adiabatic processes and two isobaric processes. The components performing the cycle's processes are compressor, turbine, heat source, cold source and recuperator. Therefore, the working fluid's mass flow runs the thermal cycle that converts thermal energy into electrical energy, able to use in spaces and land devices. The objective is numerically to model the Brayton thermal cycle gas on nominal operation with one turbomachine composed for a radial-inflow compressor and turbine of a 40.8 kWe Brayton Rotating Unit (BRU). The Brayton cycle numerical modeling is being performed with the program RELAP5-3D, version 4.3.4. The nominal operation uses as working fluid a mixture 40 g/mole He-Xe with a flow rate of 1.85 kg/s, shaft rotational speed of 45 krpm, compressor and turbine inlet temperature of 400 K and 1149 K, respectively, and compressor exit pressure 0.931 MPa. Then, the aim is to get physical corresponding data to operate each cycle component and the general cycle on this nominal operation. (author)

Brayton Cycle Numerical Modeling using the RELAP5-3D code, version 4.3.4

International Nuclear Information System (INIS)

Longhini, Eduardo P.; Lobo, Paulo D.C.; Guimarães, Lamartine N.F.; Filho, Francisco A.B.; Ribeiro, Guilherme B.

2017-01-01

This work contributes to enable and develop technologies to mount fast micro reactors, to generate heat and electric energy, for the purpose to warm and to supply electrically spacecraft equipment and, also, the production of nuclear space propulsion effect. So, for this purpose, the Brayton Cycle demonstrates to be an optimum approach for space nuclear power. The Brayton thermal cycle gas has as characteristic to be a closed cycle, with two adiabatic processes and two isobaric processes. The components performing the cycle's processes are compressor, turbine, heat source, cold source and recuperator. Therefore, the working fluid's mass flow runs the thermal cycle that converts thermal energy into electrical energy, able to use in spaces and land devices. The objective is numerically to model the Brayton thermal cycle gas on nominal operation with one turbomachine composed for a radial-inflow compressor and turbine of a 40.8 kWe Brayton Rotating Unit (BRU). The Brayton cycle numerical modeling is being performed with the program RELAP5-3D, version 4.3.4. The nominal operation uses as working fluid a mixture 40 g/mole He-Xe with a flow rate of 1.85 kg/s, shaft rotational speed of 45 krpm, compressor and turbine inlet temperature of 400 K and 1149 K, respectively, and compressor exit pressure 0.931 MPa. Then, the aim is to get physical corresponding data to operate each cycle component and the general cycle on this nominal operation. (author)

Electro-thermal analysis and optimisation of edge termination of power diode supported by 2-D/3-D numerical modelling and simulation

International Nuclear Information System (INIS)

Príbytný, P; Donoval, D; Chvála, A; Marek, J; Molnár, M

2014-01-01

Numerical modelling and simulation provide an efficient tool for analysis and optimization of device structure design. In this paper we present the analysis and the geometry optimization of the power module with high power pin diode structure supported by the advanced 2-D/3-D mixed-mode electro-thermal device simulation. The structure under investigation is P + NN + power diode device designed for high reverse voltages and very high forward currents, with a maximum forward surge current up to 2.7 kA.

Calculation of residual stresses by means of a 3D numerical weld simulation

International Nuclear Information System (INIS)

Nicak, Tomas; Huemmer, Matthias

2008-01-01

The numerical weld simulation has developed very fast in recent years. The problem complexity has increased from simple 2D models to full 3D models, which can describe the entire welding process more realistically. As recent research projects indicate, a quantitative assessment of the residual stresses by means of a 3D analysis is possible. The structure integrity can be assessed based on the weld simulation results superimposed with the operating load. Moreover, to support the qualification of welded components parametric studies for optimization of the residual stress distribution in the weld region can be performed. In this paper a full 3D numerical weld simulation for a man-hole drainage nozzle in a steam generator will be presented. The residual stresses are calculated by means of an uncoupled transient thermal and mechanical FE analysis. The paper will present a robust procedure allowing reasonable predictions of the residual stresses for complex structures in industrial practice. (authors)

Numerical and experimental investigation of the 3D free surface flow in a model Pelton turbine

International Nuclear Information System (INIS)

Fiereder, R; Riemann, S; Schilling, R

2010-01-01

This investigation focuses on the numerical and experimental analysis of the 3D free surface flow in a Pelton turbine. In particular, two typical flow conditions occurring in a full scale Pelton turbine - a configuration with a straight inlet as well as a configuration with a 90 degree elbow upstream of the nozzle - are considered. Thereby, the effect of secondary flow due to the 90 degree bending of the upstream pipe on the characteristics of the jet is explored. The hybrid flow field consists of pure liquid flow within the conduit and free surface two component flow of the liquid jet emerging out of the nozzle into air. The numerical results are validated against experimental investigations performed in the laboratory of the Institute of Fluid Mechanics (FLM). For the numerical simulation of the flow the in-house unstructured fully parallelized finite volume solver solver3D is utilized. An advanced interface capturing model based on the classic Volume of Fluid method is applied. In order to ensure sharp interface resolution an additional convection term is added to the transport equation of the volume fraction. A collocated variable arrangement is used and the set of non-linear equations, containing fluid conservation equations and model equations for turbulence and volume fraction, are solved in a segregated manner. For pressure-velocity coupling the SIMPLE and PISO algorithms are implemented. Detailed analysis of the observed flow patterns in the jet and of the jet geometry are presented.

Numerical and experimental investigation of the 3D free surface flow in a model Pelton turbine

Energy Technology Data Exchange (ETDEWEB)

Fiereder, R; Riemann, S; Schilling, R, E-mail: fiereder@lhm.mw.tum.d [Department of Fluid Mechanics, Technische Universitaet Muenchen Bolzmannstrasse 15, Garching, 85748 (Germany)

2010-08-15

This investigation focuses on the numerical and experimental analysis of the 3D free surface flow in a Pelton turbine. In particular, two typical flow conditions occurring in a full scale Pelton turbine - a configuration with a straight inlet as well as a configuration with a 90 degree elbow upstream of the nozzle - are considered. Thereby, the effect of secondary flow due to the 90 degree bending of the upstream pipe on the characteristics of the jet is explored. The hybrid flow field consists of pure liquid flow within the conduit and free surface two component flow of the liquid jet emerging out of the nozzle into air. The numerical results are validated against experimental investigations performed in the laboratory of the Institute of Fluid Mechanics (FLM). For the numerical simulation of the flow the in-house unstructured fully parallelized finite volume solver solver3D is utilized. An advanced interface capturing model based on the classic Volume of Fluid method is applied. In order to ensure sharp interface resolution an additional convection term is added to the transport equation of the volume fraction. A collocated variable arrangement is used and the set of non-linear equations, containing fluid conservation equations and model equations for turbulence and volume fraction, are solved in a segregated manner. For pressure-velocity coupling the SIMPLE and PISO algorithms are implemented. Detailed analysis of the observed flow patterns in the jet and of the jet geometry are presented.

Numerical and experimental investigation of the 3D free surface flow in a model Pelton turbine

Science.gov (United States)

Fiereder, R.; Riemann, S.; Schilling, R.

2010-08-01

This investigation focuses on the numerical and experimental analysis of the 3D free surface flow in a Pelton turbine. In particular, two typical flow conditions occurring in a full scale Pelton turbine - a configuration with a straight inlet as well as a configuration with a 90 degree elbow upstream of the nozzle - are considered. Thereby, the effect of secondary flow due to the 90 degree bending of the upstream pipe on the characteristics of the jet is explored. The hybrid flow field consists of pure liquid flow within the conduit and free surface two component flow of the liquid jet emerging out of the nozzle into air. The numerical results are validated against experimental investigations performed in the laboratory of the Institute of Fluid Mechanics (FLM). For the numerical simulation of the flow the in-house unstructured fully parallelized finite volume solver solver3D is utilized. An advanced interface capturing model based on the classic Volume of Fluid method is applied. In order to ensure sharp interface resolution an additional convection term is added to the transport equation of the volume fraction. A collocated variable arrangement is used and the set of non-linear equations, containing fluid conservation equations and model equations for turbulence and volume fraction, are solved in a segregated manner. For pressure-velocity coupling the SIMPLE and PISO algorithms are implemented. Detailed analysis of the observed flow patterns in the jet and of the jet geometry are presented.

Minimal duality breaking in the Kallen-Lehman approach to 3D Ising model: A numerical test

International Nuclear Information System (INIS)

Astorino, Marco; Canfora, Fabrizio; Martinez, Cristian; Parisi, Luca

2008-01-01

A Kallen-Lehman approach to 3D Ising model is analyzed numerically both at low and high temperatures. It is shown that, even assuming a minimal duality breaking, one can fix three parameters of the model to get a very good agreement with the Monte Carlo results at high temperatures. With the same parameters the agreement is satisfactory both at low and near critical temperatures. How to improve the agreement with Monte Carlo results by introducing a more general duality breaking is shortly discussed

Landslide/reservoir interaction: 3D numerical modelling of the Vajont rockslide and generated water wave

Science.gov (United States)

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

Magnetic materials and 3D finite element modeling

CERN Document Server

Bastos, Joao Pedro A

2014-01-01

Magnetic Materials and 3D Finite Element Modeling explores material characterization and finite element modeling (FEM) applications. This book relates to electromagnetic analysis based on Maxwell’s equations and application of the finite element (FE) method to low frequency devices. A great source for senior undergraduate and graduate students in electromagnetics, it also supports industry professionals working in magnetics, electromagnetics, ferromagnetic materials science and electrical engineering. The authors present current concepts on ferromagnetic material characterizations and losses. They provide introductory material; highlight basic electromagnetics, present experimental and numerical modeling related to losses and focus on FEM applied to 3D applications. They also explain various formulations, and discuss numerical codes.

Numerical modelling 2 D and 3 D of circulating fluidized bed: application to the realization of regime diagrams; Modelisation numerique 2D et 3D de lit fluidise circulant: application a la realisation du diagramme des regimes

Energy Technology Data Exchange (ETDEWEB)

Begis, J.; Balzer, G.

1997-02-01

The numerical modelling of internal CFB boilers flows faced with complex phenomenons due to the flows un-stationariness, the heterogeneousness of the particle size distribution, and interactions between the two phases and the walls. Our study consisted in applying numerical models to the experimental configuration of cold circulating fluidized bed studied at the Cerchar. Special attention was given to the analysis of particles - wall interactions models, stemming from Jenkins (1992) and Louge`s (1994) theories, as well as the influence of the particles on fluid turbulence. In order to realize numerical simulations, we have used Eulerian two-phases flow codes developed at NHL medolif(2D), ESTET-ASTRID(3D). From different tests we have deducted that the most appropriate model for the realization of CFB`s prediction is the model taking in account the influence of particles on fluid turbulence. Then, to evaluate the validity limits of this model, we have built the regime diagram, and we have compared it with the experimental diagram. We have concluded that the simulation allows to describe the different CFB`s working regimes, and especially transitions. We have also noticed the importance of the choice of the mean diameter of the simulated particles. In this way, making a correction of the simulated particles` diameter in comparison with Sauter mean particle diameter, we obtained numerical results in good agreement with experimental data. (authors) 13 refs.

3D numerical simulation and analysis of railgun gouging mechanism

Directory of Open Access Journals (Sweden)

Jin-guo Wu

2016-04-01

Full Text Available A gouging phenomenon with a hypervelocity sliding electrical contact in railgun not only shortens the rail lifetime but also affects the interior ballistic performance. In this paper, a 3-D numerical model was introduced to simulate and analyze the generation mechanism and evolution of the rail gouging phenomenon. The results show that a rail surface bulge is an important factor to induce gouging. High density and high pressure material flow on the contact surface, obliquely extruded into the rail when accelerating the armature to a high velocity, can produce gouging. Both controlling the bulge size to a certain range and selecting suitable materials for rail surface coating will suppress the formation of gouging. The numerical simulation had a good agreement with experiments, which validated the computing model and methodology are reliable.

3D Numerical Rift Modeling with Application to the East African Rift System

Science.gov (United States)

Glerum, A.; Brune, S.; Naliboff, J.

2017-12-01

As key components of plate tectonics, continental rifting and the formation of passive margins have been extensively studied with both analogue models and numerical techniques. Only recently however, technical advances have enabled numerical investigations into rift evolution in three dimensions, as is actually required for including those processes that cause rift-parallel variability, such as structural inheritance and oblique extension (Brune 2016). We use the massively parallel finite element code ASPECT (Kronbichler et al. 2012; Heister et al. 2017) to investigate rift evolution. ASPECT's adaptive mesh refinement enables us to focus resolution on the regions of interest (i.e. the rift center), while leaving other areas such as the asthenospheric mantle at coarse resolution, leading to kilometer-scale local mesh resolution in 3D. Furthermore, we implemented plastic and viscous strain weakening of the nonlinear viscoplastic rheology required to develop asymmetric rift geometries (e.g. Huismans and Beaumont 2003). Additionally created plugins to ASPECT allow us to specify initial temperature and composition conditions based on geophysical data (e.g. LITHO1.0, Pasyanos et al. 2014) or to prescribe more general along-strike variation in the initial strain seeding the rift. Employing the above functionality, we construct regional models of the East African Rift System (EARS), the world's largest currently active rift. As the EARS is characterized by both orthogonal and oblique rift sections, multi-phase extension histories as well as magmatic and a-magmatic branches (e.g. Chorowicz 2005; Ebinger and Scholz 2011), it constitutes an extensive natural laboratory for our research into the 3D nature of continental rifting. References:Brune, S. (2016), in Plate boundaries and natural hazards, AGU Geophysical Monograph 219, J. C. Duarte and W. P. Schellart (Eds.). Chorowicz, J. (2005). J. Afr. Earth Sci., 43, 379-410. Ebinger, C. and Scholz, C. A. (2011), in Tectonics of

Numeric 3D modelling of a geohydrothermal duplet system in malm karst; Numerische 3D-Modellierung eines geohydrothermalen Dublettenbetriebs im Malmkarst

Energy Technology Data Exchange (ETDEWEB)

Wenderoth, F.; Gropius, M. [FlowFM GbR, Berlin (Germany); Fritzer, T. [Bayerisches Geologisches Landesamt, Muenchen (Germany); Huber, B. [HydroConsult GmbH, Augsburg (Germany); Schubert, A. [ERDWERK GmbH, Muenchen (Germany)

2005-10-01

The Bavarian Geological State Office charged the HYDRO Consult GmbH, Augsburg, and the FlowFM GbR, Berlin with the task of investigating the possibilibies of numeric 3D modelling using the FEFLOW {sup registered} developed by WASY GmbH incorporating Darcy's porosity equation. Preliminary data acquisition and parallel measurements were carried out by ERDWERK GmbH on behalf of SWM-Versorgungs-GmbH, GB Energie-Erzeugung. The investigations were to provide information on whether this method can provide reliable data for licensing procedures under mining law. In particular, the following questions were to be answered: a) maximum injection and production pressures for a given operating regime; b) regional hydraulic effects of geothermal plant operation at a given volume flow and injection temperature; c) temporal and regional cold water distribution as a result of injection; d) time requirements for heat regeneration after terminaiton of geothermal duplet operation. The research project was funded by the Bavarian state minister of economics, infrastructure, transportation and technology. (orig.)

Implementation of virtual models from sheet metal forming simulation into physical 3D colour models using 3D printing

Science.gov (United States)

Junk, S.

2016-08-01

Today the methods of numerical simulation of sheet metal forming offer a great diversity of possibilities for optimization in product development and in process design. However, the results from simulation are only available as virtual models. Because there are any forming tools available during the early stages of product development, physical models that could serve to represent the virtual results are therefore lacking. Physical 3D-models can be created using 3D-printing and serve as an illustration and present a better understanding of the simulation results. In this way, the results from the simulation can be made more “comprehensible” within a development team. This paper presents the possibilities of 3D-colour printing with particular consideration of the requirements regarding the implementation of sheet metal forming simulation. Using concrete examples of sheet metal forming, the manufacturing of 3D colour models will be expounded upon on the basis of simulation results.

Numerical Study on the Tensile Behavior of 3D Four Directional Cylindrical Braided Composite Shafts

Science.gov (United States)

Zhao, Guoqi; Wang, Jiayi; Hao, Wenfeng; Liu, Yinghua; Luo, Ying

2017-10-01

The tensile behavior of 3D four directional cylindrical braided composite shafts was analyzed with the numerical method. The unit cell models for the 3D four directional cylindrical braided composite shafts with various braiding angles were constructed with ABAQUS. Hashin's failure criterion was used to analyze the tensile strength and the damage evolution of the unit cells. The influence of the braiding angle on the tensile behavior of the 3D four directional cylindrical braided composite shafts was analyzed. The numerical results showed that the tensile strength along the braiding direction increased as the braiding angle decreased. These results should play an integral role in the design of braiding composites shafts.

3-D numerical modelling of groundwater flow for scenario-based ...

African Journals Online (AJOL)

The data related to piezometric water levels, canal gauges, well logs, meteorological and lithological information were collected from Punjab Irrigation Department (PID), Water and Power Development Authority (WAPDA). Groundwater flow models for both steady and transient conditions were set-up using FEFLOW-3D.

The Effect of Numerical 2D and 3D Fem Element Modelling on Strain and Stress Distributions at Laser Weld Notches in Steel Sandwich Type Panels

Directory of Open Access Journals (Sweden)

Niklas Karol

2018-03-01

Full Text Available Like other means of transport, merchant ships face the problem of increasing requirements concerning the environment protection, which, among other issues, implies the reduction of fuel consumption by the ship. Here, the conventional approach which consists in making use of higher strength steels to decrease the mass of the ship hull can be complemented by the use of new steel structures of sandwich panel type. However, the lack of knowledge and experience concerning, among other issues, fatigue strength assessment of thin-walled sandwich structures makes their use limited. Untypical welds imply the need for individual approach to the fatigue analysis. The article presents the effect of numerical FEM modelling with the aid of two-dimensional (2D and three-dimensional (3D elements on the results of strain and stress distributions in the areas of toe and root notches of the analysed laser weld. The presented results of computer simulation reveal that modelling of strain and stress states in 2D (instead of full 3D affects only the results in close vicinity of the notch, and the observed differences rapidly disappear at a distance of 0.05 mm from the bottom of the notch. The obtained results confirm the possibility of use of numerically effective 2D strain and stress state models for analysing the fatigue strength of laser weld according to local approach.

Numerical Simulation of Effective Properties of 3D Piezoelectric Composites

Directory of Open Access Journals (Sweden)

Ri-Song Qin

2014-01-01

Full Text Available The prediction of the overall effective properties of fibre-reinforced piezocomposites has drawn much interest from investigators recently. In this work, an algorithm used in two-dimensional (2D analysis for calculating transversely isotropic material properties is developed. Since the finite element (FE meshing patterns on the opposite areas are the same, constraint equations can be applied directly to generate appropriate load. The numerical results derived using this model have found a good agreement with those in the literature. The 2D algorithm is then modified and improved in such a way that it is valid for three-dimensional (3D analysis in the case of random distributed shorts and inclusions. Linear interpolation of displacement field is employed to establish constraint equations of nodal displacements between two adjacent elements.

3D Geothermal Modelling of the Mount Amiata Hydrothermal System in Italy

Directory of Open Access Journals (Sweden)

Paolo Fulignati

2014-11-01

Full Text Available In this paper we build a subsurface model that helps in visualizing and understanding the structural framework, geology and their interactions with the Mt. Amiata geothermal system. Modelling in 3D provides the possibility to interpolate the geometry of structures and is an effective way of understanding geological features. The 3D modelling approach appears to be crucial for further progress in the reconstruction of the assessment of the geothermal model of Mt. Amiata. Furthermore, this model is used as the basis of a 3D numerical thermo-fluid-dynamic model of the existing reservoir(s. The integration between borehole data and numerical modelling results allows reconstructing the temperature distribution in the subsoil of the Mt. Amiata area.

AxiSEM3D: broadband seismic wavefields in 3-D aspherical Earth models

Science.gov (United States)

Leng, K.; Nissen-Meyer, T.; Zad, K. H.; van Driel, M.; Al-Attar, D.

2017-12-01

Seismology is the primary tool for data-informed inference of Earth structure and dynamics. Simulating seismic wave propagation at a global scale is fundamental to seismology, but remains as one of most challenging problems in scientific computing, because of both the multiscale nature of Earth's interior and the observable frequency band of seismic data. We present a novel numerical method to simulate global seismic wave propagation in realistic 3-D Earth models. Our method, named AxiSEM3D, is a hybrid of spectral element method and pseudospectral method. It reduces the azimuthal dimension of wavefields by means of a global Fourier series parameterization, of which the number of terms can be locally adapted to the inherent azimuthal smoothness of the wavefields. AxiSEM3D allows not only for material heterogeneities, such as velocity, density, anisotropy and attenuation, but also for finite undulations on radial discontinuities, both solid-solid and solid-fluid, and thereby a variety of aspherical Earth features such as ellipticity, topography, variable crustal thickness, and core-mantle boundary topography. Such interface undulations are equivalently interpreted as material perturbations of the contiguous media, based on the "particle relabelling transformation". Efficiency comparisons show that AxiSEM3D can be 1 to 3 orders of magnitude faster than conventional 3-D methods, with the speedup increasing with simulation frequency and decreasing with model complexity, but for all realistic structures the speedup remains at least one order of magnitude. The observable frequency range of global seismic data (up to 1 Hz) has been covered for wavefield modelling upon a 3-D Earth model with reasonable computing resources. We show an application of surface wave modelling within a state-of-the-art global crustal model (Crust1.0), with the synthetics compared to real data. The high-performance C++ code is released at github.com/AxiSEM3D/AxiSEM3D.

3D SPH numerical simulation of the wave generated by the Vajont rockslide

Science.gov (United States)

Vacondio, R.; Mignosa, P.; Pagani, S.

2013-09-01

A 3D numerical modeling of the wave generated by the Vajont slide, one of the most destructive ever occurred, is presented in this paper. A meshless Lagrangian Smoothed Particle Hydrodynamics (SPH) technique was adopted to simulate the highly fragmented violent flow generated by the falling slide in the artificial reservoir. The speed-up achievable via General Purpose Graphic Processing Units (GP-GPU) allowed to adopt the adequate resolution to describe the phenomenon. The comparison with the data available in literature showed that the results of the numerical simulation reproduce satisfactorily the maximum run-up, also the water surface elevation in the residual lake after the event. Moreover, the 3D velocity field of the flow during the event and the discharge hydrograph which overtopped the dam, were obtained.

Numerical 3D modelling of oil dispersion in the sea due to different accident scenarios

Science.gov (United States)

Guandalini, Roberto; Agate, Giordano; Moia, Fabio

2017-04-01

The purpose of the study has been the development of a methodology, based on a numerical 3D approach, for the analysis of oil dispersion in the sea, in order to simulate with a high level of accuracy the dynamic behavior of the oil plume and its displacement in the environment. As a matter of fact, the numerical simulation is the only approach currently able to analyse in detail possible accident scenarios, even with an high degree of complexity, of different type and intensity, allowing to follow their evolution both in time and space, and to evaluate the effectiveness of suggested prevention or recovery actions. The software for these calculations is therefore an essential tool in order to simulate the impact effects in the short, medium and long period, able to account for the complexity of the sea system involved in the dispersion process and its dependency on the meteorological, marine and morphological local conditions. This software, generally based on fluid dynamic 3D simulators and modellers, is therefore extremely specialized and requires expertise for an appropriate usage, but at the same time it allows detailed scenario analyses and design verifications. It takes into account different parameters as the sea current field and its turbulence, the wind acting on the sea surface, the salinity and temperature gradients, the local coastal morphology, the seabed bathymetry and the tide. The applied methodology is based on the Integrated Fluid Dynamic Simulation System HyperSuite developed by RSE. This simulation system includes the consideration of all the parameters previously listed, in the frame of a 3D Eulerian finite element fluid dynamic model, which accuracy is guaranteed by a very detailed spatial mesh and by an automatically optimized time step management. In order to assess the methodology features, an area of more than 2500 km2 and depth of 200 m located in the middle Adriatic Sea has been modelled. The information required for the simulation in

3D Geological Model for "LUSI" - a Deep Geothermal System

Science.gov (United States)

Sohrabi, Reza; Jansen, Gunnar; Mazzini, Adriano; Galvan, Boris; Miller, Stephen A.

2016-04-01

Geothermal applications require the correct simulation of flow and heat transport processes in porous media, and many of these media, like deep volcanic hydrothermal systems, host a certain degree of fracturing. This work aims to understand the heat and fluid transport within a new-born sedimentary hosted geothermal system, termed Lusi, that began erupting in 2006 in East Java, Indonesia. Our goal is to develop conceptual and numerical models capable of simulating multiphase flow within large-scale fractured reservoirs such as the Lusi region, with fractures of arbitrary size, orientation and shape. Additionally, these models can also address a number of other applications, including Enhanced Geothermal Systems (EGS), CO2 sequestration (Carbon Capture and Storage CCS), and nuclear waste isolation. Fractured systems are ubiquitous, with a wide-range of lengths and scales, making difficult the development of a general model that can easily handle this complexity. We are developing a flexible continuum approach with an efficient, accurate numerical simulator based on an appropriate 3D geological model representing the structure of the deep geothermal reservoir. Using previous studies, borehole information and seismic data obtained in the framework of the Lusi Lab project (ERC grant n°308126), we present here the first 3D geological model of Lusi. This model is calculated using implicit 3D potential field or multi-potential fields, depending on the geological context and complexity. This method is based on geological pile containing the geological history of the area and relationship between geological bodies allowing automatic computation of intersections and volume reconstruction. Based on the 3D geological model, we developed a new mesh algorithm to create hexahedral octree meshes to transfer the structural geological information for 3D numerical simulations to quantify Thermal-Hydraulic-Mechanical-Chemical (THMC) physical processes.

3-D numerical simulations of coronal loops oscillations

Directory of Open Access Journals (Sweden)

M. Selwa

2009-10-01

Full Text Available We present numerical results of 3-D MHD model of a dipole active region field containing a loop with a higher density than its surroundings. We study different ways of excitation of vertical kink oscillations by velocity perturbation: as an initial condition, and as an impulsive excitation with a pulse of a given position, duration, and amplitude. These properties are varied in the parametric studies. We find that the amplitude of vertical kink oscillations is significantly amplified in comparison to horizontal kink oscillations for exciters located centrally (symmetrically below the loop, but not if the exciter is located a significant distance to the side of the loop. This explains why the pure vertical kink mode is so rarely observed in comparison to the horizontally polarized one. We discuss the role of curved magnetic field lines and the pulse overlapping at one of the loop's footpoints in 3-D active regions (AR's on the excitation and the damping of slow standing waves. We find that footpoint excitation becomes more efficient in 3-D curved loops than in 2-D curved arcades and that slow waves can be excited within an interval of time that is comparable to the observed one wave-period due to the combined effect of the pulse inside and outside the loop. Additionally, we study the effect of AR topology on the excitation and trapping of loop oscillations. We find that a perturbation acting directly on a single loop excites oscillations, but results in an increased leakage compared to excitation of oscillations in an AR field by an external source.

3-D numerical simulations of coronal loops oscillations

Directory of Open Access Journals (Sweden)

M. Selwa

2009-10-01

Full Text Available We present numerical results of 3-D MHD model of a dipole active region field containing a loop with a higher density than its surroundings. We study different ways of excitation of vertical kink oscillations by velocity perturbation: as an initial condition, and as an impulsive excitation with a pulse of a given position, duration, and amplitude. These properties are varied in the parametric studies. We find that the amplitude of vertical kink oscillations is significantly amplified in comparison to horizontal kink oscillations for exciters located centrally (symmetrically below the loop, but not if the exciter is located a significant distance to the side of the loop. This explains why the pure vertical kink mode is so rarely observed in comparison to the horizontally polarized one. We discuss the role of curved magnetic field lines and the pulse overlapping at one of the loop's footpoints in 3-D active regions (AR's on the excitation and the damping of slow standing waves. We find that footpoint excitation becomes more efficient in 3-D curved loops than in 2-D curved arcades and that slow waves can be excited within an interval of time that is comparable to the observed one wave-period due to the combined effect of the pulse inside and outside the loop. Additionally, we study the effect of AR topology on the excitation and trapping of loop oscillations. We find that a perturbation acting directly on a single loop excites oscillations, but results in an increased leakage compared to excitation of oscillations in an AR field by an external source.

AUTOMATIC TEXTURE MAPPING OF ARCHITECTURAL AND ARCHAEOLOGICAL 3D MODELS

Directory of Open Access Journals (Sweden)

T. P. Kersten

2012-07-01

Full Text Available Today, detailed, complete and exact 3D models with photo-realistic textures are increasingly demanded for numerous applications in architecture and archaeology. Manual texture mapping of 3D models by digital photographs with software packages, such as Maxon Cinema 4D, Autodesk 3Ds Max or Maya, still requires a complex and time-consuming workflow. So, procedures for automatic texture mapping of 3D models are in demand. In this paper two automatic procedures are presented. The first procedure generates 3D surface models with textures by web services, while the second procedure textures already existing 3D models with the software tmapper. The program tmapper is based on the Multi Layer 3D image (ML3DImage algorithm and developed in the programming language C++. The studies showing that the visibility analysis using the ML3DImage algorithm is not sufficient to obtain acceptable results of automatic texture mapping. To overcome the visibility problem the Point Cloud Painter algorithm in combination with the Z-buffer-procedure will be applied in the future.

Automatic Texture Mapping of Architectural and Archaeological 3d Models

Science.gov (United States)

Kersten, T. P.; Stallmann, D.

2012-07-01

Today, detailed, complete and exact 3D models with photo-realistic textures are increasingly demanded for numerous applications in architecture and archaeology. Manual texture mapping of 3D models by digital photographs with software packages, such as Maxon Cinema 4D, Autodesk 3Ds Max or Maya, still requires a complex and time-consuming workflow. So, procedures for automatic texture mapping of 3D models are in demand. In this paper two automatic procedures are presented. The first procedure generates 3D surface models with textures by web services, while the second procedure textures already existing 3D models with the software tmapper. The program tmapper is based on the Multi Layer 3D image (ML3DImage) algorithm and developed in the programming language C++. The studies showing that the visibility analysis using the ML3DImage algorithm is not sufficient to obtain acceptable results of automatic texture mapping. To overcome the visibility problem the Point Cloud Painter algorithm in combination with the Z-buffer-procedure will be applied in the future.

Numerical modeling techniques for flood analysis

Science.gov (United States)

Anees, Mohd Talha; Abdullah, K.; Nawawi, M. N. M.; Ab Rahman, Nik Norulaini Nik; Piah, Abd. Rahni Mt.; Zakaria, Nor Azazi; Syakir, M. I.; Mohd. Omar, A. K.

2016-12-01

Topographic and climatic changes are the main causes of abrupt flooding in tropical areas. It is the need to find out exact causes and effects of these changes. Numerical modeling techniques plays a vital role for such studies due to their use of hydrological parameters which are strongly linked with topographic changes. In this review, some of the widely used models utilizing hydrological and river modeling parameters and their estimation in data sparse region are discussed. Shortcomings of 1D and 2D numerical models and the possible improvements over these models through 3D modeling are also discussed. It is found that the HEC-RAS and FLO 2D model are best in terms of economical and accurate flood analysis for river and floodplain modeling respectively. Limitations of FLO 2D in floodplain modeling mainly such as floodplain elevation differences and its vertical roughness in grids were found which can be improve through 3D model. Therefore, 3D model was found to be more suitable than 1D and 2D models in terms of vertical accuracy in grid cells. It was also found that 3D models for open channel flows already developed recently but not for floodplain. Hence, it was suggested that a 3D model for floodplain should be developed by considering all hydrological and high resolution topographic parameter's models, discussed in this review, to enhance the findings of causes and effects of flooding.

A Numerical Study of Forbush Decreases with a 3D Cosmic-Ray Modulation Model Based on an SDE Approach

Energy Technology Data Exchange (ETDEWEB)

Luo, Xi; Feng, Xueshang [SIGMA Weather Group, State Key Laboratory of Space Weather, National Space Science Center, Chinese Academy of Sciences, Beijing 100190 (China); Potgieter, Marius S. [Centre for Space Research, North-West University, Potchefstroom 2520 (South Africa); Zhang, Ming [Department of Physics and Space Sciences, Florida Institute of Technology, 150 West University Boulevard, Melbourne, FL 32901 (United States)

2017-04-10

Based on the reduced diffusion mechanism for producing Forbush decreases (Fds) in the heliosphere, we constructed a three-dimensional (3D) diffusion barrier, and by incorporating it into a stochastic differential equation (SDE) based time-dependent, cosmic-ray transport model, a 3D numerical model for simulating Fds is built and applied to a period of relatively quiet solar activity. This SDE model generally corroborates previous Fd simulations concerning the effects of the solar magnetic polarity, the tilt angle of the heliospheric current sheet (HCS), and cosmic-ray particle energy. Because the modulation processes in this 3D model are multi-directional, the barrier’s geometrical features affect the intensity profiles of Fds differently. We find that both the latitudinal and longitudinal extent of the barrier have relatively fewer effects on these profiles than its radial extent and the level of decreased diffusion inside the disturbance. We find, with the 3D approach, that the HCS rotational motion causes the relative location from the observation point to the HCS to vary, so that a periodic pattern appears in the cosmic-ray intensity at the observing location. Correspondingly, the magnitude and recovery time of an Fd change, and the recovering intensity profile contains oscillation as well. Investigating the Fd magnitude variation with heliocentric radial distance, we find that the magnitude decreases overall and, additionally, that the Fd magnitude exhibits an oscillating pattern as the radial distance increases, which coincides well with the wavy profile of the HCS under quiet solar modulation conditions.

Numerical Investigation of T-joints with 3D Four Directional Braided Composite Fillers Under Tensile Loading

Science.gov (United States)

Li, Xiao-kang; Liu, Zhen-guo; Hu, Long; Wang, Yi-bo; Lei, Bing; Huang, Xiang

2017-02-01

Numerical studied on T-joints with three-dimensional four directional (3D4D) braided composite fillers was presented in this article. Compared with conventional unidirectional prepreg fillers, the 3D braided composite fillers have excellent ability to prevent crack from penetrating trigone fillers, which constantly occurred in the conventional fillers. Meanwhile, the 3D braided composite fillers had higher fiber volume fraction and eliminated the fiber folding problem in unidirectional prepreg fillers. The braiding technology and mechanical performance of 3D4D braided fillers were studied. The numerical model of carbon fiber T-joints with 3D4D braided composite fillers was built by finite element analysis software. The damage formation, extension and failing process of T-joints with 3D4D braided fillers under tensile load were investigated. Further investigation was extended to the effect of 3D4D braided fillers with different braiding angles on mechanical behavior of the T-joints. The study results revealed that the filling area was the weakest part of the T-joints where the damage first appeared and the crack then rapidly spread to the glue film around the filling area and the interface between over-laminate and soleplate. The 3D4D braided fillers were undamaged and the braiding angle change induced a little effect on the bearing capacity of T-joints.

Local seismic hazard assessment in explosive volcanic settings by 3D numerical analyses

Science.gov (United States)

Razzano, Roberto; Pagliaroli, Alessandro; Moscatelli, Massimiliano; Gaudiosi, Iolanda; Avalle, Alessandra; Giallini, Silvia; Marcini, Marco; Polpetta, Federica; Simionato, Maurizio; Sirianni, Pietro; Sottili, Gianluca; Vignaroli, Gianluca; Bellanova, Jessica; Calamita, Giuseppe; Perrone, Angela; Piscitelli, Sabatino

2017-04-01

This work deals with the assessment of local seismic response in the explosive volcanic settings by reconstructing the subsoil model of the Stracciacappa maar (Sabatini Volcanic District, central Italy), whose pyroclastic succession records eruptive phases ended about 0.09 Ma ago. Heterogeneous characteristics of the Stracciacappa maar (stratification, structural setting, lithotypes, and thickness variation of depositional units) make it an ideal case history for understanding mechanisms and processes leading to modifications of amplitude-frequency-duration of seismic waves generated at earthquake sources and propagating through volcanic settings. New geological map and cross sections, constrained with recently acquired geotechnical and geophysical data, illustrate the complex geometric relationships among different depositional units forming the maar. A composite interfingering between internal lacustrine sediments and epiclastic debris, sourced from the rim, fills the crater floor; a 45 meters thick continuous coring borehole was drilled in the maar with sampling of undisturbed samples. Electrical Resistivity Tomography surveys and 2D passive seismic arrays were also carried out for constraining the geological model and the velocity profile of the S-waves, respectively. Single station noise measurements were collected in order to define natural amplification frequencies. Finally, the nonlinear cyclic soil behaviour was investigated through simple shear tests on the undisturbed samples. The collected dataset was used to define the subsoil model for 3D finite difference site response numerical analyses by using FLAC 3D software (ITASCA). Moreover, 1D and 2D numerical analyses were carried out for comparison purposes. Two different scenarios were selected as input motions: a moderate magnitude (volcanic event) and a high magnitude (tectonic event). Both earthquake scenarios revealed significant ground motion amplification (up to 15 in terms of spectral acceleration

Integrated Biogeomorphological Modeling Using Delft3D

Science.gov (United States)

Ye, Q.; Jagers, B.

2011-12-01

The skill of numerical morphological models has improved significantly from the early 2D uniform, total load sediment models (with steady state or infrequent wave updates) to recent 3D hydrodynamic models with multiple suspended and bed load sediment fractions and bed stratigraphy (online coupled with waves). Although there remain many open questions within this combined field of hydro- and morphodynamics, we observe an increasing need to include biological processes in the overall dynamics. In riverine and inter-tidal environments, there is often an important influence by riparian vegetation and macrobenthos. Over the past decade more and more researchers have started to extend the simulation environment with wrapper scripts and other quick code hacks to estimate their influence on morphological development in coastal, estuarine and riverine environments. Although one can in this way quickly analyze different approaches, these research tools have generally not been designed with reuse, performance and portability in mind. We have now implemented a reusable, flexible, and efficient two-way link between the Delft3D open source framework for hydrodynamics, waves and morphology, and the water quality and ecology modules. The same link will be used for 1D, 2D and 3D modeling on networks and both structured and unstructured grids. We will describe the concepts of the overall system, and illustrate it with some first results.

Experimental validation of a numerical 3-D finite model applied to wind turbines design under vibration constraints: TREVISE platform

Science.gov (United States)

Sellami, Takwa; Jelassi, Sana; Darcherif, Abdel Moumen; Berriri, Hanen; Mimouni, Med Faouzi

2018-04-01

With the advancement of wind turbines towards complex structures, the requirement of trusty structural models has become more apparent. Hence, the vibration characteristics of the wind turbine components, like the blades and the tower, have to be extracted under vibration constraints. Although extracting the modal properties of blades is a simple task, calculating precise modal data for the whole wind turbine coupled to its tower/foundation is still a perplexing task. In this framework, this paper focuses on the investigation of the structural modeling approach of modern commercial micro-turbines. Thus, the structural model a complex designed wind turbine, which is Rutland 504, is established based on both experimental and numerical methods. A three-dimensional (3-D) numerical model of the structure was set up based on the finite volume method (FVM) using the academic finite element analysis software ANSYS. To validate the created model, experimental vibration tests were carried out using the vibration test system of TREVISE platform at ECAM-EPMI. The tests were based on the experimental modal analysis (EMA) technique, which is one of the most efficient techniques for identifying structures parameters. Indeed, the poles and residues of the frequency response functions (FRF), between input and output spectra, were calculated to extract the mode shapes and the natural frequencies of the structure. Based on the obtained modal parameters, the numerical designed model was up-dated.

The Engelbourg's ruins: from 3D TLS point cloud acquisition to 3D virtual and historic models

Science.gov (United States)

Koehl, Mathieu; Berger, Solveig; Nobile, Sylvain

2014-05-01

The Castle of Engelbourg was built at the beginning of the 13th century, at the top of the Schlossberg. It is situated on the territory of the municipality of Thann (France), at the crossroads of Alsace and Lorraine, and dominates the outlet of the valley of Thur. Its strategic position was one of the causes of its systematic destructions during the 17th century, and Louis XIV finished his fate by ordering his demolition in 1673. Today only few vestiges remain, of which a section of the main tower from about 7m of diameter and 4m of wide laying on its slice, unique characteristic in the regional castral landscape. It is visible since the valley, was named "the Eye of the witch", and became a key attraction of the region. The site, which extends over approximately one hectare, is for several years the object of numerous archaeological studies and is at the heart of a project of valuation of the vestiges today. It was indeed a key objective, among the numerous planned works, to realize a 3D model of the site in its current state, in other words, a virtual model "such as seized", exploitable as well from a cultural and tourist point of view as by scientists and in archaeological researches. The team of the ICube/INSA lab had in responsibility the realization of this model, the acquisition of the data until the delivery of the virtual model, thanks to 3D TLS and topographic surveying methods. It was also planned to integrate into this 3D model, data of 2D archives, stemming from series of former excavations. The objectives of this project were the following ones: • Acquisition of 3D digital data of the site and 3D modelling • Digitization of the 2D archaeological data and integration in the 3D model • Implementation of a database connected to the 3D model • Virtual Visit of the site The obtained results allowed us to visualize every 3D object individually, under several forms (point clouds, 3D meshed objects and models, etc.) and at several levels of detail

Numerical investigation of band gaps in 3D printed cantilever-in-mass metamaterials

Science.gov (United States)

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.

Micro-mechanics based damage mechanics for 3D Orthogonal Woven Composites: Experiment and Numerical Modelling

KAUST Repository

Saleh, Mohamed Nasr; Lubineau, Gilles; Potluri, Prasad; Withers, Philip; Soutis, Constantinos

2016-01-01

Damage initiation and evolution of three-dimensional (3D) orthogonal woven carbon fibre composite (3DOWC) is investigated experimentally and numerically. Meso-scale homogenisation of the representative volume element (RVE) is utilised to predict the elastic properties, simulate damage initiation and evolution when loaded in tension. The effect of intra-yarns transverse cracking and shear diffused damage on the in-plane transverse modulus and shear modulus is investigated while one failure criterion is introduced to simulate the matrix damage. The proposed model is based on two major assumptions. First, the effect of the binder yarns, on the in-plane properties, is neglected, so the 3DOWC unit cell can be approximated as a (0o/90o) cross-ply laminate. Second, a micro-mechanics based damage approach is used at the meso-scale, so damage indicators can be correlated, explicitly, to the density of cracks within the material. Results from the simulated RVE are validated against experimental results along the warp (0o direction) and weft (90o direction). This approach paves the road for more predictive models as damage evolution laws are obtained from micro mechanical considerations and rely on few well-defined material parameters. This largely differs from classical damage mechanics approaches in which the evolution law is obtained by retrofitting experimental observations.

Micro-mechanics based damage mechanics for 3D Orthogonal Woven Composites: Experiment and Numerical Modelling

KAUST Repository

Saleh, Mohamed Nasr

2016-01-08

Damage initiation and evolution of three-dimensional (3D) orthogonal woven carbon fibre composite (3DOWC) is investigated experimentally and numerically. Meso-scale homogenisation of the representative volume element (RVE) is utilised to predict the elastic properties, simulate damage initiation and evolution when loaded in tension. The effect of intra-yarns transverse cracking and shear diffused damage on the in-plane transverse modulus and shear modulus is investigated while one failure criterion is introduced to simulate the matrix damage. The proposed model is based on two major assumptions. First, the effect of the binder yarns, on the in-plane properties, is neglected, so the 3DOWC unit cell can be approximated as a (0o/90o) cross-ply laminate. Second, a micro-mechanics based damage approach is used at the meso-scale, so damage indicators can be correlated, explicitly, to the density of cracks within the material. Results from the simulated RVE are validated against experimental results along the warp (0o direction) and weft (90o direction). This approach paves the road for more predictive models as damage evolution laws are obtained from micro mechanical considerations and rely on few well-defined material parameters. This largely differs from classical damage mechanics approaches in which the evolution law is obtained by retrofitting experimental observations.

Numerical prediction of pressure loss in tight-lattice rod bundle by use of 3-dimensional two-fluid model simulation code ACE-3D

International Nuclear Information System (INIS)

Yoshida, Hiroyuki; Takase, Kazuyuki; Suzuki, Takayuki

2009-01-01

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 the ACE-3D. 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 the ACE-3D, numerical simulation of boiling two-phase in tight-lattice rod bundle with no lift force model was also performed. In the results, the lift force model has direct effects on void fraction concentration in gap region, and pressure distribution in horizontal plane induced by void fraction distribution cause of bubble movement from the gap region to the subchannel region. The predicted pressure loss in the section that includes no spacer accorded with experimental results with around 10% of differences. The predicted friction pressure loss was underestimated around 20% of measured values, and the effect of the turbulence model is considered as one of the causes of this underestimation. (author)

Implementation and verification of interface constitutive model in FLAC3D

Directory of Open Access Journals (Sweden)

Hai-min Wu

2011-09-01

Full Text Available Due to the complexity of soil-structure interaction, simple constitutive models typically used for interface elements in general computer programs cannot satisfy the requirements of discontinuous deformation analysis of structures that contain different interfaces. In order to simulate the strain-softening characteristics of interfaces, a nonlinear strain-softening interface constitutive model was incorporated into fast Lagrange analysis of continua in three dimensions (FLAC3D through a user-defined program in the FISH environment. A numerical simulation of a direct shear test for geosynthetic interfaces was conducted to verify that the interface model was implemented correctly. Results of the numerical tests show good agreement with the results obtained from theoretical calculations, indicating that the model incorporated into FLAC3D can simulate the nonlinear strain-softening behavior of interfaces involving geosynthetic materials. The results confirmed the validity and reliability of the improved interface model. The procedure and method of implementing an interface constitutive model into a commercial computer program also provide a reference for implementation of a new interface constitutive model in FLAC3D.

Modeling 3D Objects for Navigation Purposes Using Laser Scanning

Directory of Open Access Journals (Sweden)

Cezary Specht

2016-07-01

Full Text Available The paper discusses the creation of 3d models and their applications in navigation. It contains a review of available methods and geometric data sources, focusing mostly on terrestrial laser scanning. It presents detailed description, from field survey to numerical elaboration, how to construct accurate model of a typical few storey building as a hypothetical reference in complex building navigation. Hence, the paper presents fields where 3d models are being used and their potential new applications.

3D Modelling and Printing Technology to Produce Patient-Specific 3D Models.

Science.gov (United States)

Birbara, Nicolette S; Otton, James M; Pather, Nalini

2017-11-10

A comprehensive knowledge of mitral valve (MV) anatomy is crucial in the assessment of MV disease. While the use of three-dimensional (3D) modelling and printing in MV assessment has undergone early clinical evaluation, the precision and usefulness of this technology requires further investigation. This study aimed to assess and validate 3D modelling and printing technology to produce patient-specific 3D MV models. A prototype method for MV 3D modelling and printing was developed from computed tomography (CT) scans of a plastinated human heart. Mitral valve models were printed using four 3D printing methods and validated to assess precision. Cardiac CT and 3D echocardiography imaging data of four MV disease patients was used to produce patient-specific 3D printed models, and 40 cardiac health professionals (CHPs) were surveyed on the perceived value and potential uses of 3D models in a clinical setting. The prototype method demonstrated submillimetre precision for all four 3D printing methods used, and statistical analysis showed a significant difference (p3D printed models, particularly using multiple print materials, were considered useful by CHPs for preoperative planning, as well as other applications such as teaching and training. This study suggests that, with further advances in 3D modelling and printing technology, patient-specific 3D MV models could serve as a useful clinical tool. The findings also highlight the potential of this technology to be applied in a variety of medical areas within both clinical and educational settings. Copyright © 2017 Australian and New Zealand Society of Cardiac and Thoracic Surgeons (ANZSCTS) and the Cardiac Society of Australia and New Zealand (CSANZ). Published by Elsevier B.V. All rights reserved.

Semi-automatic registration of 3D orthodontics models from photographs

Science.gov (United States)

Destrez, Raphaël.; Treuillet, Sylvie; Lucas, Yves; Albouy-Kissi, Benjamin

2013-03-01

In orthodontics, a common practice used to diagnose and plan the treatment is the dental cast. After digitization by a CT-scan or a laser scanner, the obtained 3D surface models can feed orthodontics numerical tools for computer-aided diagnosis and treatment planning. One of the pre-processing critical steps is the 3D registration of dental arches to obtain the occlusion of these numerical models. For this task, we propose a vision based method to automatically compute the registration based on photos of patient mouth. From a set of matched singular points between two photos and the dental 3D models, the rigid transformation to apply to the mandible to be in contact with the maxillary may be computed by minimizing the reprojection errors. During a precedent study, we established the feasibility of this visual registration approach with a manual selection of singular points. This paper addresses the issue of automatic point detection. Based on a priori knowledge, histogram thresholding and edge detection are used to extract specific points in 2D images. Concurrently, curvatures information detects 3D corresponding points. To improve the quality of the final registration, we also introduce a combined optimization of the projection matrix with the 2D/3D point positions. These new developments are evaluated on real data by considering the reprojection errors and the deviation angles after registration in respect to the manual reference occlusion realized by a specialist.

3D Numerical Examination of Continental Mantle Lithosphere Response to Lower Crust Eclogitization and Nearby Slab Subduction

Science.gov (United States)

Janbakhsh, P.; Pysklywec, R.

2017-12-01

2D numerical modeling techniques have made great contribution to understanding geodynamic processes involved in crustal and lithospheric scale deformations for the past 20 years. The aim of this presentation is to expand the scope covered by previous researchers to 3 dimensions to address out-of-plane intrusion and extrusion of mantle material in and out of model space, and toroidal mantle wedge flows. In addition, 3D velocity boundary conditions can create more realistic models to replicate real case scenarios. 3D numerical experiments that will be presented are designed to investigate the density and viscosity effects of lower crustal eclogitization on the decoupling process of continental mantle lithosphere from the crust and its delamination. In addition, these models examine near-field effects of a subducting ocean lithosphere and a lithospheric scale fault zone on the evolution of the processes. The model solutions and predictions will also be compared against the Anatolian geology where subduction of Aegean and Arabian slabs, and the northern boundary with the North Anatolian Fault Zone are considered as two main contributing factors to anomalous crustal uplift, missing mantle lithosphere, and anomalous surface heat flux.

3D Finite Element Modelling of Drilling Process of Al2024-T3 Alloy with solid tooling and Experimental Validation

DEFF Research Database (Denmark)

Davoudinejad, Ali; Tosello, Guido

2017-01-01

Drilling is an indispensable process for many manufacturing industries due to the importance of the process for assembling components. This study presents a 3D finite element modeling (3D FEM) approach for drilling process of aluminum 2024-T3. The 3D model of tool for two facet HSSCo and four facet...... area were determined numerically. The results confirm the ability and advantage of 3D FE model of the drilling process....... HSS were generated base on the details geometry. The simulations were carried out for both drills in different cutting conditions. The numerically obtained thrust forces were compared against experimental results. The tool stress distribution, chip formation and temperature distribution in the chip...

3D modeling of high-Tc superconductors by finite element software

International Nuclear Information System (INIS)

Zhang Min; Coombs, T A

2012-01-01

A three-dimensional (3D) numerical model is proposed to solve the electromagnetic problems involving transport current and background field of a high-T c superconducting (HTS) system. The model is characterized by the E–J power law and H-formulation, and is successfully implemented using finite element software. We first discuss the model in detail, including the mesh methods, boundary conditions and computing time. To validate the 3D model, we calculate the ac loss and trapped field solution for a bulk material and compare the results with the previously verified 2D solutions and an analytical solution. We then apply our model to test some typical problems such as superconducting bulk array and twisted conductors, which cannot be tackled by the 2D models. The new 3D model could be a powerful tool for researchers and engineers to investigate problems with a greater level of complicity.

Numerical modeling of NI-monitored 3D infiltration experiment

Science.gov (United States)

Dohnal, Michal; Dusek, Jaromir; Snehota, Michal; Sacha, Jan; Vogel, Tomas; Votrubova, Jana

2014-05-01

It is well known that the temporal changes of saturated hydraulic conductivity caused by the occurrence of air phase discontinuities often play an important role in water flow and solute transport experiments. In the present study, a series of infiltration-outflow experiments was conducted to test several working hypotheses about the mechanism of air phase trapping. The experiments were performed on a porous sample with artificial internal structure, using three sandy materials with contrasting hydraulic properties. The sample was axially symmetric with continuous preferential pathways and separate porous matrix blocks (the sample was 3.4 cm in diameter and 8.8 cm high). The infiltration experiments were monitored by neutron imaging (NI). The NI data were then used to quantify the water content of the selected sample regions. The flow regime in the sample was studied using a three-dimensional model based on Richards' equation. The equation was solved by the finite element method. The results of the numerical simulations of the infiltration experiments were compared with the measured outflow rates and with the spatial distribution of water content determined by NI. The research was supported by the Czech Science Foundation Project No. 14-03691S.

3D modeling of buildings outstanding sites

CERN Document Server

Héno, Rapha?le

2014-01-01

Conventional topographic databases, obtained by capture on aerial or spatial images provide a simplified 3D modeling of our urban environment, answering the needs of numerous applications (development, risk prevention, mobility management, etc.). However, when we have to represent and analyze more complex sites (monuments, civil engineering works, archeological sites, etc.), these models no longer suffice and other acquisition and processing means have to be implemented. This book focuses on the study of adapted lifting means for "notable buildings". The methods tackled in this book cover las

Coolability of a 3D homogeneous debris bed, experimental and numerical investigations

International Nuclear Information System (INIS)

Atkhen, K.; Berthoud, G.

2001-01-01

Within the framework of nuclear safety analysis, we present here experimental and numerical results in the field of debris bed coolability. Experimental data are provided by the SILFIDE 3D experimental facility in which the debris bed is heated by induction, at Electricite de France (EDF). Numerical computations are obtained with MC3D-REPO which is a 3-phase and 3D code developed by the Commissariat a l'Energie Atomique (CEA). The uniform debris bed consists of 2 and 3,17 mm diameter steel beads contained in a 50 cm x 60 cm x 10 cm vessel. Water is used as a coolant and can be introduced either by the top or the bottom of the bed at a determined temperature. Due to heterogeneous power distribution within the bed, two definitions for the critical heat flux are proposed: the classical mean value and the local flux (much higher). Even in the first case, the measured dryout heat flux is higher than the Lipinsky 1-D flux. Temperature curve analyses show that the dryout phenomenon is very local, therefore one should be careful about the right flux definition to use. As the injected power is being increased stepwise, steady temperature stages above saturation temperature before dryout can be observed. A discussion is proposed. For some very high values of the induction power, some spheres melted together, leading to a bigger non-porous region. Even if the local temperature went over 1300 C, the bed was still coolable and the critical heat flux value was not impacted. Some parametric studies led to the following conclusions: bottom coolant injection leads to a twice time higher critical flux than by top injection, the influence of the height of the water pool above debris bed is negligible, a sub-cooled liquid injection has no influence on the coolability. Fluidization of surface particles is also discussed. The MC3D-REPO model assumes a thermal equilibrium between the three phases, which gives results in agreement with experiments until dryout occurs. (author)

3D object-oriented image analysis in 3D geophysical modelling

DEFF Research Database (Denmark)

Fadel, I.; van der Meijde, M.; Kerle, N.

2015-01-01

Non-uniqueness of satellite gravity interpretation has traditionally been reduced by using a priori information from seismic tomography models. This reduction in the non-uniqueness has been based on velocity-density conversion formulas or user interpretation of the 3D subsurface structures (objects......) based on the seismic tomography models and then forward modelling these objects. However, this form of object-based approach has been done without a standardized methodology on how to extract the subsurface structures from the 3D models. In this research, a 3D object-oriented image analysis (3D OOA......) approach was implemented to extract the 3D subsurface structures from geophysical data. The approach was applied on a 3D shear wave seismic tomography model of the central part of the East African Rift System. Subsequently, the extracted 3D objects from the tomography model were reconstructed in the 3D...

Toward a scalable flexible-order model for 3D nonlinear water waves

DEFF Research Database (Denmark)

Engsig-Karup, Allan Peter; Ducrozet, Guillaume; Bingham, Harry B.

For marine and coastal applications, current work are directed toward the development of a scalable numerical 3D model for fully nonlinear potential water waves over arbitrary depths. The model is high-order accurate, robust and efficient for large-scale problems, and support will be included...... for flexibility in the description of structures by the use of curvilinear boundary-fitted meshes. The mathematical equations for potential waves in the physical domain is transformed through $\\sigma$-mapping(s) to a time-invariant boundary-fitted domain which then becomes a basis for an efficient solution...... strategy on a time-invariant mesh. The 3D numerical model is based on a finite difference method as in the original works \\cite{LiFleming1997,BinghamZhang2007}. Full details and other aspects of an improved 3D solution can be found in \\cite{EBL08}. The new and improved approach for three...

Numerical simulation of a 3-D flow within a storage area hexagonal modular pavement systems

International Nuclear Information System (INIS)

Rashid, M A; Abustan, I; Hamzah, M O

2013-01-01

This numerical study has been performed to predict the flow patterns and characteristics within a storage area Hexagonal Modular Pavement Systems. Throughout the design and planning period for future construction are increasingly integrating computational fluid dynamics (CFD) into the process. A commercially known software, FLOW-3D, is applied to numerically solve the Navier-stokes equations for solution domains which are separated into three regions with overlapping boundaries to efficiently accommodate the grid resolutions, namely the honeycomb shaped modular, gravel and combined honeycomb shaped modular with gravel fill. The filtration of the fluid within the interstices of a permeable pavement is evaluated by integrating the Reynolds Averaged Navier-Stokes equations (RANS) inside the voids rather than making use of the widespread porous media approach. In conclusion, the results from numerical simulation are generally well agreed with the existing data and flow information such as flow patterns at increased flow, discharge rate and pressure is obtained to be used for engineering design purpose. Overall, the potential for FLOW-3D to model various geometries and configurations appears great. It should be noted that CFD should not be considered a complete replacement for physical modelling; however, it can definitely be used as a supplementary tool throughout the pavement design process.

Development of a 3D cell-centered Lagrangian scheme for the numerical modeling of the gas dynamics and hyper-elasticity systems

International Nuclear Information System (INIS)

Georges, Gabriel

2016-01-01

High Energy Density Physics (HEDP) flows are multi-material flows characterized by strong shock waves and large changes in the domain shape due to rare faction waves. Numerical schemes based on the Lagrangian formalism are good candidates to model this kind of flows since the computational grid follows the fluid motion. This provides accurate results around the shocks as well as a natural tracking of multi-material interfaces and free-surfaces. In particular, cell-centered Finite Volume Lagrangian schemes such as GLACE (Godunov-type Lagrangian scheme Conservative for total Energy) and EUCCLHYD (Explicit Unstructured Cell-Centered Lagrangian Hydrodynamics) provide good results on both the modeling of gas dynamics and elastic-plastic equations. The work produced during this PhD thesis is in continuity with the work of Maire and Nkonga [JCP, 2009] for the hydrodynamic part and the work of Kluth and Despres [JCP, 2010] for the hyper elasticity part. More precisely, the aim of this thesis is to develop robust and accurate methods for the 3D extension of the EUCCLHYD scheme with a second-order extension based on MUSCL (Monotonic Upstream-centered Scheme for Conservation Laws) and GRP (Generalized Riemann Problem) procedures. A particular care is taken on the preservation of symmetries and the monotonicity of the solutions. The scheme robustness and accuracy are assessed on numerous Lagrangian test cases for which the 3D extensions are very challenging. (author) [fr

Large-scale thermal convection of viscous fluids in a faulted system: 3D test case for numerical codes

Science.gov (United States)

Magri, Fabien; Cacace, Mauro; Fischer, Thomas; Kolditz, Olaf; Wang, Wenqing; Watanabe, Norihiro

2017-04-01

In contrast to simple homogeneous 1D and 2D systems, no appropriate analytical solutions exist to test onset of thermal convection against numerical models of complex 3D systems that account for variable fluid density and viscosity as well as permeability heterogeneity (e.g. presence of faults). Owing to the importance of thermal convection for the transport of energy and minerals, the development of a benchmark test for density/viscosity driven flow is crucial to ensure that the applied numerical models accurately simulate the physical processes at hands. The presented study proposes a 3D test case for the simulation of thermal convection in a faulted system that accounts for temperature dependent fluid density and viscosity. The linear stability analysis recently developed by Malkovsky and Magri (2016) is used to estimate the critical Rayleigh number above which thermal convection of viscous fluids is triggered. The numerical simulations are carried out using the finite element technique. OpenGeoSys (Kolditz et al., 2012) and Moose (Gaston et al., 2009) results are compared to those obtained using the commercial software FEFLOW (Diersch, 2014) to test the ability of widely applied codes in matching both the critical Rayleigh number and the dynamical features of convective processes. The methodology and Rayleigh expressions given in this study can be applied to any numerical model that deals with 3D geothermal processes in faulted basins as by example the Tiberas Basin (Magri et al., 2016). References Kolditz, O., Bauer, S., Bilke, L., Böttcher, N., Delfs, J. O., Fischer, T., U. J. Görke, T. Kalbacher, G. Kosakowski, McDermott, C. I., Park, C. H., Radu, F., Rink, K., Shao, H., Shao, H.B., Sun, F., Sun, Y., Sun, A., Singh, K., Taron, J., Walther, M., Wang,W., Watanabe, N., Wu, Y., Xie, M., Xu, W., Zehner, B., 2012. OpenGeoSys: an open-source initiative for numerical simulation of thermo-hydro-mechanical/chemical (THM/C) processes in porous media. Environmental

Comparing a quasi-3D to a full 3D nearshore circulation model: SHORECIRC and ROMS

Science.gov (United States)

Haas, Kevin A.; Warner, John C.

2009-01-01

Predictions of nearshore and surf zone processes are important for determining coastal circulation, impacts of storms, navigation, and recreational safety. Numerical modeling of these systems facilitates advancements in our understanding of coastal changes and can provide predictive capabilities for resource managers. There exists many nearshore coastal circulation models, however they are mostly limited or typically only applied as depth integrated models. SHORECIRC is an established surf zone circulation model that is quasi-3D to allow the effect of the variability in the vertical structure of the currents while maintaining the computational advantage of a 2DH model. Here we compare SHORECIRC to ROMS, a fully 3D ocean circulation model which now includes a three dimensional formulation for the wave-driven flows. We compare the models with three different test applications for: (i) spectral waves approaching a plane beach with an oblique angle of incidence; (ii) monochromatic waves driving longshore currents in a laboratory basin; and (iii) monochromatic waves on a barred beach with rip channels in a laboratory basin. Results identify that the models are very similar for the depth integrated flows and qualitatively consistent for the vertically varying components. The differences are primarily the result of the vertically varying radiation stress utilized by ROMS and the utilization of long wave theory for the radiation stress formulation in vertical varying momentum balance by SHORECIRC. The quasi-3D model is faster, however the applicability of the fully 3D model allows it to extend over a broader range of processes, temporal, and spatial scales.

3D for Geosciences: Interactive Tangibles and Virtual Models

Science.gov (United States)

Pippin, J. E.; Matheney, M.; Kitsch, N.; Rosado, G.; Thompson, Z.; Pierce, S. A.

2016-12-01

Point cloud processing provides a method of studying and modelling geologic features relevant to geoscience systems and processes. Here, software including Skanect, MeshLab, Blender, PDAL, and PCL are used in conjunction with 3D scanning hardware, including a Structure scanner and a Kinect camera, to create and analyze point cloud images of small scale topography, karst features, tunnels, and structures at high resolution. This project successfully scanned internal karst features ranging from small stalactites to large rooms, as well as an external waterfall feature. For comparison purposes, multiple scans of the same object were merged into single object files both automatically, using commercial software, and manually using open source libraries and code. Files with format .ply were manually converted into numeric data sets to be analyzed for similar regions between files in order to match them together. We can assume a numeric process would be more powerful and efficient than the manual method, however it could lack other useful features that GUI's may have. The digital models have applications in mining as efficient means of replacing topography functions such as measuring distances and areas. Additionally, it is possible to make simulation models such as drilling templates and calculations related to 3D spaces. Advantages of using methods described here for these procedures include the relatively quick time to obtain data and the easy transport of the equipment. With regard to openpit mining, obtaining 3D images of large surfaces and with precision would be a high value tool by georeferencing scan data to interactive maps. The digital 3D images obtained from scans may be saved as printable files to create physical 3D-printable models to create tangible objects based on scientific information, as well as digital "worlds" able to be navigated virtually. The data, models, and algorithms explored here can be used to convey complex scientific ideas to a range of

Model-based meta-analysis for comparing Vitamin D2 and D3 parent-metabolite pharmacokinetics.

Science.gov (United States)

Ocampo-Pelland, Alanna S; Gastonguay, Marc R; Riggs, Matthew M

2017-08-01

Association of Vitamin D (D3 & D2) and its 25OHD metabolite (25OHD3 & 25OHD2) exposures with various diseases is an active research area. D3 and D2 dose-equivalency and each form's ability to raise 25OHD concentrations are not well-defined. The current work describes a population pharmacokinetic (PK) model for D2 and 25OHD2 and the use of a previously developed D3-25OHD3 PK model [1] for comparing D3 and D2-related exposures. Public-source D2 and 25OHD2 PK data in healthy or osteoporotic populations, including 17 studies representing 278 individuals (15 individual-level and 18 arm-level units), were selected using search criteria in PUBMED. Data included oral, single and multiple D2 doses (400-100,000 IU/d). Nonlinear mixed effects models were developed simultaneously for D2 and 25OHD2 PK (NONMEM v7.2) by considering 1- and 2-compartment models with linear or nonlinear clearance. Unit-level random effects and residual errors were weighted by arm sample size. Model simulations compared 25OHD exposures, following repeated D2 and D3 oral administration across typical dosing and baseline ranges. D2 parent and metabolite were each described by 2-compartment models with numerous parameter estimates shared with the D3-25OHD3 model [1]. Notably, parent D2 was eliminated (converted to 25OHD) through a first-order clearance whereas the previously published D3 model [1] included a saturable non-linear clearance. Similar to 25OHD3 PK model results [1], 25OHD2 was eliminated by a first-order clearance, which was almost twice as fast as the former. Simulations at lower baselines, following lower equivalent doses, indicated that D3 was more effective than D2 at raising 25OHD concentrations. Due to saturation of D3 clearance, however, at higher doses or baselines, the probability of D2 surpassing D3's ability to raise 25OHD concentrations increased substantially. Since 25OHD concentrations generally surpassed 75 nmol/L at these higher baselines by 3 months, there would be no

3D engineered models for highway construction : the Iowa experience.

Science.gov (United States)

2015-06-01

3D engineered modeling is a relatively new and developing technology that can provide numerous bene ts to owners, engineers, : contractors, and the general public. This manual is for highway agencies that are considering or are in the process of s...

Implementation of 3D models in the Monte Carlo code MCNP

International Nuclear Information System (INIS)

Lopes, Vivaldo; Millian, Felix M.; Guevara, Maria Victoria M.; Garcia, Fermin; Sena, Isaac; Menezes, Hugo

2009-01-01

On the area of numerical dosimetry Applied to medical physics, the scientific community focuses on the elaboration of new hybrids models based on 3D models. But different steps of the process of simulation with 3D models needed improvement and optimization in order to expedite the calculations and accuracy using this methodology. This project was developed with the aim of optimize the process of introduction of 3D models within the simulation code of radiation transport by Monte Carlo (MCNP). The fast implementation of these models on the simulation code allows the estimation of the dose deposited on the patient organs on a more personalized way, increasing the accuracy with this on the estimates and reducing the risks to health, caused by ionizing radiations. The introduction o these models within the MCNP was made through a input file, that was constructed through a sequence of images, bi-dimensional in the 3D model, generated using the program '3DSMAX', imported by the program 'TOMO M C' and thus, introduced as INPUT FILE of the MCNP code. (author)

A 3D heat conduction model for block-type high temperature reactors and its implementation into the code DYN3D

International Nuclear Information System (INIS)

Baier, Silvio; Kliem, Soeren; Rohde, Ulrich

2011-01-01

The gas-cooled high temperature reactor is a concept to produce energy at high temperatures with a high level of inherent safety. It gets special attraction due to e.g. high thermal efficiency and the possibility of hydrogen production. In addition to the PBMR (Pebble Bed Modular Reactor) the (V)HTR (Very high temperature reactor) concept has been established. The basic design of a prismatic HTR consists of the following elements. The fuel is coated with four layers of isotropic materials. These so-called TRISO particles are dispersed into compacts which are placed in a graphite block matrix. The graphite matrix additionally contains holes for the coolant gas. A one-dimensional model is sufficient to describe (the radial) heat transfer in LWRs. But temperature gradients in a prismatic HTR can occur in axial as well as in radial direction, since regions with different heat source release and with different coolant temperature heat up are coupled through the graphite matrix elements. Furthermore heat transfer into reflector elements is possible. DYN3D is a code system for coupled neutron and thermal hydraulics core calculations developed at the Helmholtzzentrum Dresden-Rossendorf. Concerning neutronics DYN3D consists of a two-group and multi-group diffusion approach based on nodal expansion methods. Furthermore a 1D thermal-hydraulics model for parallel coolant flow channels is included. The DYN3D code was extensively verified and validated via numerous numerical and experimental benchmark problems. That includes the NEA CRP benchmarks for PWR and BWR, the Three-Miles-Island-1 main steam line break and the Peach Bottom Turbine Trip benchmarks, as well as measurements carried out in an original-size VVER-1000 mock-up. An overview of the verification and validation activities can be found. Presently a DYN3D-HTR version is under development. It involves a 3D heat conduction model to deal with higher-(than one)-dimensional effects of heat transfer and heat conduction in

Radon transport modelling: User's guide to RnMod3d

International Nuclear Information System (INIS)

Andersen, C.E.

2000-08-01

RnMod3d is a numerical computer model of soil-gas and radon transport in porous media. It can be used, for example, to study radon entry from soil into houses in response to indoor-outdoor pressure differences or changes in atmospheric pressure. It can also be used for flux calculations of radon from the soil surface or to model radon exhalation from building materials such as concrete. The finite-volume model is a technical research tool, and it cannot be used meaningfully without good understanding of the involved physical equations. Some understanding of numerical mathematics and the programming language Pascal is also required. Originally, the code was developed for internal use at Risoe only. With this guide, however, it should be possible for others to use the model. Three-dimensional steady-state or transient problems with Darcy flow of soil gas and combined generation, radioactive decay, diffusion and advection of radon can be solved. Moisture is included in the model, and partitioning of radon between air, water and soil grains (adsorption) is taken into account. Most parameters can change in time and space, and transport parameters (diffusivity and permeability) may be anisotropic. This guide includes benchmark tests based on simple problems with known solutions. RnMod3d has also been part of an international model intercomparison exercise based on more complicated problems without known solutions. All tests show that RnMod3d gives results of good quality. (au)

3D thermal modeling of TRISO fuel coupled with neutronic simulation

Energy Technology Data Exchange (ETDEWEB)

Hu, Jianwei [Los Alamos National Laboratory; Uddin, Rizwan [UNIV OF ILLINIOS

2010-01-01

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

Development of numerical solution techniques in the KIKO3D code

International Nuclear Information System (INIS)

Panka, Istvan; Kereszturi, Andras; Hegedus, Csaba

2005-01-01

The paper describes the numerical methods applied in KIKO3D three-dimensional reactor dynamics code and present a new, more effective method (Bi-CGSTAB) for accelerating the large sparse matrix equation solution. The convergence characteristics were investigated in a given macro time step of a Control Rod Ejection transient. The results obtained by the old GMRES and new Bi-CGSTAB methods are compared. It is concluded that the real relative errors of the solutions obtained by GMRES or Bi - CGSTAB algorithms are in fact closer together than the estimated relative errors. The KIKO3D-Bi-CGSTAB method converges safely and it is 7-12 % faster than the old KIKO3D-GMRES solution (Authors)

Single-sided sheet-to-tube spot welding investigated by 3D numerical simulations

DEFF Research Database (Denmark)

Nielsen, Chris Valentin; Chergui, Azeddine; Zhang, Wenqi

The single-sided resistance spot welding process is analyzed by a 3D numerical study of sheet-to-tube joining. Finite element simulations are carried out in SORPAS® 3D. Two levels of electrode force and five levels of welding current are simulated. The overall effects of changing current and force...

Development of Numerical Grids for UZ Flow and Transport Modeling

International Nuclear Information System (INIS)

Hinds, J.

2001-01-01

This Analysis/Model Report (AMR) describes the methods used to develop numerical grids of the unsaturated hydrogeologic system beneath Yucca Mountain. Numerical grid generation is an integral part of the development of a complex, three-dimensional (3-D) model, such as the Unsaturated-Zone Flow and Transport Model (UZ Model) of Yucca Mountain. The resulting numerical grids, developed using current geologic, hydrogeologic, and mineralogic data, provide the necessary framework to: (1) develop calibrated hydrogeologic property sets and flow fields, (2) test conceptual hypotheses of flow and transport, and (3) predict flow and transport behavior under a variety of climatic and thermal loading conditions. Revision 00 of the work described herein follows the planning and work direction outlined in the ''Development of Numerical Grids for UZ Flow and Transport Modeling'' (CRWMS M and O 1999c). The technical scope, content, and management of ICN 01 of this AMR is currently controlled by the planning document, ''Technical Work Plan for Unsaturated Zone (UZ) Flow and Transport Process Model Report'' (BSC 2001a). The scope for the TBV resolution actions in this ICN is described in the ''Technical Work Plan for: Integrated Management of Technical Product Input Department'' (BSC 2001 b, Addendum B, Section 4.1). The steps involved in numerical grid development include: (1) defining the location of important calibration features, (2) determining model grid layers and fault geometry based on the Geologic Framework Model (GFM), the Integrated Site Model (ISM), and definition of hydrogeologic units (HGUs), (3) analyzing and extracting GFM and ISM data pertaining to layer contacts and property distributions, (4) discretizing and refining the two-dimensional (2-D), plan-view numerical grid, (5) generating the 3-D grid with finer resolution at the repository horizon and within the Calico Hills nonwelded (CHn) hydrogeologic unit, and (6) formulating the dual-permeability mesh. The

Integration of 3D photogrammetric outcrop models in the reservoir modelling workflow

Science.gov (United States)

Deschamps, Remy; Joseph, Philippe; Lerat, Olivier; Schmitz, Julien; Doligez, Brigitte; Jardin, Anne

2014-05-01

3D technologies are now widely used in geosciences to reconstruct outcrops in 3D. The technology used for the 3D reconstruction is usually based on Lidar, which provides very precise models. Such datasets offer the possibility to build well-constrained outcrop analogue models for reservoir study purposes. The photogrammetry is an alternate methodology which principles are based in determining the geometric properties of an object from photographic pictures taken from different angles. Outcrop data acquisition is easy, and this methodology allows constructing 3D outcrop models with many advantages such as: - light and fast acquisition, - moderate processing time (depending on the size of the area of interest), - integration of field data and 3D outcrops into the reservoir modelling tools. Whatever the method, the advantages of digital outcrop model are numerous as already highlighted by Hodgetts (2013), McCaffrey et al. (2005) and Pringle et al. (2006): collection of data from otherwise inaccessible areas, access to different angles of view, increase of the possible measurements, attributes analysis, fast rate of data collection, and of course training and communication. This paper proposes a workflow where 3D geocellular models are built by integrating all sources of information from outcrops (surface picking, sedimentological sections, structural and sedimentary dips…). The 3D geomodels that are reconstructed can be used at the reservoir scale, in order to compare the outcrop information with subsurface models: the detailed facies models of the outcrops are transferred into petrophysical and acoustic models, which are used to test different scenarios of seismic and fluid flow modelling. The detailed 3D models are also used to test new techniques of static reservoir modelling, based either on geostatistical approaches or on deterministic (process-based) simulation techniques. A modelling workflow has been designed to model reservoir geometries and properties from

Radon transport modelling: User's guide to RnMod3d

DEFF Research Database (Denmark)

Andersen, Claus Erik

2000-01-01

RnMod3d is a numerical computer model of soil-gas and radon transport in porous media. It can be used, for example, to study radon entry from soil into houses in response to indoor-outdoor pressure differences or changes in atmospheric pressure. It canalso be used for flux calculations of radon...... decay, diffusion and advection of radon can be solved. Moisture is included in the model, and partitioning ofradon between air, water and soil grains (adsorption) is taken into account. Most parameters can change in time and space, and transport parameters (diffusivity and permeability) may...... be anisotropic. This guide includes benchmark tests based on simpleproblems with known solutions. RnMod3d has also been part of an international model intercomparison exercise based on more complicated problems without known solutions. All tests show that RnMod3d gives results of good quality....

3D numerical modeling of YSO accretion shocks

Directory of Open Access Journals (Sweden)

Matsakos T.

2014-01-01

Full Text Available The dynamics of YSO accretion shocks is determined by radiative processes as well as the strength and structure of the magnetic field. A quasi-periodic emission signature is theoretically expected to be observed, but observations do not confirm any such pattern. In this work, we assume a uniform background field, in the regime of optically thin energy losses, and we study the multi-dimensional shock evolution in the presence of perturbations, i.e. clumps in the stream and an acoustic energy flux flowing at the base of the chromosphere. We perform 3D MHD simulations using the PLUTO code, modelling locally the impact of the infalling gas onto the chromosphere. We find that the structure and dynamics of the post-shock region is strongly dependent on the plasma-beta (thermal over magnetic pressure, different values of which may give distinguishable emission signatures, relevant for observations. In particular, a strong magnetic field effectively confines the plasma inside its flux tubes and leads to the formation of quasi-independent fibrils. The fibrils may oscillate out of phase and hence the sum of their contributions in the emission results in a smooth overall profile. On the contrary, a weak magnetic field is not found to have any significant effect on the shocked plasma and the turbulent hot slab that forms is found to retain its periodic signature.

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

Science.gov (United States)

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 Modelling of the Czochralski Growth of β-Ga2O3

Directory of Open Access Journals (Sweden)

Wolfram Miller

2017-01-01

Full Text Available Our numerical modelling of the Czochralski growth of single crystalline β-Ga 2 O 3 crystals (monoclinic symmetry starts at the 2D heat transport analysis within the crystal growth furnace, proceeds with the 3D heat transport and fluid flow analysis in the crystal-melt-crucible arrangement and targets the 3D thermal stress analysis within the β-Ga 2 O 3 crystal. In order to perform the stress analysis, we measured the thermal expansion coefficients and the elastic stiffness coefficients in two samples of a β-Ga 2 O 3 crystal grown at IKZ. Additionally, we analyse published data of β-Ga 2 O 3 material properties and use data from literature for comparative calculations. The computations were performed by the software packages CrysMAS, CGsim, Ansys-cfx and comsol Multiphysics. By the hand of two different thermal expansion data sets and two different crystal orientations, we analyse the elastic stresses in terms of the von-Mises stress.

3D Digital Modelling

DEFF Research Database (Denmark)

Hundebøl, Jesper

wave of new building information modelling tools demands further investigation, not least because of industry representatives' somewhat coarse parlance: Now the word is spreading -3D digital modelling is nothing less than a revolution, a shift of paradigm, a new alphabet... Research qeustions. Based...... on empirical probes (interviews, observations, written inscriptions) within the Danish construction industry this paper explores the organizational and managerial dynamics of 3D Digital Modelling. The paper intends to - Illustrate how the network of (non-)human actors engaged in the promotion (and arrest) of 3...... important to appreciate the analysis. Before turning to the presentation of preliminary findings and a discussion of 3D digital modelling, it begins, however, with an outline of industry specific ICT strategic issues. Paper type. Multi-site field study...

Mechanical Behaviour of 3D Multi-layer Braided Composites: Experimental, Numerical and Theoretical Study

Science.gov (United States)

Deng, Jian; Zhou, Guangming; Ji, Le; Wang, Xiaopei

2017-12-01

Mechanical properties and failure mechanisms of a newly designed 3D multi-layer braided composites are evaluated by experimental, numerical and theoretical studies. The microstructure of the composites is introduced. The unit cell technique is employed to address the periodic arrangement of the structure. The volume averaging method is used in theoretical solutions while FEM with reasonable periodic boundary conditions and meshing technique in numerical simulations. Experimental studies are also conducted to verify the feasibility of the proposed models. Predicted elastic properties agree well with the experimental data, indicating the feasibility of the proposed models. Numerical evaluation is more accurate than theoretical assessment. Deformations and stress distributions of the unit cell under tension shows displacement and traction continuity, guaranteeing the rationality of the applied periodic boundary conditions. Although compression and tension modulus are close, the compressive strength only reaches 70% of the tension strength. This indicates that the composites can be weakened in compressive loading. Additionally, by analysing the micrograph of fracture faces and strain-stress curves, a brittle failure mechanism is observed both in composites under tension and compression.

Calculation of Tsunami Damage and preparation of Inundation Maps by 2D and 3D numerical modeling in Göcek, Turkey

Science.gov (United States)

Ozer Sozdinler, C.; Arikawa, T.; Necmioglu, O.; Ozel, N. M.

2016-12-01

The Aegean and its surroundings form the most active part of the Africa-Eurasia collision zone responsible for the high level of seismicity in this region. It constitutes more than 60% of the expected seismicity in Europe up to Mw=8.2 (Moratto et al., 2007; Papazachos, 1990). Shaw and Jackson (2010) argued that the existing system of Hellenic Arc subduction-zone is capable of allowing very large but rare earthquakes on splay faults, such as the one occurred in 365, together with the contribution of small earthquakes. Based on an extensive earthquake generated tsunami scenario database, Necmioğlu and Özel (2015) showed that maximum wave heights in the Eastern Mediterranean for shallow earthquakes defined is >3 m in locations in, around and orthogonal to the Hellenic Arc. Considering the seismicity and the tsunami potential in Eastern Mediterranean, the investigation and monitoring of earthquake and tsunami hazard, and the preparation of mitigation strategies and national resilience plans would become inevitable in Turkey. Gocek town, as one of the Tsunami Forecast Points having a unique geography with many small bays and islands and a very popular touristic destination especially for yachtsmen, is selected in this study for the tsunami modeling by using high resolution bathymetric and topographic data with less than 4m grid size. The tsunami analyses are performed by the numerical codes NAMIDANCE (NAMIDANCE,2011) for 2D modeling and STOC-CADMAS (Arikawa,2014) for 3D modeling for the calculations of tsunami hydrodynamic parameters. Froude numbers, as one of the most important indicators for tsunami damage (Ozer, 2012) and the directions of current velocities inside marinas are solved by NAMIDANCE while STOC-CADMAS determines the tsunami pressure and force exerted onto the sea and land structures with 3D and non-hydrostatic approaches. The results are then used to determine the tsunami inundation and structural resilience and establish the tsunami preparedness and

Spectral-element Method for 3D Marine Controlled-source EM Modeling

Science.gov (United States)

Liu, L.; Yin, C.; Zhang, B., Sr.; Liu, Y.; Qiu, C.; Huang, X.; Zhu, J.

2017-12-01

As one of the predrill reservoir appraisal methods, marine controlled-source EM (MCSEM) has been widely used in mapping oil reservoirs to reduce risk of deep water exploration. With the technical development of MCSEM, the need for improved forward modeling tools has become evident. We introduce in this paper spectral element method (SEM) for 3D MCSEM modeling. It combines the flexibility of finite-element and high accuracy of spectral method. We use Galerkin weighted residual method to discretize the vector Helmholtz equation, where the curl-conforming Gauss-Lobatto-Chebyshev (GLC) polynomials are chosen as vector basis functions. As a kind of high-order complete orthogonal polynomials, the GLC have the characteristic of exponential convergence. This helps derive the matrix elements analytically and improves the modeling accuracy. Numerical 1D models using SEM with different orders show that SEM method delivers accurate results. With increasing SEM orders, the modeling accuracy improves largely. Further we compare our SEM with finite-difference (FD) method for a 3D reservoir model (Figure 1). The results show that SEM method is more effective than FD method. Only when the mesh is fine enough, can FD achieve the same accuracy of SEM. Therefore, to obtain the same precision, SEM greatly reduces the degrees of freedom and cost. Numerical experiments with different models (not shown here) demonstrate that SEM is an efficient and effective tool for MSCEM modeling that has significant advantages over traditional numerical methods.This research is supported by Key Program of National Natural Science Foundation of China (41530320), China Natural Science Foundation for Young Scientists (41404093), and Key National Research Project of China (2016YFC0303100, 2017YFC0601900).

Constraining the Timescales of Rehydration in Nominally Anhydrous Minerals Using 3D Numerical Diffusion Models

Science.gov (United States)

Lynn, K. J.; Warren, J. M.

2017-12-01

Nominally anhydrous minerals (NAMs) are important for characterizing deep-Earth water reservoirs, but the water contents of olivine (ol), orthopyroxene (opx), and clinopyroxene (cpx) in peridotites generally do not reflect mantle equilibrium conditions. Ol is typically "dry" and decoupled from H in cpx and opx, which is inconsistent with models of partial melting and/or diffusive loss of H during upwelling beneath mid-ocean ridges. The rehydration of mantle pyroxenes via late-stage re-fertilization has been invoked to explain their relatively high water contents. Here, we use sophisticated 3D diffusion models (after Shea et al., 2015, Am Min) of H in ol, opx, and cpx to investigate the timescales of rehydration across a range of conditions relevant for melt-rock interaction and serpentinization of peridotites. Numerical crystals with 1 mm c-axis lengths and realistic crystal morphologies are modeled using recent H diffusivities that account for compositional variation and diffusion anisotropy. Models were run over timescales of minutes to millions of years and temperatures from 300 to 1200°C. Our 3D models show that, at the high-T end of the range, H concentrations in the cores of NAMs are partially re-equilibrated in as little as a few minutes, and completely re-equilibrated within hours to weeks. At low-T (300°C), serpentinization can induce considerable diffusion in cpx and opx. H contents are 30% re-equilibrated after continuous exposure to hydrothermal fluids for 102 and 105 years, respectively, which is inconsistent with previous interpretations that there is no effect on H in opx under similar conditions. Ol is unaffected after 1 Myr due to the slower diffusivity of the proton-vacancy mechanism at 300°C (2-4 log units lower than for opx). In the middle of the T range (700-1000°C), rehydration of opx and cpx occurs over hours to days, while ol is somewhat slower to respond (days to weeks), potentially allowing the decoupling observed in natural samples to

Surface Simplification of 3D Animation Models Using Robust Homogeneous Coordinate Transformation

Directory of Open Access Journals (Sweden)

Juin-Ling Tseng

2014-01-01

Full Text Available The goal of 3D surface simplification is to reduce the storage cost of 3D models. A 3D animation model typically consists of several 3D models. Therefore, to ensure that animation models are realistic, numerous triangles are often required. However, animation models that have a high storage cost have a substantial computational cost. Hence, surface simplification methods are adopted to reduce the number of triangles and computational cost of 3D models. Quadric error metrics (QEM has recently been identified as one of the most effective methods for simplifying static models. To simplify animation models by using QEM, Mohr and Gleicher summed the QEM of all frames. However, homogeneous coordinate problems cannot be considered completely by using QEM. To resolve this problem, this paper proposes a robust homogeneous coordinate transformation that improves the animation simplification method proposed by Mohr and Gleicher. In this study, the root mean square errors of the proposed method were compared with those of the method proposed by Mohr and Gleicher, and the experimental results indicated that the proposed approach can preserve more contour features than Mohr’s method can at the same simplification ratio.

Geodynamical Evolution of the En echelon Basins in the Hexi Corridor: Implications From 3-D Numerical Modeling

Science.gov (United States)

Li, W.; Shi, Y.; Zhang, H.; Cheng, H.

2017-12-01

The Hexi Corridor, located between the Alax block and the Caledon fold belt in the North Qilian Mountains, is the forefront area of northward thrust of the Tibet Plateau. Most notably, this active tectonic region consists of a series of faults and western-northwest trending Cenozoic basins. Therefore, it's a pivotal part in terms of recording tectonic pattern of the Tibet Plateau and also demonstrating the northward growth of Tibetan Plateau. In order to explain the mechanism of formation and evolution of the paired basins in the Hexi Corridor and based on the visco-elasticity-plasticity constitutive relation, we construct a 3-D finite element numerical model, including the Altun Tagh fault zone, the northern Qilian Shan-Hexi corridor faults system and the Haiyuan fault zone in northeast of the Tibet Plateau.The boundary conditions are constrained by GPS observations and fault slip rate provided by field geology, with steady rate of deformation of north-south compression and lateral shear along the approximately east-west strike fault zones.In our numerical model, different blocks are given different mechanical features and major fault zones are assumed mechanical weak zones. The long-term (5Ma) accumulation of lithospheric stress, displacement and fault dislocation of the Hexi Corridor and its adjacent regions are calculated in different models for comparison. Meanwhile, we analyze analyzed how the crustal heterogeneity affecting the tectonic deformations in this region. Comparisons between the numerical results and the geological observations indicate that under compression-shear boundary conditions, heterogeneous blocks of various scales may lead to the development of en echelon faults and basins in the Hexi corridor. And the ectonic deformation of Alax and the North Qilian Mountains are almost simultaneous, which may be earlier than the initiation of en echelon basins in the Hexi Corridor and the faults between the en echelon basins. Calculated horizontal and

Strategy for a numerical Rock Mechanics Site Descriptive Model. Further development of the theoretical/numerical approach

International Nuclear Information System (INIS)

Olofsson, Isabelle; Fredriksson, Anders

2005-05-01

The Swedish Nuclear and Fuel Management Company (SKB) is conducting Preliminary Site Investigations at two different locations in Sweden in order to study the possibility of a Deep Repository for spent fuel. In the frame of these Site Investigations, Site Descriptive Models are achieved. These products are the result of an interaction of several disciplines such as geology, hydrogeology, and meteorology. The Rock Mechanics Site Descriptive Model constitutes one of these models. Before the start of the Site Investigations a numerical method using Discrete Fracture Network (DFN) models and the 2D numerical software UDEC was developed. Numerical simulations were the tool chosen for applying the theoretical approach for characterising the mechanical rock mass properties. Some shortcomings were identified when developing the methodology. Their impacts on the modelling (in term of time and quality assurance of results) were estimated to be so important that the improvement of the methodology with another numerical tool was investigated. The theoretical approach is still based on DFN models but the numerical software used is 3DEC. The main assets of the programme compared to UDEC are an optimised algorithm for the generation of fractures in the model and for the assignment of mechanical fracture properties. Due to some numerical constraints the test conditions were set-up in order to simulate 2D plane strain tests. Numerical simulations were conducted on the same data set as used previously for the UDEC modelling in order to estimate and validate the results from the new methodology. A real 3D simulation was also conducted in order to assess the effect of the '2D' conditions in the 3DEC model. Based on the quality of the results it was decided to update the theoretical model and introduce the new methodology based on DFN models and 3DEC simulations for the establishment of the Rock Mechanics Site Descriptive Model. By separating the spatial variability into two parts, one

Computational and experimental analysis of supersonic air ejector: Turbulence modeling and assessment of 3D effects

International Nuclear Information System (INIS)

Mazzelli, Federico; Little, Adrienne B.; Garimella, Srinivas; Bartosiewicz, Yann

2015-01-01

Highlights: • Computational and experimental assessment of computational techniques for ejector flows. • Comparisons to 2D/3D (k–ε, k–ε realizable, k–ω SST, and stress–ω RSM) turbulence models. • k–ω SST model performs best while ε-based models more accurate at low motive pressures. • Good on-design agreement across 2D and 3D models; off-design needs 3D simulations. - Abstract: Numerical and experimental analyses are performed on a supersonic air ejector to evaluate the effectiveness of commonly-used computational techniques when predicting ejector flow characteristics. Three series of experimental curves at different operating conditions are compared with 2D and 3D simulations using RANS, steady, wall-resolved models. Four different turbulence models are tested: k–ε, k–ε realizable, k–ω SST, and the stress–ω Reynolds Stress Model. An extensive analysis is performed to interpret the differences between numerical and experimental results. The results show that while differences between turbulence models are typically small with respect to the prediction of global parameters such as ejector inlet mass flow rates and Mass Entrainment Ratio (MER), the k–ω SST model generally performs best whereas ε-based models are more accurate at low motive pressures. Good agreement is found across all 2D and 3D models at on-design conditions. However, prediction at off-design conditions is only acceptable with 3D models, making 3D simulations mandatory to correctly predict the critical pressure and achieve reasonable results at off-design conditions. This may partly depend on the specific geometry under consideration, which in the present study has a rectangular cross section with low aspect ratio.

AxiSEM3D: a new fast method for global wave propagation in 3-D Earth models with undulating discontinuities

Science.gov (United States)

Leng, K.; Nissen-Meyer, T.; van Driel, M.; Al-Attar, D.

2016-12-01

We present a new, computationally efficient numerical method to simulate global seismic wave propagation in realistic 3-D Earth models with laterally heterogeneous media and finite boundary perturbations. Our method is a hybrid of pseudo-spectral and spectral element methods (SEM). We characterize the azimuthal dependence of 3-D wavefields in terms of Fourier series, such that the 3-D equations of motion reduce to an algebraic system of coupled 2-D meridional equations, which can be solved by a 2-D spectral element method (based on www.axisem.info). Computational efficiency of our method stems from lateral smoothness of global Earth models (with respect to wavelength) as well as axial singularity of seismic point sources, which jointly confine the Fourier modes of wavefields to a few lower orders. All boundary perturbations that violate geometric spherical symmetry, including Earth's ellipticity, topography and bathymetry, undulations of internal discontinuities such as Moho and CMB, are uniformly considered by means of a Particle Relabeling Transformation.The MPI-based high performance C++ code AxiSEM3D, is now available for forward simulations upon 3-D Earth models with fluid outer core, ellipticity, and both mantle and crustal structures. We show novel benchmarks for global wave solutions in 3-D mantle structures between our method and an independent, fully discretized 3-D SEM with remarkable agreement. Performance comparisons are carried out on three state-of-the-art tomography models, with seismic period going down to 5s. It is shown that our method runs up to two orders of magnitude faster than the 3-D SEM for such settings, and such computational advantage scales favourably with seismic frequency. By examining wavefields passing through hypothetical Gaussian plumes of varying sharpness, we identify in model-wavelength space the limits where our method may lose its advantage.

3-D numerical investigation of subsurface flow in anisotropic porous media using multipoint flux approximation method

KAUST Repository

Negara, Ardiansyah

2013-01-01

Anisotropy of hydraulic properties of subsurface geologic formations is an essential feature that has been established as a consequence of the different geologic processes that they undergo during the longer geologic time scale. With respect to petroleum reservoirs, in many cases, anisotropy plays significant role in dictating the direction of flow that becomes no longer dependent only on the pressure gradient direction but also on the principal directions of anisotropy. Furthermore, in complex systems involving the flow of multiphase fluids in which the gravity and the capillarity play an important role, anisotropy can also have important influences. Therefore, there has been great deal of motivation to consider anisotropy when solving the governing conservation laws numerically. Unfortunately, the two-point flux approximation of finite difference approach is not capable of handling full tensor permeability fields. Lately, however, it has been possible to adapt the multipoint flux approximation that can handle anisotropy to the framework of finite difference schemes. In multipoint flux approximation method, the stencil of approximation is more involved, i.e., it requires the involvement of 9-point stencil for the 2-D model and 27-point stencil for the 3-D model. This is apparently challenging and cumbersome when making the global system of equations. In this work, we apply the equation-type approach, which is the experimenting pressure field approach that enables the solution of the global problem breaks into the solution of multitude of local problems that significantly reduce the complexity without affecting the accuracy of numerical solution. This approach also leads in reducing the computational cost during the simulation. We have applied this technique to a variety of anisotropy scenarios of 3-D subsurface flow problems and the numerical results demonstrate that the experimenting pressure field technique fits very well with the multipoint flux approximation

Simulation Of Seawater Intrusion With 2D And 3D Models: Nauru Island Case Study

Science.gov (United States)

Ghassemi, F.; Jakeman, A. J.; Jacobson, G.; Howard, K. W. F.

1996-03-01

With the advent of large computing capacities during the past few decades, sophisticated models have been developed for the simulation of seawater intrusion in coastal and island aquifers. Currently, several models are commercially available for the simulation of this problem. This paper describes the mathematical basis and application of the SUTRA and HST3D models to simulate seawater intrusion in Nauru Island, in the central Pacific Ocean. A comparison of the performance and limitations of these two models in simulating a real problem indicates that three-dimensional simulation of seawater intrusion with the HST3D model has the major advantage of being able to specify natural boundary conditions as well as pumping stresses. However, HST3D requires a small grid size and short time steps in order to maintain numerical stability and accuracy. These requirements lead to solution of a large set of linear equations that requires the availability of powerful computing facilities in terms of memory and computing speed. Combined results of the two simulation models indicate a safe pumping rate of 400 m3/d for the aquifer on Nauru Island, where additional fresh water is presently needed for the rehabilitation of mined-out land.

Numerical Modeling of 3D Seismic Wave Propagation around Yogyakarta, the Southern Part of Central Java, Indonesia, Using Spectral-Element Method on MPI-GPU Cluster

Science.gov (United States)

Sudarmaji; Rudianto, Indra; Eka Nurcahya, Budi

2018-04-01

A strong tectonic earthquake with a magnitude of 5.9 Richter scale has been occurred in Yogyakarta and Central Java on May 26, 2006. The earthquake has caused severe damage in Yogyakarta and the southern part of Central Java, Indonesia. The understanding of seismic response of earthquake among ground shaking and the level of building damage is important. We present numerical modeling of 3D seismic wave propagation around Yogyakarta and the southern part of Central Java using spectral-element method on MPI-GPU (Graphics Processing Unit) computer cluster to observe its seismic response due to the earthquake. The homogeneous 3D realistic model is generated with detailed topography surface. The influences of free surface topography and layer discontinuity of the 3D model among the seismic response are observed. The seismic wave field is discretized using spectral-element method. The spectral-element method is solved on a mesh of hexahedral elements that is adapted to the free surface topography and the internal discontinuity of the model. To increase the data processing capabilities, the simulation is performed on a GPU cluster with implementation of MPI (Message Passing Interface).

CityGML - Interoperable semantic 3D city models

Science.gov (United States)

Gröger, Gerhard; Plümer, Lutz

2012-07-01

CityGML is the international standard of the Open Geospatial Consortium (OGC) for the representation and exchange of 3D city models. It defines the three-dimensional geometry, topology, semantics and appearance of the most relevant topographic objects in urban or regional contexts. These definitions are provided in different, well-defined Levels-of-Detail (multiresolution model). The focus of CityGML is on the semantical aspects of 3D city models, its structures, taxonomies and aggregations, allowing users to employ virtual 3D city models for advanced analysis and visualization tasks in a variety of application domains such as urban planning, indoor/outdoor pedestrian navigation, environmental simulations, cultural heritage, or facility management. This is in contrast to purely geometrical/graphical models such as KML, VRML, or X3D, which do not provide sufficient semantics. CityGML is based on the Geography Markup Language (GML), which provides a standardized geometry model. Due to this model and its well-defined semantics and structures, CityGML facilitates interoperable data exchange in the context of geo web services and spatial data infrastructures. Since its standardization in 2008, CityGML has become used on a worldwide scale: tools from notable companies in the geospatial field provide CityGML interfaces. Many applications and projects use this standard. CityGML is also having a strong impact on science: numerous approaches use CityGML, particularly its semantics, for disaster management, emergency responses, or energy-related applications as well as for visualizations, or they contribute to CityGML, improving its consistency and validity, or use CityGML, particularly its different Levels-of-Detail, as a source or target for generalizations. This paper gives an overview of CityGML, its underlying concepts, its Levels-of-Detail, how to extend it, its applications, its likely future development, and the role it plays in scientific research. Furthermore, its

Comparison between a 1D and a 2D numerical model of an active magnetic regenerative refrigerator

DEFF Research Database (Denmark)

Petersen, Thomas Frank; Engelbrecht, Kurt; Bahl, Christian Robert Haffenden

2008-01-01

The active magnetic regenerator (AMR) refrigeration system represents an environmentally attractive alternative to vapour-compression refrigeration. This paper compares the results of two numerical AMR models: (1) a 1D finite difference model and (2) a 2D finite element model. Both models simulate...... a reciprocating AMR and can determine the cyclical steady-state temperature profile of the system as well as performance parameters such as the refrigeration capacity, the work input and the coefficient of performance (COP). The models are used to analyse an AMR with a regenerator made of flat parallel plates...... results of overall results such as the refrigeration capacity but that a 2D model is required for a detailed analysis of the phenomena occurring inside the AMR....

3D Strain Modelling of Tear Fault Analogues

Science.gov (United States)

Hindle, D.; Vietor, T.

2005-12-01

Tear faults can be described as vertical discontinuities, with near fault parallel displacements terminating on some sort of shallow detachment. As such, they are difficult to study in "cross section" i.e. 2 dimensions as is often the case for fold-thrust systems. Hence, little attempt has been made to model the evolution of strain around tear faults and the processes of strain localisation in such structures due to the necessity of describing these systems in 3 dimensions and the problems this poses for both numerical and analogue modelling. Field studies suggest that strain in such regions can be distributed across broad zones on minor tear systems, which are often not easily mappable. Such strain is probably assumed to be due to distributed strain and to displacement gradients which are themselves necessary for the initiation of the tear itself. We present a numerical study of the effects of a sharp, basal discontinutiy parallel to the transport direction in a shortening wedge of material. The discontinuity is represented by two adjacent basal surfaces with strongly contrasting (0.5 and 0.05) friction coefficient. The material is modelled using PFC3D distinct element software for simulating granular material, whose properties are chosen to simulate upper crustal, sedimentary rock. The model geometry is a rectangular bounding box, 2km x 1km, and 0.35-0.5km deep, with a single, driving wall of constant velocity. We show the evolution of strain in the model in horizontal and vertical sections, and interpret strain localization as showing the spontaneous development of tear fault like features. The strain field in the model is asymmetrical, rotated towards the strong side of the model. Strain increments seem to oscillate in time, suggesting achievement of a steady state. We also note that our model cannot be treated as a critical wedge, since the 3rd dimension and the lateral variations of strength rule out this type of 2D approximation.

3-D electromagnetic modeling for very early time sounding of shallow targets using integral equations

International Nuclear Information System (INIS)

Xiong, Z.; Tripp, A.C.

1994-01-01

This paper presents an integral equation algorithm for 3D EM modeling at high frequencies for applications in engineering an environmental studies. The integral equation method remains the same for low and high frequencies, but the dominant roles of the displacements currents complicate both numerical treatments and interpretations. With singularity extraction technique they successively extended the application of the Hankel filtering technique to the computation of Hankel integrals occurring in high frequency EM modeling. Time domain results are calculated from frequency domain results via Fourier transforms. While frequency domain data are not obvious for interpretations, time domain data show wave-like pictures that resemble seismograms. Both 1D and 3D numerical results show clearly the layer interfaces

The 3D-city model

DEFF Research Database (Denmark)

Holmgren, Steen; Rüdiger, Bjarne; Tournay, Bruno

2001-01-01

We have worked with the construction and use of 3D city models for about ten years. This work has given us valuable experience concerning model methodology. In addition to this collection of knowledge, our perception of the concept of city models has changed radically. In order to explain...... of 3D city models....

Flow above and within granular media composed of spherical and non-spherical particles - using a 3D numerical model

Science.gov (United States)

Bartzke, Gerhard; Kuhlmann, Jannis; Huhn, Katrin

2016-04-01

The entrainment of single grains and, hence, their erosion characteristics are dependent on fluid forcing, grain size and density, but also shape variations. To quantitatively describe and capture the hydrodynamic conditions around individual grains, researchers commonly use empirical approaches such as laboratory flume tanks. Nonetheless, it is difficult with such physical experiments to measure the flow velocities in the direct vicinity or within the pore spaces of sediments, at a sufficient resolution and in a non-invasive way. As a result, the hydrodynamic conditions in the water column, at the fluid-porous interface and within pore spaces of a granular medium of various grain shapes is not yet fully understood. For that reason, there is a strong need for numerical models, since these are capable of quantifying fluid speeds within a granular medium. A 3D-SPH (Smooth Particle Hydrodynamics) numerical wave tank model was set up to provide quantitative evidence on the flow velocities in the direct vicinity and in the interior of granular beds composed of two shapes as a complementary method to the difficult task of in situ measurement. On the basis of previous successful numerical wave tank models with SPH, the model geometry was chosen in dimensions of X=2.68 [m], Y=0.48 [m], and Z=0.8 [m]. Three suites of experiments were designed with a range of particle shape models: (1) ellipsoids with the long axis oriented in the across-stream direction, (2) ellipsoids with the long axis oriented in the along-stream direction, and (3) spheres. Particle diameters ranged from 0.04 [m] to 0.08 [m]. A wave was introduced by a vertical paddle that accelerated to 0.8 [m/s] perpendicular to the granular bed. Flow measurements showed that the flow velocity values into the beds were highest when the grains were oriented across the stream direction and lowest in case when the grains were oriented parallel to the stream, indicating that the model was capable to simulate simultaneously

3D numerical simulations of negative hydrogen ion extraction using realistic plasma parameters, geometry of the extraction aperture and full 3D magnetic field map

Science.gov (United States)

Mochalskyy, S.; Wünderlich, D.; Ruf, B.; Franzen, P.; Fantz, U.; Minea, T.

2014-02-01

Decreasing the co-extracted electron current while simultaneously keeping negative ion (NI) current sufficiently high is a crucial issue on the development plasma source system for ITER Neutral Beam Injector. To support finding the best extraction conditions the 3D Particle-in-Cell Monte Carlo Collision electrostatic code ONIX (Orsay Negative Ion eXtraction) has been developed. Close collaboration with experiments and other numerical models allows performing realistic simulations with relevant input parameters: plasma properties, geometry of the extraction aperture, full 3D magnetic field map, etc. For the first time ONIX has been benchmarked with commercial positive ions tracing code KOBRA3D. A very good agreement in terms of the meniscus position and depth has been found. Simulation of NI extraction with different e/NI ratio in bulk plasma shows high relevance of the direct negative ion extraction from the surface produced NI in order to obtain extracted NI current as in the experimental results from BATMAN testbed.

3D numerical simulation of transient processes in hydraulic turbines

International Nuclear Information System (INIS)

Cherny, S; Chirkov, D; Lapin, V; Eshkunova, I; Bannikov, D; Avdushenko, A; Skorospelov, V

2010-01-01

An approach for numerical simulation of 3D hydraulic turbine flows in transient operating regimes is presented. The method is based on a coupled solution of incompressible RANS equations, runner rotation equation, and water hammer equations. The issue of setting appropriate boundary conditions is considered in detail. As an illustration, the simulation results for runaway process are presented. The evolution of vortex structure and its effect on computed runaway traces are analyzed.

3D numerical simulation of transient processes in hydraulic turbines

Science.gov (United States)

Cherny, S.; Chirkov, D.; Bannikov, D.; Lapin, V.; Skorospelov, V.; Eshkunova, I.; Avdushenko, A.

2010-08-01

An approach for numerical simulation of 3D hydraulic turbine flows in transient operating regimes is presented. The method is based on a coupled solution of incompressible RANS equations, runner rotation equation, and water hammer equations. The issue of setting appropriate boundary conditions is considered in detail. As an illustration, the simulation results for runaway process are presented. The evolution of vortex structure and its effect on computed runaway traces are analyzed.

A Unit-Cell Model for Predicting the Elastic Constants of 3D Four Directional Cylindrical Braided Composite Shafts

Science.gov (United States)

Hao, Wenfeng; Liu, Ye; Huang, Xinrong; Liu, Yinghua; Zhu, Jianguo

2018-06-01

In this work, the elastic constants of 3D four directional cylindrical braided composite shafts were predicted using analytical and numerical methods. First, the motion rule of yarn carrier of 3D four directional cylindrical braided composite shafts was analyzed, and the horizontal projection of yarn motion trajectory was obtained. Then, the geometry models of unit-cells with different braiding angles and fiber volume contents were built up, and the meso-scale models of 3D cylindrical braided composite shafts were obtained. Finally, the effects of braiding angles and fiber volume contents on the elastic constants of 3D braided composite shafts were analyzed theoretically and numerically. These results play a crucial role in investigating the mechanical properties of 3D 4-directional braided composites shafts.

Do-It-Yourself: 3D Models of Hydrogenic Orbitals through 3D Printing

Science.gov (United States)

Griffith, Kaitlyn M.; de Cataldo, Riccardo; Fogarty, Keir H.

2016-01-01

Introductory chemistry students often have difficulty visualizing the 3-dimensional shapes of the hydrogenic electron orbitals without the aid of physical 3D models. Unfortunately, commercially available models can be quite expensive. 3D printing offers a solution for producing models of hydrogenic orbitals. 3D printing technology is widely…

Numerical method for calculation of 3D viscous turbomachine flow taking into account stator/rotor unsteady interaction

Energy Technology Data Exchange (ETDEWEB)

Rusanov, A V; Yershov, S V [Institute of Mechanical Engineering Problems of National Academy of Sciences of Ukraine Kharkov (Ukraine)

1998-12-31

The numerical method is suggested for the calculation of the 3D periodically unsteady viscous cascade flow evoked by the aerodynamics interaction of blade rows. Such flow is described by the thin-layer Reynolds-averaged unsteady Navier-Stokes equations. The turbulent effects are simulated with the modified Baldwin-Lomax turbulence model. The problem statement allows to consider an unsteady flow through either a single turbo-machine stage or a multi stage turbomachine. The sliding mesh techniques and the time-space non-oscillatory square interpolation are used in axial spacings to calculate the flow in a computational domain that contains the reciprocally moving elements. The gasdynamical equations are integrated numerically with the implicit quasi-monotonous Godunov`s type ENO scheme of the second or third order of accuracy. The suggested numerical method is incorporated in the FlowER code developed by authors for calculations of the 3D viscous compressible flows through multi stage turbomachines. The numerical results are presented for unsteady turbine stage throughflows. The method suggested is shown to simulate qualitatively properly the main unsteady cascade effects in particular the periodically blade loadings, the propagation of stator wakes through rotor blade passage and the unsteady temperature flowfields for stages with cooled stator blades. (author) 21 refs.

Numerical method for calculation of 3D viscous turbomachine flow taking into account stator/rotor unsteady interaction

Energy Technology Data Exchange (ETDEWEB)

Rusanov, A.V.; Yershov, S.V. [Institute of Mechanical Engineering Problems of National Academy of Sciences of Ukraine Kharkov (Ukraine)

1997-12-31

The numerical method is suggested for the calculation of the 3D periodically unsteady viscous cascade flow evoked by the aerodynamics interaction of blade rows. Such flow is described by the thin-layer Reynolds-averaged unsteady Navier-Stokes equations. The turbulent effects are simulated with the modified Baldwin-Lomax turbulence model. The problem statement allows to consider an unsteady flow through either a single turbo-machine stage or a multi stage turbomachine. The sliding mesh techniques and the time-space non-oscillatory square interpolation are used in axial spacings to calculate the flow in a computational domain that contains the reciprocally moving elements. The gasdynamical equations are integrated numerically with the implicit quasi-monotonous Godunov`s type ENO scheme of the second or third order of accuracy. The suggested numerical method is incorporated in the FlowER code developed by authors for calculations of the 3D viscous compressible flows through multi stage turbomachines. The numerical results are presented for unsteady turbine stage throughflows. The method suggested is shown to simulate qualitatively properly the main unsteady cascade effects in particular the periodically blade loadings, the propagation of stator wakes through rotor blade passage and the unsteady temperature flowfields for stages with cooled stator blades. (author) 21 refs.

A new general purpose event horizon finder for 3D numerical spacetimes

International Nuclear Information System (INIS)

Diener, Peter

2003-01-01

I present a new general purpose event horizon finder for full 3D numerical spacetimes. It works by evolving a complete null surface backwards in time. The null surface is described as the zero-level set of a scalar function, which in principle is defined everywhere. This description of the surface allows the surface, trivially, to change topology, making this event horizon finder able to handle numerical spacetimes where two (or more) black holes merge into a single final black hole

Radiosity diffusion model in 3D

Science.gov (United States)

Riley, Jason D.; Arridge, Simon R.; Chrysanthou, Yiorgos; Dehghani, Hamid; Hillman, Elizabeth M. C.; Schweiger, Martin

2001-11-01

We present the Radiosity-Diffusion model in three dimensions(3D), as an extension to previous work in 2D. It is a method for handling non-scattering spaces in optically participating media. We present the extension of the model to 3D including an extension to the model to cope with increased complexity of the 3D domain. We show that in 3D more careful consideration must be given to the issues of meshing and visibility to model the transport of light within reasonable computational bounds. We demonstrate the model to be comparable to Monte-Carlo simulations for selected geometries, and show preliminary results of comparisons to measured time-resolved data acquired on resin phantoms.

Development of a 3D bone marrow adipose tissue model.

Science.gov (United States)

Fairfield, Heather; Falank, Carolyne; Farrell, Mariah; Vary, Calvin; Boucher, Joshua M; Driscoll, Heather; Liaw, Lucy; Rosen, Clifford J; Reagan, Michaela R

2018-01-26

targets. In addition, proteomic characterization as well as microarray data (expression of >22,000 genes) coupled with KEGG pathway analysis and gene set expression analysis (GSEA) supported our development of less-inflammatory 3D BMAT compared to 2D culture. In sum, we developed the first 3D, tissue-engineered bone marrow adipose tissue model, which is a versatile, novel model that can be used to study numerous diseases and biological processes involved with the bone marrow. Copyright © 2018. Published by Elsevier Inc.

Radon transport modelling: User's guide to RnMod3d

Energy Technology Data Exchange (ETDEWEB)

Andersen, C.E

2000-08-01

RnMod3d is a numerical computer model of soil-gas and radon transport in porous media. It can be used, for example, to study radon entry from soil into houses in response to indoor-outdoor pressure differences or changes in atmospheric pressure. It can also be used for flux calculations of radon from the soil surface or to model radon exhalation from building materials such as concrete. The finite-volume model is a technical research tool, and it cannot be used meaningfully without good understanding of the involved physical equations. Some understanding of numerical mathematics and the programming language Pascal is also required. Originally, the code was developed for internal use at Risoe only. With this guide, however, it should be possible for others to use the model. Three-dimensional steady-state or transient problems with Darcy flow of soil gas and combined generation, radioactive decay, diffusion and advection of radon can be solved. Moisture is included in the model, and partitioning of radon between air, water and soil grains (adsorption) is taken into account. Most parameters can change in time and space, and transport parameters (diffusivity and permeability) may be anisotropic. This guide includes benchmark tests based on simple problems with known solutions. RnMod3d has also been part of an international model intercomparison exercise based on more complicated problems without known solutions. All tests show that RnMod3d gives results of good quality. (au)

3D random Voronoi grain-based models for simulation of brittle rock damage and fabric-guided micro-fracturing

Directory of Open Access Journals (Sweden)

E. Ghazvinian

2014-12-01

Full Text Available A grain-based distinct element model featuring three-dimensional (3D Voronoi tessellations (random poly-crystals is proposed for simulation of crack damage development in brittle rocks. The grain boundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rock and allow for numerical replication of crack damage progression through initiation and propagation of micro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the past for brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi models has limited its application to two-dimensional (2D codes. The proposed approach is implemented in Neper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files that can be read directly into 3DEC. A series of Unconfined Compressive Strength (UCS tests are simulated in 3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate the relationship between each micro-parameter and the model's macro-response. The possibility of numerical replication of the classical U-shape strength curve for anisotropic rocks is also investigated in numerical UCS tests by using complex-shaped (elongated grains that are cemented to one another along their adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models for accurate replication of the mechanical behaviour of isotropic and anisotropic (containing a fabric rocks.

3D-radiative transfer in terrestrial atmosphere: An efficient parallel numerical procedure

Science.gov (United States)

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

Influence of crosslinking on the mechanical behavior of 3D printed alginate scaffolds: Experimental and numerical approaches.

Science.gov (United States)

Naghieh, Saman; Karamooz-Ravari, Mohammad Reza; Sarker, M D; Karki, Eva; Chen, Xiongbiao

2018-04-01

Tissue scaffolds fabricated by three-dimensional (3D) bioprinting are attracting considerable attention for tissue engineering applications. Because the mechanical properties of hydrogel scaffolds should match the damaged tissue, changing various parameters during 3D bioprinting has been studied to manipulate the mechanical behavior of the resulting scaffolds. Crosslinking scaffolds using a cation solution (such as CaCl 2 ) is also important for regulating the mechanical properties, but has not been well documented in the literature. Here, the effect of varied crosslinking agent volume and crosslinking time on the mechanical behavior of 3D bioplotted alginate scaffolds was evaluated using both experimental and numerical methods. Compression tests were used to measure the elastic modulus of each scaffold, then a finite element model was developed and a power model used to predict scaffold mechanical behavior. Results showed that crosslinking time and volume of crosslinker both play a decisive role in modulating the mechanical properties of 3D bioplotted scaffolds. Because mechanical properties of scaffolds can affect cell response, the findings of this study can be implemented to modulate the elastic modulus of scaffolds according to the intended application. Copyright © 2018 Elsevier Ltd. All rights reserved.

3D Face Apperance Model

DEFF Research Database (Denmark)

Lading, Brian; Larsen, Rasmus; Astrom, K

2006-01-01

We build a 3D face shape model, including inter- and intra-shape variations, derive the analytical Jacobian of its resulting 2D rendered image, and show example of its fitting performance with light, pose, id, expression and texture variations......We build a 3D face shape model, including inter- and intra-shape variations, derive the analytical Jacobian of its resulting 2D rendered image, and show example of its fitting performance with light, pose, id, expression and texture variations...

Hybrid 2D-3D modelling of GTA welding with filler wire addition

KAUST Repository

Traidia, Abderrazak

2012-07-01

A hybrid 2D-3D model for the numerical simulation of Gas Tungsten Arc welding is proposed in this paper. It offers the possibility to predict the temperature field as well as the shape of the solidified weld joint for different operating parameters, with relatively good accuracy and reasonable computational cost. Also, an original approach to simulate the effect of immersing a cold filler wire in the weld pool is presented. The simulation results reveal two important observations. First, the weld pool depth is locally decreased in the presence of filler metal, which is due to the energy absorption by the cold feeding wire from the hot molten pool. In addition, the weld shape, maximum temperature and thermal cycles in the workpiece are relatively well predicted even when a 2D model for the arc plasma region is used. © 2012 Elsevier Ltd. All rights reserved.

3D printing application and numerical simulations in a fracture system

Science.gov (United States)

Yoon, H.; Martinez, M. J.

2017-12-01

The hydrogeological and mechanical properties in fractured and porous media are fundamental to predicting coupled multiphysics processes in the subsurface. Recent advances in experimental methods and multi-scale imaging capabilities have revolutionized our ability to quantitatively characterize geomaterials and digital counterparts are now routinely used for numerical simulations to characterize petrophysical and mechanical properties across scales. 3D printing is a very effective and creative technique that reproduce the digital images in a controlled way. For geoscience applications, 3D printing can be co-opted to print reproducible porous and fractured structures derived from CT-imaging of actual rocks and theoretical algorithms for experimental testing. In this work we used a stereolithography (SLA) method to create a single fracture network. The fracture in shale was first scanned using a microCT system and then the digital fracture network was printed into two parts and assembled. Aperture ranges from 0.3 to 1 mm. In particular, we discuss the design of single fracture network and the progress of printing practices to reproduce the fracture network system. Printed samples at different scales are used to measure the permeability and surface roughness. Various numerical simulations including (non-)reactive transport and multiphase flow cases are performed to study fluid flow characterization. We will also discuss the innovative advancement of 3D printing techniques applicable for coupled processes in the subsurface. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Elastic oscillations of water column in the 2003 Tokachi-oki tsunami source: in-situ measurements and 3-D numerical modelling

Directory of Open Access Journals (Sweden)

M. A. Nosov

2007-01-01

Full Text Available During the 2003 Tokachi-Oki tsunamigenic earthquake the real-time JAMSTEC observatory obtained records which provided a unique opportunity to have a look deep inside the tsunami source. Considering water column as a compressible medium we processed the bottom pressure records in order to estimate amplitude, duration and velocity of bottom displacement. Spectral analysis of the records revealed a clear manifestation of the low-frequency elastic oscillations of water column. We also presented 3-D finite-difference numerical model developed in the framework of linear potential theory of ideal compressible fluid to better understand dynamical processes in the tsunami source. The model reproduces position of the main spectral maximum rather correctly. However, due to neglecting of crust elasticity and to lack of exact knowledge of spatiotemporal laws of bottom motion, there is an essential difference between in-situ observed and computed spectra.

3D modelling for multipurpose cadastre

NARCIS (Netherlands)

Abduhl Rahman, A.; Van Oosterom, P.J.M.; Hua, T.C.; Sharkawi, K.H.; Duncan, E.E.; Azri, N.; Hassan, M.I.

2012-01-01

Three-dimensional (3D) modelling of cadastral objects (such as legal spaces around buildings, around utility networks and other spaces) is one of the important aspects for a multipurpose cadastre (MPC). This paper describes the 3D modelling of the objects for MPC and its usage to the knowledge of 3D

Multi-view and 3D deformable part models.

Science.gov (United States)

Pepik, Bojan; Stark, Michael; Gehler, Peter; Schiele, Bernt

2015-11-01

As objects are inherently 3D, they have been modeled in 3D in the early days of computer vision. Due to the ambiguities arising from mapping 2D features to 3D models, 3D object representations have been neglected and 2D feature-based models are the predominant paradigm in object detection nowadays. While such models have achieved outstanding bounding box detection performance, they come with limited expressiveness, as they are clearly limited in their capability of reasoning about 3D shape or viewpoints. In this work, we bring the worlds of 3D and 2D object representations closer, by building an object detector which leverages the expressive power of 3D object representations while at the same time can be robustly matched to image evidence. To that end, we gradually extend the successful deformable part model [1] to include viewpoint information and part-level 3D geometry information, resulting in several different models with different level of expressiveness. We end up with a 3D object model, consisting of multiple object parts represented in 3D and a continuous appearance model. We experimentally verify that our models, while providing richer object hypotheses than the 2D object models, provide consistently better joint object localization and viewpoint estimation than the state-of-the-art multi-view and 3D object detectors on various benchmarks (KITTI [2] , 3D object classes [3] , Pascal3D+ [4] , Pascal VOC 2007 [5] , EPFL multi-view cars[6] ).

3D Gravity Modeling of Complex Salt Features in the Southern Gulf of Mexico

Directory of Open Access Journals (Sweden)

Mauricio Nava-Flores

2016-01-01

Full Text Available We present a three-dimensional (3D gravity modeling and inversion approach and its application to complex geological settings characterized by several allochthonous salt bodies embedded in terrigenous sediments. Synthetic gravity data were computed for 3D forward modeling of salt bodies interpreted from Prestack Depth Migration (PSDM seismic images. Density contrasts for the salt bodies surrounded by sedimentary units are derived from density-compaction curves for the northern Gulf of Mexico’s oil exploration surveys. By integrating results from different shape- and depth-source estimation algorithms, we built an initial model for the gravity anomaly inversion. We then applied a numerically optimized 3D simulated annealing gravity inversion method. The inverted 3D density model successfully retrieves the synthetic salt body ensemble. Results highlight the significance of integrating high-resolution potential field data for salt and subsalt imaging in oil exploration.

3D finite element modelling of force transmission and particle fracture of sand

Energy Technology Data Exchange (ETDEWEB)

Imseeh, Wadi H.; Alshibli, Khalid A. (Tennessee-K)

2018-02-01

Global compressive loading of granular media causes rearrangements of particles into a denser configuration. Under 1D compression, researchers observed that particles initially translate and rotate which lead to more contacts between particles and the development of force chains to resist applied loads. Particles within force chains resist most of the applied loads while neighbor particles provide lateral support to prevent particles within force chains from buckling. Several experimental and numerical models have been proposed in the literature to characterize force chains within granular materials. This paper presents a 3D finite element (FE) model that simulates 1D compression experiment on F-75 Ottawa sand. The FE mesh of particles closely matched 3D physical shape of sand particles that were acquired using 3D synchrotron micro-computed tomography (SMT) technique. The paper presents a quantitative assessment of the model, in which evolution of force chains, fracture modes, and stress-strain relationships showed an excellent agreement with experimental measurements reported by Cil et al. Alshibli (2017).

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

International Nuclear Information System (INIS)

Chacon, L.

2007-01-01

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

USING 3D MODELLING AND NUMERICAL SIMULATION TO OPTIMIZE LARGE WASTEWATER INSTALLATIONS

Directory of Open Access Journals (Sweden)

NEDELCUŢ Florin

2014-09-01

the evaluating the wastewater treatment efficiency. Another benefit of the CFD analysis was the possibility to visualize the magnitude of the turbulences and large swirls existing inside the clarifiers, as their distribution. With the help of the 3D modelling and CFD we managed to compare two operation cases, in order to better estimate the interval between the required sludge cleaning operations. This more accurate view of operation of final clarifiers, conducted by the CFD analysis presented here, will provide theoretical and technical basis for future proposals when modifying and improving the design of such clarifiers, in the benefit of the designers, engineers and users.

Full 3-D OCT-based pseudophakic custom computer eye model

Science.gov (United States)

Sun, M.; Pérez-Merino, P.; Martinez-Enriquez, E.; Velasco-Ocana, M.; Marcos, S.

2016-01-01

We compared measured wave aberrations in pseudophakic eyes implanted with aspheric intraocular lenses (IOLs) with simulated aberrations from numerical ray tracing on customized computer eye models, built using quantitative 3-D OCT-based patient-specific ocular geometry. Experimental and simulated aberrations show high correlation (R = 0.93; poptical geometrical and surgically-related factors to image quality, and are an excellent tool for characterizing and improving cataract surgery. PMID:27231608

3D modeling based on CityEngine

Science.gov (United States)

Jia, Guangyin; Liao, Kaiju

2017-03-01

Currently, there are many 3D modeling softwares, like 3DMAX, AUTOCAD, and more populous BIM softwares represented by REVIT. CityEngine modeling software introduced in this paper can fully utilize the existing GIS data and combine other built models to make 3D modeling on internal and external part of buildings in a rapid and batch manner, so as to improve the 3D modeling efficiency.

The mechanical properties of human dentin for 3-D finite element modeling: Numerical and analytical evaluation.

Science.gov (United States)

Grzebieluch, Wojciech; Będziński, Romuald; Czapliński, Tomasz; Kaczmarek, Urszula

2017-07-01

The FEM is often used in investigations of dentin loading conditions; however, its anisotropy is mostly neglected. The purpose of the study was to evaluate the anisotropy and the elastic properties of an equivalent homogenous material model of human dentin as well as to compare isotropic and anisotropic dentin FE-models. Analytical and numerical dentin homogenization according to Luciano and Barbero was performed and E-modulus (E), Poisson's ratios (v) G-modulus (G) were calculated. The E-modulus of the dentin matrix was 28.0 GPa, Poisson's ratio (v) was 0.3; finite element models of orthotropic and isotropic dentin were created, loaded and compared using Ansys® 14.5 and CodeAster® 11.2 software. Anisotropy of the dentin ranged from 6.9 to 35.2%. E-modulus and G-modulus were as follows: E1 = 22.0-26.0 GPa, E2/E3 = 15.7-23.0 GPa; G12/G13 = 6.96-9.35 GPa and G23 = 6.08-8.09 GPa (highest values in the superficial layer). In FEM analysis of the displacement values were higher in the isotropic than in the orthotropic model, reaching up to 16% by shear load, 37% by compression and 23% in the case of shear with bending. Strain values were higher in the isotropic model, up to 35% for the shear load, 31% for compression and 35% in the case of shear with bending. The decrease in the volumetric fraction and diameter of tubules increased the G and E values. Anisotropy of the dentin applied during FEM analysis decreased the displacements and strain values. The numerical and analytical homogenization of dentin showed similar results.

3D multidisciplinary numerical model of polychlorinated biphenyl dynamics on the Black Sea north-western shelf

Science.gov (United States)

Bagaiev, Andrii; Ivanov, Vitaliy

2014-05-01

The Black Sea north-western shelf plays a key role in economics of the developing countries such as Ukraine due to food supply, invaluable recreational potential and variety of the relevant maritime shipping routes. On the other hand, a shallow flat shelf is mostly affected by anthropogenic pollution, eutrophication, hypoxia and harmful algae blooms. The research is focused on modeling the transport and transformation of PCBs (PolyChlorinated Biphenyls) because they are exceedingly toxic and highly resistant to degradation, hence cumulatively affect marine ecosystems. Being lipophilic compounds, PCBs demonstrate the distinguishing sorption/desorption activity taking part in the biogeochemical fluxes via the organic matter particles and sediments. In the framework of the research, the coastal in-situ data on PCB concentration in the water column and sediments are processed, visualized and analyzed. It is concluded that the main sources of PCBs are related to the Danube discharge and resuspension from the shallow-water sediments. Developed 3D numerical model is aimed at simulation of PCB contamination of the water column and sediment. The model integrates the full physics hydrodynamic block as well as modules, which describe detritus transport and transformation and PCB dynamics. Three state variables are simulated in PCB transport module: concentration in solute, on the settling particles of detritus and in the top layer of sediments. PCB adsorption/desorption on detritus; the reversible PCB fluxes at the water-sediment boundary; destruction of detritus are taken into consideration. Formalization of PCB deposition/resuspension in the sediments is adapted from Van Rijn's model of the suspended sediment transport. The model was spun up to reconstruct the short term scenario of the instantaneous PCB release from the St. George Arm of Danube. It has been shown that PCB transport on sinking detritus represents the natural buffer mechanism damping the spreading PCB

Numerical relativity for D dimensional axially symmetric space-times: Formalism and code tests

International Nuclear Information System (INIS)

Zilhao, Miguel; Herdeiro, Carlos; Witek, Helvi; Nerozzi, Andrea; Sperhake, Ulrich; Cardoso, Vitor; Gualtieri, Leonardo

2010-01-01

The numerical evolution of Einstein's field equations in a generic background has the potential to answer a variety of important questions in physics: from applications to the gauge-gravity duality, to modeling black hole production in TeV gravity scenarios, to analysis of the stability of exact solutions, and to tests of cosmic censorship. In order to investigate these questions, we extend numerical relativity to more general space-times than those investigated hitherto, by developing a framework to study the numerical evolution of D dimensional vacuum space-times with an SO(D-2) isometry group for D≥5, or SO(D-3) for D≥6. Performing a dimensional reduction on a (D-4) sphere, the D dimensional vacuum Einstein equations are rewritten as a 3+1 dimensional system with source terms, and presented in the Baumgarte, Shapiro, Shibata, and Nakamura formulation. This allows the use of existing 3+1 dimensional numerical codes with small adaptations. Brill-Lindquist initial data are constructed in D dimensions and a procedure to match them to our 3+1 dimensional evolution equations is given. We have implemented our framework by adapting the Lean code and perform a variety of simulations of nonspinning black hole space-times. Specifically, we present a modified moving puncture gauge, which facilitates long-term stable simulations in D=5. We further demonstrate the internal consistency of the code by studying convergence and comparing numerical versus analytic results in the case of geodesic slicing for D=5, 6.

Numerical modelling techniques of soft soil improvement via stone columns: A brief review

Science.gov (United States)

Zukri, Azhani; Nazir, Ramli

2018-04-01

There are a number of numerical studies on stone column systems in the literature. Most of the studies found were involved with two-dimensional analysis of the stone column behaviour, while only a few studies used three-dimensional analysis. The most popular software utilised in those studies was Plaxis 2D and 3D. Other types of software that used for numerical analysis are DIANA, EXAMINE, ZSoil, ABAQUS, ANSYS, NISA, GEOSTUDIO, CRISP, TOCHNOG, CESAR, GEOFEM (2D & 3D), FLAC, and FLAC 3. This paper will review the methodological approaches to model stone column numerically, both in two-dimensional and three-dimensional analyses. The numerical techniques and suitable constitutive model used in the studies will also be discussed. In addition, the validation methods conducted were to verify the numerical analysis conducted will be presented. This review paper also serves as a guide for junior engineers through the applicable procedures and considerations when constructing and running a two or three-dimensional numerical analysis while also citing numerous relevant references.

Anti-3D Weapon Model Detection for Safe 3D Printing Based on Convolutional Neural Networks and D2 Shape Distribution

Directory of Open Access Journals (Sweden)

Giao N. Pham

2018-03-01

Full Text Available With the development of 3D printing, weapons are easily printed without any restriction from the production managers. Therefore, anti-3D weapon model detection is necessary issue in safe 3D printing to prevent the printing of 3D weapon models. In this paper, we would like to propose an anti-3D weapon model detection algorithm to prevent the printing of anti-3D weapon models for safe 3D printing based on the D2 shape distribution and an improved convolutional neural networks (CNNs. The purpose of the proposed algorithm is to detect anti-3D weapon models when they are used in 3D printing. The D2 shape distribution is computed from random points on the surface of a 3D weapon model and their geometric features in order to construct a D2 vector. The D2 vector is then trained by improved CNNs. The CNNs are used to detect anti-3D weapon models for safe 3D printing by training D2 vectors which have been constructed from the D2 shape distribution of 3D weapon models. Experiments with 3D weapon models proved that the D2 shape distribution of 3D weapon models in the same class is the same. Training and testing results also verified that the accuracy of the proposed algorithm is higher than the conventional works. The proposed algorithm is applied in a small application, and it could detect anti-3D weapon models for safe 3D printing.

UT simulation using a fully automated 3D hybrid model: Application to planar backwall breaking defects inspection

Science.gov (United States)

Imperiale, Alexandre; Chatillon, Sylvain; Darmon, Michel; Leymarie, Nicolas; Demaldent, Edouard

2018-04-01

The high frequency models gathered in the CIVA software allow fast computations and provide satisfactory quantitative predictions in a wide range of situations. However, the domain of validity of these models is limited since they do not accurately predict the ultrasound response in configurations involving subwavelength complex phenomena. In addition, when modelling backwall breaking defects inspection, an important challenge remains to capture the propagation of the creeping waves that are generated at the critical angle. Hybrid models combining numerical and asymptotic methods have already been shown to be an effective strategy to overcome these limitations in 2D [1]. However, 3D simulations remain a crucial issue for industrial applications because of the computational cost of the numerical solver. A dedicated three dimensional high order finite element model combined with a domain decomposition method has been recently proposed to tackle 3D limitations [2]. In this communication, we will focus on the specific case of planar backwall breaking defects, with an adapted coupling strategy in order to efficiently model the propagation of creeping waves. Numerical and experimental validations will be proposed on various configurations.

Insights from 3D numerical simulations on the dynamics of the India-Asia collision zone

Science.gov (United States)

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

Intrusion of Magmatic Bodies Into the Continental Crust: 3-D Numerical Models

Science.gov (United States)

Gorczyk, Weronika; Vogt, Katharina

2018-03-01

Magma intrusion is a major material transfer process in the Earth's continental crust. Yet the mechanical behavior of the intruding magma and its host are a matter of debate. In this study we present a series of numerical thermomechanical simulations on magma emplacement in 3-D. Our results demonstrate the response of the continental crust to magma intrusion. We observe change in intrusion geometries between dikes, cone sheets, sills, plutons, ponds, funnels, finger-shaped and stock-like intrusions, and injection time. The rheology and temperature of the host are the main controlling factors in the transition between these different modes of intrusion. Viscous deformation in the warm and deep crust favors host rock displacement and plutons at the crust-mantle boundary forming deep-seated plutons or magma ponds in the lower to middle crust. Brittle deformation in the cool and shallow crust induces cone-shaped fractures in the host rock and enables emplacement of finger- or stock-like intrusions at shallow or intermediate depth. Here the passage of magmatic and hydrothermal fluids from the intrusion through the fracture pattern may result in the formation of ore deposits. A combination of viscous and brittle deformation forms funnel-shaped intrusions in the middle crust. Intrusion of low-density magma may more over result in T-shaped intrusions in cross section with magma sheets at the surface.

Topology reconstruction for B-Rep modeling from 3D mesh in reverse engineering applications

Science.gov (United States)

Bénière, Roseline; Subsol, Gérard; Gesquière, Gilles; Le Breton, François; Puech, William

2012-03-01

Nowadays, most of the manufactured objects are designed using CAD (Computer-Aided Design) software. Nevertheless, for visualization, data exchange or manufacturing applications, the geometric model has to be discretized into a 3D mesh composed of a finite number of vertices and edges. But, in some cases, the initial model may be lost or unavailable. In other cases, the 3D discrete representation may be modified, for example after a numerical simulation, and does not correspond anymore to the initial model. A reverse engineering method is then required to reconstruct a 3D continuous representation from the discrete one. In previous work, we have presented a new approach for 3D geometric primitive extraction. In this paper, to complete our automatic and comprehensive reverse engineering process, we propose a method to construct the topology of the retrieved object. To reconstruct a B-Rep model, a new formalism is now introduced to define the adjacency relations. Then a new process is used to construct the boundaries of the object. The whole process is tested on 3D industrial meshes and bring a solution to recover B-Rep models.

From 2D to 3D turbulence through 2D3C configurations

Science.gov (United States)

Buzzicotti, Michele; Biferale, Luca; Linkmann, Moritz

2017-11-01

We study analytically and numerically the geometry of the nonlinear interactions and the resulting energy transfer directions of 2D3C flows. Through a set of suitably designed Direct Numerical Simulations we also study the coupling between several 2D3C flows, where we explore the transition between 2D and fully 3D turbulence. In particular, we find that the coupling of three 2D3C flows on mutually orthogonal planes subject to small-scale forcing leads to a stationary 3D out-of-equilibrium dynamics at the energy containing scales where the inverse cascade is directly balanced by a forward cascade carried by a different subsets of interactions. ERC AdG Grant No 339032 NewTURB.

3D Hilbert Space Filling Curves in 3D City Modeling for Faster Spatial Queries

DEFF Research Database (Denmark)

Ujang, Uznir; Antón Castro, Francesc/François; Azri, Suhaibah

2014-01-01

The advantages of three dimensional (3D) city models can be seen in various applications including photogrammetry, urban and regional planning, computer games, etc. They expand the visualization and analysis capabilities of Geographic Information Systems on cities, and they can be developed using...... method, retrieving portions of and especially searching these 3D city models, will not be done optimally. Even though current developments are based on an open data model allotted by the Open Geospatial Consortium (OGC) called CityGML, its XML-based structure makes it challenging to cluster the 3D urban...... objects. In this research, the authors propose an opponent data constellation technique of space-filling curves (3D Hilbert curves) for 3D city model data representation. Unlike previous methods, that try to project 3D or n-dimensional data down to 2D or 3D using Principal Component Analysis (PCA...

3D finite element modelling of sheet metal blanking process

Science.gov (United States)

Bohdal, Lukasz; Kukielka, Leon; Chodor, Jaroslaw; Kulakowska, Agnieszka; Patyk, Radoslaw; Kaldunski, Pawel

2018-05-01

The shearing process such as the blanking of sheet metals has been used often to prepare workpieces for subsequent forming operations. The use of FEM simulation is increasing for investigation and optimizing the blanking process. In the current literature a blanking FEM simulations for the limited capability and large computational cost of the three dimensional (3D) analysis has been largely limited to two dimensional (2D) plane axis-symmetry problems. However, a significant progress in modelling which takes into account the influence of real material (e.g. microstructure of the material), physical and technological conditions can be obtained by using 3D numerical analysis methods in this area. The objective of this paper is to present 3D finite element analysis of the ductile fracture, strain distribution and stress in blanking process with the assumption geometrical and physical nonlinearities. The physical, mathematical and computer model of the process are elaborated. Dynamic effects, mechanical coupling, constitutive damage law and contact friction are taken into account. The application in ANSYS/LS-DYNA program is elaborated. The effect of the main process parameter a blanking clearance on the deformation of 1018 steel and quality of the blank's sheared edge is analyzed. The results of computer simulations can be used to forecasting quality of the final parts optimization.

Referent 3D tumor model at cellular level in radionuclide therapy

International Nuclear Information System (INIS)

Spaic, R.; Ilic, R.D.; Petrovic, B.J.

2002-01-01

Aim Conventional internal dosimetry has a lot of limitations because of tumor dose nonuniformity. The best approach for absorbed dose at cellular level for different tumors in radionuclide therapy calculation is Monte Carlo method. The purpose of this study is to introduce referent tumor 3D model at cellular level for Monte Carlo simulation study in radionuclide therapy. Material and Methods The moment when tumor is detectable and when same therapy can start is time period in which referent 3D tumor model at cellular level was defined. In accordance with tumor growth rate at that moment he was a sphere with same radius (10 000 μm). In that tumor there are cells or cluster of cells, which are randomly distributed spheres. Distribution of cells/cluster of cells can be calculated from histology data but it was assumed that this distribution is normal with the same mean value and standard deviation (100±50 mm). Second parameter, which was selected to define referent tumor, is volume density of cells (30%). In this referent tumor there are no necroses. Stroma is defined as space between spheres with same concentration of materials as in spheres. Results: Referent tumor defined on this way have about 2,2 10 5 cells or cluster of cells random distributed. Using this referent 3D tumor model and for same concentration of radionuclides (1:100) and energy of beta emitters (1000 keV) which are homogeneously distributed in labeled cells absorbed dose for all cells was calculated. Simulations are done using FOTELP Monte Carlo code, which is modified for this purposes. Results of absorbed dose in cells are given in numerical values (1D distribution) and as the images (2D or 3D distributions). Conclusion Geometrical module for Monte Carlo simulation study can be standardized by introducing referent 3D tumor model at cellular level. This referent 3D tumor model gives most realistic presentation of different tumors at the moment of their detectability. Referent 3D tumor model at

a Proposal for Generalization of 3d Models

Science.gov (United States)

Uyar, A.; Ulugtekin, N. N.

2017-11-01

In recent years, 3D models have been created of many cities around the world. Most of the 3D city models have been introduced as completely graphic or geometric models, and the semantic and topographic aspects of the models have been neglected. In order to use 3D city models beyond the task, a generalization is necessary. CityGML is an open data model and XML-based format for the storage and exchange of virtual 3D city models. Level of Details (LoD) which is an important concept for 3D modelling, can be defined as outlined degree or prior representation of real-world objects. The paper aim is first describes some requirements of 3D model generalization, then presents problems and approaches that have been developed in recent years. In conclude the paper will be a summary and outlook on problems and future work.

NUMERICAL PREDICTION MODELS FOR AIR POLLUTION BY MOTOR VEHICLE EMISSIONS

Directory of Open Access Journals (Sweden)

M. M. Biliaiev

2016-12-01

Full Text Available Purpose. Scientific work involves: 1 development of 3D numerical models that allow calculating the process of air pollution by motor vehicles emissions; 2 creation of models which would allow predicting the air pollution level in urban areas. Methodology. To solve the problem upon assessing the level of air pollution by motor vehicles emissions fundamental equations of aerodynamics and mass transfer are used. For the solution of differential equations of aerodynamics and mass transfer finite-difference methods are used. For the numerical integration of the equation for the velocity potential the method of conditional approximations is applied. The equation for the velocity potential written in differential form, splits into two equations, where at each step of splitting an unknown value of the velocity potential is determined by an explicit scheme of running computation, while the difference scheme is implicit one. For the numerical integration of the emissions dispersion equation in the atmosphere applies the implicit alternating-triangular difference scheme of splitting. Emissions from the road are modeled by a series of point sources of given intensity. Developed numerical models form is the basis of the created software package. Findings. 3D numerical models were developed; they belong to the class of «diagnostic models». These models take into account main physical factors that influence the process of dispersion of harmful substances in the atmosphere when emissions from vehicles in the city occur. Based on the constructed numerical models the computational experiment was conducted to assess the level of air pollution in the street. Originality. Authors have developed numerical models that allow to calculate the 3D aerodynamics of the wind flow in urban areas and the process of mass transfer emissions from the highway. Calculations to determine the area of contamination, which is formed near the buildings, located along the highway were

The 3-dimensional core model DYN3D

Energy Technology Data Exchange (ETDEWEB)

Grundmann, U.; Mittag, S.; Rohde, U.

1999-01-01

Analyzing the safety margins in transients and accidents of nuclear reactors 3-dimensional models of the core were used to avoid conservative assumptions needed for point kinetics or 1-dimensional models. Therefore, the 3D code DYN3D has been developed for the analysis of reactivity initiated accidents (RIA) in thermal nuclear reactors. The power distributions are calculated with the help of nodal expansion methods (NEM) for hexagonal and Cartesian geometry. The fuel rod model and the thermohydraulic part provide fuel temperatures, coolant temperatures and densities as well as boron concentrations for the calculation of feedback effects on the basis of cross section libraries generated by cell codes. Safety relevant parameters like maximum fuel and cladding temperatures, critical heat flux and degree of cladding oxidation are estimated. DYN3D can analyze RIA initiated by moved control rods and/or perturbations of the coolant flow. Stationary and transient boundary conditions for the coolant flow, the core inlet temperatures and boron concentrations at the core inlet have to be given. For analyzing more complex transients the code DYN3D is coupled with the plant model ATHLET of the GRS. The extensive validation work accomplished for DYN3D is presented in several examples. Some applications of the code are described. (orig.) [Deutsch] Die Verwendung 3-dimensionaler Kernmodelle zur Untersuchung der Sicherheitsreserven bei Uebergangsprozessen und Stoerfaellen in Kernreaktoren vermeidet konservative Annahmen, die bei der Benutzung des Punktmodells oder 1-dimensionaler Modelle erforderlich sind. Aus diesen Gruenden wurde das 3-dimensionale Rechenprogramm DYN3D fuer die Untersuchung von Reaktivitaetsstoerfaellen in thermischen Reaktoren entwickelt. Die Leistungsverteilung wird mit nodalen Methoden fuer hexagonale oder kartesische Geometrie berechnet. Das Brennstabmodell und der thermohydraulische Teil von DYN3D liefert die Brennstofftemperaturen, Kuehlmitteltemperaturen

A Validation Approach for Quasistatic Numerical/Experimental Indentation Analysis in Soft Materials Using 3D Digital Image Correlation.

Science.gov (United States)

Felipe-Sesé, Luis; López-Alba, Elías; Hannemann, Benedikt; Schmeer, Sebastian; Diaz, Francisco A

2017-06-28

A quasistatic indentation numerical analysis in a round section specimen made of soft material has been performed and validated with a full field experimental technique, i.e., Digital Image Correlation 3D. The contact experiment specifically consisted of loading a 25 mm diameter rubber cylinder of up to a 5 mm indentation and then unloading. Experimental strains fields measured at the surface of the specimen during the experiment were compared with those obtained by performing two numerical analyses employing two different hyperplastic material models. The comparison was performed using an Image Decomposition new methodology that makes a direct comparison of full-field data independently of their scale or orientation possible. Numerical results show a good level of agreement with those measured during the experiments. However, since image decomposition allows for the differences to be quantified, it was observed that one of the adopted material models reproduces lower differences compared to experimental results.

3D Printing of Molecular Models

Science.gov (United States)

Gardner, Adam; Olson, Arthur

2016-01-01

Physical molecular models have played a valuable role in our understanding of the invisible nano-scale world. We discuss 3D printing and its use in producing models of the molecules of life. Complex biomolecular models, produced from 3D printed parts, can demonstrate characteristics of molecular structure and function, such as viral self-assembly,…

Estimation of thawing cryolithic area with numerical modeling in 3D geometry while exploiting underground small nuclear power plant

Directory of Open Access Journals (Sweden)

Melnikov N. N.

2016-03-01

Full Text Available The paper presents results on 3D numerical calculation of a thermal task related to assessing a thawing area when placing modules with reactor and steam-turbine facility of a small nuclear power plant in thickness of permafrost rocks. The paper discusses influence of the coefficient of thermal conductivity for large-scaled underground excavations lining and cryolithic area porosity on thawing depth and front movement velocity under different spatial directions

NUMERICAL AND MECHANICAL ANALYSES OF A 3D-PRINTED TITANIUM TRABECULAR DENTAL IMPLANT

Directory of Open Access Journals (Sweden)

Luboš Řehounek

2017-06-01

Full Text Available The main focus of this paper is to investigate and describe a novel biomaterial structure. The trabecular structure has only recently been recognized as a viable alternative for prostheses and implants and seems to have very promising biocompatibility and mechanical properties. The 3D printing technique was used to create test specimens. These specimens were then tested by nanoindentation and tensile and compression tests. A numerical model was created and curve-fitted to represent the mechanical behavior of the trabecular structure. A significant reduction in the values of Young’s modulus E was observed. The values of E for conventional implant materials are approximately 110–120GPa and the trabecular structure reached a value just below 1GPa. The next effort will be to apply the model onto a real implant. It is the “four leaf clover” implant variant by authors F. Denk Jr., A. Jíra and F. Denk Sr.

Vascular stents: Coupling full 3-D with reduced-order structural models

International Nuclear Information System (INIS)

Avdeev, I; Shams, M

2010-01-01

Self-expanding nitinol stents are used to treat peripheral arterial disease. The peripheral arteries are subjected to a combination of mechanical forces such as compression, torsion, bending, and contraction. Most commercially available peripheral self-expanding stents are composed of a series of sub-millimeter V-shaped struts, which are laser-cut from a nitinol tube and surface-treated for better fatigue performance. The numerical stent models must accurately predict location and distribution of local stresses and strains caused by large arterial deformations. Full 3-D finite element non-linear analysis of an entire stent is computationally expensive to the point of being prohibitive, especially for longer stents. Reduced-order models based on beam or shell elements are fairly accurate in capturing global deformations, but are not very helpful in predicting stent failure. We propose a mixed approach that combines the full 3-D model and reduced-order models. Several global-local, full 3-D/reduced-order finite element models of a peripheral self-expanding stent were validated and compared with experimental data. The kinematic constraint method used to couple various elements together was found to be very efficient and easily applicable to commercial FEA codes. The proposed mixed models can be used to accurately predict stent failure based on realistic (patient-specific), non-linear kinematic behavior of peripheral arteries.

Virtual 3d City Modeling: Techniques and Applications

Science.gov (United States)

Singh, S. P.; Jain, K.; Mandla, V. R.

2013-08-01

3D city model is a digital representation of the Earth's surface and it's related objects such as Building, Tree, Vegetation, and some manmade feature belonging to urban area. There are various terms used for 3D city models such as "Cybertown", "Cybercity", "Virtual City", or "Digital City". 3D city models are basically a computerized or digital model of a city contains the graphic representation of buildings and other objects in 2.5 or 3D. Generally three main Geomatics approach are using for Virtual 3-D City models generation, in first approach, researcher are using Conventional techniques such as Vector Map data, DEM, Aerial images, second approach are based on High resolution satellite images with LASER scanning, In third method, many researcher are using Terrestrial images by using Close Range Photogrammetry with DSM & Texture mapping. We start this paper from the introduction of various Geomatics techniques for 3D City modeling. These techniques divided in to two main categories: one is based on Automation (Automatic, Semi-automatic and Manual methods), and another is Based on Data input techniques (one is Photogrammetry, another is Laser Techniques). After details study of this, finally in short, we are trying to give the conclusions of this study. In the last, we are trying to give the conclusions of this research paper and also giving a short view for justification and analysis, and present trend for 3D City modeling. This paper gives an overview about the Techniques related with "Generation of Virtual 3-D City models using Geomatics Techniques" and the Applications of Virtual 3D City models. Photogrammetry, (Close range, Aerial, Satellite), Lasergrammetry, GPS, or combination of these modern Geomatics techniques play a major role to create a virtual 3-D City model. Each and every techniques and method has some advantages and some drawbacks. Point cloud model is a modern trend for virtual 3-D city model. Photo-realistic, Scalable, Geo-referenced virtual 3

3D mathematical modelling of scour around a circular pile in current

DEFF Research Database (Denmark)

Roulund, Andreas; Sumer, B. Mutlu; Fredsøe, Jørgen

1999-01-01

This paper deals with scour around a circular pile exposed to a steady current. A 3D numerical model incorporated with the k-w,SST closure coupled with the sediment-continuity equation and a bedload sediment transport formula has been used to predict the scour. 3D calculations have also been...... carried out for a plane rigid bottom for reference purpose. The predicted flow features are apparently in fairly good agreement with the experimental data. Early calculations indicate that the model is able to predict the scour properties satisfactorily in the initial stages of the scour process, up...... to scour depth of 0.6-0.7 times the pile diameter. Calculations that describe the entire scour process (including the equilibrium stage) are underway....

The Esri 3D city information model

International Nuclear Information System (INIS)

Reitz, T; Schubiger-Banz, S

2014-01-01

With residential and commercial space becoming increasingly scarce, cities are going vertical. Managing the urban environments in 3D is an increasingly important and complex undertaking. To help solving this problem, Esri has released the ArcGIS for 3D Cities solution. The ArcGIS for 3D Cities solution provides the information model, tools and apps for creating, analyzing and maintaining a 3D city using the ArcGIS platform. This paper presents an overview of the 3D City Information Model and some sample use cases

Application of an unstructured 3D finite volume numerical model to flows and salinity dynamics in the San Francisco Bay-Delta

Science.gov (United States)

Martyr-Koller, R.C.; Kernkamp, H.W.J.; Van Dam, Anne A.; Mick van der Wegen,; Lucas, Lisa; Knowles, N.; Jaffe, B.; Fregoso, T.A.

2017-01-01

A linked modeling approach has been undertaken to understand the impacts of climate and infrastructure on aquatic ecology and water quality in the San Francisco Bay-Delta region. The Delft3D Flexible Mesh modeling suite is used in this effort for its 3D hydrodynamics, salinity, temperature and sediment dynamics, phytoplankton and water-quality coupling infrastructure, and linkage to a habitat suitability model. The hydrodynamic model component of the suite is D-Flow FM, a new 3D unstructured finite-volume model based on the Delft3D model. In this paper, D-Flow FM is applied to the San Francisco Bay-Delta to investigate tidal, seasonal and annual dynamics of water levels, river flows and salinity under historical environmental and infrastructural conditions. The model is driven by historical winds, tides, ocean salinity, and river flows, and includes federal, state, and local freshwater withdrawals, and regional gate and barrier operations. The model is calibrated over a 9-month period, and subsequently validated for water levels, flows, and 3D salinity dynamics over a 2 year period.Model performance was quantified using several model assessment metrics and visualized through target diagrams. These metrics indicate that the model accurately estimated water levels, flows, and salinity over wide-ranging tidal and fluvial conditions, and the model can be used to investigate detailed circulation and salinity patterns throughout the Bay-Delta. The hydrodynamics produced through this effort will be used to drive affiliated sediment, phytoplankton, and contaminant hindcast efforts and habitat suitability assessments for fish and bivalves. The modeling framework applied here will serve as a baseline to ultimately shed light on potential ecosystem change over the current century.

Application of an unstructured 3D finite volume numerical model to flows and salinity dynamics in the San Francisco Bay-Delta

Science.gov (United States)

Martyr-Koller, R. C.; Kernkamp, H. W. J.; van Dam, A.; van der Wegen, M.; Lucas, L. V.; Knowles, N.; Jaffe, B.; Fregoso, T. A.

2017-06-01

A linked modeling approach has been undertaken to understand the impacts of climate and infrastructure on aquatic ecology and water quality in the San Francisco Bay-Delta region. The Delft3D Flexible Mesh modeling suite is used in this effort for its 3D hydrodynamics, salinity, temperature and sediment dynamics, phytoplankton and water-quality coupling infrastructure, and linkage to a habitat suitability model. The hydrodynamic model component of the suite is D-Flow FM, a new 3D unstructured finite-volume model based on the Delft3D model. In this paper, D-Flow FM is applied to the San Francisco Bay-Delta to investigate tidal, seasonal and annual dynamics of water levels, river flows and salinity under historical environmental and infrastructural conditions. The model is driven by historical winds, tides, ocean salinity, and river flows, and includes federal, state, and local freshwater withdrawals, and regional gate and barrier operations. The model is calibrated over a 9-month period, and subsequently validated for water levels, flows, and 3D salinity dynamics over a 2 year period. Model performance was quantified using several model assessment metrics and visualized through target diagrams. These metrics indicate that the model accurately estimated water levels, flows, and salinity over wide-ranging tidal and fluvial conditions, and the model can be used to investigate detailed circulation and salinity patterns throughout the Bay-Delta. The hydrodynamics produced through this effort will be used to drive affiliated sediment, phytoplankton, and contaminant hindcast efforts and habitat suitability assessments for fish and bivalves. The modeling framework applied here will serve as a baseline to ultimately shed light on potential ecosystem change over the current century.

3D Printing device adaptable to Computer Numerical Control (CNC)

OpenAIRE

GARDAN , Julien; Danesi , F.; Roucoules , Lionel; Schneider , A.

2014-01-01

This article presents the development of a 3D printing device for the additive manufacturing adapted to a CNC machining. The application involves the integration of a specific printing head. Additive manufacturing technology is most commonly used for modeling, prototyping, tooling through an exclusive machine or 3D printer. A global review and analysis of technologies show the additive manufacturing presents little independent solutions [6][9]. The problem studied especially the additive manu...

Using high-order polynomial basis in 3-D EM forward modeling based on volume integral equation method

Science.gov (United States)

Kruglyakov, Mikhail; Kuvshinov, Alexey

2018-05-01

3-D interpretation of electromagnetic (EM) data of different origin and scale becomes a common practice worldwide. However, 3-D EM numerical simulations (modeling)—a key part of any 3-D EM data analysis—with realistic levels of complexity, accuracy and spatial detail still remains challenging from the computational point of view. We present a novel, efficient 3-D numerical solver based on a volume integral equation (IE) method. The efficiency is achieved by using a high-order polynomial (HOP) basis instead of the zero-order (piecewise constant) basis that is invoked in all routinely used IE-based solvers. We demonstrate that usage of the HOP basis allows us to decrease substantially the number of unknowns (preserving the same accuracy), with corresponding speed increase and memory saving.

Development of Numerical Grids for UZ Flow and Transport Modeling

International Nuclear Information System (INIS)

P. Dobson

2003-01-01

This Scientific Analysis report describes the methods used to develop numerical grids of the unsaturated hydrogeologic system beneath Yucca Mountain. Numerical grid generation is an integral part of the development of the Unsaturated Zone Flow and Transport Model (UZ Model), a complex, three-dimensional (3-D) model of Yucca Mountain. This revision incorporates changes made to both the geologic framework model and the proposed repository layout. The resulting numerical grids, developed using current geologic, hydrogeologic, and mineralogic data, provide the necessary framework to: (1) develop calibrated hydrogeologic property sets and flow fields, (2) test conceptual hypotheses of flow and transport, and (3) predict flow and transport behavior under a variety of climatic and thermal-loading conditions. The technical scope, content, and management of this Scientific Analysis report was initially controlled by the planning document, ''Technical Work Plan (TWP) for: Unsaturated Zone Sections of License Application Chapters 8 and 12'' (BSC 2002 [159051], Section 1.6.4). This TWP was later superseded by ''Technical Work Plan for: Performance Assessment Unsaturated Zone'' (BSC 2002 [160819]), which contains the Data Qualification Plan used to qualify the DTN: MO0212GWLSSPAX.000 [161271] (See Attachment IV). Grids generated and documented in this report supersede those documented in previous versions of this report (BSC 2001 [159356]). The constraints, assumptions, and limitations associated with this report are discussed in the appropriate sections that follow. There were no deviations from the TWP scope of work in this report. Two software packages not listed in Table IV-2 of the TWP (BSC 2002 [159051]), ARCINFO V7.2.1 (CRWMS M and O 2000 [157019]; USGS 2000 [148304]) and 2kgrid8.for V1.0 (LBNL 2002 [154787]), were utilized in the development of the numerical grids; the use of additional software is accounted for in the TWP (BSC 2002 [159051], Section 13). The use of

Mechano-sensing and cell migration: a 3D model approach

International Nuclear Information System (INIS)

Borau, C; García-Aznar, J M; Kamm, R D

2011-01-01

Cell migration is essential for tissue development in different physiological and pathological conditions. It is a complex process orchestrated by chemistry, biological factors, microstructure and surrounding mechanical properties. Focusing on the mechanical interactions, cells do not only exert forces on the matrix that surrounds them, but they also sense and react to mechanical cues in a process called mechano-sensing. Here, we hypothesize the involvement of mechano-sensing in the regulation of directional cell migration through a three-dimensional (3D) matrix. For this purpose, we develop a 3D numerical model of individual cell migration, which incorporates the mechano-sensing process of the cell as the main mechanism regulating its movement. Consistent with this hypothesis, we found that factors, such as substrate stiffness, boundary conditions and external forces, regulate specific and distinct cell movements

3D-modeling and 3D-printing explorations on Japanese tea ceremony utensils

NARCIS (Netherlands)

Levy, P.D.; Yamada, Shigeru

2017-01-01

In this paper, we inquire aesthetical aspects of the Japanese tea ceremony, described as the aesthetics in the imperfection, based on novel fabrication technologies: 3D-modeling and 3D-printing. To do so, 3D-printed utensils (chashaku and chasen) were iteratively designed for the ceremony and were

Realistic deformable 3D numeric phantom for transcutaneous ultrasound

Energy Technology Data Exchange (ETDEWEB)

Cardoso, Fernando Mitsuyama; Moraes, Matheus Cardoso; Furuie, Sergio Shiguemi, E-mail: fernando.okara@gmail.com [Universidade de Sao Paulo (USP), SP (Brazil). Escola de Engenharia

2017-01-15

Introduction: Numerical phantoms are important tools to design, calibrate and evaluate several methods in various image-processing applications, such as echocardiography and mammography. We present a framework for creating ultrasound numerical deformable phantoms based on Finite Element Method (FEM), Linear Isomorphism and Field II. The proposed method considers that the scatterers map is a property of the tissue; therefore, the scatterers should move according to the tissue strain. Methods: First, a volume representing the target tissue is loaded. Second, parameter values, such as Young's Modulus, scatterers density, attenuation and scattering amplitudes are inserted for each different regions of the phantom. Then, other parameters related to the ultrasound equipment, such as ultrasound frequency and number of transducer elements, are also defined in order to perform the ultrasound acquisition using Field II. Third, the size and position of the transducer and the pressures that are applied against the tissue are defined. Subsequently, FEM is executed and deformation is computed. Next, 3D linear isomorphism is performed to displace the scatterers according to the deformation. Finally, Field II is carried out to generate the non-deformed and deformed ultrasound data. Results: The framework is evaluated by comparing strain values obtained the numerical simulation and from the physical phantom from CIRS. The mean difference between both phantoms is lesser than 10%. Conclusion: The acoustic and deformation outcomes are similar to those obtained using a physical phantom. This framework led to a tool, which is available online and free of charges for educational and research purposes. (author)

A Hybrid 3D Indoor Space Model

Directory of Open Access Journals (Sweden)

A. Jamali

2016-10-01

Full Text Available GIS integrates spatial information and spatial analysis. An important example of such integration is for emergency response which requires route planning inside and outside of a building. Route planning requires detailed information related to indoor and outdoor environment. Indoor navigation network models including Geometric Network Model (GNM, Navigable Space Model, sub-division model and regular-grid model lack indoor data sources and abstraction methods. In this paper, a hybrid indoor space model is proposed. In the proposed method, 3D modeling of indoor navigation network is based on surveying control points and it is less dependent on the 3D geometrical building model. This research proposes a method of indoor space modeling for the buildings which do not have proper 2D/3D geometrical models or they lack semantic or topological information. The proposed hybrid model consists of topological, geometrical and semantical space.

BEAMS3D Neutral Beam Injection Model

Energy Technology Data Exchange (ETDEWEB)

Lazerson, Samuel

2014-04-14

With the advent of applied 3D fi elds in Tokamaks and modern high performance stellarators, a need has arisen to address non-axisymmetric effects on neutral beam heating and fueling. We report on the development of a fully 3D neutral beam injection (NBI) model, BEAMS3D, which addresses this need by coupling 3D equilibria to a guiding center code capable of modeling neutral and charged particle trajectories across the separatrix and into the plasma core. Ionization, neutralization, charge-exchange, viscous velocity reduction, and pitch angle scattering are modeled with the ADAS atomic physics database [1]. Benchmark calculations are presented to validate the collisionless particle orbits, neutral beam injection model, frictional drag, and pitch angle scattering effects. A calculation of neutral beam heating in the NCSX device is performed, highlighting the capability of the code to handle 3D magnetic fields.

Finite-Source Inversion for the 2004 Parkfield Earthquake using 3D Velocity Model Green's Functions

Science.gov (United States)

Kim, A.; Dreger, D.; Larsen, S.

2008-12-01

We determine finite fault models of the 2004 Parkfield earthquake using 3D Green's functions. Because of the dense station coverage and detailed 3D velocity structure model in this region, this earthquake provides an excellent opportunity to examine how the 3D velocity structure affects the finite fault inverse solutions. Various studies (e.g. Michaels and Eberhart-Phillips, 1991; Thurber et al., 2006) indicate that there is a pronounced velocity contrast across the San Andreas Fault along the Parkfield segment. Also the fault zone at Parkfield is wide as evidenced by mapped surface faults and where surface slip and creep occurred in the 1966 and the 2004 Parkfield earthquakes. For high resolution images of the rupture process"Ait is necessary to include the accurate 3D velocity structure for the finite source inversion. Liu and Aurchuleta (2004) performed finite fault inversions using both 1D and 3D Green's functions for 1989 Loma Prieta earthquake using the same source paramerization and data but different Green's functions and found that the models were quite different. This indicates that the choice of the velocity model significantly affects the waveform modeling at near-fault stations. In this study, we used the P-wave velocity model developed by Thurber et al (2006) to construct the 3D Green's functions. P-wave speeds are converted to S-wave speeds and density using by the empirical relationships of Brocher (2005). Using a finite difference method, E3D (Larsen and Schultz, 1995), we computed the 3D Green's functions numerically by inserting body forces at each station. Using reciprocity, these Green's functions are recombined to represent the ground motion at each station due to the slip on the fault plane. First we modeled the waveforms of small earthquakes to validate the 3D velocity model and the reciprocity of the Green"fs function. In the numerical tests we found that the 3D velocity model predicted the individual phases well at frequencies lower than 0

Creating safer coastal and port infrastructure with innovative physical and numerical modelling

CSIR Research Space (South Africa)

Tulsi, K

2015-10-01

Full Text Available Infrastructure with Innovative Physical and Numerical Modelling Kishan Tulsi  Physical and Numerical modelling  Breakwater Monitoring  Armour track  Vessel manoeuvring simulations for safe port design and operations  Simflex software... – Integrated Port Operations Support System  Virtual Buoy Physical modelling Numerical modelling Armour Track Armour Track Armour Track Armour Track Armour track using 3D data points Ship manoeuvring simulations: Ship Manoeuvring simulations Port...

3D numerical model of the spherical particle saltation in a channel with a rough fixed bed

Czech Academy of Sciences Publication Activity Database

Lukerchenko, Nikolay; Piatsevich, Siarhei; Chára, Zdeněk; Vlasák, Pavel

2009-01-01

Roč. 57, č. 2 (2009), s. 100-112 ISSN 0042-790X R&D Projects: GA ČR GA103/06/1487 Institutional research plan: CEZ:AV0Z20600510 Keywords : 3D Saltation Model * Bed-Load Transport * Particle-Bed Collision * Particle Rotation * Particle Lateral Dispersion Subject RIV: BK - Fluid Dynamics Impact factor: 1.000, year: 2009

DKE: a fast numerical solver for the 3-D relativistic bounce-averaged electron drift kinetic equation

Energy Technology Data Exchange (ETDEWEB)

Decker, J.; Peysson, Y

2004-12-01

A new original code for solving the 3-D relativistic and bounce-averaged electron drift kinetic equation is presented. It designed for the current drive problem in tokamak with an arbitrary magnetic equilibrium. This tool allows self-consistent calculations of the bootstrap current in presence of other external current sources. RF current drive for arbitrary type of waves may be used. Several moments of the electron distribution function are determined, like the exact and effective fractions of trapped electrons, the plasma current, absorbed RF power, runaway and magnetic ripple loss rates and non-thermal Bremsstrahlung. Advanced numerical techniques have been used to make it the first fully implicit (reverse time) 3-D solver, particularly well designed for implementation in a chain of code for realistic current drive calculations in high {beta}{sub p} plasmas. All the details of the physics background and the numerical scheme are presented, as well a some examples to illustrate main code capabilities. Several important numerical points are addressed concerning code stability and potential numerical and physical limitations. (authors)

DKE: a fast numerical solver for the 3-D relativistic bounce-averaged electron drift kinetic equation

International Nuclear Information System (INIS)

Decker, J.; Peysson, Y.

2004-12-01

A new original code for solving the 3-D relativistic and bounce-averaged electron drift kinetic equation is presented. It designed for the current drive problem in tokamak with an arbitrary magnetic equilibrium. This tool allows self-consistent calculations of the bootstrap current in presence of other external current sources. RF current drive for arbitrary type of waves may be used. Several moments of the electron distribution function are determined, like the exact and effective fractions of trapped electrons, the plasma current, absorbed RF power, runaway and magnetic ripple loss rates and non-thermal Bremsstrahlung. Advanced numerical techniques have been used to make it the first fully implicit (reverse time) 3-D solver, particularly well designed for implementation in a chain of code for realistic current drive calculations in high β p plasmas. All the details of the physics background and the numerical scheme are presented, as well a some examples to illustrate main code capabilities. Several important numerical points are addressed concerning code stability and potential numerical and physical limitations. (authors)

Development of a stereolithography (STL input and computer numerical control (CNC output algorithm for an entry-level 3-D printer

Directory of Open Access Journals (Sweden)

Brown, Andrew

2014-08-01

Full Text Available This paper presents a prototype Stereolithography (STL file format slicing and tool-path generation algorithm, which serves as a data front-end for a Rapid Prototyping (RP entry- level three-dimensional (3-D printer. Used mainly in Additive Manufacturing (AM, 3-D printers are devices that apply plastic, ceramic, and metal, layer by layer, in all three dimensions on a flat surface (X, Y, and Z axis. 3-D printers, unfortunately, cannot print an object without a special algorithm that is required to create the Computer Numerical Control (CNC instructions for printing. An STL algorithm therefore forms a critical component for Layered Manufacturing (LM, also referred to as RP. The purpose of this study was to develop an algorithm that is capable of processing and slicing an STL file or multiple files, resulting in a tool-path, and finally compiling a CNC file for an entry-level 3- D printer. The prototype algorithm was implemented for an entry-level 3-D printer that utilises the Fused Deposition Modelling (FDM process or Solid Freeform Fabrication (SFF process; an AM technology. Following an experimental method, the full data flow path for the prototype algorithm was developed, starting with STL data files, and then processing the STL data file into a G-code file format by slicing the model and creating a tool-path. This layering method is used by most 3-D printers to turn a 2-D object into a 3-D object. The STL algorithm developed in this study presents innovative opportunities for LM, since it allows engineers and architects to transform their ideas easily into a solid model in a fast, simple, and cheap way. This is accomplished by allowing STL models to be sliced rapidly, effectively, and without error, and finally to be processed and prepared into a G-code print file.

Validation of a 2-D semi-coupled numerical model for fluid-structure-seabed interaction

Science.gov (United States)

Ye, Jianhong; Jeng, Dongsheng; Wang, Ren; Zhu, Changqi

2013-10-01

A 2-D semi-coupled model PORO-WSSI 2D (also be referred as FSSI-CAS 2D) for the Fluid-Structure-Seabed Interaction (FSSI) has been developed by employing RANS equations for wave motion in fluid domain, VARANS equations for porous flow in porous structures; and taking the dynamic Biot's equations (known as "u - p" approximation) for soil as the governing equations. The finite difference two-step projection method and the forward time difference method are adopted to solve the RANS, VARANS equations; and the finite element method is adopted to solve the "u - p" approximation. A data exchange port is developed to couple the RANS, VARANS equations and the dynamic Biot's equations together. The analytical solution proposed by Hsu and Jeng (1994) and some experiments conducted in wave flume or geotechnical centrifuge in which various waves involved are used to validate the developed semi-coupled numerical model. The sandy bed involved in these experiments is poro-elastic or poro-elastoplastic. The inclusion of the interaction between fluid, marine structures and poro-elastoplastic seabed foundation is a special point and highlight in this paper, which is essentially different with other previous coupled models The excellent agreement between the numerical results and the experiment data indicates that the developed coupled model is highly reliablefor the FSSI problem.

Numerical simulation of runaway electrons: 3-D effects on synchrotron radiation and impurity-based runaway current dissipation

Science.gov (United States)

del-Castillo-Negrete, D.; Carbajal, L.; Spong, D.; Izzo, V.

2018-05-01

Numerical simulations of runaway electrons (REs) with a particular emphasis on orbit dependent effects in 3-D magnetic fields are presented. The simulations were performed using the recently developed Kinetic Orbit Runaway electron Code (KORC) that computes the full-orbit relativistic dynamics in prescribed electric and magnetic fields including radiation damping and collisions. The two main problems of interest are synchrotron radiation and impurity-based RE dissipation. Synchrotron radiation is studied in axisymmetric fields and in 3-D magnetic configurations exhibiting magnetic islands and stochasticity. For passing particles in axisymmetric fields, neglecting orbit effects might underestimate or overestimate the total radiation power depending on the direction of the radial shift of the drift orbits. For trapped particles, the spatial distribution of synchrotron radiation exhibits localized "hot" spots at the tips of the banana orbits. In general, the radiation power per particle for trapped particles is higher than the power emitted by passing particles. The spatial distribution of synchrotron radiation in stochastic magnetic fields, obtained using the MHD code NIMROD, is strongly influenced by the presence of magnetic islands. 3-D magnetic fields also introduce a toroidal dependence on the synchrotron spectra, and neglecting orbit effects underestimates the total radiation power. In the presence of magnetic islands, the radiation damping of trapped particles is larger than the radiation damping of passing particles. Results modeling synchrotron emission by RE in DIII-D quiescent plasmas are also presented. The computation uses EFIT reconstructed magnetic fields and RE energy distributions fitted to the experimental measurements. Qualitative agreement is observed between the numerical simulations and the experiments for simplified RE pitch angle distributions. However, it is noted that to achieve quantitative agreement, it is necessary to use pitch angle

3D numerical simulation on heat transfer performance of a cylindrical liquid immersion solar receiver

International Nuclear Information System (INIS)

Xiang Haijun; Wang Yiping; Zhu Li; Han Xinyue; Sun Yong; Zhao Zhengjian

2012-01-01

Highlights: ► Establishment of a three-dimensional numerical simulation model of a cylindrical liquid immersion solar receiver. ► Determination of model parameters and validation of the model by using the real-collected data. ► Optimization of liquid flow rate and fin’s structure for better heat transfer performance. - Abstract: Liquid immersion cooling for a cylindrical solar receiver in a dish concentrator photovoltaic system has been experimentally verified to be a promising method of removing surplus heat from densely packed solar cells. In the present study, a three-dimensional (3D) numerical simulation model of the prototype was established for better understanding the mechanism of the direct-contact heat transfer process. With the selection of standard k–ε turbulent model, the detailed simulation results of velocity field and temperature characteristics were obtained. The heat transfer performance of two structural modules (bare module and finned module) under actual weather conditions was simulated. It was found that the predicted temperature distribution of the two structural modules at the axial and lateral direction was in good agreement with the experimental data. Based on the validated simulation model, the influence of liquid flow rate and module geometric parameters on the cell temperature was then investigated. The simulated results indicated that the cell module with fin height of 4 mm and fin number of 11 has the best heat transfer performance and will be used in further works.

Image based 3D city modeling : Comparative study

Directory of Open Access Journals (Sweden)

S. P. Singh

2014-06-01

Full Text Available 3D city model is a digital representation of the Earth’s surface and it’s related objects such as building, tree, vegetation, and some manmade feature belonging to urban area. The demand of 3D city modeling is increasing rapidly for various engineering and non-engineering applications. Generally four main image based approaches were used for virtual 3D city models generation. In first approach, researchers were used Sketch based modeling, second method is Procedural grammar based modeling, third approach is Close range photogrammetry based modeling and fourth approach is mainly based on Computer Vision techniques. SketchUp, CityEngine, Photomodeler and Agisoft Photoscan are the main softwares to represent these approaches respectively. These softwares have different approaches & methods suitable for image based 3D city modeling. Literature study shows that till date, there is no complete such type of comparative study available to create complete 3D city model by using images. This paper gives a comparative assessment of these four image based 3D modeling approaches. This comparative study is mainly based on data acquisition methods, data processing techniques and output 3D model products. For this research work, study area is the campus of civil engineering department, Indian Institute of Technology, Roorkee (India. This 3D campus acts as a prototype for city. This study also explains various governing parameters, factors and work experiences. This research work also gives a brief introduction, strengths and weakness of these four image based techniques. Some personal comment is also given as what can do or what can’t do from these softwares. At the last, this study shows; it concluded that, each and every software has some advantages and limitations. Choice of software depends on user requirements of 3D project. For normal visualization project, SketchUp software is a good option. For 3D documentation record, Photomodeler gives good

3D Numerical Investigation of the Flow Pattern in Settling Basins by SSIIM2 Software

Directory of Open Access Journals (Sweden)

Edris Merufinia

2017-01-01

Full Text Available Settling basins are considered as a major and important component in conventional water treatment processes. The high cost of constructing these basins accounts for approximately 30% of the total capital investment of water treatment plants. Hence, it is essential to model and optimize their performance beforehand. In settling and sedimentation basins, different areas of flow including secondary and rotational flows occur due to the velocity gradient. Such phenomena cause short paths, increase flow stationary and dead zones, and change the mixing rates of the flow, which collectively prevent laminar conditions to be created for the sedimentation process, and thereby reduce process efficiency. The remedy is to reduce as far as possible the dead zones in the flow. The first step to optimize a settling basin is to calculate accurately the velocity field and the volume of rotation zones. The present study presents the numerical simulation of a flow in a rectangular basin. Continuity and Navier-Stokes equations are solved using finite volume method. A 3D flow simulation is performed using the standard k-ε turbulence model for settling basins with and without baffles. Finally, the numerical results obtained are compared with experimental results reported elsewhere.

A 3D Printing Model Watermarking Algorithm Based on 3D Slicing and Feature Points

Directory of Open Access Journals (Sweden)

Giao N. Pham

2018-02-01

Full Text Available With the increase of three-dimensional (3D printing applications in many areas of life, a large amount of 3D printing data is copied, shared, and used several times without any permission from the original providers. Therefore, copyright protection and ownership identification for 3D printing data in communications or commercial transactions are practical issues. This paper presents a novel watermarking algorithm for 3D printing models based on embedding watermark data into the feature points of a 3D printing model. Feature points are determined and computed by the 3D slicing process along the Z axis of a 3D printing model. The watermark data is embedded into a feature point of a 3D printing model by changing the vector length of the feature point in OXY space based on the reference length. The x and y coordinates of the feature point will be then changed according to the changed vector length that has been embedded with a watermark. Experimental results verified that the proposed algorithm is invisible and robust to geometric attacks, such as rotation, scaling, and translation. The proposed algorithm provides a better method than the conventional works, and the accuracy of the proposed algorithm is much higher than previous methods.

From 2D to 3D: Using Illumination Cones to Build 3d Face Model

International Nuclear Information System (INIS)

Xiao, S S; Jin, M

2006-01-01

To solve the problem derivate by lighting condition and position of the camera, a new method using illumination cones to build 3d face model has been proposed. Due to illumination variability, the same object can show dramatic difference even as being viewed in fixed pose. To handle this variability, an object recognition system must employ a representation that is either invariant to, or can model this variability. The proposed technique presents an appearance-based method for modeling the variability due to illumination in the images of objects. The method differs from past appearance-based methods. Evenmore, a small set of training images is used to generate a representation that the illumination cone models the complete set of images of an object with Lambertian reflectance surface under a combination of arbitrary point light sources at infinity. After building up the illumination cones, researches focus on how to present the 3d model of the face. Combining illumination and texture feature to build up 3d model of the face make it easy solving the problem in recognition of face under different pose

An Embedded 3D Fracture Modeling Approach for Simulating Fracture-Dominated Fluid Flow and Heat Transfer in Geothermal Reservoirs

Energy Technology Data Exchange (ETDEWEB)

Johnston, Henry [National Renewable Energy Laboratory (NREL), Golden, CO (United States); Wang, Cong [Colorado School of Mines; Winterfeld, Philip [Colorado School of Mines; Wu, Yu-Shu [Colorado School of Mines

2018-02-14

An efficient modeling approach is described for incorporating arbitrary 3D, discrete fractures, such as hydraulic fractures or faults, into modeling fracture-dominated fluid flow and heat transfer in fractured geothermal reservoirs. This technique allows 3D discrete fractures to be discretized independently from surrounding rock volume and inserted explicitly into a primary fracture/matrix grid, generated without including 3D discrete fractures in prior. An effective computational algorithm is developed to discretize these 3D discrete fractures and construct local connections between 3D fractures and fracture/matrix grid blocks of representing the surrounding rock volume. The constructed gridding information on 3D fractures is then added to the primary grid. This embedded fracture modeling approach can be directly implemented into a developed geothermal reservoir simulator via the integral finite difference (IFD) method or with TOUGH2 technology This embedded fracture modeling approach is very promising and computationally efficient to handle realistic 3D discrete fractures with complicated geometries, connections, and spatial distributions. Compared with other fracture modeling approaches, it avoids cumbersome 3D unstructured, local refining procedures, and increases computational efficiency by simplifying Jacobian matrix size and sparsity, while keeps sufficient accuracy. Several numeral simulations are present to demonstrate the utility and robustness of the proposed technique. Our numerical experiments show that this approach captures all the key patterns about fluid flow and heat transfer dominated by fractures in these cases. Thus, this approach is readily available to simulation of fractured geothermal reservoirs with both artificial and natural fractures.

Modeling variably saturated multispecies reactive groundwater solute transport with MODFLOW-UZF and RT3D

Science.gov (United States)

Bailey, Ryan T.; Morway, Eric D.; Niswonger, Richard G.; Gates, Timothy K.

2013-01-01

A numerical model was developed that is capable of simulating multispecies reactive solute transport in variably saturated porous media. This model consists of a modified version of the reactive transport model RT3D (Reactive Transport in 3 Dimensions) that is linked to the Unsaturated-Zone Flow (UZF1) package and MODFLOW. Referred to as UZF-RT3D, the model is tested against published analytical benchmarks as well as other published contaminant transport models, including HYDRUS-1D, VS2DT, and SUTRA, and the coupled flow and transport modeling system of CATHY and TRAN3D. Comparisons in one-dimensional, two-dimensional, and three-dimensional variably saturated systems are explored. While several test cases are included to verify the correct implementation of variably saturated transport in UZF-RT3D, other cases are included to demonstrate the usefulness of the code in terms of model run-time and handling the reaction kinetics of multiple interacting species in variably saturated subsurface systems. As UZF1 relies on a kinematic-wave approximation for unsaturated flow that neglects the diffusive terms in Richards equation, UZF-RT3D can be used for large-scale aquifer systems for which the UZF1 formulation is reasonable, that is, capillary-pressure gradients can be neglected and soil parameters can be treated as homogeneous. Decreased model run-time and the ability to include site-specific chemical species and chemical reactions make UZF-RT3D an attractive model for efficient simulation of multispecies reactive transport in variably saturated large-scale subsurface systems.

Improving the seismic small-scale modelling by comparison with numerical methods

Science.gov (United States)

Pageot, Damien; Leparoux, Donatienne; Le Feuvre, Mathieu; Durand, Olivier; Côte, Philippe; Capdeville, Yann

2017-10-01

The potential of experimental seismic modelling at reduced scale provides an intermediate step between numerical tests and geophysical campaigns on field sites. Recent technologies such as laser interferometers offer the opportunity to get data without any coupling effects. This kind of device is used in the Mesures Ultrasonores Sans Contact (MUSC) measurement bench for which an automated support system makes possible to generate multisource and multireceivers seismic data at laboratory scale. Experimental seismic modelling would become a great tool providing a value-added stage in the imaging process validation if (1) the experimental measurement chain is perfectly mastered, and thus if the experimental data are perfectly reproducible with a numerical tool, as well as if (2) the effective source is reproducible along the measurement setup. These aspects for a quantitative validation concerning devices with piezoelectrical sources and a laser interferometer have not been yet quantitatively studied in published studies. Thus, as a new stage for the experimental modelling approach, these two key issues are tackled in the proposed paper in order to precisely define the quality of the experimental small-scale data provided by the bench MUSC, which are available in the scientific community. These two steps of quantitative validation are dealt apart any imaging techniques in order to offer the opportunity to geophysicists who want to use such data (delivered as free data) of precisely knowing their quality before testing any imaging technique. First, in order to overcome the 2-D-3-D correction usually done in seismic processing when comparing 2-D numerical data with 3-D experimental measurement, we quantitatively refined the comparison between numerical and experimental data by generating accurate experimental line sources, avoiding the necessity of geometrical spreading correction for 3-D point-source data. The comparison with 2-D and 3-D numerical modelling is based on

Aerosol numerical modelling at local scale

International Nuclear Information System (INIS)

Albriet, Bastien

2007-01-01

At local scale and in urban areas, an important part of particulate pollution is due to traffic. It contributes largely to the high number concentrations observed. Two aerosol sources are mainly linked to traffic. Primary emission of soot particles and secondary nanoparticle formation by nucleation. The emissions and mechanisms leading to the formation of such bimodal distribution are still badly understood nowadays. In this thesis, we try to provide an answer to this problematic by numerical modelling. The Modal Aerosol Model MAM is used, coupled with two 3D-codes: a CFD (Mercure Saturne) and a CTM (Polair3D). A sensitivity analysis is performed, at the border of a road but also in the first meters of an exhaust plume, to identify the role of each process involved and the sensitivity of different parameters used in the modelling. (author) [fr

Numerical Simulation of 3D Solid-Liquid Turbulent Flow in a Low Specific Speed Centrifugal Pump: Flow Field Analysis

Directory of Open Access Journals (Sweden)

Baocheng Shi

2014-06-01

Full Text Available For numerically simulating 3D solid-liquid turbulent flow in low specific speed centrifugal pumps, the iteration convergence problem caused by complex internal structure and high rotational speed of pump is always a problem for numeral simulation researchers. To solve this problem, the combination of three measures of dynamic underrelaxation factor adjustment, step method, and rotational velocity control means according to residual curves trends of operating parameters was used to improve the numerical convergence. Numeral simulation of 3D turbulent flow in a low specific speed solid-liquid centrifugal pump was performed, and the results showed that the improved solution strategy is greatly helpful to the numerical convergence. Moreover, the 3D turbulent flow fields in pumps have been simulated for the bottom ash-particles with the volume fraction of 10%, 20%, and 30% at the same particle diameter of 0.1 mm. The two-phase calculation results are compared with those of single-phase clean water flow. The calculated results gave the main region of the abrasion of the impeller and volute casing and improve the hydraulic design of the impeller in order to decrease the abrasion and increase the service life of the pump.

1D and 2D Numerical Modeling for Solving Dam-Break Flow Problems Using Finite Volume Method

Directory of Open Access Journals (Sweden)

Szu-Hsien Peng

2012-01-01

Full Text Available The purpose of this study is to model the flow movement in an idealized dam-break configuration. One-dimensional and two-dimensional motion of a shallow flow over a rigid inclined bed is considered. The resulting shallow water equations are solved by finite volumes using the Roe and HLL schemes. At first, the one-dimensional model is considered in the development process. With conservative finite volume method, splitting is applied to manage the combination of hyperbolic term and source term of the shallow water equation and then to promote 1D to 2D. The simulations are validated by the comparison with flume experiments. Unsteady dam-break flow movement is found to be reasonably well captured by the model. The proposed concept could be further developed to the numerical calculation of non-Newtonian fluid or multilayers fluid flow.

Improving Semantic Updating Method on 3d City Models Using Hybrid Semantic-Geometric 3d Segmentation Technique

Science.gov (United States)

Sharkawi, K.-H.; Abdul-Rahman, A.

2013-09-01

Cities and urban areas entities such as building structures are becoming more complex as the modern human civilizations continue to evolve. The ability to plan and manage every territory especially the urban areas is very important to every government in the world. Planning and managing cities and urban areas based on printed maps and 2D data are getting insufficient and inefficient to cope with the complexity of the new developments in big cities. The emergence of 3D city models have boosted the efficiency in analysing and managing urban areas as the 3D data are proven to represent the real world object more accurately. It has since been adopted as the new trend in buildings and urban management and planning applications. Nowadays, many countries around the world have been generating virtual 3D representation of their major cities. The growing interest in improving the usability of 3D city models has resulted in the development of various tools for analysis based on the 3D city models. Today, 3D city models are generated for various purposes such as for tourism, location-based services, disaster management and urban planning. Meanwhile, modelling 3D objects are getting easier with the emergence of the user-friendly tools for 3D modelling available in the market. Generating 3D buildings with high accuracy also has become easier with the availability of airborne Lidar and terrestrial laser scanning equipments. The availability and accessibility to this technology makes it more sensible to analyse buildings in urban areas using 3D data as it accurately represent the real world objects. The Open Geospatial Consortium (OGC) has accepted CityGML specifications as one of the international standards for representing and exchanging spatial data, making it easier to visualize, store and manage 3D city models data efficiently. CityGML able to represents the semantics, geometry, topology and appearance of 3D city models in five well-defined Level-of-Details (LoD), namely LoD0

Modelling river bank erosion processes and mass failure mechanisms using 2-D depth averaged numerical model

Science.gov (United States)

Die Moran, Andres; El kadi Abderrezzak, Kamal; Tassi, Pablo; Herouvet, Jean-Michel

2014-05-01

Bank erosion is a key process that may cause a large number of economic and environmental problems (e.g. land loss, damage to structures and aquatic habitat). Stream bank erosion (toe erosion and mass failure) represents an important form of channel morphology changes and a significant source of sediment. With the advances made in computational techniques, two-dimensional (2-D) numerical models have become valuable tools for investigating flow and sediment transport in open channels at large temporal and spatial scales. However, the implementation of mass failure process in 2D numerical models is still a challenging task. In this paper, a simple, innovative algorithm is implemented in the Telemac-Mascaret modeling platform to handle bank failure: failure occurs whether the actual slope of one given bed element is higher than the internal friction angle. The unstable bed elements are rotated around an appropriate axis, ensuring mass conservation. Mass failure of a bank due to slope instability is applied at the end of each sediment transport evolution iteration, once the bed evolution due to bed load (and/or suspended load) has been computed, but before the global sediment mass balance is verified. This bank failure algorithm is successfully tested using two laboratory experimental cases. Then, bank failure in a 1:40 scale physical model of the Rhine River composed of non-uniform material is simulated. The main features of the bank erosion and failure are correctly reproduced in the numerical simulations, namely the mass wasting at the bank toe, followed by failure at the bank head, and subsequent transport of the mobilised material in an aggradation front. Volumes of eroded material obtained are of the same order of magnitude as the volumes measured during the laboratory tests.

Development of a 3D numerical model to evaluate the Stromboli NW flank instability in relation to magma intrusion

Science.gov (United States)

Apuani, T.; Merri, A.

2009-04-01

A stress-strain analysis of the Stromboli volcano was performed using a three-dimensional explicit finite difference numerical code (FLAC 3D, ITASCA, 2005), to evaluate the effects associated to the presence of magma pressure in magmatic conduit and to foresee the evolution of the magmatic feeding complex. The simulations considered both the ordinary state for the Stromboli, characterized by a partial fill of the active dyke with regular emission of gas and lava fountains and the paroxysmal conditions observed during the March 2007's eruptive crisis, with the magma level in the active dyke reaching the topographic surface along the Sciara del Fuoco depression. The modeling contributes to identify the most probable directions of propagation of new dikes, and the effects of their propagation on the stability of the volcano edifice. The numerical model extends 6 x 6 x 2.6 km3, with a mesh resolution of 100 m, adjusting the grid to fit the shape of the object to be modeled. An elasto-plastic constitutive law was adopted and an homogeneous Mohr-Coulomb strength criterion was chosen for the volcanic cone, assuming one lithotechnical unit (alternation of lava and breccia layers "lava-breccia unit"- Apuani et al 2005). The dykes are represented as discontinuities of the grid, and are modeled by means of interfaces. The magmatic pressure is imposed to the model as normal pressure applied on both sides of the interfaces. The magmastatic pressure was calculated as Pm=d•h, where d is the magma unit weight assumed equal to 25 KN/m3, and h (m) is the height of the magma column. Values of overpressure between 0 and 1 MPa were added to simulate the paroxysmal eruption. The simulation was implemented in successive stages, assuming the results of the previous stages as condition for the next one. A progressive propagation of the dike was simulated, in accordance with the stress conditions identified step by step, and in accordance with the evidences detected by in situ survey, and

Methane Fluxes in West Siberia: 3-D Regional Model Simulation

International Nuclear Information System (INIS)

Jagovkina, S. V.; Karol, I. L.; Zubov, V. A.; Lagun, V. E.; Reshetnikov, A. I.; Rozanov, E. V.

2001-01-01

The West Siberian region is one of the main contributors of the atmospheric greenhouse gas methane due to the large areas of wetlands, rivers, lakes and numerous gas deposits situated there.But there are no reliable estimations of integral methane flux from this area into the atmosphere. For assessment of methane fluxes in West Siberia the specially constructed 3-D regional chemical transport model was applied. The 3-D distribution of methane is calculated on the basis of the current meteorological data fields(wind, temperature, geopotential) updated 4 times a day. The methane concentrations measured near the main gas fields of West Siberia in the summer season of 1999, were used for correction of methane flux intensity estimates obtained previously by comparison of measurements carried out in summer 1993 and 1996 with modelled methane mixing ratio distribution. This set of field and model experiments confirmed the preliminary conclusion about low leakage intensity: anthropogenic methane flux does not exceed 5-15% of total summer methane flux, estimated as 11-12 Mt CH 4 in summer from this region, in spite of the large areas of gas deposits located there

3D Modeling Techniques for Print and Digital Media

Science.gov (United States)

Stephens, Megan Ashley

In developing my thesis, I looked to gain skills using ZBrush to create 3D models, 3D scanning, and 3D printing. The models created compared the hearts of several vertebrates and were intended for students attending Comparative Vertebrate Anatomy. I used several resources to create a model of the human heart and was able to work from life while creating heart models from other vertebrates. I successfully learned ZBrush and 3D scanning, and successfully printed 3D heart models. ZBrush allowed me to create several intricate models for use in both animation and print media. The 3D scanning technique did not fit my needs for the project, but may be of use for later projects. I was able to 3D print using two different techniques as well.

Which Fault Orientations Occur during Oblique Rifting? Combining Analog and Numerical 3d Models with Observations from the Gulf of Aden

Science.gov (United States)

Autin, J.; Brune, S.

2013-12-01

Oblique rift systems like the Gulf of Aden are intrinsically three-dimensional. In order to understand the evolution of these systems, one has to decode the fundamental mechanical similarities of oblique rifts. One way to accomplish this, is to strip away the complexity that is generated by inherited fault structures. In doing so, we assume a laterally homogeneous segment of Earth's lithosphere and ask how many different fault populations are generated during oblique extension inbetween initial deformation and final break-up. We combine results of an analog and a numerical model that feature a 3D segment of a layered lithosphere. In both cases, rift evolution is recorded quantitatively in terms of crustal fault geometries. For the numerical model, we adopt a novel post-processing method that allows to infer small-scale crustal fault orientation from the surface stress tensor. Both models involve an angle of 40 degrees between the rift normal and the extensional direction which allows comparison to the Gulf of Aden rift system. The resulting spatio-temporal fault pattern of our models shows three normal fault orientations: rift-parallel, extension-orthogonal, and intermediate, i.e. with a direction inbetween the two previous orientations. The rift evolution involves three distinct phases: (i) During the initial rift phase, wide-spread faulting with intermediate orientation occurs. (ii) Advanced lithospheric necking enables rift-parallel normal faulting at the rift flanks, while strike-slip faulting in the central part of the rift system indicates strain partitioning. (iii) During continental break-up, displacement-orthogonal as well as intermediate faults occur. We compare our results to the structural evolution of the Eastern Gulf of Aden. External parts of the rift exhibit intermediate and displacement-orthogonal faults while rift-parallel faults are present at the rift borders. The ocean-continent transition mainly features intermediate and displacement

Modeling radiation belt dynamics using a 3-D layer method code

Science.gov (United States)

Wang, C.; Ma, Q.; Tao, X.; Zhang, Y.; Teng, S.; Albert, J. M.; Chan, A. A.; Li, W.; Ni, B.; Lu, Q.; Wang, S.

2017-08-01

A new 3-D diffusion code using a recently published layer method has been developed to analyze radiation belt electron dynamics. The code guarantees the positivity of the solution even when mixed diffusion terms are included. Unlike most of the previous codes, our 3-D code is developed directly in equatorial pitch angle (α0), momentum (p), and L shell coordinates; this eliminates the need to transform back and forth between (α0,p) coordinates and adiabatic invariant coordinates. Using (α0,p,L) is also convenient for direct comparison with satellite data. The new code has been validated by various numerical tests, and we apply the 3-D code to model the rapid electron flux enhancement following the geomagnetic storm on 17 March 2013, which is one of the Geospace Environment Modeling Focus Group challenge events. An event-specific global chorus wave model, an AL-dependent statistical plasmaspheric hiss wave model, and a recently published radial diffusion coefficient formula from Time History of Events and Macroscale Interactions during Substorms (THEMIS) statistics are used. The simulation results show good agreement with satellite observations, in general, supporting the scenario that the rapid enhancement of radiation belt electron flux for this event results from an increased level of the seed population by radial diffusion, with subsequent acceleration by chorus waves. Our results prove that the layer method can be readily used to model global radiation belt dynamics in three dimensions.

AUTOMATIC MESH GENERATION OF 3-D GEOMETRIC MODELS

Institute of Scientific and Technical Information of China (English)

刘剑飞

2003-01-01

In this paper the presentation of the ball-packing method is reviewed,and a scheme to generate mesh for complex 3-D geometric models is given,which consists of 4 steps:(1)create nodes in 3-D models by ball-packing method,(2)connect nodes to generate mesh by 3-D Delaunay triangulation,(3)retrieve the boundary of the model after Delaunay triangulation,(4)improve the mesh.

Comparison of 3D Echocardiogram-Derived 3D Printed Valve Models to Molded Models for Simulated Repair of Pediatric Atrioventricular Valves.

Science.gov (United States)

Scanlan, Adam B; Nguyen, Alex V; Ilina, Anna; Lasso, Andras; Cripe, Linnea; Jegatheeswaran, Anusha; Silvestro, Elizabeth; McGowan, Francis X; Mascio, Christopher E; Fuller, Stephanie; Spray, Thomas L; Cohen, Meryl S; Fichtinger, Gabor; Jolley, Matthew A

2018-03-01

Mastering the technical skills required to perform pediatric cardiac valve surgery is challenging in part due to limited opportunity for practice. Transformation of 3D echocardiographic (echo) images of congenitally abnormal heart valves to realistic physical models could allow patient-specific simulation of surgical valve repair. We compared materials, processes, and costs for 3D printing and molding of patient-specific models for visualization and surgical simulation of congenitally abnormal heart valves. Pediatric atrioventricular valves (mitral, tricuspid, and common atrioventricular valve) were modeled from transthoracic 3D echo images using semi-automated methods implemented as custom modules in 3D Slicer. Valve models were then both 3D printed in soft materials and molded in silicone using 3D printed "negative" molds. Using pre-defined assessment criteria, valve models were evaluated by congenital cardiac surgeons to determine suitability for simulation. Surgeon assessment indicated that the molded valves had superior material properties for the purposes of simulation compared to directly printed valves (p 3D echo-derived molded valves are a step toward realistic simulation of complex valve repairs but require more time and labor to create than directly printed models. Patient-specific simulation of valve repair in children using such models may be useful for surgical training and simulation of complex congenital cases.

Improving 3d Spatial Queries Search: Newfangled Technique of Space Filling Curves in 3d City Modeling

Science.gov (United States)

Uznir, U.; Anton, F.; Suhaibah, A.; Rahman, A. A.; Mioc, D.

2013-09-01

The advantages of three dimensional (3D) city models can be seen in various applications including photogrammetry, urban and regional planning, computer games, etc.. They expand the visualization and analysis capabilities of Geographic Information Systems on cities, and they can be developed using web standards. However, these 3D city models consume much more storage compared to two dimensional (2D) spatial data. They involve extra geometrical and topological information together with semantic data. Without a proper spatial data clustering method and its corresponding spatial data access method, retrieving portions of and especially searching these 3D city models, will not be done optimally. Even though current developments are based on an open data model allotted by the Open Geospatial Consortium (OGC) called CityGML, its XML-based structure makes it challenging to cluster the 3D urban objects. In this research, we propose an opponent data constellation technique of space-filling curves (3D Hilbert curves) for 3D city model data representation. Unlike previous methods, that try to project 3D or n-dimensional data down to 2D or 3D using Principal Component Analysis (PCA) or Hilbert mappings, in this research, we extend the Hilbert space-filling curve to one higher dimension for 3D city model data implementations. The query performance was tested using a CityGML dataset of 1,000 building blocks and the results are presented in this paper. The advantages of implementing space-filling curves in 3D city modeling will improve data retrieval time by means of optimized 3D adjacency, nearest neighbor information and 3D indexing. The Hilbert mapping, which maps a subinterval of the [0, 1] interval to the corresponding portion of the d-dimensional Hilbert's curve, preserves the Lebesgue measure and is Lipschitz continuous. Depending on the applications, several alternatives are possible in order to cluster spatial data together in the third dimension compared to its

Illustrating the disassembly of 3D models

KAUST Repository

Guo, Jianwei; Yan, Dongming; Li, Er; Dong, Weiming; Wonka, Peter; Zhang, Xiaopeng

2013-01-01

We present a framework for the automatic disassembly of 3D man-made models and the illustration of the disassembly process. Given an assembled 3D model, we first analyze the individual parts using sharp edge loops and extract the contact faces

3D physical modeling for patterning process development

Science.gov (United States)

Sarma, Chandra; Abdo, Amr; Bailey, Todd; Conley, Will; Dunn, Derren; Marokkey, Sajan; Talbi, Mohamed

2010-03-01

In this paper we will demonstrate how a 3D physical patterning model can act as a forensic tool for OPC and ground-rule development. We discuss examples where the 2D modeling shows no issues in printing gate lines but 3D modeling shows severe resist loss in the middle. In absence of corrective measure, there is a high likelihood of line discontinuity post etch. Such early insight into process limitations of prospective ground rules can be invaluable for early technology development. We will also demonstrate how the root cause of broken poly-line after etch could be traced to resist necking in the region of STI step with the help of 3D models. We discuss different cases of metal and contact layouts where 3D modeling gives an early insight in to technology limitations. In addition such a 3D physical model could be used for early resist evaluation and selection for required ground-rule challenges, which can substantially reduce the cycle time for process development.

DYNA3D2000*, Explicit 3-D Hydrodynamic FEM Program

International Nuclear Information System (INIS)

Lin, J.

2002-01-01

1 - Description of program or function: DYNA3D2000 is a nonlinear explicit finite element code for analyzing 3-D structures and solid continuum. The code is vectorized and available on several computer platforms. The element library includes continuum, shell, beam, truss and spring/damper elements to allow maximum flexibility in modeling physical problems. Many materials are available to represent a wide range of material behavior, including elasticity, plasticity, composites, thermal effects and rate dependence. In addition, DYNA3D has a sophisticated contact interface capability, including frictional sliding, single surface contact and automatic contact generation. 2 - Method of solution: Discretization of a continuous model transforms partial differential equations into algebraic equations. A numerical solution is then obtained by solving these algebraic equations through a direct time marching scheme. 3 - Restrictions on the complexity of the problem: Recent software improvements have eliminated most of the user identified limitations with dynamic memory allocation and a very large format description that has pushed potential problem sizes beyond the reach of most users. The dominant restrictions remain in code execution speed and robustness, which the developers constantly strive to improve

Towards the graviton from spinfoams: the 3d toy model

International Nuclear Information System (INIS)

Speziale, Simone

2006-01-01

Recently, a proposal has appeared for the extraction of the 2-point function of linearised quantum gravity, within the spinfoam formalism. This relies on the use of a boundary state, which introduces a semi-classical flat geometry on the boundary. In this paper, we investigate this proposal considering a toy model in the (Riemannian) 3d case, where the semi-classical limit is better understood. We show that in this limit the propagation kernel of the model is the one for the harmonic oscillator. This is at the origin of the expected 1/l behaviour of the 2-point function. Furthermore, we numerically study the short scales regime, where deviations from this behaviour occur

Numerical Study of Transonic Axial Flow Rotating Cascade Aerodynamics – Part 1: 2D Case

Directory of Open Access Journals (Sweden)

Irina Carmen ANDREI

2014-06-01

Full Text Available The purpose of this paper is to present a 2D study regarding the numerical simulation of flow within a transonic highly-loaded rotating cascade from an axial compressor. In order to describe an intricate flow pattern of a complex geometry and given specific conditions of cascade’s loading and operation, an appropriate accurate flow model is a must. For such purpose, the Navier-Stokes equations system was used as flow model; from the computational point of view, the mathematical support is completed by a turbulence model. A numerical comparison has been performed for different turbulence models (e.g. KE, KO, Reynolds Stress and Spallart-Allmaras models. The convergence history was monitored in order to focus on the numerical accuracy. The force vector has been reported in order to express the aerodynamics of flow within the rotating cascade at the running regime, in terms of Lift and Drag. The numerical results, expressed by plots of the most relevant flow parameters, have been compared. It comes out that the selecting of complex flow models and appropriate turbulence models, in conjunction with CFD techniques, allows to obtain the best computational accuracy of the numerical results. This paper aims to carry on a 2D study and a prospective 3D will be intended for the same architecture.

3D visualization of numeric planetary data using JMARS

Science.gov (United States)

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.

3D Face modeling using the multi-deformable method.

Science.gov (United States)

Hwang, Jinkyu; Yu, Sunjin; Kim, Joongrock; Lee, Sangyoun

2012-09-25

In this paper, we focus on the problem of the accuracy performance of 3D face modeling techniques using corresponding features in multiple views, which is quite sensitive to feature extraction errors. To solve the problem, we adopt a statistical model-based 3D face modeling approach in a mirror system consisting of two mirrors and a camera. The overall procedure of our 3D facial modeling method has two primary steps: 3D facial shape estimation using a multiple 3D face deformable model and texture mapping using seamless cloning that is a type of gradient-domain blending. To evaluate our method's performance, we generate 3D faces of 30 individuals and then carry out two tests: accuracy test and robustness test. Our method shows not only highly accurate 3D face shape results when compared with the ground truth, but also robustness to feature extraction errors. Moreover, 3D face rendering results intuitively show that our method is more robust to feature extraction errors than other 3D face modeling methods. An additional contribution of our method is that a wide range of face textures can be acquired by the mirror system. By using this texture map, we generate realistic 3D face for individuals at the end of the paper.

3D nonlinear modeling of the coupling and phase locking of magnetic Islands in tokamaks

Science.gov (United States)

Jardin, Stephen; Ferraro, Nathaniel; Chen, Jin; Pfefferle, David

2017-10-01

Many tokamak discharges develop multiple tearing modes possessing different mode numbers. These modes are observed to phase lock to one another, resulting in a flattening of the core toroidal plasma rotation profile, which can have deleterious effects on transport and MHD stability. In order to study these phenomena with minimum assumptions, we use the M3D-C1 3D nonlinear MHD code to perform initial value simulations of the evolution of equilibria unstable to both the 2/1 and 3/2 modes, but having sheared toroidal rotation. Initial attempts to perform these simulations led to numerical instabilities developing once the islands got to a certain size. In order to study the cause of this instability, we developed a small model code that solves a pure convection equation in 1D. We find that an implicit Crank-Nicholson method in time and Hermite Cubic finite elements (as are used in the toroidal direction in the M3D-C1 code) is not a convergent algorithm. Adding a small second order diffusion term, proportional to the velocity, improves the numerical stability properties but is not convergent in the first-derivative of the solution. Instead, adding a much smaller forth-order spatial derivative term proportional to the velocity leads to an algorithm in which both the solution and the first derivative converge as 1/N2,. Adding similar toroidal forth derivative terms to the M3D-C1 code eliminated the numerical instability. This work was supported by US DOE Contract DE-AC02-09-CH11466.

3-D time-dependent numerical model of flow patterns within a large-scale Czochralski system

Science.gov (United States)

Nam, Phil-Ouk; O, Sang-Kun; Yi, Kyung-Woo

2008-04-01

Silicon single crystals grown through the Czochralski (Cz) method have increased in size to 300 mm, resulting in the use of larger crucibles. The objective of this study is to investigate the continuous Cz method in a large crucible (800 mm), which is performed by inserting a polycrystalline silicon rod into the melt. The numerical model is based on a time-dependent and three-dimensional standard k- ɛ turbulent model using the analytical software package CFD-ACE+, version 2007. Wood's metal melt, which has a low melting point ( Tm=70 °C), was used as the modeling fluid. Crystal rotation given in the clockwise direction with rotation rates varying from 0 to 15 rpm, while the crucible was rotated counter-clockwise, with rotation rates between 0 and 3 rpm. The results show that asymmetrical phenomena of fluid flow arise as results of crystal and crucible rotation, and that these phenomena move with the passage of time. Near the crystal, the flow moves towards the crucible at the pole of the asymmetrical phenomena. Away from the poles, a vortex begins to form, which is strongly pronounced in the region between the poles.

Application of the perfectly matched layer in 3-D marine controlled-source electromagnetic modelling

Science.gov (United States)

Li, Gang; Li, Yuguo; Han, Bo; Liu, Zhan

2018-01-01

In this study, the complex frequency-shifted perfectly matched layer (CFS-PML) in stretching Cartesian coordinates is successfully applied to 3-D frequency-domain marine controlled-source electromagnetic (CSEM) field modelling. The Dirichlet boundary, which is usually used within the traditional framework of EM modelling algorithms, assumes that the electric or magnetic field values are zero at the boundaries. This requires the boundaries to be sufficiently far away from the area of interest. To mitigate the boundary artefacts, a large modelling area may be necessary even though cell sizes are allowed to grow toward the boundaries due to the diffusion of the electromagnetic wave propagation. Compared with the conventional Dirichlet boundary, the PML boundary is preferred as the modelling area of interest could be restricted to the target region and only a few absorbing layers surrounding can effectively depress the artificial boundary effect without losing the numerical accuracy. Furthermore, for joint inversion of seismic and marine CSEM data, if we use the PML for CSEM field simulation instead of the conventional Dirichlet, the modelling area for these two different geophysical data collected from the same survey area could be the same, which is convenient for joint inversion grid matching. We apply the CFS-PML boundary to 3-D marine CSEM modelling by using the staggered finite-difference discretization. Numerical test indicates that the modelling algorithm using the CFS-PML also shows good accuracy compared to the Dirichlet. Furthermore, the modelling algorithm using the CFS-PML shows advantages in computational time and memory saving than that using the Dirichlet boundary. For the 3-D example in this study, the memory saving using the PML is nearly 42 per cent and the time saving is around 48 per cent compared to using the Dirichlet.

Development of a 3D electromagnetic model for eddy current tubing inspection application to steam generator tubing

Energy Technology Data Exchange (ETDEWEB)

Maillot, V. [Institut de Radioprotection et de Surete Nucleaire, IRSN, 92 - Fontenay aux Roses (France); Pichenot, G.; Premel, D.; Sollier, T. [CEA Saclay, DRT/DECS, 91 - Gif-sur-Yvette (France)

2003-10-01

In nuclear plants, the inspection of heat exchanger tubes is usually carried out by using eddy current nondestructive testing. A numerical model, based on a volume integral approach using the Green's dyadic formalism, has been developed, with support from the French Institute for Radiological Protection and Nuclear Safety, to predict the response of an eddy current bobbin coil to 3D flaws located in the tube's wall. With an aim of integrating this model into the NDE multi techniques platform CIVA, it has been validated with experimental data for 2D and 3D flaws. (authors)

Stress field sensitivity analysis within Mesozoic successions in the Swiss Alpine foreland using 3-D-geomechanical-numerical models

Science.gov (United States)

Reiter, Karsten; Hergert, Tobias; Heidbach, Oliver

2016-04-01

The in situ stress conditions are of key importance for the evaluation of radioactive waste repositories. In stage two of the Swiss site selection program, the three siting areas of high-level radioactive waste are located in the Alpine foreland in northern Switzerland. The sedimentary succession overlays the basement, consisting of variscan crystalline rocks as well as partly preserved Permo-Carboniferous deposits in graben structures. The Mesozoic sequence represents nearly the complete era and is covered by Cenozoic Molasse deposits as well as Quaternary sediments, mainly in the valleys. The target horizon (designated host rock) is an >100 m thick argillaceous Jurassic deposit (Opalinus Clay). To enlighten the impact of site-specific features on the state of stress within the sedimentary succession, 3-D-geomechanical-numerical models with elasto-plastic rock properties are set up for three potential siting areas. The lateral extent of the models ranges between 12 and 20 km, the vertical extent is up to a depth of 2.5 or 5 km below sea level. The sedimentary sequence plus the basement are separated into 10 to 14 rock mechanical units. The Mesozoic succession is intersected by regional fault zones; two or three of them are present in each model. The numerical problem is solved with the finite element method with a resolution of 100-150 m laterally and 10-30 m vertically. An initial stress state is established for all models taking into account the depth-dependent overconsolidation ratio in Opalinus Clay in northern Switzerland. The influence of topography, rock properties, friction on the faults as well as the impact of tectonic shortening on the state of stress is investigated. The tectonic stress is implemented with lateral displacement boundary conditions, calibrated on stress data that are compiled in Northern Switzerland. The model results indicate that the stress perturbation by the topography is significant to depths greater than the relief contrast. The

Automatic balancing of 3D models

DEFF Research Database (Denmark)

Christiansen, Asger Nyman; Schmidt, Ryan; Bærentzen, Jakob Andreas

2014-01-01

3D printing technologies allow for more diverse shapes than are possible with molds and the cost of making just one single object is negligible compared to traditional production methods. However, not all shapes are suitable for 3D print. One of the remaining costs is therefore human time spent......, in these cases, we will apply a rotation of the object which only deforms the shape a little near the base. No user input is required but it is possible to specify manufacturing constraints related to specific 3D print technologies. Several models have successfully been balanced and printed using both polyjet...... is solved by creating cavities of air and distributing dense materials inside the model. Consequently, the surface is not deformed. However, printing materials with significantly different densities is often not possible and adding cavities of air is often not enough to make the model balance. Consequently...

A spherical harmonics intensity model for 3D segmentation and 3D shape analysis of heterochromatin foci.

Science.gov (United States)

Eck, Simon; Wörz, Stefan; Müller-Ott, Katharina; Hahn, Matthias; Biesdorf, Andreas; Schotta, Gunnar; Rippe, Karsten; Rohr, Karl

2016-08-01

The genome is partitioned into regions of euchromatin and heterochromatin. The organization of heterochromatin is important for the regulation of cellular processes such as chromosome segregation and gene silencing, and their misregulation is linked to cancer and other diseases. We present a model-based approach for automatic 3D segmentation and 3D shape analysis of heterochromatin foci from 3D confocal light microscopy images. Our approach employs a novel 3D intensity model based on spherical harmonics, which analytically describes the shape and intensities of the foci. The model parameters are determined by fitting the model to the image intensities using least-squares minimization. To characterize the 3D shape of the foci, we exploit the computed spherical harmonics coefficients and determine a shape descriptor. We applied our approach to 3D synthetic image data as well as real 3D static and real 3D time-lapse microscopy images, and compared the performance with that of previous approaches. It turned out that our approach yields accurate 3D segmentation results and performs better than previous approaches. We also show that our approach can be used for quantifying 3D shape differences of heterochromatin foci. Copyright © 2016 Elsevier B.V. All rights reserved.

Development of 3D and 4D Bridge Models and Plans

Science.gov (United States)

2018-05-28

Since 2012, MDOT has been leading national efforts to modernize design development with 3D modeling. Early focus on roadway projects yielded streamlined plan production and digital data for construction. As MDOT pivots to 3D model-centric design, nat...

Numerical simulation of 3-D turbulent flow through entire stage in a multistage centrifugal pump

International Nuclear Information System (INIS)

Huang, S.; Islam, M.F.; Liu, P.

2005-01-01

A three-dimensional turbulent flow through a multistage centrifugal pump is numerically simulated using a commercial CFD software package. The simulation and analysis include flow fields in rotating impeller and stationary diffuser and is completed in a multiple reference frame. The standard k-ε turbulence model is applied. The analysis of the simulation reveals that the reverse flows exist in the zone near the impeller exit and diffuser entrance, resulting in flow field asymmetric and unsteady. There is a considerable interference on velocity field at impeller exit due to the interaction between impeller blades and diffuser vanes. The hydraulic performance is connected and evaluated with the 3-D computational flow field. The current computation is verified by comparing predicted and measured head. (author)

3-D model-based vehicle tracking.

Science.gov (United States)

Lou, Jianguang; Tan, Tieniu; Hu, Weiming; Yang, Hao; Maybank, Steven J

2005-10-01

This paper aims at tracking vehicles from monocular intensity image sequences and presents an efficient and robust approach to three-dimensional (3-D) model-based vehicle tracking. Under the weak perspective assumption and the ground-plane constraint, the movements of model projection in the two-dimensional image plane can be decomposed into two motions: translation and rotation. They are the results of the corresponding movements of 3-D translation on the ground plane (GP) and rotation around the normal of the GP, which can be determined separately. A new metric based on point-to-line segment distance is proposed to evaluate the similarity between an image region and an instantiation of a 3-D vehicle model under a given pose. Based on this, we provide an efficient pose refinement method to refine the vehicle's pose parameters. An improved EKF is also proposed to track and to predict vehicle motion with a precise kinematics model. Experimental results with both indoor and outdoor data show that the algorithm obtains desirable performance even under severe occlusion and clutter.

Numerical simulation of 3D boundary-driven acoustic streaming in microfluidic devices.

Science.gov (United States)

Lei, Junjun; Hill, Martyn; Glynne-Jones, Peter

2014-02-07

This article discusses three-dimensional (3D) boundary-driven streaming in acoustofluidic devices. Firstly, the 3D Rayleigh streaming pattern in a microchannel is simulated and its effect on the movement of microparticles of various sizes is demonstrated. The results obtained from this model show good comparisons with 3D experimental visualisations and demonstrate the fully 3D nature of the acoustic streaming field and the associated acoustophoretic motion of microparticles in acoustofluidic devices. This method is then applied to another acoustofluidic device in order to gain insights into an unusual in-plane streaming pattern. The origin of this streaming has not been fully described and its characteristics cannot be explained from the classical theory of Rayleigh streaming. The simulated in-plane streaming pattern was in good agreement with the experimental visualisation. The mechanism behind it is shown to be related to the active sound intensity field, which supports our previous findings on the mechanism of the in-plane acoustic streaming pattern visualised and modelled in a thin-layered capillary device.

3D Bioprinting of Tissue/Organ Models.

Science.gov (United States)

Pati, Falguni; Gantelius, Jesper; Svahn, Helene Andersson

2016-04-04

In vitro tissue/organ models are useful platforms that can facilitate systematic, repetitive, and quantitative investigations of drugs/chemicals. The primary objective when developing tissue/organ models is to reproduce physiologically relevant functions that typically require complex culture systems. Bioprinting offers exciting prospects for constructing 3D tissue/organ models, as it enables the reproducible, automated production of complex living tissues. Bioprinted tissues/organs may prove useful for screening novel compounds or predicting toxicity, as the spatial and chemical complexity inherent to native tissues/organs can be recreated. In this Review, we highlight the importance of developing 3D in vitro tissue/organ models by 3D bioprinting techniques, characterization of these models for evaluating their resemblance to native tissue, and their application in the prioritization of lead candidates, toxicity testing, and as disease/tumor models. © 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Transition of R&D into Operations at Fleet Numerical Meteorology and Oceanography Center

Science.gov (United States)

Clancy, R. M.

2006-12-01

The U.S. Navy's Fleet Numerical Meteorology and Oceanography Center (FNMOC) plays a significant role in the National capability for operational weather and ocean prediction through its operation of sophisticated global and regional meteorological and oceanographic models, extending from the top of the atmosphere to the bottom of the ocean. FNMOC uniquely satisfies the military's requirement for a global operational weather prediction capability based on software certified to DoD Information Assurance standards and operated in a secure classified computer environment protected from outside intrusion by DoD certified firewalls. FNMOC operates around-the-clock, 365 days per year and distributes products to military and civilian users around the world, both ashore and afloat, through a variety of means. FNMOC's customers include all branches of the Department of Defense, other government organizations such as the National Weather Service, private companies, a number of colleges and universities, and the general public. FNMOC employs three primary models, the Navy Operational Global Atmospheric Prediction System (NOGAPS), the Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS), and the WaveWatch III model (WW3), along with a number of specialized models and related applications. NOGAPS is a global weather model, driving nearly all other FNMOC models and applications in some fashion. COAMPS is a high- resolution regional model that has proved to be particularly valuable for forecasting weather and ocean conditions in highly complex coastal areas. WW3 is a state-of-the-art ocean wave model that is employed both globally and regionally in support of a wide variety of naval operations. Other models support and supplement the main models with predictions of ocean thermal structure, ocean currents, sea-ice characteristics, and other data. Fleet Numerical operates at the leading edge of science and technology, and benefits greatly from collocation with its supporting

Parallel 3-D numerical simulation of dielectric barrier discharge plasma actuators

Science.gov (United States)

Houba, Tomas

Dielectric barrier discharge plasma actuators have shown promise in a range of applications including flow control, sterilization and ozone generation. Developing numerical models of plasma actuators is of great importance, because a high-fidelity parallel numerical model allows new design configurations to be tested rapidly. Additionally, it provides a better understanding of the plasma actuator physics which is useful for further innovation. The physics of plasma actuators is studied numerically. A loosely coupled approach is utilized for the coupling of the plasma to the neutral fluid. The state of the art in numerical plasma modeling is advanced by the development of a parallel, three-dimensional, first-principles model with detailed air chemistry. The model incorporates 7 charged species and 18 reactions, along with a solution of the electron energy equation. To the author's knowledge, a parallel three-dimensional model of a gas discharge with a detailed air chemistry model and the solution of electron energy is unique. Three representative geometries are studied using the gas discharge model. The discharge of gas between two parallel electrodes is used to validate the air chemistry model developed for the gas discharge code. The gas discharge model is then applied to the discharge produced by placing a dc powered wire and grounded plate electrodes in a channel. Finally, a three-dimensional simulation of gas discharge produced by electrodes placed inside a riblet is carried out. The body force calculated with the gas discharge model is loosely coupled with a fluid model to predict the induced flow inside the riblet.

A Sketching Interface for Freeform 3D Modeling

Science.gov (United States)

Igarashi, Takeo

This chapter introduces Teddy, a sketch-based modeling system to quickly and easily design freeform models such as stuffed animals and other rotund objects. The user draws several 2D freeform strokes interactively on the screen and the system automatically constructs plausible 3D polygonal surfaces. Our system supports several modeling operations, including the operation to construct a 3D polygonal surface from a 2D silhouette drawn by the user: it inflates the region surrounded by the silhouette making a wide area fat, and a narrow area thin. Teddy, our prototype system, is implemented as a Java program, and the mesh construction is done in real-time on a standard PC. Our informal user study showed that a first-time user masters the operations within 10 minutes, and can construct interesting 3D models within minutes. We also report the result of a case study where a high school teacher taught various 3D concepts in geography using the system.

Using 3D Morphable Models for face recognition in video

NARCIS (Netherlands)

van Rootseler, R.T.A.; Spreeuwers, Lieuwe Jan; Veldhuis, Raymond N.J.

The 3D Morphable Face Model (3DMM) has been used for over a decade for creating 3D models from single images of faces. This model is based on a PCA model of the 3D shape and texture generated from a limited number of 3D scans. The goal of fitting a 3DMM to an image is to find the model coefficients,

DynEarthSol3D: numerical studies of basal crevasses and calving blocks

Science.gov (United States)

Logan, E.; Lavier, L. L.; Choi, E.; Tan, E.; Catania, G. A.

2014-12-01

DynEarthSol3D (DES) is a thermomechanical model for the simulation of dynamic ice flow. We present the application of DES toward two case studies - basal crevasses and calving blocks - to illustrate the potential of the model to aid in understanding calving processes. Among the advantages of using DES are: its unstructured meshes which adaptively resolve zones of high interest; its use of multiple rheologies to simulate different types of dynamic behavior; and its explicit and parallel numerical core which both make the implementation of different boundary conditions easy and the model highly scalable. We examine the initiation and development of both basal crevasses and calving blocks through time using visco-elasto-plastic rheology. Employing a brittle-to-ductile transition zone (BDTZ) based on local strain rate shows that the style and development of brittle features like crevasses differs markedly on the rheological parameters. Brittle and ductile behavior are captured by Mohr-Coulomb elastoplasticity and Maxwell viscoelasticity, respectively. We explore the parameter spaces which define these rheologies (including temperature) as well as the BDTZ threshold (shown in the literature as 10-7 Pa s), using time-to-failure as a metric for accuracy within the model. As the time it takes for a block of ice to fail can determine an iceberg's size, this work has implications for calving laws.

Continuous modelling study of numerical volumes - Applications to the visualization of anatomical structures

International Nuclear Information System (INIS)

Goret, C.

1990-12-01

Several technics of imaging (IRM, image scanners, tomoscintigraphy, echography) give numerical informations presented by means of a stack of parallel cross-sectional images. Since many years, 3-D mathematical tools have been developed and allow the 3 D images synthesis of surfaces. In first part, we give the technics of numerical volume exploitation and their medical applications to diagnosis and therapy. The second part is about a continuous modelling of the volume with a tensor product of cubic splines. We study the characteristics of this representation and its clinical validation. Finally, we treat of the problem of surface visualization of objects contained in the volume. The results show the interest of this model and allow to propose specifications for 3-D workstation realization [fr

3D Models of Immunotherapy

Science.gov (United States)

This collaborative grant is developing 3D models of both mouse and human biology to investigate aspects of therapeutic vaccination in order to answer key questions relevant to human cancer immunotherapy.

Improving 3D spatial queries search: newfangled technique of space filling curves in 3D city modeling

DEFF Research Database (Denmark)

Uznir, U.; Anton, François; Suhaibah, A.

2013-01-01

, retrieving portions of and especially searching these 3D city models, will not be done optimally. Even though current developments are based on an open data model allotted by the Open Geospatial Consortium (OGC) called CityGML, its XML-based structure makes it challenging to cluster the 3D urban objects......The advantages of three dimensional (3D) city models can be seen in various applications including photogrammetry, urban and regional planning, computer games, etc.. They expand the visualization and analysis capabilities of Geographic Information Systems on cities, and they can be developed using...... modeling will improve data retrieval time by means of optimized 3D adjacency, nearest neighbor information and 3D indexing. The Hilbert mapping, which maps a subinterval of the [0, 1] interval to the corresponding portion of the d-dimensional Hilbert’s curve, preserves the Lebesgue measure and is Lipschitz...

Numerical modeling and experimental testing of a wave energy converter: deliverable D4.2

Energy Technology Data Exchange (ETDEWEB)

Zurkinden, A.S.; Kramer, M.; Ferri, F.; Kofoed, J.P.

2013-05-15

The objective of this document is to summarize the outcome of the research which has been carried out during the period May 2011 until June 2012 i.e. during the first year of the PhD study. The work has been done in collaboration with the co-authors. The aim of the project was primarily to provide numerical values for comparison with the experimental test results which were carried out in the same time. It is for this reason why Chapter 4 does consist exclusively of numerical values. Experimental values and measured time series of wave elevations have been used throughout the report in order to a) validate the numerical model and b) preform stochastic analysis. The latter technique is introduced in order to optimize the control parameters of the power take off system. (Author)

NoSQL Based 3D City Model Management System

Science.gov (United States)

Mao, B.; Harrie, L.; Cao, J.; Wu, Z.; Shen, J.

2014-04-01

To manage increasingly complicated 3D city models, a framework based on NoSQL database is proposed in this paper. The framework supports import and export of 3D city model according to international standards such as CityGML, KML/COLLADA and X3D. We also suggest and implement 3D model analysis and visualization in the framework. For city model analysis, 3D geometry data and semantic information (such as name, height, area, price and so on) are stored and processed separately. We use a Map-Reduce method to deal with the 3D geometry data since it is more complex, while the semantic analysis is mainly based on database query operation. For visualization, a multiple 3D city representation structure CityTree is implemented within the framework to support dynamic LODs based on user viewpoint. Also, the proposed framework is easily extensible and supports geoindexes to speed up the querying. Our experimental results show that the proposed 3D city management system can efficiently fulfil the analysis and visualization requirements.

Integration of aerial oblique imagery and terrestrial imagery for optimized 3D modeling in urban areas

Science.gov (United States)

Wu, Bo; Xie, Linfu; Hu, Han; Zhu, Qing; Yau, Eric

2018-05-01

Photorealistic three-dimensional (3D) models are fundamental to the spatial data infrastructure of a digital city, and have numerous potential applications in areas such as urban planning, urban management, urban monitoring, and urban environmental studies. Recent developments in aerial oblique photogrammetry based on aircraft or unmanned aerial vehicles (UAVs) offer promising techniques for 3D modeling. However, 3D models generated from aerial oblique imagery in urban areas with densely distributed high-rise buildings may show geometric defects and blurred textures, especially on building façades, due to problems such as occlusion and large camera tilt angles. Meanwhile, mobile mapping systems (MMSs) can capture terrestrial images of close-range objects from a complementary view on the ground at a high level of detail, but do not offer full coverage. The integration of aerial oblique imagery with terrestrial imagery offers promising opportunities to optimize 3D modeling in urban areas. This paper presents a novel method of integrating these two image types through automatic feature matching and combined bundle adjustment between them, and based on the integrated results to optimize the geometry and texture of the 3D models generated from aerial oblique imagery. Experimental analyses were conducted on two datasets of aerial and terrestrial images collected in Dortmund, Germany and in Hong Kong. The results indicate that the proposed approach effectively integrates images from the two platforms and thereby improves 3D modeling in urban areas.

geomIO: A tool for geodynamicists to turn 2D cross-sections into 3D geometries

Science.gov (United States)

Baumann, Tobias; Bauville, Arthur

2016-04-01

In numerical deformation models, material properties are usually defined on elements (e.g., in body-fitted finite elements), or on a set of Lagrangian markers (Eulerian, ALE or mesh-free methods). In any case, geometrical constraints are needed to assign different material properties to the model domain. Whereas simple geometries such as spheres, layers or cuboids can easily be programmed, it quickly gets complex and time-consuming to create more complicated geometries for numerical model setups, especially in three dimensions. geomIO (geometry I/O, http://geomio.bitbucket.org/) is a MATLAB-based library that has two main functionalities. First, it can be used to create 3D volumes based on series of 2D vector drawings similar to a CAD program; and second, it uses these 3D volumes to assign material properties to the numerical model domain. The drawings can conveniently be created using the open-source vector graphics software Inkscape. Adobe Illustrator is also partially supported. The drawings represent a series of cross-sections in the 3D model domain, for example, cross-sectional interpretations of seismic tomography. geomIO is then used to read the drawings and to create 3D volumes by interpolating between the cross-sections. In the second part, the volumes are used to assign material phases to markers inside the volumes. Multiple volumes can be created at the same time and, depending on the order of assignment, unions or intersections can be built to assign additional material phases. geomIO also offers the possibility to create 3D temperature structures for geodynamic models based on depth dependent parameterisations, for example the half space cooling model. In particular, this can be applied to geometries of subducting slabs of arbitrary shape. Yet, geomIO is held very general, and can be used for a variety of applications. We present examples of setup generation from pictures of micro-scale tectonics and lithospheric scale setups of 3D present-day model

2D Numerical Modelling of the Resin Injection Pultrusion Process Including Experimental Resin Kinetics and Temperature Validation

DEFF Research Database (Denmark)

Rasmussen, Filip Salling; Sonne, Mads Rostgaard; Larsen, Martin

In the present study, a two-dimensional (2D) transient Eulerian thermo-chemical analysis of a carbon fibre epoxy thermosetting Resin Injection Pultrusion (RIP) process is carried out. The numerical model is implemented using the well known unconditionally stable Alternating Direction Implicit (ADI......) scheme. The total heat of reaction and the cure kinetics of the epoxy thermosetting are determined using Differential Scanning Calorimetry (DSC). A very good agreement is observed between the fitted cure kinetic model and the experimental measurements. The numerical steady state temperature predictions...

Prediction of hydroxyl concentrations in cement pore water using a numerical cement hydration model

NARCIS (Netherlands)

Eijk, van R.J.; Brouwers, H.J.H.

2000-01-01

In this paper, a 3D numerical cement hydration model is used for predicting alkali and hydroxyl concentrations in cement pore water. First, this numerical model is calibrated for Dutch cement employing both chemical shrinkage and calorimetric experiments. Secondly, the strength development of some

3D unified CFD to modeling of bubbles phenomena

International Nuclear Information System (INIS)

Vladimir V Chudanov; Anna E Aksenova; Valerii A Pervichko

2005-01-01

Full text of publication follows: During of the last ten years the developed numerical methods and algorithms for solving of heat and mass transfer problems in compressible/incompressible fluids were successfully tested at simulation of interaction of two immiscible liquids. Now these computing tools are extended on a case of two-phase flows, such as a liquids-gas system as follows: outside of bubbles the non-stationary incompressible Navier-Stokes equations in the primitive variables coupled with the heat transfer equation are used; inside of bubble a compressible medium model with low Mach limit is applied. To observe of an interface of liquid-gas system we use the modified level set method and three-dimensional advective schemes of TVD-type with small scheme diffusion with use of sub-grid simulation. These schemes with small diffusion were already applied by us under using of sub-grid simulation for interface transfer in case of two non-mixing liquids. For bubble phenomena a numerical technique based on the developed algorithms with a small scheme diffusion, for which the discrete approximations are constructed using the finite-volume methods and fully staggered grids is adapted. Testing of the developed approach is carried out on the set of test problems and a good agreement is obtained between numerical predictions and experimental data. The numerical technique was successfully utilized for numerical support of 3D experiment financed by Nuclear Energy Agency at the Organization economic cooperation and development within the framework of MASKA Project, where computational fluid dynamics of two non-mixing fluids such as corium and steel was investigated. In this paper there is application of developed approach for simulation of bubble flows, in particular, for study of coalescence of two drops. The developed technique has a high degree of efficiency and allows on a personal computer (3 GHz and 2 Gbytes RAM) to carry out CFD calculations on a grid with 10 7

Measuring Visual Closeness of 3-D Models

KAUST Repository

Gollaz Morales, Jose Alejandro

2012-09-01

Measuring visual closeness of 3-D models is an important issue for different problems and there is still no standardized metric or algorithm to do it. The normal of a surface plays a vital role in the shading of a 3-D object. Motivated by this, we developed two applications to measure visualcloseness, introducing normal difference as a parameter in a weighted metric in Metro’s sampling approach to obtain the maximum and mean distance between 3-D models using 3-D and 6-D correspondence search structures. A visual closeness metric should provide accurate information on what the human observers would perceive as visually close objects. We performed a validation study with a group of people to evaluate the correlation of our metrics with subjective perception. The results were positive since the metrics predicted the subjective rankings more accurately than the Hausdorff distance.

Modelling 3D spatial objects in a geo-DBMS using a 3D primitive

Science.gov (United States)

Arens, Călin; Stoter, Jantien; van Oosterom, Peter

2005-03-01

There is a growing interest in modelling the world in three dimensions, both in applications and in science. At the same time, geographical information systems are changing into integrated architecture in which administrative and spatial data are maintained in one environment. It is for this reason that mainstream Data Base Management Systems (DBMSs) have implemented spatial data types according to the 'Simple Feature Specifications for SQL', described by the OpenGeospatial Consortium. However, these specifications are 2D, as indeed are the implementations in DBMSs. At the Section GIS Technology of TU Delft, research has been carried out in which a 3D primitive was implemented in a DBMS (Oracle Spatial). To explore the possibilities and complications, a fairly simple 3D primitive was chosen to start with: a polyhedron. In the future the study will be extended with more complex primitives, the ultimate aim being to build 3D models with features closer to the real world. Besides the data structure, a validation function was developed to check the geometric accuracy of the data. Rules for validation were established and translated into prototype implementations with the aid of literature. In order to manipulate the data, a list of useful 3D functions was specified. Most of these were translated into algorithms, which were implemented in the DBMS. The algorithms for these functions were obtained from the relevant literature. The research also comprised a comparative performance test on spatial indexing in 2D and 3D, using an R-tree. Finally, existing software was used to visualize 3D objects structured with the implemented 3D primitive. This research is a first attempt to implement a true 3D primitive in a DBMS. Future research will focus on extending and improving the implementations and on optimizing maintenance and query of 3D objects in DBMSs.

Nonlinear plasma wave models in 3D fluid simulations of laser-plasma interaction

Science.gov (United States)

Chapman, Thomas; Berger, Richard; Arrighi, Bill; Langer, Steve; Banks, Jeffrey; Brunner, Stephan

2017-10-01

Simulations of laser-plasma interaction (LPI) in inertial confinement fusion (ICF) conditions require multi-mm spatial scales due to the typical laser beam size and durations of order 100 ps in order for numerical laser reflectivities to converge. To be computationally achievable, these scales necessitate a fluid-like treatment of light and plasma waves with a spatial grid size on the order of the light wave length. Plasma waves experience many nonlinear phenomena not naturally described by a fluid treatment, such as frequency shifts induced by trapping, a nonlinear (typically suppressed) Landau damping, and mode couplings leading to instabilities that can cause the plasma wave to decay rapidly. These processes affect the onset and saturation of stimulated Raman and Brillouin scattering, and are of direct interest to the modeling and prediction of deleterious LPI in ICF. It is not currently computationally feasible to simulate these Debye length-scale phenomena in 3D across experimental scales. Analytically-derived and/or numerically benchmarked models of processes occurring at scales finer than the fluid simulation grid offer a path forward. We demonstrate the impact of a range of kinetic processes on plasma reflectivity via models included in the LPI simulation code pF3D. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.

AUTOMATIC MESH GENERATION OF 3—D GEOMETRIC MODELS

Institute of Scientific and Technical Information of China (English)

刘剑飞

2003-01-01

In this paper the presentation of the ball-packing method is reviewed, and a schemeto generate mesh for complex 3-D geometric models is given, which consists of 4 steps: (1) createnodes in 3-D models by ball-packing method, (2) connect nodes to generate mesh by 3-D Delaunaytriangulation, (3) retrieve the boundary of the model after Delaunay triangulation, (4) improve themesh.

A novel alternative method for 3D visualisation in Parasitology: the construction of a 3D model of a parasite from 2D illustrations.

Science.gov (United States)

Teo, B G; Sarinder, K K S; Lim, L H S

2010-08-01

Three-dimensional (3D) models of the marginal hooks, dorsal and ventral anchors, bars and haptoral reservoirs of a parasite, Sundatrema langkawiense Lim & Gibson, 2009 (Monogenea) were developed using the polygonal modelling method in Autodesk 3ds Max (Version 9) based on two-dimensional (2D) illustrations. Maxscripts were written to rotate the modelled 3D structures. Appropriately orientated 3D haptoral hard-parts were then selected and positioned within the transparent 3D outline of the haptor and grouped together to form a complete 3D haptoral entity. This technique is an inexpensive tool for constructing 3D models from 2D illustrations for 3D visualisation of the spatial relationships between the different structural parts within organisms.

Proposal of custom made wrist orthoses based on 3D modelling and 3D printing.

Science.gov (United States)

Abreu de Souza, Mauren; Schmitz, Cristiane; Marega Pinhel, Marcelo; Palma Setti, Joao A; Nohama, Percy

2017-07-01

Accessibility to three-dimensional (3D) technologies, such as 3D scanning systems and additive manufacturing (like 3D printers), allows a variety of 3D applications. For medical applications in particular, these modalities are gaining a lot of attention enabling several opportunities for healthcare applications. The literature brings several cases applying both technologies, but none of them focus on the spreading of how this technology could benefit the health segment. This paper proposes a new methodology, which employs both 3D modelling and 3D printing for building orthoses, which could better fit the demands of different patients. Additionally, there is an opportunity for sharing expertise, as it represents a trendy in terms of the maker-movement. Therefore, as a result of the proposed approach, we present a case study based on a volunteer who needs an immobilization orthosis, which was built for exemplification of the whole process. This proposal also employs freely available 3D models and software, having a strong social impact. As a result, it enables the implementation and effective usability for a variety of built to fit solutions, hitching useful and smarter technologies for the healthcare sector.

Quantitative model for the generic 3D shape of ICMEs at 1 AU

Science.gov (United States)

Démoulin, P.; Janvier, M.; Masías-Meza, J. J.; Dasso, S.

2016-10-01

Context. Interplanetary imagers provide 2D projected views of the densest plasma parts of interplanetary coronal mass ejections (ICMEs), while in situ measurements provide magnetic field and plasma parameter measurements along the spacecraft trajectory, that is, along a 1D cut. The data therefore only give a partial view of the 3D structures of ICMEs. Aims: By studying a large number of ICMEs, crossed at different distances from their apex, we develop statistical methods to obtain a quantitative generic 3D shape of ICMEs. Methods: In a first approach we theoretically obtained the expected statistical distribution of the shock-normal orientation from assuming simple models of 3D shock shapes, including distorted profiles, and compared their compatibility with observed distributions. In a second approach we used the shock normal and the flux rope axis orientations together with the impact parameter to provide statistical information across the spacecraft trajectory. Results: The study of different 3D shock models shows that the observations are compatible with a shock that is symmetric around the Sun-apex line as well as with an asymmetry up to an aspect ratio of around 3. Moreover, flat or dipped shock surfaces near their apex can only be rare cases. Next, the sheath thickness and the ICME velocity have no global trend along the ICME front. Finally, regrouping all these new results and those of our previous articles, we provide a quantitative ICME generic 3D shape, including the global shape of the shock, the sheath, and the flux rope. Conclusions: The obtained quantitative generic ICME shape will have implications for several aims. For example, it constrains the output of typical ICME numerical simulations. It is also a base for studying the transport of high-energy solar and cosmic particles during an ICME propagation as well as for modeling and forecasting space weather conditions near Earth.

GPU-accelerated 3-D model-based tracking

International Nuclear Information System (INIS)

Brown, J Anthony; Capson, David W

2010-01-01

Model-based approaches to tracking the pose of a 3-D object in video are effective but computationally demanding. While statistical estimation techniques, such as the particle filter, are often employed to minimize the search space, real-time performance remains unachievable on current generation CPUs. Recent advances in graphics processing units (GPUs) have brought massively parallel computational power to the desktop environment and powerful developer tools, such as NVIDIA Compute Unified Device Architecture (CUDA), have provided programmers with a mechanism to exploit it. NVIDIA GPUs' single-instruction multiple-thread (SIMT) programming model is well-suited to many computer vision tasks, particularly model-based tracking, which requires several hundred 3-D model poses to be dynamically configured, rendered, and evaluated against each frame in the video sequence. Using 6 degree-of-freedom (DOF) rigid hand tracking as an example application, this work harnesses consumer-grade GPUs to achieve real-time, 3-D model-based, markerless object tracking in monocular video.

Numerical simulation of the time-dependent deformation behaviour of clay-stone rock mass at the Tournemire site with 2D and 3D models

International Nuclear Information System (INIS)

Rutenberg, M.; Lux, K. H.

2011-01-01

Clay-stone rock masses are a reasonable alternative to e.g. salt rock masses as a host rock for underground radioactive waste repositories because of their very low permeability as well as their radionuclide retention capacity. Though clay-stone has been explored for many years, there is still a need for further research on its hydro-mechanical behaviour. Convergence measurements over a 4-year period in the tunnel system of the argillaceous Tournemire site in France yielded the presence of a time-dependent deformation behaviour in indurated clay. Moreover, a mine-by test was carried out with extensometer measurements capturing the rock mass deformation during the excavation process of a new gallery in 2003.This work focuses on the validation of a constitutive model by means of a three-dimensional (3D) simulation of the mine-by test. The utilised constitutive model Hou/Lux-T is based on the viscous constitutive model Lubby2 with which time-dependent deformation behaviour of salt rock can appropriately be simulated. It has been adapted to clay-stone by considering anisotropy effects, and in addition it features a strain-dependent fracture and failure criterion. The results of the mine-by-test simulation show that the calculated stresses and deformations in the rock mass seem to behave reasonably under this constitutive model with respect to time-dependency. A comparison of the 3D results to the results of a simplified two-dimensional (2D) simulation confirms the adequacy of using a 2D model with the constitutive model Hou/Lux-T for the setting at hand, described in the text (material parameters, time scale), in order to assess load-bearing capacity and deformability of the gallery near field away from heading face and tunnel crossing. Finally, a comparison of the 3D simulation results to the extensometer measurement results yields the principal ability of the used constitutive model to describe time-dependent evolutions of stresses and deformations during a three

NUMERICAL MODELLING OF THE SOIL BEHAVIOUR BY USING NEWLY DEVELOPED ADVANCED MATERIAL MODEL

Directory of Open Access Journals (Sweden)

Jan Veselý

2017-02-01

Full Text Available This paper describes a theoretical background, implementation and validation of the newly developed Jardine plastic hardening-softening model (JPHS model, which can be used for numerical modelling of the soils behaviour. Although the JPHS model is based on the elasto-plastic theory, like the Mohr-Coulomb model that is widely used in geotechnics, it contains some improvements, which removes the main disadvantages of the MC model. The presented model is coupled with an isotopically hardening and softening law, non-linear elastic stress-strain law, non-associated elasto-plastic material description and a cap yield surface. The validation of the model is done by comparing the numerical results with real measured data from the laboratory tests and by testing of the model on the real project of the tunnel excavation. The 3D numerical analysis is performed and the comparison between the JPHS, Mohr-Coulomb, Modified Cam-Clay, Hardening small strain model and monitoring in-situ data is done.

Three-dimensional direct numerical simulation of electromagnetically driven multiscale shallow layer flows: Numerical modeling and physical properties

Science.gov (United States)

Lardeau, Sylvain; Ferrari, Simone; Rossi, Lionel

2008-12-01

Three-dimensional (3D) direct numerical simulations of a flow driven by multiscale electromagnetic forcing are performed in order to reproduce with maximum accuracy the quasi-two-dimensional (2D) flow generated by the same multiscale forcing in the laboratory. The method presented is based on a 3D description of the flow and the electromagnetic forcing. Very good agreements between our simulations and the experiments are found both on velocity and acceleration field, this last comparison being, to our knowledge, done for the first time. Such agreement requires that both experiments and simulations are carefully performed and, more importantly, that the underlying simplification to model the experiments and the multiscale electromagnetic forcing do not introduce significant errors. The results presented in this paper differ significantly from previous 2D direct numerical simulation in which a classical linear Rayleigh friction modeling term was used to mimic the effect of the wall-normal friction. Indeed, purely 2D simulations are found to underestimate the Reynolds number and, due to the dominance of nonhomogeneous bottom friction, lead to the wrong physical mechanism. For the range of conditions presented in this paper, the Reynolds number, defined by the ratio between acceleration and viscous terms, remains the order of unity, and the Hartmann number, defined by the ratio between electromagnetic force terms and viscous terms, is about 2. The main conclusion is that 3D simulations are required to model the (3D) electromagnetic forces and the wall-normal shear. Indeed, even if the flow is quasi-2D in terms of energy, a full 3D approach is required to simulate these shallow layer flows driven by multiscale electromagnetic forcing. In the range of forcing intensity investigated in this paper, these multiscale flows remain quasi-2D, with negligible energy in the wall-normal velocity component. It is also shown that the driving terms are the electromagnetic forcing and

Multifractal modelling and 3D lacunarity analysis

International Nuclear Information System (INIS)

Hanen, Akkari; Imen, Bhouri; Asma, Ben Abdallah; Patrick, Dubois; Hedi, Bedoui Mohamed

2009-01-01

This study presents a comparative evaluation of lacunarity of 3D grey level models with different types of inhomogeneity. A new method based on the 'Relative Differential Box Counting' was developed to estimate the lacunarity features of grey level volumes. To validate our method, we generated a set of 3D grey level multifractal models with random, anisotropic and hierarchical properties. Our method gives a lacunarity measurement correlated with the theoretical one and allows a better model classification compared with a classical approach.

APPLICATION OF 3D MODELING IN 3D PRINTING FOR THE LOWER JAW RECONSTRUCTION

Directory of Open Access Journals (Sweden)

Yu. Yu. Dikov

2015-01-01

Full Text Available Aim of study: improvement of functional and aesthetic results of microsurgery reconstructions of the lower jaw due to the use of the methodology of 3D modeling and 3D printing. Application of this methodology has been demonstrated on the example of treatment of 4 patients with locally distributed tumors of the mouth cavity, who underwent excision of the tumor with simultaneous reconstruction of the lower jaw with revascularized fibular graft.Before, one patient has already undergo segmental resection of the lower jaw with the defect replacement with the avascular ileac graft and a reconstruction plate. Then, a relapse of the disease and lysis of the graft has developed with him. Modeling of the graft according to the shape of the lower jaw was performed by making osteotomies of the bone part of the graft using three-dimensional virtual models created by computed tomography data. Then these 3D models were printed with a 3D printer of plastic with the scale of 1:1 with the fused deposition modeling (FDM technology and were used during the surgery in the course of modeling of the graft. Sterilizing of the plastic model was performed in the formalin chamber.This methodology allowed more specific reconstruction of the resected fragment of the lower jaw and get better functional and aesthetic results and prepare patients to further dental rehabilitation. Advantages of this methodology are the possibility of simultaneous performance of stages of reconstruction and resection and shortening of the time of surgery.

Correlation between 2D and 3D flow curve modelling of DP steels using a microstructure-based RVE approach

International Nuclear Information System (INIS)

Ramazani, A.; Mukherjee, K.; Quade, H.; Prahl, U.; Bleck, W.

2013-01-01

A microstructure-based approach by means of representative volume elements (RVEs) is employed to evaluate the flow curve of DP steels using virtual tensile tests. Microstructures with different martensite fractions and morphologies are studied in two- and three-dimensional approaches. Micro sections of DP microstructures with various amounts of martensite have been converted to 2D RVEs, while 3D RVEs were constructed statistically with randomly distributed phases. A dislocation-based model is used to describe the flow curve of each ferrite and martensite phase separately as a function of carbon partitioning and microstructural features. Numerical tensile tests of RVE were carried out using the ABAQUS/Standard code to predict the flow behaviour of DP steels. It is observed that 2D plane strain modelling gives an underpredicted flow curve for DP steels, while the 3D modelling gives a quantitatively reasonable description of flow curve in comparison to the experimental data. In this work, a von Mises stress correlation factor σ 3D /σ 2D has been identified to compare the predicted flow curves of these two dimensionalities showing a third order polynomial relation with respect to martensite fraction and a second order polynomial relation with respect to equivalent plastic strain, respectively. The quantification of this polynomial correlation factor is performed based on laboratory-annealed DP600 chemistry with varying martensite content and it is validated for industrially produced DP qualities with various chemistry, strength level and martensite fraction.

Coupling 2D/3D registration method and statistical model to perform 3D reconstruction from partial x-rays images data.

Science.gov (United States)

Cresson, T; Chav, R; Branchaud, D; Humbert, L; Godbout, B; Aubert, B; Skalli, W; De Guise, J A

2009-01-01

3D reconstructions of the spine from a frontal and sagittal radiographs is extremely challenging. The overlying features of soft tissues and air cavities interfere with image processing. It is also difficult to obtain information that is accurate enough to reconstruct complete 3D models. To overcome these problems, the proposed method efficiently combines the partial information contained in two images from a patient with a statistical 3D spine model generated from a database of scoliotic patients. The algorithm operates through two simultaneous iterating processes. The first one generates a personalized vertebra model using a 2D/3D registration process with bone boundaries extracted from radiographs, while the other one infers the position and the shape of other vertebrae from the current estimation of the registration process using a statistical 3D model. Experimental evaluations have shown good performances of the proposed approach in terms of accuracy and robustness when compared to CT-scan.

Development of 3D statistical mandible models for cephalometric measurements

International Nuclear Information System (INIS)

Kim, Sung Goo; Yi, Won Jin; Hwang, Soon Jung; Choi, Soon Chul; Lee, Sam Sun; Heo, Min Suk; Huh, Kyung Hoe; Kim, Tae Il; Hong, Helen; Yoo, Ji Hyun

2012-01-01

The aim of this study was to provide sex-matched three-dimensional (3D) statistical shape models of the mandible, which would provide cephalometric parameters for 3D treatment planning and cephalometric measurements in orthognathic surgery. The subjects used to create the 3D shape models of the mandible included 23 males and 23 females. The mandibles were segmented semi-automatically from 3D facial CT images. Each individual mandible shape was reconstructed as a 3D surface model, which was parameterized to establish correspondence between different individual surfaces. The principal component analysis (PCA) applied to all mandible shapes produced a mean model and characteristic models of variation. The cephalometric parameters were measured directly from the mean models to evaluate the 3D shape models. The means of the measured parameters were compared with those from other conventional studies. The male and female 3D statistical mean models were developed from 23 individual mandibles, respectively. The male and female characteristic shapes of variation produced by PCA showed a large variability included in the individual mandibles. The cephalometric measurements from the developed models were very close to those from some conventional studies. We described the construction of 3D mandibular shape models and presented the application of the 3D mandibular template in cephalometric measurements. Optimal reference models determined from variations produced by PCA could be used for craniofacial patients with various types of skeletal shape.

Development of 3D statistical mandible models for cephalometric measurements

Energy Technology Data Exchange (ETDEWEB)

Kim, Sung Goo; Yi, Won Jin; Hwang, Soon Jung; Choi, Soon Chul; Lee, Sam Sun; Heo, Min Suk; Huh, Kyung Hoe; Kim, Tae Il [School of Dentistry, Seoul National University, Seoul (Korea, Republic of); Hong, Helen; Yoo, Ji Hyun [Division of Multimedia Engineering, Seoul Women' s University, Seoul (Korea, Republic of)

2012-09-15

The aim of this study was to provide sex-matched three-dimensional (3D) statistical shape models of the mandible, which would provide cephalometric parameters for 3D treatment planning and cephalometric measurements in orthognathic surgery. The subjects used to create the 3D shape models of the mandible included 23 males and 23 females. The mandibles were segmented semi-automatically from 3D facial CT images. Each individual mandible shape was reconstructed as a 3D surface model, which was parameterized to establish correspondence between different individual surfaces. The principal component analysis (PCA) applied to all mandible shapes produced a mean model and characteristic models of variation. The cephalometric parameters were measured directly from the mean models to evaluate the 3D shape models. The means of the measured parameters were compared with those from other conventional studies. The male and female 3D statistical mean models were developed from 23 individual mandibles, respectively. The male and female characteristic shapes of variation produced by PCA showed a large variability included in the individual mandibles. The cephalometric measurements from the developed models were very close to those from some conventional studies. We described the construction of 3D mandibular shape models and presented the application of the 3D mandibular template in cephalometric measurements. Optimal reference models determined from variations produced by PCA could be used for craniofacial patients with various types of skeletal shape.

Conjunction of 2D and 3D modified flow solvers for simulating spatio-temporal wind induced hydrodynamics in the Caspian Sea

Science.gov (United States)

Zounemat-Kermani, Mohammad; Sabbagh-Yazdi, Saeed-Reza

2010-06-01

The main objective of this study is the simulation of flow dynamics in the deep parts of the Caspian Sea, in which the southern and middle deep regions are surrounded by considerable areas of shallow zones. To simulate spatio-temporal wind induced hydrodynamics in deep waters, a conjunctive numerical model consisting of a 2D depth average model and a 3D pseudo compressible model is proposed. The 2D model is applied to determine time dependent free surface oscillations as well as the surface velocity patterns and is conjunct to the 3D flow solver for computing three-dimensional velocity and pressure fields which coverage to steady state for the top boundary condition. The modified 2D and 3D sets of equations are conjunct considering interface shear stresses. Both sets of 2D and 3D equations are solved on unstructured triangular and tetrahedral meshes using the Galerkin Finite Volume Method. The conjunctive model is utilized to investigate the deep currents affected by wind, Coriolis forces and the river inflow conditions of the Caspian Sea. In this study, the simulation of flow field due to major winds as well as transient winds in the Caspian Sea during a period of 6 hours in the winter season has been conducted and the numerical results for water surface level are then compared to the 2D numerical results.

An interactive display system for large-scale 3D models

Science.gov (United States)

Liu, Zijian; Sun, Kun; Tao, Wenbing; Liu, Liman

2018-04-01

With the improvement of 3D reconstruction theory and the rapid development of computer hardware technology, the reconstructed 3D models are enlarging in scale and increasing in complexity. Models with tens of thousands of 3D points or triangular meshes are common in practical applications. Due to storage and computing power limitation, it is difficult to achieve real-time display and interaction with large scale 3D models for some common 3D display software, such as MeshLab. In this paper, we propose a display system for large-scale 3D scene models. We construct the LOD (Levels of Detail) model of the reconstructed 3D scene in advance, and then use an out-of-core view-dependent multi-resolution rendering scheme to realize the real-time display of the large-scale 3D model. With the proposed method, our display system is able to render in real time while roaming in the reconstructed scene and 3D camera poses can also be displayed. Furthermore, the memory consumption can be significantly decreased via internal and external memory exchange mechanism, so that it is possible to display a large scale reconstructed scene with over millions of 3D points or triangular meshes in a regular PC with only 4GB RAM.

Hydrogeological conceptual model development and numerical modelling using CONNECTFLOW, Forsmark modelling stage 2.3

Energy Technology Data Exchange (ETDEWEB)

Follin, Sven (SF GeoLogic AB, Taeby (Sweden)); Hartley, Lee; Jackson, Peter; Roberts, David (Serco TAP (United Kingdom)); Marsic, Niko (Kemakta Konsult AB, Stockholm (Sweden))

2008-05-15

Three versions of a site descriptive model (SDM) have been completed for the Forsmark area. Version 0 established the state of knowledge prior to the start of the site investigation programme. Version 1.1 was essentially a training exercise and was completed during 2004. Version 1.2 was a preliminary site description and concluded the initial site investigation work (ISI) in June 2005. Three modelling stages are planned for the complete site investigation work (CSI). These are labelled stage 2.1, 2.2 and 2.3, respectively. An important component of each of these stages is to address and continuously try to resolve discipline-specific uncertainties of importance for repository engineering and safety assessment. Stage 2.1 included an updated geological model for Forsmark and aimed to provide a feedback from the modelling working group to the site investigation team to enable completion of the site investigation work. Stage 2.2 described the conceptual understanding and the numerical modelling of the bedrock hydrogeology in the Forsmark area based on data freeze 2.2. The present report describes the modelling based on data freeze 2.3, which is the final data freeze in Forsmark. In comparison, data freeze 2.3 is considerably smaller than data freeze 2.2. Therefore, stage 2.3 deals primarily with model confirmation and uncertainty analysis, e.g. verification of important hypotheses made in stage 2.2 and the role of parameter uncertainty in the numerical modelling. On the whole, the work reported here constitutes an addendum to the work reported in stage 2.2. Two changes were made to the CONNECTFLOW code in stage 2.3. These serve to: 1) improve the representation of the hydraulic properties of the regolith, and 2) improve the conditioning of transmissivity of the deformation zones against single-hole hydraulic tests. The changes to the modelling of the regolith were made to improve the consistency with models made with the MIKE SHE code, which involved the introduction

Multifractal modelling and 3D lacunarity analysis

Energy Technology Data Exchange (ETDEWEB)

Hanen, Akkari, E-mail: bettaieb.hanen@topnet.t [Laboratoire de biophysique, TIM, Faculte de Medecine (Tunisia); Imen, Bhouri, E-mail: bhouri_imen@yahoo.f [Unite de recherche ondelettes et multifractals, Faculte des sciences (Tunisia); Asma, Ben Abdallah, E-mail: asma.babdallah@cristal.rnu.t [Laboratoire de biophysique, TIM, Faculte de Medecine (Tunisia); Patrick, Dubois, E-mail: pdubois@chru-lille.f [INSERM, U 703, Lille (France); Hedi, Bedoui Mohamed, E-mail: medhedi.bedoui@fmm.rnu.t [Laboratoire de biophysique, TIM, Faculte de Medecine (Tunisia)

2009-09-28

This study presents a comparative evaluation of lacunarity of 3D grey level models with different types of inhomogeneity. A new method based on the 'Relative Differential Box Counting' was developed to estimate the lacunarity features of grey level volumes. To validate our method, we generated a set of 3D grey level multifractal models with random, anisotropic and hierarchical properties. Our method gives a lacunarity measurement correlated with the theoretical one and allows a better model classification compared with a classical approach.

Tangible 3D Modelling

DEFF Research Database (Denmark)

Hejlesen, Aske K.; Ovesen, Nis

2012-01-01

This paper presents an experimental approach to teaching 3D modelling techniques in an Industrial Design programme. The approach includes the use of tangible free form models as tools for improving the overall learning. The paper is based on lecturer and student experiences obtained through...... facilitated discussions during the course as well as through a survey distributed to the participating students. The analysis of the experiences shows a mixed picture consisting of both benefits and limits to the experimental technique. A discussion about the applicability of the technique and about...

Numerical modelling of coupled fluid, heat, and solute transport in deformable fractured rock

International Nuclear Information System (INIS)

Chan, T.; Reid, J.A.K.

1987-01-01

This paper reports on a three-dimensional (3D) finite-element code, MOTIF (model of transport in fractured/porous media), developed to model the coupled processes of groundwater flow, heat transport, brine transport, and one-species radionuclide transport in geological media. Three types of elements are available: a 3D continuum element, a planar fracture element that can be oriented in any arbitrary direction in 3D space or pipe flow in 3D space, and a line element for simulating fracture flow in 2D space or pipe flow in 3D space. As a quality-assurance measure, the MOTIF code was verified by comparison of its results with analytical solutions and other published numerical solutions

The Three-Dimensional (3D) Numerical Stability Analysis of Hyttemalmen Open-Pit

Science.gov (United States)

Cała, Marek; Kowalski, Michał; Stopkowicz, Agnieszka

2014-10-01

The purpose of this paper was to perform the 3D numerical calculations allowing slope stability analysis of Hyttemalmen open pit (location Kirkenes, Finnmark Province, Norway). After a ramp rock slide, which took place in December 2010, as well as some other small-scale rock slope stability problems, it proved necessary to perform a serious stability analyses. The Hyttemalmen open pit was designed with a depth up to 100 m, a bench height of 24 m and a ramp width of 10 m. The rock formation in the iron mining district of Kirkenes is called the Bjornevaten Group. This is the most structurally complicated area connected with tectonic process such as folding, faults and metamorphosis. The Bjornevaten Group is a volcano-sedimentary sequence. Rock slope stability depends on the mechanical properties of the rock, hydro-geological conditions, slope topography, joint set systems and seismic activity. However, rock slope stability is mainly connected with joint sets. Joints, or general discontinuities, are regarded as weak planes within rock which have strength reducing consequences with regard to rock strength. Discontinuities within the rock mass lead to very low tensile strength. Several simulations were performed utilising the RocLab (2007) software to estimate the gneiss cohesion for slopes of different height. The RocLab code is dedicated to estimate rock mass strength using the Hoek-Brown failure criterion. Utilising both the GSI index and the Hoek-Brown strength criterion the equivalent Mohr-Coulomb parameters (cohesion and angle of internal friction) can be calculated. The results of 3D numerical calculations (with FLA3D code) show that it is necessary to redesign the slope-bench system in the Hyttemalmen open pit. Changing slope inclination for lower stages is recommended. The minimum factor of safety should be equal 1.3. At the final planned stage of excavation, the factor of safety drops to 1.06 with failure surface ranging through all of the slopes. In the case

A continuous flow microfluidic calorimeter: 3-D numerical modeling with aqueous reactants

International Nuclear Information System (INIS)

Sen, Mehmet A.; Kowalski, Gregory J.; Fiering, Jason; Larson, Dale

2015-01-01

Highlights: • A co-flow microreactor is modeled in flow, reaction/diffusion, and thermal domains. • Analysis shows how arrayed temperature sensors can provide enthalpy of reaction. • Optical plasmonic temperature sensors could be arrayed suitably for calorimetry. • The reactor studied has a volume of 25 nL. - Abstract: A computational analysis of the reacting flow field, species diffusion and heat transfer processes with thermal boundary layer effects in a microchannel reactor with a coflow configuration was performed. Two parallel adjacent streams of aqueous reactants flow along a wide, shallow, enclosed channel in contact with a substrate, which is affixed to a temperature controlled plate. The Fluent computational fluid dynamics package solved the Navier–Stokes, mass transport and energy equations. The energy model, including the enthalpy of reaction as a nonuniform heat source, was validated by calculating the energy balance at several control volumes in the microchannel. Analysis reveals that the temperature is nearly uniform across the channel thickness, in the direction normal to the substrate surface; hence, measurements made by sensors at or near the surface are representative of the average temperature. Additionally, modeling the channel with a glass substrate and a silicone cover shows that heat transfer is predominantly due to the glass substrate. Finally, using the numerical results, we suggest that a microcalorimeter could be based on this configuration, and that temperature sensors such as optical nanohole array sensors could have sufficient spatial resolution to determine enthalpy of reaction

A continuous flow microfluidic calorimeter: 3-D numerical modeling with aqueous reactants

Energy Technology Data Exchange (ETDEWEB)

Sen, Mehmet A., E-mail: mehmet.sen@mathworks.com [Northeastern University, Department of Mechanical and Industrial Engineering, 360 Hungtington Avenue, 334 Snell Engineering Center, Boston, MA 02115 (United States); Kowalski, Gregory J., E-mail: gkowal@coe.neu.edu [Northeastern University, Department of Mechanical and Industrial Engineering, 360 Hungtington Avenue, 334 Snell Engineering Center, Boston, MA 02115 (United States); Fiering, Jason, E-mail: jfiering@draper.com [Charles Stark Draper Laboratory, 555 Technology Square, Cambridge, MA 02139 (United States); Larson, Dale, E-mail: dlarson@draper.com [Charles Stark Draper Laboratory, 555 Technology Square, Cambridge, MA 02139 (United States)

2015-03-10

Highlights: • A co-flow microreactor is modeled in flow, reaction/diffusion, and thermal domains. • Analysis shows how arrayed temperature sensors can provide enthalpy of reaction. • Optical plasmonic temperature sensors could be arrayed suitably for calorimetry. • The reactor studied has a volume of 25 nL. - Abstract: A computational analysis of the reacting flow field, species diffusion and heat transfer processes with thermal boundary layer effects in a microchannel reactor with a coflow configuration was performed. Two parallel adjacent streams of aqueous reactants flow along a wide, shallow, enclosed channel in contact with a substrate, which is affixed to a temperature controlled plate. The Fluent computational fluid dynamics package solved the Navier–Stokes, mass transport and energy equations. The energy model, including the enthalpy of reaction as a nonuniform heat source, was validated by calculating the energy balance at several control volumes in the microchannel. Analysis reveals that the temperature is nearly uniform across the channel thickness, in the direction normal to the substrate surface; hence, measurements made by sensors at or near the surface are representative of the average temperature. Additionally, modeling the channel with a glass substrate and a silicone cover shows that heat transfer is predominantly due to the glass substrate. Finally, using the numerical results, we suggest that a microcalorimeter could be based on this configuration, and that temperature sensors such as optical nanohole array sensors could have sufficient spatial resolution to determine enthalpy of reaction.

Investigation of Pumped Storage Hydropower Power-Off Transient Process Using 3D Numerical Simulation Based on SP-VOF Hybrid Model

Directory of Open Access Journals (Sweden)

Daqing Zhou

2018-04-01

Full Text Available The transient characteristic of the power-off process is investigated due to its close relation to hydraulic facilities’ safety in a pumped storage hydropower (PSH. In this paper, power-off transient characteristics of a PSH station in pump mode was studied using a three-dimensional (3D unsteady numerical method based on a single-phase and volume of fluid (SP-VOF coupled model. The computational domain covered the entire flow system, including reservoirs, diversion tunnel, surge tank, pump-turbine unit, and tailrace tunnel. The fast changing flow fields and dynamic characteristic parameters, such as unit flow rate, runner rotate speed, pumping lift, and static pressure at measuring points were simulated, and agreed well with experimental results. During the power-off transient process, the PSH station underwent pump mode, braking mode, and turbine mode, with the dynamic characteristics and inner flow configurations changing significantly. Intense pressure fluctuation occurred in the region between the runner and guide vanes, and its frequency and amplitude were closely related to the runner’s rotation speed and pressure gradient, respectively. While the reversed flow rate of the PSH unit reached maximum, some parameters, such as static pressure, torque, and pumping lift would suddenly jump significantly, due to the water hammer effect. The moment these marked jumps occurred was commonly considered as the most dangerous moment during the power-off transient process, due to the blade passages being clogged by vortexes, and chaos pressure distribution on the blade surfaces. The results of this study confirm that 3D SP-VOF hybrid simulation is an effective method to reveal the hydraulic mechanism of the PSH transient process.

Impact of aircraft exhaust on the atmosphere. Box model studies and 3-D mesoscale numerical case studies of seasonal differences

Energy Technology Data Exchange (ETDEWEB)

Petry, H; Ebel, A; Franzkowiak, V; Hendricks, J; Lippert, E; Moellhoff, M [Koeln Univ. (Germany). Inst. fuer Geophysik und Meteorologie

1998-12-31

The impact of aircraft emissions released in the tropopause region on atmospheric trace gases as O{sub 3} or HNO{sub 3} is investigated by means of model studies. Special emphasis is drawn on seasonal effects. A box model is applied as well as a 3-D mesoscale chemistry transport model. These model studies show that the impact of aircraft emissions on ozone in the tropopause region is much stronger in summer than in late autumn with a difference of one order of magnitude. (author) 14 refs.

Impact of aircraft exhaust on the atmosphere. Box model studies and 3-D mesoscale numerical case studies of seasonal differences

Energy Technology Data Exchange (ETDEWEB)

Petry, H.; Ebel, A.; Franzkowiak, V.; Hendricks, J.; Lippert, E.; Moellhoff, M. [Koeln Univ. (Germany). Inst. fuer Geophysik und Meteorologie

1997-12-31

The impact of aircraft emissions released in the tropopause region on atmospheric trace gases as O{sub 3} or HNO{sub 3} is investigated by means of model studies. Special emphasis is drawn on seasonal effects. A box model is applied as well as a 3-D mesoscale chemistry transport model. These model studies show that the impact of aircraft emissions on ozone in the tropopause region is much stronger in summer than in late autumn with a difference of one order of magnitude. (author) 14 refs.

High-accuracy critical exponents for O(N) hierarchical 3D sigma models

International Nuclear Information System (INIS)

Godina, J. J.; Li, L.; Meurice, Y.; Oktay, M. B.

2006-01-01

The critical exponent γ and its subleading exponent Δ in the 3D O(N) Dyson's hierarchical model for N up to 20 are calculated with high accuracy. We calculate the critical temperatures for the measure δ(φ-vector.φ-vector-1). We extract the first coefficients of the 1/N expansion from our numerical data. We show that the leading and subleading exponents agree with Polchinski equation and the equivalent Litim equation, in the local potential approximation, with at least 4 significant digits

PDE-based geophysical modelling using finite elements: examples from 3D resistivity and 2D magnetotellurics

International Nuclear Information System (INIS)

Schaa, R; Gross, L; Du Plessis, J

2016-01-01

We present a general finite-element solver, escript, tailored to solve geophysical forward and inverse modeling problems in terms of partial differential equations (PDEs) with suitable boundary conditions. Escript’s abstract interface allows geoscientists to focus on solving the actual problem without being experts in numerical modeling. General-purpose finite element solvers have found wide use especially in engineering fields and find increasing application in the geophysical disciplines as these offer a single interface to tackle different geophysical problems. These solvers are useful for data interpretation and for research, but can also be a useful tool in educational settings. This paper serves as an introduction into PDE-based modeling with escript where we demonstrate in detail how escript is used to solve two different forward modeling problems from applied geophysics (3D DC resistivity and 2D magnetotellurics). Based on these two different cases, other geophysical modeling work can easily be realized. The escript package is implemented as a Python library and allows the solution of coupled, linear or non-linear, time-dependent PDEs. Parallel execution for both shared and distributed memory architectures is supported and can be used without modifications to the scripts. (paper)

PDE-based geophysical modelling using finite elements: examples from 3D resistivity and 2D magnetotellurics

Science.gov (United States)

Schaa, R.; Gross, L.; du Plessis, J.

2016-04-01

We present a general finite-element solver, escript, tailored to solve geophysical forward and inverse modeling problems in terms of partial differential equations (PDEs) with suitable boundary conditions. Escript’s abstract interface allows geoscientists to focus on solving the actual problem without being experts in numerical modeling. General-purpose finite element solvers have found wide use especially in engineering fields and find increasing application in the geophysical disciplines as these offer a single interface to tackle different geophysical problems. These solvers are useful for data interpretation and for research, but can also be a useful tool in educational settings. This paper serves as an introduction into PDE-based modeling with escript where we demonstrate in detail how escript is used to solve two different forward modeling problems from applied geophysics (3D DC resistivity and 2D magnetotellurics). Based on these two different cases, other geophysical modeling work can easily be realized. The escript package is implemented as a Python library and allows the solution of coupled, linear or non-linear, time-dependent PDEs. Parallel execution for both shared and distributed memory architectures is supported and can be used without modifications to the scripts.

On modeling of statistical properties of classical 3D spin glasses

International Nuclear Information System (INIS)

Gevorkyan, A.S.; Abajyan, H.G.; Ayryan, E.A.

2011-01-01

We study statistical properties of 3D classical spin glass layer of certain width and infinite length. The 3D spin glass is represented as an ensemble of disordered 1D spatial spin chains (SSC) where interactions are random between spin chains (nonideal ensemble of 1D SSCs). It is proved that in the limit of Birkhoff's ergodic hypothesis performance, 3D spin glasses can be generated by Hamiltonian of disordered 1D SSC with random environment. Disordered 1D SSC is defined on a regular lattice where one randomly oriented spin is put on each node of lattice. Also, it is supposed that each spin randomly interacts with six nearest-neighboring spins (two spins on lattice and four in the environment). The recurrent transcendental equations are obtained on the nodes of spin-chain lattice. These equations, combined with the Silvester conditions, allow step-by-step construction of spin chain in the ground state of energy where all spins are in the minimal energy of a classical Hamiltonian. On the basis of these equations an original high-performance parallel algorithm is developed for 3D spin glasses simulation. Distributions of different parameters of unperturbed spin glass are calculated. In particular, it is analytically proved and numerical calculations show that the distribution of spin-spin interaction constant in Heisenberg nearest-neighboring Hamiltonian model, as opposed to widely used Gauss-Edwards-Anderson distribution, satisfies the Levy alpha-stable distribution law which does not have variance. A new formula is proposed for construction of partition function in the form of a one-dimensional integral on the energy distribution of 1D SSCs

3D thermal model of laser surface glazing for H13 tool steel

Science.gov (United States)

Kabir, I. R.; Yin, D.; Naher, S.

2017-10-01

In this work a three dimensional (3D) finite element model of laser surface glazing (LSG) process has been developed. The purpose of the 3D thermal model of LSG was to achieve maximum accuracy towards the predicted outcome for optimizing the process. A cylindrical geometry of 10mm diameter and 1mm length was used in ANSYS 15 software. Temperature distribution, depth of modified zone and cooling rates were analysed from the thermal model. Parametric study was carried out varying the laser power from 200W-300W with constant beam diameter and residence time which were 0.2mm and 0.15ms respectively. The maximum surface temperature 2554°K was obtained for power 300W and minimum surface temperature 1668°K for power 200W. Heating and cooling rates increased with increasing laser power. The depth of the laser modified zone attained for 300W power was 37.5µm and for 200W power was 30µm. No molten zone was observed at 200W power. Maximum surface temperatures obtained from 3D model increased 4% than 2D model presented in author's previous work. In order to verify simulation results an analytical solution of temperature distribution for laser surface modification was used. The surface temperature after heating was calculated for similar laser parameters which is 1689°K. The difference in maximum surface temperature is around 20.7°K between analytical and numerical analysis of LSG for power 200W.

Animation of 3D Model of Human Head

Directory of Open Access Journals (Sweden)

V. Michalcin

2007-04-01

Full Text Available The paper deals with the new algorithm of animation of 3D model of the human head in combination with its global motion. The designed algorithm is very fast and with low calculation requirements, because it does not need the synthesis of the input videosequence for estimation of the animation parameters as well as the parameters of global motion. The used 3D model Candide generates different expressions using its animation units which are controlled by the animation parameters. These ones are estimated on the basis of optical flow without the need of extracting of the feature points in the frames of the input videosequence because they are given by the selected vertices of the animation units of the calibrated 3D model Candide. The established multiple iterations inside the designed animation algorithm of 3D model of the human head between two successive frames significantly improved its accuracy above all for the large motion.

UCVM: An Open Source Framework for 3D Velocity Model Research

Science.gov (United States)

Gill, D.; Maechling, P. J.; Jordan, T. H.; Plesch, A.; Taborda, R.; Callaghan, S.; Small, P.

2013-12-01

Three-dimensional (3D) seismic velocity models provide fundamental input data to ground motion simulations, in the form of structured or unstructured meshes or grids. Numerous models are available for California, as well as for other parts of the United States and Europe, but models do not share a common interface. Being able to interact with these models in a standardized way is critical in order to configure and run 3D ground motion simulations. The Unified Community Velocity Model (UCVM) software, developed by researchers at the Southern California Earthquake Center (SCEC), is an open source framework designed to provide a cohesive way to interact with seismic velocity models. We describe the several ways in which we have improved the UCVM software over the last year. We have simplified the UCVM installation process by automating the installation of various community codebases, improving the ease of use.. We discuss how UCVM software was used to build velocity meshes for high-frequency (4Hz) deterministic 3D wave propagation simulations, and how the UCVM framework interacts with other open source resources, such as NetCDF file formats for visualization. The UCVM software uses a layered software architecture that transparently converts geographic coordinates to the coordinate systems used by the underlying velocity models and supports inclusion of a configurable near-surface geotechnical layer, while interacting with the velocity model codes through their existing software interfaces. No changes to the velocity model codes are required. Our recent UCVM installation improvements bundle UCVM with a setup script, written in Python, which guides users through the process that installs the UCVM software along with all the user-selectable velocity models. Each velocity model is converted into a standardized (configure, make, make install) format that is easily downloaded and installed via the script. UCVM is often run in specialized high performance computing (HPC

Comparison of PIV measurements and a discrete particle model in a rectangular 3D spout-fluid bed

NARCIS (Netherlands)

Link, J.M.; Deen, N.G.; Kuipers, J.A.M.

2004-01-01

Particle image velocimetry and a 3D hard sphere discrete particle model were applied to determine particle velocity profiles in the plane around a spout in a spoutfluid bed for various initial bed heights, spout and background fluidization velocities. Comparison between experimental and numerical

Understanding the seismic wave propagation inside and around an underground cavity from a 3D numerical survey

Science.gov (United States)

Esterhazy, Sofi; Schneider, Felix; Perugia, Ilaria; Bokelmann, Götz

2017-04-01

Motivated by the need to detect an underground cavity within the procedure of an On-Site-Inspection (OSI) of the Comprehensive Nuclear Test Ban Treaty Organization (CTBTO), which might be caused by a nuclear explosion/weapon testing, we aim to provide a basic numerical study of the wave propagation around and inside such an underground cavity. One method to investigate the geophysical properties of an underground cavity allowed by the Comprehensive Nuclear-test Ban Treaty is referred to as "resonance seismometry" - a resonance method that uses passive or active seismic techniques, relying on seismic cavity vibrations. This method is in fact not yet entirely determined by the Treaty and so far, there are only very few experimental examples that have been suitably documented to build a proper scientific groundwork. This motivates to investigate this problem on a purely numerical level and to simulate these events based on recent advances in numerical modeling of wave propagation problems. Our numerical study includes the full elastic wave field in three dimensions. We consider the effects from an incoming plane wave as well as point source located in the surrounding of the cavity at the surface. While the former can be considered as passive source like a tele-seismic earthquake, the latter represents a man-made explosion or a viborseis as used for/in active seismic techniques. Further we want to demonstrate the specific characteristics of the scattered wave field from a P-waves and S-wave separately. For our simulations in 3D we use the discontinuous Galerkin Spectral Element Code SPEED developed by MOX (The Laboratory for Modeling and Scientific Computing, Department of Mathematics) and DICA (Department of Civil and Environmental Engineering) at the Politecnico di Milano. The computations are carried out on the Vienna Scientific Cluster (VSC). The accurate numerical modeling can facilitate the development of proper analysis techniques to detect the remnants of an

3D MODELING OF THE ARCHAIC AMPHORAS OF IONIA

Directory of Open Access Journals (Sweden)

A. Denker

2015-04-01

Full Text Available Few other regions offer such a rich collection of amphoras than the cities of Ionia. Throughout history amphoras of these cities had been spread all over the Mediterranean. Despite their common characteristics, amphora manufacturing cities of Ionia had their own distinctive styles that can be identified. They differed in details of shape and decoration. Each city produced an authentic type of amphora which served as a trademark of itself and enabled its attribution to where it originated from. That’s why, amphoras provide important insight into commerce of old ages and yield evidence into ancient sailing routes. Owing to this our knowledge of the ancient trade is profoundly enriched. The following is based on the finds of amphoras which originated from the Ionian cities of Chios, Clazomenai, Lesbos, Miletus, and Samos. Starting from city-specific forms which offer interpretative advantages in provenancing, this article surveys the salient features of the regional forms and styles of the those Ionian cities. 3D modeling is utilized with the aim of bringing fresh glimpses of the investigated amphoras by showing how they originally looked. Due to their virtual indestructibility these models offer interpretative advantages by enabling experimental testing of hypotheses upon the finds without risking them. The 3D models in the following sections were reconstructed from numerous fragments of necks, handles, body sherds and bases. They convey in color- unlike the monochrome drawings which we were accustomed to-the texture, decoration, tint and the vitality of the amphoras of Ionia.

MULTI SENSOR DATA INTEGRATION FOR AN ACCURATE 3D MODEL GENERATION

Directory of Open Access Journals (Sweden)

S. Chhatkuli

2015-05-01

Full Text Available The aim of this paper is to introduce a novel technique of data integration between two different data sets, i.e. laser scanned RGB point cloud and oblique imageries derived 3D model, to create a 3D model with more details and better accuracy. In general, aerial imageries are used to create a 3D city model. Aerial imageries produce an overall decent 3D city models and generally suit to generate 3D model of building roof and some non-complex terrain. However, the automatically generated 3D model, from aerial imageries, generally suffers from the lack of accuracy in deriving the 3D model of road under the bridges, details under tree canopy, isolated trees, etc. Moreover, the automatically generated 3D model from aerial imageries also suffers from undulated road surfaces, non-conforming building shapes, loss of minute details like street furniture, etc. in many cases. On the other hand, laser scanned data and images taken from mobile vehicle platform can produce more detailed 3D road model, street furniture model, 3D model of details under bridge, etc. However, laser scanned data and images from mobile vehicle are not suitable to acquire detailed 3D model of tall buildings, roof tops, and so forth. Our proposed approach to integrate multi sensor data compensated each other’s weakness and helped to create a very detailed 3D model with better accuracy. Moreover, the additional details like isolated trees, street furniture, etc. which were missing in the original 3D model derived from aerial imageries could also be integrated in the final model automatically. During the process, the noise in the laser scanned data for example people, vehicles etc. on the road were also automatically removed. Hence, even though the two dataset were acquired in different time period the integrated data set or the final 3D model was generally noise free and without unnecessary details.

Numerical study of heat and mass transfer optimization in a 3D inclined solar distiller

Directory of Open Access Journals (Sweden)

Ghachem Kaouther

2017-01-01

Full Text Available A numerical study of the 3-D double-diffusive natural convection in an inclined solar distiller was established. The flow is considered laminar and caused by the interaction of thermal energy and the chemical species diffusions. The governing equations of the problem, are formulated using vector potential-vorticity formalism in its 3-D form, then solved by the finite volumes method. The Rayleigh number is fixed at Ra = 105 and effects of the buoyancy ratio and inclination are studied for opposed temperature and concentration gradients. The main purpose of the study is to find the optimum inclination angle of the distiller which promotes the maximum mass and heat transfer.

Computer numerical control (CNC) lithography: light-motion synchronized UV-LED lithography for 3D microfabrication

International Nuclear Information System (INIS)

Kim, Jungkwun; Allen, Mark G; Yoon, Yong-Kyu

2016-01-01

This paper presents a computer-numerical-controlled ultraviolet light-emitting diode (CNC UV-LED) lithography scheme for three-dimensional (3D) microfabrication. The CNC lithography scheme utilizes sequential multi-angled UV light exposures along with a synchronized switchable UV light source to create arbitrary 3D light traces, which are transferred into the photosensitive resist. The system comprises a switchable, movable UV-LED array as a light source, a motorized tilt-rotational sample holder, and a computer-control unit. System operation is such that the tilt-rotational sample holder moves in a pre-programmed routine, and the UV-LED is illuminated only at desired positions of the sample holder during the desired time period, enabling the formation of complex 3D microstructures. This facilitates easy fabrication of complex 3D structures, which otherwise would have required multiple manual exposure steps as in the previous multidirectional 3D UV lithography approach. Since it is batch processed, processing time is far less than that of the 3D printing approach at the expense of some reduction in the degree of achievable 3D structure complexity. In order to produce uniform light intensity from the arrayed LED light source, the UV-LED array stage has been kept rotating during exposure. UV-LED 3D fabrication capability was demonstrated through a plurality of complex structures such as V-shaped micropillars, micropanels, a micro-‘hi’ structure, a micro-‘cat’s claw,’ a micro-‘horn,’ a micro-‘calla lily,’ a micro-‘cowboy’s hat,’ and a micro-‘table napkin’ array. (paper)

3D Numerical Model of Continental Breakup via Plume Lithosphere Interaction Near Cratonic Blocks: Implications for the Tanzanian Craton

Science.gov (United States)

Koptev, A.; Calais, E.; Burov, E. B.; Leroy, S. D.; Gerya, T.

2014-12-01

Although many continental rift basins and their successfully rifted counterparts at passive continental margins are magmatic, some are not. This dichotomy prompted end-member views of the mechanism driving continental rifting, deep-seated and mantle plume-driven for some, owing to shallow lithospheric stretching for others. In that regard, the East African Rift (EAR), the 3000 km-long divergent boundary between the Nubian and Somalian plates, provides a unique setting with the juxtaposition of the eastern, magma-rich, and western, magma-poor, branches on either sides of the 250-km thick Tanzanian craton. Here we implement high-resolution rheologically realistic 3D numerical model of plume-lithosphere interactions in extensional far-field settings to explain this contrasted behaviour in a unified framework starting from simple, symmetrical initial conditions with an isolated mantle plume rising beneath a craton in an east-west tensional far field stress. The upwelling mantle plume is deflected by the cratonic keel and preferentially channelled along one of its sides. This leads to the coeval development of a magma-rich branch above the plume head and a magma-poor one along the opposite side of the craton, the formation of a rotating microplate between the two rift branches, and the feeding of melt to both branches form a single mantle source. The model bears strong similarities with the evolution of the eastern and western branches of the central EAR and the geodetically observed rotation of the Victoria microplate. This result reconciles the passive (plume-activated) versus active (far-field tectonic stresses) rift models as our experiments shows both processes in action and demonstrate the possibility of developing both magmatic and amagmatic rifts in identical geotectonic environments.

TLS for generating multi-LOD of 3D building model

International Nuclear Information System (INIS)

Akmalia, R; Setan, H; Majid, Z; Suwardhi, D; Chong, A

2014-01-01

The popularity of Terrestrial Laser Scanners (TLS) to capture three dimensional (3D) objects has been used widely for various applications. Development in 3D models has also led people to visualize the environment in 3D. Visualization of objects in a city environment in 3D can be useful for many applications. However, different applications require different kind of 3D models. Since a building is an important object, CityGML has defined a standard for 3D building models at four different levels of detail (LOD). In this research, the advantages of TLS for capturing buildings and the modelling process of the point cloud can be explored. TLS will be used to capture all the building details to generate multi-LOD. This task, in previous works, involves usually the integration of several sensors. However, in this research, point cloud from TLS will be processed to generate the LOD3 model. LOD2 and LOD1 will then be generalized from the resulting LOD3 model. Result from this research is a guiding process to generate the multi-LOD of 3D building starting from LOD3 using TLS. Lastly, the visualization for multi-LOD model will also be shown

TLS for generating multi-LOD of 3D building model

Science.gov (United States)

Akmalia, R.; Setan, H.; Majid, Z.; Suwardhi, D.; Chong, A.

2014-02-01

The popularity of Terrestrial Laser Scanners (TLS) to capture three dimensional (3D) objects has been used widely for various applications. Development in 3D models has also led people to visualize the environment in 3D. Visualization of objects in a city environment in 3D can be useful for many applications. However, different applications require different kind of 3D models. Since a building is an important object, CityGML has defined a standard for 3D building models at four different levels of detail (LOD). In this research, the advantages of TLS for capturing buildings and the modelling process of the point cloud can be explored. TLS will be used to capture all the building details to generate multi-LOD. This task, in previous works, involves usually the integration of several sensors. However, in this research, point cloud from TLS will be processed to generate the LOD3 model. LOD2 and LOD1 will then be generalized from the resulting LOD3 model. Result from this research is a guiding process to generate the multi-LOD of 3D building starting from LOD3 using TLS. Lastly, the visualization for multi-LOD model will also be shown.

Modelling reinforced concrete structures in DYNA3D

International Nuclear Information System (INIS)

Broadhouse, B.J.; Neilson, A.J.

1987-10-01

A material model for reinforced concrete has been implemented in the transient structural dynamics code DYNA3D. This paper outlines the constitutive material model, and presents comparisons of DYNA3D calculations and experiments on impulsively loaded panels, covering the full range of panel damage states from light cracking through to panel collapse or perforation. The results are presented using the post-processor code TAURUS, which has also been modified to provide mesh diagrams with superimposed crack patterns from the DYNA3D predictions. (author)

3D MHD Flux emergence experiments

DEFF Research Database (Denmark)

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

2012-01-01

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