Finite-difference solutions of the 3-D eikonal equation
Fei, Tong; Fehler, M.C.; Hildebrand, S.T.
1995-12-31
Prestack Kirchhoff depth migration requires the computation of traveltimes from surface source and receiver locations to subsurface image locations. In 3-D problems, computational efficiency becomes important. Finite-difference solutions of the eikonal equation provide computationally efficient methods for generating the traveltime information. Here, a novel finite-difference solutions of the eikonal equation provide computationally efficient methods for generating the traveltime information. Here, a novel finite-difference method for computing the first arrival traveltime by solving the eikonal equation has been developed in Cartesian coordinates. The method, which is unconditionally stable and computationally efficient, can handle instabilities due to caustics and provide information about head waves. The comparison of finite-difference solutions of the acoustic wave equation with the traveltime solutions from the eikonal equation in various structure models demonstrate that the method developed here can provide correct first arrival traveltime information even in areas of complex velocity structure.
3D finite-difference seismic migration with parallel computers
Ober, C.C.; Gjertsen, R.; Minkoff, S.; Womble, D.E.
1998-11-01
The ability to image complex geologies such as salt domes in the Gulf of Mexico and thrusts in mountainous regions is essential for reducing the risk associated with oil exploration. Imaging these structures, however, is computationally expensive as datasets can be terabytes in size. Traditional ray-tracing migration methods cannot handle complex velocity variations commonly found near such salt structures. Instead the authors use the full 3D acoustic wave equation, discretized via a finite difference algorithm. They reduce the cost of solving the apraxial wave equation by a number of numerical techniques including the method of fractional steps and pipelining the tridiagonal solves. The imaging code, Salvo, uses both frequency parallelism (generally 90% efficient) and spatial parallelism (65% efficient). Salvo has been tested on synthetic and real data and produces clear images of the subsurface even beneath complicated salt structures.
3D Finite Difference Modelling of Basaltic Region
NASA Astrophysics Data System (ADS)
Engell-Sørensen, L.
2003-04-01
The main purpose of the work was to generate realistic data to be applied for testing of processing and migration tools for basaltic regions. The project is based on the three - dimensional finite difference code (FD), TIGER, made by Sintef. The FD code was optimized (parallelized) by the author, to run on parallel computers. The parallel code enables us to model large-scale realistic geological models and to apply traditional seismic and micro seismic sources. The parallel code uses multiple processors in order to manipulate subsets of large amounts of data simultaneously. The general anisotropic code uses 21 elastic coefficients. Eight independent coefficients are needed as input parameters for the general TI medium. In the FD code, the elastic wave field computation is implemented by a higher order FD solution to the elastic wave equation and the wave fields are computed on a staggered grid, shifted half a node in one or two directions. The geological model is a gridded basalt model, which covers from 24 km to 37 km of a real shot line in horizontal direction and from the water surface to the depth of 3.5 km. The 2frac {1}{2}D model has been constructed using the compound modeling software from Norsk Hydro. The vertical parameter distribution is obtained from observations in two wells. At The depth of between 1100 m to 1500 m, a basalt horizon covers the whole sub surface layers. We have shown that it is possible to simulate a line survey in realistic (3D) geological models in reasonable time by using high performance computers. The author would like to thank Norsk Hydro, Statoil, GEUS, and SINTEF for very helpful discussions and Parallab for being helpful with the new IBM, p690 Regatta system.
Spatial parallelism of a 3D finite difference, velocity-stress elastic wave propagation code
Minkoff, S.E.
1999-12-01
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately, finite difference simulations for 3D elastic wave propagation are expensive. The authors model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MPI library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speedup. Because I/O is handled largely outside of the time-step loop (the most expensive part of the simulation) the authors have opted for straight-forward broadcast and reduce operations to handle I/O. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ghost cells. When this communication is balanced against computation by allocating subdomains of reasonable size, they observe excellent scaled speedup. Allocating subdomains of size 25 x 25 x 25 on each node, they achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
Spatial Parallelism of a 3D Finite Difference, Velocity-Stress Elastic Wave Propagation Code
MINKOFF,SUSAN E.
1999-12-09
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately. finite difference simulations for 3D elastic wave propagation are expensive. We model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MP1 library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speed up. Because i/o is handled largely outside of the time-step loop (the most expensive part of the simulation) we have opted for straight-forward broadcast and reduce operations to handle i/o. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ''ghost cells''. When this communication is balanced against computation by allocating subdomains of reasonable size, we observe excellent scaled speed up. Allocating subdomains of size 25 x 25 x 25 on each node, we achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
Ground motion simulations in Marmara (Turkey) region from 3D finite difference method
NASA Astrophysics Data System (ADS)
Aochi, Hideo; Ulrich, Thomas; Douglas, John
2016-04-01
In the framework of the European project MARSite (2012-2016), one of the main contributions from our research team was to provide ground-motion simulations for the Marmara region from various earthquake source scenarios. We adopted a 3D finite difference code, taking into account the 3D structure around the Sea of Marmara (including the bathymetry) and the sea layer. We simulated two moderate earthquakes (about Mw4.5) and found that the 3D structure improves significantly the waveforms compared to the 1D layer model. Simulations were carried out for different earthquakes (moderate point sources and large finite sources) in order to provide shake maps (Aochi and Ulrich, BSSA, 2015), to study the variability of ground-motion parameters (Douglas & Aochi, BSSA, 2016) as well as to provide synthetic seismograms for the blind inversion tests (Diao et al., GJI, 2016). The results are also planned to be integrated in broadband ground-motion simulations, tsunamis generation and simulations of triggered landslides (in progress by different partners). The simulations are freely shared among the partners via the internet and the visualization of the results is diffused on the project's homepage. All these simulations should be seen as a reference for this region, as they are based on the latest knowledge that obtained during the MARSite project, although their refinement and validation of the model parameters and the simulations are a continuing research task relying on continuing observations. The numerical code used, the models and the simulations are available on demand.
GPU-accelerated 3D neutron diffusion code based on finite difference method
Xu, Q.; Yu, G.; Wang, K.
2012-07-01
Finite difference method, as a traditional numerical solution to neutron diffusion equation, although considered simpler and more precise than the coarse mesh nodal methods, has a bottle neck to be widely applied caused by the huge memory and unendurable computation time it requires. In recent years, the concept of General-Purpose computation on GPUs has provided us with a powerful computational engine for scientific research. In this study, a GPU-Accelerated multi-group 3D neutron diffusion code based on finite difference method was developed. First, a clean-sheet neutron diffusion code (3DFD-CPU) was written in C++ on the CPU architecture, and later ported to GPUs under NVIDIA's CUDA platform (3DFD-GPU). The IAEA 3D PWR benchmark problem was calculated in the numerical test, where three different codes, including the original CPU-based sequential code, the HYPRE (High Performance Pre-conditioners)-based diffusion code and CITATION, were used as counterpoints to test the efficiency and accuracy of the GPU-based program. The results demonstrate both high efficiency and adequate accuracy of the GPU implementation for neutron diffusion equation. A speedup factor of about 46 times was obtained, using NVIDIA's Geforce GTX470 GPU card against a 2.50 GHz Intel Quad Q9300 CPU processor. Compared with the HYPRE-based code performing in parallel on an 8-core tower server, the speedup of about 2 still could be observed. More encouragingly, without any mathematical acceleration technology, the GPU implementation ran about 5 times faster than CITATION which was speeded up by using the SOR method and Chebyshev extrapolation technique. (authors)
3D frequency-domain finite-difference modeling of acoustic wave propagation
NASA Astrophysics Data System (ADS)
Operto, S.; Virieux, J.
2006-12-01
We present a 3D frequency-domain finite-difference method for acoustic wave propagation modeling. This method is developed as a tool to perform 3D frequency-domain full-waveform inversion of wide-angle seismic data. For wide-angle data, frequency-domain full-waveform inversion can be applied only to few discrete frequencies to develop reliable velocity model. Frequency-domain finite-difference (FD) modeling of wave propagation requires resolution of a huge sparse system of linear equations. If this system can be solved with a direct method, solutions for multiple sources can be computed efficiently once the underlying matrix has been factorized. The drawback of the direct method is the memory requirement resulting from the fill-in of the matrix during factorization. We assess in this study whether representative problems can be addressed in 3D geometry with such approach. We start from the velocity-stress formulation of the 3D acoustic wave equation. The spatial derivatives are discretized with second-order accurate staggered-grid stencil on different coordinate systems such that the axis span over as many directions as possible. Once the discrete equations were developed on each coordinate system, the particle velocity fields are eliminated from the first-order hyperbolic system (following the so-called parsimonious staggered-grid method) leading to second-order elliptic wave equations in pressure. The second-order wave equations discretized on each coordinate system are combined linearly to mitigate the numerical anisotropy. Secondly, grid dispersion is minimized by replacing the mass term at the collocation point by its weighted averaging over all the grid points of the stencil. Use of second-order accurate staggered- grid stencil allows to reduce the bandwidth of the matrix to be factorized. The final stencil incorporates 27 points. Absorbing conditions are PML. The system is solved using the parallel direct solver MUMPS developed for distributed
Acceleration of 3D Finite Difference AWP-ODC for seismic simulation on GPU Fermi Architecture
NASA Astrophysics Data System (ADS)
Zhou, J.; Cui, Y.; Choi, D.
2011-12-01
AWP-ODC, a highly scalable parallel finite-difference application, enables petascale 3D earthquake calculations. This application generates realistic dynamic earthquake source description and detailed physics-based anelastic ground motions at frequencies pertinent to safe building design. In 2010, the code achieved M8, a full dynamical simulation of a magnitude-8 earthquake on the southern San Andreas fault up to 2-Hz, the largest-ever earthquake simulation. Building on the success of the previous work, we have implemented CUDA on AWP-ODC to accelerate wave propagation on GPU platform. Our CUDA development aims on aggressive parallel efficiency, optimized global and shared memory access to make the best use of GPU memory hierarchy. The benchmark on NVIDIA Tesla C2050 graphics cards demonstrated many tens of speedup in single precision compared to serial implementation at a testing problem size, while an MPI-CUDA implementation is in the progress to extend our solver to multi-GPU clusters. Our CUDA implementation has been carefully verified for accuracy.
Rigorous interpolation near tilted interfaces in 3-D finite-difference EM modelling
NASA Astrophysics Data System (ADS)
Shantsev, Daniil V.; Maaø, Frank A.
2015-02-01
We present a rigorous method for interpolation of electric and magnetic fields close to an interface with a conductivity contrast. The method takes into account not only a well-known discontinuity in the normal electric field, but also discontinuity in all the normal derivatives of electric and magnetic tangential fields. The proposed method is applied to marine 3-D controlled-source electromagnetic modelling (CSEM) where sources and receivers are located close to the seafloor separating conductive seawater and resistive formation. For the finite-difference scheme based on the Yee grid, the new interpolation is demonstrated to be much more accurate than alternative methods (interpolation using nodes on one side of the interface or interpolation using nodes on both sides, but ignoring the derivative jumps). The rigorous interpolation can handle arbitrary orientation of interface with respect to the grid, which is demonstrated on a marine CSEM example with a dipping seafloor. The interpolation coefficients are computed by minimizing a misfit between values at the nearest nodes and linear expansions of the continuous field components in the coordinate system aligned with the interface. The proposed interpolation operators can handle either uniform or non-uniform grids and can be applied to interpolation for both sources and receivers.
FDFD: A 3D Finite-Difference Frequency-Domain Code for Electromagnetic Induction Tomography
NASA Astrophysics Data System (ADS)
Champagne, Nathan J.; Berryman, James G.; Buettner, H. Michael
2001-07-01
A new 3D code for electromagnetic induction tomography with intended applications to environmental imaging problems has been developed. The approach consists of calculating the fields within a volume using an implicit finite-difference frequency-domain formulation. The volume is terminated by an anisotropic perfectly matched layer region that simulates an infinite domain by absorbing outgoing waves. Extensive validation of this code has been done using analytical and semianalytical results from other codes, and some of those results are presented in this paper. The new code is written in Fortran 90 and is designed to be easily parallelized. Finally, an adjoint field method of data inversion, developed in parallel for solving the fully nonlinear inverse problem for electrical conductivity imaging (e.g., for mapping underground conducting plumes), uses this code to provide solvers for both forward and adjoint fields. Results obtained from this inversion method for high-contrast media are encouraging and provide a significant improvement over those obtained from linearized inversion methods.
NASA Astrophysics Data System (ADS)
Bayona, Victor; Kindelan, Manuel
2013-10-01
Laminar flame propagation is an important problem in combustion modelling for which great advances have been achieved both in its theoretical understanding and in the numerical solution of the governing equations in 2D and 3D. Most of these numerical simulations use finite difference techniques on simple geometries (channels, ducts, ...) with equispaced nodes. The objective of this work is to explore the applicability of the radial basis function generated finite difference (RBF-FD) method to laminar flame propagation modelling. This method is specially well suited for the solution of problems with complex geometries and irregular boundaries. Another important advantage is that the method is independent of the dimension of the problem and, therefore, it is very easy to apply in 3D problems with complex geometries. In this work we use the RBF-FD method to compute 2D and 3D numerical results that simulate premixed laminar flames with different Lewis numbers propagating in open ducts.
NASA Astrophysics Data System (ADS)
Oprsal, I.; Faeh, D.; Giardini, D.
2002-12-01
The disastrous Basel earthquake of October 18, 1356 (I0=X, M ≈ 6.9), appeared in, today seismically modest, Basel region (Upper Rhine Graben). The lack of strong ground motion seismic data can be effectively supplied by numerical modeling. We applied the 3D finite differences (FD) to predict ground motions which can be used for microzonation and hazard assessment studies. The FD method is formulated for topography models on irregular rectangular grids. It is a 3D explicit FD formulation of the hyperbolic partial differential equation (PDE). Elastodynamic PDE is solved in the time domain. The Hooke's isotropic inhomogeneous medium contains discontinuities and a topographic free surface. The 3D elastic FD modeling is applied on a newly established P and S-wave velocities structure model. This complex structure contains main interfaces and gradients inside some layers. It is adjacent to the earth surface and includes topography (Kind, Faeh and Giardini, 2002, A 3D Reference Model for the Area of Basel, in prep.). The first attempt was done for a double-couple point source and relatively simple source function. Numerical tests are planned for several finite-extent source histories because the 1356 Basel earthquake source features have not been well determined, yet. The presumed finite-extent source is adjacent to the free surface. The results are compared to the macroseismic information of the Basel area.
Finite-difference model for 3-D flow in bays and estuaries
Smith, Peter E.; Larock, Bruce E.
1993-01-01
This paper describes a semi-implicit finite-difference model for the numerical solution of three-dimensional flow in bays and estuaries. The model treats the gravity wave and vertical diffusion terms in the governing equations implicitly, and other terms explicitly. The model achieves essentially second-order accurate and stable solutions in strongly nonlinear problems by using a three-time-level leapfrog-trapezoidal scheme for the time integration.
Finite-Difference Algorithm for Simulating 3D Electromagnetic Wavefields in Conductive Media
NASA Astrophysics Data System (ADS)
Aldridge, D. F.; Bartel, L. C.; Knox, H. A.
2013-12-01
Electromagnetic (EM) wavefields are routinely used in geophysical exploration for detection and characterization of subsurface geological formations of economic interest. Recorded EM signals depend strongly on the current conductivity of geologic media. Hence, they are particularly useful for inferring fluid content of saturated porous bodies. In order to enhance understanding of field-recorded data, we are developing a numerical algorithm for simulating three-dimensional (3D) EM wave propagation and diffusion in heterogeneous conductive materials. Maxwell's equations are combined with isotropic constitutive relations to obtain a set of six, coupled, first-order partial differential equations governing the electric and magnetic vectors. An advantage of this system is that it does not contain spatial derivatives of the three medium parameters electric permittivity, magnetic permeability, and current conductivity. Numerical solution methodology consists of explicit, time-domain finite-differencing on a 3D staggered rectangular grid. Temporal and spatial FD operators have order 2 and N, where N is user-selectable. We use an artificially-large electric permittivity to maximize the FD timestep, and thus reduce execution time. For the low frequencies typically used in geophysical exploration, accuracy is not unduly compromised. Grid boundary reflections are mitigated via convolutional perfectly matched layers (C-PMLs) imposed at the six grid flanks. A shared-memory-parallel code implementation via OpenMP directives enables rapid algorithm execution on a multi-thread computational platform. Good agreement is obtained in comparisons of numerically-generated data with reference solutions. EM wavefields are sourced via point current density and magnetic dipole vectors. Spatially-extended inductive sources (current carrying wire loops) are under development. We are particularly interested in accurate representation of high-conductivity sub-grid-scale features that are common
NASA Astrophysics Data System (ADS)
Anderson, T. S.; Miller, R.; Greenfield, R.; Fisk, D.
2002-12-01
The propagation of seismic waves through regions of complex topography is not thoroughly understood. Surface waves, are of particular interest, as they are large in amplitude and can characterize the source depth, magnitude, and frequency content. The amplitude and frequency content of seismic waves that propagate in regions with large topographical variations are affected by both the scattering and blockage of the wave energy. The ability to predict the 3-d scattering due to topography will improve the understanding of both regional scale surface wave magnitudes, and refine surface wave discriminants as well as at the local scale (<2 km ) where it will aid in the development of rule of thumb guide lines for array sensor placement for real time sensing technologies. Ideally, when validating the numerical accuracy of a propagation model against field data, the input geologic parameters would be known and thus eliminates geology as a source of error in the calculation. In March of 2001, Kansas Geological Survey (KGS) performed a detailed seismic site characterization at the Smart Weapons Test Range, Yuma Proving Ground, Arizona. The result of the KGS characterization study is a high-resolution 3-d model that is used in our seismic simulations. The velocities Vs, Vp are calculated by tomography and refraction, attenuation coefficients estimated from the surface wave and from p-waves and are provided in a model with attributes resolved in 3-d to 0.5 meters. In the present work, we present comparisons of synthetic data with seismic data collected at the Smart Weapons Test Range to benchmark the accuracy achieved in simulating 3-d wave propagation in the vicinity of a topographical anomaly (trench). Synthetic seismograms are generated using a 3-d 8th order staggered grid visco-elastic finite difference code that accounts for topography. The geologic model is based on the Yuma site characterization. The size of these calculations required use of the DoD High Performance
Accurate 3-D finite difference computation of traveltimes in strongly heterogeneous media
NASA Astrophysics Data System (ADS)
Noble, M.; Gesret, A.; Belayouni, N.
2014-12-01
Seismic traveltimes and their spatial derivatives are the basis of many imaging methods such as pre-stack depth migration and tomography. A common approach to compute these quantities is to solve the eikonal equation with a finite-difference scheme. If many recently published algorithms for resolving the eikonal equation do now yield fairly accurate traveltimes for most applications, the spatial derivatives of traveltimes remain very approximate. To address this accuracy issue, we develop a new hybrid eikonal solver that combines a spherical approximation when close to the source and a plane wave approximation when far away. This algorithm reproduces properly the spherical behaviour of wave fronts in the vicinity of the source. We implement a combination of 16 local operators that enables us to handle velocity models with sharp vertical and horizontal velocity contrasts. We associate to these local operators a global fast sweeping method to take into account all possible directions of wave propagation. Our formulation allows us to introduce a variable grid spacing in all three directions of space. We demonstrate the efficiency of this algorithm in terms of computational time and the gain in accuracy of the computed traveltimes and their derivatives on several numerical examples.
3D Finite-Difference Modeling of Acoustic Radiation from Seismic Sources
NASA Astrophysics Data System (ADS)
Chael, E. P.; Aldridge, D. F.; Jensen, R. P.
2013-12-01
Shallow seismic events, earthquakes as well as explosions, often generate acoustic waves in the atmosphere observable at local or even regional distances. Recording both the seismic and acoustic signals can provide additional constraints on source parameters such as epicenter coordinates, depth, origin time, moment, and mechanism. Recent advances in finite-difference (FD) modeling methods enable accurate numerical treatment of wave propagation across the ground surface between the (solid) elastic and (fluid) acoustic domains. Using a fourth-order, staggered-grid, velocity-stress FD algorithm, we are investigating the effects of various source parameters on the acoustic (or infrasound) signals transmitted from the solid earth into the atmosphere. Compressional (P), shear (S), and Rayleigh waves all radiate some acoustic energy into the air at the ground surface. These acoustic wavefronts are typically conical in shape, since their phase velocities along the surface exceed the sound speed in air. Another acoustic arrival with a spherical wavefront can be generated from the vicinity of the epicenter of a shallow event, due to the strong vertical ground motions directly above the buried source. Images of acoustic wavefields just above the surface reveal the radiation patterns and relative amplitudes of the various arrivals. In addition, we compare the relative effectiveness of different seismic source mechanisms for generating acoustic energy. For point sources at a fixed depth, double-couples with almost any orientation produce stronger acoustic signals than isotropic explosions, due to higher-amplitude S and Rayleigh waves. Of course, explosions tend to be shallower than most earthquakes, which can offset the differences due to mechanism. Low-velocity material in the shallow subsurface acts to increase vertical seismic motions there, enhancing the coupling to acoustic waves in air. If either type of source breaks the surface (e.g., an earthquake with surface rupture
NASA Astrophysics Data System (ADS)
Li, Y.; Han, B.; Métivier, L.; Brossier, R.
2016-09-01
We investigate an optimal fourth-order staggered-grid finite-difference scheme for 3D frequency-domain viscoelastic wave modeling. An anti-lumped mass strategy is incorporated to minimize the numerical dispersion. The optimal finite-difference coefficients and the mass weighting coefficients are obtained by minimizing the misfit between the normalized phase velocities and the unity. An iterative damped least-squares method, the Levenberg-Marquardt algorithm, is utilized for the optimization. Dispersion analysis shows that the optimal fourth-order scheme presents less grid dispersion and anisotropy than the conventional fourth-order scheme with respect to different Poisson's ratios. Moreover, only 3.7 grid-points per minimum shear wavelength are required to keep the error of the group velocities below 1%. The memory cost is then greatly reduced due to a coarser sampling. A parallel iterative method named CARP-CG is used to solve the large ill-conditioned linear system for the frequency-domain modeling. Validations are conducted with respect to both the analytic viscoacoustic and viscoelastic solutions. Compared with the conventional fourth-order scheme, the optimal scheme generates wavefields having smaller error under the same discretization setups. Profiles of the wavefields are presented to confirm better agreement between the optimal results and the analytic solutions.
Iterative methods for 3D implicit finite-difference migration using the complex Padé approximation
NASA Astrophysics Data System (ADS)
Costa, Carlos A. N.; Campos, Itamara S.; Costa, Jessé C.; Neto, Francisco A.; Schleicher, Jörg; Novais, Amélia
2013-08-01
Conventional implementations of 3D finite-difference (FD) migration use splitting techniques to accelerate performance and save computational cost. However, such techniques are plagued with numerical anisotropy that jeopardises the correct positioning of dipping reflectors in the directions not used for the operator splitting. We implement 3D downward continuation FD migration without splitting using a complex Padé approximation. In this way, the numerical anisotropy is eliminated at the expense of a computationally more intensive solution of a large-band linear system. We compare the performance of the iterative stabilized biconjugate gradient (BICGSTAB) and that of the multifrontal massively parallel direct solver (MUMPS). It turns out that the use of the complex Padé approximation not only stabilizes the solution, but also acts as an effective preconditioner for the BICGSTAB algorithm, reducing the number of iterations as compared to the implementation using the real Padé expansion. As a consequence, the iterative BICGSTAB method is more efficient than the direct MUMPS method when solving a single term in the Padé expansion. The results of both algorithms, here evaluated by computing the migration impulse response in the SEG/EAGE salt model, are of comparable quality.
Preliminary simulation of a M6.5 earthquake on the Seattle Fault using 3D finite-difference modeling
Stephenson, William J.; Frankel, Arthur D.
2000-01-01
A three-dimensional finite-difference simulation of a moderate-sized (M 6.5) thrust-faulting earthquake on the Seattle fault demonstrates the effects of the Seattle Basin on strong ground motion in the Puget lowland. The model area includes the cities of Seattle, Bremerton and Bellevue. We use a recently developed detailed 3D-velocity model of the Seattle Basin in these simulations. The model extended to 20-km depth and assumed rupture on a finite fault with random slip distribution. Preliminary results from simulations of frequencies 0.5 Hz and lower suggest amplification can occur at the surface of the Seattle Basin by the trapping of energy in the Quaternary sediments. Surface waves generated within the basin appear to contribute to amplification throughout the modeled region. Several factors apparently contribute to large ground motions in downtown Seattle: (1) radiation pattern and directivity from the rupture; (2) amplification and energy trapping within the Quaternary sediments; and (3) basin geometry and variation in depth of both Quaternary and Tertiary sediments
NASA Astrophysics Data System (ADS)
Schultz, A.
2010-12-01
describe our ongoing efforts to achieve massive parallelization on a novel hybrid GPU testbed machine currently configured with 12 Intel Westmere Xeon CPU cores (or 24 parallel computational threads) with 96 GB DDR3 system memory, 4 GPU subsystems which in aggregate contain 960 NVidia Tesla GPU cores with 16 GB dedicated DDR3 GPU memory, and a second interleved bank of 4 GPU subsystems containing in aggregate 1792 NVidia Fermi GPU cores with 12 GB dedicated DDR5 GPU memory. We are applying domain decomposition methods to a modified version of Weiss' (2001) 3D frequency domain full physics EM finite difference code, an open source GPL licensed f90 code available for download from www.OpenEM.org. This will be the core of a new hybrid 3D inversion that parallelizes frequencies across CPUs and individual forward solutions across GPUs. We describe progress made in modifying the code to use direct solvers in GPU cores dedicated to each small subdomain, iteratively improving the solution by matching adjacent subdomain boundary solutions, rather than iterative Krylov space sparse solvers as currently applied to the whole domain.
TACO3D. 3-D Finite Element Heat Transfer Code
Mason, W.E.
1992-03-04
TACO3D is a three-dimensional, finite-element program for heat transfer analysis. An extension of the two-dimensional TACO program, it can perform linear and nonlinear analyses and can be used to solve either transient or steady-state problems. The program accepts time-dependent or temperature-dependent material properties, and materials may be isotropic or orthotropic. A variety of time-dependent and temperature-dependent boundary conditions and loadings are available including temperature, flux, convection, and radiation boundary conditions and internal heat generation. Additional specialized features treat enclosure radiation, bulk nodes, and master/slave internal surface conditions (e.g., contact resistance). Data input via a free-field format is provided. A user subprogram feature allows for any type of functional representation of any independent variable. A profile (bandwidth) minimization option is available. The code is limited to implicit time integration for transient solutions. TACO3D has no general mesh generation capability. Rows of evenly-spaced nodes and rows of sequential elements may be generated, but the program relies on separate mesh generators for complex zoning. TACO3D does not have the ability to calculate view factors internally. Graphical representation of data in the form of time history and spatial plots is provided through links to the POSTACO and GRAPE postprocessor codes.
NASA Astrophysics Data System (ADS)
Liu, Yang; Sen, Mrinal K.
2011-09-01
Most conventional finite-difference methods adopt second-order temporal and (2M)th-order spatial finite-difference stencils to solve the 3D acoustic wave equation. When spatial finite-difference stencils devised from the time-space domain dispersion relation are used to replace these conventional spatial finite-difference stencils devised from the space domain dispersion relation, the accuracy of modelling can be increased from second-order along any directions to (2M)th-order along 48 directions. In addition, the conventional high-order spatial finite-difference modelling accuracy can be improved by using a truncated finite-difference scheme. In this paper, we combine the time-space domain dispersion-relation-based finite difference scheme and the truncated finite-difference scheme to obtain optimised spatial finite-difference coefficients and thus to significantly improve the modelling accuracy without increasing computational cost, compared with the conventional space domain dispersion-relation-based finite difference scheme. We developed absorbing boundary conditions for the 3D acoustic wave equation, based on predicting wavefield values in a transition area by weighing wavefield values from wave equations and one-way wave equations. Dispersion analyses demonstrate that high-order spatial finite-difference stencils have greater accuracy than low-order spatial finite-difference stencils for high frequency components of wavefields, and spatial finite-difference stencils devised in the time-space domain have greater precision than those devised in the space domain under the same discretisation. The modelling accuracy can be improved further by using the truncated spatial finite-difference stencils. Stability analyses show that spatial finite-difference stencils devised in the time-space domain have better stability condition. Numerical modelling experiments for homogeneous, horizontally layered and Society of Exploration Geophysicists/European Association of
3-D Finite Element Heat Transfer
Energy Science and Technology Software Center (ESTSC)
1992-02-01
TOPAZ3D is a three-dimensional implicit finite element computer code for heat transfer analysis. TOPAZ3D can be used to solve for the steady-state or transient temperature field on three-dimensional geometries. Material properties may be temperature-dependent and either isotropic or orthotropic. A variety of time-dependent and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation. By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functionalmore » representation of boundary conditions and internal heat generation. TOPAZ3D can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.« less
NASA Astrophysics Data System (ADS)
Borisov, Dmitry; Singh, Satish C.; Fuji, Nobuaki
2015-09-01
Seismic full waveform inversion is an objective method to estimate elastic properties of the subsurface and is an important area of research, particularly in seismic exploration community. It is a data-fitting approach, where the difference between observed and synthetic data is minimized iteratively. Due to a very high computational cost, the practical implementation of waveform inversion has so far been restricted to a 2-D geometry with different levels of physics incorporated in it (e.g. elasticity/viscoelasticity) or to a 3-D geometry but using an acoustic approximation. However, the earth is three-dimensional, elastic and heterogeneous and therefore a full 3-D elastic inversion is required in order to obtain more accurate and valuable models of the subsurface. Despite the recent increase in computing power, the application of 3-D elastic full waveform inversion to real-scale problems remains quite challenging on the current computer architecture. Here, we present an efficient method to perform 3-D elastic full waveform inversion for time-lapse seismic data using a finite-difference injection method. In this method, the wavefield is computed in the whole model and is stored on a surface above a finite volume where the model is perturbed and localized inversion is performed. Comparison of the final results using the 3-D finite-difference injection method and conventional 3-D inversion performed within the whole volume shows that our new method provides significant reductions in computational time and memory requirements without any notable loss in accuracy. Our approach shows a big potential for efficient reservoir monitoring in real time-lapse experiments.
NASA Technical Reports Server (NTRS)
Sohn, Kiho D.; Ip, Shek-Se P.
1988-01-01
Three-dimensional finite element models were generated and transferred into three-dimensional finite difference models to perform transient thermal analyses for the SSME high pressure fuel turbopump's first stage nozzles and rotor blades. STANCOOL was chosen to calculate the heat transfer characteristics (HTCs) around the airfoils, and endwall effects were included at the intersections of the airfoils and platforms for the steady-state boundary conditions. Free and forced convection due to rotation effects were also considered in hollow cores. Transient HTCs were calculated by taking ratios of the steady-state values based on the flow rates and fluid properties calculated at each time slice. Results are presented for both transient plots and three-dimensional color contour isotherm plots; they were also converted into universal files to be used for FEM stress analyses.
NASA Astrophysics Data System (ADS)
Son, Sang-Kil; Chu, Shih-I.
2008-05-01
We introduce a new computational method on unstructured grids in the three-dimensional (3D) spaces to investigate the electronic structure of polyatomic molecules. The Voronoi-cell finite difference (VFD) method realizes a simple discrete Laplacian operator on unstructured grids based on Voronoi cells and their natural neighbors. The feature of unstructured grids enables us to choose intuitive pictures for an optimal molecular grid system. The new VFD method achieves highly adaptability by the Voronoi-cell diagram and yet simplicity by the finite difference scheme. It has no limitation in local refinement of grids in the vicinity of nuclear positions and provides an explicit expression at each grid without any integration. This method augmented by unstructured molecular grids is suitable for solving the Schr"odinger equation with the realistic 3D Coulomb potentials regardless of symmetry of molecules. For numerical examples, we test accuracies for electronic structures of one-electron polyatomic systems: linear H2^+ and triangular H3^++. We also extend VFD to the density functional theory (DFT) for many-electron polyatomic molecules.
3D Finite Element Analysis of Particle-Reinforced Aluminum
NASA Technical Reports Server (NTRS)
Shen, H.; Lissenden, C. J.
2002-01-01
Deformation in particle-reinforced aluminum has been simulated using three distinct types of finite element model: a three-dimensional repeating unit cell, a three-dimensional multi-particle model, and two-dimensional multi-particle models. The repeating unit cell model represents a fictitious periodic cubic array of particles. The 3D multi-particle (3D-MP) model represents randomly placed and oriented particles. The 2D generalized plane strain multi-particle models were obtained from planar sections through the 3D-MP model. These models were used to study the tensile macroscopic stress-strain response and the associated stress and strain distributions in an elastoplastic matrix. The results indicate that the 2D model having a particle area fraction equal to the particle representative volume fraction of the 3D models predicted the same macroscopic stress-strain response as the 3D models. However, there are fluctuations in the particle area fraction in a representative volume element. As expected, predictions from 2D models having different particle area fractions do not agree with predictions from 3D models. More importantly, it was found that the microscopic stress and strain distributions from the 2D models do not agree with those from the 3D-MP model. Specifically, the plastic strain distribution predicted by the 2D model is banded along lines inclined at 45 deg from the loading axis while the 3D model prediction is not. Additionally, the triaxial stress and maximum principal stress distributions predicted by 2D and 3D models do not agree. Thus, it appears necessary to use a multi-particle 3D model to accurately predict material responses that depend on local effects, such as strain-to-failure, fracture toughness, and fatigue life.
NASA Astrophysics Data System (ADS)
Lee, H.; Min, D.; Lim, S.; Yang, J.; Kwon, B.; Yoo, H.
2009-12-01
In a conventional marine seismic data analysis, pressure data have been usually interpreted on the basis of acoustic wave equation. The acoustic wave equation, however, only deals with P-wave propagation, and it cannot correctly describe the wave propagation in acoustic-elastic (fluid-solid) coupled media. Recently, in 4C OBC survey (4-component ocean bottom cable), it is possible to acquire both pressure and 3-component displacements (measured at the sea-bottom). Combining pressure and displacement data allows us to interpret subsurface structures more accurately. In order to accurately simulate wave propagation in fluid-solid coupled media, we need an acoustic-elastic coupled modeling algorithm, which deals with displacements in elastic region and pressure in acoustic region. For waveform inversion and reverse-time migration that require a great number of forward modeling, it is essential to develop an efficient scheme that reduces computing time and computer core memory. In this study, we present a 3D time-domain acoustic-elastic coupled modeling algorithm on the basis of the cell-based finite difference method. The cell-based method has proven to properly describe the free-surface boundary, which indicates that it will also properly describe the fluid-solid interface boundaries. In the acoustic-elastic coupled modeling, we first compose cell-based finite differences individually for the 3D acoustic and elastic media, and then combine the differences using the fluid-solid interface boundary conditions. Considering that the 2D acoustic-elastic coupled modeling algorithm gives numerical solutions comparable to analytic solutions, we expect that the 3D acoustic-elastic coupled modeling will correctly describe wave propagation in the fluid-solid coupled media. We apply our algorithm to 3D horizontal two- and three-layer models. Numerical experiments show that the cell-based coupled modeling algorithm properly describes S- and converted waves as well as P-waves. The
3-D Finite Element Code Postprocessor
Energy Science and Technology Software Center (ESTSC)
1996-07-15
TAURUS is an interactive post-processing application supporting visualization of finite element analysis results on unstructured grids. TAURUS provides the ability to display deformed geometries and contours or fringes of a large number of derived results on meshes consisting of beam, plate, shell, and solid type finite elements. Time history plotting is also available.
NASA Astrophysics Data System (ADS)
Wright, G.; Flyer, N.; Yuen, D. A.; Monnereau, M.; Zhang, S.; Wang, S. M.
2009-05-01
Many numerical methods, such as finite-differences, finite-volume, their yin-yang variants, finite-elements and spectral methods have been employed to study 3-D mantle convection. All have their own strengths, but also serious weaknesses. Spectrally accurate methods do not practically allow for node refinement and often involve cumbersome algebra while finite difference, volume, or element methods are generally low-order, adding excessive numerical diffusion to the model. For the 3-D mantle convection problem, we have introduced a new mesh-free approach, using radial basis functions (RBF). This method has the advantage of being algorithmic simple, spectrally accurate for arbitrary node layouts in multi-dimensions and naturally allows for node-refinement. One virtue of the RBF scheme allows the user to use a simple Cartesian geometry, while implementing the required boundary conditions for the temperature, velocities and stress components on a spherical surface at both the planetary surface and the core-mantle boundary. We have studied time- dependent mantle convection, using both a RBF-pseudospectral code and a code which uses spherical- harmonics in the angular direction and second-order finite volume in the radial direction. We have employed a third code , which uses spherical harmonics and higher-order finite-difference method a la Fornberg in the radial coordinate.We first focus on the onset of time-dependence at Rayleigh number Ra of 70,000. We follow the development of stronger time-dependence to a Ra of one million, using high enough resolution with 120 to 200 points in the radial direction and 128 to 256 spherical harmonics.
Yang, Chun; Tang, Dalin; Atluri, Satya
2011-01-01
Previously, we introduced a computational procedure based on three-dimensional meshless generalized finite difference (MGFD) method and serial magnetic resonance imaging (MRI) data to quantify patient-specific carotid atherosclerotic plaque growth functions and simulate plaque progression. Structure-only models were used in our previous report. In this paper, fluid-stricture interaction (FSI) was added to improve on prediction accuracy. One participating patient was scanned three times (T1, T2, and T3, at intervals of about 18 months) to obtain plaque progression data. Blood flow was assumed to laminar, Newtonian, viscous and incompressible. The Navier-Stokes equations with arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing equations. Plaque material was assumed to be uniform, homogeneous, isotropic, linear, and nearly incompressible. The linear elastic model was used. The 3D FSI plaque model was discretized and solved using a meshless generalized finite difference (GFD) method. Growth functions with a) morphology alone; b) morphology and plaque wall stress (PWS); morphology and flow shear stress (FSS), and d) morphology, PWS and FSS were introduced to predict future plaque growth based on previous time point data. Starting from the T2 plaque geometry, plaque progression was simulated by solving the FSI model and adjusting plaque geometry using plaque growth functions iteratively until T3 is reached. Numerically simulated plaque progression agreed very well with the target T3 plaque geometry with errors ranging from 8.62%, 7.22%, 5.77% and 4.39%, with the growth function including morphology, plaque wall stress and flow shear stress terms giving the best predictions. Adding flow shear stress term to the growth function improved the prediction error from 7.22% to 4.39%, a 40% improvement. We believe this is the first time 3D plaque progression FSI simulation based on multi-year patient-tracking data was reported. Serial MRI-based progression
3D finite element simulations of high velocity projectile impact
NASA Astrophysics Data System (ADS)
Ožbolt, Joško; İrhan, Barış; Ruta, Daniela
2015-09-01
An explicit three-dimensional (3D) finite element (FE) code is developed for the simulation of high velocity impact and fragmentation events. The rate sensitive microplane material model, which accounts for large deformations and rate effects, is used as a constitutive law. In the code large deformation frictional contact is treated by forward incremental Lagrange multiplier method. To handle highly distorted and damaged elements the approach based on the element deletion is employed. The code is then used in 3D FE simulations of high velocity projectile impact. The results of the numerical simulations are evaluated and compared with experimental results. It is shown that it realistically predicts failure mode and exit velocities for different geometries of plain concrete slab. Moreover, the importance of some relevant parameters, such as contact friction, rate sensitivity, bulk viscosity and deletion criteria are addressed.
Finite sampling corrected 3D noise with confidence intervals.
Haefner, David P; Burks, Stephen D
2015-05-20
When evaluated with a spatially uniform irradiance, an imaging sensor exhibits both spatial and temporal variations, which can be described as a three-dimensional (3D) random process considered as noise. In the 1990s, NVESD engineers developed an approximation to the 3D power spectral density for noise in imaging systems known as 3D noise. The goal was to decompose the 3D noise process into spatial and temporal components identify potential sources of origin. To characterize a sensor in terms of its 3D noise values, a finite number of samples in each of the three dimensions (two spatial, one temporal) were performed. In this correspondence, we developed the full sampling corrected 3D noise measurement and the corresponding confidence bounds. The accuracy of these methods was demonstrated through Monte Carlo simulations. Both the sampling correction as well as the confidence intervals can be applied a posteriori to the classic 3D noise calculation. The Matlab functions associated with this work can be found on the Mathworks file exchange ["Finite sampling corrected 3D noise with confidence intervals," https://www.mathworks.com/matlabcentral/fileexchange/49657-finite-sampling-corrected-3d-noise-with-confidence-intervals.]. PMID:26192530
Ozan, Oguz; Ramoglu, Serhat
2015-06-01
Implant-supported overdentures with self-aligning attachment systems are preferred to improve the stability and retention of complete dentures. The positioning of the implant attachments is a very important aspect of two-implant overdentures in obtaining better stress distribution. Therefore, the objective of this study was to compare two different attachment systems in a two-implant overdenture by evaluating the stress distributions in peri-implant bone and stresses on the attachments with positioning at different height levels using the 3D FEA method. Six models with ball attachments and 6 models with locator attachments-totaling 12 models (including 2 controls)-with the left implant positioned unilaterally at different height levels were subjected to 3 loading conditions (anterior, right posterior, and left posterior). Data for Von Misses stresses were produced numerically, color coded, and compared among the models for attachments and peri-implant cortical bone. The configurations in which implants presented 3 mm height differences in the bone level showed the most successful results in the peri-implant bone. When stresses on the attachments were compared, greater stress values were obtained from the ball attachments. As a conclusion, the configurations with a considerable (3 mm) height difference between quadrants of the mandible in the anterior segment showed the most successful results in the peri-implant bone. On the contrary, peak stress values around the implant observed from the models with less (1 mm) bone height difference may require leveling of the bone during surgery. However, these findings should be corroborated with clinical studies. PMID:24471769
Finite volume and finite element methods applied to 3D laminar and turbulent channel flows
Louda, Petr; Příhoda, Jaromír; Sváček, Petr; Kozel, Karel
2014-12-10
The work deals with numerical simulations of incompressible flow in channels with rectangular cross section. The rectangular cross section itself leads to development of various secondary flow patterns, where accuracy of simulation is influenced by numerical viscosity of the scheme and by turbulence modeling. In this work some developments of stabilized finite element method are presented. Its results are compared with those of an implicit finite volume method also described, in laminar and turbulent flows. It is shown that numerical viscosity can cause errors of same magnitude as different turbulence models. The finite volume method is also applied to 3D turbulent flow around backward facing step and good agreement with 3D experimental results is obtained.
Advances in 3D electromagnetic finite element modeling
Nelson, E.M.
1997-08-01
Numerous advances in electromagnetic finite element analysis (FEA) have been made in recent years. The maturity of frequency domain and eigenmode calculations, and the growth of time domain applications is briefly reviewed. A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will also be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis is also discussed.
3D unstructured mesh discontinuous finite element hydro
Prasad, M.K.; Kershaw, D.S.; Shaw, M.J.
1995-07-01
The authors present detailed features of the ICF3D hydrodynamics code used for inertial fusion simulations. This code is intended to be a state-of-the-art upgrade of the well-known fluid code, LASNEX. ICF3D employs discontinuous finite elements on a discrete unstructured mesh consisting of a variety of 3D polyhedra including tetrahedra, prisms, and hexahedra. The authors discussed details of how the ROE-averaged second-order convection was applied on the discrete elements, and how the C++ coding interface has helped to simplify implementing the many physics and numerics modules within the code package. The author emphasized the virtues of object-oriented design in large scale projects such as ICF3D.
3-D dynamic rupture simulations by a finite volume method
NASA Astrophysics Data System (ADS)
Benjemaa, M.; Glinsky-Olivier, N.; Cruz-Atienza, V. M.; Virieux, J.
2009-07-01
Dynamic rupture of a 3-D spontaneous crack of arbitrary shape is investigated using a finite volume (FV) approach. The full domain is decomposed in tetrahedra whereas the surface, on which the rupture takes place, is discretized with triangles that are faces of tetrahedra. First of all, the elastodynamic equations are described into a pseudo-conservative form for an easy application of the FV discretization. Explicit boundary conditions are given using criteria based on the conservation of discrete energy through the crack surface. Using a stress-threshold criterion, these conditions specify fluxes through those triangles that have suffered rupture. On these broken surfaces, stress follows a linear slip-weakening law, although other friction laws can be implemented. For The Problem Version 3 of the dynamic-rupture code verification exercise conducted by the SCEC/USGS, numerical solutions on a planar fault exhibit a very high convergence rate and are in good agreement with the reference one provided by a finite difference (FD) technique. For a non-planar fault of parabolic shape, numerical solutions agree satisfactorily well with those obtained with a semi-analytical boundary integral method in terms of shear stress amplitudes, stopping phases arrival times and stress overshoots. Differences between solutions are attributed to the low-order interpolation of the FV approach, whose results are particularly sensitive to the mesh regularity (structured/unstructured). We expect this method, which is well adapted for multiprocessor parallel computing, to be competitive with others for solving large scale dynamic ruptures scenarios of seismic sources in the near future.
Not Available
1984-10-01
STEALTH is a family of computer codes that can be used to calculate a variety of physical processes in which the dynamic behavior of a continuum is involved. The version of STEALTH described in this volume is designed for calculations of fluid-structure interaction. This version of the program consists of a hydrodynamic version of STEALTH which has been coupled to a finite-element code, WHAMSE. STEALTH computes the transient response of the fluid continuum, while WHAMSE computes the transient response of shell and beam structures under external fluid loadings. The coupling between STEALTH and WHAMSE is performed during each cycle or step of a calculation. Separate calculations of fluid response and structure response are avoided, thereby giving a more accurate model of the dynamic coupling between fluid and structure. This volume provides the theoretical background, the finite-difference equations, the finite-element equations, a discussion of several sample problems, a listing of the input decks for the sample problems, a programmer's manual and a description of the input records for the STEALTH/WHAMSE computer program.
3D Finite Element Trajectory Code with Adaptive Meshing
NASA Astrophysics Data System (ADS)
Ives, Lawrence; Bui, Thuc; Vogler, William; Bauer, Andy; Shephard, Mark; Beal, Mark; Tran, Hien
2004-11-01
Beam Optics Analysis, a new, 3D charged particle program is available and in use for the design of complex, 3D electron guns and charged particle devices. The code reads files directly from most CAD and solid modeling programs, includes an intuitive Graphical User Interface (GUI), and a robust mesh generator that is fully automatic. Complex problems can be set up, and analysis initiated in minutes. The program includes a user-friendly post processor for displaying field and trajectory data using 3D plots and images. The electrostatic solver is based on the standard nodal finite element method. The magnetostatic field solver is based on the vector finite element method and is also called during the trajectory simulation process to solve for self magnetic fields. The user imports the geometry from essentially any commercial CAD program and uses the GUI to assign parameters (voltages, currents, dielectric constant) and designate emitters (including work function, emitter temperature, and number of trajectories). The the mesh is generated automatically and analysis is performed, including mesh adaptation to improve accuracy and optimize computational resources. This presentation will provide information on the basic structure of the code, its operation, and it's capabilities.
Zherdetsky, Aleksej; Prakonina, Alena; Malony, Allen D.
2014-01-01
The Electrical Impedance Tomography (EIT) and electroencephalography (EEG) forward problems in anisotropic inhomogeneous media like the human head belongs to the class of the three-dimensional boundary value problems for elliptic equations with mixed derivatives. We introduce and explore the performance of several new promising numerical techniques, which seem to be more suitable for solving these problems. The proposed numerical schemes combine the fictitious domain approach together with the finite-difference method and the optimally preconditioned Conjugate Gradient- (CG-) type iterative method for treatment of the discrete model. The numerical scheme includes the standard operations of summation and multiplication of sparse matrices and vector, as well as FFT, making it easy to implement and eligible for the effective parallel implementation. Some typical use cases for the EIT/EEG problems are considered demonstrating high efficiency of the proposed numerical technique. PMID:24527060
Turovets, Sergei; Volkov, Vasily; Zherdetsky, Aleksej; Prakonina, Alena; Malony, Allen D
2014-01-01
The Electrical Impedance Tomography (EIT) and electroencephalography (EEG) forward problems in anisotropic inhomogeneous media like the human head belongs to the class of the three-dimensional boundary value problems for elliptic equations with mixed derivatives. We introduce and explore the performance of several new promising numerical techniques, which seem to be more suitable for solving these problems. The proposed numerical schemes combine the fictitious domain approach together with the finite-difference method and the optimally preconditioned Conjugate Gradient- (CG-) type iterative method for treatment of the discrete model. The numerical scheme includes the standard operations of summation and multiplication of sparse matrices and vector, as well as FFT, making it easy to implement and eligible for the effective parallel implementation. Some typical use cases for the EIT/EEG problems are considered demonstrating high efficiency of the proposed numerical technique. PMID:24527060
A finite element solver for 3-D compressible viscous flows
NASA Technical Reports Server (NTRS)
Reddy, K. C.; Reddy, J. N.; Nayani, S.
1990-01-01
Computation of the flow field inside a space shuttle main engine (SSME) requires the application of state of the art computational fluid dynamic (CFD) technology. Several computer codes are under development to solve 3-D flow through the hot gas manifold. Some algorithms were designed to solve the unsteady compressible Navier-Stokes equations, either by implicit or explicit factorization methods, using several hundred or thousands of time steps to reach a steady state solution. A new iterative algorithm is being developed for the solution of the implicit finite element equations without assembling global matrices. It is an efficient iteration scheme based on a modified nonlinear Gauss-Seidel iteration with symmetric sweeps. The algorithm is analyzed for a model equation and is shown to be unconditionally stable. Results from a series of test problems are presented. The finite element code was tested for couette flow, which is flow under a pressure gradient between two parallel plates in relative motion. Another problem that was solved is viscous laminar flow over a flat plate. The general 3-D finite element code was used to compute the flow in an axisymmetric turnaround duct at low Mach numbers.
Finite element solver for 3-D compressible viscous flows
NASA Technical Reports Server (NTRS)
Reddy, K. C.; Reddy, J. N.
1986-01-01
The space shuttle main engine (SSME) has extremely complex internal flow structure. The geometry of the flow domain is three-dimensional with complicated topology. The flow is compressible, viscous, and turbulent with large gradients in flow quantities and regions of recirculations. The analysis of the flow field in SSME involves several tedious steps. One is the geometrical modeling of the particular zone of the SSME being studied. Accessing the geometry definition, digitalizing it, and developing surface interpolations suitable for an interior grid generator require considerable amount of manual labor. There are several types of grid generators available with some general-purpose finite element programs. An efficient and robust computational scheme for solving 3D Navier-Stokes equations has to be implemented. Post processing software has to be adapted to visualize and analyze the computed 3D flow field. The progress made in a project to develop software for the analysis of the flow is discussed. The technical approach to the development of the finite element scheme and the relaxation procedure are discussed. The three dimensional finite element code for the compressible Navier-Stokes equations is listed.
Higher Order Lagrange Finite Elements In M3D
J. Chen; H.R. Strauss; S.C. Jardin; W. Park; L.E. Sugiyama; G. Fu; J. Breslau
2004-12-17
The M3D code has been using linear finite elements to represent multilevel MHD on 2-D poloidal planes. Triangular higher order elements, up to third order, are constructed here in order to provide M3D the capability to solve highly anisotropic transport problems. It is found that higher order elements are essential to resolve the thin transition layer characteristic of the anisotropic transport equation, particularly when the strong anisotropic direction is not aligned with one of the Cartesian coordinates. The transition layer is measured by the profile width, which is zero for infinite anisotropy. It is shown that only higher order schemes have the ability to make this layer converge towards zero when the anisotropy gets stronger and stronger. Two cases are considered. One has the strong transport direction partially aligned with one of the element edges, the other doesn't have any alignment. Both cases have the strong transport direction misaligned with the grid line by some angles.
NASA Astrophysics Data System (ADS)
Li, Dian-Sen; Fang, Dai-Ning; Lu, Zi-Xing; Yang, Zhen-Yu; Jiang, Nan
2010-08-01
In the first part of the work, we have established a new parameterized three-dimensional (3D) finite element model (FEM) which precisely simulated the spatial configuration of the braiding yarns and considered the cross-section deformation as well as the surface contact relationship between the yarns. This paper presents a prediction of the effective elastic properties and the meso-scale mechanical response of 3D braided composites to verify the validation of the FEM. The effects of the braiding parameters on the mechanical properties are investigated in detail. By analyzing the deformation and stress nephogram of the model, a reasonable overall stress field is provided and the results well support the strength prediction. The results indicate it is convenient to predict all the elastic constants of 3D braided composites with different parameters simultaneously using the FEM. Moreover, the FEM can successfully predict the meso-scale mechanical response of 3D braided composites containing periodical structures.
Beam and Truss Finite Element Verification for DYNA3D
Rathbun, H J
2007-07-16
The explicit finite element (FE) software program DYNA3D has been developed at Lawrence Livermore National Laboratory (LLNL) to simulate the dynamic behavior of structures, systems, and components. This report focuses on verification of beam and truss element formulations in DYNA3D. An efficient protocol has been developed to verify the accuracy of these structural elements by generating a set of representative problems for which closed-form quasi-static steady-state analytical reference solutions exist. To provide as complete coverage as practically achievable, problem sets are developed for each beam and truss element formulation (and their variants) in all modes of loading and physical orientation. Analyses with loading in the elastic and elastic-plastic regimes are performed. For elastic loading, the FE results are within 1% of the reference solutions for all cases. For beam element bending and torsion loading in the plastic regime, the response is heavily dependent on the numerical integration rule chosen, with higher refinement yielding greater accuracy (agreement to within 1%). Axial loading in the plastic regime produces accurate results (agreement to within 0.01%) for all integration rules and element formulations. Truss elements are also verified to provide accurate results (within 0.01%) for elastic and elastic-plastic loading. A sample problem to verify beam element response in ParaDyn, the parallel version DYNA3D, is also presented.
3-D Finite Element Analyses of the Egan Cavern Field
Klamerus, E.W.; Ehgartner, B.L.
1999-02-01
Three-dimensional finite element analyses were performed for the two gas-filled storage caverns at the Egan field, Jennings dome, Louisiana. The effects of cavern enlargement on surface subsidence, storage loss, and cavern stability were investigated. The finite element model simulated the leaching of caverns to 6 and 8 billion cubic feet (BCF) and examined their performance at various operating conditions. Operating pressures varied from 0.15 psi/ft to 0.9 psi/ft at the bottom of the lowest cemented casing. The analysis also examined the stability of the web or pillar of salt between the caverns under differential pressure loadings. The 50-year simulations were performed using JAC3D, a three dimensional finite element analysis code for nonlinear quasistatic solids. A damage criterion based on onset of dilatancy was used to evaluate cavern instability. Dilation results from the development of microfractures in salt and, hence, potential increases in permeability onset occurs well before large scale failure. The analyses predicted stable caverns throughout the 50-year period for the range of pressures investigated. Some localized salt damage was predicted near the bottom walls of the caverns if the caverns are operated at minimum pressure for long periods of time. Volumetric cavern closures over time due to creep were moderate to excessive depending on the salt creep properties and operating pressures. However, subsidence above the cavern field was small and should pose no problem, to surface facilities.
NASA Astrophysics Data System (ADS)
Skarlatoudis, A. A.; Papazachos, C. B.; Theodoulidis, N.; Kristek, J.; Moczo, P.
2010-07-01
The site effects of seismic motion in the metropolitan area of the city of Thessaloniki (Northern Greece) are investigated using a 3-D finite-difference modelling approach. Three different seismic scenarios are assumed with two different focal mechanisms for each one. Standard spectral ratios (SSR) are calculated from 3-D synthetics and compared with the ratios from the recorded motion, as well as ratios obtained from 1-D and 2-D modelling by other researchers. The average SSR curves from the six scenarios are in good agreement with the empirical ones, whereas the SSR results from 3-D modelling are different from those from 1-D modelling, exhibiting higher fundamental frequencies and larger amplification amplitudes, in much better agreement with observed SSR ratios. Comparisons of Fourier amplitude spectra obtained for various scenarios for the broader area of Thessaloniki show considerable dependence of the site effects on the source properties (position, depth and fault-plane solution), as well as on the local structure.
3D finite element model for treatment of cleft lip
NASA Astrophysics Data System (ADS)
Jiao, Chun; Hong, Dongming; Lu, Hongbing; Wang, Jianqi; Lin, Qin; Liang, Zhengrong
2009-02-01
Cleft lip is a congenital facial deformity with high occurrence rate in China. Surgical procedure involving Millard or Tennison methods is usually employed for treatment of cleft lip. However, due to the elasticity of the soft tissues and the mechanical interaction between skin and maxillary, the occurrence rate of facial abnormality or dehisce is still high after the surgery, leading to multiple operations of the patient. In this study, a framework of constructing a realistic 3D finite element model (FEM) for the treatment of cleft lip has been established. It consists of two major steps. The first one is the reconstruction of a 3D geometrical model of the cleft lip from scanning CT data. The second step is the build-up of a FEM for cleft lip using the geometric model, where the material property of all the tetrahedrons was calculated from the CT densities directly using an empirical curve. The simulation results demonstrated (1) the deformation procedure of the model step-by-step when forces were applied, (2) the stress distribution inside the model, and (3) the displacement of all elements in the model. With the computer simulation, the minimal force of having the cleft be repaired is predicted, as well as whether a given force sufficient for the treatment of a specific individual. It indicates that the proposed framework could integrate the treatment planning with stress analysis based on a realistic patient model.
3D finite element simulation of TIG weld pool
NASA Astrophysics Data System (ADS)
Kong, X.; Asserin, O.; Gounand, S.; Gilles, P.; Bergheau, J. M.; Medale, M.
2012-07-01
The aim of this paper is to propose a three-dimensional weld pool model for the moving gas tungsten arc welding (GTAW) process, in order to understand the main factors that limit the weld quality and improve the productivity, especially with respect to the welding speed. Simulation is a very powerful tool to help in understanding the physical phenomena in the weld process. A 3D finite element model of heat and fluid flow in weld pool considering free surface of the pool and traveling speed has been developed for the GTAW process. Cast3M software is used to compute all the governing equations. The free surface of the weld pool is calculated by minimizing the total surface energy. The combined effects of surface tension gradient, buoyancy force, arc pressure, arc drag force to drive the fluid flow is included in our model. The deformation of the weld pool surface and the welding speed affect fluid flow, heat flow and thus temperature gradients and molten pool dimensions. Welding trials study is presented to compare our numerical results with macrograph of the molten pool.
A finite element analysis of a 3D auxetic textile structure for composite reinforcement
NASA Astrophysics Data System (ADS)
Ge, Zhaoyang; Hu, Hong; Liu, Yanping
2013-08-01
This paper reports the finite element analysis of an innovative 3D auxetic textile structure consisting of three yarn systems (weft, warp and stitch yarns). Different from conventional 3D textile structures, the proposed structure exhibits an auxetic behaviour under compression and can be used as a reinforcement to manufacture auxetic composites. The geometry of the structure is first described. Then a 3D finite element model is established using ANSYS software and validated by the experimental results. The deformation process of the structure at different compression strains is demonstrated, and the validated finite element model is finally used to simulate the auxetic behaviour of the structure with different structural parameters and yarn properties. The results show that the auxetic behaviour of the proposed structure increases with increasing compression strain, and all the structural parameters and yarn properties have significant effects on the auxetic behaviour of the structure. It is expected that the study could provide a better understanding of 3D auxetic textile structures and could promote their application in auxetic composites.
NASA Technical Reports Server (NTRS)
Zhao, W.; Newman, J. C., Jr.; Sutton, M. A.; Wu, X. R.; Shivakumar, K. N.
1995-01-01
Stress intensity factors for quarter-elliptical corner cracks emanating from a circular hole are determined using a 3-D weight function method combined with a 3-D finite element method. The 3-D finite element method is used to analyze uncracked configuration and provide stress distribution in the region where crack is to occur. Using this stress distribution as input, the 3-D weight function method is used to determine stress intensity factors. Three different loading conditions, i.e. remote tension, remote bending and wedge loading, are considered for a wide range in geometrical parameters. The significance in using 3-D uncracked stress distribution and the difference between single and double corner cracks are studied. Typical crack opening displacements are also provided. Comparisons are made with solutions available in the literature.
NASA Astrophysics Data System (ADS)
Popov, Anton; Kaus, Boris
2015-04-01
This software project aims at bringing the 3D lithospheric deformation modeling to a qualitatively different level. Our code LaMEM (Lithosphere and Mantle Evolution Model) is based on the following building blocks: * Massively-parallel data-distributed implementation model based on PETSc library * Light, stable and accurate staggered-grid finite difference spatial discretization * Marker-in-Cell pedictor-corector time discretization with Runge-Kutta 4-th order * Elastic stress rotation algorithm based on the time integration of the vorticity pseudo-vector * Staircase-type internal free surface boundary condition without artificial viscosity contrast * Geodynamically relevant visco-elasto-plastic rheology * Global velocity-pressure-temperature Newton-Raphson nonlinear solver * Local nonlinear solver based on FZERO algorithm * Coupled velocity-pressure geometric multigrid preconditioner with Galerkin coarsening Staggered grid finite difference, being inherently Eulerian and rather complicated discretization method, provides no natural treatment of free surface boundary condition. The solution based on the quasi-viscous sticky-air phase introduces significant viscosity contrasts and spoils the convergence of the iterative solvers. In LaMEM we are currently implementing an approximate stair-case type of the free surface boundary condition which excludes the empty cells and restores the solver convergence. Because of the mutual dependence of the stress and strain-rate tensor components, and their different spatial locations in the grid, there is no straightforward way of implementing the nonlinear rheology. In LaMEM we have developed and implemented an efficient interpolation scheme for the second invariant of the strain-rate tensor, that solves this problem. Scalable efficient linear solvers are the key components of the successful nonlinear problem solution. In LaMEM we have a range of PETSc-based preconditioning techniques that either employ a block factorization of
Vector algorithms for geometrically nonlinear 3D finite element analysis
NASA Technical Reports Server (NTRS)
Whitcomb, John D.
1989-01-01
Algorithms for geometrically nonlinear finite element analysis are presented which exploit the vector processing capability of the VPS-32, which is closely related to the CYBER 205. By manipulating vectors (which are long lists of numbers) rather than individual numbers, very high processing speeds are obtained. Long vector lengths are obtained without extensive replication or reordering by storage of intermediate results in strategic patterns at all stages of the computations. Comparisons of execution times with those from programs using either scalar or other vector programming techniques indicate that the algorithms presented are quite efficient.
Equivalent Body Force Finite Elements Method and 3-D Earth Model Applied In 2004 Sumatra Earthquake
NASA Astrophysics Data System (ADS)
Qu, W.; Cheng, H.; Shi, Y.
2015-12-01
The 26 December 2004 Sumatra-Andaman earthquake with moment magnitude (Mw) of 9.1 to 9.3 is the first great earthquake recorded by digital broadband, high-dynamic-range seismometers and global positioning system (GPS) equipment, which recorded many high-quality geophysical data sets. The spherical curvature is not negligible in far field especially for large event and the real Earth is laterally inhomogeneity and the analytical results still are difficult to explain the geodetic measurements. We use equivalent body force finite elements method Zhang et al. (2015) and mesh the whole earth, to compute global co-seismic displacements using four fault slip models of the 2004 Sumatra earthquake provided by different authors. Comparisons of calculated co-seismic displacements and GPS show that the confidences are well in near field for four models, and the confidences are according to different models. In the whole four models, the Chlieh model (Chlieh et al., 2007) is the best as this slip model not only accord well with near field data but also far field data. And then we use the best slip model, Chlieh model to explore influence of three dimensional lateral earth structure on both layered spherically symmetric (PREM) and real 3-D heterogeneous earth model (Crust 1.0 model and GyPSuM). Results show that the effects of 3-D heterogeneous earth model are not negligible and decrease concomitantly with increasing distance from the epicenter. The relative effects of 3-D crust model are 23% and 40% for horizontal and vertical displacements, respectively. The effects of the 3-D mantle model are much smaller than that of 3-D crust model but with wider impacting area.
A 3D finite element ALE method using an approximate Riemann solution
Chiravalle, V. P.; Morgan, N. R.
2016-08-09
Arbitrary Lagrangian–Eulerian finite volume methods that solve a multidimensional Riemann-like problem at the cell center in a staggered grid hydrodynamic (SGH) arrangement have been proposed. This research proposes a new 3D finite element arbitrary Lagrangian–Eulerian SGH method that incorporates a multidimensional Riemann-like problem. Here, two different Riemann jump relations are investigated. A new limiting method that greatly improves the accuracy of the SGH method on isentropic flows is investigated. A remap method that improves upon a well-known mesh relaxation and remapping technique in order to ensure total energy conservation during the remap is also presented. Numerical details and test problemmore » results are presented.« less
Parallel 3D Finite Element Numerical Modelling of DC Electron Guns
Prudencio, E.; Candel, A.; Ge, L.; Kabel, A.; Ko, K.; Lee, L.; Li, Z.; Ng, C.; Schussman, G.; /SLAC
2008-02-04
In this paper we present Gun3P, a parallel 3D finite element application that the Advanced Computations Department at the Stanford Linear Accelerator Center is developing for the analysis of beam formation in DC guns and beam transport in klystrons. Gun3P is targeted specially to complex geometries that cannot be described by 2D models and cannot be easily handled by finite difference discretizations. Its parallel capability allows simulations with more accuracy and less processing time than packages currently available. We present simulation results for the L-band Sheet Beam Klystron DC gun, in which case Gun3P is able to reduce simulation time from days to some hours.
NASA Astrophysics Data System (ADS)
Hu, Shengsun; Guo, Chaobo; Wang, Dongpo; Wang, Zhijiang
2016-07-01
The nonuniform distributions of the residual stress were simulated by a 3D finite element model to analyze the elastic-plastic dynamic ultrasonic impact treatment (UIT) process of multiple impacts on the 2024 aluminum alloy. The evolution of the stress during the impact process was discussed. The successive impacts during the UIT process improve the uniformity of the plastic deformation and decrease the maximum compressive residual stress beneath the former impact indentations. The influences of different controlled parameters, including the initial impact velocity, pin diameter, pin tip, device moving, and offset distances, on the residual stress distributions were analyzed. The influences of the controlled parameters on the residual stress distributions are apparent in the offset direction due to the different surface coverage in different directions. The influences can be used to understand the UIT process and to obtain the desired residual stress by optimizing the controlled parameters.
El-Anwar, Mohamed; Ghali, Rami; Aboelnagga, Mona
2016-01-01
AIM: This study aimed to estimate the stress patterns induced by the masticatory loads on a removable prosthesis supported and retained by bar splinted implants placed in the reconstructed mandible with two different clip materials and without clip, in the fibula-jaw bone and prosthesis using finite element analysis. METHODS: Two 3D finite element models were constructed, that models components were modeled on commercial CAD/CAM software then assembled into finite element package. Vertical loads were applied simulating the masticatory forces unilaterally in the resected site and bilaterally in the central fossa of the lower first molar as 100N (tension and compression). Analysis was based on the assumption full osseointegration between different types of bones, and between implants and fibula while fixing the top surface of the TMJ in place. RESULTS: The metallic bar connecting the three implants is insensitive to the clips material. Its supporting implants showed typical behavior with maximum stress values at the neck region. Fibula and jaw bone showed stresses within physiologic, while clips material effect seems to be very small due to its relatively small size. CONCLUSION: Switching loading force direction from tensile to compression did-not change the stresses and deformations distribution, but reversed their sign from positive to negative. PMID:27275353
Description of a parallel, 3D, finite element, hydrodynamics-diffusion code
Milovich, J L; Prasad, M K; Shestakov, A I
1999-04-11
We describe a parallel, 3D, unstructured grid finite element, hydrodynamic diffusion code for inertial confinement fusion (ICF) applications and the ancillary software used to run it. The code system is divided into two entities, a controller and a stand-alone physics code. The code system may reside on different computers; the controller on the user's workstation and the physics code on a supercomputer. The physics code is composed of separate hydrodynamic, equation-of-state, laser energy deposition, heat conduction, and radiation transport packages and is parallelized for distributed memory architectures. For parallelization, a SPMD model is adopted; the domain is decomposed into a disjoint collection of subdomains, one per processing element (PE). The PEs communicate using MPI. The code is used to simulate the hydrodynamic implosion of a spherical bubble.
An augmented Lagrangian finite element formulation for 3D contact of biphasic tissues.
Guo, Hongqiang; Spilker, Robert L
2014-01-01
Biphasic contact analysis is essential to obtain a complete understanding of soft tissue biomechanics, and the importance of physiological structure on the joint biomechanics has long been recognised; however, up to date, there are no successful developments of biphasic finite element contact analysis for three-dimensional (3D) geometries of physiological joints. The aim of this study was to develop a finite element formulation for biphasic contact of 3D physiological joints. The augmented Lagrangian method was used to enforce the continuity of contact traction and fluid pressure across the contact interface. The biphasic contact method was implemented in the commercial software COMSOL Multiphysics 4.2(®) (COMSOL, Inc., Burlington, MA). The accuracy of the implementation was verified using 3D biphasic contact problems, including indentation with a flat-ended indenter and contact of glenohumeral cartilage layers. The ability of the method to model multibody biphasic contact of physiological joints was proved by a 3D knee model. The 3D biphasic finite element contact method developed in this study can be used to study the biphasic behaviours of the physiological joints. PMID:23181617
Multiple-mode Lamb wave scattering simulations using 3D elastodynamic finite integration technique.
Leckey, Cara A C; Rogge, Matthew D; Miller, Corey A; Hinders, Mark K
2012-02-01
We have implemented three-dimensional (3D) elastodynamic finite integration technique (EFIT) simulations to model Lamb wave scattering for two flaw-types in an aircraft-grade aluminum plate, a rounded rectangle flat-bottom hole and a disbond of the same shape. The plate thickness and flaws explored in this work include frequency-thickness regions where several Lamb wave modes exist and sometimes overlap in phase and/or group velocity. For the case of the flat-bottom hole the depth was incrementally increased to explore progressive changes in multiple-mode Lamb wave scattering due to the damage. The flat-bottom hole simulation results have been compared to experimental data and are shown to provide key insight for this well-defined experimental case by explaining unexpected results in experimental waveforms. For the rounded rectangle disbond flaw, which would be difficult to implement experimentally, we found that Lamb wave behavior differed significantly from the flat-bottom hole flaw. Most of the literature in this field is restricted to low frequency-thickness regions due to difficulties in interpreting data when multiple modes exist. We found that benchmarked 3D EFIT simulations can yield an understanding of scattering behavior for these higher frequency-thickness regions and in cases that would be difficult to set up experimentally. Additionally, our results show that 2D simulations would not have been sufficient for modeling the complicated scattering that occurred. PMID:21908011
Finite Element Code For 3D-Hydraulic Fracture Propagation Equations (3-layer).
Energy Science and Technology Software Center (ESTSC)
1992-03-24
HYFRACP3D is a finite element program for simulation of a pseudo three-dimensional fracture geometries with a two-dimensional planar solution. The model predicts the height, width and winglength over time for a hydraulic fracture propagating in a three-layered system of rocks with variable rock mechanics properties.
NASA Technical Reports Server (NTRS)
Vos, R. G.; Straayer, J. W.
1975-01-01
The BOPACE 3-D is a finite element computer program, which provides a general family of three-dimensional isoparametric solid elements, and includes a new algorithm for improving the efficiency of the elastic-plastic-creep solution procedure. Theoretical, user, and programmer oriented sections are presented to describe the program.
NASA Technical Reports Server (NTRS)
Zhao, W.; Newman, J. C., Jr.; Sutton, M. A.; Shivakumar, K. N.; Wu, X. R.
1995-01-01
Parallel with the work in Part-1, stress intensity factors for semi-elliptical surface cracks emanating from a circular hole are determined. The 3-D weight function method with the 3D finite element solutions for the uncracked stress distribution as in Part-1 is used for the analysis. Two different loading conditions, i.e. remote tension and wedge loading, are considered for a wide range in geometrical parameters. Both single and double surface cracks are studied and compared with other solutions available in the literature. Typical crack opening displacements are also provided.
2D-3D hybrid stabilized finite element method for tsunami runup simulations
NASA Astrophysics Data System (ADS)
Takase, S.; Moriguchi, S.; Terada, K.; Kato, J.; Kyoya, T.; Kashiyama, K.; Kotani, T.
2016-09-01
This paper presents a two-dimensional (2D)-three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier-Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.
2D-3D hybrid stabilized finite element method for tsunami runup simulations
NASA Astrophysics Data System (ADS)
Takase, S.; Moriguchi, S.; Terada, K.; Kato, J.; Kyoya, T.; Kashiyama, K.; Kotani, T.
2016-05-01
This paper presents a two-dimensional (2D)-three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier-Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.
Accurate, finite-volume methods for 3D MHD on unstructured Lagrangian meshes
Barnes, D.C.; Rousculp, C.L.
1998-10-01
Previous 2D methods for magnetohydrodynamics (MHD) have contributed both to development of core code capability and to physics applications relevant to AGEX pulsed-power experiments. This strategy is being extended to 3D by development of a modular extension of an ASCI code. Extension to 3D not only increases complexity by problem size, but also introduces new physics, such as magnetic helicity transport. The authors have developed a method which incorporates all known conservation properties into the difference scheme on a Lagrangian unstructured mesh. Because the method does not depend on the mesh structure, mesh refinement is possible during a calculation to prevent the well known problem of mesh tangling. Arbitrary polyhedral cells are decomposed into tetrahedrons. The action of the magnetic vector potential, A {center_dot} {delta}l, is centered on the edges of this extended mesh. For ideal flow, this maintains {del} {center_dot} B = 0 to round-off error. Vertex forces are derived by the variation of magnetic energy with respect to vertex positions, F = {minus}{partial_derivative}W{sub B}/{partial_derivative}r. This assures symmetry as well as magnetic flux, momentum, and energy conservation. The method is local so that parallelization by domain decomposition is natural for large meshes. In addition, a simple, ideal-gas, finite pressure term has been included. The resistive diffusion part is calculated using the support operator method, to obtain an energy conservative, symmetric method on an arbitrary mesh. Implicit time difference equations are solved by preconditioned, conjugate gradient methods. Results of convergence tests are presented. Initial results of an annular Z-pinch implosion problem illustrate the application of these methods to multi-material problems.
3D finite element analysis of porous Ti-based alloy prostheses.
Mircheski, Ile; Gradišar, Marko
2016-11-01
In this paper, novel designs of porous acetabular cups are created and tested with 3D finite element analysis (FEA). The aim is to develop a porous acetabular cup with low effective radial stiffness of the structure, which will be near to the architectural and mechanical behavior of the natural bone. For the realization of this research, a 3D-scanner technology was used for obtaining a 3D-CAD model of the pelvis bone, a 3D-CAD software for creating a porous acetabular cup, and a 3D-FEA software for virtual testing of a novel design of the porous acetabular cup. The results obtained from this research reveal that a porous acetabular cup from Ti-based alloys with 60 ± 5% porosity has the mechanical behavior and effective radial stiffness (Young's modulus in radial direction) that meet and exceed the required properties of the natural bone. The virtual testing with 3D-FEA of a novel design with porous structure during the very early stage of the design and the development of orthopedic implants, enables obtaining a new or improved biomedical implant for a relatively short time and reduced price. PMID:27015664
2D to 3D conversion implemented in different hardware
NASA Astrophysics Data System (ADS)
Ramos-Diaz, Eduardo; Gonzalez-Huitron, Victor; Ponomaryov, Volodymyr I.; Hernandez-Fragoso, Araceli
2015-02-01
Conversion of available 2D data for release in 3D content is a hot topic for providers and for success of the 3D applications, in general. It naturally completely relies on virtual view synthesis of a second view given by original 2D video. Disparity map (DM) estimation is a central task in 3D generation but still follows a very difficult problem for rendering novel images precisely. There exist different approaches in DM reconstruction, among them manually and semiautomatic methods that can produce high quality DMs but they demonstrate hard time consuming and are computationally expensive. In this paper, several hardware implementations of designed frameworks for an automatic 3D color video generation based on 2D real video sequence are proposed. The novel framework includes simultaneous processing of stereo pairs using the following blocks: CIE L*a*b* color space conversions, stereo matching via pyramidal scheme, color segmentation by k-means on an a*b* color plane, and adaptive post-filtering, DM estimation using stereo matching between left and right images (or neighboring frames in a video), adaptive post-filtering, and finally, the anaglyph 3D scene generation. Novel technique has been implemented on DSP TMS320DM648, Matlab's Simulink module over a PC with Windows 7, and using graphic card (NVIDIA Quadro K2000) demonstrating that the proposed approach can be applied in real-time processing mode. The time values needed, mean Similarity Structural Index Measure (SSIM) and Bad Matching Pixels (B) values for different hardware implementations (GPU, Single CPU, and DSP) are exposed in this paper.
Accurate Finite Difference Algorithms
NASA Technical Reports Server (NTRS)
Goodrich, John W.
1996-01-01
Two families of finite difference algorithms for computational aeroacoustics are presented and compared. All of the algorithms are single step explicit methods, they have the same order of accuracy in both space and time, with examples up to eleventh order, and they have multidimensional extensions. One of the algorithm families has spectral like high resolution. Propagation with high order and high resolution algorithms can produce accurate results after O(10(exp 6)) periods of propagation with eight grid points per wavelength.
Application of edge-based finite elements and vector ABCs in 3D scattering
NASA Technical Reports Server (NTRS)
Chatterjee, A.; Jin, J. M.; Volakis, John L.
1992-01-01
A finite element absorbing boundary condition (FE-ABC) solution of the scattering by arbitrary 3-D structures is considered. The computational domain is discretized using edge-based tetrahedral elements. In contrast to the node-based elements, edge elements can treat geometries with sharp edges, are divergence-less, and easily satisfy the field continuity condition across dielectric interfaces. They do, however, lead to a higher unknown count but this is balanced by the greater sparsity of the resulting finite element matrix. Thus, the computation time required to solve such a system iteratively with a given degree of accuracy is less than the traditional node-based approach. The purpose is to examine the derivation and performance of the ABC's when applied to 2-D and 3-D problems and to discuss the specifics of our FE-ABC implementation.
Numerical solution of 3-D magnetotelluric using vector finite element method
NASA Astrophysics Data System (ADS)
Prihantoro, Rudy; Sutarno, Doddy; Nurhasan
2015-09-01
Magnetotelluric (MT) is a passive electromagnetic (EM) method which measure natural variations of electric and magnetic vector fields at the Earth surface to map subsurface electrical conductivity/resistivity structure. In this study, we obtained numerical solution of three-dimensional (3-D) MT using vector finite element method by solving second order Maxwell differential equation describing diffusion of plane wave through the conductive earth. Rather than the nodes of the element, the edges of the element is used as a vector basis to overcome the occurrence of nonphysical solutions that usually faced by scalar (node based) finite element method. Electric vector fields formulation was used and the resulting system of equation was solved using direct solution method to obtain the electric vector field distribution throughout the earth resistivity model structure. The resulting MT response functions was verified with 1-D layered Earth and 3-D2 COMMEMI outcropping structure. Good agreement is achieved for both structure models.
3D modeling of high-Tc superconductors by finite element software
NASA Astrophysics Data System (ADS)
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-Tc 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.
Application of 3D X-ray CT data sets to finite element analysis
Bossart, P.L.; Martz, H.E.; Brand, H.R.; Hollerbach, K.
1995-08-31
Finite Element Modeling (FEM) is becoming more important as industry drives toward concurrent engineering. A fundamental hindrance to fully exploiting the power of FEM is the human effort required to acquire complex part geometry, particularly as-built geometry, as a FEM mesh. Many Quantitative Non Destructive Evaluation (QNDE) techniques that produce three-dimensional (3D) data sets provide a substantial reduction in the effort required to apply FEM to as-built parts. This paper describes progress at LLNL on the application of 3D X-ray computed tomography (CT) data sets to more rapidly produce high-quality FEM meshes of complex, as-built geometries. Issues related to the volume segmentation of the 3D CT data as well as the use of this segmented data to tailor generic hexahedral FEM meshes to part specific geometries are discussed. The application of these techniques to FEM analysis in the medical field is reported here.
The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications
NASA Technical Reports Server (NTRS)
Bravo, Ramiro H.; Chen, Ching-Jen
1992-01-01
In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.
A Multi-Compartment 3-D Finite Element Model of Rectocele and Its Interaction with Cystocele
Luo, Jiajia; Chen, Luyun; Fenner, Dee E.; Ashton-Miller, James A.; DeLancey, John O. L.
2015-01-01
We developed a subject-specific 3-D finite element model to understand the mechanics underlying formation of female pelvic organ prolapse, specifically a rectocele and its interaction with a cystocele. The model was created from MRI 3-D geometry of a healthy 45 year-old multiparous woman. It included anterior and posterior vaginal walls, levator ani muscle, cardinal and uterosacral ligaments, anterior and posterior arcus tendineus fascia pelvis, arcus tendineus levator ani, perineal body, perineal membrane and anal sphincter. Material properties were mostly from the literature. Tissue impairment was modeled as decreased tissue stiffness based on previous clinical studies. Model equations were solved using Abaqus v 6.11. The sensitivity of anterior and posterior vaginal wall geometry was calculated for different combinations tissue impairments under increasing intraabdominal pressure. Prolapse size was reported as POP-Q point at point Bp for rectocele and point Ba for cystocele. Results show that a rectocele resulted from impairments of the levator ani and posterior compartment support. For 20% levator and 85% posterior support impairments, simulated rectocele size (at POP-Q point: Bp) increased 0.29 mm/cm H2O without apical impairment and 0.36 mm/cm H2O with 60% apical impairment, as intraabdominal pressures increased from 0 to 150 cm H2O. Apical support impairment could result in the development of either a cystocele or rectocele. Simulated repair of posterior compartment support decreased rectocele but increased a preexisting cystocele. We conclude that development of rectocele and cystocele depend on the presence of anterior, posterior, levator and/or or apical support impairments, as well as the interaction of the prolapse with the opposing compartment. PMID:25757664
Development of a 3D finite element model of lens microcirculation
2012-01-01
Background It has been proposed that in the absence of a blood supply, the ocular lens operates an internal microcirculation system. This system delivers nutrients, removes waste products and maintains ionic homeostasis in the lens. The microcirculation is generated by spatial differences in membrane transport properties; and previously has been modelled by an equivalent electrical circuit and solved analytically. While effective, this approach did not fully account for all the anatomical and functional complexities of the lens. To encapsulate these complexities we have created a 3D finite element computer model of the lens. Methods Initially, we created an anatomically-correct representative mesh of the lens. We then implemented the Stokes and advective Nernst-Plank equations, in order to model the water and ion fluxes respectively. Next we complemented the model with experimentally-measured surface ionic concentrations as boundary conditions and solved it. Results Our model calculated the standing ionic concentrations and electrical potential gradients in the lens. Furthermore, it generated vector maps of intra- and extracellular space ion and water fluxes that are proposed to circulate throughout the lens. These fields have only been measured on the surface of the lens and our calculations are the first 3D representation of their direction and magnitude in the lens. Conclusion Values for steady state standing fields for concentration and electrical potential plus ionic and fluid fluxes calculated by our model exhibited broad agreement with observed experimental values. Our model of lens function represents a platform to integrate new experimental data as they emerge and assist us to understand how the integrated structure and function of the lens contributes to the maintenance of its transparency. PMID:22992294
A 3-D implicit finite-volume model of shallow water flows
NASA Astrophysics Data System (ADS)
Wu, Weiming; Lin, Qianru
2015-09-01
A three-dimensional (3-D) model has been developed to simulate shallow water flows in large water bodies, such as coastal and estuarine waters. The eddy viscosity is determined using a newly modified mixing length model that uses different mixing length functions for the horizontal and vertical shear strain rates. The 3-D shallow water flow equations with the hydrostatic pressure assumption are solved using an implicit finite-volume method based on a quadtree (telescoping) rectangular mesh on the horizontal plane and the sigma coordinate in the vertical direction. The quadtree technique can locally refine the mesh around structures or in high-gradient regions by splitting a coarse cell into four child cells. The grid nodes are numbered with a one-dimensional index system that has unstructured grid feature for better grid flexibility. All the primary variables are arranged in a non-staggered grid system. Fluxes at cell faces are determined using a Rhie and Chow-type momentum interpolation, to avoid the possible spurious checkerboard oscillations caused by linear interpolation. Each of the discretized governing equations is solved iteratively using the flexible GMRES method with ILUT preconditioning, and coupling of water level and velocity among these equations is achieved by using the SIMPLEC algorithm with under-relaxation. The model has been tested in four cases, including steady flow near a spur-dyke, tidal flows in San Francisco Bay and Gironde Estuary, and wind-induced current in a flume. The calculated water levels and velocities are in good agreement with the measured values.
NASA Astrophysics Data System (ADS)
Komatitsch, Dimitri; Michéa, David; Erlebacher, Gordon; Göddeke, Dominik
2010-05-01
We first accelerate a three-dimensional finite-difference in the time domain (FDTD) wave propagation code by a factor of about 50 using Graphics Processing Unit (GPU) computing on a cheap NVIDIA graphics card with the CUDA programming language. We implement the code in CUDA in the case of the fully heterogeneous elastic wave equation. We also implement Convolution Perfectly Matched Layers (CPMLs) on the graphics card to efficiently absorb outgoing waves on the fictitious edges of the grid. We show that the code that runs on the graphics card gives the expected results by comparing our results to those obtained by running the same simulation on a classical processor core. The methodology that we present can be used for Maxwell's equations as well because their form is similar to that of the seismic wave equation written in velocity vector and stress tensor. We then implement a high-order finite-element (spectral-element) application, which performs the numerical simulation of seismic wave propagation resulting for instance from earthquakes at the scale of a continent or from active seismic acquisition experiments in the oil industry, on a cluster of NVIDIA Tesla graphics cards using the CUDA programming language and non blocking message passing based on MPI. We compare it to the implementation in C language and MPI on a classical cluster of CPU nodes. We use mesh coloring to efficiently handle summation operations over degrees of freedom on an unstructured mesh, and we exchange information between nodes using non blocking MPI messages. Using non-blocking communications allows us to overlap the communications across the network and the data transfer between the GPU card and the CPU node on which it is installed with calculations on that GPU card. We perform a number of numerical tests to validate the single-precision CUDA and MPI implementation and assess its accuracy. We then analyze performance measurements and in average we obtain a speedup of 20x to 25x.
Validation Studies of the Finite Orbit Width version of the CQL3D code
NASA Astrophysics Data System (ADS)
Petrov, Yu. V.; Harvey, R. W.
2014-10-01
The Finite-Orbit-Width (FOW) version of the CQL3D bounce-averaged Fokker-Planck (FP) code has been further developed and tested. The neoclassical radial transport appears naturally in this version by averaging the local collision coefficients along guiding center orbits, with a proper transformation matrix from local (R,Z) coordinates to the midplane computational coordinates, where the FP equation is solved. In a similar way, the local quasilinear rf diffusion terms give rise to additional radial transport of orbits. The main challenge is the internal boundary conditions (IBC) which add many elements into the matrix of coefficients for the solution of FPE on the computational grid, effectively making it a non-banded matrix (but still sparse). Steady state runs have been achieved at NERSC supercomputers in typically 10 time steps. Validation tests are performed for NSTX conditions, but using different scaling factors of equilibrium magnetic field, from 0.5 to 8.0. The bootstrap current calculations for ions show a reasonable agreement of current density profiles with Sauter et al. model equations which are based on 1st order expansion, although the magnitudes of currents may differ by up to 30%. Supported by USDOE grants SC0006614, ER54744, and ER44649.
Finite Element Analysis of Thermo-Mechanical Properties of 3D Braided Composites
NASA Astrophysics Data System (ADS)
Jiang, Li-li; Xu, Guo-dong; Cheng, Su; Lu, Xia-mei; Zeng, Tao
2014-04-01
This paper presents a modified finite element model (FEM) to investigate the thermo-mechanical properties of three-dimensional (3D) braided composite. The effective coefficients of thermal expansion (CTE) and the meso-scale mechanical response of 3D braided composites are predicted. The effects of the braiding angle and fiber volume fraction on the effective CTE are evaluated. The results are compared to the experimental data available in the literature to demonstrate the accuracy and reliability of the present method. The tensile stress distributions of the representative volume element (RVE) are also outlined. It is found that the stress of the braiding yarn has a significant increase with temperature rise; on the other hand, the temperature change has an insignificant effect on the stress of the matrix. In addition, a rapid decrease in the tensile strength of 3D braided composites is observed with the increase in temperature. It is revealed that the thermal conditions have a significant effect on the strength of 3D braided composites. The present method provides an effective tool to predict the stresses of 3D braided composites under thermo-mechanical loading.
NASA Technical Reports Server (NTRS)
Raju, I. S.; Newman, J. C., Jr.
1993-01-01
A computer program, surf3d, that uses the 3D finite-element method to calculate the stress-intensity factors for surface, corner, and embedded cracks in finite-thickness plates with and without circular holes, was developed. The cracks are assumed to be either elliptic or part eliptic in shape. The computer program uses eight-noded hexahedral elements to model the solid. The program uses a skyline storage and solver. The stress-intensity factors are evaluated using the force method, the crack-opening displacement method, and the 3-D virtual crack closure methods. In the manual the input to and the output of the surf3d program are described. This manual also demonstrates the use of the program and describes the calculation of the stress-intensity factors. Several examples with sample data files are included with the manual. To facilitate modeling of the user's crack configuration and loading, a companion program (a preprocessor program) that generates the data for the surf3d called gensurf was also developed. The gensurf program is a three dimensional mesh generator program that requires minimal input and that builds a complete data file for surf3d. The program surf3d is operational on Unix machines such as CRAY Y-MP, CRAY-2, and Convex C-220.
Charged-particle Gun Design with 3D Finite-element Methods
NASA Astrophysics Data System (ADS)
Humphries, Stanley
2002-04-01
The DARHT second-axis injector poses a major challenge for computer simulation. The relativistic electrons are subject to strong beam-generated electric and magnetic forces. The beam and applied fields are fully three-dimensional. Furthermore, accurate field calculations at surfaces are critical to model Child-law emission. Although several 2D relativistic beam codes are available, there is presently no 3D tool that can address all important processes in the DARHT injector. As a result, we created the OmniTrak 3D finite-element code suite. This talk gives a basic tutorial on finite-element methods with emphasis on electron gun design via the ray-tracing technique. Four main areas are covered: 1) the mesh as a tool to organize space, 2) transformation of the Poisson equation through the minimum residual principle, 3) orbit tracking in a complex environment and 4) handling self-consistent beam-generated fields. The components of a volume mesh (elements, nodes and facets) are reviewed. We consider motivations for choosing a 3D mesh style: structured versus unstructured, tetrahedrons versus hexahedrons. We discuss methods for taking volume integrals over arbitrary hexahedrons through normal coordinates and shape functions, leading to the fundamental field equations. The special problems of 3D magnetic field solutions and the advantages of the reduced potential method are outlined. Accurate field interpolations for orbit calculations require fast identification of occupied elements. A method for fast element identification that also yields the orbit penetration point on the element surface is described. The final topics are the assignment of charge and current to meshes from calculated orbits and techniques for space-charge-limited emission from multiple arbitrary 3D surfaces.
ATHENA 3D: A finite element code for ultrasonic wave propagation
NASA Astrophysics Data System (ADS)
Rose, C.; Rupin, F.; Fouquet, T.; Chassignole, B.
2014-04-01
The understanding of wave propagation phenomena requires use of robust numerical models. 3D finite element (FE) models are generally prohibitively time consuming. However, advances in computing processor speed and memory allow them to be more and more competitive. In this context, EDF R&D developed the 3D version of the well-validated FE code ATHENA2D. The code is dedicated to the simulation of wave propagation in all kinds of elastic media and in particular, heterogeneous and anisotropic materials like welds. It is based on solving elastodynamic equations in the calculation zone expressed in terms of stress and particle velocities. The particularity of the code relies on the fact that the discretization of the calculation domain uses a Cartesian regular 3D mesh while the defect of complex geometry can be described using a separate (2D) mesh using the fictitious domains method. This allows combining the rapidity of regular meshes computation with the capability of modelling arbitrary shaped defects. Furthermore, the calculation domain is discretized with a quasi-explicit time evolution scheme. Thereby only local linear systems of small size have to be solved. The final step to reduce the computation time relies on the fact that ATHENA3D has been parallelized and adapted to the use of HPC resources. In this paper, the validation of the 3D FE model is discussed. A cross-validation of ATHENA 3D and CIVA is proposed for several inspection configurations. The performances in terms of calculation time are also presented in the cases of both local computer and computation cluster use.
Finite-element 3D simulation tools for high-current relativistic electron beams
NASA Astrophysics Data System (ADS)
Humphries, Stanley; Ekdahl, Carl
2002-08-01
The DARHT second-axis injector is a challenge for computer simulations. Electrons are subject to strong beam-generated forces. The fields are fully three-dimensional and accurate calculations at surfaces are critical. We describe methods applied in OmniTrak, a 3D finite-element code suite that can address DARHT and the full range of charged-particle devices. The system handles mesh generation, electrostatics, magnetostatics and self-consistent particle orbits. The MetaMesh program generates meshes of conformal hexahedrons to fit any user geometry. The code has the unique ability to create structured conformal meshes with cubic logic. Organized meshes offer advantages in speed and memory utilization in the orbit and field solutions. OmniTrak is a versatile charged-particle code that handles 3D electric and magnetic field solutions on independent meshes. The program can update both 3D field solutions from the calculated beam space-charge and current-density. We shall describe numerical methods for orbit tracking on a hexahedron mesh. Topics include: 1) identification of elements along the particle trajectory, 2) fast searches and adaptive field calculations, 3) interpolation methods to terminate orbits on material surfaces, 4) automatic particle generation on multiple emission surfaces to model space-charge-limited emission and field emission, 5) flexible Child law algorithms, 6) implementation of the dual potential model for 3D magnetostatics, and 7) assignment of charge and current from model particle orbits for self-consistent fields.
Calculation by the finite element method of 3-D turbulent flow in a centrifugal pump
NASA Astrophysics Data System (ADS)
Combes, J. F.
1992-02-01
In order to solve industrial flow problems in complex geometries, a finite element code, N3S, was developed. It allows the computation of a wide variety of 2-D or 3-D unsteady incompressible flows, by solving the Reynolds averaged Navier-Stokes equations together with a k-epsilon turbulence model. Some recent developments of this code concern turbomachinery flows, where one has to take into account periodic boundary conditions, as well as Coriolis and centrifugal forces. The numerical treatment is based on a fractional step method: at each time step, an advection step is solved successively by means of a characteristic method; a diffusion step for the scalar terms; and finally, a Generalized Stokes Problem by using a preconditioned Uzawa algorithm. The space discretization uses a standard Galerkin finite element method with a mixed formulation for the velocity and pressure. An application is presented of this code to the flow inside a centrifugal pump which was extensively tested on several air and water test rigs, and for which many quasi-3-D or Euler calculations were reported. The present N3S calculation is made on a finite element mesh comprising about 28000 tetrahedrons and 43000 nodes.
NASA Technical Reports Server (NTRS)
Raju, I. S.
1992-01-01
A computer program that generates three-dimensional (3D) finite element models for cracked 3D solids was written. This computer program, gensurf, uses minimal input data to generate 3D finite element models for isotropic solids with elliptic or part-elliptic cracks. These models can be used with a 3D finite element program called surf3d. This report documents this mesh generator. In this manual the capabilities, limitations, and organization of gensurf are described. The procedures used to develop 3D finite element models and the input for and the output of gensurf are explained. Several examples are included to illustrate the use of this program. Several input data files are included with this manual so that the users can edit these files to conform to their crack configuration and use them with gensurf.
Justification for a 2D versus 3D fingertip finite element model during static contact simulations.
Harih, Gregor; Tada, Mitsunori; Dolšak, Bojan
2016-10-01
The biomechanical response of a human hand during contact with various products has not been investigated in details yet. It has been shown that excessive contact pressure on the soft tissue can result in discomfort, pain and also cumulative traumatic disorders. This manuscript explores the benefits and limitations of a simplified two-dimensional vs. an anatomically correct three-dimensional finite element model of a human fingertip. Most authors still use 2D FE fingertip models due to their simplicity and reduced computational costs. However we show that an anatomically correct 3D FE fingertip model can provide additional insight into the biomechanical behaviour. The use of 2D fingertip FE models is justified when observing peak contact pressure values as well as displacement during the contact for the given studied cross-section. On the other hand, an anatomically correct 3D FE fingertip model provides a contact pressure distribution, which reflects the fingertip's anatomy. PMID:26856769
NASA Technical Reports Server (NTRS)
Nakazawa, S.
1988-01-01
This annual status report presents the results of work performed during the fourth year of the 3-D Inelastic Analysis Methods for Hot Section Components program (NASA Contract NAS3-23697). The objective of the program is to produce a series of new computer codes permitting more accurate and efficient 3-D analysis of selected hot section components, i.e., combustor liners, turbine blades and turbine vanes. The computer codes embody a progression of math models and are streamlined to take advantage of geometrical features, loading conditions, and forms of material response that distinguish each group of selected components. Volume 1 of this report discusses the special finite element models developed during the fourth year of the contract.
3D imaging of soil pore network: two different approaches
NASA Astrophysics Data System (ADS)
Matrecano, M.; Di Matteo, B.; Mele, G.; Terribile, F.
2009-04-01
Pore geometry imaging and its quantitative description is a key factor for advances in the knowledge of physical, chemical and biological soil processes. For many years photos from flattened surfaces of undisturbed soil samples impregnated with fluorescent resin and from soil thin sections under microscope have been the only way available for exploring pore architecture at different scales. Earlier 3D representations of the internal structure of the soil based on not destructive methods have been obtained using medical tomographic systems (NMR and X-ray CT). However, images provided using such equipments, show strong limitations in terms of spatial resolution. In the last decade very good results have then been obtained using imaging from very expensive systems based on synchrotron radiation. More recently, X-ray Micro-Tomography has resulted the most widely applied being the technique showing the best compromise between costs, resolution and size of the images. Conversely, the conceptually simpler but destructive method of "serial sectioning" has been progressively neglected for technical problems in sample preparation and time consumption needed to obtain an adequate number of serial sections for correct 3D reconstruction of soil pore geometry. In this work a comparison between the two methods above has been carried out in order to define advantages, shortcomings and to point out their different potential. A cylindrical undisturbed soil sample 6.5cm in diameter and 6.5cm height of an Ap horizon of an alluvial soil showing vertic characteristics, has been reconstructed using both a desktop X-ray micro-tomograph Skyscan 1172 and the new automatic serial sectioning system SSAT (Sequential Section Automatic Tomography) set up at CNR ISAFOM in Ercolano (Italy) with the aim to overcome most of the typical limitations of such a technique. Image best resolution of 7.5 µm per voxel resulted using X-ray Micro CT while 20 µm was the best value using the serial sectioning
Dynamic Analysis of 2D Electromagnetic Resonant Optical Scanner Using 3D Finite Element Method
NASA Astrophysics Data System (ADS)
Hirata, Katsuhiro; Hong, Sara; Maeda, Kengo
The optical scanner is a scanning device in which a laser beam is reflected by a mirror that can be rotated or oscillated. In this paper, we propose a new 2D electromagnetic resonant optical scanner that employs electromagnets and leaf springs. Torque characteristics and resonance characteristics of the scanner are analyzed using the 3D finite element method. The validity of the analysis is shown by comparing the characteristics inferred from the analysis with the characteristics of the prototype. Further, 2D resonance is investigated by introducing a superimposed-frequency current in a single coil.
Parallel goal-oriented adaptive finite element modeling for 3D electromagnetic exploration
NASA Astrophysics Data System (ADS)
Zhang, Y.; Key, K.; Ovall, J.; Holst, M.
2014-12-01
We present a parallel goal-oriented adaptive finite element method for accurate and efficient electromagnetic (EM) modeling of complex 3D structures. An unstructured tetrahedral mesh allows this approach to accommodate arbitrarily complex 3D conductivity variations and a priori known boundaries. The total electric field is approximated by the lowest order linear curl-conforming shape functions and the discretized finite element equations are solved by a sparse LU factorization. Accuracy of the finite element solution is achieved through adaptive mesh refinement that is performed iteratively until the solution converges to the desired accuracy tolerance. Refinement is guided by a goal-oriented error estimator that uses a dual-weighted residual method to optimize the mesh for accurate EM responses at the locations of the EM receivers. As a result, the mesh refinement is highly efficient since it only targets the elements where the inaccuracy of the solution corrupts the response at the possibly distant locations of the EM receivers. We compare the accuracy and efficiency of two approaches for estimating the primary residual error required at the core of this method: one uses local element and inter-element residuals and the other relies on solving a global residual system using a hierarchical basis. For computational efficiency our method follows the Bank-Holst algorithm for parallelization, where solutions are computed in subdomains of the original model. To resolve the load-balancing problem, this approach applies a spectral bisection method to divide the entire model into subdomains that have approximately equal error and the same number of receivers. The finite element solutions are then computed in parallel with each subdomain carrying out goal-oriented adaptive mesh refinement independently. We validate the newly developed algorithm by comparison with controlled-source EM solutions for 1D layered models and with 2D results from our earlier 2D goal oriented
Isoparametric 3-D Finite Element Mesh Generation Using Interactive Computer Graphics
NASA Technical Reports Server (NTRS)
Kayrak, C.; Ozsoy, T.
1985-01-01
An isoparametric 3-D finite element mesh generator was developed with direct interface to an interactive geometric modeler program called POLYGON. POLYGON defines the model geometry in terms of boundaries and mesh regions for the mesh generator. The mesh generator controls the mesh flow through the 2-dimensional spans of regions by using the topological data and defines the connectivity between regions. The program is menu driven and the user has a control of element density and biasing through the spans and can also apply boundary conditions, loads interactively.
Nonstandard finite difference schemes
NASA Technical Reports Server (NTRS)
Mickens, Ronald E.
1995-01-01
The major research activities of this proposal center on the construction and analysis of nonstandard finite-difference schemes for ordinary and partial differential equations. In particular, we investigate schemes that either have zero truncation errors (exact schemes) or possess other significant features of importance for numerical integration. Our eventual goal is to bring these methods to bear on problems that arise in the modeling of various physical, engineering, and technological systems. At present, these efforts are extended in the direction of understanding the exact nature of these nonstandard procedures and extending their use to more complicated model equations. Our presentation will give a listing (obtained to date) of the nonstandard rules, their application to a number of linear and nonlinear, ordinary and partial differential equations. In certain cases, numerical results will be presented.
The numerical integration and 3-D finite element formulation of a viscoelastic model of glass
Chambers, R.S.
1994-08-01
The use of glasses is widespread in making hermetic, insulating seals for many electronic components. Flat panel displays and fiber optic connectors are other products utilizing glass as a structural element. When glass is cooled from sealing temperatures, residual stresses are generated due to mismatches in thermal shrinkage created by the dissimilar material properties of the adjoining materials. Because glass is such a brittle material at room temperature, tensile residual stresses must be kept small to ensure durability and avoid cracking. Although production designs and the required manufacturing process development can be deduced empirically, this is an expensive and time consuming process that does not necessarily lead to an optimal design. Agile manufacturing demands that analyses be used to reduce development costs and schedules by providing insight and guiding the design process through the development cycle. To make these gains, however, viscoelastic models of glass must be available along with the right tool to use them. A viscoelastic model of glass can be used to simulate the stress and volume relaxation that occurs at elevated temperatures as the molecular structure of the glass seeks to equilibrate to the state of the supercooled liquid. The substance of the numerical treatment needed to support the implementation of the model in a 3-D finite element program is presented herein. An accurate second-order, central difference integrator is proposed for the constitutive equations, and numerical solutions are compared to those obtained with other integrators. Inherent convergence problems are reviewed and fixes are described. The resulting algorithms are generally applicable to the broad class of viscoelastic material models. First-order error estimates are used as a basis for developing a scheme for automatic time step controls, and several demonstration problems are presented to illustrate the performance of the methodology.
Glacial isostatic adjustment on 3-D Earth models: a finite-volume formulation
NASA Astrophysics Data System (ADS)
Latychev, Konstantin; Mitrovica, Jerry X.; Tromp, Jeroen; Tamisiea, Mark E.; Komatitsch, Dimitri; Christara, Christina C.
2005-05-01
We describe and present results from a finite-volume (FV) parallel computer code for forward modelling the Maxwell viscoelastic response of a 3-D, self-gravitating, elastically compressible Earth to an arbitrary surface load. We implement a conservative, control volume discretization of the governing equations using a tetrahedral grid in Cartesian geometry and a low-order, linear interpolation. The basic starting grid honours all major radial discontinuities in the Preliminary Reference Earth Model (PREM), and the models are permitted arbitrary spatial variations in viscosity and elastic parameters. These variations may be either continuous or discontinuous at a set of grid nodes forming a 3-D surface within the (regional or global) modelling domain. In the second part of the paper, we adopt the FV methodology and a spherically symmetric Earth model to generate a suite of predictions sampling a broad class of glacial isostatic adjustment (GIA) data types (3-D crustal motions, long-wavelength gravity anomalies). These calculations, based on either a simple disc load history or a global Late Pleistocene ice load reconstruction (ICE-3G), are benchmarked against predictions generated using the traditional normal-mode approach to GIA. The detailed comparison provides a guide for future analyses (e.g. what grid resolution is required to obtain a specific accuracy?) and it indicates that discrepancies in predictions of 3-D crustal velocities less than 0.1 mm yr-1 are generally obtainable for global grids with ~3 × 106 nodes; however, grids of higher resolution are required to predict large-amplitude (>1 cm yr-1) radial velocities in zones of peak post-glacial uplift (e.g. James bay) to the same level of absolute accuracy. We conclude the paper with a first application of the new formulation to a 3-D problem. Specifically, we consider the impact of mantle viscosity heterogeneity on predictions of present-day 3-D crustal motions in North America. In these tests, the
Mimetic finite difference method
NASA Astrophysics Data System (ADS)
Lipnikov, Konstantin; Manzini, Gianmarco; Shashkov, Mikhail
2014-01-01
The mimetic finite difference (MFD) method mimics fundamental properties of mathematical and physical systems including conservation laws, symmetry and positivity of solutions, duality and self-adjointness of differential operators, and exact mathematical identities of the vector and tensor calculus. This article is the first comprehensive review of the 50-year long history of the mimetic methodology and describes in a systematic way the major mimetic ideas and their relevance to academic and real-life problems. The supporting applications include diffusion, electromagnetics, fluid flow, and Lagrangian hydrodynamics problems. The article provides enough details to build various discrete operators on unstructured polygonal and polyhedral meshes and summarizes the major convergence results for the mimetic approximations. Most of these theoretical results, which are presented here as lemmas, propositions and theorems, are either original or an extension of existing results to a more general formulation using polyhedral meshes. Finally, flexibility and extensibility of the mimetic methodology are shown by deriving higher-order approximations, enforcing discrete maximum principles for diffusion problems, and ensuring the numerical stability for saddle-point systems.
Comparison of different precondtioners for nonsymmtric finite volume element methods
Mishev, I.D.
1996-12-31
We consider a few different preconditioners for the linear systems arising from the discretization of 3-D convection-diffusion problems with the finite volume element method. Their theoretical and computational convergence rates are compared and discussed.
A least-squares finite element method for 3D incompressible Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Jiang, Bo-Nan; Lin, T. L.; Hou, Lin-Jun; Povinelli, Louis A.
1993-01-01
The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations, and results in symmetric, positive definite algebraic system. An additional compatibility equation, i.e., the divergence of vorticity vector should be zero, is included to make the first-order system elliptic. The Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. The flow in a half of 3D cubic cavity is calculated at Re = 100, 400, and 1,000 with 50 x 52 x 25 trilinear elements. The Taylor-Gortler-like vortices are observed at Re = 1,000.
An accurate quadrature technique for the contact boundary in 3D finite element computations
NASA Astrophysics Data System (ADS)
Duong, Thang X.; Sauer, Roger A.
2015-01-01
This paper presents a new numerical integration technique for 3D contact finite element implementations, focusing on a remedy for the inaccurate integration due to discontinuities at the boundary of contact surfaces. The method is based on the adaptive refinement of the integration domain along the boundary of the contact surface, and is accordingly denoted RBQ for refined boundary quadrature. It can be used for common element types of any order, e.g. Lagrange, NURBS, or T-Spline elements. In terms of both computational speed and accuracy, RBQ exhibits great advantages over a naive increase of the number of quadrature points. Also, the RBQ method is shown to remain accurate for large deformations. Furthermore, since the sharp boundary of the contact surface is determined, it can be used for various purposes like the accurate post-processing of the contact pressure. Several examples are presented to illustrate the new technique.
NASA Technical Reports Server (NTRS)
Nakazawa, S.
1987-01-01
This Annual Status Report presents the results of work performed during the third year of the 3-D Inelastic Analysis Methods for Hot Section Components program (NASA Contract NAS3-23697). The objective of the program is to produce a series of new computer codes that permit more accurate and efficient three-dimensional analysis of selected hot section components, i.e., combustor liners, turbine blades, and turbine vanes. The computer codes embody a progression of mathematical models and are streamlined to take advantage of geometrical features, loading conditions, and forms of material response that distinguish each group of selected components. This report is presented in two volumes. Volume 1 describes effort performed under Task 4B, Special Finite Element Special Function Models, while Volume 2 concentrates on Task 4C, Advanced Special Functions Models.
Ash3d: A finite-volume, conservative numerical model for ash transport and tephra deposition
Schwaiger, Hans F.; Denlinger, Roger P.; Mastin, Larry G.
2012-01-01
We develop a transient, 3-D Eulerian model (Ash3d) to predict airborne volcanic ash concentration and tephra deposition during volcanic eruptions. This model simulates downwind advection, turbulent diffusion, and settling of ash injected into the atmosphere by a volcanic eruption column. Ash advection is calculated using time-varying pre-existing wind data and a robust, high-order, finite-volume method. Our routine is mass-conservative and uses the coordinate system of the wind data, either a Cartesian system local to the volcano or a global spherical system for the Earth. Volcanic ash is specified with an arbitrary number of grain sizes, which affects the fall velocity, distribution and duration of transport. Above the source volcano, the vertical mass distribution with elevation is calculated using a Suzuki distribution for a given plume height, eruptive volume, and eruption duration. Multiple eruptions separated in time may be included in a single simulation. We test the model using analytical solutions for transport. Comparisons of the predicted and observed ash distributions for the 18 August 1992 eruption of Mt. Spurr in Alaska demonstrate to the efficacy and efficiency of the routine.
NASA Technical Reports Server (NTRS)
Nakazawa, Shohei
1989-01-01
The user options available for running the MHOST finite element analysis package is described. MHOST is a solid and structural analysis program based on the mixed finite element technology, and is specifically designed for 3-D inelastic analysis. A family of 2- and 3-D continuum elements along with beam and shell structural elements can be utilized, many options are available in the constitutive equation library, the solution algorithms and the analysis capabilities. The outline of solution algorithms is discussed along with the data input and output, analysis options including the user subroutines and the definition of the finite elements implemented in the program package.
Kolotilina, L.; Nikishin, A.; Yeremin, A.
1994-12-31
The solution of large systems of linear equations is a crucial bottleneck when performing 3D finite element analysis of structures. Also, in many cases the reliability and robustness of iterative solution strategies, and their efficiency when exploiting hardware resources, fully determine the scope of industrial applications which can be solved on a particular computer platform. This is especially true for modern vector/parallel supercomputers with large vector length and for modern massively parallel supercomputers. Preconditioned iterative methods have been successfully applied to industrial class finite element analysis of structures. The construction and application of high quality preconditioners constitutes a high percentage of the total solution time. Parallel implementation of high quality preconditioners on such architectures is a formidable challenge. Two common types of existing preconditioners are the implicit preconditioners and the explicit preconditioners. The implicit preconditioners (e.g. incomplete factorizations of several types) are generally high quality but require solution of lower and upper triangular systems of equations per iteration which are difficult to parallelize without deteriorating the convergence rate. The explicit type of preconditionings (e.g. polynomial preconditioners or Jacobi-like preconditioners) require sparse matrix-vector multiplications and can be parallelized but their preconditioning qualities are less than desirable. The authors present results of numerical experiments with Factorized Sparse Approximate Inverses (FSAI) for symmetric positive definite linear systems. These are high quality preconditioners that possess a large resource of parallelism by construction without increasing the serial complexity.
FERM3D: A finite element R-matrix electron molecule scattering code
NASA Astrophysics Data System (ADS)
Tonzani, Stefano
2007-01-01
FERM3D is a three-dimensional finite element program, for the elastic scattering of a low energy electron from a general polyatomic molecule, which is converted to a potential scattering problem. The code is based on tricubic polynomials in spherical coordinates. The electron-molecule interaction is treated as a sum of three terms: electrostatic, exchange, and polarization. The electrostatic term can be extracted directly from ab initio codes ( GAUSSIAN 98 in the work described here), while the exchange term is approximated using a local density functional. A local polarization potential based on density functional theory [C. Lee, W. Yang, R.G. Parr, Phys. Rev. B 37 (1988) 785] describes the long range attraction to the molecular target induced by the scattering electron. Photoionization calculations are also possible and illustrated in the present work. The generality and simplicity of the approach is important in extending electron-scattering calculations to more complex targets than it is possible with other methods. Program summaryTitle of program:FERM3D Catalogue identifier:ADYL_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYL_v1_0 Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Computer for which the program is designed and others on which it has been tested:Intel Xeon, AMD Opteron 64 bit, Compaq Alpha Operating systems or monitors under which the program has been tested:HP Tru64 Unix v5.1, Red Hat Linux Enterprise 3 Programming language used:Fortran 90 Memory required to execute with typical data:900 MB (neutral CO 2), 2.3 GB (ionic CO 2), 1.4 GB (benzene) No. of bits in a word:32 No. of processors used:1 Has the code been vectorized?:No No. of lines in distributed program, including test data, etc.:58 383 No. of bytes in distributed program, including test data, etc.:561 653 Distribution format:tar.gzip file CPC Program library subprograms used:ADDA, ACDP Nature of physical problem:Scattering of an
NASA Astrophysics Data System (ADS)
Moortgat, Joachim; Firoozabadi, Abbas
2016-06-01
Problems of interest in hydrogeology and hydrocarbon resources involve complex heterogeneous geological formations. Such domains are most accurately represented in reservoir simulations by unstructured computational grids. Finite element methods accurately describe flow on unstructured meshes with complex geometries, and their flexible formulation allows implementation on different grid types. In this work, we consider for the first time the challenging problem of fully compositional three-phase flow in 3D unstructured grids, discretized by any combination of tetrahedra, prisms, and hexahedra. We employ a mass conserving mixed hybrid finite element (MHFE) method to solve for the pressure and flux fields. The transport equations are approximated with a higher-order vertex-based discontinuous Galerkin (DG) discretization. We show that this approach outperforms a face-based implementation of the same polynomial order. These methods are well suited for heterogeneous and fractured reservoirs, because they provide globally continuous pressure and flux fields, while allowing for sharp discontinuities in compositions and saturations. The higher-order accuracy improves the modeling of strongly non-linear flow, such as gravitational and viscous fingering. We review the literature on unstructured reservoir simulation models, and present many examples that consider gravity depletion, water flooding, and gas injection in oil saturated reservoirs. We study convergence rates, mesh sensitivity, and demonstrate the wide applicability of our chosen finite element methods for challenging multiphase flow problems in geometrically complex subsurface media.
A NURBS-based generalized finite element scheme for 3D simulation of heterogeneous materials
NASA Astrophysics Data System (ADS)
Safdari, Masoud; Najafi, Ahmad R.; Sottos, Nancy R.; Geubelle, Philippe H.
2016-08-01
A 3D NURBS-based interface-enriched generalized finite element method (NIGFEM) is introduced to solve problems with complex discontinuous gradient fields observed in the analysis of heterogeneous materials. The method utilizes simple structured meshes of hexahedral elements that do not necessarily conform to the material interfaces in heterogeneous materials. By avoiding the creation of conforming meshes used in conventional FEM, the NIGFEM leads to significant simplification of the mesh generation process. To achieve an accurate solution in elements that are crossed by material interfaces, the NIGFEM utilizes Non-Uniform Rational B-Splines (NURBS) to enrich the solution field locally. The accuracy and convergence of the NIGFEM are tested by solving a benchmark problem. We observe that the NIGFEM preserves an optimal rate of convergence, and provides additional advantages including the accurate capture of the solution fields in the vicinity of material interfaces and the built-in capability for hierarchical mesh refinement. Finally, the use of the NIGFEM in the computational analysis of heterogeneous materials is discussed.
Finite element modeling of a 3D coupled foot-boot model.
Qiu, Tian-Xia; Teo, Ee-Chon; Yan, Ya-Bo; Lei, Wei
2011-12-01
Increasingly, musculoskeletal models of the human body are used as powerful tools to study biological structures. The lower limb, and in particular the foot, is of interest because it is the primary physical interaction between the body and the environment during locomotion. The goal of this paper is to adopt the finite element (FE) modeling and analysis approaches to create a state-of-the-art 3D coupled foot-boot model for future studies on biomechanical investigation of stress injury mechanism, foot wear design and parachute landing fall simulation. In the modeling process, the foot-ankle model with lower leg was developed based on Computed Tomography (CT) images using ScanIP, Surfacer and ANSYS. Then, the boot was represented by assembling the FE models of upper, insole, midsole and outsole built based on the FE model of the foot-ankle, and finally the coupled foot-boot model was generated by putting together the models of the lower limb and boot. In this study, the FE model of foot and ankle was validated during balance standing. There was a good agreement in the overall patterns of predicted and measured plantar pressure distribution published in literature. The coupled foot-boot model will be fully validated in the subsequent works under both static and dynamic loading conditions for further studies on injuries investigation in military and sports, foot wear design and characteristics of parachute landing impact in military. PMID:21676642
A new 3D finite element model of the IEC 60318-1 artificial ear
NASA Astrophysics Data System (ADS)
Bravo, Agustín; Barham, Richard; Ruiz, Mariano; López, Juan Manuel; DeArcas, Guillermo; Recuero, Manuel
2008-08-01
The artificial ear specified in IEC 60318-1 is used for the measurement of headphones and has been designed to present an acoustic load equivalent to that of normal human ears. In this respect it is specified in terms of an acoustical impedance, and modelled by a lumped parameter approach. However, this has some inherent frequency limitations and becomes less valid as the acoustic wavelength approaches the characteristic dimensions within the device. In addition, when sound propagates through structures such as narrow tubes, annular slits or over sharp corners, noticeable thermal and viscous effects take place causing further departure from the lumped parameter model. A new numerical model has therefore been developed, which gives proper consideration to the aforementioned effects. Both kinds of losses can be simulated by means of the LMS Virtual Lab acoustic software which facilitates finite and boundary element modelling of the whole artificial ear. A full 3D model of the artificial ear has therefore been developed based on key dimensional data found in IEC 60318-1. The model has been used to calculate the acoustical impedance, and the results compared with the corresponding data determined from the lumped parameter model. The numerical simulation of the artificial ear has been shown to provide realistic results, and is a powerful tool for developing a detailed understanding of the device. It is also proving valuable in the revision of IEC 60318-1 that is currently in progress.
3D finite element simulation of effects of deflection rate on energy absorption for TRIP steel
NASA Astrophysics Data System (ADS)
Hayashi, Asuka; Pham, Hang; Iwamoto, Takeshi
2015-09-01
Recently, with the requirement of lighter weight and more safety for a design of automobile, energy absorption capability of structural materials has become important. TRIP (Transformation-induced Plasticity) steel is expected to apply to safety members because of excellent energy absorption capability and ductility. Past studies proved that such excellent characteristics in TRIP steel are dominated by strain-induced martensitic transformation (SIMT) during plastic deformation. Because SIMT strongly depends on deformation rate and temperature, an investigation of the effects of deformation rate and temperature on energy absorption in TRIP is essential. Although energy absorption capability of material can be estimated by J-integral experimentally by using pre-cracked specimen, it is difficult to determine volume fraction of martensite and temperature rise during the crack extension. In addition, their effects on J-integral, especially at high deformation rate in experiment might be quite hard. Thus, a computational prediction needs to be performed. In this study, bending deformation behavior of pre-cracked specimen until the onset point of crack extension are predicted by 3D finite element simulation based on the transformation kinetics model proposed by Iwamoto et al. (1998). It is challenged to take effects of temperature, volume fraction of martensite and deformation rate into account. Then, the mechanism for higher energy absorption characteristic will be discussed.
A 2-D Interface Element for Coupled Analysis of Independently Modeled 3-D Finite Element Subdomains
NASA Technical Reports Server (NTRS)
Kandil, Osama A.
1998-01-01
Over the past few years, the development of the interface technology has provided an analysis framework for embedding detailed finite element models within finite element models which are less refined. This development has enabled the use of cascading substructure domains without the constraint of coincident nodes along substructure boundaries. The approach used for the interface element is based on an alternate variational principle often used in deriving hybrid finite elements. The resulting system of equations exhibits a high degree of sparsity but gives rise to a non-positive definite system which causes difficulties with many of the equation solvers in general-purpose finite element codes. Hence the global system of equations is generally solved using, a decomposition procedure with pivoting. The research reported to-date for the interface element includes the one-dimensional line interface element and two-dimensional surface interface element. Several large-scale simulations, including geometrically nonlinear problems, have been reported using the one-dimensional interface element technology; however, only limited applications are available for the surface interface element. In the applications reported to-date, the geometry of the interfaced domains exactly match each other even though the spatial discretization within each domain may be different. As such, the spatial modeling of each domain, the interface elements and the assembled system is still laborious. The present research is focused on developing a rapid modeling procedure based on a parametric interface representation of independently defined subdomains which are also independently discretized.
NASA Astrophysics Data System (ADS)
Sakai, Hirotaka; Urakawa, Fumihiro; Aikawa, Akira; Namura, Akira
The vibration of concrete sleepers is an important factor engendering track deterioration. In this paper, we created a three-dimensional finite element model to reproduce a prestressed concrete (PC) sleeper in detail, expressing influence of ballast layers with a 3D spring series and dampers to reproduce their vibration and dynamic characteristics. Determination of these parameters bases on the experimental modal analysis using an impact excitation technique for PC sleepers by adjusting the accelerance between the analytical results and experimental results. Furthermore, we compared the difference of these characteristics between normal sleepers and those with some structural modifications. Analytical results clarified that such means as sleeper width extension and increased sleeper thickness will influence the reduction of ballasted track vibration as improvements of PC sleepers.
NASA Astrophysics Data System (ADS)
Giasin, Khaled; Ayvar-Soberanis, Sabino; French, Toby; Phadnis, Vaibhav
2016-07-01
Machining Glass fibre aluminium reinforced epoxy (GLARE) is cumbersome due to distinctively different mechanical and thermal properties of its constituents, which makes it challenging to achieve damage-free holes with the acceptable surface quality. The proposed work focuses on the study of the machinability of thin (~2.5 mm) GLARE laminate. Drilling trials were conducted to analyse the effect of feed rate and spindle speed on the cutting forces and hole quality. The resulting hole quality metrics (surface roughness, hole size, circularity error, burr formation and delamination) were assessed using surface profilometry and optical scanning techniques. A three dimensional (3D) finite-element (FE) model of drilling GLARE laminate was also developed using ABAQUS/Explicit to help understand the mechanism of drilling GLARE. The homogenised ply-level response of GLARE laminate was considered in the FE model to predict cutting forces in the drilling process.
NASA Technical Reports Server (NTRS)
Nakazawa, Shohei
1989-01-01
The internal structure is discussed of the MHOST finite element program designed for 3-D inelastic analysis of gas turbine hot section components. The computer code is the first implementation of the mixed iterative solution strategy for improved efficiency and accuracy over the conventional finite element method. The control structure of the program is covered along with the data storage scheme and the memory allocation procedure and the file handling facilities including the read and/or write sequences.
NASA Astrophysics Data System (ADS)
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.
A 3D, finite element model for baroclinic circulation on the Vancouver Island continental shelf
Walters, R.A.; Foreman, M.G.G.
1992-01-01
This paper describes the development and application of a 3-dimensional model of the barotropic and baroclinic circulation on the continental shelf west of Vancouver Island, Canada. A previous study with a 2D barotropic model and field data revealed that several tidal constituents have a significant baroclinic component (the K1 in particular). Thus we embarked on another study with a 3D model to study the baroclinic effects on the residual and several selected tidal constituents. The 3D model uses a harmonic expansion in time and a finite element discretization in space. All nonlinear terms are retained, including quadratic bottom stress, advection and wave transport (continuity nonlinearity). The equations are solved as a global and a local problem, where the global problem is the solution of the wave equation formulation of the shallow water equations, and the local problem is the solution of the momentum equation for the vertical velocity profile. These equations are coupled to the advection-diffusion equation for density so that density gradient forcing is included in the momentum equations. However, the study presented here describes diagnostic calculations for the baroclinic residual circulation only. The model is sufficiently efficient that it encourages sensitivity testing with a large number of model runs. In this sense, the model is akin to an extension of analytical solutions to the domain of irregular geometry and bottom topography where this parameter space can be explored in some detail. In particular, the consequences of the sigma coordinate system used by the model are explored. Test cases using an idealized representation of the continental shelf, shelf break and shelf slope, lead to an estimation of the velocity errors caused by interpolation errors inherent in the sigma coordinate system. On the basis of these estimates, the computational grid used in the 2D model is found to have inadequate resolution. Thus a new grid is generated with increased
Robust and scalable 3-D geo-electromagnetic modelling approach using the finite element method
NASA Astrophysics Data System (ADS)
Grayver, Alexander V.; Bürg, Markus
2014-07-01
We present a robust and scalable solver for time-harmonic Maxwell's equations for problems with large conductivity contrasts, wide range of frequencies, stretched grids and locally refined meshes. The solver is part of the fully distributed adaptive 3-D electromagnetic modelling scheme which employs the finite element method and unstructured non-conforming hexahedral meshes for spatial discretization using the open-source software deal.II. We use the complex-valued electric field formulation and split it into two real-valued equations for which we utilize an optimal block-diagonal pre-conditioner. Application of this pre-conditioner requires the solution of two smaller real-valued symmetric problems. We solve them by using either a direct solver or the conjugate gradient method pre-conditioned with the recently introduced auxiliary space technique. The auxiliary space pre-conditioner reformulates the original problem in form of several simpler ones, which are then solved using highly efficient algebraic multigrid methods. In this paper, we consider the magnetotelluric case and verify our numerical scheme by using COMMEMI 3-D models. Afterwards, we run a series of numerical experiments and demonstrate that the solver converges in a small number of iterations for a wide frequency range and variable problem sizes. The number of iterations is independent of the problem size, but exhibits a mild dependency on frequency. To test the stability of the method on locally refined meshes, we have implemented a residual-based a posteriori error estimator and compared it with uniform mesh refinement for problems up to 200 million unknowns. We test the scalability of the most time consuming parts of our code and show that they fulfill the strong scaling assumption as long as each MPI process possesses enough degrees of freedom to alleviate communication overburden. Finally, we refer back to a direct solver-based pre-conditioner and analyse its complexity in time. The results show
3D Finite Element Analysis of Spider Non-isothermal Forging Process
NASA Astrophysics Data System (ADS)
Niu, Ling; Wei, Wei; Wei, Kun Xia; Alexandrov, Igor V.; Hu, Jing
2016-05-01
The differences of effective stress, effective strain, velocity field, and the load-time curves between the spider isothermal and non-isothermal forging processes are investigated by making full use of 3D FEA, and verified by the production experiment of spider forging. Effective stress is mainly concentrated on the pin, and becomes lower closer to the front of the pin. The maximum effective strain in the non-isothermal forging is lower than that in the isothermal. The great majority of strain in the non-isothermal forging process is 1.76, which is larger than the strain of 1.31 in the isothermal forging. The maximum load required in the isothermal forging is higher than that in the non-isothermal. The maximum experimental load and deformation temperature in the spider production are in good agreement with those in the non-isothermal FEA. The results indicate that the non-isothermal 3D FEA results can guide the design of the spider forging process.
3D Finite Element Analysis of Spider Non-isothermal Forging Process
NASA Astrophysics Data System (ADS)
Niu, Ling; Wei, Wei; Wei, Kun Xia; Alexandrov, Igor V.; Hu, Jing
2016-06-01
The differences of effective stress, effective strain, velocity field, and the load-time curves between the spider isothermal and non-isothermal forging processes are investigated by making full use of 3D FEA, and verified by the production experiment of spider forging. Effective stress is mainly concentrated on the pin, and becomes lower closer to the front of the pin. The maximum effective strain in the non-isothermal forging is lower than that in the isothermal. The great majority of strain in the non-isothermal forging process is 1.76, which is larger than the strain of 1.31 in the isothermal forging. The maximum load required in the isothermal forging is higher than that in the non-isothermal. The maximum experimental load and deformation temperature in the spider production are in good agreement with those in the non-isothermal FEA. The results indicate that the non-isothermal 3D FEA results can guide the design of the spider forging process.
SIMULATIONS OF 2D AND 3D THERMOCAPILLARY FLOWS BY A LEAST-SQUARES FINITE ELEMENT METHOD. (R825200)
Numerical results for time-dependent 2D and 3D thermocapillary flows are presented in this work. The numerical algorithm is based on the Crank-Nicolson scheme for time integration, Newton's method for linearization, and a least-squares finite element method, together with a matri...
Eraslan, Oğuz; Inan, Ozgür
2010-08-01
The biomechanical behavior of implant thread plays an important role on stresses at implant-bone interface. Information about the effect of different thread profiles upon the bone stresses is limited. The purpose of this study was to evaluate the effects of different implant thread designs on stress distribution characteristics at supporting structures. In this study, three-dimensional (3D) finite element (FE) stress-analysis method was used. Four types of 3D mathematical models simulating four different thread-form configurations for a solid screw implant was prepared with supporting bone structure. V-thread (1), buttress (2), reverse buttress (3), and square thread designs were simulated. A 100-N static axial occlusal load was applied to occlusal surface of abutment to calculate the stress distributions. Solidworks/Cosmosworks structural analysis programs were used for FE modeling/analysis. The analysis of the von Mises stress values revealed that maximum stress concentrations were located at loading areas of implant abutments and cervical cortical bone regions for all models. Stress concentration at cortical bone (18.3 MPa) was higher than spongious bone (13.3 MPa), and concentration of first thread (18 MPa) was higher than other threads (13.3 MPa). It was seen that, while the von Mises stress distribution patterns at different implant thread models were similar, the concentration of compressive stresses were different. The present study showed that the use of different thread form designs did not affect the von Mises concentration at supporting bone structure. However, the compressive stress concentrations differ by various thread profiles. PMID:19543925
NASA Astrophysics Data System (ADS)
Mulder, W. A.; Zhebel, E.; Minisini, S.
2014-02-01
We analyse the time-stepping stability for the 3-D acoustic wave equation, discretized on tetrahedral meshes. Two types of methods are considered: mass-lumped continuous finite elements and the symmetric interior-penalty discontinuous Galerkin method. Combining the spatial discretization with the leap-frog time-stepping scheme, which is second-order accurate and conditionally stable, leads to a fully explicit scheme. We provide estimates of its stability limit for simple cases, namely, the reference element with Neumann boundary conditions, its distorted version of arbitrary shape, the unit cube that can be partitioned into six tetrahedra with periodic boundary conditions and its distortions. The Courant-Friedrichs-Lewy stability limit contains an element diameter for which we considered different options. The one based on the sum of the eigenvalues of the spatial operator for the first-degree mass-lumped element gives the best results. It resembles the diameter of the inscribed sphere but is slightly easier to compute. The stability estimates show that the mass-lumped continuous and the discontinuous Galerkin finite elements of degree 2 have comparable stability conditions, whereas the mass-lumped elements of degree one and three allow for larger time steps.
NASA Astrophysics Data System (ADS)
Koldan, Jelena; Puzyrev, Vladimir; de la Puente, Josep; Houzeaux, Guillaume; Cela, José María
2014-06-01
We present an elaborate preconditioning scheme for Krylov subspace methods which has been developed to improve the performance and reduce the execution time of parallel node-based finite-element (FE) solvers for 3-D electromagnetic (EM) numerical modelling in exploration geophysics. This new preconditioner is based on algebraic multigrid (AMG) that uses different basic relaxation methods, such as Jacobi, symmetric successive over-relaxation (SSOR) and Gauss-Seidel, as smoothers and the wave front algorithm to create groups, which are used for a coarse-level generation. We have implemented and tested this new preconditioner within our parallel nodal FE solver for 3-D forward problems in EM induction geophysics. We have performed series of experiments for several models with different conductivity structures and characteristics to test the performance of our AMG preconditioning technique when combined with biconjugate gradient stabilized method. The results have shown that, the more challenging the problem is in terms of conductivity contrasts, ratio between the sizes of grid elements and/or frequency, the more benefit is obtained by using this preconditioner. Compared to other preconditioning schemes, such as diagonal, SSOR and truncated approximate inverse, the AMG preconditioner greatly improves the convergence of the iterative solver for all tested models. Also, when it comes to cases in which other preconditioners succeed to converge to a desired precision, AMG is able to considerably reduce the total execution time of the forward-problem code-up to an order of magnitude. Furthermore, the tests have confirmed that our AMG scheme ensures grid-independent rate of convergence, as well as improvement in convergence regardless of how big local mesh refinements are. In addition, AMG is designed to be a black-box preconditioner, which makes it easy to use and combine with different iterative methods. Finally, it has proved to be very practical and efficient in the
Exponential Finite-Difference Technique
NASA Technical Reports Server (NTRS)
Handschuh, Robert F.
1989-01-01
Report discusses use of explicit exponential finite-difference technique to solve various diffusion-type partial differential equations. Study extends technique to transient-heat-transfer problems in one dimensional cylindrical coordinates and two and three dimensional Cartesian coordinates and to some nonlinear problems in one or two Cartesian coordinates.
A Simulation of crustal deformation around sourthwest Japan using 3D Finite Element Method
NASA Astrophysics Data System (ADS)
Oma, T.; Ito, T.; Sasajima, R.
2015-12-01
In southwest Japan, the Philippine Sea plate is subducting beneath the Amurian plate at the Nankai Trough. Megathrust earthquakes have been occurred with recurrence intervals of about 100-150 years. Previous studies have estimated co-seismic slip distribution at the 1944 Tokankai and the 1946 Nankai earthquakes and interplate plate coupling along the Nankai Trough. Many of previous studies employed a homogeneous elastic half space or elastic and viscoelastic layers structure. However, these assumptions as mentioned above are inadequate, since inhomogeneous structure is exceled in the real earth result from subducting plate. Therefore, in order to estimate the effect of inhomogeneous structure on the crustal deformation, we calculate crustal deformation due to Megathrust earthquake using 3-dimensional Finite Element Method (FEM). We use FEM software PyLith v2.1. In this study, we construct a finite element mesh with the region of 3000km(SW) × 2300km(NS) × 400km(depth) cover Japanese Islands, using Cubit 13.0. This mesh is considered topography, the Philippine Sea plate, the Pacific plate, Moho discontinuity, and curvature of the earth. In order to examine differences of surface displacement between inhomogeneous and homogeneous structures, we use co-seismic slip distribution of the 1944 and 1946 earthquakes estimated by Sagiya and Thatcher (1999). In result, surface elastic response under inhomogeneous structure becomes 30% larger than it's homogeneous structure at the Muroto cape. This difference indicates that co-seismic slip or plate coupling distribution estimated from Green's function under an assumption of homogeneous structure is overestimated. Then, we calculate viscoelastic response assuming Maxwell rheology model and viscosity as 1×1019. As a result, predicted horizontal velocity of viscoelastic response due to the events corresponds to 10 % of observed present deformation. It suggest that spatial pattern of plate coupling might be change when we
2013-01-01
PURPOSE This study was accomplished to assess the biomechanical state of different retaining methods of bar implant-overdenture. MATERIALS AND METHODS Two 3D finite element models were designed. The first model included implant overdenture retained by Hader-clip attachment, while the second model included two extracoronal resilient attachment (ERA) studs added distally to Hader splint bar. A non-linear frictional contact type was assumed between overdentures and mucosa to represent sliding and rotational movements among different attachment components. A 200 N was applied at the molar region unilaterally and perpendicular to the occlusal plane. Additionally, the mandible was restrained at their ramus ends. The maximum equivalent stress and strain (von Mises) were recorded and analyzed at the bone-implant interface level. RESULTS The values of von Mises stress and strain of the first model at bone-implant interface were higher than their counterparts of the second model. Stress concentration and high value of strain were recognized surrounding implant of the unloaded side in both models. CONCLUSION There were different patterns of stress-strain distribution at bone-implant interface between the studied attachment designs. Hader bar-clip attachment showed better biomechanical behavior than adding ERA studs distal to hader bar. PMID:24049576
Hsu, Christina M. L.; Palmeri, Mark L.; Segars, W. Paul; Veress, Alexander I.; Dobbins, James T.
2011-01-01
Purpose: The authors previously introduced a methodology to generate a realistic three-dimensional (3D), high-resolution, computer-simulated breast phantom based on empirical data. One of the key components of such a phantom is that it provides a means to produce a realistic simulation of clinical breast compression. In the current study, they have evaluated a finite element (FE) model of compression and have demonstrated the effect of a variety of mechanical properties on the model using a dense mesh generated from empirical breast data. While several groups have demonstrated an effective compression simulation with lower density finite element meshes, the presented study offers a mesh density that is able to model the morphology of the inner breast structures more realistically than lower density meshes. This approach may prove beneficial for multimodality breast imaging research, since it provides a high level of anatomical detail throughout the simulation study. Methods: In this paper, the authors describe methods to improve the high-resolution performance of a FE compression model. In order to create the compressible breast phantom, dedicated breast CT data was segmented and a mesh was generated with 4-noded tetrahedral elements. Using an explicit FE solver to simulate breast compression, several properties were analyzed to evaluate their effect on the compression model including: mesh density, element type, density, and stiffness of various tissue types, friction between the skin and the compression plates, and breast density. Following compression, a simulated projection was generated to demonstrate the ability of the compressible breast phantom to produce realistic simulated mammographic images. Results: Small alterations in the properties of the breast model can change the final distribution of the tissue under compression by more than 1 cm; which ultimately results in different representations of the breast model in the simulated images. The model
NASA Astrophysics Data System (ADS)
Zhang, Zhi-Qian; Liu, G. R.; Khoo, Boo Cheong
2013-02-01
A three-dimensional immersed smoothed finite element method (3D IS-FEM) using four-node tetrahedral element is proposed to solve 3D fluid-structure interaction (FSI) problems. The 3D IS-FEM is able to determine accurately the physical deformation of the nonlinear solids placed within the incompressible viscous fluid governed by Navier-Stokes equations. The method employs the semi-implicit characteristic-based split scheme to solve the fluid flows and smoothed finite element methods to calculate the transient dynamics responses of the nonlinear solids based on explicit time integration. To impose the FSI conditions, a novel, effective and sufficiently general technique via simple linear interpolation is presented based on Lagrangian fictitious fluid meshes coinciding with the moving and deforming solid meshes. In the comparisons to the referenced works including experiments, it is clear that the proposed 3D IS-FEM ensures stability of the scheme with the second order spatial convergence property; and the IS-FEM is fairly independent of a wide range of mesh size ratio.
NASA Astrophysics Data System (ADS)
Wendling, A.; Daniel, J. L.; Hivet, G.; Vidal-Sallé, E.; Boisse, P.
2015-12-01
Numerical simulation is a powerful tool to predict the mechanical behavior and the feasibility of composite parts. Among the available numerical approaches, as far as woven reinforced composites are concerned, 3D finite element simulation at the mesoscopic scale leads to a good compromise between realism and complexity. At this scale, the fibrous reinforcement is modeled by an interlacement of yarns assumed to be homogeneous that have to be accurately represented. Among the numerous issues induced by these simulations, the first one consists in providing a representative meshed geometrical model of the unit cell at the mesoscopic scale. The second one consists in enabling a fast data input in the finite element software (contacts definition, boundary conditions, elements reorientation, etc.) so as to obtain results within reasonable time. Based on parameterized 3D CAD modeling tool of unit-cells of dry fabrics already developed, this paper presents an efficient strategy which permits an automated meshing of the models with 3D hexahedral elements and to accelerate of several orders of magnitude the simulation data input. Finally, the overall modeling strategy is illustrated by examples of finite element simulation of the mechanical behavior of fabrics.
A feasibility study of a 3-D finite element solution scheme for aeroengine duct acoustics
NASA Technical Reports Server (NTRS)
Abrahamson, A. L.
1980-01-01
The advantage from development of a 3-D model of aeroengine duct acoustics is the ability to analyze axial and circumferential liner segmentation simultaneously. The feasibility of a 3-D duct acoustics model was investigated using Galerkin or least squares element formulations combined with Gaussian elimination, successive over-relaxation, or conjugate gradient solution algorithms on conventional scalar computers and on a vector machine. A least squares element formulation combined with a conjugate gradient solver on a CDC Star vector computer initially appeared to have great promise, but severe difficulties were encountered with matrix ill-conditioning. These difficulties in conditioning rendered this technique impractical for realistic problems.
Hallquist, J.O.
1981-01-01
A user's manual is provided for NIKE3D, a fully implicit three-dimensional finite element code for analyzing the large deformation static and dynamic response of inelastic solids. A contact-impact algorithm permits gaps and sliding along material interfaces. By a specialization of this algorithm, such interfaces can be rigidly tied to admit variable zoning without the need of transition regions. Spatial discretization is achieved by the use of 8-node constant pressure solid elements. Bandwidth minimization is optional. Post-processors for NIKE3D include GRAPE for plotting deformed shapes and stress contours and DYNAP for plotting time histories.
Advances in 3D-Printed Pediatric Prostheses for Upper Extremity Differences.
Tanaka, Kara S; Lightdale-Miric, Nina
2016-08-01
➤The prohibitive cost of cutting-edge prostheses prevents many children with a limb difference from obtaining them; however, new developments in 3-dimensional (3D) printing have the potential to increase the accessibility, customization, and procurement of such devices.➤Children with upper limb differences are ideal candidates for currently available 3D-printed devices because they quickly damage and outgrow prostheses, and the low cost of 3D printing makes repairs and upgrades substantially more affordable.➤Physicians and medical practitioners should become familiar with the possibilities of 3D-printed devices in order to determine the benefits and utility for their patients. PMID:27489324
Study of Multi Pass Equal Channel Angular Pressing Using 3D Finite Element Analysis
NASA Astrophysics Data System (ADS)
Setia, Rajat; Sharma, Rahul Swarup; Sharma, Shanti Swarup; Raj, K. Hans
2011-01-01
Equal Channel Angular Pressing (ECAP) has emerged as most prominent Severe Plastic Deformation (SPD) technique used to produce an ultrafine grained (UFG) structure in metals in order to improve their mechanical and physical properties. In this work Finite Element modeling of ECAP is attempted in FORGE 2007 environment. Four passes of the ECAP process of 10mm square shaped AL 6061 billet were carried out for routes A, BA and C for different channel angles and values of coefficient of friction to investigate their influence on the billet. The models were developed assuming a range of friction conditions at the billet-die contact region considering eight distinct friction coefficient (μ) values of 0.0, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35 and 0.40, respectively. The simulations are carried out using three distinct situations of die channel angles (Φ), 90°, 105°, and 120° respectively. Route `BA' emerged as a better method among the three routes studied and 90° channel angle appeared to be optimal in terms of producing high equivalent strain.
3-d finite element model development for biomechanics: a software demonstration
Hollerbach, K.; Hollister, A.M.; Ashby, E.
1997-03-01
Finite element analysis is becoming an increasingly important part of biomechanics and orthopedic research, as computational resources become more powerful, and data handling algorithms become more sophisticated. Until recently, tools with sufficient power did not exist or were not accessible to adequately model complicated, three-dimensional, nonlinear biomechanical systems. In the past, finite element analyses in biomechanics have often been limited to two-dimensional approaches, linear analyses, or simulations of single tissue types. Today, we have the resources to model fully three-dimensional, nonlinear, multi-tissue, and even multi-joint systems. The authors will present the process of developing these kinds of finite element models, using human hand and knee examples, and will demonstrate their software tools.
Visual Short-Term Memory Benefit for Objects on Different 3-D Surfaces
ERIC Educational Resources Information Center
Xu, Yaoda; Nakayama, Ken
2007-01-01
Visual short-term memory (VSTM) plays an important role in visual cognition. Although objects are located on different 3-dimensional (3-D) surfaces in the real world, how VSTM capacity may be influenced by the presence of multiple 3-D surfaces has never been examined. By manipulating binocular disparities of visual displays, the authors found that…
Lin, Jie; Zheng, Zhiqiang; Shinya, Akikazu; Matinlinna, Jukka Pekka; Botelho, Michael George; Shinya, Akiyoshi
2015-09-01
The purpose of this in vitro study was to compare the stress distribution and natural frequency of different shape and thickness retainer designs for maxillary posterior resin-bonded prostheses using finite element (FE) method. A 3D FE model of a three unit posterior resin-bonded prosthesis analysis model was generated. Three different shaped retainer designs, viz. C-shaped (three axial surface wraparounds), D-shaped (three axial surface wraparounds with central groove) and O-shaped (360° wraparounds), and three different thicknesses, viz., 0.4, 0.8, and 1.2 mm, resin-bonded prostheses were used in this study. The resin-bonded prosthesis analysis model was imported into an FE analysis software (ANSYS 10.0, ANSYS, USA) and attribution of material properties. The nodes at the bottom surface of the roots were assigned fixed zero displacement in the three spatial dimensions. A simulated angle of 45° loading of a 100 N force was applied to the node of the pontic lingual cusp surface. The stress distributions and corresponding natural frequencies were analyzed and resolved. The C-shaped retainer for 0.4 mm thickness recorded the greatest von Mises stresses of 71.4 MPa for all three groups. C-shaped, D-shaped and O-shaped retainer presented natural frequencies 3,988, 7,754, and 10,494 Hz, respectively. D-shaped retainer and O-shaped retainer increased natural frequencies and structural rigidity over the traditional C-shaped retainer. The maximum von Mises stresses values of the remaining tooth and prosthesis decreased with greater retainer thickness. D-shaped retainer and O-shaped retainer increased natural frequencies and structural rigidity over the traditional C-shaped retainer. PMID:25200313
Edge-based finite elements and vector ABCs applied to 3D scattering
NASA Technical Reports Server (NTRS)
Chatterjee, A.; Jin, J. M.; Volakis, John L.
1992-01-01
An edge based finite element formulation with vector absorbing boundary conditions is presented for scattering by composite structures having boundaries satisfying impedance and/or transition conditions. Remarkably accurate results are obtained by placing the mesh a small fraction of a wavelength away from the scatterer.
Väänänen, Sami P; Grassi, Lorenzo; Flivik, Gunnar; Jurvelin, Jukka S; Isaksson, Hanna
2015-08-01
Areal bone mineral density (aBMD), as measured by dual-energy X-ray absorptiometry (DXA), predicts hip fracture risk only moderately. Simulation of bone mechanics based on DXA imaging of the proximal femur, may help to improve the prediction accuracy. Therefore, we collected three (1-3) image sets, including CT images and DXA images of 34 proximal cadaver femurs (set 1, including 30 males, 4 females), 35 clinical patient CT images of the hip (set 2, including 27 males, 8 females) and both CT and DXA images of clinical patients (set 3, including 12 female patients). All CT images were segmented manually and landmarks were placed on both femurs and pelvises. Two separate statistical appearance models (SAMs) were built using the CT images of the femurs and pelvises in sets 1 and 2, respectively. The 3D shape of the femur was reconstructed from the DXA image by matching the SAMs with the DXA images. The orientation and modes of variation of the SAMs were adjusted to minimize the sum of the absolute differences between the projection of the SAMs and a DXA image. The mesh quality and the location of the SAMs with respect to the manually placed control points on the DXA image were used as additional constraints. Then, finite element (FE) models were built from the reconstructed shapes. Mean point-to-surface distance between the reconstructed shape and CT image was 1.0 mm for cadaver femurs in set 1 (leave-one-out test) and 1.4 mm for clinical subjects in set 3. The reconstructed volumetric BMD showed a mean absolute difference of 140 and 185 mg/cm(3) for set 1 and set 3 respectively. The generation of the SAM and the limitation of using only one 2D image were found to be the most significant sources of errors in the shape reconstruction. The noise in the DXA images had only small effect on the accuracy of the shape reconstruction. DXA-based FE simulation was able to explain 85% of the CT-predicted strength of the femur in stance loading. The present method can be used to
NASA Astrophysics Data System (ADS)
Kis, M.; Detzky, G.; Koppán, A.
2012-04-01
phenomenon in general. Authors calculated the deformations of a simple-geometry 3D cavity, which is caused by variable gravity loads. Dependence of the cavity effect on changing of distinct elastic properties in categorized models has been investigated. Authors introduced qualifying parameter fields calculated using the results of the FE modelling (nodal displacements as a model answer for the gravity load), in order to characterize the effect. Modelling results can be used as an estimation not only for the absolute cavity effect rate of the intended arrangement, furthermore the sensitivity of the given system against a particular geometric property. As an application example finite element modelling were carried out in order to estimate the influence of the complicated cavity system surrounding the "Budapest-Matyashegy" Gravity and Geodynamical Observatory of the Eotvos Lorand Geophysical Institute of Hungary.
The Complex-Step-Finite-Difference method
NASA Astrophysics Data System (ADS)
Abreu, Rafael; Stich, Daniel; Morales, Jose
2015-07-01
We introduce the Complex-Step-Finite-Difference method (CSFDM) as a generalization of the well-known Finite-Difference method (FDM) for solving the acoustic and elastic wave equations. We have found a direct relationship between modelling the second-order wave equation by the FDM and the first-order wave equation by the CSFDM in 1-D, 2-D and 3-D acoustic media. We present the numerical methodology in order to apply the introduced CSFDM and show an example for wave propagation in simple homogeneous and heterogeneous models. The CSFDM may be implemented as an extension into pre-existing numerical techniques in order to obtain fourth- or sixth-order accurate results with compact three time-level stencils. We compare advantages of imposing various types of initial motion conditions of the CSFDM and demonstrate its higher-order accuracy under the same computational cost and dispersion-dissipation properties. The introduced method can be naturally extended to solve different partial differential equations arising in other fields of science and engineering.
Simulating hydroplaning of submarine landslides by quasi 3D depth averaged finite element method
NASA Astrophysics Data System (ADS)
De Blasio, Fabio; Battista Crosta, Giovanni
2014-05-01
G.B. Crosta, H. J. Chen, and F.V. De Blasio Dept. Of Earth and Environmental Sciences, Università degli Studi di Milano Bicocca, Milano, Italy Klohn Crippen Berger, Calgary, Canada Subaqueous debris flows/submarine landslides, both in the open ocean as well as in fresh waters, exhibit extremely high mobility, quantified by a ratio between vertical to horizontal displacement of the order 0.01 or even much less. It is possible to simulate subaqueous debris flows with small-scale experiments along a flume or a pool using a cohesive mixture of clay and sand. The results have shown a strong enhancement of runout and velocity compared to the case in which the same debris flow travels without water, and have indicated hydroplaning as a possible explanation (Mohrig et al. 1998). Hydroplaning is started when the snout of the debris flow travels sufficiently fast. This generates lift forces on the front of the debris flow exceeding the self-weight of the sediment, which so begins to travel detached from the bed, literally hovering instead of flowing. Clearly, the resistance to flow plummets because drag stress against water is much smaller than the shear strength of the material. The consequence is a dramatic increase of the debris flow speed and runout. Does the process occur also for subaqueous landslides and debris flows in the ocean, something twelve orders of magnitude larger than the experimental ones? Obviously, no experiment will ever be capable to replicate this size, one needs to rely on numerical simulations. Results extending a depth-integrated numerical model for debris flows (Imran et al., 2001) indicate that hydroplaning is possible (De Blasio et al., 2004), but more should be done especially with alternative numerical methodologies. In this work, finite element methods are used to simulate hydroplaning using the code MADflow (Chen, 2014) adopting a depth averaged solution. We ran some simulations on the small scale of the laboratory experiments, and secondly
Simulation of 3D tumor cell growth using nonlinear finite element method.
Dong, Shoubing; Yan, Yannan; Tang, Liqun; Meng, Junping; Jiang, Yi
2016-06-01
We propose a novel parallel computing framework for a nonlinear finite element method (FEM)-based cell model and apply it to simulate avascular tumor growth. We derive computation formulas to simplify the simulation and design the basic algorithms. With the increment of the proliferation generations of tumor cells, the FEM elements may become larger and more distorted. Then, we describe a remesh and refinement processing of the distorted or over large finite elements and the parallel implementation based on Message Passing Interface to improve the accuracy and efficiency of the simulation. We demonstrate the feasibility and effectiveness of the FEM model and the parallelization methods in simulations of early tumor growth. PMID:26213205
A 3-D Finite-Volume Non-hydrostatic Icosahedral Model (NIM)
NASA Astrophysics Data System (ADS)
Lee, Jin
2014-05-01
The Nonhydrostatic Icosahedral Model (NIM) formulates the latest numerical innovation of the three-dimensional finite-volume control volume on the quasi-uniform icosahedral grid suitable for ultra-high resolution simulations. NIM's modeling goal is to improve numerical accuracy for weather and climate simulations as well as to utilize the state-of-art computing architecture such as massive parallel CPUs and GPUs to deliver routine high-resolution forecasts in timely manner. NIM dynamic corel innovations include: * A local coordinate system remapped spherical surface to plane for numerical accuracy (Lee and MacDonald, 2009), * Grid points in a table-driven horizontal loop that allow any horizontal point sequence (A.E. MacDonald, et al., 2010), * Flux-Corrected Transport formulated on finite-volume operators to maintain conservative positive definite transport (J.-L, Lee, ET. Al., 2010), *Icosahedral grid optimization (Wang and Lee, 2011), * All differentials evaluated as three-dimensional finite-volume integrals around the control volume. The three-dimensional finite-volume solver in NIM is designed to improve pressure gradient calculation and orographic precipitation over complex terrain. NIM dynamical core has been successfully verified with various non-hydrostatic benchmark test cases such as internal gravity wave, and mountain waves in Dynamical Cores Model Inter-comparisons Projects (DCMIP). Physical parameterizations suitable for NWP are incorporated into NIM dynamical core and successfully tested with multimonth aqua-planet simulations. Recently, NIM has started real data simulations using GFS initial conditions. Results from the idealized tests as well as real-data simulations will be shown in the conference.
Binary 3-D Markov Chain Random Fields: Finite-size Scaling Analysis of Percolation Properties
NASA Astrophysics Data System (ADS)
Harter, T.
2004-12-01
Percolation phenomena in random media have been extensively studied in a wide variety of fields in physics, chemistry, engineering, bio-, earth-, and environmental sciences. Most work has focused on uncorrelated random fields. The critical behavior in media with short-range correlations is thought to be identical to that in uncorrelated systems. However, the percolation threshold, pc, which is 0.3116 in uncorrelated media, has been observed to vary with the correlation scale and also with the random field type. Here, we present percolation properties and finite-size scaling effects in three-dimensional binary cubic lattices represented by correlated Markov-chain random fields and compare them to those in sequential Gaussian and sequential indicator random fields. We find that the computed percolation threshold in correlated random fields is significantly lower than in the uncorrelated lattice and decreases with increasing correlation scale. The rate of decrease rapidly flattens out for correlation lengths larger than 2-3 grid-blocks. At correlation scales of 5-6 grid blocks, pc is found to be 0.126 for the Markov chain random fields and slightly higher for sequential Gaussian and indicator random fields. The universal scaling constants for mean cluster size, backbone fraction, and connectivity are found to be consistent with results on uncorrelated lattices. For numerical studies, it is critical to understand finite-size effects on the percolation and associated phase connectivity properties of lattices. We present detailed statistical results on the percolation properties in finite sized lattice and their dependence on correlation scale. We show that appropriate grid resolution and choice of simulation boundaries is critical to properly simulate correlated natural geologic systems, which may display significant finite-size effects.
A 3D finite element simulation model for TBM tunnelling in soft ground
NASA Astrophysics Data System (ADS)
Kasper, Thomas; Meschke, Günther
2004-12-01
A three-dimensional finite element simulation model for shield-driven tunnel excavation is presented. The model takes into account all relevant components of the construction process (the soil and the ground water, the tunnel boring machine with frictional contact to the soil, the hydraulic jacks, the tunnel lining and the tail void grouting). The paper gives a detailed description of the model components and the stepwise procedure to simulate the construction process. The soil and the grout material are modelled as saturated porous media using a two-field finite element formulation. This allows to take into account the groundwater, the grouting pressure and the fluid interaction between the soil and slurry at the cutting face and between the soil and grout around the tail void. A Cam-Clay plasticity model is used to describe the material behaviour of cohesive soils. The cementitious grouting material in the tail void is modelled as an ageing elastic material with time-dependent stiffness and permeability. To allow for an automated computation of arbitrarily long and also curvilinear driving paths with suitable finite element meshes, the simulation procedure has been fully automated. The simulation of a tunnel advance in soft cohesive soil below the ground water table is presented and the results are compared with measurements taken from the literature. Copyright
NASA Astrophysics Data System (ADS)
Uhrig, Matthias P.; Kim, Jin-Yeon; Jacobs, Laurence J.
2016-02-01
This research presents a 3D numerical finite element (FE) model which, previously developed, precisely simulates non-contact, air-coupled measurements of nonlinear Rayleigh wave propagation. The commercial FE-solver ABAQUS is used to perform the simulations. First, frequency dependent pressure wave attenuation is investigated numerically to reconstruct the sound pressure distribution along the active surface of the non-contact receiver. Second, constitutive law and excitation source properties are optimized to match nonlinear ultrasonic experimental data. Finally, the FE-model data are fit with analytical solutions showing a good agreement and thus, indicating the significance of the study performed.
Boda-Heggemann, Judit . E-mail: judit.boda-heggemann@radonk.ma.uni-heidelberg.de; Walter, Cornelia; Rahn, Angelika; Wertz, Hansjoerg; Loeb, Iris; Lohr, Frank; Wenz, Frederik
2006-12-01
Purpose: The repositioning accuracy of mask-based fixation systems has been assessed with two-dimensional/two-dimensional or two-dimensional/three-dimensional (3D) matching. We analyzed the accuracy of commercially available head mask systems, using true 3D/3D matching, with X-ray volume imaging and cone-beam CT. Methods and Materials: Twenty-one patients receiving radiotherapy (intracranial/head-and-neck tumors) were evaluated (14 patients with rigid and 7 with thermoplastic masks). X-ray volume imaging was analyzed online and offline separately for the skull and neck regions. Translation/rotation errors of the target isocenter were analyzed. Four patients were treated to neck sites. For these patients, repositioning was aided by additional body tattoos. A separate analysis of the setup error on the basis of the registration of the cervical vertebra was performed. The residual error after correction and intrafractional motility were calculated. Results: The mean length of the displacement vector for rigid masks was 0.312 {+-} 0.152 cm (intracranial) and 0.586 {+-} 0.294 cm (neck). For the thermoplastic masks, the value was 0.472 {+-} 0.174 cm (intracranial) and 0.726 {+-} 0.445 cm (neck). Rigid masks with body tattoos had a displacement vector length in the neck region of 0.35 {+-} 0.197 cm. The intracranial residual error and intrafractional motility after X-ray volume imaging correction for rigid masks was 0.188 {+-} 0.074 cm, and was 0.134 {+-} 0.14 cm for thermoplastic masks. Conclusions: The results of our study have demonstrated that rigid masks have a high intracranial repositioning accuracy per se. Given the small residual error and intrafractional movement, thermoplastic masks may also be used for high-precision treatments when combined with cone-beam CT. The neck region repositioning accuracy was worse than the intracranial accuracy in both cases. However, body tattoos and image guidance improved the accuracy. Finally, the combination of both mask
NASA Astrophysics Data System (ADS)
Pereira, J. P.; Duarte, C. A.; Jiao, X.; Guoy, D.
2009-06-01
This paper presents a study of generalized enrichment functions for 3D curved crack fronts. Two coordinate systems used in the definition of singular curved crack front enrichment functions are analyzed. In the first one, a set of Cartesian coordinate systems defined along the crack front is used. In the second case, the geometry of the crack front is approximated by a set of curvilinear coordinate systems. A description of the computation of derivatives of enrichment functions and curvilinear base vectors is presented. The coordinate systems are automatically defined using geometrical information provided by an explicit representation of the crack surface. A detailed procedure to accurately evaluate the surface normal, conormal and tangent vectors along curvilinear crack fronts in explicit crack surface representations is also presented. An accurate and robust definition of orthonormal vectors along crack fronts is crucial for the proper definition of enrichment functions. Numerical experiments illustrate the accuracy and robustness of the proposed approaches.
Image forces on 3d dislocation structures in crystals of finite volume
El-Azab, A.
1999-07-01
The present work aims at studying the image stress and image Peach-Koehler force fields for three-dimensional dislocation configurations in a single crystal of finite volume. It is shown that the image stress field is significant within the entire crystal volume, and that the image Peach-Koehler force can be of the same order of magnitude as the direct interaction force calculated from the infinite domain solution. The results demonstrate that image stress gives rise to long-range interaction forces that are important in meso-scale dynamics of dislocation structures.
Image Forces on 3-D Dislocation Structures in Crystals of Finite Volume
El-Azab, Anter ); V.V. Bulatov
1999-01-01
The present work aims at studying the image stress and image Peach-Koehler force fields for three-dimensional dislocation configurations in a single crystal of finite volume. It is shown that the image stress field is significant within the entire crystal volume, and that the image Peach-Koehler force can be of the same order of magnitude as the direct interaction force calculated from the infinite domain solution. The results demonstrate that image stress gives rise to long-range interaction forces that are important in meso-scale dynamics of dislocation structures.
A Lagrange-Galerkin hp-Finite Element Method for a 3D Nonhydrostatic Ocean Model
NASA Astrophysics Data System (ADS)
Galán del Sastre, Pedro; Bermejo, Rodolfo
2016-03-01
We introduce in this paper a Lagrange-Galerkin hp-finite element method to calculate the numerical solution of a nonhydrostatic ocean model. The Lagrange-Galerkin method yields a Stokes-like problem the solution of which is computed by a second-order rotational splitting scheme that separates the calculation of the velocity and pressure, the latter is decomposed into hydrostatic and nonhydrostatic components. We have tested the method in flows where the nonhydrostatic effects are important. The results are very encouraging.
NASA Technical Reports Server (NTRS)
Leser, Patrick E.; Hochhalter, Jacob D.; Newman, John A.; Leser, William P.; Warner, James E.; Wawrzynek, Paul A.; Yuan, Fuh-Gwo
2015-01-01
Utilizing inverse uncertainty quantification techniques, structural health monitoring can be integrated with damage progression models to form probabilistic predictions of a structure's remaining useful life. However, damage evolution in realistic structures is physically complex. Accurately representing this behavior requires high-fidelity models which are typically computationally prohibitive. In the present work, a high-fidelity finite element model is represented by a surrogate model, reducing computation times. The new approach is used with damage diagnosis data to form a probabilistic prediction of remaining useful life for a test specimen under mixed-mode conditions.
Kılıç, Emre Eibert, Thomas F.
2015-05-01
An approach combining boundary integral and finite element methods is introduced for the solution of three-dimensional inverse electromagnetic medium scattering problems. Based on the equivalence principle, unknown equivalent electric and magnetic surface current densities on a closed surface are utilized to decompose the inverse medium problem into two parts: a linear radiation problem and a nonlinear cavity problem. The first problem is formulated by a boundary integral equation, the computational burden of which is reduced by employing the multilevel fast multipole method (MLFMM). Reconstructed Cauchy data on the surface allows the utilization of the Lorentz reciprocity and the Poynting's theorems. Exploiting these theorems, the noise level and an initial guess are estimated for the cavity problem. Moreover, it is possible to determine whether the material is lossy or not. In the second problem, the estimated surface currents form inhomogeneous boundary conditions of the cavity problem. The cavity problem is formulated by the finite element technique and solved iteratively by the Gauss–Newton method to reconstruct the properties of the object. Regularization for both the first and the second problems is achieved by a Krylov subspace method. The proposed method is tested against both synthetic and experimental data and promising reconstruction results are obtained.
Finite element and finite difference methods in electromagnetic scattering
NASA Astrophysics Data System (ADS)
Morgan, Michael A.
Finite-difference and finite-element methods for the computational analysis of EM scattering phenomena are examined in chapters contributed by leading experts. Topics addressed include an FEM for composite scatterers, coupled finite- and boundary-element methods for EM scattering, absorbing boundary conditions for the direct solution PDEs arising in EM scattering problems, application of the control-region approximation to two-dimensional EM scattering, coupled potentials for EM fields in inhomogeneous media, the method of conforming boundary elements for transient electromagnetics, and the finite-difference time-domain method for numerical modeling of EM wave interactions with arbitrary structures. Extensive diagrams and graphs of typical results are provided.
Finite Orbit Width versions of the CQL3D code: Hybrid-FOW and Full-FOW
NASA Astrophysics Data System (ADS)
Petrov, Yu. V.; Harvey, R. W.
2012-10-01
Finite-Orbit-Width (FOW) effects are being added into the CQL3D bounce-averaged Fokker-Planck code [1] using two main options. In the Hybrid-FOW option, partial FOW capabilities are implemented which add FOW features into the particle source (NB) operator, RF quasilinear operator, diagnostics, and guiding center orbit losses with gyro-radius correction. Collisions remain Zero-Orbit-Width (ZOW). The Hybrid-FOW version provides a greatly improved agreement with signals measured by the NSTX Fast Ion Diagnostic [2]. The advantage of the Hybrid-FOW version is that run time increases by only a factor of two compared to ZOW runs. The Full-FOW option further adds neoclassical radial transport features into the FP coding. The collisional coefficients are averaged along guiding center orbits, with a proper transformation matrix from local coordinates to the midplane coordinates, where the FP equation is solved. All radial terms are included. The computations are parallelized in velocity-grid index, typically using 128 CPU cores. We emphasize that this theory includes nonthermal and full-orbit, not first order correction, neoclassical theory. [4pt] [1] R.W. Harvey and M. McCoy, ``The CQL3D Fokker Planck Code,'' www.compxco.com/cql3d [0pt] [2] R.W. Harvey, Yu. Petrov, D. Liu, W. Heidbrink, P. Bonoli, this mtg (2012)
The Combined Finite-Discrete Element Method applied to the Study of Rock Fracturing Behavior in 3D
Rougier, Esteban; Bradley, Christopher R.; Broom, Scott T.; Knight, Earl E.; Munjiza, Ante; Sussman, Aviva J.; Swift, Robert P.
2011-01-01
Since its introduction the combined finite-discrete element method (FEM/DEM), has become an excellent tool to address a wide range of problems involving fracturing and fragmentation of solids. Within the context of rock mechanics, the FEM/DEM method has been applied to many complex industrial problems such as block caving, deep mining techniques, rock blasting, seismic waves, packing problems, rock crushing problems, etc. In the real world most of the problems involving fracture and fragmentation of solids are three dimensional problems. With the aim of addressing these problems an improved 2D/3D FEM/DEM capability has been developed at Los Alamos National Laboratory (LANL). These capabilities include state of the art 3D contact detection, contact interaction, constitutive material models, and fracture models. In this paper, Split Hopkinson Pressure Bar (SHPB) Brazilian experiments are simulated using this improved 2D/3D FEM/DEM approach which is implemented in LANL's MUNROU (Munjiza-Rougier) code. The results presented in this work show excellent agreement with both the SHPB experiments and previous 2D numerical simulations performed by other FEM/DEM research groups.
NASA Astrophysics Data System (ADS)
Haddag, B.; Kagnaya, T.; Nouari, M.; Cutard, T.
2013-01-01
Modelling machining operations allows estimating cutting parameters which are difficult to obtain experimentally and in particular, include quantities characterizing the tool-workpiece interface. Temperature is one of these quantities which has an impact on the tool wear, thus its estimation is important. This study deals with a new modelling strategy, based on two steps of calculation, for analysis of the heat transfer into the cutting tool. Unlike the classical methods, considering only the cutting tool with application of an approximate heat flux at the cutting face, estimated from experimental data (e.g. measured cutting force, cutting power), the proposed approach consists of two successive 3D Finite Element calculations and fully independent on the experimental measurements; only the definition of the behaviour of the tool-workpiece couple is necessary. The first one is a 3D thermomechanical modelling of the chip formation process, which allows estimating cutting forces, chip morphology and its flow direction. The second calculation is a 3D thermal modelling of the heat diffusion into the cutting tool, by using an adequate thermal loading (applied uniform or non-uniform heat flux). This loading is estimated using some quantities obtained from the first step calculation, such as contact pressure, sliding velocity distributions and contact area. Comparisons in one hand between experimental data and the first calculation and at the other hand between measured temperatures with embedded thermocouples and the second calculation show a good agreement in terms of cutting forces, chip morphology and cutting temperature.
NASA Astrophysics Data System (ADS)
Gu, Xuejun; Jelen, Urszula; Li, Jinsheng; Jia, Xun; Jiang, Steve B.
2011-06-01
Targeting at the development of an accurate and efficient dose calculation engine for online adaptive radiotherapy, we have implemented a finite-size pencil beam (FSPB) algorithm with a 3D-density correction method on graphics processing unit (GPU). This new GPU-based dose engine is built on our previously published ultrafast FSPB computational framework (Gu et al 2009 Phys. Med. Biol. 54 6287-97). Dosimetric evaluations against Monte Carlo dose calculations are conducted on ten IMRT treatment plans (five head-and-neck cases and five lung cases). For all cases, there is improvement with the 3D-density correction over the conventional FSPB algorithm and for most cases the improvement is significant. Regarding the efficiency, because of the appropriate arrangement of memory access and the usage of GPU intrinsic functions, the dose calculation for an IMRT plan can be accomplished well within 1 s (except for one case) with this new GPU-based FSPB algorithm. Compared to the previous GPU-based FSPB algorithm without 3D-density correction, this new algorithm, though slightly sacrificing the computational efficiency (~5-15% lower), has significantly improved the dose calculation accuracy, making it more suitable for online IMRT replanning.
Gu, Xuejun; Jelen, Urszula; Li, Jinsheng; Jia, Xun; Jiang, Steve B
2011-06-01
Targeting at the development of an accurate and efficient dose calculation engine for online adaptive radiotherapy, we have implemented a finite-size pencil beam (FSPB) algorithm with a 3D-density correction method on graphics processing unit (GPU). This new GPU-based dose engine is built on our previously published ultrafast FSPB computational framework (Gu et al 2009 Phys. Med. Biol. 54 6287-97). Dosimetric evaluations against Monte Carlo dose calculations are conducted on ten IMRT treatment plans (five head-and-neck cases and five lung cases). For all cases, there is improvement with the 3D-density correction over the conventional FSPB algorithm and for most cases the improvement is significant. Regarding the efficiency, because of the appropriate arrangement of memory access and the usage of GPU intrinsic functions, the dose calculation for an IMRT plan can be accomplished well within 1 s (except for one case) with this new GPU-based FSPB algorithm. Compared to the previous GPU-based FSPB algorithm without 3D-density correction, this new algorithm, though slightly sacrificing the computational efficiency (∼5-15% lower), has significantly improved the dose calculation accuracy, making it more suitable for online IMRT replanning. PMID:21558589
Orthodontic intrusion of maxillary incisors: a 3D finite element method study
Saga, Armando Yukio; Maruo, Hiroshi; Argenta, Marco André; Maruo, Ivan Toshio; Tanaka, Orlando Motohiro
2016-01-01
Objective: In orthodontic treatment, intrusion movement of maxillary incisors is often necessary. Therefore, the objective of this investigation is to evaluate the initial distribution patterns and magnitude of compressive stress in the periodontal ligament (PDL) in a simulation of orthodontic intrusion of maxillary incisors, considering the points of force application. Methods: Anatomic 3D models reconstructed from cone-beam computed tomography scans were used to simulate maxillary incisors intrusion loading. The points of force application selected were: centered between central incisors brackets (LOAD 1); bilaterally between the brackets of central and lateral incisors (LOAD 2); bilaterally distal to the brackets of lateral incisors (LOAD 3); bilaterally 7 mm distal to the center of brackets of lateral incisors (LOAD 4). Results and Conclusions: Stress concentrated at the PDL apex region, irrespective of the point of orthodontic force application. The four load models showed distinct contour plots and compressive stress values over the midsagittal reference line. The contour plots of central and lateral incisors were not similar in the same load model. LOAD 3 resulted in more balanced compressive stress distribution. PMID:27007765
3D finite element modelling of guided wave scattering at delaminations in composites
NASA Astrophysics Data System (ADS)
Murat, Bibi Intan Suraya; Fromme, Paul
2016-02-01
Carbon fiber laminate composites are increasingly used for aerospace structures as they offer a number of advantages including a good strength to weight ratio. However, impact during the operation and servicing of the aircraft can lead to barely visible and difficult to detect damage. Depending on the severity of the impact, delaminations can occur, reducing the load carrying capacity of the structure. Efficient nondestructive testing of composite panels can be achieved using guided ultrasonic waves propagating along the structure. The guided wave (A0 Lamb wave mode) scattering at delaminations was modeled using full three-dimensional Finite Element (FE) simulations. The influence of the delamination size was systematically investigated from a parameter study. A significant influence of the delamination width on the guided wave scattering was found, especially on the angular dependency of the scattered guided wave amplitude. The sensitivity of guided ultrasonic waves for the detection of delamination damage in composite panels is discussed.
Active tectonics in Taiwan: insights from a 3-D viscous finite element model
NASA Astrophysics Data System (ADS)
Sun, Yujun; Liu, Mian; Dong, Shuwen; Zhang, Huai; Shi, Yaolin
2015-12-01
Taiwan is a young orogenic belt with complex spatial distributions of deformation and earthquakes. We have constructed a three-dimensional finite element model to explore how the interplays between lithospheric structure and plate boundary processes control the distribution of stress and strain rates in the Taiwan region. The model assumes a liberalized power-law rheology and incorporates main lithospheric structures; the model domain is loaded by the present-day crustal velocity applied at its boundaries. The model successfully reproduces the main features of the GPS-measured strain rate patterns and the earthquake-indicated stress states in the Taiwan region. The best fitting model requires the viscosity of the lower crust to be two orders of magnitude lower than that of the upper crust and lithospheric mantle. The calculated deviatoric stress is high in regions of thrust faulting and low in regions of extensional and strike-slip faulting, consistent with the spatial pattern of seismic intensity in Taiwan.
Finite element methods of analysis for 3D inviscid compressible flows
NASA Technical Reports Server (NTRS)
Peraire, Jaime
1990-01-01
The applicants have developed a finite element based approach for the solution of three-dimensional compressible flows. The procedure enables flow solutions to be obtained on tetrahedral discretizations of computational domains of complex form. A further development was the incorporation of a solution adaptive mesh strategy in which the adaptivity is achieved by complete remeshing of the solution domain. During the previous year, the applicants were working with the Advanced Aerodynamics Concepts Branch at NASA Ames Research Center with an implementation of the basic meshing and solution procedure. The objective of the work to be performed over this twelve month period was the transfer of the adaptive mesh technology and also the undertaking of basic research into alternative flow algorithms for the Euler equations on unstructured meshes.
Simulation of dielectrophoretic assembly of carbon nanotubes using 3D finite element analysis.
Berger, S D; McGruer, N E; Adams, G G
2015-04-17
One of the most important methods for selective and repeatable assembly of carbon nanotubes (CNTs) is alternating current dielectrophoresis (DEP). This method has been demonstrated experimentally as a viable technique for nano-scale manufacturing of novel CNT based devices. Previous numerical analyses have studied the motion of nanotubes, the volume from which they are assembled, and the rate of assembly, but have been restricted by various simplifying assumptions. In this paper we present a method for simulating the motion and behavior of CNTs subjected to dielectrophoresis using a three-dimensional electrostatic finite element analysis. By including the CNT in the finite element model, we can accurately predict the effect of the CNT on the electric field and the resulting force distribution across the CNT can be determined. We have used this information to calculate the motion of CNTs assembling onto the electrodes, and show how they tend to move towards the center of an electrode and come into contact at highly skewed angles. Our analysis suggests that the CNTs move to the electrode gap only after initially contacting the electrodes. We have also developed a model of the elastic deformation of CNTs as they approach the electrodes demonstrating how the induced forces can significantly alter the CNT shape during assembly. These results show that the CNT does not behave as a rigid body when in close proximity to the electrodes. In the future this method can be applied to a variety of real electrode geometries on a case-by-case basis and will provide more detailed insight into the specific motion and assembly parameters necessary for effective DEP assembly. PMID:25804394
Simulation of dielectrophoretic assembly of carbon nanotubes using 3D finite element analysis
NASA Astrophysics Data System (ADS)
Berger, S. D.; McGruer, N. E.; Adams, G. G.
2015-04-01
One of the most important methods for selective and repeatable assembly of carbon nanotubes (CNTs) is alternating current dielectrophoresis (DEP). This method has been demonstrated experimentally as a viable technique for nano-scale manufacturing of novel CNT based devices. Previous numerical analyses have studied the motion of nanotubes, the volume from which they are assembled, and the rate of assembly, but have been restricted by various simplifying assumptions. In this paper we present a method for simulating the motion and behavior of CNTs subjected to dielectrophoresis using a three-dimensional electrostatic finite element analysis. By including the CNT in the finite element model, we can accurately predict the effect of the CNT on the electric field and the resulting force distribution across the CNT can be determined. We have used this information to calculate the motion of CNTs assembling onto the electrodes, and show how they tend to move towards the center of an electrode and come into contact at highly skewed angles. Our analysis suggests that the CNTs move to the electrode gap only after initially contacting the electrodes. We have also developed a model of the elastic deformation of CNTs as they approach the electrodes demonstrating how the induced forces can significantly alter the CNT shape during assembly. These results show that the CNT does not behave as a rigid body when in close proximity to the electrodes. In the future this method can be applied to a variety of real electrode geometries on a case-by-case basis and will provide more detailed insight into the specific motion and assembly parameters necessary for effective DEP assembly.
NASA Astrophysics Data System (ADS)
Moortgat, J.; Firoozabadi, A.
2013-12-01
Most problems of interest in hydrogeology and subsurface energy resources involve complex heterogeneous geological formations. Such domains are most naturally represented in numerical reservoir simulations by unstructured computational grids. Finite element methods are a natural choice to describe fluid flow on unstructured meshes, because the governing equations can be readily discretized for any grid-element geometry. In this work, we consider the challenging problem of fully compositional three-phase flow in 3D unstructured grids, discretized by tetrahedra, prisms, or hexahedra, and compare to simulations on 3D structured grids. We employ a combination of mixed hybrid finite element methods to solve for the pressure and flux fields in a fractional flow formulation, and higher-order discontinuous Galerkin methods for the mass transport equations. These methods are well suited to simulate flow in heterogeneous and fractured reservoirs, because they provide a globally continuous pressure and flux field, while allowing for sharp discontinuities in the phase properties, such as compositions and saturations. The increased accuracy from using higher-order methods improves the modeling of highly non-linear flow, such as gravitational and viscous fingering. We present several numerical examples to study convergence rates and the (lack of) sensitivity to gridding/mesh orientation, and mesh quality. These examples consider gravity depletion, water and gas injection in oil saturated subsurface reservoirs with species exchange between up to three fluid phases. The examples demonstrate the wide applicability of our chosen finite element methods in the study of challenging multiphase flow problems in porous, geometrically complex, subsurface media.
Sutradhar, Alok; Park, Jaejong; Carrau, Diana; Miller, Michael J
2014-09-01
With the dawn of 3D printing technology, patient-specific implant designs are set to have a paradigm shift. A topology optimization method in designing patient-specific craniofacial implants has been developed to ensure adequate load transfer mechanism and restore the form and function of the mid-face. Patient-specific finite element models are used to design these implants and to validate whether they are viable for physiological loading such as mastication. Validation of these topology optimized finite element models using mechanical testing is a critical step. Instead of inserting the implants into a cadaver or patient, we embed the implants into the computer-aided skull model of a patient and, fuse them together to 3D print the complete skull model with the implant. Masticatory forces are applied in the molar region to simulate chewing and measure the stress-strain trajectory. Until recently, strain gages have been used to measure strains for validation. Digital Image Correlation (DIC) method is a relatively new technique for full-field strain measurement which provides a continuous deformation field data. The main objective of this study is to validate the finite element model of patient-specific craniofacial implants against the strain data from the DIC obtained during the mastication simulation and show that the optimized shapes provide adequate load-transfer mechanism. Patient-specific models are obtained from CT scans. The principal maximum and minimum strains are compared. The computational and experimental approach to designing patient-specific implants proved to be a viable technique for mid-face craniofacial reconstruction. PMID:24992729
Bone stress and strain modification in diastema closure: 3D analysis using finite element method.
Geramy, Allahyar; Bouserhal, Joseph; Martin, Domingo; Baghaeian, Pedram
2015-09-01
The aim of this study was to analyse the stress and strain distribution in the alveolar bone between two central incisors in the process of diastema closure with a constant force. A 3-dimensional computer modeling based on finite element techniques was used for this purpose. A model of an anterior segment of the mandible containing cortical bone, spongy bone, gingivae, PDL and two central incisors with a bracket in the labial surface of each tooth were designed. The von Mises stress and strain was evaluated in alveolar bone along a path of nodes defined in a cresto-apical direction in the midline between two teeth. It was observed that stress and strain of alveolar bone increased in midline with a constant force to close the diastema regardless of the type of movement in gradual steps of diastema closure, however the stress was higher in the tipping movement than the bodily so it can be suggested that a protocol of force system modification should be introduced to compensate for the stress and strain changes caused by the reduced distance to avoid the unwanted stress alteration during the diastema closure. PMID:26277458
Electron scattering from large molecules: a 3d finite element R-matrix approach
NASA Astrophysics Data System (ADS)
Tonzani, Stefano; Greene, Chris H.
2005-05-01
To solve the Schr"odinger equation for scattering of a low energy electron from a molecule, we present a three-dimensional finite element R-matrix method [S. Tonzani and C. H. Greene, J. Chem. Phys. 122 01411, (2005)]. Using the static exchange and local density approximations, we can use directly the molecular potentials extracted from ab initio codes (GAUSSIAN 98 in the work described here). A local polarization potential based on density functional theory [F. A. Gianturco and A. Rodriguez-Ruiz, Phys. Rev. A 47, 1075 (1993)] approximately describes the long range attraction to the molecular target induced by the scattering electron without adjustable parameters. We have used this approach successfully in calculations of cross sections for small and medium sized molecules (like SF6, XeF6, C60 and Uracil). This method will be useful to treat the electron-induced dynamics of extended molecular systems, possibly of biological interest, where oth er more complex ab initio methods are difficult to apply.
NASA Astrophysics Data System (ADS)
Korneev, V. G.
2012-09-01
BPS is a well known an efficient and rather general domain decomposition Dirichlet-Dirichlet type preconditioner, suggested in the famous series of papers Bramble, Pasciak and Schatz (1986-1989). Since then, it has been serving as the origin for the whole family of domain decomposition Dirichlet-Dirichlet type preconditioners-solvers as for h so hp discretizations of elliptic problems. For its original version, designed for h discretizations, the named authors proved the bound O(1 + log2 H/ h) for the relative condition number under some restricting conditions on the domain decomposition and finite element discretization. Here H/ h is the maximal relation of the characteristic size H of a decomposition subdomain to the mesh parameter h of its discretization. It was assumed that subdomains are images of the reference unite cube by trilinear mappings. Later similar bounds related to h discretizations were proved for more general domain decompositions, defined by means of coarse tetrahedral meshes. These results, accompanied by the development of some special tools of analysis aimed at such type of decompositions, were summarized in the book of Toselli and Widlund (2005). This paper is also confined to h discretizations. We further expand the range of admissible domain decompositions for constructing BPS preconditioners, in which decomposition subdomains can be convex polyhedrons, satisfying some conditions of shape regularity. We prove the bound for the relative condition number with the same dependence on H/ h as in the bound given above. Along the way to this result, we simplify the proof of the so called abstract bound for the relative condition number of the domain decomposition preconditioner. In the part, related to the analysis of the interface sub-problem preconditioning, our technical tools are generalization of those used by Bramble, Pasciak and Schatz.
NASA Astrophysics Data System (ADS)
Hafiz Mahayudin, Mohd; Che Mat, Ruzinoor
2016-06-01
The main objective of this paper is to discuss on the effectiveness of visualising terrain draped with Unmanned Aerial Vehicle (UAV) images generated from different contour intervals using Unity 3D game engine in online environment. The study area that was tested in this project was oil palm plantation at Sintok, Kedah. The contour data used for this study are divided into three different intervals which are 1m, 3m and 5m. ArcGIS software were used to clip the contour data and also UAV images data to be similar size for the overlaying process. The Unity 3D game engine was used as the main platform for developing the system due to its capabilities which can be launch in different platform. The clipped contour data and UAV images data were process and exported into the web format using Unity 3D. Then process continue by publishing it into the web server for comparing the effectiveness of different 3D terrain data (contour data) draped with UAV images. The effectiveness is compared based on the data size, loading time (office and out-of-office hours), response time, visualisation quality, and frame per second (fps). The results were suggest which contour interval is better for developing an effective online 3D terrain visualisation draped with UAV images using Unity 3D game engine. It therefore benefits decision maker and planner related to this field decide on which contour is applicable for their task.
NASA Astrophysics Data System (ADS)
Vattré, A.; Devincre, B.; Feyel, F.; Gatti, R.; Groh, S.; Jamond, O.; Roos, A.
2014-02-01
A unified model coupling 3D dislocation dynamics (DD) simulations with the finite element (FE) method is revisited. The so-called Discrete-Continuous Model (DCM) aims to predict plastic flow at the (sub-)micron length scale of materials with complex boundary conditions. The evolution of the dislocation microstructure and the short-range dislocation-dislocation interactions are calculated with a DD code. The long-range mechanical fields due to the dislocations are calculated by a FE code, taking into account the boundary conditions. The coupling procedure is based on eigenstrain theory, and the precise manner in which the plastic slip, i.e. the dislocation glide as calculated by the DD code, is transferred to the integration points of the FE mesh is described in full detail. Several test cases are presented, and the DCM is applied to plastic flow in a single-crystal Nickel-based superalloy.
Finite-Orbit-Width version of the CQL3D for description of RF-enhanced neoclassical transport
NASA Astrophysics Data System (ADS)
Petrov, Yu. V.; Harvey, R. W.
2015-12-01
The CQL3D bounce-averaged Fokker-Planck (FP) code [l] has been upgraded to include Finite-Orbit-Width (FOW) effects. The calculations can be done either with a fast Hybrid-FOW option or with a slower but neoclassically complete full-FOW option. The neoclassical radial transport appears naturally in the full-FOW version by averaging the local collision coefficients along guiding center orbits, with a proper transformation matrix from local (R, Z) coordinates to the midplane computational coordinates, where the FP equation is solved. In a similar way, the local quasilinear rf diffusion terms give rise to additional radial transport of orbits. The results of validation tests for the full-FOW version are presented.
Edge-based finite element approach to the simulation of geoelectromagnetic induction in a 3-D sphere
NASA Astrophysics Data System (ADS)
Yoshimura, Ryokei; Oshiman, Naoto
2002-02-01
We present a new simulator based on an edge-based finite element method (FEM) for computing the global-scale electromagnetic (EM) induction responses in a 3-D conducting sphere excited by an external source current for a variety of frequencies. The formulation is in terms of the magnetic vector potential. The edge-element approach assigns the degrees of freedom to the edges rather than to the nodes of the element. This edge-element strictly satisfies the discontinuity of the normal boundary conditions without considering the enforced normal boundary conditions that are usually practiced in a node-based FEM. To verify our simulation code, we compare our results with those of other solvers for two test computations, corresponding to azimuthally symmetric and asymmetric models. The results are in good agreement with one another.
NASA Astrophysics Data System (ADS)
Umar Alkali, Adam; Lenggo Ginta, Turnad; Majdi Abdul-Rani, Ahmad
2015-04-01
This paper presents a 3D transient finite element modelling of the workpiece temperature field produced during the travelling heat sourced from oxyacetylene flame. The proposed model was given in terms of preheat-only test applicable during thermally enhanced machining using the oxyacetylene flame as a heat source. The FEA model as well as the experimental test investigated the surface temperature distribution on 316L stainless steel at scanning speed of 100mm/min, 125mm/min 160mm/min, 200mm/min and 250mm/min. The parametric properties of the heat source maintained constant are; lead distance Ld =10mm, focus height Fh=7.5mm, oxygen gas pressure Poxy=15psi and acetylene gas pressure Pacty=25psi. An experimental validation of the temperature field induced on type 316L stainless steel reveal that temperature distribution increases when the travelling speed decreases.
ZIP3D: An elastic and elastic-plastic finite-element analysis program for cracked bodies
NASA Technical Reports Server (NTRS)
Shivakumar, K. N.; Newman, J. C., Jr.
1990-01-01
ZIP3D is an elastic and an elastic-plastic finite element program to analyze cracks in three dimensional solids. The program may also be used to analyze uncracked bodies or multi-body problems involving contacting surfaces. For crack problems, the program has several unique features including the calculation of mixed-mode strain energy release rates using the three dimensional virtual crack closure technique, the calculation of the J integral using the equivalent domain integral method, the capability to extend the crack front under monotonic or cyclic loading, and the capability to close or open the crack surfaces during cyclic loading. The theories behind the various aspects of the program are explained briefly. Line-by-line data preparation is presented. Input data and results for an elastic analysis of a surface crack in a plate and for an elastic-plastic analysis of a single-edge-crack-tension specimen are also presented.
3D Finite Element Model for Writing Long-Period Fiber Gratings by CO2 Laser Radiation
Coelho, João M. P.; Nespereira, Marta; Abreu, Manuel; Rebordão, José
2013-01-01
In the last years, mid-infrared radiation emitted by CO2 lasers has become increasing popular as a tool in the development of long-period fiber gratings. However, although the development and characterization of the resulting sensing devices have progressed quickly, further research is still necessary to consolidate functional models, especially regarding the interaction between laser radiation and the fiber's material. In this paper, a 3D finite element model is presented to simulate the interaction between laser radiation and an optical fiber and to determine the resulting refractive index change. Dependence with temperature of the main parameters of the optical fiber materials (with special focus on the absorption of incident laser radiation) is considered, as well as convection and radiation losses. Thermal and residual stress analyses are made for a standard single mode fiber, and experimental results are presented. PMID:23941908
Characterization of 3D and planar Si diodes with different neutron converter materials
NASA Astrophysics Data System (ADS)
Mendicino, R.; Boscardin, M.; Carturan, S.; Betta, G.-F. Dalla; Palma, M. Dalla; Maggioni, G.; Quaranta, A.; Ronchin, S.
2015-10-01
In this paper, we report on the characterization of silicon 3D and planar sensors, coupled with different neutron converter materials, such as 10B, B104 C and 6LiF, with different deposition thickness. Selected results from the electrical and functional characterization of the devices are shown and comparatively discussed with the aid of SRIM and Geant4 simulations. The limited neutron detection efficiency, on the order of ≃ 1% (planar) and ≃ 8% (3D) from simulations, is understood, and hints for the optimization of the devices have been derived.
On the Finite-Time Splash and Splat Singularities for the 3-D Free-Surface Euler Equations
NASA Astrophysics Data System (ADS)
Coutand, Daniel; Shkoller, Steve
2014-01-01
We prove that the 3-D free-surface incompressible Euler equations with regular initial geometries and velocity fields have solutions which can form a finite-time "splash" (or "splat") singularity first introduced in Castro et al. (Splash singularity for water waves, http://arxiv.org/abs/1106.2120v2, 2011), wherein the evolving 2-D hypersurface, the moving boundary of the fluid domain, self-intersects at a point (or on surface). Such singularities can occur when the crest of a breaking wave falls unto its trough, or in the study of drop impact upon liquid surfaces. Our approach is founded upon the Lagrangian description of the free-boundary problem, combined with a novel approximation scheme of a finite collection of local coordinate charts; as such we are able to analyze a rather general set of geometries for the evolving 2-D free-surface of the fluid. We do not assume the fluid is irrotational, and as such, our method can be used for a number of other fluid interface problems, including compressible flows, plasmas, as well as the inclusion of surface tension effects.
Desai, Shrikar R.; Karthikeyan, I.; Gaddale, Reetika
2013-01-01
Purpose: The purpose of this finite element study was to compare the stresses, strains, and displacements of double versus single implant in immediate loading for replacing mandibular molar. Materials and Methods: Two 3D FEM (finite element method) models were made to simulate implant designs. The first model used 5-mm-wide diameter implant to support a single molar crown. The second model used 3.75-3.75 double implant design. Anisotropic properties were assigned to bone model. Each model was analyzed with single force magnitude (100 N) in vertical axis. Results: This FEM study suggested that micromotion can be controlled better for double implants compared to single wide-diameter implants. The Von Mises stress for double implant showed 74.44% stress reduction compared to that of 5-mm implant. The Von Mises elastic strain was reduced by 61% for double implant compared to 5-mm implant. Conclusion: Within the limitations of the study, when the mesiodistal space for artificial tooth is more than 12.5 mm, under immediate loading, the double implant support should be considered. PMID:24554890
NASA Astrophysics Data System (ADS)
Gómez-Gutiérrez, Álvaro; Susanne, Schnabel; Conoscenti, Christian; Caraballo-Arias, Nathalie A.; Ferro, Vito; di Stefano, Constanza; Juan de Sanjosé, José; Berenguer-Sempere, Fernando; de Matías, Javier
2014-05-01
Recent developments made in tri-dimensional photo-reconstruction techniques (3D-PR), such as the use of Structure from Motion (SfM) and MultiView Stereo (MVS) techniques together, have allowed obtaining high resolution 3D point clouds. In order to achieve final point clouds with these techniques, only oblique images from consumer un-calibrated and non-metric cameras are needed. Here, these techniques are used in order to measure, monitor and quantify geomorphological features and processes. Three different applications through a range of scales and landforms are presented here. Firstly, five small gully headcuts located in a small catchment in SW Spain were monitored with the aim of estimating headcut retreat rates. During this field work, 3D models obtained by means of a Terrestrial Laser Scanner (TLS) were captured and used as benchmarks to analyze 3D-PR method accuracy. Results of this analysis showed centimeter-level accuracies with average distances between the 3D-PR model and the TLS model ranging from 0.009 to 0.025 m. Estimated soil loss ranged from -0.246 m3 to 0.114 m3 for a wet period (289 mm) of 54 days in 2013. Secondly, a calanchi type badland in Sicily (Italy) was photo-reconstructed and the quality of the 3D-PR model was analyzed using a Digital Elevation Model produced by classic digital photogrammetry with photos captured by an Unmanned Aerial Vehicle (UAV). In this case, sub-meter calculated accuracies (0.30) showed that it is possible to describe badland morphology using 3D-PR models but it is not feasible to use these models to quantify annual rates of soil erosion in badlands (10 mm eroded per year). Finally, a high-resolution model of the Veleta rock glacier (in SE Spain) was elaborated with 3D-PR techniques and compared with a 3D model obtained by means of a TLS. Results indicated that 3D-PR method can be applied to the micro-scale study of glacier morphologies and processes with average distances to the TLS point cloud of 0.21 m.
Shurbaji Mozayek, Rami; Allaf, Mirza; B. Abuharb, Mohammad
2016-01-01
Background. Long span is seen in many clinical situations. Treatmentplanning options of these cases are difficult and may require FPD, RPD or ISP. Each option has its own disadvantages, including mechanical problems, patient comfort and cost. This article will evaluate the stress distribution of a different treatment option, which consists of adding a single sup-porting implant to the FPD by using 3D finite element analysis. Methods. Three models, each consisting of 5 units, were created as follows: 1. Tooth Pontic Pontic Pontic Tooth; 2. Tooth Pontic Implant Pontic Tooth; 3. Tooth Pontic Pontic Implant Tooth. An axial force was applied to the prostheses by using 3D finite element method and stresses were evaluated. Results. The maximum stress was found in the prostheses in all the models; the highest stress values in all the shared components of the models were almost similar. Stress in implants was lower in the second model than the third one. Conclusion. Adding a supporting implant in long-span FPD has no advantages while it has the disadvantages of complicating treatment and the complications that may occur to the implant and surrounding bone. PMID:27429723
Shurbaji Mozayek, Rami; Allaf, Mirza; B Abuharb, Mohammad
2016-01-01
Background. Long span is seen in many clinical situations. Treatmentplanning options of these cases are difficult and may require FPD, RPD or ISP. Each option has its own disadvantages, including mechanical problems, patient comfort and cost. This article will evaluate the stress distribution of a different treatment option, which consists of adding a single sup-porting implant to the FPD by using 3D finite element analysis. Methods. Three models, each consisting of 5 units, were created as follows: 1. Tooth Pontic Pontic Pontic Tooth; 2. Tooth Pontic Implant Pontic Tooth; 3. Tooth Pontic Pontic Implant Tooth. An axial force was applied to the prostheses by using 3D finite element method and stresses were evaluated. Results. The maximum stress was found in the prostheses in all the models; the highest stress values in all the shared components of the models were almost similar. Stress in implants was lower in the second model than the third one. Conclusion. Adding a supporting implant in long-span FPD has no advantages while it has the disadvantages of complicating treatment and the complications that may occur to the implant and surrounding bone. PMID:27429723
Hybrid finite element-finite difference method for thermal analysis of blood vessels.
Blanchard, C H; Gutierrez, G; White, J A; Roemer, R B
2000-01-01
A hybrid finite-difference/finite-element technique for the thermal analysis of blood vessels embedded in perfused tissue has been developed and evaluated. This method provides efficient and accurate solutions to the conjugated heat transfer problem of convection by blood coupled to conduction in the tissue. The technique uses a previously developed 3D automatic meshing method for creating a finite element mesh in the tissue surrounding the vessels, coupled iteratively with a 1-D marching finite difference method for the interior of the vessels. This hybrid technique retains the flexibility and ease of automated finite-element meshing techniques for modelling the complex geometry of blood vessels and irregularly shaped tissues, and speeds the solution time by using a simple finite-difference method to calculate the bulk mean temperatures within all blood vessels. The use of the 1D finite-difference technique in the blood vessels also eliminates the large computer memory requirements needed to accurately solve large vessel network problems when fine FE meshes are used in the interior of vessels. The accuracy of the hybrid technique has been verified against previously verified numerical solutions. In summary, the hybrid technique combines the accuracy and flexibility found in automated finite-element techniques, with the speed and reduction of computational memory requirements associated with the 1D finite-difference technique, something which has not been done before. This method, thus, has the potential to provide accurate, flexible and relatively fast solutions for the thermal analysis of coupled perfusion/blood vessel problems, and large vessel network problems. PMID:10949130
Fessler, H.; Edwards, C.D.
1983-05-01
Combined strip and rosette gauge measurements and results from three-dimensional, finite element calculations are in excellent agreement with frozen stress photoelastic results for an efficient shape of cast-steel node under axial, brace loading. Three different meshes showed that two layers of elements through the thickness are needed.
NASA Technical Reports Server (NTRS)
Ponchak, George E.; Dalton, Edan; Tentzeris, Emmanouil M.; Papapolymerou, John; Williams, W. Dan (Technical Monitor)
2001-01-01
Three-dimensional circuits built upon multiple layers of polyimide are required for constructing Si/SiGe monolithic microwave/millimeter-wave integrated circuits on complementary metal oxide semiconductor (CMOS) (low resistivity) Si wafers. Thin film microstrip lines (TFMS) with finite width ground planes embedded in the polyimide are often used. However, the closely spaced TFMS lines are susceptible to high levels of coupling, which degrades circuit performance. In this paper, Finite Difference Time Domain (FDTD) analysis and experimental measurements are used to show that the ground planes must be connected by via holes to reduce coupling in both the forward and backward directions.
Pavarino, E.; Neves, L. A.; Machado, J. M.; de Godoy, M. F.; Shiyou, Y.; Momente, J. C.; Zafalon, G. F. D.; Pinto, A. R.; Valêncio, C. R.
2013-01-01
The Finite Element Method is a well-known technique, being extensively applied in different areas. Studies using the Finite Element Method (FEM) are targeted to improve cardiac ablation procedures. For such simulations, the finite element meshes should consider the size and histological features of the target structures. However, it is possible to verify that some methods or tools used to generate meshes of human body structures are still limited, due to nondetailed models, nontrivial preprocessing, or mainly limitation in the use condition. In this paper, alternatives are demonstrated to solid modeling and automatic generation of highly refined tetrahedral meshes, with quality compatible with other studies focused on mesh generation. The innovations presented here are strategies to integrate Open Source Software (OSS). The chosen techniques and strategies are presented and discussed, considering cardiac structures as a first application context. PMID:23762031
NASA Astrophysics Data System (ADS)
Simpson, D. G.; Lipatov, A. S.; Sittler, E. C.; Cooper, J. F.; Hartle, R. E.; Sarantos, M.
2012-12-01
In this report we discuss the results of a 3D hybrid modeling of the interaction between Saturn's magnetosphere and Titan's atmosphere/ionosphere for the T5 encounter. The T5 flyby is the only encounter when the two main ionizing sources of Titan's atmosphere, solar radiation and corotating plasma, align quasi-anti-parallel. The model is based on recent analysis of the Cassini Plasma Spectrometer (CAPS) and the Cassini Ion and Neutral Mass Spectrometer (INMS) measurements during the T5 flyby through Titan's ram-side and polar ionosphere [1,2]. Magnetic field data was used from the MAG instrument [3]. In our model the background ions (O+, H+), all pickup ions, and ionospheric ions are considered as a particles, whereas the electrons are described as a fluid (see e.g. [4]). Inhomogeneous photoionization (in the dayside ionosphere), electron-impact ionization, and charge exchange are included in our model. The temperature of the background electrons and pickup electrons was also incorporated into the generalized Ohm's law. We also take into account collisions between ions and neutrals. In our hybrid simulations we use Chamberlain profiles for the exosphere's components. The moon is considered as a weakly conducting body. The first results of our hybrid modeling show a strong asymmetry in the background (H+, O+) and pickup (H2+, N2+, CH4+) ion density profiles. Such strong asymmetry cannot be explained by a single-fluid multi-species 3D MHD model [5], which includes complex chemistry but does not produce finite gyroradius and kinetic effects. References [1] Sittler, et al., Energy Deposition Processes in Titan's Atmosphere and Its Induced Magnetosphere. In: Titan from Cassini-Huygens, Brown, R.H., Lebreton, J.P., Waite, J.H., Eds., Springer, (Dordrecht, Heidelberg, London, New York), pp. 393-455, 2010. [2] Agren, K., et al., On magnetosphere electron impact ionization and dynamics in Titan's ram-side and polar ionosphere -- a Cassini case study, Ann. Geophys., 25, 2359
NASA Astrophysics Data System (ADS)
Yonetsu, Daigo; Tanaka, Kazufumi; Hara, Takehisa
In recent years, induction-heating (IH) cookers that can be used to heat nonmagnetic metals such as aluminum have been produced. Occasionally, a light pan moves on a glass plate due to buoyancy when heated by an IH cooker. In some IH cookers, an aluminum plate is mounted between the glass plate and the coil in order to reduce the buoyancy effect. The objective of this research is to evaluate the buoyancy-reduction effect and the heating effect of buoyancy-reduction plates. Eddy current analysis is carried out by 3D finite element method, and the electromagnetic force and the heat distribution on the heating plate are calculated. After this calculation is performed, the temperature distribution of the heating plate is calculated by heat transfer analysis. It is found that the shape, area, and the position of the buoyancy reduction plate strongly affect the buoyancy and the heat distribution. The impact of the shape, area, and position of the buoyancy reduction plate was quantified. The phenomena in the heating were elucidated qualitatively.
Application of the Finite Orbit Width Version of the CQL3D Code to NBI +RF Heating of NSTX Plasma
NASA Astrophysics Data System (ADS)
Petrov, Yu. V.; Harvey, R. W.
2015-11-01
The CQL3D bounce-averaged Fokker-Planck (FP) code has been upgraded to include Finite-Orbit-Width (FOW) effects. The calculations can be done either with a fast Hybrid-FOW option or with a slower but neoclassically complete full-FOW option. The banana regime neoclassical radial transport appears naturally in the full-FOW version by averaging the local collision coefficients along guiding center orbits, with a proper transformation matrix from local (R, Z) coordinates to the midplane computational coordinates, where the FP equation is solved. In a similar way, the local quasilinear rf diffusion terms give rise to additional radial transport of orbits. The full-FOW version is applied to simulation of ion heating in NSTX plasma. It is demonstrated that it can describe the physics of transport phenomena in plasma with auxiliary heating, in particular, the enhancement of the radial transport of ions by RF heating and the occurrence of the bootstrap current. Because of the bounce-averaging on the FPE, the results are obtained in a relatively short computational time. A typical full-FOW run time is 30 min using 140 MPI cores. Due to an implicit solver, calculations with a large time step (tested up to dt = 0.5 sec) remain stable. Supported by USDOE grants SC0006614, ER54744, and ER44649.
NASA Astrophysics Data System (ADS)
Zhang, Qi-Hua
2015-10-01
Finite element generation of complicated fracture networks is the core issue and source of technical difficulty in three-dimensional (3-D) discrete fracture network (DFN) flow models. Due to the randomness and uncertainty in the configuration of a DFN, the intersection lines (traces) are arbitrarily distributed in each face (fracture and other surfaces). Hence, subdivision of the fractures is an issue relating to subdivision of two-dimensional (2-D) domains with arbitrarily-distributed constraints. When the DFN configuration is very complicated, the well-known approaches (e.g. Voronoi Delaunay-based methods and advancing-front techniques) cannot operate properly. This paper proposes an algorithm to implement end-to-end connection between traces to subdivide 2-D domains into closed loops. The compositions of the vertices in the common edges between adjacent loops (which may belong to a single fracture or two connected fractures) are thus ensured to be topologically identical. The paper then proposes an approach for triangulating arbitrary loops which does not add any nodes to ensure consistency of the meshes at the common edges. In addition, several techniques relating to tolerance control and improving code robustness are discussed. Finally, the equivalent permeability of the rock mass is calculated for some very complicated DFNs (the DFN may contain 1272 fractures, 633 connected fractures, and 16,270 closed loops). The results are compared with other approaches to demonstrate the veracity and efficiency of the approach proposed in this paper.
NASA Technical Reports Server (NTRS)
Nakazawa, Shohei
1991-01-01
Formulations and algorithms implemented in the MHOST finite element program are discussed. The code uses a novel concept of the mixed iterative solution technique for the efficient 3-D computations of turbine engine hot section components. The general framework of variational formulation and solution algorithms are discussed which were derived from the mixed three field Hu-Washizu principle. This formulation enables the use of nodal interpolation for coordinates, displacements, strains, and stresses. Algorithmic description of the mixed iterative method includes variations for the quasi static, transient dynamic and buckling analyses. The global-local analysis procedure referred to as the subelement refinement is developed in the framework of the mixed iterative solution, of which the detail is presented. The numerically integrated isoparametric elements implemented in the framework is discussed. Methods to filter certain parts of strain and project the element discontinuous quantities to the nodes are developed for a family of linear elements. Integration algorithms are described for linear and nonlinear equations included in MHOST program.
Stochastic finite-difference time-domain
NASA Astrophysics Data System (ADS)
Smith, Steven Michael
2011-12-01
This dissertation presents the derivation of an approximate method to determine the mean and the variance of electro-magnetic fields in the body using the Finite-Difference Time-Domain (FDTD) method. Unlike Monte Carlo analysis, which requires repeated FDTD simulations, this method directly computes the variance of the fields at every point in space at every sample of time in the simulation. This Stochastic FDTD simulation (S-FDTD) has at its root a new wave called the Variance wave, which is computed in the time domain along with the mean properties of the model space in the FDTD simulation. The Variance wave depends on the electro-magnetic fields, the reflections and transmission though the different dielectrics, and the variances of the electrical properties of the surrounding materials. Like the electro-magnetic fields, the Variance wave begins at zero (there is no variance before the source is turned on) and is computed in the time domain until all fields reach steady state. This process is performed in a fraction of the time of a Monte Carlo simulation and yields the first two statistical parameters (mean and variance). The mean of the field is computed using the traditional FDTD equations. Variance is computed by approximating the correlation coefficients between the constituitive properties and the use of the S-FDTD equations. The impetus for this work was the simulation time it takes to perform 3D Specific Absorption Rate (SAR) FDTD analysis of the human head model for cell phone power absorption in the human head due to the proximity of a cell phone being used. In many instances, Monte Carlo analysis is not performed due to the lengthy simulation times required. With the development of S-FDTD, these statistical analyses could be performed providing valuable statistical information with this information being provided in a small fraction of the time it would take to perform a Monte Carlo analysis.
Ultra-Wideband Time-Difference-of-Arrival High Resolution 3D Proximity Tracking System
NASA Technical Reports Server (NTRS)
Ni, Jianjun; Arndt, Dickey; Ngo, Phong; Phan, Chau; Dekome, Kent; Dusl, John
2010-01-01
This paper describes a research and development effort for a prototype ultra-wideband (UWB) tracking system that is currently under development at NASA Johnson Space Center (JSC). The system is being studied for use in tracking of lunar./Mars rovers and astronauts during early exploration missions when satellite navigation systems are not available. U IATB impulse radio (UWB-IR) technology is exploited in the design and implementation of the prototype location and tracking system. A three-dimensional (3D) proximity tracking prototype design using commercially available UWB products is proposed to implement the Time-Difference- Of-Arrival (TDOA) tracking methodology in this research effort. The TDOA tracking algorithm is utilized for location estimation in the prototype system, not only to exploit the precise time resolution possible with UWB signals, but also to eliminate the need for synchronization between the transmitter and the receiver. Simulations show that the TDOA algorithm can achieve the fine tracking resolution with low noise TDOA estimates for close-in tracking. Field tests demonstrated that this prototype UWB TDOA High Resolution 3D Proximity Tracking System is feasible for providing positioning-awareness information in a 3D space to a robotic control system. This 3D tracking system is developed for a robotic control system in a facility called "Moonyard" at Honeywell Defense & System in Arizona under a Space Act Agreement.
Ausiello, P; Apicella, A; Davidson, C L; Rengo, S
2001-10-01
The combination of diverse materials and complex geometry makes stress distribution analysis in teeth very complicated. Simulation in a computerized model might enable a study of the simultaneous interaction of the many variables. A 3D solid model of a human maxillary premolar was prepared and exported into a 3D-finite element model (FEM). Additionally, a generic class II MOD cavity preparation and restoration was simulated in the FEM model by a proper choice of the mesh volumes. A validation procedure of the FEM model was executed based on a comparison of theoretical calculations and experimental data. Different rigidities were assigned to the adhesive system and restorative materials. Two different stress conditions were simulated: (a) stresses arising from the polymerization shrinkage and (b) stresses resulting from shrinkage stress in combination with vertical occlusal loading. Three different cases were analyzed: a sound tooth, a tooth with a class II MOD cavity, adhesively restored with a high (25 GPa) and one with a low (12.5GPa) elastic modulus composite. The cusp movements induced by polymerization stress and (over)-functional occlusal loading were evaluated. While cusp displacement was higher for the more rigid composites due to the pre-stressing from polymerization shrinkage, cusp movements turned out to be lower for the more flexible composites in case the restored tooth which was stressed by the occlusal loading. This preliminary study by 3D FEA on adhesively restored teeth with a class II MOD cavity indicated that Young's modulus values of the restorative materials play an essential role in the success of the restoration. Premature failure due to stresses arising from polymerization shrinkage and occlusal loading can be prevented by proper selection and combination of materials. PMID:11522306
NASA Astrophysics Data System (ADS)
Kajzer, A.; Pozorski, J.; Szewc, K.
2014-08-01
In the paper we present Large-eddy simulation (LES) results of 3D Taylor- Green vortex obtained by the three different computational approaches: Smoothed Particle Hydrodynamics (SPH), Lattice Boltzmann Method (LBM) and Finite Volume Method (FVM). The Smagorinsky model was chosen as a subgrid-scale closure in LES for all considered methods and a selection of spatial resolutions have been investigated. The SPH and LBM computations have been carried out with the use of the in-house codes executed on GPU and compared, for validation purposes, with the FVM results obtained using the open-source CFD software OpenFOAM. A comparative study in terms of one-point statistics and turbulent energy spectra shows a good agreement of LES results for all methods. An analysis of the GPU code efficiency and implementation difficulties has been made. It is shown that both SPH and LBM may offer a significant advantage over mesh-based CFD methods.
Perceptual integration for qualitatively different 3-D cues in the human brain.
Dövencioğlu, Dicle; Ban, Hiroshi; Schofield, Andrew J; Welchman, Andrew E
2013-09-01
The visual system's flexibility in estimating depth is remarkable: We readily perceive 3-D structure under diverse conditions from the seemingly random dots of a "magic eye" stereogram to the aesthetically beautiful, but obviously flat, canvasses of the Old Masters. Yet, 3-D perception is often enhanced when different cues specify the same depth. This perceptual process is understood as Bayesian inference that improves sensory estimates. Despite considerable behavioral support for this theory, insights into the cortical circuits involved are limited. Moreover, extant work tested quantitatively similar cues, reducing some of the challenges associated with integrating computationally and qualitatively different signals. Here we address this challenge by measuring fMRI responses to depth structures defined by shading, binocular disparity, and their combination. We quantified information about depth configurations (convex "bumps" vs. concave "dimples") in different visual cortical areas using pattern classification analysis. We found that fMRI responses in dorsal visual area V3B/KO were more discriminable when disparity and shading concurrently signaled depth, in line with the predictions of cue integration. Importantly, by relating fMRI and psychophysical tests of integration, we observed a close association between depth judgments and activity in this area. Finally, using a cross-cue transfer test, we found that fMRI responses evoked by one cue afford classification of responses evoked by the other. This reveals a generalized depth representation in dorsal visual cortex that combines qualitatively different information in line with 3-D perception. PMID:23647559
Development and application of 3-D foot-shape measurement system under different loads
NASA Astrophysics Data System (ADS)
Liu, Guozhong; Wang, Boxiong; Shi, Hui; Luo, Xiuzhi
2008-03-01
The 3-D foot-shape measurement system under different loads based on laser-line-scanning principle was designed and the model of the measurement system was developed. 3-D foot-shape measurements without blind areas under different loads and the automatic extraction of foot-parameter are achieved with the system. A global calibration method for CCD cameras using a one-axis motion unit in the measurement system and the specialized calibration kits is presented. Errors caused by the nonlinearity of CCD cameras and other devices and caused by the installation of the one axis motion platform, the laser plane and the toughened glass plane can be eliminated by using the nonlinear coordinate mapping function and the Powell optimized method in calibration. Foot measurements under different loads for 170 participants were conducted and the statistic foot parameter measurement results for male and female participants under non-weight condition and changes of foot parameters under half-body-weight condition, full-body-weight condition and over-body-weight condition compared with non-weight condition are presented. 3-D foot-shape measurement under different loads makes it possible to realize custom-made shoe-making and shows great prosperity in shoe design, foot orthopaedic treatment, shoe size standardization, and establishment of a feet database for consumers and athletes.
Sotelo, Julio; Urbina, Jesus; Valverde, Israel; Tejos, Cristian; Irarrazaval, Pablo; Andia, Marcelo E; Uribe, Sergio; Hurtado, Daniel E
2016-06-01
Several 2D methods have been proposed to estimate WSS and OSI from PC-MRI, neglecting the longitudinal velocity gradients that typically arise in cardiovascular flow, particularly on vessel geometries whose cross section and centerline orientation strongly vary in the axial direction. Thus, the contribution of longitudinal velocity gradients remains understudied. In this work, we propose a 3D finite-element method for the quantification of WSS and OSI from 3D-CINE PC-MRI that accounts for both in-plane and longitudinal velocity gradients. We demonstrate the convergence and robustness of the method on cylindrical geometries using a synthetic phantom based on the Poiseuille flow equation. We also show that, in the presence of noise, the method is both stable and accurate. Using computational fluid dynamics simulations, we show that the proposed 3D method results in more accurate WSS estimates than those obtained from a 2D analysis not considering out-of-plane velocity gradients. Further, we conclude that for irregular geometries the accurate prediction of WSS requires the consideration of longitudinal gradients in the velocity field. Additionally, we compute 3D maps of WSS and OSI for 3D-CINE PC-MRI data sets from an aortic phantom and sixteen healthy volunteers and two patients. The OSI values show a greater dispersion than WSS, which is strongly dependent on the PC-MRI resolution. We envision that the proposed 3D method will improve the estimation of WSS and OSI from 3D-CINE PC-MRI images, allowing for more accurate estimates in vessels with pathologies that induce high longitudinal velocity gradients, such as coarctations and aneurisms. PMID:26780787
NASA Astrophysics Data System (ADS)
Yin, Y.; Hung, S.
2007-12-01
Seismic tomographic imaging has played a key component to unravel the deep processes that caused the surface morphology and rift magmatism in the southwest United States. Several studies used teleseismic body- wave arrivals recorded by the La Ristra experiment, a dense broadband array of 950-km in length deployed during 1999-2001 and run through the Great Plains, the Rio Grande Rift, and the Colorado Plateau, to construct a 2-D tomographic image of the upper mantle structure beneath this linear array (e.g., Gao et al., 2004). However, because of the inevitable smoothing and damping imposed in the tomographic model, the resulting velocity contrast is too weak to explain distinct P and S waveform changes across the array (Song and Helmberger, 2007). In this study, we include all the data from the La Ristra and available nearby arrays and reexamine finite- frequency travel time delays measured by inter-station cross correlation of waveforms at both high- (0.3-2 Hz for P and 0.1-0.5 Hz for S) and low-frequencies (0.03-0.125 Hz for P and 0.03-0.1 Hz for S). Differing from the previous models that rely on classical ray theory and simple grid parameterization, our inversion considers more realistic 3-D sensitivity kernels for relative travel-time delays and a wavelet-based, multi-scale parameterization that enables to yield robust features with spatially-varying resolutions. Our preliminary P-wave model reveals a prominent low-velocity zone extending from near surface to the depth of 300 km beneath the Rio Grande Rift, while the upper mantle which underlies the Great Plains and the Colorado Plateau is seismically fast. We will demonstrate the difference and improvement of 3-D tomographic models through the use of finite-frequency kernels and multi-scale parameterization.
Statistical 3D shape analysis of gender differences in lateral ventricles
NASA Astrophysics Data System (ADS)
He, Qing; Karpman, Dmitriy; Duan, Ye
2010-03-01
This paper aims at analyzing gender differences in the 3D shapes of lateral ventricles, which will provide reference for the analysis of brain abnormalities related to neurological disorders. Previous studies mostly focused on volume analysis, and the main challenge in shape analysis is the required step of establishing shape correspondence among individual shapes. We developed a simple and efficient method based on anatomical landmarks. 14 females and 10 males with matching ages participated in this study. 3D ventricle models were segmented from MR images by a semiautomatic method. Six anatomically meaningful landmarks were identified by detecting the maximum curvature point in a small neighborhood of a manually clicked point on the 3D model. Thin-plate spline was used to transform a randomly selected template shape to each of the rest shape instances, and the point correspondence was established according to Euclidean distance and surface normal. All shapes were spatially aligned by Generalized Procrustes Analysis. Hotelling T2 twosample metric was used to compare the ventricle shapes between males and females, and False Discovery Rate estimation was used to correct for the multiple comparison. The results revealed significant differences in the anterior horn of the right ventricle.
NASA Astrophysics Data System (ADS)
Kohout, B.; Pirinen, J.; Ruiter, N. V.
2012-03-01
The established standard screening method to detect breast cancer is X-ray mammography. However X-ray mammography often has low contrast for tumors located within glandular tissue. A new approach is 3D Ultrasound Computer Tomography (USCT), which is expected to detect small tumors at an early stage. This paper describes the development, improvement and the results of Finite Element Method (FEM) simulations of the Transducer Array System (TAS) used in our 3D USCT. The focus of this work is on researching the influence of meshing and material parameters on the electrical impedance curves. Thereafter, these findings are used to optimize the simulation model. The quality of the simulation was evaluated by comparing simulated impedance characteristics with measured data of the real TAS. The resulting FEM simulation model is a powerful tool to analyze and optimize transducer array systems applied for USCT. With this simulation model, the behavior of TAS for different geometry modifications was researched. It provides a means to understand the acoustical performances inside of any ultrasound transducer represented by its electrical impedance characteristic.
ERIC Educational Resources Information Center
Nazari, Mohammad Ali; Perrier, Pascal; Payan, Yohan
2013-01-01
Purpose: The authors aimed to design a distributed lambda model (DLM), which is well adapted to implement three-dimensional (3-D), finite-element descriptions of muscles. Method: A muscle element model was designed. Its stress-strain relationships included the active force-length characteristics of the ? model along the muscle fibers, together…
NASA Astrophysics Data System (ADS)
Kanellopoulos, V. N.; Webb, J. P.
1993-03-01
A 3D vector analysis of plane wave scattering by a metallic sphere using finite elements and Absorbing Boundary Conditions (ABCs) is presented. The ABCs are applied on the outer surface that truncates the infinitely extending domain. Mixed order curvilinear covariantprojection elements are used to avoid spurious corruptions. The second order ABC is superior to the first at no extra computational cost. The errors due to incomplete absorption decrease as the outer surface is moved further away from the scatterer. An error of about 1% in near-field values was obtained with the second order ABC, when the outer surface was less than half a wavelength from the scatterer. Une analyse tridimensionnelle vectorielle de la diffusion d'onde plane sur une sphère métallique utilisant des éléments finis et des Conditions aux Limites Absorbantes (CLA) est présentée. Les CLA sont appliquées sur la surface exteme tronquant le domaine s'étendant à l'infini. Des éléments curvilignes mixtes utilisant des projections covariantes sont utilisés pour éviter des solutions parasites. La CLA de second ordre est supérieure à celle de premier ordre sans effort de calcul additionnel. Les erreurs dues à l'absorption incomplète décroissent à mesure que l'on déplace la surface externe à une distance croissante du diffuseur. Un taux d'erreur d'environ 1 % dans les valeurs du champ proche a été obtenu avec les CLA de second ordre lorsque la surface externe était placée à une distance inférieure à une demi-longueur de la source de diffusion.
3D Scaffolds with Different Stiffness but the Same Microstructure for Bone Tissue Engineering.
Chen, Guobao; Dong, Chanjuan; Yang, Li; Lv, Yonggang
2015-07-29
A growing body of evidence has shown that extracellular matrix (ECM) stiffness can modulate stem cell adhesion, proliferation, migration, differentiation, and signaling. Stem cells can feel and respond sensitively to the mechanical microenvironment of the ECM. However, most studies have focused on classical two-dimensional (2D) or quasi-three-dimensional environments, which cannot represent the real situation in vivo. Furthermore, most of the current methods used to generate different mechanical properties invariably change the fundamental structural properties of the scaffolds (such as morphology, porosity, pore size, and pore interconnectivity). In this study, we have developed novel three-dimensional (3D) scaffolds with different degrees of stiffness but the same 3D microstructure that was maintained by using decellularized cancellous bone. Mixtures of collagen and hydroxyapatite [HA: Ca10(PO4)6(OH)2] with different proportions were coated on decellularized cancellous bone to vary the stiffness (local stiffness, 13.00 ± 5.55 kPa, 13.87 ± 1.51 kPa, and 37.7 ± 19.6 kPa; bulk stiffness, 6.74 ± 1.16 kPa, 8.82 ± 2.12 kPa, and 23.61 ± 8.06 kPa). Microcomputed tomography (μ-CT) assay proved that there was no statistically significant difference in the architecture of the scaffolds before or after coating. Cell viability, osteogenic differentiation, cell recruitment, and angiogenesis were determined to characterize the scaffolds and evaluate their biological responses in vitro and in vivo. The in vitro results indicate that the scaffolds developed in this study could sustain adhesion and growth of rat mesenchymal stem cells (MSCs) and promote their osteogenic differentiation. The in vivo results further demonstrated that these scaffolds could help to recruit MSCs from subcutaneous tissue, induce them to differentiate into osteoblasts, and provide the 3D environment for angiogenesis. These findings showed that the method we developed can build scaffolds with
Holford, D.J.
1994-01-01
This document is a user`s manual for the Rn3D finite element code. Rn3D was developed to simulate gas flow and radon transport in variably saturated, nonisothermal porous media. The Rn3D model is applicable to a wide range of problems involving radon transport in soil because it can simulate either steady-state or transient flow and transport in one-, two- or three-dimensions (including radially symmetric two-dimensional problems). The porous materials may be heterogeneous and anisotropic. This manual describes all pertinent mathematics related to the governing, boundary, and constitutive equations of the model, as well as the development of the finite element equations used in the code. Instructions are given for constructing Rn3D input files and executing the code, as well as a description of all output files generated by the code. Five verification problems are given that test various aspects of code operation, complete with example input files, FORTRAN programs for the respective analytical solutions, and plots of model results. An example simulation is presented to illustrate the type of problem Rn3D is designed to solve. Finally, instructions are given on how to convert Rn3D to simulate systems other than radon, air, and water.
Evaluation of the marginal seal of CEREC 3D restorations using two different luting agents.
Jahangiri, Leila; Agosta, Claudine; Estafan, Denise
2007-01-01
This in vitro study used two resin cements to evaluate the microleakage of CEREC 3D restorations. An mesio-occluso-distal-lingual (MODL) onlay preparation was placed on an ivorine premolar tooth. Sixty elastomeric impressions of this preparation were used to fabricate resin dies. Imaging propellant was applied to the resin dies. Optical impressions were taken with the CEREC 3 and the CEREC 3D software and used to design and mill 60 onlays. The reliability coefficient was alpha = 0.76 with an examiner agreement rate of 96.4%. Restorations cemented with posterior composite demonstrated a mean microleakage value of 1.46 mm (SD+/- 1.64); by contrast, those restored with dual-cure cement reported a mean microleakage value of 0.21 mm (SD+/-0.57). The difference in microleakage was very highly significant (alpha = 0.001). Posterior composite CEREC 3D restorations cemented with the dual-cure cement showed significantly less microleakage compared to those cemented with posterior composite. PMID:17333982
Finite-difference computations of rotor loads
NASA Technical Reports Server (NTRS)
Caradonna, F. X.; Tung, C.
1985-01-01
This paper demonstrates the current and future potential of finite-difference methods for solving real rotor problems which now rely largely on empiricism. The demonstration consists of a simple means of combining existing finite-difference, integral, and comprehensive loads codes to predict real transonic rotor flows. These computations are performed for hover and high-advance-ratio flight. Comparisons are made with experimental pressure data.
Finite-difference computations of rotor loads
NASA Technical Reports Server (NTRS)
Caradonna, F. X.; Tung, C.
1985-01-01
The current and future potential of finite difference methods for solving real rotor problems which now rely largely on empiricism are demonstrated. The demonstration consists of a simple means of combining existing finite-difference, integral, and comprehensive loads codes to predict real transonic rotor flows. These computations are performed for hover and high-advanced-ratio flight. Comparisons are made with experimental pressure data.
Verri, Fellippo Ramos; Cruz, Ronaldo Silva; de Souza Batista, Victor Eduardo; Almeida, Daniel Augusto de Faria; Verri, Ana Caroline Gonçales; Lemos, Cleidiel Aparecido de Araújo; Santiago Júnior, Joel Ferreira; Pellizzer, Eduardo Piza
2016-11-01
The aim of this study was to assess stress/strain of different implant modeling simplifications by 3D-FEA. Three variation of external hexagon implant (Ø3.75 × 10 mm) supporting one molar crown were simulated: A (no threads); B (slightly threads simplification); C (original design). 200 N (axial) and 100 N (oblique) were applied. Cortical bone was evaluated by maximum principal stress and microstrain qualitatively and quantitatively (ANOVA and Tukey post hoc (p < 0.05)). Higher stress levels (p < 0.05) were observed in model A. Models B and C presented similar stress transmission. It was possible to conclude that slightly simplification should be used for studies evaluating stress transferring for bone tissue. PMID:27082041
The Performance Evaluation of Multi-Image 3d Reconstruction Software with Different Sensors
NASA Astrophysics Data System (ADS)
Mousavi, V.; Khosravi, M.; Ahmadi, M.; Noori, N.; Naveh, A. Hosseini; Varshosaz, M.
2015-12-01
Today, multi-image 3D reconstruction is an active research field and generating three dimensional model of the objects is one the most discussed issues in Photogrammetry and Computer Vision that can be accomplished using range-based or image-based methods. Very accurate and dense point clouds generated by range-based methods such as structured light systems and laser scanners has introduced them as reliable tools in the industry. Image-based 3D digitization methodologies offer the option of reconstructing an object by a set of unordered images that depict it from different viewpoints. As their hardware requirements are narrowed down to a digital camera and a computer system, they compose an attractive 3D digitization approach, consequently, although range-based methods are generally very accurate, image-based methods are low-cost and can be easily used by non-professional users. One of the factors affecting the accuracy of the obtained model in image-based methods is the software and algorithm used to generate three dimensional model. These algorithms are provided in the form of commercial software, open source and web-based services. Another important factor in the accuracy of the obtained model is the type of sensor used. Due to availability of mobile sensors to the public, popularity of professional sensors and the advent of stereo sensors, a comparison of these three sensors plays an effective role in evaluating and finding the optimized method to generate three-dimensional models. Lots of research has been accomplished to identify a suitable software and algorithm to achieve an accurate and complete model, however little attention is paid to the type of sensors used and its effects on the quality of the final model. The purpose of this paper is deliberation and the introduction of an appropriate combination of a sensor and software to provide a complete model with the highest accuracy. To do this, different software, used in previous studies, were compared and
Whirley, R.G.; Engelmann, B.E.
1993-11-01
This report is the User Manual for the 1993 version of DYNA3D, and also serves as a User Guide. DYNA3D is a nonlinear, explicit, finite element code for analyzing the transient dynamic response of three-dimensional solids and structures. The code is fully vectorized and is available on several computer platforms. DYNA3D includes solid, shell, beam, and truss elements to allow maximum flexibility in modeling physical problems. Many material models 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 and single surface contact. Rigid materials provide added modeling flexibility. A material model driver with interactive graphics display is incorporated into DYNA3D to permit accurate modeling of complex material response based on experimental data. Along with the DYNA3D Example Problem Manual, this document provides the information necessary to apply DYNA3D to solve a wide range of engineering analysis problems.
Saha, Debashish; Babler, Matthaus U; Holzner, Markus; Soos, Miroslav; Lüthi, Beat; Liberzon, Alex; Kinzelbach, Wolfgang
2016-01-12
Aggregates grown in mild shear flow are released, one at a time, into homogeneous isotropic turbulence, where their motion and intermittent breakup is recorded by three-dimensional particle tracking velocimetry (3D-PTV). The aggregates have an open structure with a fractal dimension of ∼2.2, and their size is 1.4 ± 0.4 mm, which is large, compared to the Kolmogorov length scale (η = 0.15 mm). 3D-PTV of flow tracers allows for the simultaneous measurement of aggregate trajectories and the full velocity gradient tensor along their pathlines, which enables us to access the Lagrangian stress history of individual breakup events. From this data, we found no consistent pattern that relates breakup to the local flow properties at the point of breakup. Also, the correlation between the aggregate size and both shear stress and normal stress at the location of breakage is found to be weaker, when compared with the correlation between size and drag stress. The analysis suggests that the aggregates are mostly broken due to the accumulation of the drag stress over a time lag on the order of the Kolmogorov time scale. This finding is explained by the fact that the aggregates are large, which gives their motion inertia and increases the time for stress propagation inside the aggregate. Furthermore, it is found that the scaling of the largest fragment and the accumulated stress at breakup follows an earlier established power law, i.e., dfrag ∼ σ(-0.6) obtained from laminar nozzle experiments. This indicates that, despite the large size and the different type of hydrodynamic stress, the microscopic mechanism causing breakup is consistent over a wide range of aggregate size and stress magnitude. PMID:26646289
NASA Technical Reports Server (NTRS)
Parikh, Paresh; Pirzadeh, Shahyar; Loehner, Rainald
1990-01-01
A set of computer programs for 3-D unstructured grid generation, fluid flow calculations, and flow field visualization was developed. The grid generation program, called VGRID3D, generates grids over complex configurations using the advancing front method. In this method, the point and element generation is accomplished simultaneously, VPLOT3D is an interactive, menudriven pre- and post-processor graphics program for interpolation and display of unstructured grid data. The flow solver, VFLOW3D, is an Euler equation solver based on an explicit, two-step, Taylor-Galerkin algorithm which uses the Flux Corrected Transport (FCT) concept for a wriggle-free solution. Using these programs, increasingly complex 3-D configurations of interest to aerospace community were gridded including a complete Space Transportation System comprised of the space-shuttle orbitor, the solid-rocket boosters, and the external tank. Flow solutions were obtained on various configurations in subsonic, transonic, and supersonic flow regimes.
FDIPS: Finite Difference Iterative Potential-field Solver
NASA Astrophysics Data System (ADS)
Toth, Gabor; van der Holst, Bartholomeus; Huang, Zhenguang
2016-06-01
FDIPS is a finite difference iterative potential-field solver that can generate the 3D potential magnetic field solution based on a magnetogram. It is offered as an alternative to the spherical harmonics approach, as when the number of spherical harmonics is increased, using the raw magnetogram data given on a grid that is uniform in the sine of the latitude coordinate can result in inaccurate and unreliable results, especially in the polar regions close to the Sun. FDIPS is written in Fortran 90 and uses the MPI library for parallel execution.
Inter-point procrustes: identifying regional and large differences in 3D anatomical shapes.
Lekadir, Karim; Frangi, Alejandro F; Yang, Guang-Zhong
2012-01-01
This paper presents a new approach for the robust alignment and interpretation of 3D anatomical structures with large and localized shape differences. In such situations, existing techniques based on the well-known Procrustes analysis can be significantly affected due to the introduced non-Gaussian distribution of the residuals. In the proposed technique, influential points that induce large dissimilarities are identified and displaced with the aim to obtain an intermediate template with an improved distribution of the residuals. The key element of the algorithm is the use of pose invariant shape variables to robustly guide both the influential point detection and displacement steps. The intermediate template is then used as the basis for the estimation of the final pose parameters between the source and destination shapes, enabling to effectively highlight the regional differences of interest. The validation using synthetic and real datasets of different morphologies demonstrates robustness up-to 50% regional differences and potential for shape classification. PMID:23286119
NASA Astrophysics Data System (ADS)
Vattré, A.; Denoual, C.
2016-07-01
A thermodynamically consistent framework for combining nonlinear elastoplasticity and multivariant phase-field theory is formulated at large strains. In accordance with the Clausius-Duhem inequality, the Helmholtz free energy and time-dependent constitutive relations give rise to displacive driving forces for pressure-induced martensitic phase transitions in materials. Inelastic forces are obtained by using a representation of the energy landscape that involves the concept of reaction pathways with respect to the point group symmetry operations of crystal lattices. On the other hand, additional elastic forces are derived for the most general case of large strains and rotations, as well as nonlinear, anisotropic, and different elastic pressure-dependent properties of phases. The phase-field formalism coupled with finite elastoplastic deformations is implemented into a three-dimensional Lagrangian finite element approach and is applied to analyze the iron body-centered cubic (α-Fe) into hexagonal close-packed (ɛ-Fe) phase transitions under high hydrostatic compression. The simulations exhibit the major role played by the plastic deformation in the morphological and microstructure evolution processes. Due to the strong long-range elastic interactions between variants without plasticity, a forward α → ɛ transition is energetically unfavorable and remains incomplete. However, plastic dissipation releases considerably the stored strain energy, leading to the α ↔ ɛ ↔α‧ (forward and reverse) polymorphic phase transformations with an unexpected selection of variants.
NASA Technical Reports Server (NTRS)
Krueger, Ronald; Paris, Isbelle L.; OBrien, T. Kevin; Minguet, Pierre J.
2004-01-01
The influence of two-dimensional finite element modeling assumptions on the debonding prediction for skin-stiffener specimens was investigated. Geometrically nonlinear finite element analyses using two-dimensional plane-stress and plane-strain elements as well as three different generalized plane strain type approaches were performed. The computed skin and flange strains, transverse tensile stresses and energy release rates were compared to results obtained from three-dimensional simulations. The study showed that for strains and energy release rate computations the generalized plane strain assumptions yielded results closest to the full three-dimensional analysis. For computed transverse tensile stresses the plane stress assumption gave the best agreement. Based on this study it is recommended that results from plane stress and plane strain models be used as upper and lower bounds. The results from generalized plane strain models fall between the results obtained from plane stress and plane strain models. Two-dimensional models may also be used to qualitatively evaluate the stress distribution in a ply and the variation of energy release rates and mixed mode ratios with delamination length. For more accurate predictions, however, a three-dimensional analysis is required.
NH3D(+) ions in ammonium perchlorate - The N-D stretching bands of different sites
NASA Astrophysics Data System (ADS)
Weier, Jane E.; Strauss, Herbert L.
1993-03-01
The N-D stretching region of the infrared spectrum of doped NH4ClO4 has been reexamined at low temperatures. The spectra are found to be still changing at 2 K, indicating the existence of energy differences of this order. The N-H spectrum in ND4ClO4 was also examined. Comparison of the two sets of spectra led to an assignment of the small energy differences as due to tunneling among the distinct positions the NX3Y(+) ions can assume in the lattice. The distinct orientation of NH3D(+) and the ND3H(+) ions have slightly different energies and this, plus the tunneling, accounts for all the eight observed lines. A complete assignment of the tunneling levels will, however, require more information.
Numerical computation of transonic flows by finite-element and finite-difference methods
NASA Technical Reports Server (NTRS)
Hafez, M. M.; Wellford, L. C.; Merkle, C. L.; Murman, E. M.
1978-01-01
Studies on applications of the finite element approach to transonic flow calculations are reported. Different discretization techniques of the differential equations and boundary conditions are compared. Finite element analogs of Murman's mixed type finite difference operators for small disturbance formulations were constructed and the time dependent approach (using finite differences in time and finite elements in space) was examined.
Finite-difference modelling of wavefield constituents
NASA Astrophysics Data System (ADS)
Robertsson, Johan O. A.; van Manen, Dirk-Jan; Schmelzbach, Cedric; Van Renterghem, Cederic; Amundsen, Lasse
2015-11-01
The finite-difference method is among the most popular methods for modelling seismic wave propagation. Although the method has enjoyed huge success for its ability to produce full wavefield seismograms in complex models, it has one major limitation which is of critical importance for many modelling applications; to naturally output up- and downgoing and P- and S-wave constituents of synthesized seismograms. In this paper, we show how such wavefield constituents can be isolated in finite-difference-computed synthetics in complex models with high numerical precision by means of a simple algorithm. The description focuses on up- and downgoing and P- and S-wave separation of data generated using an isotropic elastic finite-difference modelling method. However, the same principles can also be applied to acoustic, electromagnetic and other wave equations.
NASA Astrophysics Data System (ADS)
Da Viá, C.; Bolle, E.; Einsweiler, K.; Garcia-Sciveres, M.; Hasi, J.; Kenney, C.; Linhart, V.; Parker, Sherwood; Pospisil, S.; Rohne, O.; Slavicek, T.; Watts, S.; Wermes, N.
2009-06-01
3D detectors, with electrodes penetrating the entire silicon wafer and active edges, were fabricated at the Stanford Nano Fabrication Facility (SNF), California, USA, with different electrode configurations. After irradiation with neutrons up to a fluence of 8.8×10 15 n eq cm -2, they were characterised using an infrared laser tuned to inject ˜2 minimum ionising particles showing signal efficiencies as high as 66% for the configuration with the shortest (56 μm) inter-electrode spacing. Sensors from the same wafer were also bump-bonded to the ATLAS FE-I3 pixel readout chip and their noise characterised. Most probable signal-to-noise ratios were calculated before and after irradiation to be as good as 38:1 after the highest irradiation level with a substrate thickness of 210 μm. These devices are promising candidates for application at the LHC such as the very forward detectors at ATLAS and CMS, the ATLAS B-Layer replacement and the general pixel upgrade. Moreover, 3D sensors could play a role in applications where high speed, high-resolution detectors are required, such as the vertex locators at the proposed Compact Linear Collider (CLIC) at CERN.
Comparison of different 3D navigation systems by a clinical "user".
Cartellieri, M; Kremser, J; Vorbeck, F
2001-01-01
Three-dimensional navigation systems are routinely used in endoscopic skull base surgery, neurosurgery, maxillo-facial and endoscopic sinus surgery. Their precision can, however, change in the course of one experiment. We have compared five different 3D navigation systems and discuss here possible reasons for the limits of system precision. A plexiglass cube on which test points were marked served as a test-model. Two well-trained system users measured the distances between the test points in each of the five systems. The results were compared with reference data provided by the NUMEREX device at the Technical University of Vienna. The accuracy data shown by all these 3D navigation systems ranged from 0.0 mm to 6.67 mm. The accuracy data of a system calculated in advance did not always correspond with the system precision on the screen. The system precision in the center of the cube was higher than on its surface, which made us conclude that the angle between the tracker system and the pointing device touching the test point may be critical for system precision. Applying an automatic registration step did not result in greater system precision. Slice thickness and the angle of the pointing device seem to be responsible for system precision. PMID:11271433
NASA Astrophysics Data System (ADS)
Bada, Adedayo; Alcaraz-Calero, Jose M.; Wang, Qi; Grecos, Christos
2014-05-01
This paper describes a comprehensive empirical performance evaluation of 3D video processing employing the physical/virtual architecture implemented in a cloud environment. Different virtualization technologies, virtual video cards and various 3D benchmarks tools have been utilized in order to analyse the optimal performance in the context of 3D online gaming applications. This study highlights 3D video rendering performance under each type of hypervisors, and other factors including network I/O, disk I/O and memory usage. Comparisons of these factors under well-known virtual display technologies such as VNC, Spice and Virtual 3D adaptors reveal the strengths and weaknesses of the various hypervisors with respect to 3D video rendering and streaming.
Deriving 3d Point Clouds from Terrestrial Photographs - Comparison of Different Sensors and Software
NASA Astrophysics Data System (ADS)
Niederheiser, Robert; Mokroš, Martin; Lange, Julia; Petschko, Helene; Prasicek, Günther; Oude Elberink, Sander
2016-06-01
Terrestrial photogrammetry nowadays offers a reasonably cheap, intuitive and effective approach to 3D-modelling. However, the important choice, which sensor and which software to use is not straight forward and needs consideration as the choice will have effects on the resulting 3D point cloud and its derivatives. We compare five different sensors as well as four different state-of-the-art software packages for a single application, the modelling of a vegetated rock face. The five sensors represent different resolutions, sensor sizes and price segments of the cameras. The software packages used are: (1) Agisoft PhotoScan Pro (1.16), (2) Pix4D (2.0.89), (3) a combination of Visual SFM (V0.5.22) and SURE (1.2.0.286), and (4) MicMac (1.0). We took photos of a vegetated rock face from identical positions with all sensors. Then we compared the results of the different software packages regarding the ease of the workflow, visual appeal, similarity and quality of the point cloud. While PhotoScan and Pix4D offer the user-friendliest workflows, they are also "black-box" programmes giving only little insight into their processing. Unsatisfying results may only be changed by modifying settings within a module. The combined workflow of Visual SFM, SURE and CloudCompare is just as simple but requires more user interaction. MicMac turned out to be the most challenging software as it is less user-friendly. However, MicMac offers the most possibilities to influence the processing workflow. The resulting point-clouds of PhotoScan and MicMac are the most appealing.
Applications of an exponential finite difference technique
NASA Technical Reports Server (NTRS)
Handschuh, Robert F.; Keith, Theo G., Jr.
1988-01-01
An exponential finite difference scheme first presented by Bhattacharya for one dimensional unsteady heat conduction problems in Cartesian coordinates was extended. The finite difference algorithm developed was used to solve the unsteady diffusion equation in one dimensional cylindrical coordinates and was applied to two and three dimensional conduction problems in Cartesian coordinates. Heat conduction involving variable thermal conductivity was also investigated. The method was used to solve nonlinear partial differential equations in one and two dimensional Cartesian coordinates. Predicted results are compared to exact solutions where available or to results obtained by other numerical methods.
Djoudi, Farid
2013-01-01
Two separate themes are presented in this paper. Aims The first theme is to present a graphical modeling approach of human anatomical structures namely, the femur and the tibia. The second theme involves making a finite element analysis of stresses, displacements and deformations in prosthetic implants (the femoral implant and the polyethylene insert). Objectives The graphical modeling approach comes in two parts. The first is the segmentation of MRI scanned images, retrieved in DICOM format for edge detection. In the second part, 3D-CAD models are generated from the results of the segmentation stage. The finite element analysis is done by first extracting the prosthetic implants from the reconstructed 3D-CAD model, then do a finite element analysis of these implants under objectively determined conditions such as; forces, allowed displacements, the materials composing implant, and the coefficient of friction. Conclusion The objective of this work is to implement an interface for exchanging data between 2D MRI images obtained from a medical diagnosis of a patient and the 3D-CAD model used in various applications, such as; the extraction of the implants, stress analysis at the knee joint and can serve as an aid to surgery, also predict the behavior of the prosthetic implants vis-a-vis the forces acting on the knee joints. PMID:24396234
NASA Astrophysics Data System (ADS)
Yao, Hepeng; Luan, Tian; Li, Chen; Zhang, Yin; Ma, Zhaoyuan; Chen, Xuzong
2016-01-01
Pursuing ultralow temperature 3D atom gas under microgravity conditions is one of the popular topics in the field of ultracold research. Many groups around the world are using, or are planning to use, delta-kick cooling (DKC) in microgravity. Our group has also proposed a two-stage crossed beam cooling (TSCBC) method that also provides a path to picokelvin temperatures. In this paper, we compare the characteristics of TSCBC and DKC for producing a picokelvin system in microgravity. Using a direct simulation Monte Carlo (DSMC) method, we simulate the cooling process of 87Rb using the two different cooling techniques. Under the same initial conditions, 87Rb can reach 7 pK in 15 s using TSCBC and 75 pK in 5.1 s with DKC. The simulation results show that TSCBC can reach lower temperatures compared with DKC, but needs more time and a more stable laser.
Maiti, Abhik; Chakravarty, Debashish
2016-01-01
3D reconstruction of geo-objects from their digital images is a time-efficient and convenient way of studying the structural features of the object being modelled. This paper presents a 3D reconstruction methodology which can be used to generate photo-realistic 3D watertight surface of different irregular shaped objects, from digital image sequences of the objects. The 3D reconstruction approach described here is robust, simplistic and can be readily used in reconstructing watertight 3D surface of any object from its digital image sequence. Here, digital images of different objects are used to build sparse, followed by dense 3D point clouds of the objects. These image-obtained point clouds are then used for generation of photo-realistic 3D surfaces, using different surface reconstruction algorithms such as Poisson reconstruction and Ball-pivoting algorithm. Different control parameters of these algorithms are identified, which affect the quality and computation time of the reconstructed 3D surface. The effects of these control parameters in generation of 3D surface from point clouds of different density are studied. It is shown that the reconstructed surface quality of Poisson reconstruction depends on Samples per node (SN) significantly, greater SN values resulting in better quality surfaces. Also, the quality of the 3D surface generated using Ball-Pivoting algorithm is found to be highly depend upon Clustering radius and Angle threshold values. The results obtained from this study give the readers of the article a valuable insight into the effects of different control parameters on determining the reconstructed surface quality. PMID:27386376
Comparison of 3-D finite element model of ashlar masonry with 2-D numerical models of ashlar masonry
NASA Astrophysics Data System (ADS)
Beran, Pavel
2016-06-01
3-D state of stress in heterogeneous ashlar masonry can be also computed by several suitable chosen 2-D numerical models of ashlar masonry. The results obtained from 2-D numerical models well correspond to the results obtained from 3-D numerical model. The character of thermal stress is the same. While using 2-D models the computational time is reduced more than hundredfold and therefore this method could be used for computation of thermal stresses during long time periods with 10 000 of steps.
StrainModeler: A MATHEMATICA™-based program for 3D analysis of finite and progressive strain
NASA Astrophysics Data System (ADS)
Bobillo-Ares, Nilo C.; Aller, Jesús; Bastida, Fernando; Menéndez, Omar; Lisle, Richard J.
2015-05-01
StrainModeler is a program constructed in the MATHEMATICA™ environment that performs 3D progressive strain calculations for lines and planes undergoing any sequence of homogeneous deformations. The main inputs to the system define the initial line or plane to be deformed and the deformation sequence to be applied, including combinations of simple shear, pure shear and volume change. For the deformation of lines, the output of the program is the change of attitude of the initial line, which can be represented by graphics or plotted in an equal-area projection. For the deformation of planes, the program has several outputs: (i) change of attitude of the initial plane; (ii) magnitudes and ratio of the semi-axes of the strain ellipse on the deformed plane; (iii) orientation of the major and minor axes of the strain ellipse on the deformed plane; (iv) orientations of the axial planes of the folds formed on the deformed plane, and (v) area change on the deformed plane. The variation of any of these parameters can be shown against a linear parameter only linked to the number of steps involved in the deformation, as a kind of "time" line, or it can be shown against the variation of a parameter of the strain ellipsoid (e. g.: major axis/minor axis ratio). A sequence of directions can be also visualized as a curve in an equal-area plot. Three applications of the program are presented. In the first, the deformation by simple shear of a plane with any orientation is analyzed. In the second, we explore the formation of recumbent folds in layers with different initial orientations for simple shear and pure shear deformations. In the third, we use StrainModeler to analyze the deformation of a set of folds located in a ductile shear zone in the Variscan Belt of NW Spain.
Coupled 2D-3D finite element method for analysis of a skin panel with a discontinuous stiffener
NASA Technical Reports Server (NTRS)
Wang, J. T.; Lotts, C. G.; Davis, D. D., Jr.; Krishnamurthy, T.
1992-01-01
This paper describes a computationally efficient analysis method which was used to predict detailed stress states in a typical composite compression panel with a discontinuous hat stiffener. A global-local approach was used. The global model incorporated both 2D shell and 3D brick elements connected by newly developed transition elements. Most of the panel was modeled with 2D elements, while 3D elements were employed to model the stiffener flange and the adjacent skin. Both linear and geometrically nonlinear analyses were performed on the global model. The effect of geometric nonlinearity induced by the eccentric load path due to the discontinuous hat stiffener was significant. The local model used a fine mesh of 3D brick elements to model the region at the end of the stiffener. Boundary conditions of the local 3D model were obtained by spline interpolation of the nodal displacements from the global analysis. Detailed in-plane and through-the-thickness stresses were calculated in the flange-skin interface near the end of the stiffener.
3D-CAD Effects on Creative Design Performance of Different Spatial Abilities Students
ERIC Educational Resources Information Center
Chang, Y.
2014-01-01
Students' creativity is an important focus globally and is interrelated with students' spatial abilities. Additionally, three-dimensional computer-assisted drawing (3D-CAD) overcomes barriers to spatial expression during the creative design process. Does 3D-CAD affect students' creative abilities? The purpose of this study was to…
Moreno, Karen; Wroe, Stephen; Clausen, Philip; McHenry, Colin; D’Amore, Domenic C; Rayfield, Emily J; Cunningham, Eleanor
2008-01-01
The Komodo dragon (Varanus komodoensis) displays a unique hold and pull-feeding technique. Its delicate ‘space-frame’ skull morphology differs greatly from that apparent in most living large prey specialists and is suggestive of a high degree of optimization, wherein use of materials is minimized. Here, using high-resolution finite element modelling based on dissection and in vivo bite and pull data, we present results detailing the mechanical performance of the giant lizard's skull. Unlike most modern predators, V. komodoensis applies minimal input from the jaw muscles when butchering prey. Instead it uses series of actions controlled by postcranial muscles. A particularly interesting feature of the performance of the skull is that it reveals considerably lower overall stress when these additional extrinsic forces are added to those of the jaw adductors. This remarkable reduction in stress in response to additional force is facilitated by both internal and external bone anatomy. Functional correlations obtained from these analyses also provide a solid basis for the interpretation of feeding ecology in extinct species, including dinosaurs and sabre-tooth cats, with which V. komodoensis shares various cranial and dental characteristics. PMID:18510503
NASA Astrophysics Data System (ADS)
Wu, Guangxi; Yu, Xiong
2015-06-01
Thermoelectric power generator has potential for small-scale and distributed power generation because of its high durability and scalability. It is very important to realize that the transient behavior of thermoelectric modules (TEM) affects a thermoelectric generator's response to dynamic working environments. Traditionally, researchers have used simplified models to describe the behavior of thermoelectric modules. In this paper we propose a comprehensive mathematical model that considers the effect of variations of chemical potential and carrier density, which are ignored by traditional models. Finite element models based on this new model are used to simulate the transient behavior of a thermoelectric module subjected to rapid changes in boundary temperature or working load. Simulation results show that transition times of thermoelectric modules affected by temperature change are much longer than those of modules affected by changes in electrical load resistance. Sudden changes in working temperature cause voltage overshoot of the TEM output, which, however, is not observed in responses to sudden changes of load resistance. Comparisons also show there are significant differences between the behavior of TEM predicted by use of this new comprehensive model and that predicted by use of traditional models, particularly for the high-temperature intrinsic ionization region and the low-temperature weak ionization region. This implies that chemical potential and carrier density variations, which are taken into account by this new model but ignored by traditional models, have major effects on the performance of TEM.
Moreno, Karen; Wroe, Stephen; Clausen, Philip; McHenry, Colin; D'Amore, Domenic C; Rayfield, Emily J; Cunningham, Eleanor
2008-06-01
The Komodo dragon (Varanus komodoensis) displays a unique hold and pull-feeding technique. Its delicate 'space-frame' skull morphology differs greatly from that apparent in most living large prey specialists and is suggestive of a high degree of optimization, wherein use of materials is minimized. Here, using high-resolution finite element modelling based on dissection and in vivo bite and pull data, we present results detailing the mechanical performance of the giant lizard's skull. Unlike most modern predators, V. komodoensis applies minimal input from the jaw muscles when butchering prey. Instead it uses series of actions controlled by postcranial muscles. A particularly interesting feature of the performance of the skull is that it reveals considerably lower overall stress when these additional extrinsic forces are added to those of the jaw adductors. This remarkable reduction in stress in response to additional force is facilitated by both internal and external bone anatomy. Functional correlations obtained from these analyses also provide a solid basis for the interpretation of feeding ecology in extinct species, including dinosaurs and sabre-tooth cats, with which V. komodoensis shares various cranial and dental characteristics. PMID:18510503
Seismic imaging using finite-differences and parallel computers
Ober, C.C.
1997-12-31
A key to reducing the risks and costs of associated with oil and gas exploration is the fast, accurate imaging of complex geologies, such as salt domes in the Gulf of Mexico and overthrust regions in US onshore regions. Prestack depth migration generally yields the most accurate images, and one approach to this is to solve the scalar wave equation using finite differences. As part of an ongoing ACTI project funded by the US Department of Energy, a finite difference, 3-D prestack, depth migration code has been developed. The goal of this work is to demonstrate that massively parallel computers can be used efficiently for seismic imaging, and that sufficient computing power exists (or soon will exist) to make finite difference, prestack, depth migration practical for oil and gas exploration. Several problems had to be addressed to get an efficient code for the Intel Paragon. These include efficient I/O, efficient parallel tridiagonal solves, and high single-node performance. Furthermore, to provide portable code the author has been restricted to the use of high-level programming languages (C and Fortran) and interprocessor communications using MPI. He has been using the SUNMOS operating system, which has affected many of his programming decisions. He will present images created from two verification datasets (the Marmousi Model and the SEG/EAEG 3D Salt Model). Also, he will show recent images from real datasets, and point out locations of improved imaging. Finally, he will discuss areas of current research which will hopefully improve the image quality and reduce computational costs.
Magnetization reversal in 3D nano-structures of different shapes
NASA Astrophysics Data System (ADS)
Blachowicz, T.; Ehrmann, A.
2016-04-01
Magnetic nano-particles have been intensively studied during the last decade due to their potential utilization in various applications. An important topic is the dependence of magnetic properties on the exact samples shape. After demonstrating the influence of shape distortions in magnetic nano half-spheres on magnetization reversal processes and hysteresis shapes, a series of different 3D nano-objects from permalloy with shape modifications has been examined with respect to their magnetic properties. Modifications are performed by cutting parts of diverse samples between the extrema of a cuboid and a half-sphere. Simulations of these samples have been performed by Magpar, using external magnetic fields along two different axes, swept with two different speeds. Depending on the original particle shape and its modifications, several phenomena can be found: Cutting a hole in a cuboid can switch the hard axis from out-of-plane to the in-plane direction. In some nano-particles, strong oscillations occur which can be suppressed by appropriate shape modifications. In some of the nano-objects, the magnetization reversal mechanism is completely altered by a change in the field sweeping speed. The article gives an overview of the different possibilities to tailor magnetic properties of nano-systems.
NASA Astrophysics Data System (ADS)
Luo, Jiao; Wu, Bin; Li, Miao-Quan
2012-02-01
The physically-based internal state variable (ISV) models were used to describe the changes of dislocation density, grain size, and flow stress in the high temperature deformation of titanium alloys in this study. The constants of the present models could be identified based on experimental results, which were conducted at deformation temperatures ranging from 1093 K to 1303 K, height reductions ranging from 20% to 60%, and the strain rates of 0.001, 0.01, 0.1, 1.0, and 10.0 s-1. The physically-based internal state variable models were implemented into the commercial finite element (FE) code. Then, a three-dimensional (3D) FE simulation system coupling of deformation, heat transfer, and microstructure evolution was developed for the blade forging of Ti-6Al-4V alloy. FE analysis was carried out to simulate the microstructure evolution in the blade forging of Ti-6Al-4V alloy. Finally, the blade forging tests of Ti-6Al-4V alloy were performed to validate the results of FE simulation. According to the tensile tests, it is seen that the mechanical properties, such as tensile strength and elongation, satisfy the application requirements well. The maximum and minimum differences between the calculated and experimental grain size of primary α phase are 11.71% and 4.23%, respectively. Thus, the industrial trials show a good agreement with FE simulation of blade forging.
3D reconstruction for sinusoidal motion based on different feature detection algorithms
NASA Astrophysics Data System (ADS)
Zhang, Peng; Zhang, Jin; Deng, Huaxia; Yu, Liandong
2015-02-01
The dynamic testing of structures and components is an important area of research. Extensive researches on the methods of using sensors for vibration parameters have been studied for years. With the rapid development of industrial high-speed camera and computer hardware, the method of using stereo vision for dynamic testing has been the focus of the research since the advantages of non-contact, full-field, high resolution and high accuracy. But in the country there is not much research about the dynamic testing based on stereo vision, and yet few people publish articles about the three-dimensional (3D) reconstruction of feature points in the case of dynamic. It is essential to the following analysis whether it can obtain accurate movement of target objects. In this paper, an object with sinusoidal motion is detected by stereo vision and the accuracy with different feature detection algorithms is investigated. Three different marks including dot, square and circle are stuck on the object and the object is doing sinusoidal motion by vibration table. Then use feature detection algorithm speed-up robust feature (SURF) to detect point, detect square corners by Harris and position the center by Hough transform. After obtaining the pixel coordinate values of the feature point, the stereo calibration parameters are used to achieve three-dimensional reconstruction through triangulation principle. The trajectories of the specific direction according to the vibration frequency and the frequency camera acquisition are obtained. At last, the reconstruction accuracy of different feature detection algorithms is compared.
Differences in 3D kinematics between volleyball and beach volleyball spike movements.
Tilp, Markus; Wagner, Herbert; Müller, Erich
2008-09-01
The purpose of this study was to identify differences between volleyball and beach volleyball spike jump movements performed on an indoor surface and sand surface respectively. Eight elite male volleyball players performed spike jump movements on both surfaces. An eight-camera motion capturing system (250Hz) was used to generate 3D kinematic data. Seven groups of variables representing the kinematics of the centre of mass, the countermovement, the approach phase, and the angular amplitudes and maximal velocities of the lower and upper limbs were examined using Hotelling's T2(2). Significant differences were observed in the movement of the centre of mass (P < 0.05), the countermovement, the kinematics of the approach phase, and the angular amplitudes of the lower limbs. However, no significant differences were observed either in the maximal angular velocities of the lower and upper limbs, or in the amplitudes of the upper limb motion. In conclusion, the participants showed significant adaptation to changed movement conditions. As a result of the compliance of the sand surface, the participants slowed down their movements, especially during the phase of transition from knee flexion to extension and during the extension phase. Furthermore, the participants demonstrated changes in foot position to reach the greatest height possible. PMID:18972887
On the wavelet optimized finite difference method
NASA Technical Reports Server (NTRS)
Jameson, Leland
1994-01-01
When one considers the effect in the physical space, Daubechies-based wavelet methods are equivalent to finite difference methods with grid refinement in regions of the domain where small scale structure exists. Adding a wavelet basis function at a given scale and location where one has a correspondingly large wavelet coefficient is, essentially, equivalent to adding a grid point, or two, at the same location and at a grid density which corresponds to the wavelet scale. This paper introduces a wavelet optimized finite difference method which is equivalent to a wavelet method in its multiresolution approach but which does not suffer from difficulties with nonlinear terms and boundary conditions, since all calculations are done in the physical space. With this method one can obtain an arbitrarily good approximation to a conservative difference method for solving nonlinear conservation laws.
FEMFLOW3D; a finite-element program for the simulation of three-dimensional aquifers; version 1.0
Durbin, Timothy J.; Bond, Linda D.
1998-01-01
This document also includes model validation, source code, and example input and output files. Model validation was performed using four test problems. For each test problem, the results of a model simulation with FEMFLOW3D were compared with either an analytic solution or the results of an independent numerical approach. The source code, written in the ANSI x3.9-1978 FORTRAN standard, and the complete input and output of an example problem are listed in the appendixes.
Pressure-selective modulation of NMDA receptor subtypes may reflect 3D structural differences.
Mor, Amir; Kuttner, Yosef Y; Levy, Shiri; Mor, Merav; Hollmann, Michael; Grossman, Yoram
2012-01-01
Professional deep-water divers exposed to high pressure (HP) above 1.1 MPa suffer from High Pressure Neurological Syndrome (HPNS), which is associated with CNS hyperexcitability. We have previously reported that HP augments N-methyl-D-aspartate receptor (NMDAR) synaptic responses, increases neuronal excitability, and potentially causes irreversible neuronal damage. We now report that HP (10.1 MPa) differentially affects eight specific NMDAR subtypes. GluN1(1a or 1b) was co-expressed with one of the four GluN2(A-D) subunits in Xenopus laevis oocytes. HP increased ionic currents (measured by two electrode voltage clamps) of one subtype, reduced the current in four others, and did not affect the current in the remaining three. 3D theoretical modeling was aimed at revealing specific receptor domains involved with HP selectivity. In light of the information on the CNS spatial distribution of the different NMDAR subtypes, we conclude that the NMDAR's diverse responses to HP may lead to selective HP effects on different brain regions. These discoveries call for further and more specific investigation of deleterious HP effects and suggest the need for a re-evaluation of deep-diving safety guidelines. PMID:22973194
NASA Technical Reports Server (NTRS)
Ponchak, G. E.; Bushyager, N.; Papapolymerou, J.; Tentzeris, E. M.; Laskar, J.
2002-01-01
Three-dimensional circuits built upon multiple layers of polyimide are required for constructing Si/SiGe monolithic microwave/mm-wave integrated circuits on CMOS (low resistivity) Si wafers. It is expected that these circuits will replace the ones fabricated on GaAs and reduce the overall system cost. However, the closely spaced transmission lines that are required for a high-density circuit environment are susceptible to high levels of cross-coupling, which degrades the overall circuit performance. In this paper, theoretical and experimental results on coupling and ways to reduce it are presented for two types of transmission lines: a) the microstrip line and b) the Finite Ground Coplanar (FGC) line. For microstrip lines it is shown that a fence of metalized via-holes can significantly reduce coupling, especially in the case when both lines are on the same polyimide layer or when the shielding structure extends through several polyimide layers. For closely spaced microstrip lines, coupling is lower for a metal filled trench shield than a via-hole fence. Coupling amongst microstrip lines is dependent on the ratio of line separation to polyimide thickness and is primarily due to magnetic fields. For FGC lines it is shown that they have in general low coupling that can be reduced significantly when there is even a small gap between the ground planes of each line. FGC lines have approximately 8 dB lower coupling than coupled coplanar waveguides (CPW). In addition, forward and backward characteristics of the FGC lines do not resemble those of other transmission lines such as microstrip. Therefore, the coupling mechanism of the FGC lines is different compared to thin film microstrip lines.
Puso, M; Maker, B N; Ferencz, R M; Hallquist, J O
2000-03-24
This report provides the NIKE3D user's manual update summary for changes made from version 3.0.0 April 24, 1995 to version 3.3.6 March 24,2000. The updates are excerpted directly from the code printed output file (hence the Courier font and formatting), are presented in chronological order and delineated by NIKE3D version number. NIKE3D is a fully implicit three-dimensional finite element code for analyzing the finite strain static and dynamic response of inelastic solids, shells, and beams. Spatial discretization is achieved by the use of 8-node solid elements, 2-node truss and beam elements, and 4-node membrane and shell elements. Thirty constitutive models are available for representing a wide range of elastic, plastic, viscous, and thermally dependent material behavior. Contact-impact algorithms permit gaps, frictional sliding, and mesh discontinuities along material interfaces. Several nonlinear solution strategies are available, including Full-, Modified-, and Quasi-Newton methods. The resulting system of simultaneous linear equations is either solved iteratively by an element-by-element method, or directly by a direct factorization method.