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)
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.
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.
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.
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.
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 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.
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.
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.
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
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.
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.
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.
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.
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.
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.
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
Software suite for finite difference method models.
Arola, T; Hannula, M; Narra, N; Malmivuo, J; Hyttinen, J
2006-01-01
We have developed a software suite for finite difference method (FDM) model construction, visualization and quasi-static simulation to be used in bioelectric field modeling. The aim of the software is to provide a full path from medical image data to simulation of bioelectric phenomena and results visualization. It is written in Java and can be run on various platforms while still supporting all features included. The software can be distributed across a network utilizing dedicated servers for calculation intensive tasks. Supported visualization modes are both two- and three-dimensional modes. PMID:17946057
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.
Efficient discretization in finite difference method
NASA Astrophysics Data System (ADS)
Rozos, Evangelos; Koussis, Antonis; Koutsoyiannis, Demetris
2015-04-01
Finite difference method (FDM) is a plausible and simple method for solving partial differential equations. The standard practice is to use an orthogonal discretization to form algebraic approximate formulations of the derivatives of the unknown function and a grid, much like raster maps, to represent the properties of the function domain. For example, for the solution of the groundwater flow equation, a raster map is required for the characterization of the discretization cells (flow cell, no-flow cell, boundary cell, etc.), and two raster maps are required for the hydraulic conductivity and the storage coefficient. Unfortunately, this simple approach to describe the topology comes along with the known disadvantages of the FDM (rough representation of the geometry of the boundaries, wasted computational resources in the unavoidable expansion of the grid refinement in all cells of the same column and row, etc.). To overcome these disadvantages, Hunt has suggested an alternative approach to describe the topology, the use of an array of neighbours. This limits the need for discretization nodes only for the representation of the boundary conditions and the flow domain. Furthermore, the geometry of the boundaries is described more accurately using a vector representation. Most importantly, graded meshes can be employed, which are capable of restricting grid refinement only in the areas of interest (e.g. regions where hydraulic head varies rapidly, locations of pumping wells, etc.). In this study, we test the Hunt approach against MODFLOW, a well established finite difference model, and the Finite Volume Method with Simplified Integration (FVMSI). The results of this comparison are examined and critically discussed.
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.
Comparison of finite-difference and analytic microwave calculation methods
Friedlander, F.I.; Jackson, H.W.; Barmatz, M.; Wagner, P.
1996-12-31
Normal modes and power absorption distributions in microwave cavities containing lossy dielectric samples were calculated for problems of interest in materials processing. The calculations were performed both using a commercially available finite-difference electromagnetic solver and by numerical evaluation of exact analytic expressions. Results obtained by the two methods applied to identical physical situations were compared. The studies validate the accuracy of the finite-difference electromagnetic solver. Relative advantages of the analytic and finite-difference methods are discussed.
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.
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
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.
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.
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.
Compact finite difference method for American option pricing
NASA Astrophysics Data System (ADS)
Zhao, Jichao; Davison, Matt; Corless, Robert M.
2007-09-01
A compact finite difference method is designed to obtain quick and accurate solutions to partial differential equation problems. The problem of pricing an American option can be cast as a partial differential equation. Using the compact finite difference method this problem can be recast as an ordinary differential equation initial value problem. The complicating factor for American options is the existence of an optimal exercise boundary which is jointly determined with the value of the option. In this article we develop three ways of combining compact finite difference methods for American option price on a single asset with methods for dealing with this optimal exercise boundary. Compact finite difference method one uses the implicit condition that solutions of the transformed partial differential equation be nonnegative to detect the optimal exercise value. This method is very fast and accurate even when the spatial step size h is large (h[greater-or-equal, slanted]0.1). Compact difference method two must solve an algebraic nonlinear equation obtained by Pantazopoulos (1998) at every time step. This method can obtain second order accuracy for space x and requires a moderate amount of time comparable with that required by the Crank Nicolson projected successive over relaxation method. Compact finite difference method three refines the free boundary value by a method developed by Barone-Adesi and Lugano [The saga of the American put, 2003], and this method can obtain high accuracy for space x. The last two of these three methods are convergent, moreover all the three methods work for both short term and long term options. Through comparison with existing popular methods by numerical experiments, our work shows that compact finite difference methods provide an exciting new tool for American option pricing.
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.
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
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 difference methods for the solution of unsteady potential flows
NASA Technical Reports Server (NTRS)
Caradonna, F. X.
1982-01-01
Various problems which are confronted in the development of an unsteady finite difference potential code are reviewed mainly in the context of what is done for a typical small disturbance and full potential method. The issues discussed include choice of equations, linearization and conservation, differencing schemes, and algorithm development. A number of applications, including unsteady three dimensional rotor calculations, are demonstrated.
Test of two methods for faulting on finite-difference calculations
Andrews, D.J.
1999-01-01
Tests of two fault boundary conditions show that each converges with second order accuracy as the finite-difference grid is refined. The first method uses split nodes so that there are disjoint grids that interact via surface traction. The 3D version described here is a generalization of a method I have used extensively in 2D; it is as accurate as the 2D version. The second method represents fault slip as inelastic strain in a fault zone. Offset of stress from its elastic value is seismic moment density. Implementation of this method is quite simple in a finite-difference scheme using velocity and stress as dependent variables.
Finite-difference lattice-Boltzmann methods for binary fluids.
Xu, Aiguo
2005-06-01
We investigate two-fluid Bhatnagar-Gross-Krook (BGK) kinetic methods for binary fluids. The developed theory works for asymmetric as well as symmetric systems. For symmetric systems it recovers Sirovich's theory and is summarized in models A and B. For asymmetric systems it contributes models C, D, and E which are especially useful when the total masses and/or local temperatures of the two components are greatly different. The kinetic models are discretized based on an octagonal discrete velocity model. The discrete-velocity kinetic models and the continuous ones are required to describe the same hydrodynamic equations. The combination of a discrete-velocity kinetic model and an appropriate finite-difference scheme composes a finite-difference lattice Boltzmann method. The validity of the formulated methods is verified by investigating (i) uniform relaxation processes, (ii) isothermal Couette flow, and (iii) diffusion behavior. PMID:16089910
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
Finite difference methods for the solution of unsteady potential flows
NASA Technical Reports Server (NTRS)
Caradonna, F. X.
1985-01-01
A brief review is presented of various problems which are confronted in the development of an unsteady finite difference potential code. This review is conducted mainly in the context of what is done for a typical small disturbance and full potential methods. The issues discussed include choice of equation, linearization and conservation, differencing schemes, and algorithm development. A number of applications including unsteady three-dimensional rotor calculation, are demonstrated.
High Order Finite Difference Methods for Multiscale Complex Compressible Flows
NASA Technical Reports Server (NTRS)
Sjoegreen, Bjoern; Yee, H. C.
2002-01-01
The classical way of analyzing finite difference schemes for hyperbolic problems is to investigate as many as possible of the following points: (1) Linear stability for constant coefficients; (2) Linear stability for variable coefficients; (3) Non-linear stability; and (4) Stability at discontinuities. We will build a new numerical method, which satisfies all types of stability, by dealing with each of the points above step by step.
Introduction to finite-difference methods for numerical fluid dynamics
Scannapieco, E.; Harlow, F.H.
1995-09-01
This work is intended to be a beginner`s exercise book for the study of basic finite-difference techniques in computational fluid dynamics. It is written for a student level ranging from high-school senior to university senior. Equations are derived from basic principles using algebra. Some discussion of partial-differential equations is included, but knowledge of calculus is not essential. The student is expected, however, to have some familiarity with the FORTRAN computer language, as the syntax of the computer codes themselves is not discussed. Topics examined in this work include: one-dimensional heat flow, one-dimensional compressible fluid flow, two-dimensional compressible fluid flow, and two-dimensional incompressible fluid flow with additions of the equations of heat flow and the {Kappa}-{epsilon} model for turbulence transport. Emphasis is placed on numerical instabilities and methods by which they can be avoided, techniques that can be used to evaluate the accuracy of finite-difference approximations, and the writing of the finite-difference codes themselves. Concepts introduced in this work include: flux and conservation, implicit and explicit methods, Lagrangian and Eulerian methods, shocks and rarefactions, donor-cell and cell-centered advective fluxes, compressible and incompressible fluids, the Boussinesq approximation for heat flow, Cartesian tensor notation, the Boussinesq approximation for the Reynolds stress tensor, and the modeling of transport equations. A glossary is provided which defines these and other terms.
Explicit Finite Difference Methods for the Delay Pseudoparabolic Equations
Amirali, I.; Amiraliyev, G. M.; Cakir, M.; Cimen, E.
2014-01-01
Finite difference technique is applied to numerical solution of the initial-boundary value problem for the semilinear delay Sobolev or pseudoparabolic equation. By the method of integral identities two-level difference scheme is constructed. For the time integration the implicit rule is being used. Based on the method of energy estimates the fully discrete scheme is shown to be absolutely stable and convergent of order two in space and of order one in time. The error estimates are obtained in the discrete norm. Some numerical results confirming the expected behavior of the method are shown. PMID:24688392
Explicit finite difference methods for the delay pseudoparabolic equations.
Amirali, I; Amiraliyev, G M; Cakir, M; Cimen, E
2014-01-01
Finite difference technique is applied to numerical solution of the initial-boundary value problem for the semilinear delay Sobolev or pseudoparabolic equation. By the method of integral identities two-level difference scheme is constructed. For the time integration the implicit rule is being used. Based on the method of energy estimates the fully discrete scheme is shown to be absolutely stable and convergent of order two in space and of order one in time. The error estimates are obtained in the discrete norm. Some numerical results confirming the expected behavior of the method are shown. PMID:24688392
Macroscopic traffic modeling with the finite difference method
Mughabghab, S.; Azarm, A.; Stock, D.
1996-03-15
A traffic congestion forecasting model (ATOP), developed in the present investigation, is described briefly. Several macroscopic models, based on the solution of the partial differential equation of conservation of vehicles by the finite difference method, were tested using actual traffic data. The functional form, as well as the parameters, of the equation of state which describes the relation between traffic speed and traffic density, were determined for a section of the Long Island Expressway. The Lax method and the forward difference technique were applied. The results of extensive tests showed that the Lax method, in addition to giving very good agreement with the traffic data, produces stable solutions.
Arrayed waveguide grating using the finite difference beam propagation method
NASA Astrophysics Data System (ADS)
Toledo, M. C. F.; Alayo, M. I.
2013-03-01
The purpose of this work is to analyze by simulation the coupling effects occurring in Arrayed Waveguide Grating (AWG) using the finite difference beam propagation method (FD-BPM). Conventional FD-BPM techniques do not immediately lend themselves to the analysis of large structures such as AWG. Cooper et al.1 introduced a description of the coupling between the interface of arrayed waveguides and star couplers using the numerically-assisted coupled-mode theory. However, when the arrayed waveguides are spatially close, such that, there is strong coupling between them, and coupled-mode theory is not adequate. On the other hand, Payne2 developed an exact eigenvalue equation for the super modes of a straight arrayed waveguide which involve a computational overhead. In this work, an integration of both methods is accomplished in order to describe the behavior of the propagation of light in guided curves. This new method is expected to reduce the necessary effort for simulation while also enabling the simulation of large and curved arrayed waveguides using a fully vectorial finite difference technique.
A finite-difference method for transonic airfoil design.
NASA Technical Reports Server (NTRS)
Steger, J. L.; Klineberg, J. M.
1972-01-01
This paper describes an inverse method for designing transonic airfoil sections or for modifying existing profiles. Mixed finite-difference procedures are applied to the equations of transonic small disturbance theory to determine the airfoil shape corresponding to a given surface pressure distribution. The equations are solved for the velocity components in the physical domain and flows with embedded shock waves can be calculated. To facilitate airfoil design, the method allows alternating between inverse and direct calculations to obtain a profile shape that satisfies given geometric constraints. Examples are shown of the application of the technique to improve the performance of several lifting airfoil sections. The extension of the method to three dimensions for designing supercritical wings is also indicated.
Accurate finite difference methods for time-harmonic wave propagation
NASA Technical Reports Server (NTRS)
Harari, Isaac; Turkel, Eli
1994-01-01
Finite difference methods for solving problems of time-harmonic acoustics are developed and analyzed. Multidimensional inhomogeneous problems with variable, possibly discontinuous, coefficients are considered, accounting for the effects of employing nonuniform grids. A weighted-average representation is less sensitive to transition in wave resolution (due to variable wave numbers or nonuniform grids) than the standard pointwise representation. Further enhancement in method performance is obtained by basing the stencils on generalizations of Pade approximation, or generalized definitions of the derivative, reducing spurious dispersion, anisotropy and reflection, and by improving the representation of source terms. The resulting schemes have fourth-order accurate local truncation error on uniform grids and third order in the nonuniform case. Guidelines for discretization pertaining to grid orientation and resolution are presented.
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
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
Elastic finite-difference method for irregular grids
Oprsal, I.; Zahradnik, J.
1999-01-01
Finite-difference (FD) modeling of complicated structures requires simple algorithms. This paper presents a new elastic FD method for spatially irregular grids that is simple and, at the same time, saves considerable memory and computing time. Features like faults, low-velocity layers, cavities, and/or nonplanar surfaces are treated on a fine grid, while the remaining parts of the model are, with equal accuracy, represented on a coarse grid. No interpolation is needed between the fine and coarse parts due to the rectangular grid cells. Relatively abrupt transitions between the small and large grid steps produce no numerical artifacts in the present method. Planar or nonplanar free surfaces, including underground cavities, are treated in a way similar to internal grid points but with consideration of the zero-valued elastic parameters and density outside the free surface (vacuum formalism). A theoretical proof that vacuum formalism fulfills the free-surface conditions is given. Numerical validation is performed through comparison with independent methods, comparing FD with explicitly prescribed boundary conditions and finite elements. Memory and computing time needed in the studied models was only about 10 to 40% of that employing regular square grids of equal accuracy. A practical example of a synthetic seismic section, showing clear signatures of a coal seam and cavity, is presented. The method can be extended to three dimensions.
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
A finite-difference contrast source inversion method
NASA Astrophysics Data System (ADS)
Abubakar, A.; Hu, W.; van den Berg, P. M.; Habashy, T. M.
2008-12-01
We present a contrast source inversion (CSI) algorithm using a finite-difference (FD) approach as its backbone for reconstructing the unknown material properties of inhomogeneous objects embedded in a known inhomogeneous background medium. Unlike the CSI method using the integral equation (IE) approach, the FD-CSI method can readily employ an arbitrary inhomogeneous medium as its background. The ability to use an inhomogeneous background medium has made this algorithm very suitable to be used in through-wall imaging and time-lapse inversion applications. Similar to the IE-CSI algorithm the unknown contrast sources and contrast function are updated alternately to reconstruct the unknown objects without requiring the solution of the full forward problem at each iteration step in the optimization process. The FD solver is formulated in the frequency domain and it is equipped with a perfectly matched layer (PML) absorbing boundary condition. The FD operator used in the FD-CSI method is only dependent on the background medium and the frequency of operation, thus it does not change throughout the inversion process. Therefore, at least for the two-dimensional (2D) configurations, where the size of the stiffness matrix is manageable, the FD stiffness matrix can be inverted using a non-iterative inversion matrix approach such as a Gauss elimination method for the sparse matrix. In this case, an LU decomposition needs to be done only once and can then be reused for multiple source positions and in successive iterations of the inversion. Numerical experiments show that this FD-CSI algorithm has an excellent performance for inverting inhomogeneous objects embedded in an inhomogeneous background medium.
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
NASA Technical Reports Server (NTRS)
Liu, C.; Liu, Z.
1993-01-01
The high order finite difference and multigrid methods have been successfully applied to direct numerical simulation (DNS) for flow transition in 3D channels and 3D boundary layers with 2D and 3D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semicoarsening multigrid method associated with line distributive relaxation scheme, and a new treatment of the outflow boundary condition, which needs only a very short buffer domain to damp all wave reflection, are developed. These approaches make the multigrid DNS code very accurate and efficient. This makes us not only able to do spatial DNS for the 3D channel and flat plate at low computational costs, but also able to do spatial DNS for transition in the 3D boundary layer with 3D single and multiple roughness elements. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments.
Optimization of a finite difference method for nonlinear wave equations
NASA Astrophysics Data System (ADS)
Chen, Miaochao
2013-07-01
Wave equations have important fluid dynamics background, which are extensively used in many fields, such as aviation, meteorology, maritime, water conservancy, etc. This paper is devoted to the explicit difference method for nonlinear wave equations. Firstly, a three-level and explicit difference scheme is derived. It is shown that the explicit difference scheme is uniquely solvable and convergent. Moreover, a numerical experiment is conducted to illustrate the theoretical results of the presented method.
Kim, S.
1994-12-31
Parallel iterative procedures based on domain decomposition techniques are defined and analyzed for the numerical solution of wave propagation by finite element and finite difference methods. For finite element methods, in a Lagrangian framework, an efficient way for choosing the algorithm parameter as well as the algorithm convergence are indicated. Some heuristic arguments for finding the algorithm parameter for finite difference schemes are addressed. Numerical results are presented to indicate the effectiveness of the methods.
A non-linear constrained optimization technique for the mimetic finite difference method
Manzini, Gianmarco; Svyatskiy, Daniil; Bertolazzi, Enrico; Frego, Marco
2014-09-30
This is a strategy for the construction of monotone schemes in the framework of the mimetic finite difference method for the approximation of diffusion problems on unstructured polygonal and polyhedral meshes.
NASA Technical Reports Server (NTRS)
Fix, G. J.; Rose, M. E.
1983-01-01
A least squares formulation of the system divu = rho, curlu = zeta is surveyed from the viewpoint of both finite element and finite difference methods. Closely related arguments are shown to establish convergence estimates.
NASA Astrophysics Data System (ADS)
Beilina, Larisa
2016-08-01
We present domain decomposition finite element/finite difference method for the solution of hyperbolic equation. The domain decomposition is performed such that finite elements and finite differences are used in different subdomains of the computational domain: finite difference method is used on the structured part of the computational domain and finite elements on the unstructured part of the domain. Explicit discretizations for both methods are constructed such that the finite element and the finite difference schemes coincide on the common structured overlapping layer between computational subdomains. Then the resulting approach can be considered as a pure finite element scheme which avoids instabilities at the interfaces. We derive an energy estimate for the underlying hyperbolic equation with absorbing boundary conditions and illustrate efficiency of the domain decomposition method on the reconstruction of the conductivity function in three dimensions.
NASA Astrophysics Data System (ADS)
Lisitsa, Vadim; Tcheverda, Vladimir; Botter, Charlotte
2016-04-01
We present an algorithm for the numerical simulation of seismic wave propagation in models with a complex near surface part and free surface topography. The approach is based on the combination of finite differences with the discontinuous Galerkin method. The discontinuous Galerkin method can be used on polyhedral meshes; thus, it is easy to handle the complex surfaces in the models. However, this approach is computationally intense in comparison with finite differences. Finite differences are computationally efficient, but in general, they require rectangular grids, leading to the stair-step approximation of the interfaces, which causes strong diffraction of the wavefield. In this research we present a hybrid algorithm where the discontinuous Galerkin method is used in a relatively small upper part of the model and finite differences are applied to the main part of the model.
NASA Astrophysics Data System (ADS)
Hallez, Hans; Vanrumste, Bart; Van Hese, Peter; D'Asseler, Yves; Lemahieu, Ignace; Van de Walle, Rik
2005-08-01
Many implementations of electroencephalogram (EEG) dipole source localization neglect the anisotropical conductivities inherent to brain tissues, such as the skull and white matter anisotropy. An examination of dipole localization errors is made in EEG source analysis, due to not incorporating the anisotropic properties of the conductivity of the skull and white matter. First, simulations were performed in a 5 shell spherical head model using the analytical formula. Test dipoles were placed in three orthogonal planes in the spherical head model. Neglecting the skull anisotropy results in a dipole localization error of, on average, 13.73 mm with a maximum of 24.51 mm. For white matter anisotropy these values are 11.21 mm and 26.3 mm, respectively. Next, a finite difference method (FDM), presented by Saleheen and Kwong (1997 IEEE Trans. Biomed. Eng. 44 800-9), is used to incorporate the anisotropy of the skull and white matter. The FDM method has been validated for EEG dipole source localization in head models with all compartments isotropic as well as in a head model with white matter anisotropy. In a head model with skull anisotropy the numerical method could only be validated if the 3D lattice was chosen very fine (grid size <=2 mm).
NASA Astrophysics Data System (ADS)
Sun, Y.; Zhang, W.; Chen, X.
2014-12-01
This paper presents a curvilinear grid finite difference method for modeling seismic wave propagation with topographic fluid (acoustic) and solid (elastic) interface. The curvilinear grid finite difference method has been successfully used for seismic wave simulation with free surface topography and earthquake dynamics with complex falut geometry. For seismic wave simulation with topographic sea bottom, we use the curvilinear grid to conform the grid to the sea bottom to avoid artifical scatterings due to stair-case approximation. We solve the acoustic wave equation in the water layer and the elastic wave equation in the solid below the sea bottom. The fluid-solid interface condition is implemented by decomposing velocity and stress components to normal and parallel directions of the sea bottom. The results exhibit high accuracy by comparsion with analytical solutions for flat interfaces and also work very well when the fluid-solid interface is topographic. The scheme can be easily extended to 3-D situation.
NASA Astrophysics Data System (ADS)
Choi, S. J.; Kim, J.; Shin, S.
2014-12-01
In this presentation, a new non-hydrostatic (NH) dynamical core using the spectral element method (SEM) in the horizontal discretization and the finite difference method (FDM) in the vertical discretization will be presented. By using horizontal SEM, which decomposes the physical domain into smaller pieces with a small communication stencil, we can achieve a high level of scalability. Also by using vertical FDM, we provide an easy way for coupling the dynamics and existing physics packages. The Euler equations used here are in a flux form based on the hybrid sigma hydrostatic pressure vertical coordinate, which are similar to those used in the Weather Research and Forecasting (WRF) model. Within these Euler equations, we use a time-split third-order Runge-Kutta (RK3) for the time discretization. In order to establish robustness, firstly the NH dynamical core is verified in a simplified two dimensional (2D) slice framework by conducting widely used standard benchmark tests, and then we verify the global three dimensional (3D) dynamical core on the cubed-sphere grid with several test cases introduced by Dynamical Core Model Intercomparison Project (DCMIP).
Kauppinen, P; Hyttinen, J; Laarne, P; Malmivuo, J
1999-02-01
There is an evolving need for new information available by employing patient tailored anatomically accurate computer models of the electrical properties of the human body. Because construction of a computer model can be difficult and laborious to perform sufficiently well, devised models have varied greatly in the level of anatomical accuracy incorporated in them. This has restricted the validity of conducted simulations. In the present study, a versatile software package was developed to transform anatomic voxel data into accurate finite difference method volume conductor models conveniently and in a short time. The package includes components for model construction, simulation, visualisation and detailed analysis of simulation output based on volume conductor theory. Due to the methods developed, models can comprise more anatomical details than the prior computer models. Several models have been constructed, for example, a highly detailed 3-D anatomically accurate computer model of the human thorax as a volume conductor utilising the US National Library of Medicine's (NLM) Visible Human Man (VHM) digital anatomy data. Based on the validation runs the developed software package is readily applicable in analysis of a wide range of bioelectric field problems. PMID:10092033
SOME NEW FINITE DIFFERENCE METHODS FOR HELMHOLTZ EQUATIONS ON IRREGULAR DOMAINS OR WITH INTERFACES.
Wan, Xiaohai; Li, Zhilin
2012-06-01
Solving a Helmholtz equation Δu + λu = f efficiently is a challenge for many applications. For example, the core part of many efficient solvers for the incompressible Navier-Stokes equations is to solve one or several Helmholtz equations. In this paper, two new finite difference methods are proposed for solving Helmholtz equations on irregular domains, or with interfaces. For Helmholtz equations on irregular domains, the accuracy of the numerical solution obtained using the existing augmented immersed interface method (AIIM) may deteriorate when the magnitude of λ is large. In our new method, we use a level set function to extend the source term and the PDE to a larger domain before we apply the AIIM. For Helmholtz equations with interfaces, a new maximum principle preserving finite difference method is developed. The new method still uses the standard five-point stencil with modifications of the finite difference scheme at irregular grid points. The resulting coefficient matrix of the linear system of finite difference equations satisfies the sign property of the discrete maximum principle and can be solved efficiently using a multigrid solver. The finite difference method is also extended to handle temporal discretized equations where the solution coefficient λ is inversely proportional to the mesh size. PMID:22701346
A FINITE-DIFFERENCE, DISCRETE-WAVENUMBER METHOD FOR CALCULATING RADAR TRACES
A hybrid of the finite-difference method and the discrete-wavenumber method is developed to calculate radar traces. The method is based on a three-dimensional model defined in the Cartesian coordinate system; the electromagnetic properties of the model are symmetric with respect ...
A FINITE-DIFFERENCE, DISCRETE-WAVENUMBER METHOD FOR CALCULATING RADAR TRACES
A hybrid of the finite-difference method and the discrete-wavenumber method is developed to calculate radar traces. The method is based on a three-dimensional model defined in the Cartesian coordinate system; the electromag-netic properties of the model are symmetric with respect...
Numerical solution of a diffusion problem by exponentially fitted finite difference methods.
D'Ambrosio, Raffaele; Paternoster, Beatrice
2014-01-01
This paper is focused on the accurate and efficient solution of partial differential differential equations modelling a diffusion problem by means of exponentially fitted finite difference numerical methods. After constructing and analysing special purpose finite differences for the approximation of second order partial derivatives, we employed them in the numerical solution of a diffusion equation with mixed boundary conditions. Numerical experiments reveal that a special purpose integration, both in space and in time, is more accurate and efficient than that gained by employing a general purpose solver. PMID:26034665
Finite-difference, spectral and Galerkin methods for time-dependent problems
NASA Technical Reports Server (NTRS)
Tadmor, E.
1983-01-01
Finite difference, spectral and Galerkin methods for the approximate solution of time dependent problems are surveyed. A unified discussion on their accuracy, stability and convergence is given. In particular, the dilemma of high accuracy versus stability is studied in some detail.
High Order Finite Difference Methods, Multidimensional Linear Problems and Curvilinear Coordinates
NASA Technical Reports Server (NTRS)
Nordstrom, Jan; Carpenter, Mark H.
1999-01-01
Boundary and interface conditions are derived for high order finite difference methods applied to multidimensional linear problems in curvilinear coordinates. The boundary and interface conditions lead to conservative schemes and strict and strong stability provided that certain metric conditions are met.
Relative and Absolute Error Control in a Finite-Difference Method Solution of Poisson's Equation
ERIC Educational Resources Information Center
Prentice, J. S. C.
2012-01-01
An algorithm for error control (absolute and relative) in the five-point finite-difference method applied to Poisson's equation is described. The algorithm is based on discretization of the domain of the problem by means of three rectilinear grids, each of different resolution. We discuss some hardware limitations associated with the algorithm,…
A coarse-mesh nodal method-diffusive-mesh finite difference method
Joo, H.; Nichols, W.R.
1994-05-01
Modern nodal methods have been successfully used for conventional light water reactor core analyses where the homogenized, node average cross sections (XSs) and the flux discontinuity factors (DFs) based on equivalence theory can reliably predict core behavior. For other types of cores and other geometries characterized by tightly-coupled, heterogeneous core configurations, the intranodal flux shapes obtained from a homogenized nodal problem may not accurately portray steep flux gradients near fuel assembly interfaces or various reactivity control elements. This may require extreme values of DFs (either very large, very small, or even negative) to achieve a desired solution accuracy. Extreme values of DFs, however, can disrupt the convergence of the iterative methods used to solve for the node average fluxes, and can lead to a difficulty in interpolating adjacent DF values. Several attempts to remedy the problem have been made, but nothing has been satisfactory. A new coarse-mesh nodal scheme called the Diffusive-Mesh Finite Difference (DMFD) technique, as contrasted with the coarse-mesh finite difference (CMFD) technique, has been developed to resolve this problem. This new technique and the development of a few-group, multidimensional kinetics computer program are described in this paper.
Finite difference methods for transient signal propagation in stratified dispersive media
NASA Technical Reports Server (NTRS)
Lam, D. H.
1975-01-01
Explicit difference equations are presented for the solution of a signal of arbitrary waveform propagating in an ohmic dielectric, a cold plasma, a Debye model dielectric, and a Lorentz model dielectric. These difference equations are derived from the governing time-dependent integro-differential equations for the electric fields by a finite difference method. A special difference equation is derived for the grid point at the boundary of two different media. Employing this difference equation, transient signal propagation in an inhomogeneous media can be solved provided that the medium is approximated in a step-wise fashion. The solutions are generated simply by marching on in time. It is concluded that while the classical transform methods will remain useful in certain cases, with the development of the finite difference methods described, an extensive class of problems of transient signal propagating in stratified dispersive media can be effectively solved by numerical methods.
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
Projection methods for incompressible flow problems with WENO finite difference schemes
NASA Astrophysics Data System (ADS)
de Frutos, Javier; John, Volker; Novo, Julia
2016-03-01
Weighted essentially non-oscillatory (WENO) finite difference schemes have been recommended in a competitive study of discretizations for scalar evolutionary convection-diffusion equations [20]. This paper explores the applicability of these schemes for the simulation of incompressible flows. To this end, WENO schemes are used in several non-incremental and incremental projection methods for the incompressible Navier-Stokes equations. Velocity and pressure are discretized on the same grid. A pressure stabilization Petrov-Galerkin (PSPG) type of stabilization is introduced in the incremental schemes to account for the violation of the discrete inf-sup condition. Algorithmic aspects of the proposed schemes are discussed. The schemes are studied on several examples with different features. It is shown that the WENO finite difference idea can be transferred to the simulation of incompressible flows. Some shortcomings of the methods, which are due to the splitting in projection schemes, become also obvious.
Transient analysis of printed lines using finite-difference time-domain method
Ahmed, Shahid
2012-03-29
Comprehensive studies of ultra-wideband pulses and electromagnetic coupling on printed coupled lines have been performed using full-wave 3D finite-difference time-domain analysis. Effects of unequal phase velocities of coupled modes, coupling between line traces, and the frequency dispersion on the waveform fidelity and crosstalk have been investigated in detail. To discriminate the contributions of different mechanisms into pulse evolution, single and coupled microstrip lines without (ϵ_{r} = 1) and with (ϵ_{r} > 1) dielectric substrates have been examined. To consistently compare the performance of the coupled lines with substrates of different permittivities and transients of different characteristic times, a generic metric similar to the electrical wavelength has been introduced. The features of pulse propagation on coupled lines with layered and pedestal substrates and on the irregular traces have been explored. Finally, physical interpretations of the simulation results are discussed in the paper.
A comparison of the finite difference and finite element methods for heat transfer calculations
NASA Technical Reports Server (NTRS)
Emery, A. F.; Mortazavi, H. R.
1982-01-01
The finite difference method and finite element method for heat transfer calculations are compared by describing their bases and their application to some common heat transfer problems. In general it is noted that neither method is clearly superior, and in many instances, the choice is quite arbitrary and depends more upon the codes available and upon the personal preference of the analyst than upon any well defined advantages of one method. Classes of problems for which one method or the other is better suited are defined.
NASA Astrophysics Data System (ADS)
Min, Xiaoyi
This thesis first presents the study of the interaction of electromagnetic waves with three-dimensional heterogeneous, dielectric, magnetic, and lossy bodies by surface integral equation modeling. Based on the equivalence principle, a set of coupled surface integral equations is formulated and then solved numerically by the method of moments. Triangular elements are used to model the interfaces of the heterogeneous body, and vector basis functions are defined to expand the unknown current in the formulation. The validity of this formulation is verified by applying it to concentric spheres for which an exact solution exists. The potential applications of this formulation to a partially coated sphere and a homogeneous human body are discussed. Next, this thesis also introduces an efficient new set of integral equations for treating the scattering problem of a perfectly conducting body coated with a thin magnetically lossy layer. These electric field integral equations and magnetic field integral equations are numerically solved by the method of moments (MoM). To validate the derived integral equations, an alternative method to solve the scattering problem of an infinite circular cylinder coated with a thin magnetic lossy layer has also been developed, based on the eigenmode expansion. Results for the radar cross section and current densities via the MoM and the eigenmode expansion method are compared. The agreement is excellent. The finite difference time domain method is subsequently implemented to solve a metallic object coated with a magnetic thin layer and numerical results are compared with that by the MoM. Finally, this thesis presents an application of the finite-difference time-domain approach to the problem of electromagnetic receiving and scattering by a cavity -backed antenna situated on an infinite conducting plane. This application involves modifications of Yee's model, which applies the difference approximations of field derivatives to differential
NASA Astrophysics Data System (ADS)
Lin, M. C.; Nieter, C.; Stoltz, P. H.; Smithe, D. N.
2009-05-01
This work introduces a conformal finite difference time domain (CFDTD) method to accurately determine the dispersion relation of an A6 relativistic magnetron. The accuracy is measured by comparing with accurate SUPERFISH calculations based on finite element method. The results show that an accuracy of 99.4% can be achieved by using only 10,000 mesh points with Dey-Mittra algorithm. By comparison, a mesh number of 360,000 is needed to preserve 99% accuracy using conventional FDTD method. This suggests one can efficiently and accurately study the hot tests of microwave tubes using CFDTD particle-in-cell method instead of conventional FDTD one.
Finite Differences and Collocation Methods for the Solution of the Two Dimensional Heat Equation
NASA Technical Reports Server (NTRS)
Kouatchou, Jules
1999-01-01
In this paper we combine finite difference approximations (for spatial derivatives) and collocation techniques (for the time component) to numerically solve the two dimensional heat equation. We employ respectively a second-order and a fourth-order schemes for the spatial derivatives and the discretization method gives rise to a linear system of equations. We show that the matrix of the system is non-singular. Numerical experiments carried out on serial computers, show the unconditional stability of the proposed method and the high accuracy achieved by the fourth-order scheme.
Finite difference methods with non-uniform meshes for nonlinear fractional differential equations
NASA Astrophysics Data System (ADS)
Li, Changpin; Yi, Qian; Chen, An
2016-07-01
In this article, finite difference methods with non-uniform meshes for solving nonlinear fractional differential equations are presented, where the non-equidistant stepsize is non-decreasing. The rectangle formula and trapezoid formula are proposed based on the non-uniform meshes. Combining the above two methods, we then establish the predictor-corrector scheme. The error and stability analysis are carefully investigated. At last, numerical examples are carried out to verify the theoretical analysis. Besides, the comparisons between non-uniform and uniform meshes are given, where the non-uniform meshes show the better performance when dealing with the less smooth problems.
An adaptive-mesh finite-difference solution method for the Navier-Stokes equations
NASA Astrophysics Data System (ADS)
Luchini, Paolo
1987-02-01
An adjustable variable-spacing grid is presented which permits the addition or deletion of single points during iterative solutions of the Navier-Stokes equations by finite difference methods. The grid is designed for application to two-dimensional steady-flow problems which can be described by partial differential equations whose second derivatives are constrained to the Laplacian operator. An explicit Navier-Stokes equations solution technique defined for use with the grid incorporates a hybrid form of the convective terms. Three methods are developed for automatic modifications of the mesh during calculations.
NASA Astrophysics Data System (ADS)
Kumar, Vivek; Raghurama Rao, S. V.
2008-04-01
Non-standard finite difference methods (NSFDM) introduced by Mickens [ Non-standard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994] are interesting alternatives to the traditional finite difference and finite volume methods. When applied to linear hyperbolic conservation laws, these methods reproduce exact solutions. In this paper, the NSFDM is first extended to hyperbolic systems of conservation laws, by a novel utilization of the decoupled equations using characteristic variables. In the second part of this paper, the NSFDM is studied for its efficacy in application to nonlinear scalar hyperbolic conservation laws. The original NSFDMs introduced by Mickens (1994) were not in conservation form, which is an important feature in capturing discontinuities at the right locations. Mickens [Construction and analysis of a non-standard finite difference scheme for the Burgers-Fisher equations, Journal of Sound and Vibration 257 (4) (2002) 791-797] recently introduced a NSFDM in conservative form. This method captures the shock waves exactly, without any numerical dissipation. In this paper, this algorithm is tested for the case of expansion waves with sonic points and is found to generate unphysical expansion shocks. As a remedy to this defect, we use the strategy of composite schemes [R. Liska, B. Wendroff, Composite schemes for conservation laws, SIAM Journal of Numerical Analysis 35 (6) (1998) 2250-2271] in which the accurate NSFDM is used as the basic scheme and localized relaxation NSFDM is used as the supporting scheme which acts like a filter. Relaxation schemes introduced by Jin and Xin [The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Communications in Pure and Applied Mathematics 48 (1995) 235-276] are based on relaxation systems which replace the nonlinear hyperbolic conservation laws by a semi-linear system with a stiff relaxation term. The relaxation parameter ( λ) is chosen locally
Serpentine: Finite Difference Methods for Wave Propagation in Second Order Formulation
Petersson, N A; Sjogreen, B
2012-03-26
second order system is significantly smaller. Another issue with re-writing a second order system into first order form is that compatibility conditions often must be imposed on the first order form. These (Saint-Venant) conditions ensure that the solution of the first order system also satisfies the original second order system. However, such conditions can be difficult to enforce on the discretized equations, without introducing additional modeling errors. This project has previously developed robust and memory efficient algorithms for wave propagation including effects of curved boundaries, heterogeneous isotropic, and viscoelastic materials. Partially supported by internal funding from Lawrence Livermore National Laboratory, many of these methods have been implemented in the open source software WPP, which is geared towards 3-D seismic wave propagation applications. This code has shown excellent scaling on up to 32,768 processors and has enabled seismic wave calculations with up to 26 Billion grid points. TheWPP calculations have resulted in several publications in the field of computational seismology, e.g.. All of our current methods are second order accurate in both space and time. The benefits of higher order accurate schemes for wave propagation have been known for a long time, but have mostly been developed for first order hyperbolic systems. For second order hyperbolic systems, it has not been known how to make finite difference schemes stable with free surface boundary conditions, heterogeneous material properties, and curvilinear coordinates. The importance of higher order accurate methods is not necessarily to make the numerical solution more accurate, but to reduce the computational cost for obtaining a solution within an acceptable error tolerance. This is because the accuracy in the solution can always be improved by reducing the grid size h. However, in practice, the available computational resources might not be large enough to solve the problem with a
Kudryavtsev, Oleg
2013-01-01
In the paper, we consider the problem of pricing options in wide classes of Lévy processes. We propose a general approach to the numerical methods based on a finite difference approximation for the generalized Black-Scholes equation. The goal of the paper is to incorporate the Wiener-Hopf factorization into finite difference methods for pricing options in Lévy models with jumps. The method is applicable for pricing barrier and American options. The pricing problem is reduced to the sequence of linear algebraic systems with a dense Toeplitz matrix; then the Wiener-Hopf factorization method is applied. We give an important probabilistic interpretation based on the infinitely divisible distributions theory to the Laurent operators in the correspondent factorization identity. Notice that our algorithm has the same complexity as the ones which use the explicit-implicit scheme, with a tridiagonal matrix. However, our method is more accurate. We support the advantage of the new method in terms of accuracy and convergence by using numerical experiments. PMID:24489518
2013-01-01
In the paper, we consider the problem of pricing options in wide classes of Lévy processes. We propose a general approach to the numerical methods based on a finite difference approximation for the generalized Black-Scholes equation. The goal of the paper is to incorporate the Wiener-Hopf factorization into finite difference methods for pricing options in Lévy models with jumps. The method is applicable for pricing barrier and American options. The pricing problem is reduced to the sequence of linear algebraic systems with a dense Toeplitz matrix; then the Wiener-Hopf factorization method is applied. We give an important probabilistic interpretation based on the infinitely divisible distributions theory to the Laurent operators in the correspondent factorization identity. Notice that our algorithm has the same complexity as the ones which use the explicit-implicit scheme, with a tridiagonal matrix. However, our method is more accurate. We support the advantage of the new method in terms of accuracy and convergence by using numerical experiments. PMID:24489518
Wang, Wei; Shu, Chi-Wang; Yee, H.C.; Sjögreen, Björn
2012-01-01
A new high order finite-difference method utilizing the idea of Harten ENO subcell resolution method is proposed for chemical reactive flows and combustion. In reaction problems, when the reaction time scale is very small, e.g., orders of magnitude smaller than the fluid dynamics time scales, the governing equations will become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present proposed method is a modified fractional step method which solves the convection step and reaction step separately. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with the computed flow variables in the shock region modified by the Harten subcell resolution idea. For numerical experiments, a fifth-order finite-difference WENO scheme and its anti-diffusion WENO variant are considered. A wide range of 1D and 2D scalar and Euler system test cases are investigated. Studies indicate that for the considered test cases, the new method maintains high order accuracy in space for smooth flows, and for stiff source terms with discontinuities, it can capture the correct propagation speed of discontinuities in very coarse meshes with reasonable CFL numbers.
A 3-dimensional finite-difference method for calculating the dynamic coefficients of seals
NASA Technical Reports Server (NTRS)
Dietzen, F. J.; Nordmann, R.
1989-01-01
A method to calculate the dynamic coefficients of seals with arbitrary geometry is presented. The Navier-Stokes equations are used in conjunction with the k-e turbulence model to describe the turbulent flow. These equations are solved by a full 3-dimensional finite-difference procedure instead of the normally used perturbation analysis. The time dependence of the equations is introduced by working with a coordinate system rotating with the precession frequency of the shaft. The results of this theory are compared with coefficients calculated by a perturbation analysis and with experimental results.
NASA Technical Reports Server (NTRS)
Liu, C.; Liu, Z.
1993-01-01
The fourth-order finite-difference scheme with fully implicit time-marching presently used to computationally study the spatial instability of planar Poiseuille flow incorporates a novel treatment for outflow boundary conditions that renders the buffer area as short as one wavelength. A semicoarsening multigrid method accelerates convergence for the implicit scheme at each time step; a line-distributive relaxation is developed as a robust fast solver that is efficient for anisotropic grids. Computational cost is no greater than that of explicit schemes, and excellent agreement with linear theory is obtained.
Experiments with explicit filtering for LES using a finite-difference method
NASA Technical Reports Server (NTRS)
Lund, T. S.; Kaltenbach, H. J.
1995-01-01
The equations for large-eddy simulation (LES) are derived formally by applying a spatial filter to the Navier-Stokes equations. The filter width as well as the details of the filter shape are free parameters in LES, and these can be used both to control the effective resolution of the simulation and to establish the relative importance of different portions of the resolved spectrum. An analogous, but less well justified, approach to filtering is more or less universally used in conjunction with LES using finite-difference methods. In this approach, the finite support provided by the computational mesh as well as the wavenumber-dependent truncation errors associated with the finite-difference operators are assumed to define the filter operation. This approach has the advantage that it is also 'automatic' in the sense that no explicit filtering: operations need to be performed. While it is certainly convenient to avoid the explicit filtering operation, there are some practical considerations associated with finite-difference methods that favor the use of an explicit filter. Foremost among these considerations is the issue of truncation error. All finite-difference approximations have an associated truncation error that increases with increasing wavenumber. These errors can be quite severe for the smallest resolved scales, and these errors will interfere with the dynamics of the small eddies if no corrective action is taken. Years of experience at CTR with a second-order finite-difference scheme for high Reynolds number LES has repeatedly indicated that truncation errors must be minimized in order to obtain acceptable simulation results. While the potential advantages of explicit filtering are rather clear, there is a significant cost associated with its implementation. In particular, explicit filtering reduces the effective resolution of the simulation compared with that afforded by the mesh. The resolution requirements for LES are usually set by the need to capture
A moving mesh finite difference method for equilibrium radiation diffusion equations
Yang, Xiaobo; Huang, Weizhang; Qiu, Jianxian
2015-10-01
An efficient moving mesh finite difference method is developed for the numerical solution of equilibrium radiation diffusion equations in two dimensions. The method is based on the moving mesh partial differential equation approach and moves the mesh continuously in time using a system of meshing partial differential equations. The mesh adaptation is controlled through a Hessian-based monitor function and the so-called equidistribution and alignment principles. Several challenging issues in the numerical solution are addressed. Particularly, the radiation diffusion coefficient depends on the energy density highly nonlinearly. This nonlinearity is treated using a predictor–corrector and lagged diffusion strategy. Moreover, the nonnegativity of the energy density is maintained using a cutoff method which has been known in literature to retain the accuracy and convergence order of finite difference approximation for parabolic equations. Numerical examples with multi-material, multiple spot concentration situations are presented. Numerical results show that the method works well for radiation diffusion equations and can produce numerical solutions of good accuracy. It is also shown that a two-level mesh movement strategy can significantly improve the efficiency of the computation.
A moving mesh finite difference method for equilibrium radiation diffusion equations
NASA Astrophysics Data System (ADS)
Yang, Xiaobo; Huang, Weizhang; Qiu, Jianxian
2015-10-01
An efficient moving mesh finite difference method is developed for the numerical solution of equilibrium radiation diffusion equations in two dimensions. The method is based on the moving mesh partial differential equation approach and moves the mesh continuously in time using a system of meshing partial differential equations. The mesh adaptation is controlled through a Hessian-based monitor function and the so-called equidistribution and alignment principles. Several challenging issues in the numerical solution are addressed. Particularly, the radiation diffusion coefficient depends on the energy density highly nonlinearly. This nonlinearity is treated using a predictor-corrector and lagged diffusion strategy. Moreover, the nonnegativity of the energy density is maintained using a cutoff method which has been known in literature to retain the accuracy and convergence order of finite difference approximation for parabolic equations. Numerical examples with multi-material, multiple spot concentration situations are presented. Numerical results show that the method works well for radiation diffusion equations and can produce numerical solutions of good accuracy. It is also shown that a two-level mesh movement strategy can significantly improve the efficiency of the computation.
NASA Astrophysics Data System (ADS)
Jia, Jinhong; Wang, Hong
2015-07-01
Numerical methods for space-fractional diffusion equations often generate dense or even full stiffness matrices. Traditionally, these methods were solved via Gaussian type direct solvers, which requires O (N3) of computational work per time step and O (N2) of memory to store where N is the number of spatial grid points in the discretization. In this paper we develop a preconditioned fast Krylov subspace iterative method for the efficient and faithful solution of finite difference methods (both steady-state and time-dependent) space-fractional diffusion equations with fractional derivative boundary conditions in one space dimension. The method requires O (N) of memory and O (Nlog N) of operations per iteration. Due to the application of effective preconditioners, significantly reduced numbers of iterations were achieved that further reduces the computational cost of the fast method. Numerical results are presented to show the utility of the method.
Samak, M. Mosleh E. Abu; Bakar, A. Ashrif A.; Kashif, Muhammad; Zan, Mohd Saiful Dzulkifly
2016-01-01
This paper discusses numerical analysis methods for different geometrical features that have limited interval values for typically used sensor wavelengths. Compared with existing Finite Difference Time Domain (FDTD) methods, the alternating direction implicit (ADI)-FDTD method reduces the number of sub-steps by a factor of two to three, which represents a 33% time savings in each single run. The local one-dimensional (LOD)-FDTD method has similar numerical equation properties, which should be calculated as in the previous method. Generally, a small number of arithmetic processes, which result in a shorter simulation time, are desired. The alternating direction implicit technique can be considered a significant step forward for improving the efficiency of unconditionally stable FDTD schemes. This comparative study shows that the local one-dimensional method had minimum relative error ranges of less than 40% for analytical frequencies above 42.85 GHz, and the same accuracy was generated by both methods.
New method of 3-D object recognition
NASA Astrophysics Data System (ADS)
He, An-Zhi; Li, Qun Z.; Miao, Peng C.
1991-12-01
In this paper, a new method of 3-D object recognition using optical techniques and a computer is presented. We perform 3-D object recognition using moire contour to obtain the object's 3- D coordinates, projecting drawings of the object in three coordinate planes to describe it and using a method of inquiring library of judgement to match objects. The recognition of a simple geometrical entity is simulated by computer and studied experimentally. The recognition of an object which is composed of a few simple geometrical entities is discussed.
NASA Astrophysics Data System (ADS)
Ying, Jinyong; Xie, Dexuan
2015-10-01
The Poisson-Boltzmann equation (PBE) is one widely-used implicit solvent continuum model for calculating electrostatics of ionic solvated biomolecule. In this paper, a new finite element and finite difference hybrid method is presented to solve PBE efficiently based on a special seven-overlapped box partition with one central box containing the solute region and surrounded by six neighboring boxes. In particular, an efficient finite element solver is applied to the central box while a fast preconditioned conjugate gradient method using a multigrid V-cycle preconditioning is constructed for solving a system of finite difference equations defined on a uniform mesh of each neighboring box. Moreover, the PBE domain, the box partition, and an interface fitted tetrahedral mesh of the central box can be generated adaptively for a given PQR file of a biomolecule. This new hybrid PBE solver is programmed in C, Fortran, and Python as a software tool for predicting electrostatics of a biomolecule in a symmetric 1:1 ionic solvent. Numerical results on two test models with analytical solutions and 12 proteins validate this new software tool, and demonstrate its high performance in terms of CPU time and memory usage.
Computationally efficient finite-difference modal method for the solution of Maxwell's equations.
Semenikhin, Igor; Zanuccoli, Mauro
2013-12-01
In this work, a new implementation of the finite-difference (FD) modal method (FDMM) based on an iterative approach to calculate the eigenvalues and corresponding eigenfunctions of the Helmholtz equation is presented. Two relevant enhancements that significantly increase the speed and accuracy of the method are introduced. First of all, the solution of the complete eigenvalue problem is avoided in favor of finding only the meaningful part of eigenmodes by using iterative methods. Second, a multigrid algorithm and Richardson extrapolation are implemented. Simultaneous use of these techniques leads to an enhancement in terms of accuracy, which allows a simple method such as the FDMM with a typical three-point difference scheme to be significantly competitive with an analytical modal method. PMID:24323014
The arbitrary order mixed mimetic finite difference method for the diffusion equation
Gyrya, Vitaliy; Lipnikov, Konstantin; Manzini, Gianmarco
2016-05-01
Here, we propose an arbitrary-order accurate mimetic finite difference (MFD) method for the approximation of diffusion problems in mixed form on unstructured polygonal and polyhedral meshes. As usual in the mimetic numerical technology, the method satisfies local consistency and stability conditions, which determines the accuracy and the well-posedness of the resulting approximation. The method also requires the definition of a high-order discrete divergence operator that is the discrete analog of the divergence operator and is acting on the degrees of freedom. The new family of mimetic methods is proved theoretically to be convergent and optimal error estimates for flux andmore » scalar variable are derived from the convergence analysis. A numerical experiment confirms the high-order accuracy of the method in solving diffusion problems with variable diffusion tensor. It is worth mentioning that the approximation of the scalar variable presents a superconvergence effect.« less
A Finite Difference Method for Modeling Migration of Impurities in Multilayer Systems
NASA Astrophysics Data System (ADS)
Tosa, V.; Kovacs, Katalin; Mercea, P.; Piringer, O.
2008-09-01
A finite difference method to solve the one-dimensional diffusion of impurities in a multilayer system was developed for the special case in which a partition coefficient K impose a ratio of the concentrations at the interface between two adiacent layers. The fictitious point method was applied to derive the algebraic equations for the mesh points at the interface, while for the non-uniform mesh points within the layers a combined method was used. The method was tested and then applied to calculate migration of impurities from multilayer systems into liquids or solids samples, in migration experiments performed for quality testing purposes. An application was developed in the field of impurities migrations from multilayer plastic packagings into food, a problem of increasing importance in food industry.
NASA Technical Reports Server (NTRS)
Strong, Stuart L.; Meade, Andrew J., Jr.
1992-01-01
Preliminary results are presented of a finite element/finite difference method (semidiscrete Galerkin method) used to calculate compressible boundary layer flow about airfoils, in which the group finite element scheme is applied to the Dorodnitsyn formulation of the boundary layer equations. The semidiscrete Galerkin (SDG) method promises to be fast, accurate and computationally efficient. The SDG method can also be applied to any smoothly connected airfoil shape without modification and possesses the potential capability of calculating boundary layer solutions beyond flow separation. Results are presented for low speed laminar flow past a circular cylinder and past a NACA 0012 airfoil at zero angle of attack at a Mach number of 0.5. Also shown are results for compressible flow past a flat plate for a Mach number range of 0 to 10 and results for incompressible turbulent flow past a flat plate. All numerical solutions assume an attached boundary layer.
A mimetic finite difference method for the Stokes problem with elected edge bubbles
Lipnikov, K; Berirao, L
2009-01-01
A new mimetic finite difference method for the Stokes problem is proposed and analyzed. The unstable P{sub 1}-P{sub 0} discretization is stabilized by adding a small number of bubble functions to selected mesh edges. A simple strategy for selecting such edges is proposed and verified with numerical experiments. The discretizations schemes for Stokes and Navier-Stokes equations must satisfy the celebrated inf-sup (or the LBB) stability condition. The stability condition implies a balance between discrete spaces for velocity and pressure. In finite elements, this balance is frequently achieved by adding bubble functions to the velocity space. The goal of this article is to show that the stabilizing edge bubble functions can be added only to a small set of mesh edges. This results in a smaller algebraic system and potentially in a faster calculations. We employ the mimetic finite difference (MFD) discretization technique that works for general polyhedral meshes and can accomodate non-uniform distribution of stabilizing bubbles.
Modelling the core convection using finite element and finite difference methods
NASA Astrophysics Data System (ADS)
Chan, K. H.; Li, Ligang; Liao, Xinhao
2006-08-01
Applications of both parallel finite element and finite difference methods to thermal convection in a rotating spherical shell modelling the fluid dynamics of the Earth's outer core are presented. The numerical schemes are verified by reproducing the convection benchmark test by Christensen et al. [Christensen, U.R., Aubert, J., Cardin, P., Dormy, E., Gibbons, S., Glatzmaier, G.A., Grote, E., Honkura, Y., Jones, C., Kono, M., Matsushima, M., Sakuraba, A., Takahashi, F., Tilgner, A., Wilcht, J., Zhang, K., 2001. A numerical dynamo benchmark. Phys. Earth Planet. Interiors 128, 25-34.]. Both global average and local characteristics agree satisfactorily with the benchmark solution. With the element-by-element (EBE) parallelization technique, the finite element code demonstrates nearly optimal linear scalability in computational speed. The finite difference code is also efficient and scalable by utilizing a parallel library Aztec [Tuminaro, R.S., Heroux, M., Hutchinson, S.A., Shadid, J.N., 1999. Official AZTEC User's Guide: Version 2.1.].
Methods for comparing 3D surface attributes
NASA Astrophysics Data System (ADS)
Pang, Alex; Freeman, Adam
1996-03-01
A common task in data analysis is to compare two or more sets of data, statistics, presentations, etc. A predominant method in use is side-by-side visual comparison of images. While straightforward, it burdens the user with the task of discerning the differences between the two images. The user if further taxed when the images are of 3D scenes. This paper presents several methods for analyzing the extent, magnitude, and manner in which surfaces in 3D differ in their attributes. The surface geometry are assumed to be identical and only the surface attributes (color, texture, etc.) are variable. As a case in point, we examine the differences obtained when a 3D scene is rendered progressively using radiosity with different form factor calculation methods. The comparison methods include extensions of simple methods such as mapping difference information to color or transparency, and more recent methods including the use of surface texture, perturbation, and adaptive placements of error glyphs.
Overview of finite difference Hartree-Fock method algorithm, implementation and application
NASA Astrophysics Data System (ADS)
Kobus, J.
2012-12-01
Two-dimensional, finite difference Hartree-Fock method has been in constant usage and development over the last two decades. The method has proved stable and efficient enough to be applied to dozens of diatomic molecules, even to systems as large as the thorium fluoride. Its latest version is presented and the dependence of its accuracy on the grid size and efficiency on the overrelaxation parameters are discussed. The method has been mainly used to develop and calibrate sequences of universal even-tempered and polarization-consistent basis sets and assess basis set truncation and superposition errors. Its modified version has proved useful in testing various exchange-correlation potentials within the density functional theory. The method has turned out to be a valuable source of reference values of total energies, multipole moments, static polarizabilities and hyperpolarizabilities (αzz, βzzz, γzzzz, Az,zz and Bzz,zz) for atoms, diatomic molecules and their ions. Recently, it has been modified to allow to calculate the electrical properties of homonuclear molecules and the results for the Li2, N2, F2 and O2 systems are presented. Electrical properties of the AlF, CS, KCl diatomics and of highly ionized krypton atom (Kr+32) are reported as well. Accuracy of both the matrix Hartree-Fock employing universal even-tempered basis sets and the finite difference Hartree-Fock methods is discussed and the basis set superposition errors of the dipole polarizability and the first hyperpolarizability of the FH molecule is reexamined. Basis set superposition errors are also discussed in case of the dipole polarizability and the second hyperpolarizability of the F2 system.
NASA Technical Reports Server (NTRS)
Schroeter, Jens; Wunsch, Carl
1986-01-01
The paper studies with finite difference nonlinear circulation models the uncertainties in interesting flow properties, such as western boundary current transport, potential and kinetic energy, owing to the uncertainty in the driving surface boundary condition. The procedure is based upon nonlinear optimization methods. The same calculations permit quantitative study of the importance of new information as a function of type, region of measurement and accuracy, providing a method to study various observing strategies. Uncertainty in a model parameter, the bottom friction coefficient, is studied in conjunction with uncertain measurements. The model is free to adjust the bottom friction coefficient such that an objective function is minimized while fitting a set of data to within prescribed bounds. The relative importance of the accuracy of the knowledge about the friction coefficient with respect to various kinds of observations is then quantified, and the possible range of the friction coefficients is calculated.
NASA Technical Reports Server (NTRS)
Shu, Chi-Wang
1998-01-01
This project is about the development of high order, non-oscillatory type schemes for computational fluid dynamics. Algorithm analysis, implementation, and applications are performed. Collaborations with NASA scientists have been carried out to ensure that the research is relevant to NASA objectives. The combination of ENO finite difference method with spectral method in two space dimension is considered, jointly with Cai [3]. The resulting scheme behaves nicely for the two dimensional test problems with or without shocks. Jointly with Cai and Gottlieb, we have also considered one-sided filters for spectral approximations to discontinuous functions [2]. We proved theoretically the existence of filters to recover spectral accuracy up to the discontinuity. We also constructed such filters for practical calculations.
Finite difference method to find period-one gait cycles of simple passive walkers
NASA Astrophysics Data System (ADS)
Dardel, Morteza; Safartoobi, Masoumeh; Pashaei, Mohammad Hadi; Ghasemi, Mohammad Hassan; Navaei, Mostafa Kazemi
2015-01-01
Passive dynamic walking refers to a class of bipedal robots that can walk down an incline with no actuation or control input. These bipeds are sensitive to initial conditions due to their style of walking. According to small basin of attraction of passive limit cycles, it is important to start with an initial condition in the basin of attraction of stable walking (limit cycle). This paper presents a study of the simplest passive walker with point and curved feet. A new approach is proposed to find proper initial conditions for a pair of stable and unstable period-one gait limit cycles. This methodology is based on finite difference method which can solve the nonlinear differential equations of motion on a discrete time. Also, to investigate the physical configurations of the walkers and the environmental influence such as the slope angle, the parameter analysis is applied. Numerical simulations reveal the performance of the presented method in finding two stable and unstable gait patterns.
NASA Astrophysics Data System (ADS)
Arias-Ramirez, Walter; Olson, Britton; Wolf, William; Lawrence Livermore National Laboratory Team; University of Campinas Team
2015-11-01
The suitability of a continuing forcing immersed boundary method (IBM) combined with a high-order finite difference method is examined on several acoustic scattering problems. A suite of two-dimensional numerical simulations of canonical cases are conducted with the aim of analyzing the error behavior associated with the IBM, through wave reflection, wave diffraction, and the shock-boundary layer interaction phenomena. The compressible Navier-Stokes equations are solved using the Miranda code developed at Lawrence Livermore National Laboratory. Comparison of analytical solution against numerical results is shown for different flow parameters. Preliminary results indicate that the continuing forcing approach has the largest error in wave reflection compared to analytical solution. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344.
NASA Astrophysics Data System (ADS)
Guda, A. A.; Guda, S. A.; Soldatov, M. A.; Lomachenko, K. A.; Bugaev, A. L.; Lamberti, C.; Gawelda, W.; Bressler, C.; Smolentsev, G.; Soldatov, A. V.; Joly, Y.
2016-05-01
Finite difference method (FDM) implemented in the FDMNES software [Phys. Rev. B, 2001, 63, 125120] was revised. Thorough analysis shows, that the calculated diagonal in the FDM matrix consists of about 96% zero elements. Thus a sparse solver would be more suitable for the problem instead of traditional Gaussian elimination for the diagonal neighbourhood. We have tried several iterative sparse solvers and the direct one MUMPS solver with METIS ordering turned out to be the best. Compared to the Gaussian solver present method is up to 40 times faster and allows XANES simulations for complex systems already on personal computers. We show applicability of the software for metal-organic [Fe(bpy)3]2+ complex both for low spin and high spin states populated after laser excitation.
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.
Numerical analysis of polarization gratings using the finite-difference time-domain method
Oh, Chulwoo; Escuti, Michael J.
2007-10-15
We report the first full numerical analysis of polarization gratings (PGs), and study their most general properties and limits by using the finite-difference time-domain (FDTD) method. In this way, we avoid limiting assumptions on material properties or grating dimensions (e.g., no paraxial approximations) and provide a more complete understanding of PG diffraction behavior. We identify the fundamental delineation between diffraction regimes (thin versus thick) for anisotropic gratings and determine the conditions for {approx_equal}100% diffraction efficiency in the framework of the coupled-wave {rho} and Q parameters. Diffraction characteristics including the efficiency, spectral response, and polarization sensitivity are investigated for the two primary types of PGs with linear and circular birefringence. The angular response and finite-grating behavior (i.e., pixelation) are also examined. Comparisons with previous analytic approximations, where applicable, show good agreement.
NASA Technical Reports Server (NTRS)
Kishoni, Doron; Taasan, Shlomo
1994-01-01
Solution of the wave equation using techniques such as finite difference or finite element methods can model elastic wave propagation in solids. This requires mapping the physical geometry into a computational domain whose size is governed by the size of the physical domain of interest and by the required resolution. This computational domain, in turn, dictates the computer memory requirements as well as the calculation time. Quite often, the physical region of interest is only a part of the whole physical body, and does not necessarily include all the physical boundaries. Reduction of the calculation domain requires positioning an artificial boundary or region where a physical boundary does not exist. It is important however that such a boundary, or region, will not affect the internal domain, i.e., it should not cause reflections that propagate back into the material. This paper concentrates on the issue of constructing such a boundary region.
Low-dispersion finite difference methods for acoustic waves in a pipe
NASA Technical Reports Server (NTRS)
Davis, Sanford
1991-01-01
A new algorithm for computing one-dimensional acoustic waves in a pipe is demonstrated by solving the acoustic equations as an initial-boundary-value problem. Conventional dissipation-free second-order finite difference methods suffer severe phase distortion for grids with less that about ten mesh points per wavelength. Using the signal generation by a piston in a duct as an example, transient acoustic computations are presented using a new compact three-point algorithm which allows about 60 percent fewer mesh points per wavelength. Both pulse and harmonic excitation are considered. Coupling of the acoustic signal with the pipe resonant modes is shown to generate a complex transient wave with rich harmonic content.
Free transverse vibration of a wrinkled annular thin film by using finite difference method
NASA Astrophysics Data System (ADS)
Wang, C. G.; Liu, Y. P.; Lan, L.; Tan, H. F.
2016-02-01
This paper investigates the free transverse vibration of a wrinkled annular thin film. The non-dimensional Hamilton motion equation of the wrinkled annular thin film is established, which is solved by using the finite difference method to acquire the vibration frequency and mode. The predicted vibration characteristics are verified by the experimental measurements based on the digital image correlation (DIC) technique. The results show that wrinkles have great effects on the vibration of the annular thin film. Especially for the heavily wrinkled cases, the local-global interactive mode dominates the vibration of the annular thin film. The frequency increases as the wrinkling level increases which is mainly due to the increased nonlinear geometric stiffness. The results provide favorable supports for understanding the role of nonlinear wrinkling on the vibration of thin films.
NASA Astrophysics Data System (ADS)
Putri, Selmi; Arif, Idam; Khotimah, Siti Nurul
2015-04-01
In this study, peritoneal dialysis transport system was numerically simulated using finite difference method. The increase in the intraperitoneal pressure due to coughing has a high value outside the working area of the void volume fraction of the hydrostatic pressure θ(P). Therefore to illustrate the effects of the pressure increment, the pressure of working area is chosen between 1 and 3 mmHg. The effects of increased pressure in peritoneal tissue cause more fluid to flow into the blood vessels and lymph. Furthermore, the increased pressure in peritoneal tissue makes the volumetric flux jv and solute flux js across the tissue also increase. The more fluid flow into the blood vessels and lymph causes the fluid to flow into tissue qv and the glucose flow qs to have more negative value and also decreases the glucose concentration CG in the tissue.
CUDA Fortran acceleration for the finite-difference time-domain method
NASA Astrophysics Data System (ADS)
Hadi, Mohammed F.; Esmaeili, Seyed A.
2013-05-01
A detailed description of programming the three-dimensional finite-difference time-domain (FDTD) method to run on graphical processing units (GPUs) using CUDA Fortran is presented. Two FDTD-to-CUDA thread-block mapping designs are investigated and their performances compared. Comparative assessment of trade-offs between GPU's shared memory and L1 cache is also discussed. This presentation is for the benefit of FDTD programmers who work exclusively with Fortran and are reluctant to port their codes to C in order to utilize GPU computing. The derived CUDA Fortran code is compared with an optimized CPU version that runs on a workstation-class CPU to present a realistic GPU to CPU run time comparison and thus help in making better informed investment decisions on FDTD code redesigns and equipment upgrades. All analyses are mirrored with CUDA C simulations to put in perspective the present state of CUDA Fortran development.
NASA Technical Reports Server (NTRS)
Mccoy, M. J.
1980-01-01
Various finite difference techniques used to solve Laplace's equation are compared. Curvilinear coordinate systems are used on two dimensional regions with irregular boundaries, specifically, regions around circles and airfoils. Truncation errors are analyzed for three different finite difference methods. The false boundary method and two point and three point extrapolation schemes, used when having the Neumann boundary condition are considered and the effects of spacing and nonorthogonality in the coordinate systems are studied.
NASA Astrophysics Data System (ADS)
Gao, Junhui
2013-05-01
Overlap grid is usually used in numerical simulation of flow with complex geometry by high order finite difference scheme. It is difficult to generate overlap grid and the connectivity information between adjacent blocks, especially when interpolation is required for non-coincident overlap grids. In this study, an interface flux reconstruction (IFR) method is proposed for numerical simulation using high order finite difference scheme with multi-block structured grids. In this method the neighboring blocks share a common face, and the fluxes on each block are matched to set the boundary conditions for each interior block. Therefore this method has the promise of allowing discontinuous grids on either side of an interior block interface. The proposed method is proven to be stable for 7-point central DRP scheme coupled with 4-point and 5-point boundary closure schemes, as well as the 4th order compact scheme coupled with 3rd order boundary closure scheme. Four problems are numerically solved with the developed code to validate the interface flux reconstruction method in this study. The IFR method coupled with the 4th order DRP scheme or compact scheme is validated to be 4th order accuracy with one and two dimensional waves propagation problems. Two dimensional pulse propagation in mean flow is computed with wavy mesh to demonstrate the ability of the proposed method for non-uniform grid. To demonstrate the ability of the proposed method for complex geometry, sound scattering by two cylinders is simulated and the numerical results are compared with the analytical data. It is shown that the numerical results agree well with the analytical data. Finally the IFR method is applied to simulate viscous flow pass a cylinder at Reynolds number 150 to show its capability for viscous problem. The computed pressure coefficient on the cylinder surface, the frequency of vortex shedding, the lift and drag coefficients are presented. The numerical results are compared with the data
Mimetic finite difference method for the stokes problem on polygonal meshes
Lipnikov, K; Beirao Da Veiga, L; Gyrya, V; Manzini, G
2009-01-01
Various approaches to extend the finite element methods to non-traditional elements (pyramids, polyhedra, etc.) have been developed over the last decade. Building of basis functions for such elements is a challenging task and may require extensive geometry analysis. The mimetic finite difference (MFD) method has many similarities with low-order finite element methods. Both methods try to preserve fundamental properties of physical and mathematical models. The essential difference is that the MFD method uses only the surface representation of discrete unknowns to build stiffness and mass matrices. Since no extension inside the mesh element is required, practical implementation of the MFD method is simple for polygonal meshes that may include degenerate and non-convex elements. In this article, we develop a MFD method for the Stokes problem on arbitrary polygonal meshes. The method is constructed for tensor coefficients, which will allow to apply it to the linear elasticity problem. The numerical experiments show the second-order convergence for the velocity variable and the first-order for the pressure.
Recognition methods for 3D textured surfaces
NASA Astrophysics Data System (ADS)
Cula, Oana G.; Dana, Kristin J.
2001-06-01
Texture as a surface representation is the subject of a wide body of computer vision and computer graphics literature. While texture is always associated with a form of repetition in the image, the repeating quantity may vary. The texture may be a color or albedo variation as in a checkerboard, a paisley print or zebra stripes. Very often in real-world scenes, texture is instead due to a surface height variation, e.g. pebbles, gravel, foliage and any rough surface. Such surfaces are referred to here as 3D textured surfaces. Standard texture recognition algorithms are not appropriate for 3D textured surfaces because the appearance of these surfaces changes in a complex manner with viewing direction and illumination direction. Recent methods have been developed for recognition of 3D textured surfaces using a database of surfaces observed under varied imaging parameters. One of these methods is based on 3D textons obtained using K-means clustering of multiscale feature vectors. Another method uses eigen-analysis originally developed for appearance-based object recognition. In this work we develop a hybrid approach that employs both feature grouping and dimensionality reduction. The method is tested using the Columbia-Utrecht texture database and provides excellent recognition rates. The method is compared with existing recognition methods for 3D textured surfaces. A direct comparison is facilitated by empirical recognition rates from the same texture data set. The current method has key advantages over existing methods including requiring less prior information on both the training and novel images.
On a finite-difference method for solving transient viscous flow problems
NASA Technical Reports Server (NTRS)
Li, C. P.
1983-01-01
A method has been developed to solve the unsteady, compressible Navier-Stokes equation with the property of consistency and the ability of minimizing the equation stiffness. It relies on innovative extensions of the state-of-the-art finite-difference techniques and is composed of: (1) the upwind scheme for split-flux and the central scheme for conventional flux terms in the inviscid and viscous regions, respectively; (2) the characteristic treatment of both inviscid and viscous boundaries; (3) an ADI procedure compatible with interior and boundary points; and (4) a scalar matrix coefficient including viscous terms. The performance of this method is assessed with four sample problems; namely, a standing shock in the Laval duct, a shock reflected from the wall, the shock-induced boundary-layer separation, and a transient internal nozzle flow. The results from the present method, an existing hybrid block method, and a well-known two-step explicit method are compared and discussed. It is concluded that this method has an optimal trade-off between the solution accuracy and computational economy, and other desirable properties for analyzing transient viscous flow problems.
NASA Technical Reports Server (NTRS)
Nordstrom, Jan; Carpenter, Mark H.
1998-01-01
Boundary and interface conditions for high order finite difference methods applied to the constant coefficient Euler and Navier-Stokes equations are derived. The boundary conditions lead to strict and strong stability. The interface conditions are stable and conservative even if the finite difference operators and mesh sizes vary from domain to domain. Numerical experiments show that the new conditions also lead to good results for the corresponding nonlinear problems.
New Finite Difference Methods Based on IIM for Inextensible Interfaces in Incompressible Flows
Li, Zhilin; Lai, Ming-Chih
2012-01-01
In this paper, new finite difference methods based on the augmented immersed interface method (IIM) are proposed for simulating an inextensible moving interface in an incompressible two-dimensional flow. The mathematical models arise from studying the deformation of red blood cells in mathematical biology. The governing equations are incompressible Stokes or Navier-Stokes equations with an unknown surface tension, which should be determined in such a way that the surface divergence of the velocity is zero along the interface. Thus, the area enclosed by the interface and the total length of the interface should be conserved during the evolution process. Because of the nonlinear and coupling nature of the problem, direct discretization by applying the immersed boundary or immersed interface method yields complex nonlinear systems to be solved. In our new methods, we treat the unknown surface tension as an augmented variable so that the augmented IIM can be applied. Since finding the unknown surface tension is essentially an inverse problem that is sensitive to perturbations, our regularization strategy is to introduce a controlled tangential force along the interface, which leads to a least squares problem. For Stokes equations, the forward solver at one time level involves solving three Poisson equations with an interface. For Navier-Stokes equations, we propose a modified projection method that can enforce the pressure jump condition corresponding directly to the unknown surface tension. Several numerical experiments show good agreement with other results in the literature and reveal some interesting phenomena. PMID:23795308
Sun, Qiang; Wu, Guo Xiong
2013-03-01
A mathematical model and a numerical solution procedure are developed to simulate flow field through a 3D permeable vessel with multibranches embedded in a solid tumour. The model is based on Poisseuille's law for the description of the flow through the vessels, Darcy's law for the fluid field inside the tumour interstitium, and Starling's law for the flux transmitted across the vascular walls. The solution procedure is based on a coupled method, in which the finite difference method is used for the flow in the vessels and the boundary element method is used for the flow in the tumour. When vessels meet each other at a junction, the pressure continuity and mass conservation are imposed at the junction. Three typical representative structures within the tumour vasculature, symmetrical dichotomous branching, asymmetrical bifurcation with uneven radius of daughter vessels and trifurcation, are investigated in detail as case studies. These results have demonstrated the features of tumour flow environment by the pressure distributions and flow velocity field. PMID:23345121
Simulations of SH wave scattering due to cracks by the 2-D finite difference method
NASA Astrophysics Data System (ADS)
Suzuki, Y.; Kawahara, J.; Okamoto, T.; Miyashita, K.
2006-05-01
We simulate SH wave scattering by 2-D parallel cracks using the finite difference method (FDM), instead of the popularly used boundary integral equation method (BIEM). Here special emphasis is put on simplicity; we apply a standard FDM (fourth-order velocity-stress scheme with a staggered grid) to media in cluding traction-freecracks, which are expressed by arrays of grid points with zero traction. Two types of accuracy tests based oncomparison with a reliable BIEM, suggest that the present method gives practically sufficient accuracy, except for the wavefields in the vicinity of cracks, which can be well handled if the second-order FDM is used instead. As an application of this method, we also simulate wave propagation in media with randomly distributed cracks of the same length. We experimentally determine the attenuation and velocity dispersion induced by scattering from the synthetic seismograms, using a waveform averaging technique. It is shown that the results are well explained by a theory based on the Foldy approximation for crack densities of up to about 01. The presence of a free surface does not affect the validity of the theory. A preliminary experiment also suggests that the validity will not change even for multi-scale cracks.
An exploratory study of a finite difference method for calculating unsteady transonic potential flow
NASA Technical Reports Server (NTRS)
Bennett, R. M.; Bland, S. R.
1979-01-01
A method for calculating transonic flow over steady and oscillating airfoils was developed by Isogai. The full potential equation is solved with a semi-implicit, time-marching, finite difference technique. Steady flow solutions are obtained from time asymptotic solutions for a steady airfoil. Corresponding oscillatory solutions are obtained by initiating an oscillation and marching in time for several cycles until a converged periodic solution is achieved. The method is described in general terms and results for the case of an airfoil with an oscillating flap are presented for Mach numbers 0.500 and 0.875. Although satisfactory results are obtained for some reduced frequencies, it is found that the numerical technique generates spurious oscillations in the indicial response functions and in the variation of the aerodynamic coefficients with reduced frequency. These oscillations are examined with a dynamic data reduction method to evaluate their effects and trends with reduced frequency and Mach number. Further development of the numerical method is needed to eliminate these oscillations.
Semi-implicit finite difference methods for three-dimensional shallow water flow
Casulli, Vincenzo; Cheng, Ralph T.
1992-01-01
A semi-implicit finite difference method for the numerical solution of three-dimensional shallow water flows is presented and discussed. The governing equations are the primitive three-dimensional turbulent mean flow equations where the pressure distribution in the vertical has been assumed to be hydrostatic. In the method of solution a minimal degree of implicitness has been adopted in such a fashion that the resulting algorithm is stable and gives a maximal computational efficiency at a minimal computational cost. At each time step the numerical method requires the solution of one large linear system which can be formally decomposed into a set of small three-diagonal systems coupled with one five-diagonal system. All these linear systems are symmetric and positive definite. Thus the existence and uniquencess of the numerical solution are assured. When only one vertical layer is specified, this method reduces as a special case to a semi-implicit scheme for solving the corresponding two-dimensional shallow water equations. The resulting two- and three-dimensional algorithm has been shown to be fast, accurate and mass-conservative and can also be applied to simulate flooding and drying of tidal mud-flats in conjunction with three-dimensional flows. Furthermore, the resulting algorithm is fully vectorizable for an efficient implementation on modern vector computers.
NASA Astrophysics Data System (ADS)
Liao, Fei; Ye, Zhengyin
2015-12-01
Despite significant progress in recent computational techniques, the accurate numerical simulations, such as direct-numerical simulation and large-eddy simulation, are still challenging. For accurate calculations, the high-order finite difference method (FDM) is usually adopted with coordinate transformation from body-fitted grid to Cartesian grid. But this transformation might lead to failure in freestream preservation with the geometric conservation law (GCL) violated, particularly in high-order computations. GCL identities, including surface conservation law (SCL) and volume conservation law (VCL), are very important in discretization of high-order FDM. To satisfy GCL, various efforts have been made. An early and successful approach was developed by Thomas and Lombard [6] who used the conservative form of metrics to cancel out metric terms to further satisfy SCL. Visbal and Gaitonde [7] adopted this conservative form of metrics for SCL identities and satisfied VCL identity through invoking VCL equation to acquire the derivative of Jacobian in computation on moving and deforming grids with central compact schemes derived by Lele [5]. Later, using the metric technique from Visbal and Gaitonde [7], Nonomura et al. [8] investigated the freestream and vortex preservation properties of high-order WENO and WCNS on stationary curvilinear grids. A conservative metric method (CMM) was further developed by Deng et al. [9] with stationary grids, and detailed discussion about the innermost difference operator of CMM was shown with proof and corresponding numerical test cases. Noticing that metrics of CMM is asymmetrical without coordinate-invariant property, Deng et al. proposed a symmetrical CMM (SCMM) [12] by using the symmetric forms of metrics derived by Vinokur and Yee [10] to further eliminate asymmetric metric errors with stationary grids considered only. The research from Abe et al. [11] presented new asymmetric and symmetric conservative forms of time metrics and
Conservative high-order-accurate finite-difference methods for curvilinear grids
NASA Technical Reports Server (NTRS)
Rai, Man M.; Chakrvarthy, Sukumar
1993-01-01
Two fourth-order-accurate finite-difference methods for numerically solving hyperbolic systems of conservation equations on smooth curvilinear grids are presented. The first method uses the differential form of the conservation equations; the second method uses the integral form of the conservation equations. Modifications to these schemes, which are required near boundaries to maintain overall high-order accuracy, are discussed. An analysis that demonstrates the stability of the modified schemes is also provided. Modifications to one of the schemes to make it total variation diminishing (TVD) are also discussed. Results that demonstrate the high-order accuracy of both schemes are included in the paper. In particular, a Ringleb-flow computation demonstrates the high-order accuracy and the stability of the boundary and near-boundary procedures. A second computation of supersonic flow over a cylinder demonstrates the shock-capturing capability of the TVD methodology. An important contribution of this paper is the dear demonstration that higher order accuracy leads to increased computational efficiency.
Simulation of optical devices using parallel finite-difference time-domain method
NASA Astrophysics Data System (ADS)
Li, Kang; Kong, Fanmin; Mei, Liangmo; Liu, Xin
2005-11-01
This paper presents a new parallel finite-difference time-domain (FDTD) numerical method in a low-cost network environment to stimulate optical waveguide characteristics. The PC motherboard based cluster is used, as it is relatively low-cost, reliable and has high computing performance. Four clusters are networked by fast Ethernet technology. Due to the simplicity nature of FDTD algorithm, a native Ethernet packet communication mechanism is used to reduce the overhead of the communication between the adjacent clusters. To validate the method, a microcavity ring resonator based on semiconductor waveguides is chosen as an instance of FDTD parallel computation. Speed-up rate under different division density is calculated. From the result we can conclude that when the decomposing size reaches a certain point, a good parallel computing speed up will be maintained. This simulation shows that through the overlapping of computation and communication method and controlling the decomposing size, the overhead of the communication of the shared data will be conquered. The result indicates that the implementation can achieve significant speed up for the FDTD algorithm. This will enable us to tackle the larger real electromagnetic problem by the low-cost PC clusters.
Use of the finite-difference time-domain method in electromagnetic dosimetry
Sullivan, D.M.
1987-01-01
Although there are acceptable methods for calculating whole body electromagnetic absorption, no completely acceptable method for calculating the local specific absorption rate (SAR) at points within the body has been developed. Frequency domain methods, such as the method of moments (MoM) have achieved some success; however, the MoM requires computer storage on the order of (3N)/sup 2/, and computation time on the order of (3N)/sup 3/ where N is the number of cells. The finite-difference time-domain (FDTD) method has been employed extensively in calculating the scattering from metallic objects, and recently is seeing some use in calculating the interaction of EM fields with complex, lossy dielectric bodies. Since the FDTD method has storage and time requirements proportional to N, it presents an attractive alternative to calculating SAR distribution in large bodies. This dissertation describes the FDTD method and evaluates it by comparing its results with analytic solutions in 2 and 3 dimensions. The results obtained demonstrate that the FDTD method is capable of calculating internal SAR distribution with acceptable accuracy. The construction of a data base to provide detailed, inhomogeneous man models for use with the FDTD method is described. Using this construction method, a model of 40,000 1.31 cm. cells is developed for use at 350 MHz, and another model consisting of 5000 2.62 cm. cells is developed for use at 100 MHz. To add more realism to the problem, a ground plane is added to the FDTD software. The needed changes to the software are described, along with a test which confirms its accuracy. Using the CRAY II supercomputer, SAR distributions in human models are calculated using incident frequencies of 100 MHz and 350 MHz for three different cases: (1) A homogeneous man model in free space, (2) an inhomogeneous man model in free space, and (3) an inhomogeneous man model standing on a ground plane.
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.
Stochastic finite difference lattice Boltzmann method for steady incompressible viscous flows
Fu, S.C.; So, R.M.C.; Leung, W.W.F.
2010-08-20
With the advent of state-of-the-art computers and their rapid availability, the time is ripe for the development of efficient uncertainty quantification (UQ) methods to reduce the complexity of numerical models used to simulate complicated systems with incomplete knowledge and data. The spectral stochastic finite element method (SSFEM) which is one of the widely used UQ methods, regards uncertainty as generating a new dimension and the solution as dependent on this dimension. A convergent expansion along the new dimension is then sought in terms of the polynomial chaos system, and the coefficients in this representation are determined through a Galerkin approach. This approach provides an accurate representation even when only a small number of terms are used in the spectral expansion; consequently, saving in computational resource can be realized compared to the Monte Carlo (MC) scheme. Recent development of a finite difference lattice Boltzmann method (FDLBM) that provides a convenient algorithm for setting the boundary condition allows the flow of Newtonian and non-Newtonian fluids, with and without external body forces to be simulated with ease. Also, the inherent compressibility effect in the conventional lattice Boltzmann method, which might produce significant errors in some incompressible flow simulations, is eliminated. As such, the FDLBM together with an efficient UQ method can be used to treat incompressible flows with built in uncertainty, such as blood flow in stenosed arteries. The objective of this paper is to develop a stochastic numerical solver for steady incompressible viscous flows by combining the FDLBM with a SSFEM. Validation against MC solutions of channel/Couette, driven cavity, and sudden expansion flows are carried out.
Bringuier, Jonathan N.; Mittra, Raj
2012-01-01
A rigorous full-wave solution, via the Finite-Difference-Time-Domain (FDTD) method, is performed in an attempt to obtain realistic communication channel models for on-body wireless transmission in Body-Area-Networks (BANs), which are local data networks using the human body as a propagation medium. The problem of modeling the coupling between body mounted antennas is often not amenable to attack by hybrid techniques owing to the complex nature of the human body. For instance, the time-domain Green's function approach becomes more involved when the antennas are not conformal. Furthermore, the human body is irregular in shape and has dispersion properties that are unique. One consequence of this is that we must resort to modeling the antenna network mounted on the body in its entirety, and the number of degrees of freedom (DoFs) can be on the order of billions. Even so, this type of problem can still be modeled by employing a parallel version of the FDTD algorithm running on a cluster. Lastly, we note that the results of rigorous simulation of BANs can serve as benchmarks for comparison with the abundance of measurement data. PMID:23012575
NASA Astrophysics Data System (ADS)
Nikkar, Samira; Nordström, Jan
2015-06-01
A time-dependent coordinate transformation of a constant coefficient hyperbolic system of equations which results in a variable coefficient system of equations is considered. By applying the energy method, well-posed boundary conditions for the continuous problem are derived. Summation-by-Parts (SBP) operators for the space and time discretization, together with a weak imposition of boundary and initial conditions using Simultaneously Approximation Terms (SATs) lead to a provable fully-discrete energy-stable conservative finite difference scheme. We show how to construct a time-dependent SAT formulation that automatically imposes boundary conditions, when and where they are required. We also prove that a uniform flow field is preserved, i.e. the Numerical Geometric Conservation Law (NGCL) holds automatically by using SBP-SAT in time and space. The developed technique is illustrated by considering an application using the linearized Euler equations: the sound generated by moving boundaries. Numerical calculations corroborate the stability and accuracy of the new fully discrete approximations.
NASA Technical Reports Server (NTRS)
Ransom, Jonathan B.
2002-01-01
A multifunctional interface method with capabilities for variable-fidelity modeling and multiple method analysis is presented. The methodology provides an effective capability by which domains with diverse idealizations can be modeled independently to exploit the advantages of one approach over another. The multifunctional method is used to couple independently discretized subdomains, and it is used to couple the finite element and the finite difference methods. The method is based on a weighted residual variational method and is presented for two-dimensional scalar-field problems. A verification test problem and a benchmark application are presented, and the computational implications are discussed.
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.
An Adaptive Finite Difference Method for Hyperbolic Systems in OneSpace Dimension
Bolstad, John H.
1982-06-01
Many problems of physical interest have solutions which are generally quite smooth in a large portion of the region of interest, but have local phenomena such as shocks, discontinuities or large gradients which require much more accurate approximations or finer grids for reasonable accuracy. Examples are atmospheric fronts, ocean currents, and geological discontinuities. In this thesis we develop and partially analyze an adaptive finite difference mesh refinement algorithm for the initial boundary value problem for hyperbolic systems in one space dimension. The method uses clusters of uniform grids which can ''move'' along with pulses or steep gradients appearing in the calculation, and which are superimposed over a uniform coarse grid. Such refinements are created, destroyed, merged, separated, recursively nested or moved based on estimates of the local truncation error. We use a four-way linked tree and sequentially allocated deques (double-ended queues) to perform these operations efficiently. The local truncation error in the interior of the region is estimated using a three-step Richardson extrapolation procedure, which can also be considered a deferred correction method. At the boundaries we employ differences to estimate the error. Our algorithm was implemented using a portable, extensible Fortran preprocessor, to which we added records and pointers. The method is applied to three model problems: the first order wave equation, the second order wave equation, and the inviscid Burgers equation. For the first two model problems our algorithm is shown to be three to five times more efficient (in computing time) than the use of a uniform coarse mesh, for the same accuracy. Furthermore, to our knowledge, our algorithm is the only one which adaptively treats time-dependent boundary conditions for hyperbolic systems.
Finite difference and lead field methods in designing implantable ECG monitor.
Väisänen, Juho; Hyttinen, Jari; Malmivuo, Jaakko
2006-10-01
To minimize time-consuming and expensive in vitro and in vivo testing, information regarding the effects of implantation and the implants on measurements should be available during the designing of active implantable devices measuring bioelectric signals such as electrocardiograms (ECG). Modeling offers a fairly inexpensive and effective means of studying and demonstrating the effects of implantation on ECG measurements prior to any in vivo tests, and can thus provide the designer with valuable information. Finite difference model (FDM) and lead field approaches offer straightforward and effective modeling methods supporting the designing of active implantable ECG devices. The present study demonstrates such methods in developing and studying ECG implants. They were applied in demonstrating the effects of implant dimensions and of electrode implantation on the measurement sensitivity of the ECG device. The results of the simulations indicated that the interelectrode distance is the factor of the implant design determining the lead sensitivity. Other parameters related implant dimensions and shape have minor effect on the morphology of the ECG or on the average sensitivity of the measurement. This is shown for example when the interelectrode distance was reduced to 1/3 of original the average lead sensitivity decreased by 69.1% while larger relative changes in other dimensions produced clearly smaller changes. It was also observed here that implanting the electrodes deeper under the skin has major effects on the local sensitivities in heart muscle and thus affect to the morphology of the ECG. The study indicated also that non-conducting medium (i.e. implant insulated body) between the electrodes increases the sensitivity on heart muscle compared to cases where only electrodes are implanted. PMID:17031715
NASA Astrophysics Data System (ADS)
Hochgraf, Kelsey
Auralization methods have been used for a long time to simulate the acoustics of a concert hall for different seat positions. The goal of this thesis was to apply the concept of auralization to a larger audience area that the listener could walk through to compare differences in acoustics for a wide range of seat positions. For this purpose, the acoustics of Rensselaer's Experimental Media and Performing Arts Center (EMPAC) Concert Hall were simulated to create signals for a 136 channel wave field synthesis (WFS) system located at Rensselaer's Collaborative Research Augmented Immersive Virtual Environment (CRAIVE) Laboratory. By allowing multiple people to dynamically experience the concert hall's acoustics at the same time, this research gained perspective on what is important for achieving objective accuracy and subjective plausibility in an auralization. A finite difference time domain (FDTD) simulation on a three-dimensional face-centered cubic grid, combined at a crossover frequency of 800 Hz with a CATT-Acoustic(TM) simulation, was found to have a reverberation time, direct to reverberant sound energy ratio, and early reflection pattern that more closely matched measured data from the hall compared to a CATT-Acoustic(TM) simulation and other hybrid simulations. In the CRAIVE lab, nine experienced listeners found all hybrid auralizations (with varying source location, grid resolution, crossover frequency, and number of loudspeakers) to be more perceptually plausible than the CATT-Acoustic(TM) auralization. The FDTD simulation required two days to compute, while the CATT-Acoustic(TM) simulation required three separate TUCT(TM) computations, each taking four hours, to accommodate the large number of receivers. Given the perceptual advantages realized with WFS for auralization of a large, inhomogeneous sound field, it is recommended that hybrid simulations be used in the future to achieve more accurate and plausible auralizations. Predictions are made for a
NASA Technical Reports Server (NTRS)
Butler, T. D.; Weatherill, W. H.; Sebastian, J. D.; Ehlers, F. E.
1977-01-01
The design and usage of a pilot program using a finite difference method for calculating the pressure distributions over harmonically oscillating wings in transonic flow are discussed. The procedure used is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The steady velocity potential which must be obtained from some other program, is required for input. The unsteady differential equation is linear, complex in form with spatially varying coefficients. Because sinusoidal motion is assumed, time is not a variable. The numerical solution is obtained through a finite difference formulation and a line relaxation solution method.
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
NASA Astrophysics Data System (ADS)
Wang, Yue
A new variable grid-size and time-step finite-difference (FD) method is developed and applied to three different geophysical problems: simulation of tube waves in boreholes, three-dimensional (3-D) ground-motion simulation in sedimentary basin models, and reverse-time migration of multicomponent data. Unlike the conventional FD method, which uses a fixed grid-size and time-step for the entire model region, spatially variable grid-sizes and time-steps are used to achieve the optimal computational efficiency. For tube wave simulations, a fine grid-spacing is used for simulation inside the borehole region, while a coarse grid is used in the exterior region. While the stability condition requires a very fine time step for the fine grid, a variable time-step method provides coarse time steps for simulation in the coarse grid. Variable grid-size and time-step changes are used to achieve both accuracy and efficiency in the simulations. Numerical tests are performed for the Bayou Choctaw salt-flank model with different borehole models. The results show the important borehole effects on the seismic wavefield for a realistic source bandwidth. The combination of variable grid-size and time-step methods reduces computational costs by several orders of magnitude for the borehole models. Viscoelastic 3-D simulations are performed for a three-layer Salt Lake basin model. The near-surface unconsolidated layer is modeled with a fine grid, and the deep part of the model is modeled by a coarse grid. Simulation results show that the 3-D basin features and the shallow layer significantly affect the amplitude and duration time of the ground motion. In the elastic case, the approximation by 2-D modeling is insufficient to simulate the 3-D ground motion response. A basin model without a shallow low-velocity layer underestimates the ground motion duration and cumulative kinetic energy by 50% or more. The simulation of a Bingham Mine blast suggests that a lower S-velocity should be used to
NASA Technical Reports Server (NTRS)
Chen, G.; Zheng, Q.; Coleman, M.; Weerakoon, S.
1983-01-01
This paper briefly reviews convergent finite difference schemes for hyperbolic initial boundary value problems and their applications to boundary control systems of hyperbolic type which arise in the modelling of vibrations. These difference schemes are combined with the primal and the dual approaches to compute the optimal control in the unconstrained case, as well as the case when the control is subject to inequality constraints. Some of the preliminary numerical results are also presented.
Simulation model of stratified thermal energy storage tank using finite difference method
NASA Astrophysics Data System (ADS)
Waluyo, Joko
2016-06-01
Stratified TES tank is normally used in the cogeneration plant. The stratified TES tanks are simple, low cost, and equal or superior in thermal performance. The advantage of TES tank is that it enables shifting of energy usage from off-peak demand for on-peak demand requirement. To increase energy utilization in a stratified TES tank, it is required to build a simulation model which capable to simulate the charging phenomenon in the stratified TES tank precisely. This paper is aimed to develop a novel model in addressing the aforementioned problem. The model incorporated chiller into the charging of stratified TES tank system in a closed system. The model was developed in one-dimensional type involve with heat transfer aspect. The model covers the main factors affect to degradation of temperature distribution namely conduction through the tank wall, conduction between cool and warm water, mixing effect on the initial flow of the charging as well as heat loss to surrounding. The simulation model is developed based on finite difference method utilizing buffer concept theory and solved in explicit method. Validation of the simulation model is carried out using observed data obtained from operating stratified TES tank in cogeneration plant. The temperature distribution of the model capable of representing S-curve pattern as well as simulating decreased charging temperature after reaching full condition. The coefficient of determination values between the observed data and model obtained higher than 0.88. Meaning that the model has capability in simulating the charging phenomenon in the stratified TES tank. The model is not only capable of generating temperature distribution but also can be enhanced for representing transient condition during the charging of stratified TES tank. This successful model can be addressed for solving the limitation temperature occurs in charging of the stratified TES tank with the absorption chiller. Further, the stratified TES tank can be
NASA Astrophysics Data System (ADS)
Pedrozo, Héctor A.; Rosenberger, Mario R.; Schvezov, Carlos E.
2016-06-01
The solution by the Finite Difference Method of the Richards equation written as a function of the degree of saturation of the domain and the matrix potential is obtained and the convergence of the solutions is analyzed. The necessary time and spatial sizes for convergence are obtained and established.
Appelo, D; Petersson, N A
2007-12-17
The isotropic elastic wave equation governs the propagation of seismic waves caused by earthquakes and other seismic events. It also governs the propagation of waves in solid material structures and devices, such as gas pipes, wave guides, railroad rails and disc brakes. In the vast majority of wave propagation problems arising in seismology and solid mechanics there are free surfaces. These free surfaces have, in general, complicated shapes and are rarely flat. Another feature, characterizing problems arising in these areas, is the strong heterogeneity of the media, in which the problems are posed. For example, on the characteristic length scales of seismological problems, the geological structures of the earth can be considered piecewise constant, leading to models where the values of the elastic properties are also piecewise constant. Large spatial contrasts are also found in solid mechanics devices composed of different materials welded together. The presence of curved free surfaces, together with the typical strong material heterogeneity, makes the design of stable, efficient and accurate numerical methods for the elastic wave equation challenging. Today, many different classes of numerical methods are used for the simulation of elastic waves. Early on, most of the methods were based on finite difference approximations of space and time derivatives of the equations in second order differential form (displacement formulation), see for example [1, 2]. The main problem with these early discretizations were their inability to approximate free surface boundary conditions in a stable and fully explicit manner, see e.g. [10, 11, 18, 20]. The instabilities of these early methods were especially bad for problems with materials with high ratios between the P-wave (C{sub p}) and S-wave (C{sub s}) velocities. For rectangular domains, a stable and explicit discretization of the free surface boundary conditions is presented in the paper [17] by Nilsson et al. In summary
Puso, M A; Laursen, T A
2002-05-02
Smoothing of contact surfaces can be used to eliminate the chatter typically seen with node on facet contact and give a better representation of the actual contact surface. The latter affect is well demonstrated for problems with interference fits. In this work we present two methods for the smoothing of contact surfaces for 3D finite element contact. In the first method, we employ Gregory patches to smooth the faceted surface in a node on facet implementation. In the second method, we employ a Bezier interpolation of the faceted surface in a mortar method implementation of contact. As is well known, node on facet approaches can exhibit locking due to the failure of the Babuska-Brezzi condition and in some instances fail the patch test. The mortar method implementation is stable and provides optimal convergence in the energy of error. In the this work we demonstrate the superiority of the smoothed versus the non-smoothed node on facet implementations. We also show where the node on facet method fails and some results from the smoothed mortar method implementation.
Flux vector splitting of the inviscid equations with application to finite difference methods
NASA Technical Reports Server (NTRS)
Steger, J. L.; Warming, R. F.
1979-01-01
The conservation-law form of the inviscid gasdynamic equations has the remarkable property that the nonlinear flux vectors are homogeneous functions of degree one. This property readily permits the splitting of flux vectors into subvectors by similarity transformations so that each subvector has associated with it a specified eigenvalue spectrum. As a consequence of flux vector splitting, new explicit and implicit dissipative finite-difference schemes are developed for first-order hyperbolic systems of equations. Appropriate one-sided spatial differences for each split flux vector are used throughout the computational field even if the flow is locally subsonic. The results of some preliminary numerical computations are included.
Memory cost of absorbing conditions for the finite-difference time-domain method.
Chobeau, Pierre; Savioja, Lauri
2016-07-01
Three absorbing layers are investigated using standard rectilinear finite-difference schemes. The perfectly matched layer (PML) is compared with basic lossy layers terminated by two types of absorbing boundary conditions, all simulated using equivalent memory consumption. Lossy layers present the advantage of being scalar schemes, whereas the PML relies on a staggered scheme where both velocity and pressure are split. Although the PML gives the lowest reflection magnitudes over all frequencies and incidence angles, the most efficient lossy layer gives reflection magnitudes of the same order as the PML from mid- to high-frequency and for restricted incidence angles. PMID:27475200
NASA Technical Reports Server (NTRS)
Bailey, Harry E.; Beam, Richard M.
1991-01-01
Finite-difference approximations for steady-state compressible Navier-Stokes equations, whose two spatial dimensions are written in generalized curvilinear coordinates and strong conservation-law form, are presently solved by means of Newton's method in order to obtain a lifting-airfoil flow field under subsonic and transonnic conditions. In addition to ascertaining the computational requirements of an initial guess ensuring convergence and the degree of computational efficiency obtainable via the approximate Newton method's freezing of the Jacobian matrices, attention is given to the need for auxiliary methods assessing the temporal stability of steady-state solutions. It is demonstrated that nonunique solutions of the finite-difference equations are obtainable by Newton's method in conjunction with a continuation method.
NASA Technical Reports Server (NTRS)
Bauld, N. R., Jr.; Goree, J. G.; Tzeng, L.-S.
1985-01-01
It is pointed out that edge delamination is a serious failure mechanism for laminated composite materials. Various numerical methods have been utilized in attempts to calculate the interlaminar stress components which precede delamination in a laminate. There are, however, discrepancies regarding the results provided by different methods, taking into account a finite-difference procedure, a perturbation procedure, and finite element approaches. The present investigation has the objective to assess the capacity of a finite difference method to predict the character and magnitude of the interlaminar stress distributions near an interface corner. A second purpose of the investigation is to determine if predictions by finite element method in-plane, interlaminar stress components near an interface corner represent actual laminate behavior.
NASA Astrophysics Data System (ADS)
Russell, Greg; Harkins, Kevin D.; Secomb, Timothy W.; Galons, Jean-Philippe; Trouard, Theodore P.
2012-02-01
A new finite difference (FD) method for calculating the time evolution of complex transverse magnetization in diffusion-weighted magnetic resonance imaging and spectroscopy experiments is described that incorporates periodic boundary conditions. The new FD method relaxes restrictions on the allowable time step size employed in modeling which can significantly reduce computation time for simulations of large physical extent and allow for more complex, physiologically relevant, geometries to be simulated.
NASA Astrophysics Data System (ADS)
Li, Gang; Zhang, Lili; Hao, Tianyao
2016-02-01
An effective solver for the large complex system of linear equations is critical for improving the accuracy of numerical solutions in three-dimensional (3D) magnetotelluric (MT) modeling using the staggered finite-difference (SFD) method. In electromagnetic modeling, the formed system of linear equations is commonly solved using preconditioned iterative relaxation methods. We present 3D MT modeling using the SFD method, based on former work. The multigrid solver and three solvers preconditioned by incomplete Cholesky decomposition—the minimum residual method, the generalized product bi-conjugate gradient method and the bi-conjugate gradient stabilized method—are used to solve the formed system of linear equations. Divergence correction for the magnetic field is applied. We also present a comparison of the stability and convergence of these iterative solvers if divergence correction is used. Model tests show that divergence correction improves the convergence of iterative solvers and the accuracy of numerical results. Divergence correction can also decrease the number of iterations for fast convergence without changing the stability of linear solvers. For consideration of the computation time and memory requirements, the multigrid solver combined with divergence correction is preferred for 3D MT field simulation.
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.
A support-operator method for 3-D rupture dynamics
NASA Astrophysics Data System (ADS)
Ely, Geoffrey P.; Day, Steven M.; Minster, Jean-Bernard
2009-06-01
We present a numerical method to simulate spontaneous shear crack propagation within a heterogeneous, 3-D, viscoelastic medium. Wave motions are computed on a logically rectangular hexahedral mesh, using the generalized finite-difference method of Support Operators (SOM). This approach enables modelling of non-planar surfaces and non-planar fault ruptures. Our implementation, the Support Operator Rupture Dynamics (SORD) code, is highly scalable, enabling large-scale, multiprocessors calculations. The fault surface is modelled by coupled double nodes, where rupture occurs as dictated by the local stress conditions and a frictional failure law. The method successfully performs test problems developed for the Southern California Earthquake Center (SCEC)/U.S. Geological Survey (USGS) dynamic earthquake rupture code validation exercise, showing good agreement with semi-analytical boundary integral method results. We undertake further dynamic rupture tests to quantify numerical errors introduced by shear deformations to the hexahedral mesh. We generate a family of meshes distorted by simple shearing, in the along-strike direction, up to a maximum of 73°. For SCEC/USGS validation problem number 3, grid-induced errors increase with mesh shear angle, with the logarithm of error approximately proportional to angle over the range tested. At 73°, rms misfits are about 10 per cent for peak slip rate, and 0.5 per cent for both rupture time and total slip, indicating that the method (which, up to now, we have applied mainly to near-vertical strike-slip faulting) is also capable of handling geometries appropriate to low-angle surface-rupturing thrust earthquakes. Additionally, we demonstrate non-planar rupture effects, by modifying the test geometry to include, respectively, cylindrical curvature and sharp kinks.
Finite-difference time-domain methods to analyze ytterbium-doped Q-switched fiber lasers.
Hattori, Haroldo T; Khaleque, Abdul
2016-03-01
Q-switched lasers are widely used in material processing, laser ranging, medicine, and nonlinear optics--in particular, Q-switched lasers in optical fibers are important since they cannot only generate high peak powers but can also concentrate high peak powers in small areas. In this paper, we present new finite-difference time-domain methods that analyze the dynamics of Q-switched fiber lasers, which are more flexible and robust than previous methods. We extend the method to analyze fiber ring lasers and compare the results with our experiments. PMID:26974625
A Review of High-Order and Optimized Finite-Difference Methods for Simulating Linear Wave Phenomena
NASA Technical Reports Server (NTRS)
Zingg, David W.
1996-01-01
This paper presents a review of high-order and optimized finite-difference methods for numerically simulating the propagation and scattering of linear waves, such as electromagnetic, acoustic, or elastic waves. The spatial operators reviewed include compact schemes, non-compact schemes, schemes on staggered grids, and schemes which are optimized to produce specific characteristics. The time-marching methods discussed include Runge-Kutta methods, Adams-Bashforth methods, and the leapfrog method. In addition, the following fourth-order fully-discrete finite-difference methods are considered: a one-step implicit scheme with a three-point spatial stencil, a one-step explicit scheme with a five-point spatial stencil, and a two-step explicit scheme with a five-point spatial stencil. For each method studied, the number of grid points per wavelength required for accurate simulation of wave propagation over large distances is presented. Recommendations are made with respect to the suitability of the methods for specific problems and practical aspects of their use, such as appropriate Courant numbers and grid densities. Avenues for future research are suggested.
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.
Petersson, N. Anders; Sjogreen, Bjorn
2015-07-20
We develop a fourth order accurate finite difference method for solving the three-dimensional elastic wave equation in general heterogeneous anisotropic materials on curvilinear grids. The proposed method is an extension of the method for isotropic materials, previously described in the paper by Sjögreen and Petersson (2012) [11]. The method we proposed discretizes the anisotropic elastic wave equation in second order formulation, using a node centered finite difference method that satisfies the principle of summation by parts. The summation by parts technique results in a provably stable numerical method that is energy conserving. Also, we generalize and evaluate the super-grid far-field technique for truncating unbounded domains. Unlike the commonly used perfectly matched layers (PML), the super-grid technique is stable for general anisotropic material, because it is based on a coordinate stretching combined with an artificial dissipation. Moreover, the discretization satisfies an energy estimate, proving that the numerical approximation is stable. We demonstrate by numerical experiments that sufficiently wide super-grid layers result in very small artificial reflections. Applications of the proposed method are demonstrated by three-dimensional simulations of anisotropic wave propagation in crystals.
Petersson, N. Anders; Sjogreen, Bjorn
2015-07-20
We develop a fourth order accurate finite difference method for solving the three-dimensional elastic wave equation in general heterogeneous anisotropic materials on curvilinear grids. The proposed method is an extension of the method for isotropic materials, previously described in the paper by Sjögreen and Petersson (2012) [11]. The method we proposed discretizes the anisotropic elastic wave equation in second order formulation, using a node centered finite difference method that satisfies the principle of summation by parts. The summation by parts technique results in a provably stable numerical method that is energy conserving. Also, we generalize and evaluate the super-grid far-fieldmore » technique for truncating unbounded domains. Unlike the commonly used perfectly matched layers (PML), the super-grid technique is stable for general anisotropic material, because it is based on a coordinate stretching combined with an artificial dissipation. Moreover, the discretization satisfies an energy estimate, proving that the numerical approximation is stable. We demonstrate by numerical experiments that sufficiently wide super-grid layers result in very small artificial reflections. Applications of the proposed method are demonstrated by three-dimensional simulations of anisotropic wave propagation in crystals.« less
NASA Astrophysics Data System (ADS)
Tan, Sirui; Huang, Lianjie
2014-05-01
For modelling large-scale 3-D scalar-wave propagation, the finite-difference (FD) method with high-order accuracy in space but second-order accuracy in time is widely used because of its relatively low requirements of computer memory. We develop a novel staggered-grid (SG) FD method with high-order accuracy not only in space, but also in time, for solving 2- and 3-D scalar-wave equations. We determine the coefficients of the FD operator in the joint time-space domain to achieve high-order accuracy in time while preserving high-order accuracy in space. Our new FD scheme is based on a stencil that contains a few more grid points than the standard stencil. It is 2M-th-order accurate in space and fourth-order accurate in time when using 2M grid points along each axis and wavefields at one time step as the standard SGFD method. We validate the accuracy and efficiency of our new FD scheme using dispersion analysis and numerical modelling of scalar-wave propagation in 2- and 3-D complex models with a wide range of velocity contrasts. For media with a velocity contrast up to five, our new FD scheme is approximately two times more computationally efficient than the standard SGFD scheme with almost the same computer-memory requirement as the latter. Further numerical experiments demonstrate that our new FD scheme loses its advantages over the standard SGFD scheme if the velocity contrast is 10. However, for most large-scale geophysical applications, the velocity contrasts often range approximately from 1 to 3. Our new method is thus particularly useful for large-scale 3-D scalar-wave modelling and full-waveform inversion.
NASA Astrophysics Data System (ADS)
Maloney, James G.; Smith, Glenn S.; Scott, Waymond R., Jr.
1990-07-01
Two antennas are considered, a cylindrical monopole and a conical monopole. Both are driven through an image plane from a coaxial transmission line. Each of these antennas corresponds to a well-posed theoretical electromagnetic boundary value problem and a realizable experimental model. These antennas are analyzed by a straightforward application of the time-domain finite-difference method. The computed results for these antennas are shown to be in excellent agreement with accurate experimental measurements for both the time domain and the frequency domain. The graphical displays presented for the transient near-zone and far-zone radiation from these antennas provide physical insight into the radiation process.
NASA Technical Reports Server (NTRS)
Weatherill, W. H.; Ehlers, F. E.
1979-01-01
The design and usage of a pilot program for calculating the pressure distributions over harmonically oscillating airfoils in transonic flow are described. The procedure used is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equations for small disturbances. The steady velocity potential which must be obtained from some other program, was required for input. The unsteady equation, as solved, is linear with spatially varying coefficients. Since sinusoidal motion was assumed, time was not a variable. The numerical solution was obtained through a finite difference formulation and either a line relaxation or an out of core direct solution method.
Efficient solution on solving 3D Maxwell equations using stable semi-implicit splitting method
NASA Astrophysics Data System (ADS)
Cen, Wei; Gu, Ning
2016-05-01
In this paper, we propose an efficient solution on solving 3-dimensional (3D) time-domain Maxwell equations using the semi-implicit Crank-Nicholson (CN) method for time domain discretization with advantage of unconditional time stability. By applying the idea of fractional steps method (FSM) to the CN scheme, the proposed method provides a much simpler and efficient implementation than a direct implementation of the CN scheme. Compared with the alternating-direction implicit (ADI) method and explicit finite-difference time-domain approach (FDTD), it significantly saves the computational resource like memory and CPU time while remains similar numerical accuracy.
Beyond Euler's Method: Implicit Finite Differences in an Introductory ODE Course
ERIC Educational Resources Information Center
Kull, Trent C.
2011-01-01
A typical introductory course in ordinary differential equations (ODEs) exposes students to exact solution methods. However, many differential equations must be approximated with numerical methods. Textbooks commonly include explicit methods such as Euler's and Improved Euler's. Implicit methods are typically introduced in more advanced courses…
Second-order explicit finite-difference methods for transient-flow analysis
NASA Technical Reports Server (NTRS)
Chaudhry, M. H.; Hussaini, M. Y.
1983-01-01
Three second-order accurate numerical methods - MacCormack's method, Lambda scheme and Gabutti scheme - are introduced to solve the quasi-linear, hyperbolic partial differential equations describing transient flows in closed conduits. The details of these methods and the treatment of boundary conditions are presented and the results computed by using these methods for a typical piping system are compared. It is shown that for the same accuracy, second-order methods require considerably lesser number of computational nodes and computer time as compared to those required by the first-order methods.
Son, Sang-Kil
2011-03-01
We introduce a new numerical grid-based method on unstructured grids in the three-dimensional real-space to investigate the electronic structure of polyatomic molecules. The Voronoi-cell finite difference (VFD) method realizes a discrete Laplacian operator based on Voronoi cells and their natural neighbors, featuring high adaptivity and simplicity. To resolve multicenter Coulomb singularity in all-electron calculations of polyatomic molecules, this method utilizes highly adaptive molecular grids which consist of spherical atomic grids. It provides accurate and efficient solutions for the Schroedinger equation and the Poisson equation with the all-electron Coulomb potentials regardless of the coordinate system and the molecular symmetry. For numerical examples, we assess accuracy of the VFD method for electronic structures of one-electron polyatomic systems, and apply the method to the density-functional theory for many-electron polyatomic molecules.
A comparative study of finite-difference methods for radiative transfer problems.
NASA Astrophysics Data System (ADS)
Mohan Rao, D.; Varghese, B. A.; Srinivasa Rao, M.
1995-06-01
The authors have compared the numerical results of two widely used difference methods for the radiative transfer equation in plane-parallel medium. The Discrete Space theory (DS) is based on the direct first-order differential equation for the specific intensity whereas Auer's Hermitian (AH) method used the second order form for the mean-intensity and flux-like variables. The numerical results of these two methods are compared with analytical solutions under the two-stream approximation in a semi-infinite atmosphere. For the multi-stream case, the numerical errors are estimated using the solution of Chandrasekhar's discrete ordinate method.
Convergency analysis of the high-order mimetic finite difference method
Lipnikov, Konstantin; Veiga Da Beirao, L; Manzini, G
2008-01-01
We prove second-order convergence of the conservative variable and its flux in the high-order MFD method. The convergence results are proved for unstructured polyhedral meshes and full tensor diffusion coefficients. For the case of non-constant coefficients, we also develop a new family of high-order MFD methods. Theoretical result are confirmed through numerical experiments.
NASA Astrophysics Data System (ADS)
Kang, Hyun-Ju; Choi, Hyeok; Kim, Jae-Hyun; Lee, Se-Hun; Yoo, Tae-Ho; Chung, Chin-Wook
2016-06-01
A modified central difference method (MCDM) is proposed to obtain the electron energy distribution functions (EEDFs) in single Langmuir probes. Numerical calculation of the EEDF with MCDM is simple and has less noise. This method provides the second derivatives at a given point as the weighted average of second order central difference derivatives calculated at different voltage intervals, weighting each by the square of the interval. In this paper, the EEDFs obtained from MCDM are compared to those calculated via the averaged central difference method. It is found that MCDM effectively suppresses the noises in the EEDF, while the same number of points are used to calculate of the second derivative.
3D Face Modeling Using the Multi-Deformable Method
Hwang, Jinkyu; Yu, Sunjin; Kim, Joongrock; Lee, Sangyoun
2012-01-01
In this paper, we focus on the problem of the accuracy performance of 3D face modeling techniques using corresponding features in multiple views, which is quite sensitive to feature extraction errors. To solve the problem, we adopt a statistical model-based 3D face modeling approach in a mirror system consisting of two mirrors and a camera. The overall procedure of our 3D facial modeling method has two primary steps: 3D facial shape estimation using a multiple 3D face deformable model and texture mapping using seamless cloning that is a type of gradient-domain blending. To evaluate our method's performance, we generate 3D faces of 30 individuals and then carry out two tests: accuracy test and robustness test. Our method shows not only highly accurate 3D face shape results when compared with the ground truth, but also robustness to feature extraction errors. Moreover, 3D face rendering results intuitively show that our method is more robust to feature extraction errors than other 3D face modeling methods. An additional contribution of our method is that a wide range of face textures can be acquired by the mirror system. By using this texture map, we generate realistic 3D face for individuals at the end of the paper. PMID:23201976
Eulerian adaptive finite-difference method for high-velocity impact and penetration problems
Barton, Philip T.; Deiterding, Ralf; Meiron, Daniel I.; Pullin, Dale I
2013-01-01
Owing to the complex processes involved, faithful prediction of high-velocity impact events demands a simulation method delivering efficient calculations based on comprehensively formulated constitutive models. Such an approach is presented herein, employing a weighted essentially non-oscillatory (WENO) method within an adaptive mesh refinement (AMR) framework for the numerical solution of hyperbolic partial differential equations. Applied widely in computational fluid dynamics, these methods are well suited to the involved locally non-smooth finite deformations, circumventing any requirement for artificial viscosity functions for shock capturing. Application of the methods is facilitated through using a model of solid dynamics based upon hyper-elastic theory comprising kinematic evolution equations for the elastic distortion tensor. The model for finite inelastic deformations is phenomenologically equivalent to Maxwell s model of tangential stress relaxation. Closure relations tailored to the expected high-pressure states are proposed and calibrated for the materials of interest. Sharp interface resolution is achieved by employing level-set functions to track boundary motion, along with a ghost material method to capture the necessary internal boundary conditions for material interactions and stress-free surfaces. The approach is demonstrated for the simulation of high velocity impacts of steel projectiles on aluminium target plates in two and three dimensions.
Xiao, Jinbiao; Sun, Xiaohan
2012-09-10
A vector mode solver for bending waveguides by using a modified finite-difference (FD) method is developed in a local cylindrical coordinate system, where the perfectly matched layer absorbing boundary conditions are incorporated. Utilizing Taylor series expansion technique and continuity condition of the longitudinal field components, a standard matrix eigenvalue equation without the averaged index approximation approach for dealing with the discrete points neighboring the dielectric interfaces is obtained. Complex effective indexes and field distributions of leaky modes for a typical rib bending waveguide and a silicon wire bend are presented, and solutions accord well with those from the film mode matching method, which shows the validity and utility of the established method. PMID:23037277
NASA Technical Reports Server (NTRS)
Panczak, Tim; Ring, Steve; Welch, Mark
1999-01-01
Thermal engineering has long been left out of the concurrent engineering environment dominated by CAD (computer aided design) and FEM (finite element method) software. Current tools attempt to force the thermal design process into an environment primarily created to support structural analysis, which results in inappropriate thermal models. As a result, many thermal engineers either build models "by hand" or use geometric user interfaces that are separate from and have little useful connection, if any, to CAD and FEM systems. This paper describes the development of a new thermal design environment called the Thermal Desktop. This system, while fully integrated into a neutral, low cost CAD system, and which utilizes both FEM and FD methods, does not compromise the needs of the thermal engineer. Rather, the features needed for concurrent thermal analysis are specifically addressed by combining traditional parametric surface based radiation and FD based conduction modeling with CAD and FEM methods. The use of flexible and familiar temperature solvers such as SINDA/FLUINT (Systems Improved Numerical Differencing Analyzer/Fluid Integrator) is retained.
NASA Technical Reports Server (NTRS)
DeBonis, James R.
2013-01-01
A computational fluid dynamics code that solves the compressible Navier-Stokes equations was applied to the Taylor-Green vortex problem to examine the code s ability to accurately simulate the vortex decay and subsequent turbulence. The code, WRLES (Wave Resolving Large-Eddy Simulation), uses explicit central-differencing to compute the spatial derivatives and explicit Low Dispersion Runge-Kutta methods for the temporal discretization. The flow was first studied and characterized using Bogey & Bailley s 13-point dispersion relation preserving (DRP) scheme. The kinetic energy dissipation rate, computed both directly and from the enstrophy field, vorticity contours, and the energy spectra are examined. Results are in excellent agreement with a reference solution obtained using a spectral method and provide insight into computations of turbulent flows. In addition the following studies were performed: a comparison of 4th-, 8th-, 12th- and DRP spatial differencing schemes, the effect of the solution filtering on the results, the effect of large-eddy simulation sub-grid scale models, and the effect of high-order discretization of the viscous terms.
NASA Astrophysics Data System (ADS)
Steich, David James
1995-01-01
The Finite Difference Time Domain (FDTD) method is a simple yet powerful method for numerically solving electromagnetic wave phenomenon on computers. The FDTD technique discretizes Maxwell's equations with finite difference equations. These finite difference equations, which approximate local field behavior, are applied to large spatial lattices allowing calculation of a vast array of electromagnetical phenomenon. The greatest strengths of the FDTD method are in its simplicity, efficiency, and diversity. FDTD is capable of modeling the scattering and coupling to lossy dielectrics, lossy magnetics, anisotropic media, dispersive media, and nonlinear materials for general geometric shapes. Wideband frequency information can be obtained using FDTD for both near and far field observation points in a single computational run. However, along with all of its benefits, the FDTD algorithm has some deficiencies. For most problems of interest, poor accuracy at geometry interfaces of differing media and at outer problem space boundarys where the spatial lattice must be truncated are the two largest error sources of the FDTD algorithm. Although most accuracy issues can be circumvented by expending large amounts of computer memory and cpu time, using excessive computer resources is not always possible and is never appealing. The purpose of this thesis is to generalize, analyze, and test various mainstream local Outer Radiating Boundary Conditions (ORBCs) for the FDTD method applied to Maxwell's equations in order to help gain a better understanding of present ORBC limitations. A common mathematical model is presented for the boundary conditions. Boundary conditions shown to fit the model include Mur, Superabsorption, Liao, Higdon, and Lindman ORBCs of varying orders. Simple operators are defined and then used to generate the final discretized equations for each of the boundary conditions, automatically, without requiring complicated high order equations. The procedure also allows
A support-operator method for viscoelastic wave modelling in 3-D heterogeneous media
NASA Astrophysics Data System (ADS)
Ely, Geoffrey P.; Day, Steven M.; Minster, Jean-Bernard
2008-01-01
We apply the method of support operators (SOM) to solve the 3-D, viscoelastic equations of motion for use in earthquake simulations. SOM is a generalized finite-difference method that can utilize meshes of arbitrary structure and incorporate irregular geometry. Our implementation uses a 3-D, logically rectangular, hexahedral mesh. Calculations are second-order in space and time. A correction term is employed for suppression of spurious zero-energy modes (hourglass oscillations). We develop a free surface boundary condition, and an absorbing boundary condition using the method of perfectly matched layers (PML). Numerical tests using a layered material model in a highly deformed mesh show good agreement with the frequency-wavenumber method, for resolutions greater than 10 nodes per wavelength. We also test a vertically incident P wave on a semi-circular canyon, for which results match boundary integral solutions at resolutions greater that 20 nodes per wavelength. We also demonstrate excellent parallel scalability of our code.
NASA Technical Reports Server (NTRS)
Taflove, A.; Umashankar, K. R.
1987-01-01
The formulation and recent applications of the finite-difference time-domain (FD-TD) method for the numerical modeling of electromagnetic scattering and interaction problems are considered. It is shown that improvements in FD-TD modeling concepts and software implementation often make it a preferable choice for structures which cannot be easily treated by conventional integral equations and asymptotic approaches. Recent FD-TD modeling validations in research areas including coupling to wires and wire bundles in free space and cavities, scattering from surfaces in relativistic motion, inverse scattering, and radiation condition theory, are reviewed. Finally, the advantages and disadvantages of FD-TD, and guidelines concerning when FD-TD should and should not be used in high-frequency electromagnetic modeling problems, are summarized.
NASA Technical Reports Server (NTRS)
Abarbanel, S.; Gottlieb, D.
1976-01-01
The paper considers the leap-frog finite-difference method (Kreiss and Oliger, 1973) for systems of partial differential equations of the form du/dt = dF/dx + dG/dy + dH/dz, where d denotes partial derivative, u is a q-component vector and a function of x, y, z, and t, and the vectors F, G, and H are functions of u only. The original leap-frog algorithm is shown to admit a modification that improves on the stability conditions for two and three dimensions by factors of 2 and 2.8, respectively, thereby permitting larger time steps. The scheme for three dimensions is considered optimal in the sense that it combines simple averaging and large time steps.
NASA Astrophysics Data System (ADS)
Ding, H.; Shu, C.; Yeo, K. S.; Xu, D.
2007-01-01
In this paper, the mesh-free least square-based finite difference (MLSFD) method is applied to numerically study the flow field around two circular cylinders arranged in side-by-side and tandem configurations. For each configuration, various geometrical arrangements are considered, in order to reveal the different flow regimes characterized by the gap between the two cylinders. In this work, the flow simulations are carried out in the low Reynolds number range, that is, Re=100 and 200. Instantaneous vorticity contours and streamlines around the two cylinders are used as the visualization aids. Some flow parameters such as Strouhal number, drag and lift coefficients calculated from the solution are provided and quantitatively compared with those provided by other researchers.
3D modelling of the electromagnetic response of geophysical targets using the FDTD method
Debroux, P.S.
1996-05-01
A publicly available and maintained electromagnetic finite-difference time domain (FDTD) code has been applied to the forward modelling of the response of 1D, 2D and 3D geophysical targets to a vertical magnetic dipole excitation. The FDTD method is used to analyze target responses in the 1 MHz to 100 MHz range, where either conduction or displacement currents may have the controlling role. The response of the geophysical target to the excitation is presented as changes in the magnetic field ellipticity. The results of the FDTD code compare favorably with previously published integral equation solutions of the response of 1D targets, and FDTD models calculated with different finite-difference cell sizes are compared to find the effect of model discretization on the solution. The discretization errors, calculated as absolute error in ellipticity, are presented for the different ground geometry models considered, and are, for the most part, below 10% of the integral equation solutions. Finally, the FDTD code is used to calculate the magnetic ellipticity response of a 2D survey and a 3D sounding of complicated geophysical targets. The response of these 2D and 3D targets are too complicated to be verified with integral equation solutions, but show the proper low- and high-frequency responses.
a Fast Method for Measuring the Similarity Between 3d Model and 3d Point Cloud
NASA Astrophysics Data System (ADS)
Zhang, Zongliang; Li, Jonathan; Li, Xin; Lin, Yangbin; Zhang, Shanxin; Wang, Cheng
2016-06-01
This paper proposes a fast method for measuring the partial Similarity between 3D Model and 3D point Cloud (SimMC). It is crucial to measure SimMC for many point cloud-related applications such as 3D object retrieval and inverse procedural modelling. In our proposed method, the surface area of model and the Distance from Model to point Cloud (DistMC) are exploited as measurements to calculate SimMC. Here, DistMC is defined as the weighted distance of the distances between points sampled from model and point cloud. Similarly, Distance from point Cloud to Model (DistCM) is defined as the average distance of the distances between points in point cloud and model. In order to reduce huge computational burdens brought by calculation of DistCM in some traditional methods, we define SimMC as the ratio of weighted surface area of model to DistMC. Compared to those traditional SimMC measuring methods that are only able to measure global similarity, our method is capable of measuring partial similarity by employing distance-weighted strategy. Moreover, our method is able to be faster than other partial similarity assessment methods. We demonstrate the superiority of our method both on synthetic data and laser scanning data.
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.
3D scanning modeling method application in ancient city reconstruction
NASA Astrophysics Data System (ADS)
Ren, Pu; Zhou, Mingquan; Du, Guoguang; Shui, Wuyang; Zhou, Pengbo
2015-07-01
With the development of optical engineering technology, the precision of 3D scanning equipment becomes higher, and its role in 3D modeling is getting more distinctive. This paper proposed a 3D scanning modeling method that has been successfully applied in Chinese ancient city reconstruction. On one hand, for the existing architectures, an improved algorithm based on multiple scanning is adopted. Firstly, two pieces of scanning data were rough rigid registered using spherical displacers and vertex clustering method. Secondly, a global weighted ICP (iterative closest points) method is used to achieve a fine rigid registration. On the other hand, for the buildings which have already disappeared, an exemplar-driven algorithm for rapid modeling was proposed. Based on the 3D scanning technology and the historical data, a system approach was proposed for 3D modeling and virtual display of ancient city.
Simulations of P-SV wave scattering due to cracks by the 2-D finite difference method
NASA Astrophysics Data System (ADS)
Suzuki, Yuji; Shiina, Takahiro; Kawahara, Jun; Okamoto, Taro; Miyashita, Kaoru
2013-12-01
We simulate P-SV wave scattering by 2-D parallel cracks using the finite difference method (FDM). Here, special emphasis is put on simplicity; we apply a standard FDM (second-order velocity-stress scheme with a staggered grid) to media including traction-free, infinitesimally thin cracks, which are expressed in a simple manner. As an accuracy test of the present method, we calculate the displacement discontinuity along an isolated crack caused by harmonic waves using the method, which is compared with the corresponding results based on a reliable boundary integral equation method. The test resultantly indicates that the present method yields sufficient accuracy. As an application of this method, we also simulate wave propagation in media with randomly distributed cracks. We experimentally determine the attenuation and velocity dispersion induced by scattering from the synthetic seismograms, using a waveform averaging technique. It is shown that the results are well explained by a theory based on the Foldy approximation, if the crack density is sufficiently low. The theory appears valid with a crack density up to at least 0.1 for SV wave incidence, whereas the validity limit appears lower for P wave incidence.
NASA Astrophysics Data System (ADS)
Wang, Cheng; Dong, XinZhuang; Shu, Chi-Wang
2015-10-01
For numerical simulation of detonation, computational cost using uniform meshes is large due to the vast separation in both time and space scales. Adaptive mesh refinement (AMR) is advantageous for problems with vastly different scales. This paper aims to propose an AMR method with high order accuracy for numerical investigation of multi-dimensional detonation. A well-designed AMR method based on finite difference weighted essentially non-oscillatory (WENO) scheme, named as AMR&WENO is proposed. A new cell-based data structure is used to organize the adaptive meshes. The new data structure makes it possible for cells to communicate with each other quickly and easily. In order to develop an AMR method with high order accuracy, high order prolongations in both space and time are utilized in the data prolongation procedure. Based on the message passing interface (MPI) platform, we have developed a workload balancing parallel AMR&WENO code using the Hilbert space-filling curve algorithm. Our numerical experiments with detonation simulations indicate that the AMR&WENO is accurate and has a high resolution. Moreover, we evaluate and compare the performance of the uniform mesh WENO scheme and the parallel AMR&WENO method. The comparison results provide us further insight into the high performance of the parallel AMR&WENO method.
NASA Astrophysics Data System (ADS)
Wei, Xiao-Kun; Shao, Wei; Shi, Sheng-Bing; Zhang, Yong; Wang, Bing-Zhong
2015-07-01
An efficient conformal locally one-dimensional finite-difference time-domain (LOD-CFDTD) method is presented for solving two-dimensional (2D) electromagnetic (EM) scattering problems. The formulation for the 2D transverse-electric (TE) case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit (ADI) CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field (TF/SF) boundary and the perfectly matched layer (PML), the radar cross section (RCS) of two 2D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method. Project supported by the National Natural Science Foundation of China (Grant Nos. 61331007 and 61471105).
NASA Astrophysics Data System (ADS)
Ahn, Jai Seok
2014-01-01
A finite difference method (FDM) applicable to a two dimensional (2D) quantum dot was developed as a non-conventional approach to the theoretical understandings of quantum devices. This method can be applied to a realistic potential with an arbitrary shape. Using this method, the Hamiltonian in a tri-diagonal matrix could be obtained from any 2D potential, and the Hamiltonian could be diagonalized numerically for the eigenvalues. The legitimacy of this method was first checked by comparing the results with a finite round well with the analytic solutions. Two truncated harmonic wells were examined as a realistic model potential for lateral double quantum dots (DQDs) and for triple quantum dots (TQDs). The successful applications of the 2D FDM were observed with the entanglements in the DQDs. The level-splitting and anticrossing behaviors of the DQDs could be obtained by varying the distance between the dots and by introducing asymmetry in the well-depths. The 2D FDM results for linear/triangular TQDs were compared with the tight binding approximations.
NASA Astrophysics Data System (ADS)
MacKinnon, Robert J.; Carey, Graham F.
2003-01-01
A new class of positivity-preserving, flux-limited finite-difference and Petrov-Galerkin (PG) finite-element methods are devised for reactive transport problems.The methods are similar to classical TVD flux-limited schemes with the main difference being that the flux-limiter constraint is designed to preserve positivity for problems involving diffusion and reaction. In the finite-element formulation, we also consider the effect of numerical quadrature in the lumped and consistent mass matrix forms on the positivity-preserving property. Analysis of the latter scheme shows that positivity-preserving solutions of the resulting difference equations can only be guaranteed if the flux-limited scheme is both implicit and satisfies an additional lower-bound condition on time-step size. We show that this condition also applies to standard Galerkin linear finite-element approximations to the linear diffusion equation. Numerical experiments are provided to demonstrate the behavior of the methods and confirm the theoretical conditions on time-step size, mesh spacing, and flux limiting for transport problems with and without nonlinear reaction.
Efficient fabrication method of nano-grating for 3D holographic display with full parallax views.
Wan, Wenqiang; Qiao, Wen; Huang, Wenbin; Zhu, Ming; Fang, Zongbao; Pu, Donglin; Ye, Yan; Liu, Yanhua; Chen, Linsen
2016-03-21
Without any special glasses, multiview 3D displays based on the diffractive optics can present high resolution, full-parallax 3D images in an ultra-wide viewing angle. The enabling optical component, namely the phase plate, can produce arbitrarily distributed view zones by carefully designing the orientation and the period of each nano-grating pixel. However, such 3D display screen is restricted to a limited size due to the time-consuming fabricating process of nano-gratings on the phase plate. In this paper, we proposed and developed a lithography system that can fabricate the phase plate efficiently. Here we made two phase plates with full nano-grating pixel coverage at a speed of 20 mm^{2}/mins, a 500 fold increment in the efficiency when compared to the method of E-beam lithography. One 2.5-inch phase plate generated 9-view 3D images with horizontal-parallax, while the other 6-inch phase plate produced 64-view 3D images with full-parallax. The angular divergence in horizontal axis and vertical axis was 1.5 degrees, and 1.25 degrees, respectively, slightly larger than the simulated value of 1.2 degrees by Finite Difference Time Domain (FDTD). The intensity variation was less than 10% for each viewpoint, in consistency with the simulation results. On top of each phase plate, a high-resolution binary masking pattern containing amplitude information of all viewing zone was well aligned. We achieved a resolution of 400 pixels/inch and a viewing angle of 40 degrees for 9-view 3D images with horizontal parallax. In another prototype, the resolution of each view was 160 pixels/inch and the view angle was 50 degrees for 64-view 3D images with full parallax. As demonstrated in the experiments, the homemade lithography system provided the key fabricating technology for multiview 3D holographic display. PMID:27136814
NASA Technical Reports Server (NTRS)
Chang, Sin-Chung
1987-01-01
The validity of the modified equation stability analysis introduced by Warming and Hyett was investigated. It is shown that the procedure used in the derivation of the modified equation is flawed and generally leads to invalid results. Moreover, the interpretation of the modified equation as the exact partial differential equation solved by a finite-difference method generally cannot be justified even if spatial periodicity is assumed. For a two-level scheme, due to a series of mathematical quirks, the connection between the modified equation approach and the von Neuman method established by Warming and Hyett turns out to be correct despite its questionable original derivation. However, this connection is only partially valid for a scheme involving more than two time levels. In the von Neumann analysis, the complex error multiplication factor associated with a wave number generally has (L-1) roots for an L-level scheme. It is shown that the modified equation provides information about only one of these roots.
Garcia-Herranz, Nuria; Cabellos, Oscar; Aragones, Jose M.; Ahnert, Carol
2003-05-15
In order to take into account in a more effective and accurate way the intranodal heterogeneities in coarse-mesh finite-difference (CMFD) methods, a new equivalent parameter generation methodology has been developed and tested. This methodology accounts for the dependence of the nodal homogeneized two-group cross sections and nodal coupling factors, with interface flux discontinuity (IFD) factors that account for heterogeneities on the flux-spectrum and burnup intranodal distributions as well as on neighbor effects.The methodology has been implemented in an analytic CMFD method, rigorously obtained for homogeneous nodes with transverse leakage and generalized now for heterogeneous nodes by including IFD heterogeneity factors. When intranodal mesh node heterogeneity vanishes, the heterogeneous solution tends to the analytic homogeneous nodal solution. On the other hand, when intranodal heterogeneity increases, a high accuracy is maintained since the linear and nonlinear feedbacks on equivalent parameters have been shown to be as a very effective way of accounting for heterogeneity effects in two-group multidimensional coarse-mesh diffusion calculations.
Light Attenuation Method for 3D data acquisition (LAM3D) of bottom particle deposits
NASA Astrophysics Data System (ADS)
Er, Jenn Wei; Law, Adrian W. K.; Adams, E. Eric; Yang, Yang
2015-11-01
We have developed a novel experimental technique, Light Attenuation Method for 3D data acquisition (LAM3D), to acquire three-dimensional spatial characteristics and temporal development of bottom particle deposits. The new technique performs data acquisition with higher spatial and temporal resolution than existing approaches with laser and ultrasonic 3D profilers, and is therefore ideal for laboratory investigations with fast varying changes in the sediment bed, such as the developing deposition profile from sediment clouds commonly formed during dredging or land reclamation projects and the dynamic evolution in movable bed processes in rivers. The principle of the technique is based on the analysis of the light attenuation due to multiple light scattering through the particle deposits layer compared to the clear water column. With appropriate calibration, the particles size and distribution thickness can be quantified by the transmitted light spectrum. In the presentation, we will first show our measurement setup with a light panel for calibrated illumination and a system of DSLR cameras for the light capturing. Subsequently, we shall present the experimental results of fast evolving deposition profile of a barge-disposed sediment cloud upon its bottom impact on the sea bed.
A 3D Level Set Method for Microwave Breast Imaging
Colgan, Timothy J.; Hagness, Susan C.; Van Veen, Barry D.
2015-01-01
Objective Conventional inverse-scattering algorithms for microwave breast imaging result in moderate resolution images with blurred boundaries between tissues. Recent 2D numerical microwave imaging studies demonstrate that the use of a level set method preserves dielectric boundaries, resulting in a more accurate, higher resolution reconstruction of the dielectric properties distribution. Previously proposed level set algorithms are computationally expensive and thus impractical in 3D. In this paper we present a computationally tractable 3D microwave imaging algorithm based on level sets. Methods We reduce the computational cost of the level set method using a Jacobian matrix, rather than an adjoint method, to calculate Frechet derivatives. We demonstrate the feasibility of 3D imaging using simulated array measurements from 3D numerical breast phantoms. We evaluate performance by comparing full 3D reconstructions to those from a conventional microwave imaging technique. We also quantitatively assess the efficacy of our algorithm in evaluating breast density. Results Our reconstructions of 3D numerical breast phantoms improve upon those of a conventional microwave imaging technique. The density estimates from our level set algorithm are more accurate than those of conventional microwave imaging, and the accuracy is greater than that reported for mammographic density estimation. Conclusion Our level set method leads to a feasible level of computational complexity for full 3D imaging, and reconstructs the heterogeneous dielectric properties distribution of the breast more accurately than conventional microwave imaging methods. Significance 3D microwave breast imaging using a level set method is a promising low-cost, non-ionizing alternative to current breast imaging techniques. PMID:26011863
3D face recognition based on a modified ICP method
NASA Astrophysics Data System (ADS)
Zhao, Kankan; Xi, Jiangtao; Yu, Yanguang; Chicharo, Joe F.
2011-11-01
3D face recognition technique has gained much more attention recently, and it is widely used in security system, identification system, and access control system, etc. The core technique in 3D face recognition is to find out the corresponding points in different 3D face images. The classic partial Iterative Closest Point (ICP) method is iteratively align the two point sets based on repetitively calculate the closest points as the corresponding points in each iteration. After several iterations, the corresponding points can be obtained accurately. However, if two 3D face images with different scale are from the same person, the classic partial ICP does not work. In this paper we propose a modified partial Iterative Closest Point (ICP) method in which the scaling effect is considered to achieve 3D face recognition. We design a 3x3 diagonal matrix as the scale matrix in each iteration of the classic partial ICP. The probing face image which is multiplied by the scale matrix will keep the similar scale with the reference face image. Therefore, we can accurately determine the corresponding points even the scales of probing image and reference image are different. 3D face images in our experiments are acquired by a 3D data acquisition system based on Digital Fringe Projection Profilometry (DFPP). A 3D database consists of 30 group images, three images with the same scale, which are from the same person with different views, are included in each group. And in different groups, the scale of the 3 images may be different from other groups. The experiment results show that our proposed method can achieve 3D face recognition, especially in the case that the scales of probing image and referent image are different.
NASA Astrophysics Data System (ADS)
Moriyama, Eduardo H.; Zangaro, Renato A.; Lobo, Paulo D. d. C.; Villaverde, Antonio G. J. B.; Watanabe-Sei, Ii; Pacheco, Marcos T. T.; Otsuka, Daniel K.
2002-06-01
Thermal damage in dental pulp during Nd:YAG laser irradiation have been studied by several researchers; but due to dentin inhomogeneous structure, laser interaction with dentin in the hypersensitivity treatment are not fully understood. In this work, heat distribution profile on human dentine samples irradiated with Nd:YAG laser was simulated at surface and subjacent layers. Calculations were carried out using the Crank-Nicolson's finite difference method. Sixteen dentin samples with 1,5 mm of thickness were evenly distributed into four groups and irradiated with Nd:YAG laser pulses, according to the following scheme: (I) 1 pulse of 900 mJ, (II) 2 pulses of 450 mJ, (III) 3 pulses of 300 mJ, (IV) 6 pulses of 150 mJ; corresponding to a total laser energy of 900 mJ. The pulse interval was 300ms, the pulse duration of 900 ms and irradiated surface area of 0,005 mm2. Laser induced morphological changes in dentin were observed for all the irradiated samples. The heat distribution throughout the dentin layer, from the external dentin surface to the pulpal chamber wall, was calculated for each case, in order to obtain further information about the pulsed Nd:YAG laser-oral hard tissue interaction. The simulation showed significant differences in the final temperature at the pulpal chamber, depending on the exposition time and the energy contained in the laser pulse.
3-D UNSTRUCTURED HEXAHEDRAL-MESH Sn TRANSPORT METHODS
J. MOREL; J. MCGHEE; ET AL
2000-11-01
This is the final report of a three-year, Laboratory-Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). We have developed a method for solving the neutral-particle transport equation on 3-D unstructured hexahedral meshes using a S{sub n} discretization in angle in conjunction with a discontinuous finite-element discretization in space and a multigroup discretization in energy. Previous methods for solving this equation in 3-D have been limited to rectangular meshes. The unstructured-mesh method that we have developed is far more efficient for solving problems with complex 3-D geometric features than rectangular-mesh methods. In spite of having to make several compromises in our spatial discretization technique and our iterative solution technique, our method has been found to be both accurate and efficient for a broad class of problems.
[An integrated segmentation method for 3D ultrasound carotid artery].
Yang, Xin; Wu, Huihui; Liu, Yang; Xu, Hongwei; Liang, Huageng; Cai, Wenjuan; Fang, Mengjie; Wang, Yujie
2013-07-01
An integrated segmentation method for 3D ultrasound carotid artery was proposed. 3D ultrasound image was sliced into transverse, coronal and sagittal 2D images on the carotid bifurcation point. Then, the three images were processed respectively, and the carotid artery contours and thickness were obtained finally. This paper tries to overcome the disadvantages of current computer aided diagnosis method, such as high computational complexity, easily introduced subjective errors et al. The proposed method could get the carotid artery overall information rapidly, accurately and completely. It could be transplanted into clinical usage for atherosclerosis diagnosis and prevention. PMID:24195385
Improving automated 3D reconstruction methods via vision metrology
NASA Astrophysics Data System (ADS)
Toschi, Isabella; Nocerino, Erica; Hess, Mona; Menna, Fabio; Sargeant, Ben; MacDonald, Lindsay; Remondino, Fabio; Robson, Stuart
2015-05-01
This paper aims to provide a procedure for improving automated 3D reconstruction methods via vision metrology. The 3D reconstruction problem is generally addressed using two different approaches. On the one hand, vision metrology (VM) systems try to accurately derive 3D coordinates of few sparse object points for industrial measurement and inspection applications; on the other, recent dense image matching (DIM) algorithms are designed to produce dense point clouds for surface representations and analyses. This paper strives to demonstrate a step towards narrowing the gap between traditional VM and DIM approaches. Efforts are therefore intended to (i) test the metric performance of the automated photogrammetric 3D reconstruction procedure, (ii) enhance the accuracy of the final results and (iii) obtain statistical indicators of the quality achieved in the orientation step. VM tools are exploited to integrate their main functionalities (centroid measurement, photogrammetric network adjustment, precision assessment, etc.) into the pipeline of 3D dense reconstruction. Finally, geometric analyses and accuracy evaluations are performed on the raw output of the matching (i.e. the point clouds) by adopting a metrological approach. The latter is based on the use of known geometric shapes and quality parameters derived from VDI/VDE guidelines. Tests are carried out by imaging the calibrated Portable Metric Test Object, designed and built at University College London (UCL), UK. It allows assessment of the performance of the image orientation and matching procedures within a typical industrial scenario, characterised by poor texture and known 3D/2D shapes.
NASA Astrophysics Data System (ADS)
Wang, Enjiang; Liu, Yang; Sen, Mrinal K.
2016-07-01
The 2D acoustic wave equation is commonly solved numerically by finite-difference (FD) methods in which the accuracy of solution is significantly affected by the FD stencils. The commonly used cross stencil can reach either only second-order accuracy for space domain dispersion-relation-based FD method or (2 M)th-order accuracy along eight specific propagation directions for time-space domain dispersion-relation-based FD method, if the conventional (2 M)th-order spatial FD and second-order temporal FD are used to discretize the equation. One other newly developed rhombus stencil can reach arbitrary even-order accuracy. However, this stencil adds significantly computational cost when the operator length is large. To achieve a balance between the solution accuracy and efficiency, we develop a new FD stencil to solve the 2D acoustic wave equation. This stencil is a combination of the cross stencil and rhombus stencil. A cross stencil with an operator length parameter M is used to approximate the spatial partial derivatives while a rhombus stencil with an operator length parameter N together with the conventional 2nd-order temporal FD is employed in approximating the temporal partial derivatives. Using this stencil, a new FD scheme is developed; we demonstrate that this scheme can reach (2 M)th-order accuracy in space and (2 N)th-order accuracy in time when spatial FD coefficients and temporal FD coefficients are derived from respective dispersion relation using Taylor-series expansion (TE) method. To further increase the accuracy, we derive the FD coefficients by employing the time-space domain dispersion relation of this FD scheme using TE. We also use least-squares (LS) optimization method to reduce dispersion at high wavenumbers. Dispersion analysis, stability analysis and modelling examples demonstrate that our new scheme has greater accuracy and better stability than conventional FD schemes, and thus can adopt large time steps. To reduce the extra computational
A three-point backward finite-difference method has been derived for a system of mixed hyperbolic_{¯¯}parabolic (convection_{¯¯}diffusion) partial differential equations (mixed PDEs). The method resorts to the three-point backward differenci...
Klose, Alexander D.; Beattie, Bradley J.; Dehghani, Hamid; Vider, Lena; Le, Carl; Ponomarev, Vladimir; Blasberg, Ronald
2010-01-01
Purpose: Bioluminescence imaging is a research tool for studying gene expression levels in small animal models of human disease. Bioluminescence light, however, is strongly scattered in biological tissue and no direct image of the light-emitting reporter probe’s location can be obtained. Therefore, the authors have developed a linear image reconstruction method for bioluminescence tomography (BLT) that recovers the three-dimensional spatial bioluminescent source distribution in small animals. Methods: The proposed reconstruction method uses third-order simplified spherical harmonics (SP3) solutions to the equation of radiative transfer for modeling the bioluminescence light propagation in optically nonuniform tissue. The SP3 equations and boundary conditions are solved with a finite-difference (FD) technique on a regular grid. The curved geometry of the animal surface was taken into account with a blocking-off region method for regular grids. Coregistered computed tomography (CT) and magnetic resonance (MR) images provide information regarding the geometry of the skin surface and internal organs. The inverse source problem is defined as an algebraic system of linear equations for the unknown source distribution and is iteratively solved given multiview and multispectral boundary measurements. The average tissue absorption parameters, which are used for the image reconstruction process, were calculated with an evolution strategy (ES) from in vivo measurements using an implanted pointlike source of known location and spectrum. Moreover, anatomical information regarding the location of the internal organs and other tissue structures within the animal’s body are provided by coregistered MR images. Results: First, the authors recovered the wavelength-dependent absorption coefficients (average error of 14%) with the ES under ideal conditions by using a numerical mouse model. Next, they reconstructed the average absorption coefficient of a small animal by using an
MR image denoising method for brain surface 3D modeling
NASA Astrophysics Data System (ADS)
Zhao, De-xin; Liu, Peng-jie; Zhang, De-gan
2014-11-01
Three-dimensional (3D) modeling of medical images is a critical part of surgical simulation. In this paper, we focus on the magnetic resonance (MR) images denoising for brain modeling reconstruction, and exploit a practical solution. We attempt to remove the noise existing in the MR imaging signal and preserve the image characteristics. A wavelet-based adaptive curve shrinkage function is presented in spherical coordinates system. The comparative experiments show that the denoising method can preserve better image details and enhance the coefficients of contours. Using these denoised images, the brain 3D visualization is given through surface triangle mesh model, which demonstrates the effectiveness of the proposed method.
Shao, Yan-Lin Faltinsen, Odd M.
2014-10-01
We propose a new efficient and accurate numerical method based on harmonic polynomials to solve boundary value problems governed by 3D Laplace equation. The computational domain is discretized by overlapping cells. Within each cell, the velocity potential is represented by the linear superposition of a complete set of harmonic polynomials, which are the elementary solutions of Laplace equation. By its definition, the method is named as Harmonic Polynomial Cell (HPC) method. The characteristics of the accuracy and efficiency of the HPC method are demonstrated by studying analytical cases. Comparisons will be made with some other existing boundary element based methods, e.g. Quadratic Boundary Element Method (QBEM) and the Fast Multipole Accelerated QBEM (FMA-QBEM) and a fourth order Finite Difference Method (FDM). To demonstrate the applications of the method, it is applied to some studies relevant for marine hydrodynamics. Sloshing in 3D rectangular tanks, a fully-nonlinear numerical wave tank, fully-nonlinear wave focusing on a semi-circular shoal, and the nonlinear wave diffraction of a bottom-mounted cylinder in regular waves are studied. The comparisons with the experimental results and other numerical results are all in satisfactory agreement, indicating that the present HPC method is a promising method in solving potential-flow problems. The underlying procedure of the HPC method could also be useful in other fields than marine hydrodynamics involved with solving Laplace equation.
A method to fabricate disconnected silver nanostructures in 3D.
Vora, Kevin; Kang, SeungYeon; Mazur, Eric
2012-01-01
The standard nanofabrication toolkit includes techniques primarily aimed at creating 2D patterns in dielectric media. Creating metal patterns on a submicron scale requires a combination of nanofabrication tools and several material processing steps. For example, steps to create planar metal structures using ultraviolet photolithography and electron-beam lithography can include sample exposure, sample development, metal deposition, and metal liftoff. To create 3D metal structures, the sequence is repeated multiple times. The complexity and difficulty of stacking and aligning multiple layers limits practical implementations of 3D metal structuring using standard nanofabrication tools. Femtosecond-laser direct-writing has emerged as a pre-eminent technique for 3D nanofabrication.(1,2) Femtosecond lasers are frequently used to create 3D patterns in polymers and glasses.(3-7) However, 3D metal direct-writing remains a challenge. Here, we describe a method to fabricate silver nanostructures embedded inside a polymer matrix using a femtosecond laser centered at 800 nm. The method enables the fabrication of patterns not feasible using other techniques, such as 3D arrays of disconnected silver voxels.(8) Disconnected 3D metal patterns are useful for metamaterials where unit cells are not in contact with each other,(9) such as coupled metal dot(10,11)or coupled metal rod(12,13) resonators. Potential applications include negative index metamaterials, invisibility cloaks, and perfect lenses. In femtosecond-laser direct-writing, the laser wavelength is chosen such that photons are not linearly absorbed in the target medium. When the laser pulse duration is compressed to the femtosecond time scale and the radiation is tightly focused inside the target, the extremely high intensity induces nonlinear absorption. Multiple photons are absorbed simultaneously to cause electronic transitions that lead to material modification within the focused region. Using this approach, one can
NASA Astrophysics Data System (ADS)
Panayappan, Kadappan
With the advent of sub-micron technologies and increasing awareness of Electromagnetic Interference and Compatibility (EMI/EMC) issues, designers are often interested in full- wave solutions of complete systems, taking to account a variety of environments in which the system operates. However, attempts to do this substantially increase the complexities involved in computing full-wave solutions, especially when the problems involve multi- scale geometries with very fine features. For such problems, even the well-established numerical methods, such as the time domain technique FDTD and the frequency domain methods FEM and MoM, are often challenged to the limits of their capabilities. In an attempt to address such challenges, three novel techniques have been introduced in this work, namely Dipole Moment (DM) Approach, Recursive Update in Frequency Domain (RUFD) and New Finite Difference Time Domain ( vFDTD). Furthermore, the efficacy of the above techniques has been illustrated, via several examples, and the results obtained by proposed techniques have been compared with other existing numerical methods for the purpose of validation. The DM method is a new physics-based approach for formulating MoM problems, which is based on the use of dipole moments (DMs), as opposed to the conventional Green's functions. The absence of the Green's functions, as well as those of the vector and scalar potentials, helps to eliminate two of the key sources of difficulties in the conventional MoM formulation, namely the singularity and low-frequency problems. Specifically, we show that there are no singularities that we need to be concerned with in the DM formulation; hence, this obviates the need for special techniques for integrating these singularities. Yet another salutary feature of the DM approach is its ability to handle thin and lossy structures, or whether they are metallic, dielectric-type, or even combinations thereof. We have found that the DM formulation can handle these
Novel 3D Compression Methods for Geometry, Connectivity and Texture
NASA Astrophysics Data System (ADS)
Siddeq, M. M.; Rodrigues, M. A.
2016-06-01
A large number of applications in medical visualization, games, engineering design, entertainment, heritage, e-commerce and so on require the transmission of 3D models over the Internet or over local networks. 3D data compression is an important requirement for fast data storage, access and transmission within bandwidth limitations. The Wavefront OBJ (object) file format is commonly used to share models due to its clear simple design. Normally each OBJ file contains a large amount of data (e.g. vertices and triangulated faces, normals, texture coordinates and other parameters) describing the mesh surface. In this paper we introduce a new method to compress geometry, connectivity and texture coordinates by a novel Geometry Minimization Algorithm (GM-Algorithm) in connection with arithmetic coding. First, each vertex ( x, y, z) coordinates are encoded to a single value by the GM-Algorithm. Second, triangle faces are encoded by computing the differences between two adjacent vertex locations, which are compressed by arithmetic coding together with texture coordinates. We demonstrate the method on large data sets achieving compression ratios between 87 and 99 % without reduction in the number of reconstructed vertices and triangle faces. The decompression step is based on a Parallel Fast Matching Search Algorithm (Parallel-FMS) to recover the structure of the 3D mesh. A comparative analysis of compression ratios is provided with a number of commonly used 3D file formats such as VRML, OpenCTM and STL highlighting the performance and effectiveness of the proposed method.
NASA Astrophysics Data System (ADS)
Thomas, Justin W.
2006-12-01
The Numerical Nuclear Reactor (NNR) is a code suite that is being developed to provide high-fidelity multi-physics capability for the analysis of light water nuclear reactors. The focus of the work here is to extend the capability of the NNR by incorporation of the neutronics module, DeCART, for Boiling Water Reactor (BWR) applications. The DeCART code has been coupled to the NNR fluid mechanics and heat transfer module STAR-CD for light water reactor applications. The coupling has been accomplished via an interface program, which is responsible for mapping the STAR-CD and DeCART meshes, managing communication, and monitoring convergence. DeCART obtains the solution of the 3-D Boltzmann transport equation by performing a series of 2-D modular ray tracing-based method of characteristics problems that are coupled within the framework of 3-D coarse-mesh finite difference. The relatively complex geometry and increased axial heterogeneity found in BWRs are beyond the modeling capability of the original version of DeCART. In this work, DeCART is extended in three primary areas. First, the geometric capability is generalized by extending the modular ray tracing scheme and permitting an unstructured mesh in the global finite difference kernel. Second, numerical instabilities, which arose as a result of the severe axial heterogeneity found in BWR cores, have been resolved. Third, an advanced nodal method has been implemented to improve the accuracy of the axial flux distribution. In this semi-analytic nodal method, the analytic solution to the transverse-integrated neutron diffusion equation is obtained, where the nonhomogeneous neutron source was first approximated by a quartic polynomial. The successful completion of these three tasks has allowed the application of the coupled DeCART/STAR-CD code to practical BWR problems.
A Kosloff/Basal method, 3D migration program implemented on the CYBER 205 supercomputer
NASA Technical Reports Server (NTRS)
Pyle, L. D.; Wheat, S. R.
1984-01-01
Conventional finite difference migration has relied on approximations to the acoustic wave equation which allow energy to propagate only downwards. Although generally reliable, such approaches usually do not yield an accurate migration for geological structures with strong lateral velocity variations or with steeply dipping reflectors. An earlier study by D. Kosloff and E. Baysal (Migration with the Full Acoustic Wave Equation) examined an alternative approach based on the full acoustic wave equation. The 2D, Fourier type algorithm which was developed was tested by Kosloff and Baysal against synthetic data and against physical model data. The results indicated that such a scheme gives accurate migration for complicated structures. This paper describes the development and testing of a vectorized, 3D migration program for the CYBER 205 using the Kosloff/Baysal method. The program can accept as many as 65,536 zero offset (stacked) traces.
Color dithering methods for LEGO-like 3D printing
NASA Astrophysics Data System (ADS)
Sun, Pei-Li; Sie, Yuping
2015-01-01
Color dithering methods for LEGO-like 3D printing are proposed in this study. The first method is work for opaque color brick building. It is a modification of classic error diffusion. Many color primaries can be chosen. However, RGBYKW is recommended as its image quality is good and the number of color primary is limited. For translucent color bricks, multi-layer color building can enhance the image quality significantly. A LUT-based method is proposed to speed the dithering proceeding and make the color distribution even smoother. Simulation results show the proposed multi-layer dithering method can really improve the image quality of LEGO-like 3D printing.
SAMA: A Method for 3D Morphological Analysis
Cerruti, Florent; Sonnenschein, Carlos; Soto, Ana M.
2016-01-01
Three-dimensional (3D) culture models are critical tools for understanding tissue morphogenesis. A key requirement for their analysis is the ability to reconstruct the tissue into computational models that allow quantitative evaluation of the formed structures. Here, we present Software for Automated Morphological Analysis (SAMA), a method by which epithelial structures grown in 3D cultures can be imaged, reconstructed and analyzed with minimum human intervention. SAMA allows quantitative analysis of key features of epithelial morphogenesis such as ductal elongation, branching and lumen formation that distinguish different hormonal treatments. SAMA is a user-friendly set of customized macros operated via FIJI (http://fiji.sc/Fiji), an open-source image analysis platform in combination with a set of functions in R (http://www.r-project.org/), an open-source program for statistical analysis. SAMA enables a rapid, exhaustive and quantitative 3D analysis of the shape of a population of structures in a 3D image. SAMA is cross-platform, licensed under the GPLv3 and available at http://montevil.theobio.org/content/sama. PMID:27035711
NASA Technical Reports Server (NTRS)
Weatherill, W. H.; Sebastian, J. D.; Ehlers, F. E.
1978-01-01
Separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady equations for small disturbances was performed. The steady velocity potential was obtained first from the well known nonlinear equation for steady transonic flow. The unsteady velocity potential was then obtained from a linear differential equation in complex form with spatially varying coefficients. Since sinusoidal motion is assumed, the unsteady equation is independent of time. The results of an investigation into the relaxation-solution-instability problem was discussed. Concepts examined include variations in outer boundary conditions, a coordinate transformation so that the boundary condition at infinity may be applied to the outer boundaries of the finite difference region, and overlapping subregions. The general conclusion was that only a full direct solution in which all unknowns are obtained at the same time will avoid the solution instabilities of relaxation. An analysis of the one-dimensional form of the unsteady transonic equation was studied to evaluate errors between exact and finite difference solutions. Pressure distributions were presented for a low-aspect-ratio clipped delta wing at Mach number of 0.9 and for a moderate-aspect-ratio rectangular wing at a Mach number of 0.875.
A method for building 3D models of barchan dunes
NASA Astrophysics Data System (ADS)
Nai, Yang; Li-lan, Su; Lin, Wan; Jie, Yang; Shi-yi, Chen; Wei-lu, Hu
2016-01-01
The distributions of barchan dunes are usually represented by digital terrain models (DTMs) overlaid with digital orthophoto maps. Given that most regions with barchan dues have low relief, a 3D map obtained from a DTM may ineffectively show the stereoscopic shape of each dune. The method of building 3D models of barchan dunes using existing modeling software seldom considers the geographical environment. As a result, barchan dune models are often inconsistent with actual DTMs and incompletely express the morphological characteristics of dunes. Manual construction of barchan dune models is also costly and time consuming. Considering these problems, the morphological characteristics of barchan dunes and the mathematical relationships between the morphological parameters of the dunes, such as length, height, and width, are analyzed in this study. The methods of extracting the morphological feature points of barchan dunes, calculating their morphological parameters and building dune outlines and skeleton lines based on the medial axes, are also presented. The dune outlines, skeleton lines, and part of the medial axes of dunes are used to construct a constrained triangulated irregular network. C# and ArcEngine are employed to build 3D models of barchan dunes automatically. Experimental results of a study conducted in Tengger Desert show that the method can be used to approximate the morphological characteristics of barchan dunes and is less time consuming than manual methods.
Breast tumour visualization using 3D quantitative ultrasound methods
NASA Astrophysics Data System (ADS)
Gangeh, Mehrdad J.; Raheem, Abdul; Tadayyon, Hadi; Liu, Simon; Hadizad, Farnoosh; Czarnota, Gregory J.
2016-04-01
Breast cancer is one of the most common cancer types accounting for 29% of all cancer cases. Early detection and treatment has a crucial impact on improving the survival of affected patients. Ultrasound (US) is non-ionizing, portable, inexpensive, and real-time imaging modality for screening and quantifying breast cancer. Due to these attractive attributes, the last decade has witnessed many studies on using quantitative ultrasound (QUS) methods in tissue characterization. However, these studies have mainly been limited to 2-D QUS methods using hand-held US (HHUS) scanners. With the availability of automated breast ultrasound (ABUS) technology, this study is the first to develop 3-D QUS methods for the ABUS visualization of breast tumours. Using an ABUS system, unlike the manual 2-D HHUS device, the whole patient's breast was scanned in an automated manner. The acquired frames were subsequently examined and a region of interest (ROI) was selected in each frame where tumour was identified. Standard 2-D QUS methods were used to compute spectral and backscatter coefficient (BSC) parametric maps on the selected ROIs. Next, the computed 2-D parameters were mapped to a Cartesian 3-D space, interpolated, and rendered to provide a transparent color-coded visualization of the entire breast tumour. Such 3-D visualization can potentially be used for further analysis of the breast tumours in terms of their size and extension. Moreover, the 3-D volumetric scans can be used for tissue characterization and the categorization of breast tumours as benign or malignant by quantifying the computed parametric maps over the whole tumour volume.
Optical Sensors and Methods for Underwater 3D Reconstruction.
Massot-Campos, Miquel; Oliver-Codina, Gabriel
2015-01-01
This paper presents a survey on optical sensors and methods for 3D reconstruction in underwater environments. The techniques to obtain range data have been listed and explained, together with the different sensor hardware that makes them possible. The literature has been reviewed, and a classification has been proposed for the existing solutions. New developments, commercial solutions and previous reviews in this topic have also been gathered and considered. PMID:26694389
Optical Sensors and Methods for Underwater 3D Reconstruction
Massot-Campos, Miquel; Oliver-Codina, Gabriel
2015-01-01
This paper presents a survey on optical sensors and methods for 3D reconstruction in underwater environments. The techniques to obtain range data have been listed and explained, together with the different sensor hardware that makes them possible. The literature has been reviewed, and a classification has been proposed for the existing solutions. New developments, commercial solutions and previous reviews in this topic have also been gathered and considered. PMID:26694389
Discrete Method of Images for 3D Radio Propagation Modeling
NASA Astrophysics Data System (ADS)
Novak, Roman
2016-09-01
Discretization by rasterization is introduced into the method of images (MI) in the context of 3D deterministic radio propagation modeling as a way to exploit spatial coherence of electromagnetic propagation for fine-grained parallelism. Traditional algebraic treatment of bounding regions and surfaces is replaced by computer graphics rendering of 3D reflections and double refractions while building the image tree. The visibility of reception points and surfaces is also resolved by shader programs. The proposed rasterization is shown to be of comparable run time to that of the fundamentally parallel shooting and bouncing rays. The rasterization does not affect the signal evaluation backtracking step, thus preserving its advantage over the brute force ray-tracing methods in terms of accuracy. Moreover, the rendering resolution may be scaled back for a given level of scenario detail with only marginal impact on the image tree size. This allows selection of scene optimized execution parameters for faster execution, giving the method a competitive edge. The proposed variant of MI can be run on any GPU that supports real-time 3D graphics.
NASA Technical Reports Server (NTRS)
Collier, Richard S.
1997-01-01
This report describes finite difference computer calculations for the Space Shuttle Launch Pad which predict lightning induced electric currents and electric and magnetic fields caused by a lightning strike to the Lightning Protection System caternary wire. Description of possible lightning threats to Shuttle Payload components together with specifications for protection of these components, result from the calculation of lightning induced electric and magnetic fields inside and outside the during a lightning event. These fields also induce currents and voltages on cables and circuits which may be connected to, or a part of, shuttle payload components. These currents and voltages are also calculated. These threat levels are intended as a guide for designers of payload equipment to specify any shielding and/or lightning protection mitigation which may be required for payload components which are in the process of preparation or being transferred into the Shuttle Orbiter.
Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes
NASA Technical Reports Server (NTRS)
Feng, Huiyu; Mavriplis, Catherine; VanderWijngaart, Rob; Biswas, Rupak
2004-01-01
High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel
NASA Astrophysics Data System (ADS)
Sharkawi, K.-H.; Abdul-Rahman, A.
2013-09-01
to LoD4. The accuracy and structural complexity of the 3D objects increases with the LoD level where LoD0 is the simplest LoD (2.5D; Digital Terrain Model (DTM) + building or roof print) while LoD4 is the most complex LoD (architectural details with interior structures). Semantic information is one of the main components in CityGML and 3D City Models, and provides important information for any analyses. However, more often than not, the semantic information is not available for the 3D city model due to the unstandardized modelling process. One of the examples is where a building is normally generated as one object (without specific feature layers such as Roof, Ground floor, Level 1, Level 2, Block A, Block B, etc). This research attempts to develop a method to improve the semantic data updating process by segmenting the 3D building into simpler parts which will make it easier for the users to select and update the semantic information. The methodology is implemented for 3D buildings in LoD2 where the buildings are generated without architectural details but with distinct roof structures. This paper also introduces hybrid semantic-geometric 3D segmentation method that deals with hierarchical segmentation of a 3D building based on its semantic value and surface characteristics, fitted by one of the predefined primitives. For future work, the segmentation method will be implemented as part of the change detection module that can detect any changes on the 3D buildings, store and retrieve semantic information of the changed structure, automatically updates the 3D models and visualize the results in a userfriendly graphical user interface (GUI).
NASA Astrophysics Data System (ADS)
Choi, S.-J.; Giraldo, F. X.; Kim, J.; Shin, S.
2014-11-01
The non-hydrostatic (NH) compressible Euler equations for dry atmosphere were solved in a simplified two-dimensional (2-D) slice framework employing a spectral element method (SEM) for the horizontal discretization and a finite difference method (FDM) for the vertical discretization. By using horizontal SEM, which decomposes the physical domain into smaller pieces with a small communication stencil, a high level of scalability can be achieved. By using vertical FDM, an easy method for coupling the dynamics and existing physics packages can be provided. The SEM uses high-order nodal basis functions associated with Lagrange polynomials based on Gauss-Lobatto-Legendre (GLL) quadrature points. The FDM employs a third-order upwind-biased scheme for the vertical flux terms and a centered finite difference scheme for the vertical derivative and integral terms. For temporal integration, a time-split, third-order Runge-Kutta (RK3) integration technique was applied. The Euler equations that were used here are in flux form based on the hydrostatic pressure vertical coordinate. The equations are the same as those used in the Weather Research and Forecasting (WRF) model, but a hybrid sigma-pressure vertical coordinate was implemented in this model. We validated the model by conducting the widely used standard tests: linear hydrostatic mountain wave, tracer advection, and gravity wave over the Schär-type mountain, as well as density current, inertia-gravity wave, and rising thermal bubble. The results from these tests demonstrated that the model using the horizontal SEM and the vertical FDM is accurate and robust provided sufficient diffusion is applied. The results with various horizontal resolutions also showed convergence of second-order accuracy due to the accuracy of the time integration scheme and that of the vertical direction, although high-order basis functions were used in the horizontal. By using the 2-D slice model, we effectively showed that the combined spatial
System and method for 3D printing of aerogels
Worsley, Marcus A.; Duoss, Eric; Kuntz, Joshua; Spadaccini, Christopher; Zhu, Cheng
2016-03-08
A method of forming an aerogel. The method may involve providing a graphene oxide powder and mixing the graphene oxide powder with a solution to form an ink. A 3D printing technique may be used to write the ink into a catalytic solution that is contained in a fluid containment member to form a wet part. The wet part may then be cured in a sealed container for a predetermined period of time at a predetermined temperature. The cured wet part may then be dried to form a finished aerogel part.
Method and simulation to study 3D crosstalk perception
NASA Astrophysics Data System (ADS)
Khaustova, Dar'ya; Blondé, Laurent; Huynh-Thu, Quan; Vienne, Cyril; Doyen, Didier
2012-03-01
To various degrees, all modern 3DTV displays suffer from crosstalk, which can lead to a decrease of both visual quality and visual comfort, and also affect perception of depth. In the absence of a perfect 3D display technology, crosstalk has to be taken into account when studying perception of 3D stereoscopic content. In order to improve 3D presentation systems and understand how to efficiently eliminate crosstalk, it is necessary to understand its impact on human perception. In this paper, we present a practical method to study the perception of crosstalk. The approach consists of four steps: (1) physical measurements of a 3DTV, (2) building of a crosstalk surface based on those measurements and representing specifically the behavior of that 3TV, (3) manipulation of the crosstalk function and application on reference images to produce test images degraded by crosstalk in various ways, and (4) psychophysical tests. Our approach allows both a realistic representation of the behavior of a 3DTV and the easy manipulation of its resulting crosstalk in order to conduct psycho-visual experiments. Our approach can be used in all studies requiring the understanding of how crosstalk affects perception of stereoscopic content and how it can be corrected efficiently.
Reconstruction of 3D structure using stochastic methods: morphology and transport properties
NASA Astrophysics Data System (ADS)
Karsanina, Marina; Gerke, Kirill; Čapek, Pavel; Vasilyev, Roman; Korost, Dmitry; Skvortsova, Elena
2013-04-01
One of the main factors defining numerous flow phenomena in rocks, soils and other porous media, including fluid and solute movements, is pore structure, e.g., pore sizes and their connectivity. Numerous numerical methods were developed to quantify single and multi-phase flow in such media on microscale. Among most popular ones are: 1) a wide range of finite difference/element/volume solutions of Navier-Stokes equations and its simplifications; 2) lattice-Boltzmann method; 3) pore-network models, among others. Each method has some advantages and shortcomings, so that different research teams usually utilize more than one, depending on the study case. Recent progress in 3D imaging of internal structure, e.g., X-ray tomography, FIB-SEM and confocal microscopy, made it possible to obtain digitized input pore parameters for such models, however, a trade-off between resolution and sample size is usually unavoidable. There are situations then only standard two-dimensional information of porous structure is known due to tomography high cost or resolution limitations. However, physical modeling on microscale requires 3D information. There are three main approaches to reconstruct (using 2D cut(s) or some other limited information/properties) porous media: 1) statistical methods (correlation functions and simulated annealing, multi-point statistics, entropy methods), 2) sequential methods (sphere or other granular packs) and 3) morphological methods. Stochastic reconstructions using correlation functions possess some important advantage - they provide a statistical description of the structure, which is known to have relationships with all physical properties. In addition, this method is more flexible for other applications to characterize porous media. Taking different 3D scans of natural and artificial porous materials (sandstones, soils, shales, ceramics) we choose some 2D cut/s as sources of input correlation functions. Based on different types of correlation functions
Single-camera fixed perspective 360-deg 3D method
NASA Astrophysics Data System (ADS)
Harding, Kevin G.; Fergan, Robert K.
1997-01-01
The use of 3D methods for such applications as feature locations within a wide field-of-view, such as for automated guided vehicles or large assembly work, offers some distinct challenges. The use of stereo viewing has often been the method of choice due to the wide area coverage and hardware simplicity. However, stereo based methods suffer from a loss of spatial position resolution for more distant object as compared to close objects due to the high demagnification needed to cover large fields-of-view. A long depth-of-field in such systems may also degrade the general ability to perform correlations due to poor focus. In addition, stereo looses distance resolution for features nearing the line of the two cameras, typically requiring movement of the cameras. The paper presents a novel method of obtaining 3D scene information as seen from the center of a cylindrical field. The method described uses a single camera with a view that is rotated through 360 degrees by means of a continuously rotating mirror. The viewing systems uses a constant field of view optical system that provides a constant X-Y resolution of features in the scene over depths of several meters. Comparing successive images with the readout from an encoder on the rotating mirror generates all locations of objects within a limited height cylinder. This paper will discuss the sources of errors and typical capabilities of this approach in light of a real-time part location tracking application useful in assembly systems.
Quantitative validation of the 3D SAR profile of hyperthermia applicators using the gamma method.
de Bruijne, Maarten; Samaras, Theodoros; Chavannes, Nicolas; van Rhoon, Gerard C
2007-06-01
For quality assurance of hyperthermia treatment planning systems, quantitative validation of the electromagnetic model of an applicator is essential. The objective of this study was to validate a finite-difference time-domain (FDTD) model implementation of the Lucite cone applicator (LCA) for superficial hyperthermia. The validation involved (i) the assessment of the match between the predicted and measured 3D specific absorption rate (SAR) distribution, and (ii) the assessment of the ratio between model power and real-world power. The 3D SAR distribution of seven LCAs was scanned in a phantom bath using the DASY4 dosimetric measurement system. The same set-up was modelled in SEMCAD X. The match between the predicted and the measured SAR distribution was quantified with the gamma method, which combines distance-to-agreement and dose difference criteria. Good quantitative agreement was observed: more than 95% of the measurement points met the acceptance criteria 2 mm/2% for all applicators. The ratio between measured and predicted power absorption ranged from 0.75 to 0.92 (mean 0.85). This study shows that quantitative validation of hyperthermia applicator models is feasible and is worth considering as a part of hyperthermia quality assurance procedures. PMID:17505090
3D reconstruction methods of coronal structures by radio observations
NASA Astrophysics Data System (ADS)
Aschwanden, Markus J.; Bastian, T. S.; White, Stephen M.
1992-11-01
The ability to carry out the three dimensional (3D) reconstruction of structures in the solar corona would represent a major advance in the study of the physical properties in active regions and in flares. Methods which allow a geometric reconstruction of quasistationary coronal structures (for example active region loops) or dynamic structures (for example flaring loops) are described: stereoscopy of multi-day imaging observations by the VLA (Very Large Array); tomography of optically thin emission (in radio or soft x-rays); multifrequency band imaging by the VLA; and tracing of magnetic field lines by propagating electron beams.
3D reconstruction methods of coronal structures by radio observations
NASA Technical Reports Server (NTRS)
Aschwanden, Markus J.; Bastian, T. S.; White, Stephen M.
1992-01-01
The ability to carry out the three dimensional (3D) reconstruction of structures in the solar corona would represent a major advance in the study of the physical properties in active regions and in flares. Methods which allow a geometric reconstruction of quasistationary coronal structures (for example active region loops) or dynamic structures (for example flaring loops) are described: stereoscopy of multi-day imaging observations by the VLA (Very Large Array); tomography of optically thin emission (in radio or soft x-rays); multifrequency band imaging by the VLA; and tracing of magnetic field lines by propagating electron beams.
NASA Astrophysics Data System (ADS)
Bogdanovich, B. Yu.; Nesterovich, A. V.; Sukhanova, L. A.; Khlestkov, Yu. A.
2014-08-01
Results of modeling of a high-current relativistic beam by the finite-difference method are compared with results obtained for a beam with the same parameters using the well-known KARAT code, which is based on the large-particle method. These two methods give similar results, which justifies the use of the finite-difference method for the numerical solution of the equations of motion describing the motion of the beam in its own and an external electromagnetic field.
NASA Astrophysics Data System (ADS)
Wang, Zhi-liang; Li, Yong-chi; Wang, J. G.
2006-12-01
The propagation and attenuation of blast-induced stress waves differs between geomedia such as rock or soil mass. This paper numerically studies the propagation and attenuation of blast-induced elastoplastic waves in deep geomedia by using a one-dimensional (1-D) finite-difference code. Firstly, the elastoplastic Cap models for rock and soil masses are introduced into the governing equations of spherical wave motion and a FORTRAN code based on the finite difference method is developed. Secondly, an underground spherical blast is simulated with this code and verified by software, RENEWTO. The propagation of stress-waves in rock and soil masses is numerically investigated, respectively. Finally, the effect of a soil cover layer on the attenuation of stress waves in the rear rock mass is studied. It is determined that large plastic deformation of geomedia can effectively dissipate the energy of stress-waves inward and the developed 1-D finite difference code coupled with elastoplastic Cap models is convenient and effective in the numerical simulations for underground spherical explosion.
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.
A perceptual preprocess method for 3D-HEVC
NASA Astrophysics Data System (ADS)
Shi, Yawen; Wang, Yongfang; Wang, Yubing
2015-08-01
A perceptual preprocessing method for 3D-HEVC coding is proposed in the paper. Firstly we proposed a new JND model, which accounts for luminance contrast masking effect, spatial masking effect, and temporal masking effect, saliency characteristic as well as depth information. We utilize spectral residual approach to obtain the saliency map and built a visual saliency factor based on saliency map. In order to distinguish the sensitivity of objects in different depth. We segment each texture frame into foreground and background by a automatic threshold selection algorithm using corresponding depth information, and then built a depth weighting factor. A JND modulation factor is built with a linear combined with visual saliency factor and depth weighting factor to adjust the JND threshold. Then, we applied the proposed JND model to 3D-HEVC for residual filtering and distortion coefficient processing. The filtering process is that the residual value will be set to zero if the JND threshold is greater than residual value, or directly subtract the JND threshold from residual value if JND threshold is less than residual value. Experiment results demonstrate that the proposed method can achieve average bit rate reduction of 15.11%, compared to the original coding scheme with HTM12.1, while maintains the same subjective quality.
Lee, D. W.; Joo, H. G.
2013-07-01
The methods and performance of a three-dimensional S{sub n} transport code employing the Discontinuous Finite Element Method (DFEM) and the Coarse Mesh Finite Difference (CMFD) formulation are presented. The mesh generator GMSH and a post processing visualization tool Visit are combined with the code for flexible geometry processing and versatile visualization. The CMFD method for DFEM Sn applications is formulated and the performance of the CMFD acceleration of eigenvalue calculations is demonstrated for a simple set of neutron transport problems. (authors)
NASA Astrophysics Data System (ADS)
Bettaibi, Soufiene; Kuznik, Frédéric; Sediki, Ezeddine
2016-02-01
This paper presents a numerical study of thermosolutal mixed convection in rectangular enclosure with sliding top lid. The fluid flow is solved by the multiple relaxation time (MRT) lattice Boltzmann method (LBM), whereas the temperature and concentration fields are computed by finite difference method (FDM). The main objective of this study is to investigate the accuracy and the effectiveness of such model to predict thermodynamics for heat and mass transfer in a driven cavity. This model is validated with different numerical methods in the current literature. A good agreement is obtained between our results and previous works. The different comparisons demonstrate the robustness and the accuracy of the proposed approach.
3D Hot Test Simulations of a 220 GHz Folded Waveguide Traveling Wave Tube Using a CFDTD PIC Method
NASA Astrophysics Data System (ADS)
Lin, Ming-Chieh; Song, Heather
2015-11-01
Millimeter or sub-THz wave sources centered at 220 GHz is of interest due to the potential for its commercial and military applications including high resolution radar, remote sensing, and high-data-rate communications. It has been demonstrated via 3D cold test finite element method (FEM) simulations that a folded waveguide traveling wave tube (FWTWT) can be designed and optimized at this frequency range with a small signal gain of 18 dB over a comparatively broad (-3 dB) bandwidth of ~ 10%. On the other hand, 3D hot test simulations of a V-band ladder TWT have been successfully demonstrated using a conformal finite-difference time-domain (CFDTD) particle-in-cell (PIC) method for center frequency of 50 GHz. In the present work, the 220 GHz FWTWT designs have been reviewed and studied. 3D Cold test simulations using both the CFDTD and FEM methods have been carried out and compared with each other as basis for 3D hot test CFDTD PIC simulations. The preliminary simulation result shows that the gain-bandwidth features at 220 GHz are achievable while carefully avoiding beam interceptions. Our study shows that the interaction characteristics are very sensitive to the operating beam parameters. Detail simulation results and discussions will be presented.
A method of PSF generation for 3D brightfield deconvolution.
Tadrous, P J
2010-02-01
This paper addresses the problem of 3D deconvolution of through focus widefield microscope datasets (Z-stacks). One of the most difficult stages in brightfield deconvolution is finding the point spread function. A theoretically calculated point spread function (called a 'synthetic PSF' in this paper) requires foreknowledge of many system parameters and still gives only approximate results. A point spread function measured from a sub-resolution bead suffers from low signal-to-noise ratio, compounded in the brightfield setting (by contrast to fluorescence) by absorptive, refractive and dispersal effects. This paper describes a method of point spread function estimation based on measurements of a Z-stack through a thin sample. This Z-stack is deconvolved by an idealized point spread function derived from the same Z-stack to yield a point spread function of high signal-to-noise ratio that is also inherently tailored to the imaging system. The theory is validated by a practical experiment comparing the non-blind 3D deconvolution of the yeast Saccharomyces cerevisiae with the point spread function generated using the method presented in this paper (called the 'extracted PSF') to a synthetic point spread function. Restoration of both high- and low-contrast brightfield structures is achieved with fewer artefacts using the extracted point spread function obtained with this method. Furthermore the deconvolution progresses further (more iterations are allowed before the error function reaches its nadir) with the extracted point spread function compared to the synthetic point spread function indicating that the extracted point spread function is a better fit to the brightfield deconvolution model than the synthetic point spread function. PMID:20096049
The COMET method in 3-D hexagonal geometry
Connolly, K. J.; Rahnema, F.
2012-07-01
The hybrid stochastic-deterministic coarse mesh radiation transport (COMET) method developed at Georgia Tech now solves reactor core problems in 3-D hexagonal geometry. In this paper, the method is used to solve three preliminary test problems designed to challenge the method with steep flux gradients, high leakage, and strong asymmetry and heterogeneity in the core. The test problems are composed of blocks taken from a high temperature test reactor benchmark problem. As the method is still in development, these problems and their results are strictly preliminary. Results are compared to whole core Monte Carlo reference solutions in order to verify the method. Relative errors are on the order of 50 pcm in core eigenvalue, and mean relative error in pin fission density calculations is less than 1% in these difficult test cores. The method requires the one-time pre-computation of a response expansion coefficient library, which may be compiled in a comparable amount of time to a single whole core Monte Carlo calculation. After the library has been computed, COMET may solve any number of core configurations on the order of an hour, representing a significant gain in efficiency over other methods for whole core transport calculations. (authors)
NASA Astrophysics Data System (ADS)
Miranda, D. D.; Howard, A. Q.
2012-12-01
Computational modelling of geophysical data is an important step in the process of hydrocarbon exploration. It consists in simulating the exploratory procedure and realistic geological environments. It allows a preliminary evaluation of the exploration feasibility of a particular terrain or geological model, indicating the best conditions for geophysical surveys. In this paper, we assess the Finite Difference frequency domain method for modelling the electromagnetic response of a horizontal electric dipole in 1D and 2.5D geometries. The non-uniform grid is refined in regions where the electromagnetic fields vary rapidly, namely the regions where we have variation in conductivity distribution and near the source dipole. We chose the horizontal electromagnetic dipole because it is the source normally used in the marine controlled-source electromagnetic surveys (mCSEM), which is the next step in our research. The mCSEM, also known as Sea Bed Logging, is a method for detection and characterization of thin resistive structures, like hydrocarbon reservoirs, often located in regions of deep water. It consists of a mobile electric dipole or a magnetic loop as a source, positioned near the sea floor where an array of electric and magnetic receivers are deployed. The source transmitter uses a low frequency signal on the order of 1Hz, that diffuses both in the ocean and in the sediments beneath it and is captured by the receivers . Amplitude and phase of this signal depend on the electrical conductivity of the seabed environment. The complexity of the environments and the large dimensions of the geological domains that we want to investigate make the modelling procedure extremely demanding, since the Finite Difference method requires a total discretization of the studied domain, resulting in large systems of linear equations, which can make the procedure long and expensive. Non-uniform grids and exploitation of the sparse property of the Finite Difference matrices are example
The 3D inelastic analysis methods for hot section components
NASA Technical Reports Server (NTRS)
Dame, L. T.; Mcknight, R. L.
1983-01-01
The objective of this research is to develop an analytical tool capable of economically evaluating the cyclic time dependent plasticity which occurs in hot section engine components in areas of strain concentration resulting from the combination of both mechanical and thermal stresses. The techniques developed must be capable of accommodating large excursions in temperatures with the associated variations in material properties including plasticity and creep. The overall objective of this proposed program is to develop advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate and yet more cost effective analysis of combustors, turbine blades, and vanes. The approach will be to develop four different theories, one linear and three higher order with increasing complexities including embedded singularities.
On 3D inelastic analysis methods for hot section components
NASA Technical Reports Server (NTRS)
Mcknight, R. L.; Chen, P. C.; Dame, L. T.; Holt, R. V.; Huang, H.; Hartle, M.; Gellin, S.; Allen, D. H.; Haisler, W. E.
1986-01-01
Accomplishments are described for the 2-year program, to develop advanced 3-D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades and vanes. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulations models were developed; an eight-noded mid-surface shell element, a nine-noded mid-surface shell element and a twenty-noded isoparametric solid element. A separate computer program was developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.
The 3D inelastic analysis methods for hot section components
NASA Technical Reports Server (NTRS)
Mcknight, R. L.; Maffeo, R. J.; Tipton, M. T.; Weber, G.
1992-01-01
A two-year program to develop advanced 3D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades, and vanes is described. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulation models were developed: an eight-noded midsurface shell element; a nine-noded midsurface shell element; and a twenty-noded isoparametric solid element. A separate computer program has been developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.
Improved time-space method for 3-D heat transfer problems including global warming
Saitoh, T.S.; Wakashima, Shinichiro
1999-07-01
In this paper, the Time-Space Method (TSM) which has been proposed for solving general heat transfer and fluid flow problems was improved in order to cover global and urban warming. The TSM is effective in almost all-transient heat transfer and fluid flow problems, and has been already applied to the 2-D melting problems (or moving boundary problems). The computer running time will be reduced to only 1/100th--1/1000th of the existing schemes for 2-D and 3-D problems. However, in order to apply to much larger-scale problems, for example, global warming, urban warming and general ocean circulation, the SOR method (or other iterative methods) in four dimensions is somewhat tedious and provokingly slow. Motivated by the above situation, the authors improved the speed of iteration of the previous TSM by introducing the following ideas: (1) Timewise chopping: Time domain is chopped into small peaches to save memory requirement; (2) Adaptive iteration: Converged region is eliminated for further iteration; (3) Internal selective iteration: Equation with slow iteration speed in iterative procedure is selectively iterated to accelerate entire convergence; and (4) False transient integration: False transient term is added to the Poisson-type equation and the relevant solution is regarded as a parabolic equation. By adopting the above improvements, the higher-order finite different schemes and the hybrid mesh, the computer running time for the TSM is reduced to some 1/4600th of the conventional explicit method for a typical 3-D natural convection problem in a closed cavity. The proposed TSM will be more efficacious for large-scale environmental problems, such as global warming, urban warming and general ocean circulation, in which a tremendous computing time would be required.
Lattice Boltzmann Method for 3-D Flows with Curved Boundary
NASA Technical Reports Server (NTRS)
Mei, Renwei; Shyy, Wei; Yu, Dazhi; Luo, Li-Shi
2002-01-01
In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.
3D Wavelet-Based Filter and Method
Moss, William C.; Haase, Sebastian; Sedat, John W.
2008-08-12
A 3D wavelet-based filter for visualizing and locating structural features of a user-specified linear size in 2D or 3D image data. The only input parameter is a characteristic linear size of the feature of interest, and the filter output contains only those regions that are correlated with the characteristic size, thus denoising the image.
Merritt, M.L.
1993-01-01
The simulation of the transport of injected freshwater in a thin brackish aquifer, overlain and underlain by confining layers containing more saline water, is shown to be influenced by the choice of the finite-difference approximation method, the algorithm for representing vertical advective and dispersive fluxes, and the values assigned to parametric coefficients that specify the degree of vertical dispersion and molecular diffusion that occurs. Computed potable water recovery efficiencies will differ depending upon the choice of algorithm and approximation method, as will dispersion coefficients estimated based on the calibration of simulations to match measured data. A comparison of centered and backward finite-difference approximation methods shows that substantially different transition zones between injected and native waters are depicted by the different methods, and computed recovery efficiencies vary greatly. Standard and experimental algorithms and a variety of values for molecular diffusivity, transverse dispersivity, and vertical scaling factor were compared in simulations of freshwater storage in a thin brackish aquifer. Computed recovery efficiencies vary considerably, and appreciable differences are observed in the distribution of injected freshwater in the various cases tested. The results demonstrate both a qualitatively different description of transport using the experimental algorithms and the interrelated influences of molecular diffusion and transverse dispersion on simulated recovery efficiency. When simulating natural aquifer flow in cross-section, flushing of the aquifer occurred for all tested coefficient choices using both standard and experimental algorithms. ?? 1993.
Study on Low-Frequency Oscillations in a Gyrotron Using a 3D CFDTD PIC Method
NASA Astrophysics Data System (ADS)
Lin, M. C.; Smithe, D. N.
2010-11-01
Low-frequency oscillations (LFOs) have been observed in a high average power gyrotron and the trapped electron population contributing to the oscillation has been measured. As high average power gyrotrons are the most promising millimeter wave source for thermonuclear fusion research, it is important to get a better understanding of this parasitic phenomenon to avoid any deterioration of the electron beam quality thus reducing the gyrotron efficiency. 2D Particle-in-cell simulations quasi-statically model the development of oscillations of the space charge in the adiabatic trap, but the physics of the electron dynamics in the adiabatic trap is only partially understood. Therefore, understanding of the LFOs remains incomplete and a full picture of this parasitic phenomenon has not been seen yet. In this work, we use a 3D conformal finite-difference time-domain (CFDTD) particle-in-cell (PIC) method to accurately and efficiently study the LFOs in a high average power gyrotron. As the CFDTD method exhibits a second order accuracy, complicated structures, such as a magnetron injection gun, can be well described. Employing a highly parallelized computation, the model can be simulated in time domain more realistically.
NASA Astrophysics Data System (ADS)
De Basabe, Jonás D.; Sen, Mrinal K.
2015-01-01
The numerical simulation of wave propagation in media with solid and fluid layers is essential for marine seismic exploration data analysis. The numerical methods for wave propagation that are applicable to this physical settings can be broadly classified as partitioned or monolithic: The partitioned methods use separate simulations in the fluid and solid regions and explicitly satisfy the interface conditions, whereas the monolithic methods use the same method in all the domain without any special treatment of the fluid-solid interface. Despite the accuracy of the partitioned methods, the monolithic methods are more common in practice because of their convenience. In this paper, we analyse the accuracy of several monolithic methods for wave propagation in the presence of a fluid-solid interface. The analysis is based on grid-dispersion criteria and numerical examples. The methods studied here include: the classical finite-difference method (FDM) based on the second-order displacement formulation of the elastic wave equation (DFDM), the staggered-grid finite difference method (SGFDM), the velocity-stress FDM with a standard grid (VSFDM) and the spectral-element method (SEM). We observe that among these, DFDM and the first-order SEM have a large amount of grid dispersion in the fluid region which renders them impractical for this application. On the other hand, SGFDM, VSFDM and SEM of order greater or equal to 2 yield accurate results for the body waves in the fluid and solid regions if a sufficient number of nodes per wavelength is used. All of the considered methods yield limited accuracy for the surface waves because the proper boundary conditions are not incorporated into the numerical scheme. Overall, we demonstrate both by analytic treatment and numerical experiments, that a first-order velocity-stress formulation can, in general, be used in dealing with fluid-solid interfaces without using staggered grids necessarily.
Methods for Geometric Data Validation of 3d City Models
NASA Astrophysics Data System (ADS)
Wagner, D.; Alam, N.; Wewetzer, M.; Pries, M.; Coors, V.
2015-12-01
Geometric quality of 3D city models is crucial for data analysis and simulation tasks, which are part of modern applications of the data (e.g. potential heating energy consumption of city quarters, solar potential, etc.). Geometric quality in these contexts is however a different concept as it is for 2D maps. In the latter case, aspects such as positional or temporal accuracy and correctness represent typical quality metrics of the data. They are defined in ISO 19157 and should be mentioned as part of the metadata. 3D data has a far wider range of aspects which influence their quality, plus the idea of quality itself is application dependent. Thus, concepts for definition of quality are needed, including methods to validate these definitions. Quality on this sense means internal validation and detection of inconsistent or wrong geometry according to a predefined set of rules. A useful starting point would be to have correct geometry in accordance with ISO 19107. A valid solid should consist of planar faces which touch their neighbours exclusively in defined corner points and edges. No gaps between them are allowed, and the whole feature must be 2-manifold. In this paper, we present methods to validate common geometric requirements for building geometry. Different checks based on several algorithms have been implemented to validate a set of rules derived from the solid definition mentioned above (e.g. water tightness of the solid or planarity of its polygons), as they were developed for the software tool CityDoctor. The method of each check is specified, with a special focus on the discussion of tolerance values where they are necessary. The checks include polygon level checks to validate the correctness of each polygon, i.e. closeness of the bounding linear ring and planarity. On the solid level, which is only validated if the polygons have passed validation, correct polygon orientation is checked, after self-intersections outside of defined corner points and edges
3D simulation of seismic wave propagation around a tunnel using the spectral element method
NASA Astrophysics Data System (ADS)
Lambrecht, L.; Friederich, W.
2010-05-01
We model seismic wave propagation in the environment of a tunnel for later application to reconnaissance. Elastic wave propagation can be simulated by different numerical techniques such as finite differences and pseudospectral methods. Their disadvantage is the lack of accuracy on free surfaces, numerical dispersion and inflexibility of the mesh. Here we use the software package SPECFEM3D_SESAME in an svn development version, which is based on the spectral element method (SEM) and can handle complex mesh geometries. A weak form of the elastic wave equation leads to a linear system of equations with a diagonal mass matrix, where the free surface boundary of the tunnel can be treated under realistic conditions and can be effectively implemented in parallel. We have designed a 3D external mesh including a tunnel and realistic features such as layers and holes to simulate elastic wave propagation in the zone around the tunnel. The source is acting at the tunnel surface so that we excite Rayleigh waves which propagate to the front face of the tunnel. A conversion takes place and a high amplitude S-wave is radiated in the direction of the tunnel axis. Reflections from perturbations in front of the tunnel can be measured by receivers implemented on the tunnel face. For a shallow tunnel the land surface has high influence on the wave propagation. By implementing additional receivers at this surface we intent to improve the prediction. It shows that the SEM is very capable to handle the complex geometry of the model and especially incorporates the free surfaces of the model.