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.
Benchmarks of 3D Laplace Equation Solvers in a Cubic Configuration for Streamer Simulation
NASA Astrophysics Data System (ADS)
Joseph-Marie, Plewa; Olivier, Ducasse; Philippe, Dessante; Carolyn, Jacobs; Olivier, Eichwald; Nicolas, Renon; Mohammed, Yousfi
2016-05-01
The aim of this paper is to test a developed SOR R&B method using the Chebyshev accelerator algorithm to solve the Laplace equation in a cubic 3D configuration. Comparisons are made in terms of precision and computing time with other elliptic equation solvers proposed in the open source LIS library. The first results, obtained by using a single core on a HPC, show that the developed SOR R&B method is efficient when the spectral radius needed for the Chebyshev acceleration is carefully pre-estimated. Preliminary results obtained with a parallelized code using the MPI library are also discussed when the calculation is distributed over one hundred cores.
Blumberg, L.N.
1992-03-01
The authors have analyzed simulated magnetic measurements data for the SXLS bending magnet in a plane perpendicular to the reference axis at the magnet midpoint by fitting the data to an expansion solution of the 3-dimensional Laplace equation in curvilinear coordinates as proposed by Brown and Servranckx. The method of least squares is used to evaluate the expansion coefficients and their uncertainties, and compared to results from an FFT fit of 128 simulated data points on a 12-mm radius circle about the reference axis. They find that the FFT method gives smaller coefficient uncertainties that the Least Squares method when the data are within similar areas. The Least Squares method compares more favorably when a larger number of data points are used within a rectangular area of 30-mm vertical by 60-mm horizontal--perhaps the largest area within the 35-mm x 75-mm vacuum chamber for which data could be obtained. For a grid with 0.5-mm spacing within the 30 x 60 mm area the Least Squares fit gives much smaller uncertainties than the FFT. They are therefore in the favorable position of having two methods which can determine the multipole coefficients to much better accuracy than the tolerances specified to General Dynamics. The FFT method may be preferable since it requires only one Hall probe rather than the four envisioned for the least squares grid data. However least squares can attain better accuracy with fewer probe movements. The time factor in acquiring the data will likely be the determining factor in choice of method. They should further explore least squares analysis of a Fourier expansion of data on a circle or arc of a circle since that method gives coefficient uncertainties without need for multiple independent sets of data as needed by the FFT method.
NASA Astrophysics Data System (ADS)
Brüning, J.; Dobrokhotov, S. Yu.; Minenkov, D. S.
2011-12-01
The aim of this paper is to construct solutions of the Dirichlet problem for the 3D Laplace equation in a layer with highly oscillating boundary. The boundary simulates the surface of a nanotube array, and the solutions are applied to compute the cold field electron emission. We suggest a family of exact solutions that solve the problem for a boundary with appropriate geometry. These solutions, along with the Fowler-Nordheim formula, allow one to present explicit asymptotic formulas for the electric field and the emission current. In this part of the paper, we consider the main mathematical aspects, restricting ourselves to the analysis of properties of the potential created by a single tube and a regular array of tubes. In the next part, we shall consider some cases corresponding to nonregular arrays of tubes and concrete physical examples.
Laplace-domain waveform modeling and inversion for the 3D acoustic-elastic coupled media
NASA Astrophysics Data System (ADS)
Shin, Jungkyun; Shin, Changsoo; Calandra, Henri
2016-06-01
Laplace-domain waveform inversion reconstructs long-wavelength subsurface models by using the zero-frequency component of damped seismic signals. Despite the computational advantages of Laplace-domain waveform inversion over conventional frequency-domain waveform inversion, an acoustic assumption and an iterative matrix solver have been used to invert 3D marine datasets to mitigate the intensive computing cost. In this study, we develop a Laplace-domain waveform modeling and inversion algorithm for 3D acoustic-elastic coupled media by using a parallel sparse direct solver library (MUltifrontal Massively Parallel Solver, MUMPS). We precisely simulate a real marine environment by coupling the 3D acoustic and elastic wave equations with the proper boundary condition at the fluid-solid interface. In addition, we can extract the elastic properties of the Earth below the sea bottom from the recorded acoustic pressure datasets. As a matrix solver, the parallel sparse direct solver is used to factorize the non-symmetric impedance matrix in a distributed memory architecture and rapidly solve the wave field for a number of shots by using the lower and upper matrix factors. Using both synthetic datasets and real datasets obtained by a 3D wide azimuth survey, the long-wavelength component of the P-wave and S-wave velocity models is reconstructed and the proposed modeling and inversion algorithm are verified. A cluster of 80 CPU cores is used for this study.
Laplace's equation and Faraday's lines of force
Narasimhan, T.N.
2007-06-01
Boundary-value problems involve two dependent variables: a potential function, and a stream function. They can be approached in two mutually independent ways. The first, introduced by Laplace, involves spatial gradients at a point. Inspired by Faraday, Maxwell introduced the other, visualizing the flow domain as a collection of flow tubes and isopotential surfaces. Boundary-value problems intrinsically entail coupled treatment (or, equivalently, optimization) of potential and stream functions Historically, potential theory avoided the cumbersome optimization task through ingenious techniques such as conformal mapping and Green's functions. Laplace's point-based approach, and Maxwell's global approach, each provides its own unique insights into boundary-value problems. Commonly, Laplace's equation is solved either algebraically, or with approximate numerical methods. Maxwell's geometry-based approach opens up novel possibilities of direct optimization, providing an independent logical basis for numerical models, rather than treating them as approximate solvers of the differential equation. Whereas points, gradients, and Darcy's law are central to posing problems on the basis of Laplace's approach, flow tubes, potential differences, and the mathematical form of Ohm's law are central to posing them in natural coordinates oriented along flow paths. Besides being of philosophical interest, optimization algorithms can provide advantages that complement the power of classical numerical models. In the spirit of Maxwell, who eloquently spoke for a balance between abstract mathematical symbolism and observable attributes of concrete objects, this paper is an examination of the central ideas of the two approaches, and a reflection on how Maxwell's integral visualization may be practically put to use in a world of digital computers.
The Geometric Solution of Laplace's Equation
NASA Astrophysics Data System (ADS)
Bakhoum, Ezzat Gamal
In 1891, J.J. Thomson--the discoverer of the electron--stated a formula that relates the first derivative of the electric field intensity to the mean curvature of an equipotential surface. That formula was later proved by others, but remained unexploited in any practical purpose to this date. This dissertation presents a numerical method based on Thomson's formula for the rapid solution of Laplace's equation, the governing equation of field theory. The presented method is based on geometric construction principles. Specifically, the method uses the concept of representing equipotential surfaces by polynomials for the rapid tracing of these surfaces; and is therefore fundamentally different from previously-known techniques which are based on discretizing the domain or the boundary of the problem. The new method is especially suited for problems which have complicated or irregular boundaries as well as problems in exterior domains. Previously, such types of problems have required a number of computations of O(N.M), where N is the number of points taken on the boundary of the problem and M is the number of points inside the domain at which the solution is to be computed. The new method requires an O(M) computations only; and is therefore significantly faster than the previous techniques. Applications include problems of electrostatics, cosmology, biomedical engineering, nuclear and particle physics, etc.
3D Laplace-domain full waveform inversion using a single GPU card
NASA Astrophysics Data System (ADS)
Shin, Jungkyun; Ha, Wansoo; Jun, Hyunggu; Min, Dong-Joo; Shin, Changsoo
2014-06-01
The Laplace-domain full waveform inversion is an efficient long-wavelength velocity estimation method for seismic datasets lacking low-frequency components. However, to invert a 3D velocity model, a large cluster of CPU cores have commonly been required to overcome the extremely long computing time caused by a large impedance matrix and a number of source positions. In this study, a workstation with a single GPU card (NVIDIA GTX 580) is successfully used for the 3D Laplace-domain full waveform inversion rather than a large cluster of CPU cores. To exploit a GPU for our inversion algorithm, the routine for the iterative matrix solver is ported to the CUDA programming language for forward and backward modeling parts with minimized modification of the remaining parts, which were originally written in Fortran 90. Using a uniformly structured grid set, nonzero values in the sparse impedance matrix can be arranged according to certain rules, which efficiently parallelize the preconditioned conjugate gradient method for a number of threads contained in the GPU card. We perform a numerical experiment to verify the accuracy of a floating point operation performed by a GPU to calculate the Laplace-domain wavefield. We also measure the efficiencies of the original CPU and modified GPU programs using a cluster of CPU cores and a workstation with a GPU card, respectively. Through the analysis, the parallelized inversion code for a GPU achieves the speedup of 14.7-24.6x compared to a CPU-based serial code depending on the degrees of freedom of the impedance matrix. Finally, the practicality of the proposed algorithm is examined by inverting a 3D long-wavelength velocity model using wide azimuth real datasets in 3.7 days.
Laplace's equation on convex polyhedra via the unified method
Ashton, A. C. L.
2015-01-01
We provide a new method to study the classical Dirichlet problem for Laplace's equation on a convex polyhedron. This new approach was motivated by Fokas’ unified method for boundary value problems. The central object in this approach is the global relation: an integral equation which couples the known boundary data and the unknown boundary values. This integral equation depends holomorphically on two complex parameters, and the resulting analysis takes place on a Banach space of complex analytic functions closely related to the classical Paley–Wiener space. We write the global relation in the form of an operator equation and prove that the relevant operator is bounded below using some novel integral identities. We give a new integral representation to the solution to the underlying boundary value problem which serves as a concrete realization of the fundamental principle of Ehrenpreis.
Open active cloaking and illusion devices for the Laplace equation
NASA Astrophysics Data System (ADS)
Ma, Qian; Yang, Fan; Jin, Tian Yu; Lei Mei, Zhong; Cui, Tie Jun
2016-04-01
We propose open active cloaking and illusion devices for the Laplace equation. Compared with the closed configurations of active cloaking and illusion devices, we focus on improving the distribution schemes for the controlled sources, which do not have to surround the protected object strictly. Instead, the controlled sources can be placed in several small discrete clusters, and produce the desired voltages along the controlled boundary, to actively hide or disguise the protected object. Numerical simulations are performed with satisfactory results, which are further validated by experimental measurements. The open cloaking and illusion devices have many advantages over the closed configurations in various potential applications.
Decoupling of the DGLAP evolution equations by Laplace method
NASA Astrophysics Data System (ADS)
Boroun, G. R.; Zarrin, S.; Teimoury, F.
2015-10-01
In this paper we derive two second-order differential equations for the gluon and singlet distribution functions by using the Laplace transform method. We decoupled the solutions of the singlet and gluon distributions into the initial conditions (function and derivative of the function) at the virtuality Q 0 2 separately as these solutions are defined by F 2 s ( x, Q 2) = F( F 0 ,∂ F s0 and G( x, Q 2)= G( G 0, ∂ G 0. We compared our results with the MSTW parameterization and the experimental measurements of F 2 p ( x, Q 2.
Effectiveness of the Young-Laplace equation at nanoscale
Liu, Hailong; Cao, Guoxin
2016-01-01
Using molecular dynamics (MD) simulations, a new approach based on the behavior of pressurized water out of a nanopore (1.3–2.7 nm) in a flat plate is developed to calculate the relationship between the water surface curvature and the pressure difference across water surface. It is found that the water surface curvature is inversely proportional to the pressure difference across surface at nanoscale, and this relationship will be effective for different pore size, temperature, and even for electrolyte solutions. Based on the present results, we cannot only effectively determine the surface tension of water and the effects of temperature or electrolyte ions on the surface tension, but also show that the Young-Laplace (Y-L) equation is valid at nanoscale. In addition, the contact angle of water with the hydrophilic material can be further calculated by the relationship between the critical instable pressure of water surface (burst pressure) and nanopore size. Combining with the infiltration behavior of water into hydrophobic microchannels, the contact angle of water at nanoscale can be more accurately determined by measuring the critical pressure causing the instability of water surface, based on which the uncertainty of measuring the contact angle of water at nanoscale is highly reduced. PMID:27033874
Effectiveness of the Young-Laplace equation at nanoscale
NASA Astrophysics Data System (ADS)
Liu, Hailong; Cao, Guoxin
2016-04-01
Using molecular dynamics (MD) simulations, a new approach based on the behavior of pressurized water out of a nanopore (1.3–2.7 nm) in a flat plate is developed to calculate the relationship between the water surface curvature and the pressure difference across water surface. It is found that the water surface curvature is inversely proportional to the pressure difference across surface at nanoscale, and this relationship will be effective for different pore size, temperature, and even for electrolyte solutions. Based on the present results, we cannot only effectively determine the surface tension of water and the effects of temperature or electrolyte ions on the surface tension, but also show that the Young-Laplace (Y-L) equation is valid at nanoscale. In addition, the contact angle of water with the hydrophilic material can be further calculated by the relationship between the critical instable pressure of water surface (burst pressure) and nanopore size. Combining with the infiltration behavior of water into hydrophobic microchannels, the contact angle of water at nanoscale can be more accurately determined by measuring the critical pressure causing the instability of water surface, based on which the uncertainty of measuring the contact angle of water at nanoscale is highly reduced.
Laplace homotopy perturbation method for Burgers equation with space- and time-fractional order
NASA Astrophysics Data System (ADS)
Johnston, S. J.; Jafari, H.; Moshokoa, S. P.; Ariyan, V. M.; Baleanu, D.
2016-07-01
The fractional Burgers equation describes the physical processes of unidirectional propagation of weakly nonlinear acoustic waves through a gas-filled pipe. The Laplace homotopy perturbation method is discussed to obtain the approximate analytical solution of space-fractional and time-fractional Burgers equations. The method used combines the Laplace transform and the homotopy perturbation method. Numerical results show that the approach is easy to implement and accurate when applied to partial differential equations of fractional orders.
NASA Technical Reports Server (NTRS)
Truong, K. V.; Unal, Aynur; Tobak, M.
1989-01-01
Various features of the solutions of Duffing's equation are described using a representation of the solutions in the Laplace-Borel transform domain. An application of this technique is illustrated for the symmetry-breaking bifurcation of a hard spring.
Discovering the Laplace Transform in Undergraduate Differential Equations
ERIC Educational Resources Information Center
Quinn, Terrance J.; Rai, Sanjay
2008-01-01
The Laplace Transform is an object of fundamental importance in pure and applied mathematics. In addition, it has special pedagogical value in that it can provide a natural and concrete setting for a student to begin thinking about the modern concepts of "operator" and "functional". Most undergraduate textbooks, however, merely define the…
NASA Astrophysics Data System (ADS)
Petrov, P.; Newman, G. A.
2010-12-01
-Fourier domain we had developed 3D code for full-wave field simulation in the elastic media which take into account nonlinearity introduced by free-surface effects. Our approach is based on the velocity-stress formulation. In the contrast to conventional formulation we defined the material properties such as density and Lame constants not at nodal points but within cells. This second order finite differences method formulated in the cell-based grid, generate numerical solutions compatible with analytical ones within the range errors determinate by dispersion analysis. Our simulator will be embedded in an inversion scheme for joint seismic- electromagnetic imaging. It also offers possibilities for preconditioning the seismic wave propagation problems in the frequency domain. References. Shin, C. & Cha, Y. (2009), Waveform inversion in the Laplace-Fourier domain, Geophys. J. Int. 177(3), 1067- 1079. Shin, C. & Cha, Y. H. (2008), Waveform inversion in the Laplace domain, Geophys. J. Int. 173(3), 922-931. Commer, M. & Newman, G. (2008), New advances in three-dimensional controlled-source electromagnetic inversion, Geophys. J. Int. 172(2), 513-535. Newman, G. A., Commer, M. & Carazzone, J. J. (2010), Imaging CSEM data in the presence of electrical anisotropy, Geophysics, in press.
Design method for electromagnetic cloak with arbitrary shapes based on Laplace's equation.
Hu, Jin; Zhou, Xiaoming; Hu, Gengkai
2009-02-01
In transformation optics, the space transformation is viewed as the deformation of a material. The permittivity and permeability tensors in the transformed space are found to correlate with the deformation field of the material. By solving the Laplace's equation, which describes how the material will deform during a transformation, we can design electromagnetic cloaks with arbitrary shapes if the boundary conditions of the cloak are considered. As examples, the material parameters of the spherical and elliptical cylindrical cloaks are derived based on the analytical solutions of the Laplace's equation. For cloaks with irregular shapes, the material parameters of the transformation medium are determined numerically by solving the Laplace's equation. Full-wave simulations based on the Maxwell's equations validate the designed cloaks. The proposed method can be easily extended to design other transformation materials for electromagnetic and acoustic wave phenomena. PMID:19188959
Solution of the Time-Dependent Schrödinger Equation by the Laplace Transform Method
Lin, S. H.; Eyring, H.
1971-01-01
The time-dependent Schrödinger equation for two quite general types of perturbation has been solved by introducing the Laplace transforms to eliminate the time variable. The resulting time-independent differential equation can then be solved by the perturbation method, the variation method, the variation-perturbation method, and other methods. PMID:16591898
Derivation of the generalized Young-Laplace equation of curved interfaces in nanoscaled solids
NASA Astrophysics Data System (ADS)
Chen, Tungyang; Chiu, Min-Sen; Weng, Chung-Ning
2006-10-01
In nanoscaled solids, the mathematical behavior of a curved interface between two different phases with interface stress effects can be described by the generalized Young-Laplace equations [T. Young, Philos. Trans. R. Soc. London 95, 65 (1805); P. S. Laplace, Traite de Mechanique Celeste (Gauthier-Villars, Paris, 1805), Vol. 4, Supplements au Livre X]. Here we present a geometric illustration to prove the equations. By considering a small element of the curved thin interface, we model the interface stresses as in-plane stresses acting along its edges, while on the top and bottom faces of the interface the tractions are contributed from its three-dimensional bulk neighborhood. With this schematic illustration, simple force balance considerations will give the Young-Laplace equations across the interface. Similar procedures can be applied to conduction phenomena. This will allow us to reconstruct one type of imperfect interfaces, referred to as highly conducting interfaces.
On the Implementation of 3D Galerkin Boundary Integral Equations
Nintcheu Fata, Sylvain; Gray, Leonard J
2010-01-01
In this article, a reverse contribution technique is proposed to accelerate the construction of the dense influence matrices associated with a Galerkin approximation of singular and hypersingular boundary integral equations of mixed-type in potential theory. In addition, a general-purpose sparse preconditioner for boundary element methods has also been developed to successfully deal with ill-conditioned linear systems arising from the discretization of mixed boundary-value problems on non-smooth surfaces. The proposed preconditioner, which originates from the precorrected-FFT method, is sparse, easy to generate and apply in a Krylov subspace iterative solution of discretized boundary integral equations. Moreover, an approximate inverse of the preconditioner is implicitly built by employing an incomplete LU factorization. Numerical experiments involving mixed boundary-value problems for the Laplace equation are included to illustrate the performance and validity of the proposed techniques.
Concentration of Laplace Eigenfunctions and Stabilization of Weakly Damped Wave Equation
NASA Astrophysics Data System (ADS)
Burq, N.; Zuily, C.
2016-08-01
In this article, we prove some universal bounds on the speed of concentration on small (frequency-dependent) neighbourhoods of sub-manifolds of L 2-norms of quasi modes for Laplace operators on compact manifolds. We deduce new results on the rate of decay of weakly damped wave equations.
Derivation of new 3D discrete ordinate equations
Ahrens, C. D.
2012-07-01
The Sn equations have been the workhorse of deterministic radiation transport calculations for many years. Here we derive two new angular discretizations of the 3D transport equation. The first set of equations, derived using Lagrange interpolation and collocation, retains the classical Sn structure, with the main difference being how the scattering source is calculated. Because of the formal similarity with the classical S n equations, it should be possible to modify existing computer codes to take advantage of the new formulation. In addition, the new S n-like equations correctly capture delta function scattering. The second set of equations, derived using a Galerkin technique, does not retain the classical Sn structure because the streaming term is not diagonal. However, these equations can be cast into a form similar to existing methods developed to reduce ray effects. Numerical investigation of both sets of equations is under way. (authors)
Numerical dispersion analysis for three-dimensional Laplace-Fourier-domain scalar wave equation
NASA Astrophysics Data System (ADS)
Chen, Jing-Bo
2016-06-01
Based on the phase velocity and attenuation propagation velocity, a method for performing numerical dispersion analysis of three-dimensional Laplace-Fourier-domain scalar wave equation is presented. This method is applied to a 27-point average-derivative optimal scheme and a 27-point finite-element scheme. Within the relative error of 1%, the 27-point average-derivative optimal scheme requires seven grid points per wavelength and pseudo-wavelength while the 27-point finite-element scheme requires 23 grid points per wavelength and pseudo-wavelength for equal and unequal directional sampling intervals. Numerical examples show that the 27-point Laplace-Fourier-domain average-derivative optimal scheme is more accurate than the 27-point Laplace-Fourier-domain finite-element scheme for the same computational cost. By using larger directional sampling intervals while maintaining accuracy, the 27-point Laplace-Fourier-domain average-derivative optimal scheme can greatly reduce the computational cost of three-dimensional Laplace-Fourier-domain modelling.
Nonlinear Laplace equation, de Sitter vacua, and information geometry
Loran, Farhang
2005-06-15
Three exact solutions say {phi}{sub 0} of massless scalar theories on Euclidean space, i.e. D=6 {phi}{sup 3}, D=4 {phi}{sup 4} and D=3 {phi}{sup 6} models are obtained which share similar properties. The information geometry of their moduli spaces coincide with the Euclidean AdS{sub 7}, AdS{sub 5} and AdS{sub 4} respectively on which {phi}{sub 0} can be described as a stable tachyon. In D=4 we recognize that the SU(2) instanton density is proportional to {phi}{sub 0}{sup 4}. The original action S[{phi}] written in terms of new scalars {phi}-tilde={phi}-{phi}{sub 0} is shown to be equivalent to an interacting scalar theory on D-dimensional de Sitter background.
A note on singularities of the 3-D Euler equation
NASA Technical Reports Server (NTRS)
Tanveer, S.
1994-01-01
In this paper, we consider analytic initial conditions with finite energy, whose complex spatial continuation is a superposition of a smooth background flow and a singular field. Through explicit calculation in the complex plane, we show that under some assumptions, the solution to the 3-D Euler equation ceases to be analytic in the real domain in finite time.
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)
Shigeta, Takemi; Young, D. L.; Liu, Chein-Shan
2012-08-01
The mixed boundary value problem of the Laplace equation is considered. The method of fundamental solutions (MFS) approximates the exact solution to the Laplace equation by a linear combination of independent fundamental solutions with different source points. The accuracy of the numerical solution depends on the distribution of source points. In this paper, a weighted greedy QR decomposition (GQRD) is proposed to choose significant source points by introducing a weighting parameter. An index called an average degree of approximation is defined to show the efficiency of the proposed method. From numerical experiments, it is concluded that the numerical solution tends to be more accurate when the average degree of approximation is larger, and that the proposed method can yield more accurate solutions with a less number of source points than the conventional GQRD.
A parallel algorithm for solving the 3d Schroedinger equation
Strickland, Michael; Yager-Elorriaga, David
2010-08-20
We describe a parallel algorithm for solving the time-independent 3d Schroedinger equation using the finite difference time domain (FDTD) method. We introduce an optimized parallelization scheme that reduces communication overhead between computational nodes. We demonstrate that the compute time, t, scales inversely with the number of computational nodes as t {proportional_to} (N{sub nodes}){sup -0.95} {sup {+-} 0.04}. This makes it possible to solve the 3d Schroedinger equation on extremely large spatial lattices using a small computing cluster. In addition, we present a new method for precisely determining the energy eigenvalues and wavefunctions of quantum states based on a symmetry constraint on the FDTD initial condition. Finally, we discuss the usage of multi-resolution techniques in order to speed up convergence on extremely large lattices.
BV supersolutions to equations of 1-Laplace and minimal surface type
NASA Astrophysics Data System (ADS)
Scheven, Christoph; Schmidt, Thomas
2016-08-01
We propose notions of BV supersolutions to (the Dirichlet problem for) the 1-Laplace equation, the minimal surface equation, and equations of similar type. We then establish some related compactness and consistency results. Our main technical tool is a generalized product of L∞ divergence-measure fields and gradient measures of BV functions. This product crucially depends on the choice of a representative of the BV function, and the proofs of its basic properties involve results on one-sided approximation and fine (semi)continuity in the BV context.
Numerical simulation of vortex breakdown via 3-D Euler equations
NASA Astrophysics Data System (ADS)
Le, T. H.; Mege, P.; Morchoisne, Y.
1990-06-01
The long term goal is the modeling of vortex breakdown that occurs in some aerodynamic configurations at high angle of attack, (i.e., fighters with highly swept delta wings or missiles). A numerical simulation was made based on solving the 3-D Euler equations for an usteady incompressible flow. Preliminary results were obtained using a pressure-velocity formulation with periodic boundary conditions, the Euler equations being discretized by 2nd order finite difference schemes. The continuation to this work by implementing more realistic boundary conditions and 4th order finite difference discretization schemes are presented.
On the applicability of Young-Laplace equation for nanoscale liquid drops
NASA Astrophysics Data System (ADS)
Yan, Hong; Wei, Jiuan; Cui, Shuwen; Xu, Shenghua; Sun, Zhiwei; Zhu, Ruzeng
2016-03-01
Debates continue on the applicability of the Young-Laplace equation for droplets, vapor bubbles and gas bubbles in nanoscale. It is more meaningful to find the error range of the Young-Laplace equation in nanoscale instead of making the judgement of its applicability. To do this, for seven liquid argon drops (containing 800, 1000, 1200, 1400, 1600, 1800, or 2000 particles, respectively) at T = 78 K we determined the radius of surface of tension R s and the corresponding surface tension γ s by molecular dynamics simulation based on the expressions of R s and γ s in terms of the pressure distribution for droplets. Compared with the two-phase pressure difference directly obtained by MD simulation, the results show that the absolute values of relative error of two-phase pressure difference given by the Young-Laplace equation are between 0.0008 and 0.027, and the surface tension of the argon droplet increases with increasing radius of surface of tension, which supports that the Tolman length of Lennard-Jones droplets is positive and that Lennard-Jones vapor bubbles is negative. Besides, the logic error in the deduction of the expressions of the radius and the surface tension of surface of tension, and in terms of the pressure distribution for liquid drops in a certain literature is corrected.
Two-equation turbulence modeling for 3-D hypersonic flows
NASA Technical Reports Server (NTRS)
Bardina, J. E.; Coakley, T. J.; Marvin, J. G.
1992-01-01
An investigation to verify, incorporate and develop two-equation turbulence models for three-dimensional high speed flows is presented. The current design effort of hypersonic vehicles has led to an intensive study of turbulence models for compressible hypersonic flows. This research complements an extensive review of experimental data and the current development of 2D turbulence models. The review of experimental data on 2D and 3D flows includes complex hypersonic flows with pressure profiles, skin friction, wall heat transfer, and turbulence statistics data. In a parallel effort, turbulence models for high speed flows have been tested against flat plate boundary layers, and are being tested against the 2D database. In the present paper, we present the results of 3D Navier-Stokes numerical simulations with an improved k-omega two-equation turbulence model against experimental data and empirical correlations of an adiabatic flat plate boundary layer, a cold wall flat plate boundary layer, and a 3D database flow, the interaction of an oblique shock wave and a thick turbulent boundary layer with a free stream Mach number = 8.18 and Reynolds number = 5 x 10 to the 6th.
NASA Astrophysics Data System (ADS)
Liang, Yingjie; Chen, Wen; Magin, Richard L.
2016-07-01
Analytical solutions to the fractional diffusion equation are often obtained by using Laplace and Fourier transforms, which conveniently encode the order of the time and the space derivatives (α and β) as non-integer powers of the conjugate transform variables (s, and k) for the spectral and the spatial frequencies, respectively. This study presents a new solution to the fractional diffusion equation obtained using the Laplace transform and expressed as a Fox's H-function. This result clearly illustrates the kinetics of the underlying stochastic process in terms of the Laplace spectral frequency and entropy. The spectral entropy is numerically calculated by using the direct integration method and the adaptive Gauss-Kronrod quadrature algorithm. Here, the properties of spectral entropy are investigated for the cases of sub-diffusion and super-diffusion. We find that the overall spectral entropy decreases with the increasing α and β, and that the normal or Gaussian case with α = 1 and β = 2, has the lowest spectral entropy (i.e., less information is needed to describe the state of a Gaussian process). In addition, as the neighborhood over which the entropy is calculated increases, the spectral entropy decreases, which implies a spatial averaging or coarse graining of the material properties. Consequently, the spectral entropy is shown to provide a new way to characterize the temporal correlation of anomalous diffusion. Future studies should be designed to examine changes of spectral entropy in physical, chemical and biological systems undergoing phase changes, chemical reactions and tissue regeneration.
Potentially singular solutions of the 3D axisymmetric Euler equations
Luo, Guo; Hou, Thomas Y.
2014-01-01
The question of finite-time blowup of the 3D incompressible Euler equations is numerically investigated in a periodic cylinder with solid boundaries. Using rotational symmetry, the equations are discretized in the (2D) meridian plane on an adaptive (moving) mesh and is integrated in time with adaptively chosen time steps. The vorticity is observed to develop a ring-singularity on the solid boundary with a growth proportional to ∼(ts − t)−2.46, where ts ∼ 0.0035056 is the estimated singularity time. A local analysis also suggests the existence of a self-similar blowup. The simulations stop at τ2 = 0.003505 at which time the vorticity amplifies by more than (3 × 108)-fold and the maximum mesh resolution exceeds (3 × 1012)2. The vorticity vector is observed to maintain four significant digits throughout the computations. PMID:25157172
Modeling tree crown dynamics with 3D partial differential equations.
Beyer, Robert; Letort, Véronique; Cournède, Paul-Henry
2014-01-01
We characterize a tree's spatial foliage distribution by the local leaf area density. Considering this spatially continuous variable allows to describe the spatiotemporal evolution of the tree crown by means of 3D partial differential equations. These offer a framework to rigorously take locally and adaptively acting effects into account, notably the growth toward light. Biomass production through photosynthesis and the allocation to foliage and wood are readily included in this model framework. The system of equations stands out due to its inherent dynamic property of self-organization and spontaneous adaptation, generating complex behavior from even only a few parameters. The density-based approach yields spatially structured tree crowns without relying on detailed geometry. We present the methodological fundamentals of such a modeling approach and discuss further prospects and applications. PMID:25101095
On the normal modes of Laplace's tidal equations for zonal wavenumber zero
NASA Technical Reports Server (NTRS)
Tanaka, H. L.; Kasahara, Akira
1992-01-01
The characteristic differences between two different rotational modes of Laplace's tidal equations for wavenumber m = 0, called the K- and the S-modes, are compared in their energy ratio and structures. It is shown that the K-mode representation captures most of the observed zonal energy with a few terms, whereas the S-mode representation requires many terms. For small vertical scale components, the K-mode series converges faster than the S-mode series. Attention is also given to the differences between the energy spectra projected upon the K- and S-modes and the merits of each set as expansion functions for the zonal atmospheric motions.
Padé approximation of Laplace transforms of some special functions in terms of Painlevé equations
NASA Astrophysics Data System (ADS)
Nakamura, Yoshimasa; Ohira, Norihiro
2004-10-01
Only in a limiting case has a closed form of the Laplace transforms of special functions been known (Abramowitz M and Stegun I A 1970 Handbook of Mathematical Functions (New York: Dover)). It is shown that Padé approximations of Laplace transforms of the Airy, Hermite-Weber, Bessel functions are computed by using Bäcklund transformations of the Painlevé equations in a purely algebraic manner.
Equations on knot polynomials and 3d/5d duality
Mironov, A.; Morozov, A.
2012-09-24
We briefly review the current situation with various relations between knot/braid polynomials (Chern-Simons correlation functions), ordinary and extended, considered as functions of the representation and of the knot topology. These include linear skein relations, quadratic Plucker relations, as well as 'differential' and (quantum) A-polynomial structures. We pay a special attention to identity between the A-polynomial equations for knots and Baxter equations for quantum relativistic integrable systems, related through Seiberg-Witten theory to 5d super-Yang-Mills models and through the AGT relation to the q-Virasoro algebra. This identity is an important ingredient of emerging a 3d- 5d generalization of the AGT relation. The shape of the Baxter equation (including the values of coefficients) depend on the choice of the knot/braid. Thus, like the case of KP integrability, where (some, so far torus) knots parameterize particular points of the Universal Grassmannian, in this relation they parameterize particular points in the moduli space of many-body integrable systems of relativistic type.
Numerical solution for Laplace equation with mixed boundary condition for ship problem in the sea
NASA Astrophysics Data System (ADS)
Silalahi, Fitriani Tupa R.; Budhi, Wono Setya; Adytia, Didit; van Groesen, E.
2015-09-01
One interesting phenomena is investigating the movement of ships at the sea. To start with the investigation in modelling of this problem, we will assume that the ship is only a one-dimensional object that is floating on the sea surface. Similarly, we assume that the water flow is uniform in parallel directions to the ship. Therefore, we simply use the two-dimensional Laplace equation in this problem. In the section that describes the surface of sea, Neumann boundary condition is imposed in part related to the ship and the Dirichlet boundary condition for others. Then on the other three boundaries, we imposed the Neumann boundary condition by assuming that the water does not flow on the bottom, and both end. The model is solved by numerical solution using the finite element method. Velocity potential solution on the whole domain is demonstrated as a result of the implementation of the finite element method. In this paper, we initiate an investigation with assuming that the ship is on the water so that the domain of the Laplace equation is rectangular. Then we assume the drift ship. Furthermore, we also study the dependence of width and depth of the domain to the velocity potential.
A parallel multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation
Riyanti, C.D. . E-mail: C.D.Riyanti@tudelft.nl; Kononov, A.; Erlangga, Y.A.; Vuik, C.; Oosterlee, C.W.; Plessix, R.-E.; Mulder, W.A.
2007-05-20
We investigate the parallel performance of an iterative solver for 3D heterogeneous Helmholtz problems related to applications in seismic wave propagation. For large 3D problems, the computation is no longer feasible on a single processor, and the memory requirements increase rapidly. Therefore, parallelization of the solver is needed. We employ a complex shifted-Laplace preconditioner combined with the Bi-CGSTAB iterative method and use a multigrid method to approximate the inverse of the resulting preconditioning operator. A 3D multigrid method with 2D semi-coarsening is employed. We show numerical results for large problems arising in geophysical applications.
Non-rigid registration of breast surfaces using the laplace and diffusion equations
2010-01-01
A semi-automated, non-rigid breast surface registration method is presented that involves solving the Laplace or diffusion equations over undeformed and deformed breast surfaces. The resulting potential energy fields and isocontours are used to establish surface correspondence. This novel surface-based method, which does not require intensity images, anatomical landmarks, or fiducials, is compared to a gold standard of thin-plate spline (TPS) interpolation. Realistic finite element simulations of breast compression and further testing against a tissue-mimicking phantom demonstrate that this method is capable of registering surfaces experiencing 6 - 36 mm compression to within a mean error of 0.5 - 5.7 mm. PMID:20149261
Lu, Dingjie; Xie, Yi Min; Huang, Xiaodong; Zhou, Shiwei; Li, Qing
2015-11-28
Analytical studies on the size effects of a simply-shaped beam fixed at both ends have successfully explained the sudden changes of effective Young's modulus as its diameter decreases below 100 nm. Yet they are invalid for complex nanostructures ubiquitously existing in nature. In accordance with a generalized Young-Laplace equation, one of the representative size effects is transferred to non-uniformly distributed pressure against an external surface due to the imbalance of inward and outward loads. Because the magnitude of pressure depends on the principal curvatures, iterative steps have to be adopted to gradually stabilize the structure in finite element analysis. Computational results are in good agreement with both experiment data and theoretical prediction. Furthermore, the investigation on strengthened and softened Young's modulus for two complex nanostructures demonstrates that the proposed computational method provides a general and effective approach to analyze the size effects for nanostructures in arbitrary shape.
Laplace's equation and the Dirichlet-Neumann map in multiply connected domains
NASA Technical Reports Server (NTRS)
Greenbaum, A.; Greengard, L.; Mcfadden, G. B.
1993-01-01
A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.
Waveform inversion in the Laplace and Laplace-Fourier domains
NASA Astrophysics Data System (ADS)
Shin, Changsoo; Ha, Wansoo
2010-05-01
adjusting the real and imaginary part of these complex frequencies, we can numerically calculate Laplace-Fourier transformed wavefields for several damped wavefields in the time domain. The complex-valued wavefield allows us to use both amplitude and phase information, while in the Laplace domain, only the real-valued wavefield is used for the inversion. Through simultaneous inversion of the amplitude and phase of complex-valued logarithmic wavefields, medium- and short-wavelength velocity structures can be recovered in addition to long-wavelength structures, and the penetration depth of the inversion can be enhanced. Since the Laplace-domain techniques are based on the damped wavefields, Laplace-domain inversion is very sensitive to noise appearing before the first arrival. Therefore, careful muting of this noise is needed for successful inversion to real data. To make the Laplace transform stable, the maximum recording time should be long enough so that the amplitude of damped signals after the maximum recording time can be ignored within the noise level. In particular, Laplace domain inversion also requires larger offset data than the frequency-domain inversion because the logarithmic residual tends to form a long-wavelength shape. Inversion results of Laplace and Laplace-Fourier domain inversion can be used as a starting velocity model for subsequent frequency-domain full waveform inversion. This two-step strategy may produce correct high-resolution velocity models from several sources of synthetic and field data. Laplace-domain inversions with a coarse grid make 3-D seismic inversion feasible. Elastic and acoustic-elastic coupled inversions in the Laplace and Laplace-Fourier domain have the potential to delineate elastic parameters from land and marine data. Data with a large offset and long recording time can produce good inversion results in either the Laplace or Laplace-Fourier domain. The selection of the optimum combination of Laplace-Fourier frequencies, as well
van Pelt, J
1992-01-01
Transient potentials in dendritic trees can be calculated by approximating the dendrite by a set of connected cylinders. The profiles for the currents and potentials in the whole system can then be obtained by imposing the proper boundary conditions and calculating these profiles along each individual cylinder. An elegant implementation of this method has been described by Holmes (1986), and is based on the Laplace transform of the cable equation. By calculating the currents and potentials only at the ends of the cylinders, the whole system of connected cylinders can be described by a set of n equations, where n denotes the number of internal and external nodes (points of connection and endpoints of the cylinders). The present study shows that the set of equations can be formulated by a simple vector equation which is essentially a generalization of Ohm's law for the whole system. The current and potential n-vectors are coupled by a n x n conductance matrix whose structure immediately reflects the connectivity pattern of the connected cylinders. The vector equation accounts for conductances, associated with driving potentials, which may be local or distributed over the membrane. It is shown that the vector equation can easily be adapted for the calculation of transients over a period in which stepwise changes in system parameters have occurred. In this adaptation it is assumed that the initial conditions for the potential profiles at the start of a new period after a stepwise change can be approximated by steady-state solutions.(ABSTRACT TRUNCATED AT 250 WORDS) PMID:1486128
Massopust, P.R.
1997-08-01
All solutions of an in its angular coordinates continuously perturbed Laplace-Beltrami equation in the open unit ball IB{sup n+2} {contained_in} IR{sup n+2}, n {ge} 1, are characterized. Moreover, it is shown that such pertubations yield distributional boundary values which are different from, but algebraically and topologically equivalent to, the hyperfunctions of Lions & Magenes. This is different from the case of radially perturbed Laplace-Beltrami operators (cf. [7]) where one has stability of distributional boundary values under such perturbations.
A quasi-spectral method for Cauchy problem of 2/D Laplace equation on an annulus
NASA Astrophysics Data System (ADS)
Saito, Katsuyoshi; Nakada, Manabu; Iijima, Kentaro; Onishi, Kazuei
2005-01-01
Real numbers are usually represented in the computer as a finite number of digits hexa-decimal floating point numbers. Accordingly the numerical analysis is often suffered from rounding errors. The rounding errors particularly deteriorate the precision of numerical solution in inverse and ill-posed problems. We attempt to use a multi-precision arithmetic for reducing the rounding error evil. The use of the multi-precision arithmetic system is by the courtesy of Dr Fujiwara of Kyoto University. In this paper we try to show effectiveness of the multi-precision arithmetic by taking two typical examples; the Cauchy problem of the Laplace equation in two dimensions and the shape identification problem by inverse scattering in three dimensions. It is concluded from a few numerical examples that the multi-precision arithmetic works well on the resolution of those numerical solutions, as it is combined with the high order finite difference method for the Cauchy problem and with the eigenfunction expansion method for the inverse scattering problem.
Use of the augmented Young-Laplace equation to model equilibrium and evaporating extended menisci
DasGupta, S.; Schonberg, J.A.; Kim, I.Y.; Wayner, P.C.Jr. )
1993-05-01
The generic importance of fluid flow and change-of-phase heat transfer in the contact line region of an extended meniscus has led to theoretical and experimental research on the details of these transport processes. Numerical solutions of equilibrium and nonequilibrium models based on the augmented Young-Laplace equation were successfully used to evaluate experimental data for an extended meniscus. The data for the equilibrium and nonequilibrium meniscus profiles were obtained optically using ellipsometry and image processing interferometry. A Taylor series expansion of the fourth-order nonlinear transport model was used to obtain the extremely sensitive initial conditions at the interline. The solid-liquid-vapor Hamaker constants for the systems were obtained from the experimental data. The consistency of the data was demonstrated by using the combining rules to calculate the unknown value of the Hamaker constant for the experimental substrate. The sensitivity of the meniscus profile to small changes in the environment was demonstrated. Both temperature and intermolecular forces need to be included in modeling transport processes in the contact line region because the chemical potential is a function of both temperature and pressure.
Exponential Mixing of the 3D Stochastic Navier-Stokes Equations Driven by Mildly Degenerate Noises
Albeverio, Sergio; Debussche, Arnaud; Xu Lihu
2012-10-15
We prove the strong Feller property and exponential mixing for 3D stochastic Navier-Stokes equation driven by mildly degenerate noises (i.e. all but finitely many Fourier modes being forced) via a Kolmogorov equation approach.
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.
Chatterjee, Jaideep
2007-10-01
This paper shows how 2 coupled Young-Laplace equations can be solved to predict the shapes of two coupled menisci formed in a capillary system. Experiments are performed, which demonstrate that the equilibrium volume of liquid retained in a vertical capillary, can be variable, even when all the properties of the system are invariant. This variability in liquid retention also leads to different equilibrium shapes of the top and bottom menisci. A coupled form of the Young-Laplace equation is solved to predict the two coupled menisci shapes. The curvature of the top meniscus is fitted to the experimentally recorded meniscus shape. The coupled Young-Laplace equation solution is used to predict the shape of the bottom meniscus. The shape of the bottom meniscus thus obtained, is shown to match the experimentally recorded bottom meniscus shape reasonably well. This observed coupling of the menisci has a significant impact on some porosimetric techniques which are based on liquid extrusion and explains why the volume of liquid that can be retained in a capillary can vary, under invariant conditions. Retention of liquids in capillaries is of interest in several applications like fabric wash. PMID:17568603
2D/1D approximations to the 3D neutron transport equation. I: Theory
Kelley, B. W.; Larsen, E. W.
2013-07-01
A new class of '2D/1D' approximations is proposed for the 3D linear Boltzmann equation. These approximate equations preserve the exact transport physics in the radial directions x and y and diffusion physics in the axial direction z. Thus, the 2D/1D equations are more accurate approximations of the 3D Boltzmann equation than the conventional 3D diffusion equation. The 2D/1D equations can be systematically discretized, to yield accurate simulation methods for 3D reactor core problems. The resulting solutions will be more accurate than 3D diffusion solutions, and less expensive to generate than standard 3D transport solutions. In this paper, we (i) show that the simplest 2D/1D equation has certain desirable properties, (ii) systematically discretize this equation, and (iii) derive a stable iteration scheme for solving the discrete system of equations. In a companion paper [1], we give numerical results that confirm the theoretical predictions of accuracy and iterative stability. (authors)
2D/1D approximations to the 3D neutron transport equation. II: Numerical comparisons
Kelley, B. W.; Collins, B.; Larsen, E. W.
2013-07-01
In a companion paper [1], (i) several new '2D/1D equations' are introduced as accurate approximations to the 3D Boltzmann transport equation, (ii) the simplest of these approximate equations is systematically discretized, and (iii) a theoretically stable iteration scheme is developed to solve the discrete equations. In this paper, numerical results are presented that confirm the theoretical predictions made in [1]. (authors)
On the Dynamic Programming Approach for the 3D Navier-Stokes Equations
Manca, Luigi
2008-06-15
The dynamic programming approach for the control of a 3D flow governed by the stochastic Navier-Stokes equations for incompressible fluid in a bounded domain is studied. By a compactness argument, existence of solutions for the associated Hamilton-Jacobi-Bellman equation is proved. Finally, existence of an optimal control through the feedback formula and of an optimal state is discussed.
A Laplacian Equation Method for Numerical Generation of Boundary-Fitted 3D Orthogonal Grids
NASA Astrophysics Data System (ADS)
Theodoropoulos, T.; Bergeles, G. C.
1989-06-01
A sethod for generating boundary fitted orthogonal curvilinear grids in 3-dimensional space is described. The mapping between the curvilinear coordinates and the Cartesian coordinates is provided by a set of Laplace equations which, expressed in curvilinear coordinates, involve the components of the metric tensor and are therefore non-linear and coupled. An iterative algorithm is described, which achieves a numerical solution. Grids appropriate for the calculation of flow fields over complex topography or in complex flow passages as those found in turbomachinery, and for other engineering applications can be constructed using the proposed method. Various examples are presented and plotted in perspective, and data for the assessment of the properties of the resulting meshes is provided.
Implementation of Advanced Two Equation Turbulence Models in the USM3D Unstructured Flow Solver
NASA Technical Reports Server (NTRS)
Wang, Qun-Zhen; Massey, Steven J.; Abdol-Hamid, Khaled S.
2000-01-01
USM3D is a widely-used unstructured flow solver for simulating inviscid and viscous flows over complex geometries. The current version (version 5.0) of USM3D, however, does not have advanced turbulence models to accurately simulate complicated flow. We have implemented two modified versions of the original Jones and Launder k-epsilon "two-equation" turbulence model and the Girimaji algebraic Reynolds stress model in USM3D. Tests have been conducted for three flat plate boundary layer cases, a RAE2822 airfoil and an ONERA M6 wing. The results are compared with those from direct numerical simulation, empirical formulae, theoretical results, and the existing Spalart-Allmaras one-equation model.
Global regular solutions for the 3D Kawahara equation posed on unbounded domains
NASA Astrophysics Data System (ADS)
Larkin, Nikolai A.; Simões, Márcio Hiran
2016-08-01
An initial boundary value problem for the 3D Kawahara equation posed on a channel-type domain was considered. The existence and uniqueness results for global regular solutions as well as exponential decay of small solutions in the H 2-norm were established.
Global regular solutions for the 3D Zakharov-Kuznetsov equation posed on unbounded domains
NASA Astrophysics Data System (ADS)
Larkin, N. A.
2015-09-01
An initial-boundary value problem for the 3D Zakharov-Kuznetsov equation posed on unbounded domains is considered. Existence and uniqueness of a global regular solution as well as exponential decay of the H2-norm for small initial data are proven.
Quasi-regular solutions to a class of 3D degenerating hyperbolic equations
NASA Astrophysics Data System (ADS)
Hristov, T. D.; Popivanov, N. I.; Schneider, M.
2012-11-01
In the fifties M. Protter stated new three-dimensional (3D) boundary value problems (BVP) for mixed type equations of first kind. For hyperbolic-elliptic equations they are multidimensional analogue of the classical two-dimensional (2D) Morawetz-Guderley transonic problem. Up to now, in this case, not a single example of nontrivial solution to the new problem, neither a general existence result is known. The difficulties appear even for BVP in the hyperbolic part of the domain, that were formulated by Protter for weakly hyperbolic equations. In that case the Protter problems are 3D analogues of the plane Darboux or Cauchy-Goursat problems. It is interesting that in contrast to the planar problems the new 3D problems are strongly ill-posed. Some of the Protter problems for degenerating hyperbolic equation without lower order terms or even for the usual wave equation have infinite-dimensional kernels. Therefore there are infinitely many orthogonality conditions for classical solvability of their adjiont problems. So it is interesting to obtain results for uniqueness of solutions adding first order terms in the equation. In the present paper we do this and find conditions for coefficients under which we prove uniqueness of quasi-regular solutions to the Protter problems.
NASA Astrophysics Data System (ADS)
Faria, Paula
2010-09-01
For the past few years, the potential of transcranial direct current stimulation (tDCS) for the treatment of several pathologies has been investigated. Knowledge of the current density distribution is an important factor in optimizing such applications of tDCS. For this goal, we used the finite element method to solve the Laplace equation in a spherical head model in order to investigate the three dimensional distribution of the current density and the variation of its intensity with depth using different electrodes montages: the traditional one with two sponge electrodes and new electrode montages: with sponge and EEG electrodes and with EEG electrodes varying the numbers of electrodes. The simulation results confirm the effectiveness of the mixed system which may allow the use of tDCS and EEG recording concomitantly and may help to optimize this neuronal stimulation technique. The numerical results were used in a promising application of tDCS in epilepsy.
A research of 3D gravity inversion based on the recovery of sparse underdetermined linear equations
NASA Astrophysics Data System (ADS)
Zhaohai, M.
2014-12-01
Because of the properties of gravity data, it is made difficult to solve the problem of multiple solutions. There are two main types of 3D gravity inversion methods：One of two methods is based on the improvement of the instability of the sensitive matrix, solving the problem of multiple solutions and instability in 3D gravity inversion. Another is to join weight function into the 3D gravity inversion iteration. Through constant iteration, it can renewal density values and weight function to achieve the purpose to solve the multiple solutions and instability of the 3D gravity data inversion. Thanks to the sparse nature of the solutions of 3D gravity data inversions, we can transform it into a sparse equation. Then, through solving the sparse equations, we can get perfect 3D gravity inversion results. The main principle is based on zero norm of sparse matrix solution of the equation. Zero norm is mainly to solve the nonzero solution of the sparse matrix. However, the method of this article adopted is same as the principle of zero norm. But the method is the opposite of zero norm to obtain zero value solution. Through the form of a Gaussian fitting solution of the zero norm, we can find the solution by using regularization principle. Moreover, this method has been proved that it had a certain resistance to random noise in the mathematics, and it was more suitable than zero norm for the solution of the geophysical data. 3D gravity which is adopted in this article can well identify abnormal body density distribution characteristics, and it can also recognize the space position of abnormal distribution very well. We can take advantage of the density of the upper and lower limit penalty function to make each rectangular residual density within a reasonable range. Finally, this 3D gravity inversion is applied to a variety of combination model test, such as a single straight three-dimensional model, the adjacent straight three-dimensional model and Y three
Laplace and Z Transform Solutions of Differential and Difference Equations With the HP-41C.
ERIC Educational Resources Information Center
Harden, Richard C.; Simons, Fred O., Jr.
1983-01-01
A previously developed program for the HP-41C programmable calculator is extended to handle models of differential and difference equations with multiple eigenvalues. How to obtain difference equation solutions via the Z transform is described. (MNS)
The small data solutions of general 3-D quasilinear wave equations. II
NASA Astrophysics Data System (ADS)
Ding, Bingbing; Witt, Ingo; Yin, Huicheng
2016-07-01
This paper is a continuation of the work in [8], where the authors established the global existence of smooth small data solutions to the general 3-D quasilinear wave equation ∑ i , j = 0 3 gij (u , ∂ u) ∂ij2 u = 0 when the weak null condition holds. In the present paper, we show that the smooth small data solutions of equation ∑ i , j = 0 3 gij (u , ∂ u) ∂ij2 u = 0 will blow up in finite time when the weak null condition does not hold and a generic nondegenerate condition on the initial data is satisfied, moreover, a precise blowup time is completely determined. Therefore, collecting the main results in this paper and [8], we have given a basically complete study on the blowup or global existence of small data solutions to the 3-D quasilinear wave equation ∑ i , j = 0 3 gij (u , ∂ u) ∂ij2 u = 0.
Upper Semicontinuity of Pullback Attractors for the 3D Nonautonomous Benjamin-Bona-Mahony Equations
Yang, Xinguang; Wang, Xiaosong; Zhang, Lingrui
2014-01-01
We will study the upper semicontinuity of pullback attractors for the 3D nonautonomouss Benjamin-Bona-Mahony equations with external force perturbation terms. Under some regular assumptions, we can prove the pullback attractors 𝒜 ε(t) of equation ut-Δut-νΔu+∇·F→(u)=ɛg(x,t), x ∈ Ω, converge to the global attractor 𝒜 of the above-mentioned equation with ε = 0 for any t ∈ ℝ. PMID:24790585
Recasting the 3D Wigner-Liouville equation with spectral components of the force
NASA Astrophysics Data System (ADS)
van de Put, Maarten; Sorée, Bart; Magnus, Wim
The phasespace approach to many-body quantum mechanics, by means of the Wigner-function is interesting through its connection to classical mechanics. Time-evolution of any statistical distribution of states under influence of a (time-dependent) Hamiltonian is obtained through use of the Wigner-Liouville equation. The standard form of this equation contains two 3D integrals, over the entire phase space. As a result, this form emphasizes the non-locality of the interaction of the potential, but lacks simplicity and ease of understanding. Furthermore, the integrals make numerical solution of the Wigner-Liouville equation challenging. We present an alternative form to the Wigner-Liouville equation based on the force rather than the potential, in alignment with the classical Boltzmann equation. Decomposition of the force in its spectral components yields a simpler form of the Wigner-Liouville equation. This new form has only one 3D integral over the spectral force components, and is local in position, simplifying both interpretation and numerical implementation. Because of its use of the force, it straightforwardly reduces to the Boltzmann equation under classical conditions.
Equation-of-State Test Suite for the DYNA3D Code
Benjamin, Russell D.
2015-11-05
This document describes the creation and implementation of a test suite for the Equationof- State models in the DYNA3D code. A customized input deck has been created for each model, as well as a script that extracts the relevant data from the high-speed edit file created by DYNA3D. Each equation-of-state model is broken apart and individual elements of the model are tested, as well as testing the entire model. The input deck for each model is described and the results of the tests are discussed. The intent of this work is to add this test suite to the validation suite presently used for DYNA3D.
Active exterior cloaking for the 2D Laplace and Helmholtz equations.
Vasquez, Fernando Guevara; Milton, Graeme W; Onofrei, Daniel
2009-08-14
A new cloaking method is presented for 2D quasistatics and the 2D Helmholtz equation that we speculate extends to other linear wave equations. For 2D quasistatics it is proven how a single active exterior cloaking device can be used to shield an object from surrounding fields, yet produce very small scattered fields. The problem is reduced to finding a polynomial which is close to 1 in a disk and close to 0 in another disk, and such a polynomial is constructed. For the 2D Helmholtz equation it is numerically shown that three exterior cloaking devices placed around the object suffice to hide it. PMID:19792644
Meng, Da; Zheng, Bin; Lin, Guang; Sushko, Maria L.
2014-08-29
We have developed efficient numerical algorithms for the solution of 3D steady-state Poisson-Nernst-Planck equations (PNP) with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are discretized by finite difference scheme and solved iteratively by Gummel method with relaxation. The Nernst-Planck equations are transformed into Laplace equations through the Slotboom transformation. Algebraic multigrid method is then applied to efficiently solve the Poisson equation and the transformed Nernst-Planck equations. A novel strategy for calculating excess chemical potentials through fast Fourier transforms is proposed which reduces computational complexity from O(N2) to O(NlogN) where N is the number of grid points. Integrals involving Dirac delta function are evaluated directly by coordinate transformation which yields more accurate result compared to applying numerical quadrature to an approximated delta function. Numerical results for ion and electron transport in solid electrolyte for Li ion batteries are shown to be in good agreement with the experimental data and the results from previous studies.
Xiong, Z.; Tripp, A.C.
1994-12-31
This paper presents an integral equation algorithm for 3D EM modeling at high frequencies for applications in engineering an environmental studies. The integral equation method remains the same for low and high frequencies, but the dominant roles of the displacements currents complicate both numerical treatments and interpretations. With singularity extraction technique they successively extended the application of the Hankel filtering technique to the computation of Hankel integrals occurring in high frequency EM modeling. Time domain results are calculated from frequency domain results via Fourier transforms. While frequency domain data are not obvious for interpretations, time domain data show wave-like pictures that resemble seismograms. Both 1D and 3D numerical results show clearly the layer interfaces.
NuSol - Numerical solver for the 3D stationary nuclear Schrödinger equation
NASA Astrophysics Data System (ADS)
Graen, Timo; Grubmüller, Helmut
2016-01-01
The classification of short hydrogen bonds depends on several factors including the shape and energy spacing between the nuclear eigenstates of the hydrogen. Here, we describe the NuSol program in which three classes of algorithms were implemented to solve the 1D, 2D and 3D time independent nuclear Schrödinger equation. The Schrödinger equation was solved using the finite differences based Numerov's method which was extended to higher dimensions, the more accurate pseudo-spectral Chebyshev collocation method and the sinc discrete variable representation by Colbert and Miller. NuSol can be applied to solve the Schrödinger equation for arbitrary analytical or numerical potentials with focus on nuclei bound by the potential of their molecular environment. We validated the methods against literature values for the 2D Henon-Heiles potential, the 3D linearly coupled sextic oscillators and applied them to study hydrogen bonding in the malonaldehyde derivate 4-cyano-2,2,6,6-tetramethyl-3,5-heptanedione. With NuSol, the extent of nuclear delocalization in a given molecular potential can directly be calculated without relying on linear reaction coordinates in 3D molecular space.
Some Properties of the M3D-C1 Form of the 3D Magnetohydrodynamics Equations
J. Breslau, N. Ferraro, S. Jardin
2009-07-10
We introduce a set of scalar variables and projection operators for the vector momentum and magnetic field evolution equations that have several unique and desirable properties, making them a preferred system for solving the magnetohydrodynamics equations in a torus with a strong toroidal magnetic field. We derive a "weak form" of these equations that explicitly conserves energy and is suitable for a Galerkin finite element formulation provided the basis elements have C1 continuity. Systems of reduced equations are discussed, along with their energy conservation properties. An implicit time advance is presented that adds diagonally dominant self-adjoint energy terms to the mass matrix to obtain numerical stability.
Analytical and numerical aspects in solving the controlled 3D Gross-Pitaevskii equation
Fedele, R.; Jovanovic, D.; De Nicola, S.; Eliasson, B.; Shukla, P. K.
2009-11-10
The results of recently developed investigations, that have been carried out within the framework of the controlling potential method (CPM), are reviewed. This method allows one to decompose a three dimensional (3D) Gross-Pitaevskii equation (GPE) into the pair of coupled Schroedinger-type equations. Under suitable mathematical conditions, the solutions of the 3D controlled GPE can be constructed from the solutions of a 2D linear Schroedinger equation (the transverse component of the GPE) coupled with a 1D nonlinear Schroedinger equation (the longitudinal component of the GPE). Such decomposition allows one to cast the solutions in the form of the product of the solutions of the transverse and the longitudinal components of the GPE. The coupling between these two equations is the functional of both the transverse and the longitudinal profiles. It is shown that the CPM can be used to obtain a new class of three-dimensional solitary waves solutions of the GPE, which governs the dynamics of Bose-Einstein condensates. By imposing an external controlling potential, the desired time-dependent shape of the localized BECs is obtained. The stability of the exact solutions was checked with direct simulations of the time -dependent, three-dimensional GPE. Our simulations show that the localized condensates are stable with respect to perturbed initial conditions.
A least-squares finite element method for 3D incompressible Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Jiang, Bo-Nan; Lin, T. L.; Hou, Lin-Jun; Povinelli, Louis A.
1993-01-01
The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations, and results in symmetric, positive definite algebraic system. An additional compatibility equation, i.e., the divergence of vorticity vector should be zero, is included to make the first-order system elliptic. The Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. The flow in a half of 3D cubic cavity is calculated at Re = 100, 400, and 1,000 with 50 x 52 x 25 trilinear elements. The Taylor-Gortler-like vortices are observed at Re = 1,000.
A novel numerical flux for the 3D Euler equations with general equation of state
NASA Astrophysics Data System (ADS)
Toro, Eleuterio F.; Castro, Cristóbal E.; Lee, Bok Jik
2015-12-01
Here we extend the flux vector splitting approach recently proposed in E.F. Toro and M.E. Vázquez-Cendón (2012) [42]. The scheme was originally presented for the 1D Euler equations for ideal gases and its extension presented in this paper is threefold: (i) we solve the three-dimensional Euler equations on general meshes; (ii) we use a general equation of state; and (iii) we achieve high order of accuracy in both space and time through application of the semi-discrete ADER methodology on general meshes. The resulting methods are systematically assessed for accuracy, robustness and efficiency on a carefully selected suite of test problems. Formal high accuracy is assessed through convergence rates studies for schemes of up to 4th order of accuracy in both space and time on unstructured meshes.
NASA Astrophysics Data System (ADS)
DeJong, Andrew
Numerical models of fluid-structure interaction have grown in importance due to increasing interest in environmental energy harvesting, airfoil-gust interactions, and bio-inspired formation flying. Powered by increasingly powerful parallel computers, such models seek to explain the fundamental physics behind the complex, unsteady fluid-structure phenomena. To this end, a high-fidelity computational model based on the high-order spectral difference method on 3D unstructured, dynamic meshes has been developed. The spectral difference method constructs continuous solution fields within each element with a Riemann solver to compute the inviscid fluxes at the element interfaces and an averaging mechanism to compute the viscous fluxes. This method has shown promise in the past as a highly accurate, yet sufficiently fast method for solving unsteady viscous compressible flows. The solver is monolithically coupled to the equations of motion of an elastically mounted 3-degree of freedom rigid bluff body undergoing flow-induced lift, drag, and torque. The mesh is deformed using 4 methods: an analytic function, Laplace equation, biharmonic equation, and a bi-elliptic equation with variable diffusivity. This single system of equations -- fluid and structure -- is advanced through time using a 5-stage, 4th-order Runge-Kutta scheme. Message Passing Interface is used to run the coupled system in parallel on up to 240 processors. The solver is validated against previously published numerical and experimental data for an elastically mounted cylinder. The effect of adding an upstream body and inducing wake galloping is observed.
A 3D GCL compatible cell-centered Lagrangian scheme for solving gas dynamics equations
NASA Astrophysics Data System (ADS)
Georges, Gabriel; Breil, Jérôme; Maire, Pierre-Henri
2016-01-01
Solving the gas dynamics equations under the Lagrangian formalism enables to simulate complex flows with strong shock waves. This formulation is well suited to the simulation of multi-material compressible fluid flows such as those encountered in the domain of High Energy Density Physics (HEDP). These types of flows are characterized by complex 3D structures such as hydrodynamic instabilities (Richtmyer-Meshkov, Rayleigh-Taylor, etc.). Recently, the 3D extension of different Lagrangian schemes has been proposed and appears to be challenging. More precisely, the definition of the cell geometry in the 3D space through the treatment of its non-planar faces and the limiting of a reconstructed field in 3D in the case of a second-order extension are of great interest. This paper proposes two new methods to solve these problems. A systematic and symmetric geometrical decomposition of polyhedral cells is presented. This method enables to define a discrete divergence operator leading to the respect of the Geometric Conservation Law (GCL). Moreover, a multi-dimensional minmod limiter is proposed. This new limiter constructs, from nodal gradients, a cell gradient which enables to ensure the monotonicity of the numerical solution even in presence of strong discontinuity. These new ingredients are employed into a cell-centered Lagrangian scheme. Robustness and accuracy are assessed against various representative test cases.
A fast rebinning algorithm for 3D positron emission tomography using John's equation
NASA Astrophysics Data System (ADS)
Defrise, Michel; Liu, Xuan
1999-08-01
Volume imaging in positron emission tomography (PET) requires the inversion of the three-dimensional (3D) x-ray transform. The usual solution to this problem is based on 3D filtered-backprojection (FBP), but is slow. Alternative methods have been proposed which factor the 3D data into independent 2D data sets corresponding to the 2D Radon transforms of a stack of parallel slices. Each slice is then reconstructed using 2D FBP. These so-called rebinning methods are numerically efficient but are approximate. In this paper a new exact rebinning method is derived by exploiting the fact that the 3D x-ray transform of a function is the solution to the second-order partial differential equation first studied by John. The method is proposed for two sampling schemes, one corresponding to a pair of infinite plane detectors and another one corresponding to a cylindrical multi-ring PET scanner. The new FORE-J algorithm has been implemented for this latter geometry and was compared with the approximate Fourier rebinning algorithm FORE and with another exact rebinning algorithm, FOREX. Results with simulated data demonstrate a significant improvement in accuracy compared to FORE, while the reconstruction time is doubled. Compared to FOREX, the FORE-J algorithm is slightly less accurate but more than three times faster.
NASA Astrophysics Data System (ADS)
Holota, Petr; Nesvadba, Otakar
2015-04-01
In this paper the reciprocal distance is used for generating Galerkin's approximations in the weak solution of Neumann's problem that has an important role in Earth's gravity field studies. The reciprocal distance has a natural tie to the fundamental solution of Laplace's partial differential equation and in the paper it is represented by means of an expansion into a series of oblate spheroidal harmonics. Subsequently, the gradient vector of the reciprocal distance is constructed. In the computation of its components the expansion mentioned above is employed. The paper then focuses on the scalar product of reciprocal distance gradients in two different points and in particular on a series representation of a volume integral of the scalar product spread over an unbounded domain given by the exterior of an oblate spheroid (oblate ellipsoid of revolution). The integral yields the entries of Galerkin's matrix. The numerical interpretation of all the expansions used as well as the respective software implementation within the OpenCL framework is treated, which concerns also a numerical evaluation of Legendre functions of a real and an imaginary argument. In parallel an approximate closed formula expressing the entries of Galerkin's matrix (with an accuracy up to terms multiplied by the square of numerical eccentricity) is derived for convenience and comparison. The paper is added extensive numerical examples that illustrate the approach applied and demonstrate the accuracy of the derived formulas. Aspects related to practical applications are discussed.
Ravetto, P.; Sumini, M.; Ganapol, B.D.
1988-01-01
In an attempt to better understand the influence of prompt and delayed neutrons on nuclear reactor dynamics, a continuous slowing down model based on Fermi age theory was developed several years ago. This model was easily incorporated into the one-group diffusion equation and provided a realistic physical picture of how delayed and prompt neutrons slow down and simultaneously diffuse throughout a medium. The model allows for different slowing down times for each delayed neutron group as well as for prompt neutrons and for spectral differences between the two typed of neutrons. Because of its generality, this model serves not only a a useful predictive tool to anticipate reactor transients, but also as an excellent educational tool to demonstrate the effect of delayed neutrons in reactor kinetics. However, because of numerical complications, the slowing down model could not be developed to its full potential. In particular, the major limitation was the inversion of the Laplace transform, which relied on a knowledge of the poles associated with the resulting transformed flux. For this reason, only one group of delayed neutrons and times longer than the slowing down times could be considered. As is shown, the new inversion procedure removes the short time limitation as well as allows for any number of delayed neutron groups. The inversion technique is versatile and is useful in teaching numerical methods in nuclear science.
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.
Fast and Robust Sixth Order Multigrid Computation for 3D Convection Diffusion Equation.
Wang, Yin; Zhang, Jun
2010-10-15
We present a sixth order explicit compact finite difference scheme to solve the three dimensional (3D) convection diffusion equation. We first use multiscale multigrid method to solve the linear systems arising from a 19-point fourth order discretization scheme to compute the fourth order solutions on both the coarse grid and the fine grid. Then an operator based interpolation scheme combined with an extrapolation technique is used to approximate the sixth order accurate solution on the fine grid. Since the multigrid method using a standard point relaxation smoother may fail to achieve the optimal grid independent convergence rate for solving convection diffusion equation with a high Reynolds number, we implement the plane relaxation smoother in the multigrid solver to achieve better grid independency. Supporting numerical results are presented to demonstrate the efficiency and accuracy of the sixth order compact scheme (SOC), compared with the previously published fourth order compact scheme (FOC). PMID:21151737
A CNN-based approach to integrate the 3-D turbolent diffusion equation
NASA Astrophysics Data System (ADS)
Nunnari, G.
2003-04-01
The paper deals with the integration of the 3-D turbulent diffusion equation. This problem is relevant in several application fields including fluid dynamics, air/water pollution, volcanic ash emissions and industrial hazard assessment. As it is well known numerical solution of such a kind of equation is very time consuming even by using modern digital computers and this represents a short-coming for on-line applications. To overcome this drawback a Cellular Neural Network Approach is proposed in this paper. CNN's proposed by Chua and Yang in 1988 are massive parallel analog non-linear circuits with local interconnections between the computing elements that allow very fast distributed computations. Nowadays several producers of semiconductors such as SGS-Thomson are producing on chip CNN's so that their massive use for heavy computing applications is expected in the near future. In the paper the methodological background of the proposed approach will be outlined. Further some results both in terms of accuracy and computation time will be presented also in comparison with traditional three-dimensional computation schemes. Some results obtained to model 3-D pollution problems in the industrial area of Siracusa (Italy), characterised by a large concentration of petrol-chemical plants, will be presented.
A dispersion minimizing scheme for the 3-D Helmholtz equation based on ray theory
NASA Astrophysics Data System (ADS)
Stolk, Christiaan C.
2016-06-01
We develop a new dispersion minimizing compact finite difference scheme for the Helmholtz equation in 2 and 3 dimensions. The scheme is based on a newly developed ray theory for difference equations. A discrete Helmholtz operator and a discrete operator to be applied to the source and the wavefields are constructed. Their coefficients are piecewise polynomial functions of hk, chosen such that phase and amplitude errors are minimal. The phase errors of the scheme are very small, approximately as small as those of the 2-D quasi-stabilized FEM method and substantially smaller than those of alternatives in 3-D, assuming the same number of gridpoints per wavelength is used. In numerical experiments, accurate solutions are obtained in constant and smoothly varying media using meshes with only five to six points per wavelength and wave propagation over hundreds of wavelengths. When used as a coarse level discretization in a multigrid method the scheme can even be used with down to three points per wavelength. Tests on 3-D examples with up to 108 degrees of freedom show that with a recently developed hybrid solver, the use of coarser meshes can lead to corresponding savings in computation time, resulting in good simulation times compared to the literature.
New equations to calculate 3D joint centres in the lower extremities.
Sandau, Martin; Heimbürger, Rikke V; Villa, Chiara; Jensen, Karl E; Moeslund, Thomas B; Aanæs, Henrik; Alkjær, Tine; Simonsen, Erik B
2015-10-01
Biomechanical movement analysis in 3D requires estimation of joint centres in the lower extremities and this estimation is based on extrapolation from markers placed on anatomical landmarks. The purpose of the present study was to quantify the accuracy of three established set of equations and provide new improved equations to predict the joint centre locations. The 'true' joint centres of the knee and ankle joint were obtained in vivo by MRI scans on 10 male subjects whereas the 'true' hip joint centre was obtained in 10 male and 10 female cadavers by CT scans. For the hip joint the errors ranged from 26.7 (8.9) to 29.6 (7.5) mm, for the knee joint 5.8 (3.1) to 22.6 (3.3) mm and for the ankle joint 14.4 (2.2) to 27.0 (4.6) mm. This differed significantly from the improved equations by which the error for the hip joint ranged from 8.2 (3.6) to 11.6 (5.6) mm, for the knee joint from 2.9 (2.1) to 4.7 (2.5) mm and for the ankle joint from 3.4 (1.3) to 4.1 (2.0) mm. The coefficients in the new hip joint equations differed significantly between sexes. This difference depends on anatomical differences of the male and female pelvis. PMID:26320760
NASA Astrophysics Data System (ADS)
Chen, Duan; Cai, Wei; Zinser, Brian; Cho, Min Hyung
2016-09-01
In this paper, we develop an accurate and efficient Nyström volume integral equation (VIE) method for the Maxwell equations for a large number of 3-D scatterers. The Cauchy Principal Values that arise from the VIE are computed accurately using a finite size exclusion volume together with explicit correction integrals consisting of removable singularities. Also, the hyper-singular integrals are computed using interpolated quadrature formulae with tensor-product quadrature nodes for cubes, spheres and cylinders, that are frequently encountered in the design of meta-materials. The resulting Nyström VIE method is shown to have high accuracy with a small number of collocation points and demonstrates p-convergence for computing the electromagnetic scattering of these objects. Numerical calculations of multiple scatterers of cubic, spherical, and cylindrical shapes validate the efficiency and accuracy of the proposed method.
Cari, C. Suparmi, A.
2014-09-30
Dirac equation of 3D harmonics oscillator plus trigonometric Scarf non-central potential for spin symmetric case is solved using supersymmetric quantum mechanics approach. The Dirac equation for exact spin symmetry reduces to Schrodinger like equation. The relativistic energy and wave function for spin symmetric case are simply obtained using SUSY quantum mechanics method and idea of shape invariance.
An unstaggered constrained transport method for the 3D ideal magnetohydrodynamic equations
NASA Astrophysics Data System (ADS)
Helzel, Christiane; Rossmanith, James A.; Taetz, Bertram
2011-05-01
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with nonlinear numerical instabilities. One approach for controlling the discrete divergence is through a so-called constrained transport method, which is based on first predicting a magnetic field through a standard finite volume solver, and then correcting this field through the appropriate use of a magnetic vector potential. In this work we develop a constrained transport method for the 3D ideal MHD equations that is based on a high-resolution wave propagation scheme. Our proposed scheme is the 3D extension of the 2D scheme developed by Rossmanith [J.A. Rossmanith, An unstaggered, high-resolution constrained transport method for magnetohydrodynamic flows, SIAM J. Sci. Comput. 28 (2006) 1766], and is based on the high-resolution wave propagation method of Langseth and LeVeque [J.O. Langseth, R.J. LeVeque, A wave propagation method for threedimensional hyperbolic conservation laws, J. Comput. Phys. 165 (2000) 126]. In particular, in our extension we take great care to maintain the three most important properties of the 2D scheme: (1) all quantities, including all components of the magnetic field and magnetic potential, are treated as cell-centered; (2) we develop a high-resolution wave propagation scheme for evolving the magnetic potential; and (3) we develop a wave limiting approach that is applied during the vector potential evolution, which controls unphysical oscillations in the magnetic field. One of the key numerical difficulties that is novel to 3D is that the transport equation that must be solved for the magnetic vector potential is only weakly hyperbolic. In presenting our numerical algorithm we describe how to numerically handle this problem of weak hyperbolicity, as well as how to choose an appropriate gauge condition. The
NASA Astrophysics Data System (ADS)
Noppel, M.; Vehkamäki, H.; Winkler, P. M.; Kulmala, M.; Wagner, P. E.
2013-10-01
Based on the results of a previous paper [M. Noppel, H. Vehkamäki, P. M. Winkler, M. Kulmala, and P. E. Wagner, J. Chem. Phys. 139, 134107 (2013)], we derive a thermodynamically consistent expression for reversible or minimal work needed to form a dielectric liquid nucleus of a new phase on a charged insoluble conducting sphere within a uniform macroscopic one- or multicomponent mother phase. The currently available model for ion-induced nucleation assumes complete spherical symmetry of the system, implying that the seed ion is immediately surrounded by the condensing liquid from all sides. We take a step further and treat more realistic geometries, where a cap-shaped liquid cluster forms on the surface of the seed particle. We derive the equilibrium conditions for such a cluster. The equalities of chemical potentials of each species between the nucleus and the vapor represent the conditions of chemical equilibrium. The generalized Young equation that relates contact angle with surface tensions, surface excess polarizations, and line tension, also containing the electrical contribution from triple line excess polarization, expresses the condition of thermodynamic equilibrium at three-phase contact line. The generalized Laplace equation gives the condition of mechanical equilibrium at vapor-liquid dividing surface: it relates generalized pressures in neighboring bulk phases at an interface with surface tension, excess surface polarization, and dielectric displacements in neighboring phases with two principal radii of surface curvature and curvatures of equipotential surfaces in neighboring phases at that point. We also re-express the generalized Laplace equation as a partial differential equation, which, along with electrostatic Laplace equations for bulk phases, determines the shape of a nucleus. We derive expressions that are suitable for calculations of the size and composition of a critical nucleus (generalized version of the classical Kelvin-Thomson equation).
Computational time analysis of the numerical solution of 3D electrostatic Poisson's equation
NASA Astrophysics Data System (ADS)
Kamboh, Shakeel Ahmed; Labadin, Jane; Rigit, Andrew Ragai Henri; Ling, Tech Chaw; Amur, Khuda Bux; Chaudhary, Muhammad Tayyab
2015-05-01
3D Poisson's equation is solved numerically to simulate the electric potential in a prototype design of electrohydrodynamic (EHD) ion-drag micropump. Finite difference method (FDM) is employed to discretize the governing equation. The system of linear equations resulting from FDM is solved iteratively by using the sequential Jacobi (SJ) and sequential Gauss-Seidel (SGS) methods, simulation results are also compared to examine the difference between the results. The main objective was to analyze the computational time required by both the methods with respect to different grid sizes and parallelize the Jacobi method to reduce the computational time. In common, the SGS method is faster than the SJ method but the data parallelism of Jacobi method may produce good speedup over SGS method. In this study, the feasibility of using parallel Jacobi (PJ) method is attempted in relation to SGS method. MATLAB Parallel/Distributed computing environment is used and a parallel code for SJ method is implemented. It was found that for small grid size the SGS method remains dominant over SJ method and PJ method while for large grid size both the sequential methods may take nearly too much processing time to converge. Yet, the PJ method reduces computational time to some extent for large grid sizes.
NASA Technical Reports Server (NTRS)
Yokota, Jeffrey W.
1988-01-01
An LU implicit multigrid algorithm is developed to calculate 3-D compressible viscous flows. This scheme solves the full 3-D Reynolds-Averaged Navier-Stokes equation with a two-equation kappa-epsilon model of turbulence. The flow equations are integrated by an efficient, diagonally inverted, LU implicit multigrid scheme while the kappa-epsilon equations are solved, uncoupled from the flow equations, by a block LU implicit algorithm. The flow equations are solved within the framework of the multigrid method using a four-grid level W-cycle, while the kappa-epsilon equations are iterated only on the finest grid. This treatment of the Reynolds-Averaged Navier-Stokes equations proves to be an efficient method for calculating 3-D compressible viscous flows.
An iterative KP1 method for solving the transport equation in 3D domains on unstructured grids
NASA Astrophysics Data System (ADS)
Kokonkov, N. I.; Nikolaeva, O. V.
2015-10-01
A two-step iterative KP1 method for solving systems of grid equations that approximate the integro-differential transport equation in 3D domains on unstructured grids using nodal SN methods is described. Results of testing the efficiency of the proposed method in solving benchmark problems of reactor protection on tetrahedral grids are presented.
Efficient 3D/1D self-consistent integral-equation analysis of ICRH antennae
NASA Astrophysics Data System (ADS)
Maggiora, R.; Vecchi, G.; Lancellotti, V.; Kyrytsya, V.
2004-08-01
This work presents a comprehensive account of the theory and implementation of a method for the self-consistent numerical analysis of plasma-facing ion-cyclotron resonance heating (ICRH) antenna arrays. The method is based on the integral-equation formulation of the boundary-value problem, solved via a weighted-residual scheme. The antenna geometry (including Faraday shield bars and a recess box) is fairly general and three-dimensional (3D), and the plasma is in the one-dimensional (1D) 'slab' approximation; finite-Larmor radius effects, as well as plasma density and temperature gradients, are considered. Feeding via the voltages in the access coaxial lines is self-consistently accounted throughout and the impedance or scattering matrix of the antenna array obtained therefrom. The problem is formulated in both the dual space (physical) and spectral (wavenumber) domains, which allows the extraction and simple handling of the terms that slow the convergence in the spectral domain usually employed. This paper includes validation tests of the developed code against measured data, both in vacuo and in the presence of plasma. An example of application to a complex geometry is also given.
NASA Astrophysics Data System (ADS)
Abdi, Daniel S.; Giraldo, Francis X.
2016-09-01
A unified approach for the numerical solution of the 3D hyperbolic Euler equations using high order methods, namely continuous Galerkin (CG) and discontinuous Galerkin (DG) methods, is presented. First, we examine how classical CG that uses a global storage scheme can be constructed within the DG framework using constraint imposition techniques commonly used in the finite element literature. Then, we implement and test a simplified version in the Non-hydrostatic Unified Model of the Atmosphere (NUMA) for the case of explicit time integration and a diagonal mass matrix. Constructing CG within the DG framework allows CG to benefit from the desirable properties of DG such as, easier hp-refinement, better stability etc. Moreover, this representation allows for regional mixing of CG and DG depending on the flow regime in an area. The different flavors of CG and DG in the unified implementation are then tested for accuracy and performance using a suite of benchmark problems representative of cloud-resolving scale, meso-scale and global-scale atmospheric dynamics. The value of our unified approach is that we are able to show how to carry both CG and DG methods within the same code and also offer a simple recipe for modifying an existing CG code to DG and vice versa.
SATELLITE DYNAMICS ON THE LAPLACE SURFACE
Tremaine, Scott; Touma, Jihad; Namouni, Fathi E-mail: jihad.touma@gmail.com
2009-03-15
The orbital dynamics of most planetary satellites is governed by the quadrupole moment from the equatorial bulge of the host planet and the tidal field from the Sun. On the Laplace surface, the long-term orbital evolution driven by the combined effects of these forces is zero, so that orbits have a fixed orientation and shape. The 'classical' Laplace surface is defined for circular orbits, and coincides with the planet's equator at small planetocentric distances and with its orbital plane at large distances. A dissipative circumplanetary disk should settle to this surface, and hence satellites formed from such a disk are likely to orbit in or near the classical Laplace surface. This paper studies the properties of Laplace surfaces. Our principal results are: (1) if the planetary obliquity exceeds 68.{sup 0}875, there is a range of semimajor axes in which the classical Laplace surface is unstable; (2) at some obliquities and planetocentric distances, there is a distinct Laplace surface consisting of nested eccentric orbits, which bifurcates from the classical Laplace surface at the point where instability sets in; (3) there is also a 'polar' Laplace surface perpendicular to the line of nodes of the planetary equator on the planetary orbit; (4) for circular orbits, the polar Laplace surface is stable at small planetocentric distances and unstable at large distances; (5) at the onset of instability, this polar Laplace surface bifurcates into two polar Laplace surfaces composed of nested eccentric orbits.
A remark on the Beale-Kato-Majda criterion for the 3D MHD equations with zero magnetic diffusivity
NASA Astrophysics Data System (ADS)
Gala, Sadek; Ragusa, Maria Alessandra
2016-06-01
In this work, we show that a smooth solution of the 3D MHD equations with zero magnetic diffusivity in the whole space ℝ3 breaks down if and only if a certain norm of the magnetic field blows up at the same time.
Solution of the Skyrme HF + BCS equation on a 3D mesh
NASA Astrophysics Data System (ADS)
Bonche, P.; Flocard, H.; Heenen, P. H.
2005-09-01
Over the years, the ev8 code has been a very useful tool for the study of nuclear mean-field theory. Its main characteristic is that it solves the Hartree-Fock plus BCS equations for Skyrme type functionals via a discretization of the individual wave-functions on a three-dimensional Cartesian mesh. This allows maximal flexibility in the determination of the nuclear shape by the variational process. For instance, the same mesh can be used to describe the oblate deformed, spherical, prolate deformed, superdeformed and fission configurations of a given nucleus. The quadrupole constraining operator yielding the deformation energy curve covering all these configurations is included in ev8. This version of the code is restricted to even-even nuclei. Program summaryTitle of program:ev8 Catalogue identifier:ADWA Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADWA Licensing provisions: none Computers on which the program has been tested: HP-RX4640, Compaq-Digital Alpha GS140, has run on several other platforms Computer for which the program is designed and others on which is has been tested:Unix, Linux Operating systems or monitors under which the program has been tested:FORTRAN-90 Programming language used:depends on problem; example given requires 60 MB Memory required to execute with typical data:yes No. of lines in distributed program, including test data, etc.:11 524 No. of bytes in distributed program, including test data, etc.:89 949 Distribution format:tar.gzip file Nature of the physical problem:By means of the Hartree-Fock plus BCS method using Skyrme type functionals, ev8 allows a study of the evolution of the binding energy of even-even nuclei for various shapes determined by the most general quadrupole constraint. Solution method:The program expands the single-particle wave-functions on a 3D Cartesian mesh. The nonlinear mean-field equations are solved by the
On the Global Regularity of a Helical-Decimated Version of the 3D Navier-Stokes Equations
NASA Astrophysics Data System (ADS)
Biferale, Luca; Titi, Edriss S.
2013-06-01
We study the global regularity, for all time and all initial data in H 1/2, of a recently introduced decimated version of the incompressible 3D Navier-Stokes (dNS) equations. The model is based on a projection of the dynamical evolution of Navier-Stokes (NS) equations into the subspace where helicity (the L 2-scalar product of velocity and vorticity) is sign-definite. The presence of a second (beside energy) sign-definite inviscid conserved quadratic quantity, which is equivalent to the H 1/2-Sobolev norm, allows us to demonstrate global existence and uniqueness, of space-periodic solutions, together with continuity with respect to the initial conditions, for this decimated 3D model. This is achieved thanks to the establishment of two new estimates, for this 3D model, which show that the H 1/2 and the time average of the square of the H 3/2 norms of the velocity field remain finite. Such two additional bounds are known, in the spirit of the work of H. Fujita and T. Kato (Arch. Ration. Mech. Anal. 16:269-315, 1964; Rend. Semin. Mat. Univ. Padova 32:243-260, 1962), to be sufficient for showing well-posedness for the 3D NS equations. Furthermore, they are directly linked to the helicity evolution for the dNS model, and therefore with a clear physical meaning and consequences.
Finite Element Code For 3D-Hydraulic Fracture Propagation Equations (3-layer).
Energy Science and Technology Software Center (ESTSC)
1992-03-24
HYFRACP3D is a finite element program for simulation of a pseudo three-dimensional fracture geometries with a two-dimensional planar solution. The model predicts the height, width and winglength over time for a hydraulic fracture propagating in a three-layered system of rocks with variable rock mechanics properties.
NASA Astrophysics Data System (ADS)
Zhdanov, M. S.; Cuma, M.; Black, N.; Wilson, G. A.
2009-12-01
The marine controlled source electromagnetic (MCSEM) method has become widely used in offshore oil and gas exploration. Interpretation of MCSEM data is still a very challenging problem, especially if one would like to take into account the realistic 3D structure of the subsurface. The inversion of MCSEM data is complicated by the fact that the EM response of a hydrocarbon-bearing reservoir is very weak in comparison with the background EM fields generated by an electric dipole transmitter in complex geoelectrical structures formed by a conductive sea-water layer and the terranes beneath it. In this paper, we present a review of the recent developments in the area of large-scale 3D EM forward modeling and inversion. Our approach is based on using a new integral form of Maxwell’s equations allowing for an inhomogeneous background conductivity, which results in a numerically effective integral representation for 3D EM field. This representation provides an efficient tool for the solution of 3D EM inverse problems. To obtain a robust inverse model of the conductivity distribution, we apply regularization based on a focusing stabilizing functional which allows for the recovery of models with both smooth and sharp geoelectrical boundaries. The method is implemented in a fully parallel computer code, which makes it possible to run large-scale 3D inversions on grids with millions of inversion cells. This new technique can be effectively used for active EM detection and monitoring of the subsurface targets.
On the transition towards slow manifold in shallow-water and 3D Euler equations in a rotating frame
NASA Technical Reports Server (NTRS)
Mahalov, A.
1994-01-01
The long-time, asymptotic state of rotating homogeneous shallow-water equations is investigated. Our analysis is based on long-time averaged rotating shallow-water equations describing interactions of large-scale, horizontal, two-dimensional motions with surface inertial-gravity waves field for a shallow, uniformly rotating fluid layer. These equations are obtained in two steps: first by introducing a Poincare/Kelvin linear propagator directly into classical shallow-water equations, then by averaging. The averaged equations describe interaction of wave fields with large-scale motions on time scales long compared to the time scale 1/f(sub o) introduced by rotation (f(sub o)/2-angular velocity of background rotation). The present analysis is similar to the one presented by Waleffe (1991) for 3D Euler equations in a rotating frame. However, since three-wave interactions in rotating shallow-water equations are forbidden, the final equations describing the asymptotic state are simplified considerably. Special emphasis is given to a new conservation law found in the asymptotic state and decoupling of the dynamics of the divergence free part of the velocity field. The possible rising of a decoupled dynamics in the asymptotic state is also investigated for homogeneous turbulence subjected to a background rotation. In our analysis we use long-time expansion, where the velocity field is decomposed into the 'slow manifold' part (the manifold which is unaffected by the linear 'rapid' effects of rotation or the inertial waves) and a formal 3D disturbance. We derive the physical space version of the long-time averaged equations and consider an invariant, basis-free derivation. This formulation can be used to generalize Waleffe's (1991) helical decomposition to viscous inhomogeneous flows (e.g. problems in cylindrical geometry with no-slip boundary conditions on the cylinder surface and homogeneous in the vertical direction).
Xie, G.; Li, J.; Majer, E.; Zuo, D.
1998-07-01
This paper describes a new 3D parallel GILD electromagnetic (EM) modeling and nonlinear inversion algorithm. The algorithm consists of: (a) a new magnetic integral equation instead of the electric integral equation to solve the electromagnetic forward modeling and inverse problem; (b) a collocation finite element method for solving the magnetic integral and a Galerkin finite element method for the magnetic differential equations; (c) a nonlinear regularizing optimization method to make the inversion stable and of high resolution; and (d) a new parallel 3D modeling and inversion using a global integral and local differential domain decomposition technique (GILD). The new 3D nonlinear electromagnetic inversion has been tested with synthetic data and field data. The authors obtained very good imaging for the synthetic data and reasonable subsurface EM imaging for the field data. The parallel algorithm has high parallel efficiency over 90% and can be a parallel solver for elliptic, parabolic, and hyperbolic modeling and inversion. The parallel GILD algorithm can be extended to develop a high resolution and large scale seismic and hydrology modeling and inversion in the massively parallel computer.
Cauchy's almost forgotten Lagrangian formulation of the Euler equation for 3D incompressible flow
NASA Astrophysics Data System (ADS)
Frisch, Uriel; Villone, Barbara
2014-09-01
Two prized papers, one by Augustin Cauchy in 1815, presented to the French Academy and the other by Hermann Hankel in 1861, presented to Göttingen University, contain major discoveries on vorticity dynamics whose impact is now quickly increasing. Cauchy found a Lagrangian formulation of 3D ideal incompressible flow in terms of three invariants that generalize to three dimensions the now well-known law of conservation of vorticity along fluid particle trajectories for two-dimensional flow. This has very recently been used to prove analyticity in time of fluid particle trajectories for 3D incompressible Euler flow and can be extended to compressible flow, in particular to cosmological dark matter. Hankel showed that Cauchy's formulation gives a very simple Lagrangian derivation of the Helmholtz vorticity-flux invariants and, in the middle of the proof, derived an intermediate result which is the conservation of the circulation of the velocity around a closed contour moving with the fluid. This circulation theorem was to be rediscovered independently by William Thomson (Kelvin) in 1869. Cauchy's invariants were only occasionally cited in the 19th century - besides Hankel, foremost by George Stokes and Maurice Lévy - and even less so in the 20th until they were rediscovered via Emmy Noether's theorem in the late 1960, but reattributed to Cauchy only at the end of the 20th century by Russian scientists.
A lattice-Boltzmann scheme of the Navier-Stokes equations on a 3D cuboid lattice
NASA Astrophysics Data System (ADS)
Min, Haoda; Peng, Cheng; Wang, Lian-Ping
2015-11-01
The standard lattice-Boltzmann method (LBM) for fluid flow simulation is based on a square (in 2D) or cubic (in 3D) lattice grids. Recently, two new lattice Boltzmann schemes have been developed on a 2D rectangular grid using the MRT (multiple-relaxation-time) collision model, by adding a free parameter in the definition of moments or by extending the equilibrium moments. Here we developed a lattice Boltzmann model on 3D cuboid lattice, namely, a lattice grid with different grid lengths in different spatial directions. We designed our MRT-LBM model by matching the moment equations from the Chapman-Enskog expansion with the Navier-Stokes equations. The model guarantees correct hydrodynamics. A second-order term is added to the equilibrium moments in order to restore the isotropy of viscosity on a cuboid lattice. The form and the coefficients of the extended equilibrium moments are determined through an inverse design process. An additional benefit of the model is that the viscosity can be adjusted independent of the stress-moment relaxation parameter, thus improving the numerical stability of the model. The resulting cuboid MRT-LBM model is then validated through benchmark simulations using laminar channel flow, turbulent channel flow, and the 3D Taylor-Green vortex flow.
The point-source method for 3D reconstructions for the Helmholtz and Maxwell equations
NASA Astrophysics Data System (ADS)
Ben Hassen, M. F.; Erhard, K.; Potthast, R.
2006-02-01
We use the point-source method (PSM) to reconstruct a scattered field from its associated far field pattern. The reconstruction scheme is described and numerical results are presented for three-dimensional acoustic and electromagnetic scattering problems. We give new proofs of the algorithms, based on the Green and Stratton-Chu formulae, which are more general than with the former use of the reciprocity relation. This allows us to handle the case of limited aperture data and arbitrary incident fields. Both for 3D acoustics and electromagnetics, numerical reconstructions of the field for different settings and with noisy data are shown. For shape reconstruction in acoustics, we develop an appropriate strategy to identify areas with good reconstruction quality and combine different such regions into one joint function. Then, we show how shapes of unknown sound-soft scatterers are found as level curves of the total reconstructed field.
Numerical solution of 3D Navier-Stokes equations with upwind implicit schemes
NASA Technical Reports Server (NTRS)
Marx, Yves P.
1990-01-01
An upwind MUSCL type implicit scheme for the three-dimensional Navier-Stokes equations is presented. Comparison between different approximate Riemann solvers (Roe and Osher) are performed and the influence of the reconstructions schemes on the accuracy of the solution as well as on the convergence of the method is studied. A new limiter is introduced in order to remove the problems usually associated with non-linear upwind schemes. The implementation of a diagonal upwind implicit operator for the three-dimensional Navier-Stokes equations is also discussed. Finally the turbulence modeling is assessed. Good prediction of separated flows are demonstrated if a non-equilibrium turbulence model is used.
Three-dimensional transient electromagnetic modeling in the Laplace Domain
Mizunaga, H.; Lee, Ki Ha; Kim, H.J.
1998-09-01
In modeling electromagnetic responses, Maxwell's equations in the frequency domain are popular and have been widely used (Nabighian, 1994; Newman and Alumbaugh, 1995; Smith, 1996, to list a few). Recently, electromagnetic modeling in the time domain using the finite difference (FDTD) method (Wang and Hohmann, 1993) has also been used to study transient electromagnetic interactions in the conductive medium. This paper presents a new technique to compute the electromagnetic response of three-dimensional (3-D) structures. The proposed new method is based on transforming Maxwell's equations to the Laplace domain. For each discrete Laplace variable, Maxwell's equations are discretized in 3-D using the staggered grid and the finite difference method (FDM). The resulting system of equations is then solved for the fields using the incomplete Cholesky conjugate gradient (ICCG) method. The new method is particularly effective in saving computer memory since all the operations are carried out in real numbers. For the same reason, the computing speed is faster than frequency domain modeling. The proposed approach can be an extremely useful tool in developing an inversion algorithm using the time domain data.
On Energy Cascades in the Forced 3D Navier-Stokes Equations
NASA Astrophysics Data System (ADS)
Dascaliuc, R.; Grujić, Z.
2016-02-01
We show—in the framework of physical scales and (K_1,K_2) -averages—that Kolmogorov's dissipation law combined with the smallness condition on a Taylor length scale is sufficient to guarantee energy cascades in the forced Navier-Stokes equations. Moreover, in the periodic case we establish restrictive scaling laws—in terms of Grashof number—for kinetic energy, energy flux, and energy dissipation rate. These are used to improve our sufficient condition for forced cascades in physical scales.
On Energy Cascades in the Forced 3D Navier-Stokes Equations
NASA Astrophysics Data System (ADS)
Dascaliuc, R.; Grujić, Z.
2016-06-01
We show—in the framework of physical scales and (K_1,K_2)-averages—that Kolmogorov's dissipation law combined with the smallness condition on a Taylor length scale is sufficient to guarantee energy cascades in the forced Navier-Stokes equations. Moreover, in the periodic case we establish restrictive scaling laws—in terms of Grashof number—for kinetic energy, energy flux, and energy dissipation rate. These are used to improve our sufficient condition for forced cascades in physical scales.
Calculations of separated 3-D flows with a pressure-staggered Navier-Stokes equations solver
NASA Technical Reports Server (NTRS)
Kim, S.-W.
1991-01-01
A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.
Ocular surface temperature: a 3D FEM prediction using bioheat equation.
Ng, E Y K; Ooi, E H
2007-06-01
Computational and mathematical human eye models from previous studies which were constructed in two-dimensions (2D) did not give a precise representation of the actual human eye. This work is an extension from an earlier published work on the 2D model. In this paper, a 3D FEM model of the human eye is simulated for the steady state temperature distribution during normal condition and during electromagnetic (EM) wave radiation. Results show a discrepancy of 0.49% for a normal condition as opposed to 1.9% of a 2D model when compared to experimental results from open literatures. Investigations on the EM wave radiations found an average power absorption density of 15,151 and 22,145 Wm(-3) for the 750 and 1500 MHz radiation, respectively. A peak temperature of 38.18( composite function)C was predicted for the 750 MHz radiation while 41.19( composite function)C was computed for the 1500 MHz radiation. These temperatures are in reasonable agreement with the simulated results computed by another report in the past. PMID:17034781
Ott, C D; Dimmelmeier, H; Marek, A; Janka, H-T; Hawke, I; Zink, B; Schnetter, E
2007-06-29
We present 2D and 3D simulations of the collapse of rotating stellar iron cores in general relativity employing a nuclear equation of state and an approximate treatment of deleptonization. We compare fully general relativistic and conformally flat evolutions and find that the latter treatment is sufficiently accurate for the core-collapse supernova problem. We focus on gravitational wave (GW) emission from rotating collapse, bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. We also track the growth of a nonaxisymmetric instability in one model, leading to strong narrow-band GW emission. PMID:17678077
NASA Astrophysics Data System (ADS)
Green, A.; Gribenko, A.; Cuma, M.; Zhdanov, M. S.
2008-12-01
In this paper we apply 3D inversion to MT data collected in Oregon as a part of the EarthScope project. We use the integral equation method as a forward modeling engine. Quasi-analytical approximation with a variable background (QAVB) method of Frechet derivative calculation is applied. This technique allows us to simplify the inversion algorithm and to use just one forward modeling on every iteration step. The receiver footprint approach considerably reduces the computational resources needed to invert the large volumes of data covering vast areas. The data set, which was used in the inversion, was obtained through the Incorporated Research Institutions for Seismology (IRIS). The long-period MT data was collected in Eastern Oregon in 2006. The inverted electrical conductivity distribution agrees reasonably well with geological features of the region as well as with 3D MT inversion results obtained by other researchers. The geoelectrical model of the Oregon deep interior produced by 3D inversion indicates several lithospheres' electrical conductivity anomalies, including a linear zone marked by low-high conductivity transition along the Klamath Blue Mountain Lineament associated with a linear trend of gravity minima. High electrical conductivity values occur in the upper crust under the accreted terrains in the Blue Mountains region.
Local existence and Gevrey regularity of 3-D Navier-Stokes equations with ℓp initial data
NASA Astrophysics Data System (ADS)
Biswas, Animikh
We obtain local existence and Gevrey regularity of 3-D periodic Navier-Stokes equations in case the sequence of Fourier coefficients of the initial data is in ℓp (p<3/2). The ℓp norm of the sequence of Fourier coefficients of the solution and its analogous Gevrey norm remains bounded on a time interval whose length depends only on the size of the body force and the ℓp norm of the Fourier coefficient sequence of the initial data. The control on the Gevrey norm produces explicit estimates on the analyticity radius of the solution as in Foias and Temam (J. Funct. Anal. 87 (1989) 359-369). The results provide an alternate approach in estimating the space-analyticity radius of solutions to Navier-Stokes equations than the one presented by Grujić and Kukavica (J. Funct. Anal. 152 (1998) 447-466).
Recent advances in Runge-Kutta schemes for solving 3-D Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Vatsa, Veer N.; Wedan, Bruce W.; Abid, Ridha
1989-01-01
A thin-layer Navier-Stokes has been developed for solving high Reynolds number, turbulent flows past aircraft components under transonic flow conditions. The computer code has been validated through data comparisons for flow past isolated wings, wing-body configurations, prolate spheroids and wings mounted inside wind-tunnels. The basic code employs an explicit Runge-Kutta time-stepping scheme to obtain steady state solution to the unsteady governing equations. Significant gain in the efficiency of the code has been obtained by implementing a multigrid acceleration technique to achieve steady-state solutions. The improved efficiency of the code has made it feasible to conduct grid-refinement and turbulence model studies in a reasonable amount of computer time. The non-equilibrium turbulence model of Johnson and King has been extended to three-dimensional flows and excellent agreement with pressure data has been obtained for transonic separated flow over a transport type of wing.
Implicit scheme for Maxwell equations solution in case of flat 3D domains
NASA Astrophysics Data System (ADS)
Boronina, Marina; Vshivkov, Vitaly
2016-02-01
We present a new finite-difference scheme for Maxwell's equations solution for three-dimensional domains with different scales in different directions. The stability condition of the standard leap-frog scheme requires decreasing of the time-step with decreasing of the minimal spatial step, which depends on the minimal domain size. We overcome the conditional stability by modifying the standard scheme adding implicitness in the direction of the smallest size. The new scheme satisfies the Gauss law for the electric and magnetic fields in the final- differences. The approximation order, the maintenance of the wave amplitude and propagation speed, the invariance of the wave propagation on angle with the coordinate axes are analyzed.
Absence of Critical Points of Solutions to the Helmholtz Equation in 3D
NASA Astrophysics Data System (ADS)
Alberti, Giovanni S.
2016-05-01
The focus of this paper is to show the absence of critical points for the solutions to the Helmholtz equation in a bounded domain {Ωsubset{R}3} , given by div(a nabla u_{ω}g)-ω qu_{ω}g=0&quad {in Ω,} u_{ω}g=g&quad{on partialΩ.} We prove that for an admissible g there exists a finite set of frequencies K in a given interval and an open cover {overline{Ω}=\\cup_{ωin K} Ω_{ω}} such that {|nabla u_{ω}g(x)| > 0} for every {ωin K} and {xinΩ_{ω}} . The set K is explicitly constructed. If the spectrum of this problem is simple, which is true for a generic domain {Ω} , the admissibility condition on g is a generic property.
Statistical shape analysis using 3D Poisson equation-A quantitatively validated approach.
Gao, Yi; Bouix, Sylvain
2016-05-01
Statistical shape analysis has been an important area of research with applications in biology, anatomy, neuroscience, agriculture, paleontology, etc. Unfortunately, the proposed methods are rarely quantitatively evaluated, and as shown in recent studies, when they are evaluated, significant discrepancies exist in their outputs. In this work, we concentrate on the problem of finding the consistent location of deformation between two population of shapes. We propose a new shape analysis algorithm along with a framework to perform a quantitative evaluation of its performance. Specifically, the algorithm constructs a Signed Poisson Map (SPoM) by solving two Poisson equations on the volumetric shapes of arbitrary topology, and statistical analysis is then carried out on the SPoMs. The method is quantitatively evaluated on synthetic shapes and applied on real shape data sets in brain structures. PMID:26874288
On the Helicity in 3D-Periodic Navier-Stokes Equations II: The Statistical Case
NASA Astrophysics Data System (ADS)
Foias, Ciprian; Hoang, Luan; Nicolaenko, Basil
2009-09-01
We study the asymptotic behavior of the statistical solutions to the Navier-Stokes equations using the normalization map [9]. It is then applied to the study of mean energy, mean dissipation rate of energy, and mean helicity of the spatial periodic flows driven by potential body forces. The statistical distribution of the asymptotic Beltrami flows are also investigated. We connect our mathematical analysis with the empirical theory of decaying turbulence. With appropriate mathematically defined ensemble averages, the Kolmogorov universal features are shown to be transient in time. We provide an estimate for the time interval in which those features may still be present. Our collaborator and friend Basil Nicolaenko passed away in September of 2007, after this work was completed. Honoring his contribution and friendship, we dedicate this article to him.
Soligno, Giuseppe; Dijkstra, Marjolein; van Roij, René
2014-12-28
Many physical problems require explicit knowledge of the equilibrium shape of the interface between two fluid phases. Here, we present a new numerical method which is simply implementable and easily adaptable for a wide range of problems involving capillary deformations of fluid-fluid interfaces. We apply a simulated annealing algorithm to find the interface shape that minimizes the thermodynamic potential of the system. First, for completeness, we provide an analytical proof that minimizing this potential is equivalent to solving the Young-Laplace equation and the Young law. Then, we illustrate our numerical method showing two-dimensional results for fluid-fluid menisci between vertical or inclined walls and curved surfaces, capillary interactions between vertical walls, equilibrium shapes of sessile heavy droplets on a flat horizontal solid surface, and of droplets pending from flat or curved solid surfaces. Finally, we show illustrative three-dimensional results to point out the applicability of the method to micro- or nano-particles adsorbed at a fluid-fluid interface. PMID:25554170
Soligno, Giuseppe; Roij, René van; Dijkstra, Marjolein
2014-12-28
Many physical problems require explicit knowledge of the equilibrium shape of the interface between two fluid phases. Here, we present a new numerical method which is simply implementable and easily adaptable for a wide range of problems involving capillary deformations of fluid-fluid interfaces. We apply a simulated annealing algorithm to find the interface shape that minimizes the thermodynamic potential of the system. First, for completeness, we provide an analytical proof that minimizing this potential is equivalent to solving the Young-Laplace equation and the Young law. Then, we illustrate our numerical method showing two-dimensional results for fluid-fluid menisci between vertical or inclined walls and curved surfaces, capillary interactions between vertical walls, equilibrium shapes of sessile heavy droplets on a flat horizontal solid surface, and of droplets pending from flat or curved solid surfaces. Finally, we show illustrative three-dimensional results to point out the applicability of the method to micro- or nano-particles adsorbed at a fluid-fluid interface.
Calculation of a simulated 3-D high speed inlet using the Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Knight, D. D.
1983-01-01
A hybrid numerical algorithm, developed to solve the full three-dimensional Navier-Stokes equations, is applied to the computation of the flowfield in a simulated three-dimensional high speed aircraft inlet at a Mach number of 2.5 and Reynolds number of 1.4 x 10 to the 7th based on inlet length. The numerical algorithm incorporates a coordinate transformation in order to handle general flow geometries, and utilizes the algebraic turbulent eddy viscosity model of Baldwin and Lomax. The hybrid algorithm has been vectorized on the CDC CYBER 203 computer using the SL/1 vector programming language developed at NASA Langley. The computed results are compared with experimental measurements of the ramp and cowl static pressures, and boundary layer pitot profiles. The results are also compared with a previous two-dimensional Navier-Stokes computation of the same configuration. The agreement with the experimental data is generally good; however, additional improvements in turbulence modeling are needed.
NASA Astrophysics Data System (ADS)
Li, H.; Feng, X. S.; Xiang, J.; Zuo, P.
2014-12-01
In Li et al. [2013, New approach for solving the inverse boundary value problem of Laplace's equation on a circle: Technique renovation of the Grad-Shafranov (GS) reconstruction, J. Geophys. Res. Space., 118, 2876-2881], a couple of Hilbert transform relations were applied to the study of the ill-posedness for the essential GS reconstructions. In this further study, a detailed derivation for these reciprocal relations are presented in case of the plane circular region, and then the reciprocal relations are extended to apply to the plane rectangular region after a conformal mapping procedure. While for the case of plane rectangular region, it is confronted by a traditional problem of the so-called corner singularities, which divided the extended reciprocal relations into four integrals with end-point singularities. With the help of the extended Euler-Maclaurin expansion, new quadrature schemes are developed for these singular integrals. Benchmark testing with the analytic solutions on a rectangle boundary has also show the efficiency and robustness of these extensions. The new solution approach is also developed with the introduced reciprocal relations, and an iterated Tikhonov regularization scheme is applied to deal with the ill-posed linear operators appearing in the discretization of the new approach. The special case on the rectangular boundary is benchmarked with the analytic solutions. Numerical experiments highlight the efficiency and robustness of the proposed method. A robust solution approach is expected to be developed based on these new results for the GS equation on any 2D region with partial-known boundary conditions.
NASA Astrophysics Data System (ADS)
Pankratov, Oleg; Kuvshinov, Alexey
2016-01-01
second part, we summarize modern trends in the development of efficient 3-D EM forward modelling schemes with special emphasis on recent advances in the integral equation approach.
NASA Astrophysics Data System (ADS)
Shi, Jian; Zhang, Qian
2016-03-01
A uniqueness result of weak solution for the 3D viscous magneto-hydrodynamics equations in {B^1_{infty,infty}} is proved by means of the Fourier localization technique and the losing derivative estimates.
NASA Astrophysics Data System (ADS)
Ge, Liang; Sotiropoulos, Fotis
2007-08-01
A novel numerical method is developed that integrates boundary-conforming grids with a sharp interface, immersed boundary methodology. The method is intended for simulating internal flows containing complex, moving immersed boundaries such as those encountered in several cardiovascular applications. The background domain (e.g. the empty aorta) is discretized efficiently with a curvilinear boundary-fitted mesh while the complex moving immersed boundary (say a prosthetic heart valve) is treated with the sharp-interface, hybrid Cartesian/immersed-boundary approach of Gilmanov and Sotiropoulos [A. Gilmanov, F. Sotiropoulos, A hybrid cartesian/immersed boundary method for simulating flows with 3d, geometrically complex, moving bodies, Journal of Computational Physics 207 (2005) 457-492.]. To facilitate the implementation of this novel modeling paradigm in complex flow simulations, an accurate and efficient numerical method is developed for solving the unsteady, incompressible Navier-Stokes equations in generalized curvilinear coordinates. The method employs a novel, fully-curvilinear staggered grid discretization approach, which does not require either the explicit evaluation of the Christoffel symbols or the discretization of all three momentum equations at cell interfaces as done in previous formulations. The equations are integrated in time using an efficient, second-order accurate fractional step methodology coupled with a Jacobian-free, Newton-Krylov solver for the momentum equations and a GMRES solver enhanced with multigrid as preconditioner for the Poisson equation. Several numerical experiments are carried out on fine computational meshes to demonstrate the accuracy and efficiency of the proposed method for standard benchmark problems as well as for unsteady, pulsatile flow through a curved, pipe bend. To demonstrate the ability of the method to simulate flows with complex, moving immersed boundaries we apply it to calculate pulsatile, physiological flow
NASA Technical Reports Server (NTRS)
Demuren, A. O.; Ibraheem, S. O.
1993-01-01
The convergence characteristics of various approximate factorizations for the 3D Euler and Navier-Stokes equations are examined using the von-Neumann stability analysis method. Three upwind-difference based factorizations and several central-difference based factorizations are considered for the Euler equations. In the upwind factorizations both the flux-vector splitting methods of Steger and Warming and van Leer are considered. Analysis of the Navier-Stokes equations is performed only on the Beam and Warming central-difference scheme. The range of CFL numbers over which each factorization is stable is presented for one-, two-, and three-dimensional flow. Also presented for each factorization is the CFL number at which the maximum eigenvalue is minimized, for all Fourier components, as well as for the high frequency range only. The latter is useful for predicting the effectiveness of multigrid procedures with these schemes as smoothers. Further, local mode analysis is performed to test the suitability of using a uniform flow field in the stability analysis. Some inconsistencies in the results from previous analyses are resolved.
NASA Astrophysics Data System (ADS)
Newton, W. G.; Stone, J. R.; Mezzacappa, A.
2006-09-01
First results from a fully self-consistent, temperature-dependent equation of state that spans the density range of neutron stars and supernova cores above neutron drip density are presented. The equation of state (EoS) is calculated using a mean-field Hartree-Fock method in three dimensions (3D). The nuclear interaction is represented by the phenomenological Skyrme model in this work, but the EoS can be obtained in our framework for any suitable form of the nucleon-nucleon effective interaction. The scheme we employ naturally allows effects such as (i) neutron drip, which results in an external neutron gas, (ii) the variety of exotic nuclear shapes expected for extremely neutron heavy nuclei, and (iii) the subsequent dissolution of these nuclei into nuclear matter. In this way, the equation of state is calculated across phase transitions without recourse to interpolation techniques between density regimes described by different physical models. EoS tables are calculated in the wide range of densities, temperature and proton/neutron ratios on the ORNL NCCS XT3, using up to 2000 processors simultaneously.
An IPOT meshless method using DC PSE approximation for fluid flow equations in 2D and 3D geometries
NASA Astrophysics Data System (ADS)
Bourantas, G. C.; Loukopoulos, V. C.; Skouras, E. D.; Burganos, V. N.; Nikiforidis, G. C.
2016-06-01
Navier-Stokes (N-S) equations, in their primitive variable (u-v-p) formulation, are numerically solved using the Implicit Potential (IPOT) numerical scheme in the context of strong form Meshless Point Collocation (MPC) method. The unknown field functions are computed using the Discretization Correction Particle Strength Exchange (DC PSE) approximation method. The latter makes use of discrete moment conditions to derive the operator kernels, which leads to low condition number for the moment matrix compared to other meshless interpolation methods and increased stability for the numerical solution. The proposed meshless scheme is applied on 2D and 3D spatial domains, using uniform or irregular set of nodes to represent the domain. The numerical results obtained are compared against those obtained using well-established methods.
NASA Astrophysics Data System (ADS)
Simpson, J. J.; Taflove, A.
2005-12-01
We report a finite-difference time-domain (FDTD) computational solution of Maxwell's equations [1] that models the possibility of detecting and characterizing ionospheric disturbances above seismic regions. Specifically, we study anomalies in Schumann resonance spectra in the extremely low frequency (ELF) range below 30 Hz as observed in Japan caused by a hypothetical cylindrical ionospheric disturbance above Taiwan. We consider excitation of the global Earth-ionosphere waveguide by lightning in three major thunderstorm regions of the world: Southeast Asia, South America (Amazon region), and Africa. Furthermore, we investigate varying geometries and characteristics of the ionospheric disturbance above Taiwan. The FDTD technique used in this study enables a direct, full-vector, three-dimensional (3-D) time-domain Maxwell's equations calculation of round-the-world ELF propagation accounting for arbitrary horizontal as well as vertical geometrical and electrical inhomogeneities and anisotropies of the excitation, ionosphere, lithosphere, and oceans. Our entire-Earth model grids the annular lithosphere-atmosphere volume within 100 km of sea level, and contains over 6,500,000 grid-points (63 km laterally between adjacent grid points, 5 km radial resolution). We use our recently developed spherical geodesic gridding technique having a spatial discretization best described as resembling the surface of a soccer ball [2]. The grid is comprised entirely of hexagonal cells except for a small fixed number of pentagonal cells needed for completion. Grid-cell areas and locations are optimized to yield a smoothly varying area difference between adjacent cells, thereby maximizing numerical convergence. We compare our calculated results with measured data prior to the Chi-Chi earthquake in Taiwan as reported by Hayakawa et. al. [3]. Acknowledgement This work was suggested by Dr. Masashi Hayakawa, University of Electro-Communications, Chofugaoka, Chofu Tokyo. References [1] A
On the Finite-Time Splash and Splat Singularities for the 3-D Free-Surface Euler Equations
NASA Astrophysics Data System (ADS)
Coutand, Daniel; Shkoller, Steve
2014-01-01
We prove that the 3-D free-surface incompressible Euler equations with regular initial geometries and velocity fields have solutions which can form a finite-time "splash" (or "splat") singularity first introduced in Castro et al. (Splash singularity for water waves, http://arxiv.org/abs/1106.2120v2, 2011), wherein the evolving 2-D hypersurface, the moving boundary of the fluid domain, self-intersects at a point (or on surface). Such singularities can occur when the crest of a breaking wave falls unto its trough, or in the study of drop impact upon liquid surfaces. Our approach is founded upon the Lagrangian description of the free-boundary problem, combined with a novel approximation scheme of a finite collection of local coordinate charts; as such we are able to analyze a rather general set of geometries for the evolving 2-D free-surface of the fluid. We do not assume the fluid is irrotational, and as such, our method can be used for a number of other fluid interface problems, including compressible flows, plasmas, as well as the inclusion of surface tension effects.
NASA Astrophysics Data System (ADS)
Gainullin, I. K.; Sonkin, M. A.
2015-03-01
A parallelized three-dimensional (3D) time-dependent Schrodinger equation (TDSE) solver for one-electron systems is presented in this paper. The TDSE Solver is based on the finite-difference method (FDM) in Cartesian coordinates and uses a simple and explicit leap-frog numerical scheme. The simplicity of the numerical method provides very efficient parallelization and high performance of calculations using Graphics Processing Units (GPUs). For example, calculation of 106 time-steps on the 1000ṡ1000ṡ1000 numerical grid (109 points) takes only 16 hours on 16 Tesla M2090 GPUs. The TDSE Solver demonstrates scalability (parallel efficiency) close to 100% with some limitations on the problem size. The TDSE Solver is validated by calculation of energy eigenstates of the hydrogen atom (13.55 eV) and affinity level of H- ion (0.75 eV). The comparison with other TDSE solvers shows that a GPU-based TDSE Solver is 3 times faster for the problems of the same size and with the same cost of computational resources. The usage of a non-regular Cartesian grid or problem-specific non-Cartesian coordinates increases this benefit up to 10 times. The TDSE Solver was applied to the calculation of the resonant charge transfer (RCT) in nanosystems, including several related physical problems, such as electron capture during H+-H0 collision and electron tunneling between H- ion and thin metallic island film.
NASA Astrophysics Data System (ADS)
Zhai, Cuili; Zhang, Ting
2015-09-01
In this article, we consider the global well-posedness to the 3-D incompressible inhomogeneous Navier-Stokes equations with a class of large velocity. More precisely, assuming a 0 ∈ B˙ q , 1 /3 q ( R 3 ) and u 0 = ( u0 h , u0 3 ) ∈ B˙ p , 1 - 1 + /3 p ( R 3 ) for p, q ∈ (1, 6) with sup ( /1 p , /1 q ) ≤ /1 3 + inf ( /1 p , /1 q ) , we prove that if C a↑0↑ B˙q1/3 q α (↑u0 3↑ B˙ p , 1 - 1 + /3 p/μ + 1 ) ≤ 1 , /C μ (↑u0 h↑ B˙ p , 1 - 1 + /3 p + ↑u03↑ B˙ p , 1 - 1 + /3 p 1 - α ↑u0h↑ B˙ p , 1 - 1 + /3 p α) ≤ 1 , then the system has a unique global solution a ∈ C ˜ ( [ 0 , ∞ ) ; B˙ q , 1 /3 q ( R 3 ) ) , u ∈ C ˜ ( [ 0 , ∞ ) ; B˙ p , 1 - 1 + /3 p ( R 3 ) ) ∩ L 1 ( R + ; B˙ p , 1 1 + /3 p ( R 3 ) ) . It improves the recent result of M. Paicu and P. Zhang [J. Funct. Anal. 262, 3556-3584 (2012)], where the exponent form of the initial smallness condition is replaced by a polynomial form.
NASA Astrophysics Data System (ADS)
Zieniuk, Eugeniusz; Kapturczak, Marta; Sawicki, Dominik
2016-06-01
In solving of boundary value problems the shapes of the boundary can be modelled by the curves widely used in computer graphics. In parametric integral equations system (PIES) such curves are directly included into the mathematical formalism. Its simplify the way of definition and modification of the shape of the boundary. Until now in PIES the B-spline, Bézier and Hermite curves were used. Recent developments in the computer graphics paid our attention, therefore we implemented in PIES possibility of defining the shape of boundary using the NURBS curves. The curves will allow us to modeling different shapes more precisely. In this paper we will compare PIES solutions (with applied NURBS) with the solutions existing in the literature.
ERIC Educational Resources Information Center
Merchant, Zahira; Goetz, Ernest T.; Keeney-Kennicutt, Wendy; Kwok, Oi-man; Cifuentes, Lauren; Davis, Trina J.
2012-01-01
We examined a model of the impact of a 3D desktop virtual reality environment on the learner characteristics (i.e. perceptual and psychological variables) that can enhance chemistry-related learning achievements in an introductory college chemistry class. The relationships between the 3D virtual reality features and the chemistry learning test as…
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.
1990-01-01
The development and applications of multiblock/multizone and adaptive grid methodologies for solving the three-dimensional simplified Navier-Stokes equations are described. Adaptive grid and multiblock/multizone approaches are introduced and applied to external and internal flow problems. These new implementations increase the capabilities and flexibility of the PAB3D code in solving flow problems associated with complex geometry.
NASA Astrophysics Data System (ADS)
Bustamante, Miguel D.
2014-11-01
We consider 3D Euler fluids endowed with a discrete symmetry whereby the velocity field is invariant under mirror reflections about a 2D surface known as the ``symmetry plane.'' This type of flow is widely used in numerical simulations of classical/magnetic/quantum turbulence and vortex reconnection. On the 2D symmetry plane, the governing equations are best written in terms of two scalars: vorticity and stretching rate of vorticity. These determine the velocity field on the symmetry plane. However, the governing equations are not closed, because of the contribution of a single pressure term that depends on the full 3D velocity profile. By modelling this pressure term we propose a one-parameter family of sensible models for the flow along the 2D symmetry plane. We apply the method of infinitesimal Lie symmetries and solve the governing equations analytically for the two scalars as functions of time. We show how the value of the model's parameter determines if the analytical solution has a finite-time blowup and obtain explicit formulae for the blowup time. We validate the models by showing that a particular choice of the model's parameter corresponds to a well-known exact solution of 3D Euler equations [Gibbon et al., Physica D 132, 497 (1999)]. We discuss practical applications. Supported by Science Foundation Ireland (SFI) under Grant Number 12/IP/1491.
NASA Technical Reports Server (NTRS)
Kwak, D.
1994-01-01
INS3D computes steady-state solutions to the incompressible Navier-Stokes equations. The INS3D approach utilizes pseudo-compressibility combined with an approximate factorization scheme. This computational fluid dynamics (CFD) code has been verified on problems such as flow through a channel, flow over a backwardfacing step and flow over a circular cylinder. Three dimensional cases include flow over an ogive cylinder, flow through a rectangular duct, wind tunnel inlet flow, cylinder-wall juncture flow and flow through multiple posts mounted between two plates. INS3D uses a pseudo-compressibility approach in which a time derivative of pressure is added to the continuity equation, which together with the momentum equations form a set of four equations with pressure and velocity as the dependent variables. The equations' coordinates are transformed for general three dimensional applications. The equations are advanced in time by the implicit, non-iterative, approximately-factored, finite-difference scheme of Beam and Warming. The numerical stability of the scheme depends on the use of higher-order smoothing terms to damp out higher-frequency oscillations caused by second-order central differencing. The artificial compressibility introduces pressure (sound) waves of finite speed (whereas the speed of sound would be infinite in an incompressible fluid). As the solution converges, these pressure waves die out, causing the derivation of pressure with respect to time to approach zero. Thus, continuity is satisfied for the incompressible fluid in the steady state. Computational efficiency is achieved using a diagonal algorithm. A block tri-diagonal option is also available. When a steady-state solution is reached, the modified continuity equation will satisfy the divergence-free velocity field condition. INS3D is capable of handling several different types of boundaries encountered in numerical simulations, including solid-surface, inflow and outflow, and far
Kamon, M.; Phillips, J.R.
1994-12-31
In this paper techniques are presented for preconditioning equations generated by discretizing constrained vector integral equations associated with magnetoquasistatic analysis. Standard preconditioning approaches often fail on these problems. The authors present a specialized preconditioning technique and prove convergence bounds independent of the constraint equations and electromagnetic excitation frequency. Computational results from analyzing several electronic packaging examples are given to demonstrate that the new preconditioning approach can sometimes reduce the number of GMRES iterations by more than an order of magnitude.
Efficient solution of 3D Ginzburg-Landau problem for mesoscopic superconductors
NASA Astrophysics Data System (ADS)
Pereira, Paulo J.; Moshchalkov, Victor V.; Chibotaru, Liviu F.
2014-03-01
The recently proposed approach for the solution of Ginzburg-Landau (GL) problem for 2D samples of arbitrary shape is, in this article, extended over 3D samples having the shape of (i) a prism with arbitrary base and (ii) a solid of revolution with arbitrary profile. Starting from the set of Laplace operator eigenfunctions of a 2D object, we construct an approximation to or the exact eigenfunctions of the Laplace operator of a 3D structure by applying an extrusion or revolution to these solutions. This set of functions is used as the basis to construct the solutions of the linearized GL equation. These solutions are then used as basis for the non-linear GL equation much like the famous LCAO method. To solve the non-linear equation, we used the Newton-Raphson method starting from the solution of the linear equation, i.e., the nucleation distribution of superconducting condensate. The vector potential approximations typically used in 2D cases, i.e., considering it as corresponding to applied constant field, are in the 3D case harder to justify. For that reason, we use a locally corrected Nystrom method to solve the second Ginzburg-Landau equation. The complete solution of GL problem is then achieved by solving self-consistently both equations.
NASA Astrophysics Data System (ADS)
Vidal, A.; San-Blas, A. A.; Quesada-Pereira, F. D.; Pérez-Soler, J.; Gil, J.; Vicente, C.; Gimeno, B.; Boria, V. E.
2015-07-01
A novel technique for the full-wave analysis of 3-D complex waveguide devices is presented. This new formulation, based on the Boundary Integral-Resonant Mode Expansion (BI-RME) method, allows the rigorous full-wave electromagnetic characterization of 3-D arbitrarily shaped metallic structures making use of extremely low CPU resources (both time and memory). The unknown electric current density on the surface of the metallic elements is represented by means of Rao-Wilton-Glisson basis functions, and an algebraic procedure based on a singular value decomposition is applied to transform such functions into the classical solenoidal and nonsolenoidal basis functions needed by the original BI-RME technique. The developed tool also provides an accurate computation of the electromagnetic fields at an arbitrary observation point of the considered device, so it can be used for predicting high-power breakdown phenomena. In order to validate the accuracy and efficiency of this novel approach, several new designs of band-pass waveguides filters are presented. The obtained results (S-parameters and electromagnetic fields) are successfully compared both to experimental data and to numerical simulations provided by a commercial software based on the finite element technique. The results obtained show that the new technique is specially suitable for the efficient full-wave analysis of complex waveguide devices considering an integrated coaxial excitation, where the coaxial probes may be in contact with the metallic insets of the component.
On Bi-Grid Local Mode Analysis of Solution Techniques for 3-D Euler and Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Ibraheem, S. O.; Demuren, A. O.
1994-01-01
A procedure is presented for utilizing a bi-grid stability analysis as a practical tool for predicting multigrid performance in a range of numerical methods for solving Euler and Navier-Stokes equations. Model problems based on the convection, diffusion and Burger's equation are used to illustrate the superiority of the bi-grid analysis as a predictive tool for multigrid performance in comparison to the smoothing factor derived from conventional von Neumann analysis. For the Euler equations, bi-grid analysis is presented for three upwind difference based factorizations, namely Spatial, Eigenvalue and Combination splits, and two central difference based factorizations, namely LU and ADI methods. In the former, both the Steger-Warming and van Leer flux-vector splitting methods are considered. For the Navier-Stokes equations, only the Beam-Warming (ADI) central difference scheme is considered. In each case, estimates of multigrid convergence rates from the bi-grid analysis are compared to smoothing factors obtained from single-grid stability analysis. Effects of grid aspect ratio and flow skewness are examined. Both predictions are compared with practical multigrid convergence rates for 2-D Euler and Navier-Stokes solutions based on the Beam-Warming central scheme.
NASA Astrophysics Data System (ADS)
Young, D. L.; Tsai, C. H.; Wu, C. S.
2015-11-01
An alternative vector potential formulation is used to solve the Navier-Stokes (N-S) equations in 3D incompressible viscous flow problems with and without through-flow boundaries. Difficulties of the vector potential formulation include the implementation of boundary conditions for through-flow boundaries and the numerical treatment of fourth-order partial differential equations. The advantages on the other hand are the automatic satisfaction of the continuity equation; and pressure is decoupled from the velocity. The objective of this paper is to introduce the appropriate gauge and boundary conditions on the vector potential formulation by a localized meshless method. To handle the divergence-free property, a Coulomb gauge condition is enforced on the vector potential to ensure its existence and uniqueness mathematically. We further improve the algorithm to through-flow problems for the boundary conditions of vector potential by introducing the concept of Stokes' theorem. Based on this innovation, there is no need to include an additional variable to tackle the through-flow fields. This process will greatly simplify the imposition of boundary conditions by the vector potential approach. Under certain conditions, the coupled fourth-order partial differential equations can be easily solved by using this meshless local differential quadrature (LDQ) method. Due to the LDQ capability to deal with the high order differential equations, this algorithm is very attractive to solve this fourth-order vector potential formulation for the N-S equations as comparing to the conventional numerical schemes such as finite element or finite difference methods. The proposed vector potential formulation is simpler and has improved accuracy and efficiency compared to other pressure-free or pressure-coupled algorithms. This investigation can be regarded as the first complete study to obtain the N-S solutions by vector potential formulation through a LDQ method. Two classic 3D benchmark
On the Rigorous Derivation of the 3D Cubic Nonlinear Schrödinger Equation with a Quadratic Trap
NASA Astrophysics Data System (ADS)
Chen, Xuwen
2013-11-01
We consider the dynamics of the three-dimensional N-body Schrödinger equation in the presence of a quadratic trap. We assume the pair interaction potential is N 3 β-1 V( N β x). We justify the mean-field approximation and offer a rigorous derivation of the three-dimensional cubic nonlinear Schrödinger equation (NLS) with a quadratic trap. We establish the space-time bound conjectured by Klainerman and Machedon (Commun Math Phys 279:169-185, 2008) for by adapting and simplifying an argument in Chen and Pavlović (Annales Henri Poincaré, 2013) which solves the problem for in the absence of a trap.
NASA Astrophysics Data System (ADS)
Yuasa, T.; Sunaguchi, N.; Ichihara, S.; Ando, M.
2013-05-01
Refraction-contrast computed tomography (CT) can image biological soft tissues and soft materials at a high contrast, which cannot be clearly depicted by contemporary CT based on absorption contrast. It reconstructs a distribution of refractive index from projections, whose data each is an angular deviation from incident direction due to refraction by an object, and is acquired by imaging methods using an angular analyzer, e.g., DEI (diffraction enhance imaging), or DFI (dark field imaging). First, a reconstruction algorithm for refraction-contrast CT is derived from the ray equation of a fundamental equation describing refraction phenomena in geometrical optics. Then, in order to demonstrate its efficacy, we performed imaging experiment using DFI-CT imaging system. A reconstructed image of human breast cancer tissue is presented.
NASA Astrophysics Data System (ADS)
Tsuzuki, Yutaka
2015-09-01
This paper is concerned with a system of heat equations with hysteresis and Navier-Stokes equations. In Tsuzuki (J Math Anal Appl 423:877-897, 2015) an existence result is obtained for the problem in a 2-dimensional domain with the Navier-Stokes equation in a weak sense. However the result does not include uniqueness for the problem due to the low regularity for solutions. This paper establishes existence and uniqueness in 2- and 3-dimensional domains with the Navier-Stokes equation in a stronger sense. Moreover this work decides required height of regularity for the initial data by introducing the fractional power of the Stokes operator.
NASA Technical Reports Server (NTRS)
Zhang, Jun; Ge, Lixin; Kouatchou, Jules
2000-01-01
A new fourth order compact difference scheme for the three dimensional convection diffusion equation with variable coefficients is presented. The novelty of this new difference scheme is that it Only requires 15 grid points and that it can be decoupled with two colors. The entire computational grid can be updated in two parallel subsweeps with the Gauss-Seidel type iterative method. This is compared with the known 19 point fourth order compact differenCe scheme which requires four colors to decouple the computational grid. Numerical results, with multigrid methods implemented on a shared memory parallel computer, are presented to compare the 15 point and the 19 point fourth order compact schemes.
OpenMP for 3D potential boundary value problems solved by PIES
NASA Astrophysics Data System (ADS)
KuŻelewski, Andrzej; Zieniuk, Eugeniusz
2016-06-01
The main purpose of this paper is examination of an application of modern parallel computing technique OpenMP to speed up the calculation in the numerical solution of parametric integral equations systems (PIES). The authors noticed, that solving more complex boundary problems by PIES sometimes requires large computing time. This paper presents the use of OpenMP and fast C++ linear algebra library Armadillo for boundary value problems modelled by 3D Laplace's equation and solved using PIES. The testing example shows that the use of mentioned technologies significantly increases speed of calculations in PIES.
Models Ion Trajectories in 2D and 3D Electrostatic and Magnetic Fields
Energy Science and Technology Software Center (ESTSC)
2000-02-21
SIMION3D7.0REV is a C based ion optics simulation program that can model complex problems using Laplace equation solutions for potential fields. The program uses an ion optics workbench that can hold up to 200 2D and/or 3D electrostatic/magnetic potential arrays. Arrays can have up to 50,000,000 points. SIMION3D7.0''s 32 bit virtual Graphics User Interface provides a highly interactive advanced user environment. All potential arrays are visualized as 3D objects that the user can cut awaymore » to inspect ion trajectories and potential energy surfaces. User programs allow the user to customize the program for specific simulations. A geometry file option supports the definition of highly complex array geometry. Algorithm modifications have improved this version''s computational speed and accuracy.« less
Models Ion Trajectories in 2D and 3D Electrostatic and Magnetic Fields
Dahl, David
2000-02-21
SIMION3D7.0REV is a C based ion optics simulation program that can model complex problems using Laplace equation solutions for potential fields. The program uses an ion optics workbench that can hold up to 200 2D and/or 3D electrostatic/magnetic potential arrays. Arrays can have up to 50,000,000 points. SIMION3D7.0''s 32 bit virtual Graphics User Interface provides a highly interactive advanced user environment. All potential arrays are visualized as 3D objects that the user can cut away to inspect ion trajectories and potential energy surfaces. User programs allow the user to customize the program for specific simulations. A geometry file option supports the definition of highly complex array geometry. Algorithm modifications have improved this version''s computational speed and accuracy.
The duality principle and inversion of Laplace-Stielties transforms
NASA Astrophysics Data System (ADS)
Pavelyev, A. G.
2016-04-01
The fundamental relation between the Laplace transform, the Stielties transform, and the generalized integral equation of refraction is revealed, and a duality principle is formulated for the solution of inverse problems of radio physics. New formulas of the Laplace-transform inversion satisfying the duality principle are obtained. There is no necessity of contour integration in a complex plane for the relations found, which considerably simplifies the reconstruction of originals and makes it possible to control systematic errors in the experimental data.
An Operator Method for Evaluating Laplace Transforms
ERIC Educational Resources Information Center
Lanoue, B. G.; Yurekli, O.
2005-01-01
This note discusses a simple operator technique based on the differentiation and shifting properties of the Laplace transform to find Laplace transforms for various elementary functions. The method is simpler than known integration techniques to evaluate Laplace transforms.
NASA Technical Reports Server (NTRS)
Cwik, T.; Jamnejad, V.; Zuffada, C.
1993-01-01
It is often desirable to calculate the electromagnetic fields inside and about a complicated system of scattering bodies, as well as in their far-field region. The finite element method (FE) is well suited to solving the interior problem, but the domain has to be limited to a manageable size. At the truncation of the FE mesh one can either impose approximate (absorbing) boundary conditions or set up an integral equation (IE) for the fields scattered from the bodies. The latter approach is preferable since it results in higher accuracy. Hence, the two techniques can be successfully combined by introducing a surface that encloses the scatterers, applying a FE model to the inner volume and setting up an IE for the tangential fields components on the surface. Here the continuity of the tangential fields is used bo obtain a consistent solution. A few coupled FE-IE methods have recently appeared in the literature. The approach presented here has the advantage of using edge-based finite elements, a type of finite elements with degrees of freedom associated with edges of the mesh. Because of their properties, they are better suited than the conventional node based elements to represent electromagnetic fields, particularly when inhomogeneous regions are modeled, since the node based elements impose an unnatural continuity of all field components across boundaries of mesh elements. Additionally, our approach is well suited to handle large size problems and lends itself to code parallelization. We will discuss the salient features that make our approach very efficient from the standpoint of numerical computation, and the fields and RCS of a few objects are illustrated as examples.
Joukar, Amin; Nammakie, Erfan; Niroomand-Oscuii, Hanieh
2015-01-01
The application of laser in ophthalmology and eye surgery is so widespread that hardly can anyone deny its importance. On the other hand, since the human eye is an organ susceptible to external factors such as heat waves, laser radiation rapidly increases the temperature of the eye and therefore the study of temperature distribution inside the eye under laser irradiation is crucial; but the use of experimental and invasive methods for measuring the temperature inside the eye is typically high-risk and hazardous. In this paper, using the three-dimensional finite element method, the distribution of heat transfer inside the eye under transient condition was studied through three different lasers named Nd:Yag, Nd:Yap and ArF. Considering the metabolic heat and blood perfusion rate in various regions of the eye, numerical solution of space-time dependant Pennes bioheat transfer equation has been applied in this study. Lambert-Beer's law has been used to model the absorption of laser energy inside the eye tissues. It should also be mentioned that the effect of the ambient temperature, tear evaporation rate, laser power and the pupil diameter on the temperature distribution have been studied. Also, temperature distribution inside the eye after applying each laser and temperature variations of six optional regions as functions of time have been investigated. The results show that these radiations cause temperature rise in various regions, which will in turn causes serious damages to the eye tissues. Investigating the temperature distribution inside the eye under the laser irradiation can be a useful tool to study and predict the thermal effects of laser radiation on the human eye and evaluate the risk involved in performing laser surgery. PMID:25774029
The Laplace transform on time scales revisited
NASA Astrophysics Data System (ADS)
Davis, John M.; Gravagne, Ian A.; Jackson, Billy J.; Marks, Robert J., II; Ramos, Alice A.
2007-08-01
In this work, we reexamine the time scale Laplace transform as defined by Bohner and Peterson [M. Bohner, A. Peterson, Dynamic Equations on Time Scales: An Introduction with Applications, Birkhauser, Boston, 2001; M. Bohner, A. Peterson, Laplace transform and Z-transform: Unification and extension, Methods Appl. Anal. 9 (1) (2002) 155-162]. In particular, we give conditions on the class of functions which have a transform, develop an inversion formula for the transform, and further, we provide a convolution for the transform. The notion of convolution leads to considering its algebraic structure--in particular the existence of an identity element--motivating the development of the Dirac delta functional on time scales. Applications and examples of these concepts are given.
NASA Astrophysics Data System (ADS)
Kolybasova, V. V.; Krutitskii, P. A.
2010-09-01
We consider a skew derivative problem for a function which is harmonic in the exterior of open arcs in a plane. This problem models electric current in a semiconductor film from electrodes of arbitrary shapes in the presence of a magnetic field. A numerical method for solving the problem is proposed. The method is based on the boundary integral equation approach. The proposed numerical method is tested for different values of parameters and different shapes of the electrodes.
NASA Astrophysics Data System (ADS)
Porter, K.
2015-12-01
There are two common ways to create a ground-motion map for a hypothetical earthquake: using ground motion prediction equations (by far the more common of the two) and using 3-D physics-based modeling. The former is very familiar to engineers, the latter much less so, and the difference can present a problem because engineers tend to trust the familiar and distrust novelty. Maps for essentially the same hypothetical earthquake using the two different methods can look very different, while appearing to present the same information. Using one or the other can lead an engineer or disaster planner to very different estimates of damage and risk. The reasons have to do with depiction of variability, spatial correlation of shaking, the skewed distribution of real-world shaking, and the upward-curving relationship between shaking and damage. The scientists who develop the two kinds of map tend to specialize in one or the other and seem to defend their turf, which can aggravate the problem of clearly communicating with engineers.The USGS Science Application for Risk Reduction's (SAFRR) HayWired scenario has addressed the challenge of explaining to engineers the differences between the two maps, and why, in a disaster planning scenario, one might want to use the less-familiar 3-D map.
NASA Astrophysics Data System (ADS)
Li, Liang; Lanteri, Stéphane; Perrussel, Ronan
2014-01-01
A Schwarz-type domain decomposition method is presented for the solution of the system of 3d time-harmonic Maxwell's equations. We introduce a hybridizable discontinuous Galerkin (HDG) scheme for the discretization of the problem based on a tetrahedrization of the computational domain. The discrete system of the HDG method on each subdomain is solved by an optimized sparse direct (LU factorization) solver. The solution of the interface system in the domain decomposition framework is accelerated by a Krylov subspace method. The formulation and the implementation of the resulting DD-HDG (Domain Decomposed-Hybridizable Discontinuous Galerkin) method are detailed. Numerical results show that the resulting DD-HDG solution strategy has an optimal convergence rate and can save both CPU time and memory cost compared to a classical upwind flux-based DD-DG (Domain Decomposed-Discontinuous Galerkin) approach.
NASA Astrophysics Data System (ADS)
Gerke, Kirill; Vasilyev, Roman; Khirevich, Siarhei; Karsanina, Marina; Collins, Daniel; Korost, Dmitry; Mallants, Dirk
2015-04-01
In this contribution we introduce a novel free software which solves the Stokes equation to obtain velocity fields for low Reynolds-number flows within externally generated 3D pore geometries. Provided with velocity fields, one can calculate permeability for known pressure gradient boundary conditions via Darcy's equation. Finite-difference schemes of 2nd and 4th order of accuracy are used together with an artificial compressibility method to iteratively converge to a steady-state solution of Stokes' equation. This numerical approach is much faster and less computationally demanding than the majority of open-source or commercial softwares employing other algorithms (finite elements/volumes, lattice Boltzmann, etc.) The software consists of two parts: 1) a pre and post-processing graphical interface, and 2) a solver. The latter is efficiently parallelized to use any number of available cores (the speedup on 16 threads was up to 10-12 depending on hardware). Due to parallelization and memory optimization our software can be used to obtain solutions for 300x300x300 voxels geometries on modern desktop PCs. The software was successfully verified by testing it against lattice Boltzmann simulations and analytical solutions. To illustrate the software's applicability for numerous problems in Earth Sciences, a number of case studies have been developed: 1) identifying the representative elementary volume for permeability determination within a sandstone sample, 2) derivation of permeability/hydraulic conductivity values for rock and soil samples and comparing those with experimentally obtained values, 3) revealing the influence of the amount of fine-textured material such as clay on filtration properties of sandy soil. This work was partially supported by RSF grant 14-17-00658 (pore-scale modelling) and RFBR grants 13-04-00409-a and 13-05-01176-a.
2-D and 3-D computations of curved accelerator magnets
Turner, L.R.
1991-01-01
In order to save computer memory, a long accelerator magnet may be computed by treating the long central region and the end regions separately. The dipole magnets for the injector synchrotron of the Advanced Photon Source (APS), now under construction at Argonne National Laboratory (ANL), employ magnet iron consisting of parallel laminations, stacked with a uniform radius of curvature of 33.379 m. Laplace's equation for the magnetic scalar potential has a different form for a straight magnet (x-y coordinates), a magnet with surfaces curved about a common center (r-{theta} coordinates), and a magnet with parallel laminations like the APS injector dipole. Yet pseudo 2-D computations for the three geometries give basically identical results, even for a much more strongly curved magnet. Hence 2-D (x-y) computations of the central region and 3-D computations of the end regions can be combined to determine the overall magnetic behavior of the magnets. 1 ref., 6 figs.
Conformal Laplace superintegrable systems in 2D: polynomial invariant subspaces
NASA Astrophysics Data System (ADS)
Escobar-Ruiz, M. A.; Miller, Willard, Jr.
2016-07-01
2nd-order conformal superintegrable systems in n dimensions are Laplace equations on a manifold with an added scalar potential and 2n-1 independent 2nd order conformal symmetry operators. They encode all the information about Helmholtz (eigenvalue) superintegrable systems in an efficient manner: there is a 1-1 correspondence between Laplace superintegrable systems and Stäckel equivalence classes of Helmholtz superintegrable systems. In this paper we focus on superintegrable systems in two-dimensions, n = 2, where there are 44 Helmholtz systems, corresponding to 12 Laplace systems. For each Laplace equation we determine the possible two-variate polynomial subspaces that are invariant under the action of the Laplace operator, thus leading to families of polynomial eigenfunctions. We also study the behavior of the polynomial invariant subspaces under a Stäckel transform. The principal new results are the details of the polynomial variables and the conditions on parameters of the potential corresponding to polynomial solutions. The hidden gl 3-algebraic structure is exhibited for the exact and quasi-exact systems. For physically meaningful solutions, the orthogonality properties and normalizability of the polynomials are presented as well. Finally, for all Helmholtz superintegrable solvable systems we give a unified construction of one-dimensional (1D) and two-dimensional (2D) quasi-exactly solvable potentials possessing polynomial solutions, and a construction of new 2D PT-symmetric potentials is established.
Ghosh, Aryya; Vaval, Nayana; Pal, Sourav
2015-07-14
Auger decay is an efficient ultrafast relaxation process of core-shell or inner-shell excited atom or molecule. Generally, it occurs in femto-second or even atto-second time domain. Direct measurement of lifetimes of Auger process of single ionized and double ionized inner-shell state of an atom or molecule is an extremely difficult task. In this paper, we have applied the highly correlated complex absorbing potential-equation-of-motion coupled cluster (CAP-EOMCC) approach which is a combination of CAP and EOMCC approach to calculate the lifetime of the states arising from 2p inner-shell ionization of an Ar atom and 3d inner-shell ionization of Kr atom. We have also calculated the lifetime of Ar{sup 2+}(2p{sup −1}3p{sup −1}) {sup 1}D, Ar{sup 2+}(2p{sup −1}3p{sup −1}) {sup 1}S, and Ar{sup 2+}(2p{sup −1}3s{sup −1}) {sup 1}P double ionized states. The predicted results are compared with the other theoretical results as well as experimental results available in the literature.
The Laplace Planes of Uranus and Pluto
NASA Technical Reports Server (NTRS)
Dobrovolskis, Anthony R.
1993-01-01
Satellite orbits close to an oblate planet precess about its equatorial plane, while distant satellites precess around the plane of the planet's heliocentric orbit. In between, satellites in nearly circular orbits precess about a warped intermediate surface called the Laplace 'plane.' Herein we derive general formulas for locating the Laplace plane. Because Uranus and Pluto have high obliquities, their Laplace planes are severely warped. We present maps of these Laplace planes, of interest in telescopic searches for new satellites. The Laplace plane of the Solar System as a whole is similarly distorted, but comets in the inner Oort cloud precess too slowly to sense the Laplace plane.
An extension of the Laplace transform to Schwartz distributions
NASA Technical Reports Server (NTRS)
Price, D. R.
1974-01-01
A characterization of the Laplace transform is developed which extends the transform to the Schwartz distributions. The class of distributions includes the impulse functions and other singular functions which occur as solutions to ordinary and partial differential equations. The standard theorems on analyticity, uniqueness, and invertibility of the transform are proved by using the characterization as the definition of the Laplace transform. The definition uses sequences of linear transformations on the space of distributions which extends the Laplace transform to another class of generalized functions, the Mikusinski operators. It is shown that the sequential definition of the transform is equivalent to Schwartz' extension of the ordinary Laplace transform to distributions but, in contrast to Schwartz' definition, does not use the distributional Fourier transform. Several theorems concerning the particular linear transformations used to define the Laplace transforms are proved. All the results proved in one dimension are extended to the n-dimensional case, but proofs are presented only for those situations that require methods different from their one-dimensional analogs.
The PROSAIC Laplace and Fourier Transform
Smith, G.A.
1994-11-01
Integral Transform methods play an extremely important role in many branches of science and engineering. The ease with which many problems may be solved using these techniques is well known. In Electrical Engineering especially, Laplace and Fourier Transforms have been used for a long time as a way to change the solution of differential equations into trivial algebraic manipulations or to provide alternate representations of signals and data. These techniques, while seemingly overshadowed by today`s emphasis on digital analysis, still form an invaluable basis in the understanding of systems and circuits. A firm grasp of the practical aspects of these subjects provides valuable conceptual tools. This tutorial paper is a review of Laplace and Fourier Transforms from an applied perspective with an emphasis on engineering applications. The interrelationship of the time and frequency domains will be stressed, in an attempt to comfort those who, after living so much of their lives in the time domain, find thinking in the frequency domain disquieting.
Electroosmotic flow through a microparallel channel with 3D wall roughness.
Chang, Long; Jian, Yongjun; Buren, Mandula; Sun, Yanjun
2016-02-01
In this paper, a perturbation method is introduced to study the EOF in a microparallel channel with 3D wall roughness. The corrugations of the two walls are periodic sinusoidal waves of small amplitude in two directions either in phase or half-period out of phase. Based on linearized Poisson-Boltzmann equation, Laplace equation, and the Navier-Stokes equations, the perturbation solutions of velocity, electrical potential, and volume flow rate are obtained. By using numerical computation, the influences of the wall corrugations on the mean velocity are analyzed. The variations of electrical potential, velocity profile, mean velocity, and their dependences on the wave number α and β of wall corrugations in two directions, the nondimensional electrokinetic width K, and the zeta potential ratio between the lower wall and the upper wall ς are analyzed graphically. PMID:26333852
Four Poission-Laplace Theory of Gravitation (I)
NASA Astrophysics Data System (ADS)
Nyambuya, Golden Gadzirayi
2015-08-01
The Poisson-Laplace equation is a working and acceptable equation of gravitation which is mostly used or applied in its differential form in Magneto-Hydro-Dynamic (MHD) modelling of e.g. molecular clouds. From a general relativistic standpoint, it describes gravitational fields in the region of low spacetime curvature as it emerges in the weak field limit. For non-static gravitational fields, this equation is not generally covariant. On the requirements of general covariance, this equation can be extended to include a time-dependent component, in which case one is led to the Four Poisson-Laplace equation. We solve the Four Poisson-Laplace equation for radial solutions, and apart from the Newtonian gravitational component, we obtain four new solutions leading to four new gravitational components capable (in-principle) of explaining e.g. the Pioneer anomaly, the Titius-Bode Law and the formation of planetary rings. In this letter, we focus only on writing down these solutions. The task showing that these new solutions might explain the aforesaid gravitational anomalies has been left for separate future readings.
NASA Astrophysics Data System (ADS)
Wang, F.; Jordan, T. H.
2012-12-01
Seismic hazard models based on empirical ground motion prediction equations (GMPEs) employ a model-based factorization to account for source, propagation, and path effects. An alternative is to simulate these effects directly using earthquake source models combined with three-dimensional (3D) models of Earth structure. We have developed an averaging-based factorization (ABF) scheme that facilitates the geographically explicit comparison of these two types of seismic hazard models. For any fault source k with epicentral position x, slip spatial and temporal distribution f, and moment magnitude m, we calculate the excitation functions G(s, k, x, m, f) for sites s in a geographical region R, such as 5% damped spectral acceleration at a particular period. Through a sequence of weighted-averaging and normalization operations following a certain hierarchy over f, m, x, k, and s, we uniquely factorize G(s, k, x, m, f) into six components: A, B(s), C(s, k), D(s, k, x), E(s, k, x, m), and F(s, k, x, m, f). Factors for a target model can be divided by those of a reference model to obtain six corresponding factor ratios, or residual factors: a, b(s), c(s, k), d(s, k, x), e(s, k, x, m), and f(s, k, x, m, f). We show that these residual factors characterize differences in basin effects primarily through b(s), distance scaling primarily through c(s, k), and source directivity primarily through d(s, k, x). We illustrate the ABF scheme by comparing CyberShake Hazard Model (CSHM) for the Los Angeles region (Graves et. al. 2010) with the Next Generation Attenuation (NGA) GMPEs modified according to the directivity relations of Spudich and Chiou (2008). Relative to CSHM, all NGA models underestimate the directivity and basin effects. In particular, the NGA models do not account for the coupling between source directivity and basin excitation that substantially enhance the low-frequency seismic hazards in the sedimentary basins of the Los Angeles region. Assuming Cyber
NASA Astrophysics Data System (ADS)
Pletinckx, D.
2011-09-01
The current 3D hype creates a lot of interest in 3D. People go to 3D movies, but are we ready to use 3D in our homes, in our offices, in our communication? Are we ready to deliver real 3D to a general public and use interactive 3D in a meaningful way to enjoy, learn, communicate? The CARARE project is realising this for the moment in the domain of monuments and archaeology, so that real 3D of archaeological sites and European monuments will be available to the general public by 2012. There are several aspects to this endeavour. First of all is the technical aspect of flawlessly delivering 3D content over all platforms and operating systems, without installing software. We have currently a working solution in PDF, but HTML5 will probably be the future. Secondly, there is still little knowledge on how to create 3D learning objects, 3D tourist information or 3D scholarly communication. We are still in a prototype phase when it comes to integrate 3D objects in physical or virtual museums. Nevertheless, Europeana has a tremendous potential as a multi-facetted virtual museum. Finally, 3D has a large potential to act as a hub of information, linking to related 2D imagery, texts, video, sound. We describe how to create such rich, explorable 3D objects that can be used intuitively by the generic Europeana user and what metadata is needed to support the semantic linking.
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITH TURB3D)
NASA Technical Reports Server (NTRS)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITHOUT TURB3D)
NASA Technical Reports Server (NTRS)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
Inversion and approximation of Laplace transforms
NASA Technical Reports Server (NTRS)
Lear, W. M.
1980-01-01
A method of inverting Laplace transforms by using a set of orthonormal functions is reported. As a byproduct of the inversion, approximation of complicated Laplace transforms by a transform with a series of simple poles along the left half plane real axis is shown. The inversion and approximation process is simple enough to be put on a programmable hand calculator.
NASA Astrophysics Data System (ADS)
Mora, Pierric; Ducasse, Eric; Deschamps, Marc
We solve the problem of the interaction of a transient guided elastic wave by a planar crack with an indirect Boundary Element approach in the Laplace domain (t → s). The originality of this work is to use the numerical Green function of the layered plate rather than the analytical Green function of each layer. As a consequence, the BEM matrices are small. To obtain the Green function in the (x,z,s) domain we first solve the equations in the Fourier transform (k,z,s) domain with a Global Matrix approach, and then perform a numerical inverse FFT. Comparisons with finite element show excellent agreement. This approach is fast and low memory consuming for planar defects in arbitrary layered media, and can be extended to arbitrary shapes and boundary conditions for a higher computational cost. It is valid in 3D, however only the 2D case is considered in this work.
NASA Astrophysics Data System (ADS)
Ren, Dandan; Ou, Yaobin
2016-08-01
In this paper, we prove the incompressible limit of all-time strong solutions to the three-dimensional full compressible Navier-Stokes equations. Here the velocity field and temperature satisfy the Dirichlet boundary condition and convective boundary condition, respectively. The uniform estimates in both the Mach number {ɛin(0,overline{ɛ}]} and time {tin[0,∞)} are established by deriving a differential inequality with decay property, where {overline{ɛ} in(0,1]} is a constant. Based on these uniform estimates, the global solution of full compressible Navier-Stokes equations with "well-prepared" initial conditions converges to the one of isentropic incompressible Navier-Stokes equations as the Mach number goes to zero.
NASA Technical Reports Server (NTRS)
Shareef, N. H.; Amirouche, F. M. L.
1991-01-01
A computational algorithmic procedure is developed and implemented for the dynamic analysis of a multibody system with rigid/flexible interconnected bodies. The algorithm takes into consideration the large rotation/translation and small elastic deformations associated with the rigid-body degrees of freedom and the flexibility of the bodies in the system respectively. Versatile three-dimensional isoparametric brick elements are employed for the modeling of the geometric configurations of the bodies. The formulation of the recursive dynamical equations of motion is based on the recursive Kane's equations, strain energy concepts, and the techniques of component mode synthesis. In order to minimize CPU-intensive matrix multiplication operations and speed up the execution process, the concepts of indexed arrays is utilized in the formulation of the equations of motion. A spin-up maneuver of a space robot with three flexible links carrying a solar panel is used as an illustrative example.
3d-3d correspondence revisited
NASA Astrophysics Data System (ADS)
Chung, Hee-Joong; Dimofte, Tudor; Gukov, Sergei; Sułkowski, Piotr
2016-04-01
In fivebrane compactifications on 3-manifolds, we point out the importance of all flat connections in the proper definition of the effective 3d {N}=2 theory. The Lagrangians of some theories with the desired properties can be constructed with the help of homological knot invariants that categorify colored Jones polynomials. Higgsing the full 3d theories constructed this way recovers theories found previously by Dimofte-Gaiotto-Gukov. We also consider the cutting and gluing of 3-manifolds along smooth boundaries and the role played by all flat connections in this operation.
NASA Astrophysics Data System (ADS)
Zheng, Xiang; Yang, Chao; Cai, Xiao-Chuan; Keyes, David
2015-03-01
We present a numerical algorithm for simulating the spinodal decomposition described by the three dimensional Cahn-Hilliard-Cook (CHC) equation, which is a fourth-order stochastic partial differential equation with a noise term. The equation is discretized in space and time based on a fully implicit, cell-centered finite difference scheme, with an adaptive time-stepping strategy designed to accelerate the progress to equilibrium. At each time step, a parallel Newton-Krylov-Schwarz algorithm is used to solve the nonlinear system. We discuss various numerical and computational challenges associated with the method. The numerical scheme is validated by a comparison with an explicit scheme of high accuracy (and unreasonably high cost). We present steady state solutions of the CHC equation in two and three dimensions. The effect of the thermal fluctuation on the spinodal decomposition process is studied. We show that the existence of the thermal fluctuation accelerates the spinodal decomposition process and that the final steady morphology is sensitive to the stochastic noise. We also show the evolution of the energies and statistical moments. In terms of the parallel performance, it is found that the implicit domain decomposition approach scales well on supercomputers with a large number of processors.
Zheng, Xiang; Yang, Chao; Cai, Xiao-Chuan; Keyes, David
2015-03-15
We present a numerical algorithm for simulating the spinodal decomposition described by the three dimensional Cahn–Hilliard–Cook (CHC) equation, which is a fourth-order stochastic partial differential equation with a noise term. The equation is discretized in space and time based on a fully implicit, cell-centered finite difference scheme, with an adaptive time-stepping strategy designed to accelerate the progress to equilibrium. At each time step, a parallel Newton–Krylov–Schwarz algorithm is used to solve the nonlinear system. We discuss various numerical and computational challenges associated with the method. The numerical scheme is validated by a comparison with an explicit scheme of high accuracy (and unreasonably high cost). We present steady state solutions of the CHC equation in two and three dimensions. The effect of the thermal fluctuation on the spinodal decomposition process is studied. We show that the existence of the thermal fluctuation accelerates the spinodal decomposition process and that the final steady morphology is sensitive to the stochastic noise. We also show the evolution of the energies and statistical moments. In terms of the parallel performance, it is found that the implicit domain decomposition approach scales well on supercomputers with a large number of processors.
Ge, Liang; Sotiropoulos, Fotis
2008-01-01
A novel numerical method is developed that integrates boundary-conforming grids with a sharp interface, immersed boundary methodology. The method is intended for simulating internal flows containing complex, moving immersed boundaries such as those encountered in several cardiovascular applications. The background domain (e.g the empty aorta) is discretized efficiently with a curvilinear boundary-fitted mesh while the complex moving immersed boundary (say a prosthetic heart valve) is treated with the sharp-interface, hybrid Cartesian/immersed-boundary approach of Gilmanov and Sotiropoulos [1]. To facilitate the implementation of this novel modeling paradigm in complex flow simulations, an accurate and efficient numerical method is developed for solving the unsteady, incompressible Navier-Stokes equations in generalized curvilinear coordinates. The method employs a novel, fully-curvilinear staggered grid discretization approach, which does not require either the explicit evaluation of the Christoffel symbols or the discretization of all three momentum equations at cell interfaces as done in previous formulations. The equations are integrated in time using an efficient, second-order accurate fractional step methodology coupled with a Jacobian-free, Newton-Krylov solver for the momentum equations and a GMRES solver enhanced with multigrid as preconditioner for the Poisson equation. Several numerical experiments are carried out on fine computational meshes to demonstrate the accuracy and efficiency of the proposed method for standard benchmark problems as well as for unsteady, pulsatile flow through a curved, pipe bend. To demonstrate the ability of the method to simulate flows with complex, moving immersed boundaries we apply it to calculate pulsatile, physiological flow through a mechanical, bileaflet heart valve mounted in a model straight aorta with an anatomical-like triple sinus. PMID:19194533
Ge, Liang; Sotiropoulos, Fotis
2007-08-01
A novel numerical method is developed that integrates boundary-conforming grids with a sharp interface, immersed boundary methodology. The method is intended for simulating internal flows containing complex, moving immersed boundaries such as those encountered in several cardiovascular applications. The background domain (e.g the empty aorta) is discretized efficiently with a curvilinear boundary-fitted mesh while the complex moving immersed boundary (say a prosthetic heart valve) is treated with the sharp-interface, hybrid Cartesian/immersed-boundary approach of Gilmanov and Sotiropoulos [1]. To facilitate the implementation of this novel modeling paradigm in complex flow simulations, an accurate and efficient numerical method is developed for solving the unsteady, incompressible Navier-Stokes equations in generalized curvilinear coordinates. The method employs a novel, fully-curvilinear staggered grid discretization approach, which does not require either the explicit evaluation of the Christoffel symbols or the discretization of all three momentum equations at cell interfaces as done in previous formulations. The equations are integrated in time using an efficient, second-order accurate fractional step methodology coupled with a Jacobian-free, Newton-Krylov solver for the momentum equations and a GMRES solver enhanced with multigrid as preconditioner for the Poisson equation. Several numerical experiments are carried out on fine computational meshes to demonstrate the accuracy and efficiency of the proposed method for standard benchmark problems as well as for unsteady, pulsatile flow through a curved, pipe bend. To demonstrate the ability of the method to simulate flows with complex, moving immersed boundaries we apply it to calculate pulsatile, physiological flow through a mechanical, bileaflet heart valve mounted in a model straight aorta with an anatomical-like triple sinus. PMID:19194533
NASA Astrophysics Data System (ADS)
Meulien Ohlmann, Odile
2013-02-01
Today the industry offers a chain of 3D products. Learning to "read" and to "create in 3D" becomes an issue of education of primary importance. 25 years professional experience in France, the United States and Germany, Odile Meulien set up a personal method of initiation to 3D creation that entails the spatial/temporal experience of the holographic visual. She will present some different tools and techniques used for this learning, their advantages and disadvantages, programs and issues of educational policies, constraints and expectations related to the development of new techniques for 3D imaging. Although the creation of display holograms is very much reduced compared to the creation of the 90ies, the holographic concept is spreading in all scientific, social, and artistic activities of our present time. She will also raise many questions: What means 3D? Is it communication? Is it perception? How the seeing and none seeing is interferes? What else has to be taken in consideration to communicate in 3D? How to handle the non visible relations of moving objects with subjects? Does this transform our model of exchange with others? What kind of interaction this has with our everyday life? Then come more practical questions: How to learn creating 3D visualization, to learn 3D grammar, 3D language, 3D thinking? What for? At what level? In which matter? for whom?
NASA Astrophysics Data System (ADS)
Bruno, Oscar P.; Cubillos, Max
2016-02-01
This paper introduces alternating-direction implicit (ADI) solvers of higher order of time-accuracy (orders two to six) for the compressible Navier-Stokes equations in two- and three-dimensional curvilinear domains. The higher-order accuracy in time results from 1) An application of the backward differentiation formulae time-stepping algorithm (BDF) in conjunction with 2) A BDF-like extrapolation technique for certain components of the nonlinear terms (which makes use of nonlinear solves unnecessary), as well as 3) A novel application of the Douglas-Gunn splitting (which greatly facilitates handling of boundary conditions while preserving higher-order accuracy in time). As suggested by our theoretical analysis of the algorithms for a variety of special cases, an extensive set of numerical experiments clearly indicate that all of the BDF-based ADI algorithms proposed in this paper are "quasi-unconditionally stable" in the following sense: each algorithm is stable for all couples (h , Δt)of spatial and temporal mesh sizes in a problem-dependent rectangular neighborhood of the form (0 ,Mh) × (0 ,Mt). In other words, for each fixed value of Δt below a certain threshold, the Navier-Stokes solvers presented in this paper are stable for arbitrarily small spatial mesh-sizes. The second-order formulation has further been rigorously shown to be unconditionally stable for linear hyperbolic and parabolic equations in two-dimensional space. Although implicit ADI solvers for the Navier-Stokes equations with nominal second-order of temporal accuracy have been proposed in the past, the algorithms presented in this paper are the first ADI-based Navier-Stokes solvers for which second-order or better accuracy has been verified in practice under non-trivial (non-periodic) boundary conditions.
Kaltenbacher, Barbara; Kaltenbacher, Manfred; Sim, Imbo
2013-01-01
We consider the second order wave equation in an unbounded domain and propose an advanced perfectly matched layer (PML) technique for its efficient and reliable simulation. In doing so, we concentrate on the time domain case and use the finite-element (FE) method for the space discretization. Our un-split-PML formulation requires four auxiliary variables within the PML region in three space dimensions. For a reduced version (rPML), we present a long time stability proof based on an energy analysis. The numerical case studies and an application example demonstrate the good performance and long time stability of our formulation for treating open domain problems. PMID:23888085
Kaltenbacher, Barbara; Kaltenbacher, Manfred; Sim, Imbo
2013-02-15
We consider the second order wave equation in an unbounded domain and propose an advanced perfectly matched layer (PML) technique for its efficient and reliable simulation. In doing so, we concentrate on the time domain case and use the finite-element (FE) method for the space discretization. Our un-split-PML formulation requires four auxiliary variables within the PML region in three space dimensions. For a reduced version (rPML), we present a long time stability proof based on an energy analysis. The numerical case studies and an application example demonstrate the good performance and long time stability of our formulation for treating open domain problems. PMID:23888085
NASA Astrophysics Data System (ADS)
Kaltenbacher, Barbara; Kaltenbacher, Manfred; Sim, Imbo
2013-02-01
We consider the second order wave equation in an unbounded domain and propose an advanced perfectly matched layer (PML) technique for its efficient and reliable simulation. In doing so, we concentrate on the time domain case and use the finite-element (FE) method for the space discretization. Our un-split-PML formulation requires four auxiliary variables within the PML region in three space dimensions. For a reduced version (rPML), we present a long time stability proof based on an energy analysis. The numerical case studies and an application example demonstrate the good performance and long time stability of our formulation for treating open domain problems.
Generalized Laplace Transforms and Extended Heaviside Calculus
ERIC Educational Resources Information Center
Deakin, Michael A. B.
2008-01-01
An extended Heaviside calculus proposed by Peraire in a recent paper is similar to a generalization of the Laplace transform proposed by the present author. This similarity will be illustrated by analysis of an example supplied by Peraire.
NASA Astrophysics Data System (ADS)
Jia, Xuanji; Zhou, Yong
2015-09-01
We prove that a weak solution (u, b) to the MHD equations is smooth on (0, T ] if \\text{M}\\in {{L}α}≤ft(0,T;{{L}γ}≤ft({{{R}}3}\\right)\\right) with 2/α +3/γ =2 , 1≤slant α <∞ and 3/2<γ ≤slant ∞ , where \\text{M} is a 3× 3 mixture matrix (see its definition below). As we will explain later, this kind of regularity criteria is more likely to capture the nature of the coupling effects between the fluid velocity and the magnetic field in the evolution of the MHD flows. Moreover, the condition on \\text{M} is scaling invariant, i.e. it is of Ladyzhenskaya-Prodi-Serrin type.
On computing Laplace's coefficients and their derivatives.
NASA Astrophysics Data System (ADS)
Gerasimov, I. A.; Vinnikov, E. L.
The algorithm of computing Laplace's coefficients and their derivatives is proposed with application of recurrent relations. The A.G.M.-method is used for the calculation of values L0(0), L0(1). The FORTRAN-program corresponding to the algorithm is given. The precision control was provided with numerical integrating by Simpsons method. The behavior of Laplace's coefficients and their third derivatives whith varying indices K, n for fixed values of the α-parameter is presented graphically.
ERIC Educational Resources Information Center
Hastings, S. K.
2002-01-01
Discusses 3 D imaging as it relates to digital representations in virtual library collections. Highlights include X-ray computed tomography (X-ray CT); the National Science Foundation (NSF) Digital Library Initiatives; output peripherals; image retrieval systems, including metadata; and applications of 3 D imaging for libraries and museums. (LRW)
Energy Science and Technology Software Center (ESTSC)
2008-10-05
Laplace is a electric field driven flow simulation program for detailed device design support. Transport processes include electrokinesis, dielectrophoresis, and diffusion. Laplace solves for the electric field in a microfluidic system and the liquid and particle flow that is produced by the electric field for the primary purpose of microfluidic design development and simulation. Laplace allows you to visualize the flow by tracking tracer particles, viewing flow streamlines, etc. Laplace can make movies of simulatedmore » particle motion to allow you to test and share the behavior of microfuidic designs. The electric field is calculated using an iterative linear solver and particle motion is solved by finite difference, finite-displacement simulation of particle trajectories. Laplace uses a bitmapped picture or drawing of a microsystem to infer the geometry. The channel depth is everywhere proportional to the magnitude of the blue channel of the image: 0 (black) = zero depth, or no channel, 256 (saturated blue) = deepest channel, and intermediate values correspond to intermediate depths. Laplace automatically applies various boundary conditions (applied voltage or current) to ports, where channels cross the edge of the image.« less
Cummings, Eric; LaJeunesse, Tony
2008-10-05
Laplace is a electric field driven flow simulation program for detailed device design support. Transport processes include electrokinesis, dielectrophoresis, and diffusion. Laplace solves for the electric field in a microfluidic system and the liquid and particle flow that is produced by the electric field for the primary purpose of microfluidic design development and simulation. Laplace allows you to visualize the flow by tracking tracer particles, viewing flow streamlines, etc. Laplace can make movies of simulated particle motion to allow you to test and share the behavior of microfuidic designs. The electric field is calculated using an iterative linear solver and particle motion is solved by finite difference, finite-displacement simulation of particle trajectories. Laplace uses a bitmapped picture or drawing of a microsystem to infer the geometry. The channel depth is everywhere proportional to the magnitude of the blue channel of the image: 0 (black) = zero depth, or no channel, 256 (saturated blue) = deepest channel, and intermediate values correspond to intermediate depths. Laplace automatically applies various boundary conditions (applied voltage or current) to ports, where channels cross the edge of the image.
NASA Astrophysics Data System (ADS)
Wang, Haogang; Liao, Tien-Hao; Shi, Jiancheng; Yu, Zherui
2014-11-01
The forthcoming Water Cycle Observation Mission (WCOM) is to understand the water cycle system among land, atmosphere, and ocean. In both active and passive microwave remote sensing of soil moisture, the surface roughness plays an important role. Electromagnetic models of roughness provide tables of emissivities and backscattering coefficients that can be used to retrieve soil moisture. In this paper, a fast and accurate three dimensional solution of Maxwell's equations is developed and employed to solve rough soil surface scattering problem at L-band. The algorithm combines QR Pre-Ranked Multilevel UV(MLUV) factorization and Hierarchical Fast Far Field Approximation. It is implemented using OpenMP interface for fast parallel calculation. In this algorithm, 1) QR based rank predetermined algorithm is derived to further compress the UV matrix pairs obtained using coarse-coarse sampling; 2) at the finer levels, MLUV is used straightforwardly to factorize the interactions between groups, while at the coarsest level, interactions between groups in the interaction list are calculated using an elegantly derived Hierarchical Fast Far Field Approximation (HFAFFA) to accelerate the calculation of interactions between large groups while keeping the accuracy of this approximation; 3) OpenMP interface is used to parallelize this new algorithm. Numerical results including the incoherent bistatic scattering coefficients and the emissivity demonstrate the efficiency of this method.
NASA Astrophysics Data System (ADS)
Fortes, A. Dominic; Suard, Emmanuelle; Lemée-Cailleau, Marie-Hélène; Pickard, Christopher J.; Needs, Richard J.
2009-10-01
We describe the results of a neutron powder diffraction study of perdeuterated ammonia monohydrate (AMH, ND3ṡD2O) carried out in the range 102
equation of state of AMH I has parameters, V0=248.00(2) Å3, K0=7.33(3) GPa with the first pressure derivative of K0 fixed at the value obtained in ab initio calculations, (∂K0/∂P)T=K0'=5.3; the implied value of the second derivative is therefore (∂2K0/∂P2)T=K0″=-0.94(1) GPa-1. At 351 MPa, we observed that the transition from AMH I to AMH II occurred over a period of 90 min, with an associated reduction in molar volume of 4.6% and an increase in the incompressibility of 19.6%.
Zaeytijd, J. de Bogaert, I.; Franchois, A.
2008-07-01
Electromagnetic scattering problems involving inhomogeneous objects can be numerically solved by applying a Method of Moments discretization to the volume integral equation. For electrically large problems, the iterative solution of the resulting linear system is expensive, both computationally and in memory use. In this paper, a hybrid MLFMA-FFT method is presented, which combines the fast Fourier transform (FFT) method and the High Frequency Multilevel Fast Multipole Algorithm (MLFMA) in order to reduce the cost of the matrix-vector multiplications needed in the iterative solver. The method represents the scatterers within a set of possibly disjoint identical cubic subdomains, which are meshed using a uniform cubic grid. This specific mesh allows for the application of FFTs to calculate the near interactions in the MLFMA and reduces the memory cost considerably, since the aggregation and disaggregation matrices of the MLFMA can be reused. Additional improvements to the general MLFMA framework, such as an extention of the FFT interpolation scheme of Sarvas et al. from the scalar to the vectorial case in combination with a more economical representation of the radiation patterns on the lowest level in vector spherical harmonics, are proposed and the choice of the subdomain size is discussed. The hybrid method performs better in terms of speed and memory use on large sparse configurations than both the FFT method and the HF MLFMA separately and it has lower memory requirements on general large problems. This is illustrated on a number of representative numerical test cases.
NASA Astrophysics Data System (ADS)
Ando, R.
2014-12-01
The boundary integral equation method formulated in the real space and time domain (BIEM-ST) has been used as a powerful tool to analyze the earthquake rupture dynamics on non-planar faults. Generally, BIEM is more accurate than volumetric methods such as the finite difference method and the finite difference method. With the recent development of the high performance computing environment, the earthquake rupture simulation studies have been conducted considering three dimensional realistic fault geometry models. However, the utility of BIEM-ST has been limited due to its heavy computational demanding increased depending on square of time steps (N2), which was needed to evaluate the historic integration. While BIEM can be efficient with the spectral domain formulation, the applications of such a method are limited to planar fault cases. In this study, we propose a new method to reduce the calculation time of BIEM-ST to linear of time step (N) without degrading the accuracy in the 3 dimensional modeling space. We extends the method proposed earlier for the case of the 2 dimensional framework, applying the asymptotic expressions of the elasto-dynamic Green's functions. This method uses the physical nature of the stress Green's function as dividing the causality cone according to the distances from the wave-fronts. The scalability of this method is shown on the parallel computing environment of the distributed memory. We demonstrate the applicability to analyses of subduction earthquake cases, suffering long time from the numerical limitations of previously available BIEMs. We analyze the dynamic rupture processes on dipping reverse faults embed in a three dimensional elastic half space.
Transfer Functions Via Laplace- And Fourier-Borel Transforms
NASA Technical Reports Server (NTRS)
Can, Sumer; Unal, Aynur
1991-01-01
Approach to solution of nonlinear ordinary differential equations involves transfer functions based on recently-introduced Laplace-Borel and Fourier-Borel transforms. Main theorem gives transform of response of nonlinear system as Cauchy product of transfer function and transform of input function of system, together with memory effects. Used to determine responses of electrical circuits containing variable inductances or resistances. Also possibility of doing all noncommutative algebra on computers in such symbolic programming languages as Macsyma, Reduce, PL1, or Lisp. Process of solution organized and possibly simplified by algebraic manipulations reducing integrals in solutions to known or tabulated forms.
Cao, Qian; Thawait, Gaurav; Gang, Grace J.; Zbijewski, Wojciech; Reigel, Thomas; Brown, Tyler; Corner, Brian; Demehri, Shadpour; Siewerdsen, Jeffrey H.
2015-01-01
Joint space morphology can be indicative of the risk, presence, progression, and/or treatment response of disease or trauma. We describe a novel methodology of characterizing joint space morphology in high-resolution 3D images [e.g., cone-beam CT (CBCT)] using a model based on elementary electrostatics that overcomes a variety of basic limitations of existing 2D and 3D methods. The method models each surface of a joint as a conductor at fixed electrostatic potential and characterizes the intra-articular space in terms of the electric field lines resulting from the solution of Gauss’ Law and the Laplace equation. As a test case, the method was applied to discrimination of healthy and osteoarthritic subjects (N = 39) in 3D images of the knee acquired on an extremity CBCT system. The method demonstrated improved diagnostic performance (area under the receiver operating characteristic curve, AUC > 0.98) compared to simpler methods of quantitative measurement and qualitative image-based assessment by three expert musculoskeletal radiologists (AUC = 0.87, p-value = 0.007). The method is applicable to simple (e.g., the knee or elbow) or multi-axial joints (e.g., the wrist or ankle) and may provide a useful means of quantitatively assessing a variety of joint pathologies. PMID:25575100
NASA Astrophysics Data System (ADS)
Cao, Qian; Thawait, Gaurav; Gang, Grace J.; Zbijewski, Wojciech; Reigel, Thomas; Brown, Tyler; Corner, Brian; Demehri, Shadpour; Siewerdsen, Jeffrey H.
2015-02-01
Joint space morphology can be indicative of the risk, presence, progression, and/or treatment response of disease or trauma. We describe a novel methodology of characterizing joint space morphology in high-resolution 3D images (e.g. cone-beam CT (CBCT)) using a model based on elementary electrostatics that overcomes a variety of basic limitations of existing 2D and 3D methods. The method models each surface of a joint as a conductor at fixed electrostatic potential and characterizes the intra-articular space in terms of the electric field lines resulting from the solution of Gauss’ Law and the Laplace equation. As a test case, the method was applied to discrimination of healthy and osteoarthritic subjects (N = 39) in 3D images of the knee acquired on an extremity CBCT system. The method demonstrated improved diagnostic performance (area under the receiver operating characteristic curve, AUC > 0.98) compared to simpler methods of quantitative measurement and qualitative image-based assessment by three expert musculoskeletal radiologists (AUC = 0.87, p-value = 0.007). The method is applicable to simple (e.g. the knee or elbow) or multi-axial joints (e.g. the wrist or ankle) and may provide a useful means of quantitatively assessing a variety of joint pathologies.
Crandall, K.R.
1987-08-01
TRACE 3-D is an interactive beam-dynamics program that calculates the envelopes of a bunched beam, including linear space-charge forces, through a user-defined transport system. TRACE 3-D provides an immediate graphics display of the envelopes and the phase-space ellipses and allows nine types of beam-matching options. This report describes the beam-dynamics calculations and gives detailed instruction for using the code. Several examples are described in detail.
NASA Astrophysics Data System (ADS)
Petrov, P.; Newman, G. A.
2013-12-01
Full Waveform Inversion (FWI) is a promising seismic imaging method which aims to compute quantitative estimates of the subsurface parameters (bulk wave velocity, shear wave velocity, rock density) from local measurements of the seismic wavefield. It based on a modeling of wave propagation from controlled seismic sources and consists in minimizing iteratively the difference between the predicted wavefields at the receivers and the recorded data. This amounts to solving a strongly nonlinear large-scale inverse problem. We have formulated a 3D inverse solution for the elastic wave problem in Laplace-Fourier domain using the non-linear conjugate gradient methods. FWI in the Laplace-Fourier domain has seen considerable interest over the last several years as an effective approach to imaging seismic data, where conventional FWI schemes have difficulties in converging to acceptable solutions. Finite difference methods, employing staggered grids are used to compute predicted data and objective functional gradients. Typically three forward modeling applications per frequency are required to produce the model update at each iteration. Because the solution's realism and complexity are still limited by the speed and memory of serial processors the code has been implemented on a massive parallel computing platform where hundreds to thousands of processors operate on the problem simultaneously. The inversion scheme is tested by inverting data produced with a forward modelling code for a few models.
NASA Astrophysics Data System (ADS)
Oldham, Mark
2015-01-01
Radiochromic materials exhibit a colour change when exposed to ionising radiation. Radiochromic film has been used for clinical dosimetry for many years and increasingly so recently, as films of higher sensitivities have become available. The two principle advantages of radiochromic dosimetry include greater tissue equivalence (radiologically) and the lack of requirement for development of the colour change. In a radiochromic material, the colour change arises direct from ionising interactions affecting dye molecules, without requiring any latent chemical, optical or thermal development, with important implications for increased accuracy and convenience. It is only relatively recently however, that 3D radiochromic dosimetry has become possible. In this article we review recent developments and the current state-of-the-art of 3D radiochromic dosimetry, and the potential for a more comprehensive solution for the verification of complex radiation therapy treatments, and 3D dose measurement in general.
NASA Astrophysics Data System (ADS)
Iliesiu, Luca; Kos, Filip; Poland, David; Pufu, Silviu S.; Simmons-Duffin, David; Yacoby, Ran
2016-03-01
We study the conformal bootstrap for a 4-point function of fermions < ψψψψ> in 3D. We first introduce an embedding formalism for 3D spinors and compute the conformal blocks appearing in fermion 4-point functions. Using these results, we find general bounds on the dimensions of operators appearing in the ψ × ψ OPE, and also on the central charge C T . We observe features in our bounds that coincide with scaling dimensions in the GrossNeveu models at large N . We also speculate that other features could coincide with a fermionic CFT containing no relevant scalar operators.
Clement, T.P.; Jones, N.L.
1998-02-01
RT3D (Reactive Transport in 3-Dimensions) is a computer code that solves coupled partial differential equations that describe reactive-flow and transport of multiple mobile and/or immobile species in a three dimensional saturated porous media. RT3D was developed from the single-species transport code, MT3D (DoD-1.5, 1997 version). As with MT3D, RT3D also uses the USGS groundwater flow model MODFLOW for computing spatial and temporal variations in groundwater head distribution. This report presents a set of tutorial problems that are designed to illustrate how RT3D simulations can be performed within the Department of Defense Groundwater Modeling System (GMS). GMS serves as a pre- and post-processing interface for RT3D. GMS can be used to define all the input files needed by RT3D code, and later the code can be launched from within GMS and run as a separate application. Once the RT3D simulation is completed, the solution can be imported to GMS for graphical post-processing. RT3D v1.0 supports several reaction packages that can be used for simulating different types of reactive contaminants. Each of the tutorials, described below, provides training on a different RT3D reaction package. Each reaction package has different input requirements, and the tutorials are designed to describe these differences. Furthermore, the tutorials illustrate the various options available in GMS for graphical post-processing of RT3D results. Users are strongly encouraged to complete the tutorials before attempting to use RT3D and GMS on a routine basis.
Laplace and Z transforms of linear dynamical systems and conic sections
NASA Astrophysics Data System (ADS)
Rodrigo, Marianito R.
2016-06-01
We consider the solution trajectories of linear continuous and discrete dynamical systems and show that in the Laplace and Z transform spaces, respectively, they lie on the intersections of hypersurfaces described by second-degree polynomial equations. In particular, in two dimensions, these intersections are conic sections whose types are determined by the eigenvalues of the coefficient matrices.
Laplace-Gauss and Helmholtz-Gauss paraxial modes in media with quadratic refraction index.
Kiselev, Aleksei P; Plachenov, Alexandr B
2016-04-01
The scalar theory of paraxial wave propagation in an axisymmetric medium where the refraction index quadratically depends on transverse variables is addressed. Exact solutions of the corresponding parabolic equation are presented, generalizing the Laplace-Gauss and Helmholtz-Gauss modes earlier known for homogeneous media. Also, a generalization of a zero-order asymmetric Bessel-Gauss beam is given. PMID:27140777
Software for 3D radiotherapy dosimetry. Validation
NASA Astrophysics Data System (ADS)
Kozicki, Marek; Maras, Piotr; Karwowski, Andrzej C.
2014-08-01
The subject of this work is polyGeVero® software (GeVero Co., Poland), which has been developed to fill the requirements of fast calculations of 3D dosimetry data with the emphasis on polymer gel dosimetry for radiotherapy. This software comprises four workspaces that have been prepared for: (i) calculating calibration curves and calibration equations, (ii) storing the calibration characteristics of the 3D dosimeters, (iii) calculating 3D dose distributions in irradiated 3D dosimeters, and (iv) comparing 3D dose distributions obtained from measurements with the aid of 3D dosimeters and calculated with the aid of treatment planning systems (TPSs). The main features and functions of the software are described in this work. Moreover, the core algorithms were validated and the results are presented. The validation was performed using the data of the new PABIGnx polymer gel dosimeter. The polyGeVero® software simplifies and greatly accelerates the calculations of raw 3D dosimetry data. It is an effective tool for fast verification of TPS-generated plans for tumor irradiation when combined with a 3D dosimeter. Consequently, the software may facilitate calculations by the 3D dosimetry community. In this work, the calibration characteristics of the PABIGnx obtained through four calibration methods: multi vial, cross beam, depth dose, and brachytherapy, are discussed as well.
Laplace, Pierre-Simon (1749-1827)
NASA Astrophysics Data System (ADS)
Murdin, P.
2000-11-01
Celestial mechanician, born in Beaumont-en-Auge, Normandy, France, became professor of mathematics at the Ecole Militaire in Paris, examining the cadet Napoleon Bonaparte. This position made Laplace well known to people in positions of power, which he opportunistically exploited, becoming, under Napoleon, Minister of the Interior (Napoleon soon removed him from office `because he brought the spir...
An approximation for inverse Laplace transforms
NASA Technical Reports Server (NTRS)
Lear, W. M.
1981-01-01
Programmable calculator runs simple finite-series approximation for Laplace transform inversions. Utilizing family of orthonormal functions, approximation is used for wide range of transforms, including those encountered in feedback control problems. Method works well as long as F(t) decays to zero as it approaches infinity and so is appliable to most physical systems.
Evaluation of the Laplace Integral. Classroom Notes
ERIC Educational Resources Information Center
Chen, Hongwei
2004-01-01
Based on the dominated convergence theorem and parametric differentiation, two different evaluations of the Laplace integral are displayed. This article presents two different proofs of (1) which may be of interest since they are based on principles within the realm of real analysis. The first method applies the dominated convergence theorem to…
A Note on Laplace's Expansion Theorem
ERIC Educational Resources Information Center
Janji, Milan
2005-01-01
A short proof of Laplace's expansion theorem is given. The proof is elementary and can be presented at any level of undergraduate studies where determinants are taught. It is derived directly from the definition so that the theorem may be used as a starting point for further investigation of determinants.
Parseval-Type Relations for Laplace Transform and their Applications
ERIC Educational Resources Information Center
Herman, S.; Maceli, J.; Rogala, M.; Yurekli, O.
2008-01-01
In the present note, two Parseval-type relations involving the Laplace transform are given. The application of the relations is demonstrated in evaluating improper integrals and Laplace transforms of trigonometric functions.
Hong, X; Gao, H
2014-06-15
Purpose: The Linear Boltzmann Transport Equation (LBTE) solved through statistical Monte Carlo (MC) method provides the accurate dose calculation in radiotherapy. This work is to investigate the alternative way for accurately solving LBTE using deterministic numerical method due to its possible advantage in computational speed from MC. Methods: Instead of using traditional spherical harmonics to approximate angular scattering kernel, our deterministic numerical method directly computes angular scattering weights, based on a new angular discretization method that utilizes linear finite element method on the local triangulation of unit angular sphere. As a Result, our angular discretization method has the unique advantage in positivity, i.e., to maintain all scattering weights nonnegative all the time, which is physically correct. Moreover, our method is local in angular space, and therefore handles the anisotropic scattering well, such as the forward-peaking scattering. To be compatible with image-guided radiotherapy, the spatial variables are discretized on the structured grid with the standard diamond scheme. After discretization, the improved sourceiteration method is utilized for solving the linear system without saving the linear system to memory. The accuracy of our 3D solver is validated using analytic solutions and benchmarked with Geant4, a popular MC solver. Results: The differences between Geant4 solutions and our solutions were less than 1.5% for various testing cases that mimic the practical cases. More details are available in the supporting document. Conclusion: We have developed a 3D LBTE solver based on a new angular discretization method that guarantees the positivity of scattering weights for physical correctness, and it has been benchmarked with Geant4 for photon dose calculation.
NASA Astrophysics Data System (ADS)
Iizuka, Keigo
2008-02-01
In order to circumvent the fact that only one observer can view the image from a stereoscopic microscope, an attachment was devised for displaying the 3D microscopic image on a large LCD monitor for viewing by multiple observers in real time. The principle of operation, design, fabrication, and performance are presented, along with tolerance measurements relating to the properties of the cellophane half-wave plate used in the design.
Tidal friction in the Earth-Moon system and Laplace planes: Darwin redux
NASA Astrophysics Data System (ADS)
Rubincam, David Parry
2016-03-01
The dynamical evolution of the Earth-Moon system due to tidal friction is treated here. George H. Darwin used Laplace planes (also called proper planes) in his study of tidal evolution. The Laplace plane approach is adapted here to the formalisms of W.M. Kaula and P. Goldreich. Like Darwin, the approach assumes a three-body problem: Earth, Moon, and Sun, where the Moon and Sun are point-masses. The tidal potential is written in terms of the Laplace plane angles. The resulting secular equations of motion can be easily integrated numerically assuming the Moon is in a circular orbit about the Earth and the Earth is in a circular orbit about the Sun. For Earth-Moon distances greater than ∼10 Earth radii, the Earth's approximate tidal response can be characterized with a single parameter, which is a ratio: a Love number times the sine of a lag angle divided by another such product. For low parameter values it can be shown that Darwin's low-viscosity molten Earth, M. Ross's and G. Schubert's model of an Earth near melting, and Goldreich's equal tidal lag angles must all give similar histories. For higher parameter values, as perhaps has been the case at times with the ocean tides, the Earth's obliquity may have decreased slightly instead of increased once the Moon's orbit evolved further than 50 Earth radii from the Earth, with possible implications for climate. This is contrast to the other tidal friction models mentioned, which have the obliquity always increasing with time. As for the Moon, its orbit is presently tilted to its Laplace plane by 5.2°. The equations do not allow the Moon to evolve out of its Laplace plane by tidal friction alone, so that if it was originally in its Laplace plane, the tilt arose with the addition of other mechanisms, such as resonance passages.
Tidal Friction in the Earth-Moon System and Laplace Planes: Darwin Redux
NASA Technical Reports Server (NTRS)
Rubincam, David P.
2015-01-01
The dynamical evolution of the Earth-Moon system due to tidal friction is treated here. George H. Darwin used Laplace planes (also called proper planes) in his study of tidal evolution. The Laplace plane approach is adapted here to the formalisms of W.M. Kaula and P. Goldreich. Like Darwin, the approach assumes a three-body problem: Earth, Moon, and Sun, where the Moon and Sun are point-masses. The tidal potential is written in terms of the Laplace plane angles. The resulting secular equations of motion can be easily integrated numerically assuming the Moon is in a circular orbit about the Earth and the Earth is in a circular orbit about the Sun. For Earth-Moon distances greater than 10 Earth radii, the Earth's approximate tidal response can be characterized with a single parameter, which is a ratio: a Love number times the sine of a lag angle divided by another such product. For low parameter values it can be shown that Darwin's low-viscosity molten Earth, M. Ross's and G. Schubert's model of an Earth near melting, and Goldreich's equal tidal lag angles must all give similar histories. For higher parameter values, as perhaps has been the case at times with the ocean tides, the Earth's obliquity may have decreased slightly instead of increased once the Moon's orbit evolved further than 50 Earth radii from the Earth, with possible implications for climate. This is contrast to the other tidal friction models mentioned, which have the obliquity always increasing with time. As for the Moon, its orbit is presently tilted to its Laplace plane by 5.2deg. The equations do not allow the Moon to evolve out of its Laplace plane by tidal friction alone, so that if it was originally in its Laplace plane, the tilt arose with the addition of other mechanisms, such as resonance passages.
Modelling a single phase voltage controlled rectifier using Laplace transforms
NASA Technical Reports Server (NTRS)
Kraft, L. Alan; Kankam, M. David
1992-01-01
The development of a 20 kHz, AC power system by NASA for large space projects has spurred a need to develop models for the equipment which will be used on these single phase systems. To date, models for the AC source (i.e., inverters) have been developed. It is the intent of this paper to develop a method to model the single phase voltage controlled rectifiers which will be attached to the AC power grid as an interface for connected loads. A modified version of EPRI's HARMFLO program is used as the shell for these models. The results obtained from the model developed in this paper are quite adequate for the analysis of problems such as voltage resonance. The unique technique presented in this paper uses the Laplace transforms to determine the harmonic content of the load current of the rectifier rather than a curve fitting technique. Laplace transforms yield the coefficient of the differential equations which model the line current to the rectifier directly.
NASA Astrophysics Data System (ADS)
Kostrzewski, Andrew A.; Aye, Tin M.; Kim, Dai Hyun; Esterkin, Vladimir; Savant, Gajendra D.
1998-09-01
Physical Optics Corporation has developed an advanced 3-D virtual reality system for use with simulation tools for training technical and military personnel. This system avoids such drawbacks of other virtual reality (VR) systems as eye fatigue, headaches, and alignment for each viewer, all of which are due to the need to wear special VR goggles. The new system is based on direct viewing of an interactive environment. This innovative holographic multiplexed screen technology makes it unnecessary for the viewer to wear special goggles.
NASA Technical Reports Server (NTRS)
1992-01-01
Ames Research Center research into virtual reality led to the development of the Convolvotron, a high speed digital audio processing system that delivers three-dimensional sound over headphones. It consists of a two-card set designed for use with a personal computer. The Convolvotron's primary application is presentation of 3D audio signals over headphones. Four independent sound sources are filtered with large time-varying filters that compensate for motion. The perceived location of the sound remains constant. Possible applications are in air traffic control towers or airplane cockpits, hearing and perception research and virtual reality development.
Spectral Laplace-Beltrami wavelets with applications in medical images.
Tan, Mingzhen; Qiu, Anqi
2015-05-01
The spectral graph wavelet transform (SGWT) has recently been developed to compute wavelet transforms of functions defined on non-Euclidean spaces such as graphs. By capitalizing on the established framework of the SGWT, we adopt a fast and efficient computation of a discretized Laplace-Beltrami (LB) operator that allows its extension from arbitrary graphs to differentiable and closed 2-D manifolds (smooth surfaces embedded in the 3-D Euclidean space). This particular class of manifolds are widely used in bioimaging to characterize the morphology of cells, tissues, and organs. They are often discretized into triangular meshes, providing additional geometric information apart from simple nodes and weighted connections in graphs. In comparison with the SGWT, the wavelet bases constructed with the LB operator are spatially localized with a more uniform "spread" with respect to underlying curvature of the surface. In our experiments, we first use synthetic data to show that traditional applications of wavelets in smoothing and edge detectio can be done using the wavelet bases constructed with the LB operator. Second, we show that multi-resolutional capabilities of the proposed framework are applicable in the classification of Alzheimer's patients with normal subjects using hippocampal shapes. Wavelet transforms of the hippocampal shape deformations at finer resolutions registered higher sensitivity (96%) and specificity (90%) than the classification results obtained from the direct usage of hippocampal shape deformations. In addition, the Laplace-Beltrami method requires consistently a smaller number of principal components (to retain a fixed variance) at higher resolution as compared to the binary and weighted graph Laplacians, demonstrating the potential of the wavelet bases in adapting to the geometry of the underlying manifold. PMID:25343758
Cevidanes, Lucia; Tucker, Scott; Styner, Martin; Kim, Hyungmin; Chapuis, Jonas; Reyes, Mauricio; Proffit, William; Turvey, Timothy; Jaskolka, Michael
2009-01-01
This paper discusses the development of methods for computer-aided jaw surgery. Computer-aided jaw surgery allows us to incorporate the high level of precision necessary for transferring virtual plans into the operating room. We also present a complete computer-aided surgery (CAS) system developed in close collaboration with surgeons. Surgery planning and simulation include construction of 3D surface models from Cone-beam CT (CBCT), dynamic cephalometry, semi-automatic mirroring, interactive cutting of bone and bony segment repositioning. A virtual setup can be used to manufacture positioning splints for intra-operative guidance. The system provides further intra-operative assistance with the help of a computer display showing jaw positions and 3D positioning guides updated in real-time during the surgical procedure. The CAS system aids in dealing with complex cases with benefits for the patient, with surgical practice, and for orthodontic finishing. Advanced software tools for diagnosis and treatment planning allow preparation of detailed operative plans, osteotomy repositioning, bone reconstructions, surgical resident training and assessing the difficulties of the surgical procedures prior to the surgery. CAS has the potential to make the elaboration of the surgical plan a more flexible process, increase the level of detail and accuracy of the plan, yield higher operative precision and control, and enhance documentation of cases. Supported by NIDCR DE017727, and DE018962 PMID:20816308
Laplace-Runge-Lenz vector for arbitrary spin
Nikitin, A. G.
2013-12-15
A countable set of superintegrable quantum mechanical systems is presented which admit the dynamical symmetry with respect to algebra so(4). This algebra is generated by the Laplace-Runge-Lenz vector generalized to the case of arbitrary spin. The presented systems describe neutral particles with non-trivial multipole momenta. Their spectra can be found algebraically like in the case of hydrogen atom. Solutions for the systems with spins 1/2 and 1 are presented explicitly, solutions for spin 3/2 can be expressed via solutions of an ordinary differential equation of first order. A more extended version of this paper including detailed calculations is published as an e-print arXiv:1308.4279.
NASA Astrophysics Data System (ADS)
Gil, José J.; San José, Ignacio
2010-11-01
From our previous definition of the indices of polarimetric purity for 3D light beams [J.J. Gil, J.M. Correas, P.A. Melero and C. Ferreira, Monogr. Semin. Mat. G. de Galdeano 31, 161 (2004)], an analysis of their geometric and physical interpretation is presented. It is found that, in agreement with previous results, the first parameter is a measure of the degree of polarization, whereas the second parameter (called the degree of directionality) is a measure of the mean angular aperture of the direction of propagation of the corresponding light beam. This pair of invariant, non-dimensional, indices of polarimetric purity contains complete information about the polarimetric purity of a light beam. The overall degree of polarimetric purity is obtained as a weighted quadratic average of the degree of polarization and the degree of directionality.
Caspi, S.; Helm, M.; Laslett, L.J.
1991-03-30
We have developed an harmonic representation for the three dimensional field components within the windings of accelerator magnets. The form by which the field is presented is suitable for interfacing with other codes that make use of the 3D field components (particle tracking and stability). The field components can be calculated with high precision and reduced cup time at any location (r,{theta},z) inside the magnet bore. The same conductor geometry which is used to simulate line currents is also used in CAD with modifications more readily available. It is our hope that the format used here for magnetic fields can be used not only as a means of delivering fields but also as a way by which beam dynamics can suggest correction to the conductor geometry. 5 refs., 70 figs.
NASA Technical Reports Server (NTRS)
2004-01-01
The Mars Exploration Rover Spirit took this 3-D navigation camera mosaic of the crater called 'Bonneville' after driving approximately 13 meters (42.7 feet) to get a better vantage point. Spirit's current position is close enough to the edge to see the interior of the crater, but high enough and far enough back to get a view of all of the walls. Because scientists and rover controllers are so pleased with this location, they will stay here for at least two more martian days, or sols, to take high resolution panoramic camera images of 'Bonneville' in its entirety. Just above the far crater rim, on the left side, is the rover's heatshield, which is visible as a tiny reflective speck.
NASA Technical Reports Server (NTRS)
1997-01-01
Many prominent rocks near the Sagan Memorial Station are featured in this image, taken in stereo by the Imager for Mars Pathfinder (IMP) on Sol 3. 3D glasses are necessary to identify surface detail. Wedge is at lower left; Shark, Half-Dome, and Pumpkin are at center. Flat Top, about four inches high, is at lower right. The horizon in the distance is one to two kilometers away.
Mars Pathfinder is the second in NASA's Discovery program of low-cost spacecraft with highly focused science goals. The Jet Propulsion Laboratory, Pasadena, CA, developed and manages the Mars Pathfinder mission for NASA's Office of Space Science, Washington, D.C. JPL is an operating division of the California Institute of Technology (Caltech). The Imager for Mars Pathfinder (IMP) was developed by the University of Arizona Lunar and Planetary Laboratory under contract to JPL. Peter Smith is the Principal Investigator.
Click below to see the left and right views individually. [figure removed for brevity, see original site] Left [figure removed for brevity, see original site] Right
NASA Technical Reports Server (NTRS)
2004-01-01
This 3-D, microscopic imager mosaic of a target area on a rock called 'Diamond Jenness' was taken after NASA's Mars Exploration Rover Opportunity ground into the surface with its rock abrasion tool for a second time.
Opportunity has bored nearly a dozen holes into the inner walls of 'Endurance Crater.' On sols 177 and 178 (July 23 and July 24, 2004), the rover worked double-duty on Diamond Jenness. Surface debris and the bumpy shape of the rock resulted in a shallow and irregular hole, only about 2 millimeters (0.08 inch) deep. The final depth was not enough to remove all the bumps and leave a neat hole with a smooth floor. This extremely shallow depression was then examined by the rover's alpha particle X-ray spectrometer.
On Sol 178, Opportunity's 'robotic rodent' dined on Diamond Jenness once again, grinding almost an additional 5 millimeters (about 0.2 inch). The rover then applied its Moessbauer spectrometer to the deepened hole. This double dose of Diamond Jenness enabled the science team to examine the rock at varying layers. Results from those grindings are currently being analyzed.
The image mosaic is about 6 centimeters (2.4 inches) across.
Shim3d Helmholtz Solution Package
Energy Science and Technology Software Center (ESTSC)
2009-01-29
This suite of codes solves the Helmholtz Equation for the steady-state propagation of single-frequency electromagnetic radiation in an arbitrary 2D or 3D dielectric medium. Materials can be either transparent or absorptive (including metals) and are described entirely by their shape and complex dielectric constant. Dielectric boundaries are assumed to always fall on grid boundaries and the material within a single grid cell is considered to be uniform. Input to the problem is in the formmore » of a Dirichlet boundary condition on a single boundary, and may be either analytic (Gaussian) in shape, or a mode shape computed using a separate code (such as the included eigenmode solver vwave20), and written to a file. Solution is via the finite difference method using Jacobi iteration for 3D problems or direct matrix inversion for 2D problems. Note that 3D problems that include metals will require different iteration parameters than described in the above reference. For structures with curved boundaries not easily modeled on a rectangular grid, the auxillary codes helmholtz11(2D), helm3d (semivectoral), and helmv3d (full vectoral) are provided. For these codes the finite difference equations are specified on a topological regular triangular grid and solved using Jacobi iteration or direct matrix inversion as before. An automatic grid generator is supplied.« less
The first eigenvalue of the Laplace operator
NASA Astrophysics Data System (ADS)
Kanguzhin, Baltabek E.; Dauitbek, Dostilek
2016-08-01
We consider a self-adjoint differential operator in the Hilbert space. The domain of the operator is changed by the perturbation of the boundary conditions so that a given neighborhood "there are no eigenvalues on neighborhood of zero" from the points of the spectrum of the perturbed operator. For the Sturm-Liouville operator on the segment and the Laplace operator on the square such a possibility is achieved through integral perturbations of boundary conditions. These statements are given with full proofs, and with a possible extension.
On Neumann and Poincare problems for Laplace equation
NASA Astrophysics Data System (ADS)
Ryazanov, Vladimir
2016-08-01
It is proved the existence of nonclassical solutions of the Neumann problem for the harmonic functions in the Jordan rectifiable domains with arbitrary measurable boundary distributions of normal derivatives. The same is stated for a special case of the Poincare problem on directional derivatives. Moreover, it is shown that the spaces of the found solutions have the infinite dimension.
Convergent radial dispersion: a Laplace transform solution for aquifer tracer testing
Moench, A.F.
1989-01-01
A Laplace transform solution was obtained for the injection of a tracer in a well situated in a homogeneous aquifer where steady, horizontal, radially convergent flow has been established due to pumping at a second well. The standard advection-dispersion equation for mass transfer was used as the controlling equation. For boundary conditions, mass balances that account for mixing of the tracer with the fluid residing in the injection and pumping wells were used. The derived solution, which can be adapted for either resident or flux-averaged concentration, is of practical use only for the pumped well. This problem is of interest because it is easily applied to field determination of aquifer dispersivity and effective porosity. Breakthrough curves were obtained by numerical inversion of the Laplace transform solution. -from Author
Waveform Inversion of Synthetic Ocean Models in the Laplace Domain
NASA Astrophysics Data System (ADS)
Rosado, H.; Blacic, T. M.; Jun, H.; Shin, C.
2014-12-01
In seismic oceanography, the processed images show where small temperature changes (as little as 0.03°C) occur, although they do not give absolute temperatures. To get a 2-D temperature map, the data must be inverted for sound speed, which is then converted to temperature using equations of state. Full waveform inversion requires a starting model that is iteratively updated until the residuals converge. Global search algorithms such as Genetic Algorithm do not require a starting model close to the true model, but are computationally exhausting. Local search inversion is less expensive, but requires a reasonably accurate starting model. Unfortunately, most marine seismic data has little associated hydrographic data and so it is difficult to create starting models close enough to the true model for convergence throughout the target area. In addition, the band-limited nature of seismic data makes it inherently challenging to extract the long wavelength sound speed trend directly from seismic data. Laplace domain inversion (LDI) developed by Changsoo Shin and colleagues requires only a rudimentary starting model to produce smooth background sound speed models without requiring prior information about the medium. It works by transforming input data to the Laplace domain, and then examining the zero frequency component of the damped wavefield to extract a smooth sound speed model - basically, removing higher frequency fluctuations to expose background trends. This ability to use frequencies below those effectively propagated by the seismic source is what enables LDI to produce the smooth background trend from the data. We applied LDI to five synthetic data sets based on simplified models of oceanographic features. Using LDI, we were able to recover smoothed versions of our synthetic models, showing the viability of the method for creating sound speed profiles suitable for use as starting models for other methods of inversion that output more detailed models.
Fragmentation functions of neutral mesons π0 and k0 with Laplace transform approach
NASA Astrophysics Data System (ADS)
Taghavi-Shahri, F.; Tehrani, S. Atashbar; Zarei, M.
2016-06-01
With an analytical solutions of DGLAP evolution equations based on the Laplace transform method, we find the fragmentation functions (FFs) of neutral mesons, π0 and k0 at NLO approximation. We also calculated the total fragmentation functions of these mesons and compared them with experimental data and those from global fits. The results show a good agreement between our solutions and other models and they are compatible with experimental data.
3D Elastic Wavefield Tomography
NASA Astrophysics Data System (ADS)
Guasch, L.; Warner, M.; Stekl, I.; Umpleby, A.; Shah, N.
2010-12-01
Wavefield tomography, or waveform inversion, aims to extract the maximum information from seismic data by matching trace by trace the response of the solid earth to seismic waves using numerical modelling tools. Its first formulation dates from the early 80's, when Albert Tarantola developed a solid theoretical basis that is still used today with little change. Due to computational limitations, the application of the method to 3D problems has been unaffordable until a few years ago, and then only under the acoustic approximation. Although acoustic wavefield tomography is widely used, a complete solution of the seismic inversion problem requires that we account properly for the physics of wave propagation, and so must include elastic effects. We have developed a 3D tomographic wavefield inversion code that incorporates the full elastic wave equation. The bottle neck of the different implementations is the forward modelling algorithm that generates the synthetic data to be compared with the field seismograms as well as the backpropagation of the residuals needed to form the direction update of the model parameters. Furthermore, one or two extra modelling runs are needed in order to calculate the step-length. Our approach uses a FD scheme explicit time-stepping by finite differences that are 4th order in space and 2nd order in time, which is a 3D version of the one developed by Jean Virieux in 1986. We chose the time domain because an explicit time scheme is much less demanding in terms of memory than its frequency domain analogue, although the discussion of wich domain is more efficient still remains open. We calculate the parameter gradients for Vp and Vs by correlating the normal and shear stress wavefields respectively. A straightforward application would lead to the storage of the wavefield at all grid points at each time-step. We tackled this problem using two different approaches. The first one makes better use of resources for small models of dimension equal
ERIC Educational Resources Information Center
Iqbal, M.
2002-01-01
In this paper we have converted the Laplace transform into an integral equation of the first kind of convolution type, which is an ill-posed problem, and used a statistical regularization method to solve it. The method is applied to three examples. It gives a good approximation to the true solution and compares well with the method given by…
NASA Astrophysics Data System (ADS)
Mehrabi Pari, Sharareh; Javidan, Kurosh; Taghavi Shahri, Fatemeh
2016-05-01
The "Laplace Transform Method" is used to solve the Fokker-Plank equation for finding the time evolution of the heavy quarks distribution functions such as charm and bottom in quark gluon plasma. These solutions will lead us to calculation of nuclear suppression factor RAA. The results have good agreement with available experiment data from the PHENIX collaboration.
Towards Informetrics: Haitun, Laplace, Zipf, Bradford and the Alvey Programme.
ERIC Educational Resources Information Center
Brookes, B. C.
1984-01-01
Review of recent developments in statistical theories for social sciences highlights Haitun's statistical distributions, Laplace's "Law of Succession" and distribution, Laplace and Bradford analysis of book-index data, inefficiency of frequency distribution analysis, Laws of Bradford and Zipf, natural categorization, and Bradford Law and household…
Feeling Wall Tension in an Interactive Demonstration of Laplace's Law
ERIC Educational Resources Information Center
Letic, Milorad
2012-01-01
Laplace's Law plays a major role in explanations of the wall tension of structures like blood vessels, the bladder, the uterus in pregnancy, bronchioles, eyeballs, and the behavior of aneurisms or the enlarged heart. The general relation of Laplace's law, expressing that the product of the radius of curvature (r) and pressure (P) is equal to wall…
NASA Astrophysics Data System (ADS)
Mediavilla, Evencio; Arribas, Santiago; Roth, Martin; Cepa-Nogué, Jordi; Sánchez, Francisco
2011-09-01
Preface; Acknowledgements; 1. Introductory review and technical approaches Martin M. Roth; 2. Observational procedures and data reduction James E. H. Turner; 3. 3D Spectroscopy instrumentation M. A. Bershady; 4. Analysis of 3D data Pierre Ferruit; 5. Science motivation for IFS and galactic studies F. Eisenhauer; 6. Extragalactic studies and future IFS science Luis Colina; 7. Tutorials: how to handle 3D spectroscopy data Sebastian F. Sánchez, Begona García-Lorenzo and Arlette Pécontal-Rousset.
Radial fractional Laplace operators and Hessian inequalities
NASA Astrophysics Data System (ADS)
Ferrari, Fausto; Verbitsky, Igor E.
In this paper we deduce a formula for the fractional Laplace operator ( on radially symmetric functions useful for some applications. We give a criterion of subharmonicity associated with (, and apply it to a problem related to the Hessian inequality of Sobolev type: ∫Rn |(u| dx⩽C∫Rn -uFk[u] dx, where Fk is the k-Hessian operator on Rn, 1⩽k
Animal Models and Integrated Nested Laplace Approximations
Holand, Anna Marie; Steinsland, Ingelin; Martino, Sara; Jensen, Henrik
2013-01-01
Animal models are generalized linear mixed models used in evolutionary biology and animal breeding to identify the genetic part of traits. Integrated Nested Laplace Approximation (INLA) is a methodology for making fast, nonsampling-based Bayesian inference for hierarchical Gaussian Markov models. In this article, we demonstrate that the INLA methodology can be used for many versions of Bayesian animal models. We analyze animal models for both synthetic case studies and house sparrow (Passer domesticus) population case studies with Gaussian, binomial, and Poisson likelihoods using INLA. Inference results are compared with results using Markov Chain Monte Carlo methods. For model choice we use difference in deviance information criteria (DIC). We suggest and show how to evaluate differences in DIC by comparing them with sampling results from simulation studies. We also introduce an R package, AnimalINLA, for easy and fast inference for Bayesian Animal models using INLA. PMID:23708299
3D Elevation Program—Virtual USA in 3D
Lukas, Vicki; Stoker, J.M.
2016-01-01
The U.S. Geological Survey (USGS) 3D Elevation Program (3DEP) uses a laser system called ‘lidar’ (light detection and ranging) to create a virtual reality map of the Nation that is very accurate. 3D maps have many uses with new uses being discovered all the time.
Concurrent 3-D motion segmentation and 3-D interpretation of temporal sequences of monocular images.
Sekkati, Hicham; Mitiche, Amar
2006-03-01
The purpose of this study is to investigate a variational method for joint multiregion three-dimensional (3-D) motion segmentation and 3-D interpretation of temporal sequences of monocular images. Interpretation consists of dense recovery of 3-D structure and motion from the image sequence spatiotemporal variations due to short-range image motion. The method is direct insomuch as it does not require prior computation of image motion. It allows movement of both viewing system and multiple independently moving objects. The problem is formulated following a variational statement with a functional containing three terms. One term measures the conformity of the interpretation within each region of 3-D motion segmentation to the image sequence spatiotemporal variations. The second term is of regularization of depth. The assumption that environmental objects are rigid accounts automatically for the regularity of 3-D motion within each region of segmentation. The third and last term is for the regularity of segmentation boundaries. Minimization of the functional follows the corresponding Euler-Lagrange equations. This results in iterated concurrent computation of 3-D motion segmentation by curve evolution, depth by gradient descent, and 3-D motion by least squares within each region of segmentation. Curve evolution is implemented via level sets for topology independence and numerical stability. This algorithm and its implementation are verified on synthetic and real image sequences. Viewers presented with anaglyphs of stereoscopic images constructed from the algorithm's output reported a strong perception of depth. PMID:16519351
MT3D was first developed by Chunmiao Zheng in 1990 at S.S. Papadopulos & Associates, Inc. with partial support from the U.S. Environmental Protection Agency (USEPA). Starting in 1990, MT3D was released as a pubic domain code from the USEPA. Commercial versions with enhanced capab...
NASA Technical Reports Server (NTRS)
1977-01-01
A market study of a proposed version of a 3-D eyetracker for initial use at NASA's Ames Research Center was made. The commercialization potential of a simplified, less expensive 3-D eyetracker was ascertained. Primary focus on present and potential users of eyetrackers, as well as present and potential manufacturers has provided an effective means of analyzing the prospects for commercialization.
Energy Science and Technology Software Center (ESTSC)
2013-10-01
Earth3D is a computer code designed to allow fast calculation of seismic rays and travel times through a 3D model of the Earth. LLNL is using this for earthquake location and global tomography efforts and such codes are of great interest to the Earth Science community.
[3-D ultrasound in gastroenterology].
Zoller, W G; Liess, H
1994-06-01
Three-dimensional (3D) sonography represents a development of noninvasive diagnostic imaging by real-time two-dimensional (2D) sonography. The use of transparent rotating scans, comparable to a block of glass, generates a 3D effect. The objective of the present study was to optimate 3D presentation of abdominal findings. Additional investigations were made with a new volumetric program to determine the volume of selected findings of the liver. The results were compared with the estimated volumes of 2D sonography and 2D computer tomography (CT). For the processing of 3D images, typical parameter constellations were found for the different findings, which facilitated processing of 3D images. In more than 75% of the cases examined we found an optimal 3D presentation of sonographic findings with respect to the evaluation criteria developed by us for the 3D imaging of processed data. Great differences were found for the estimated volumes of the findings of the liver concerning the three different techniques applied. 3D ultrasound represents a valuable method to judge morphological appearance in abdominal findings. The possibility of volumetric measurements enlarges its potential diagnostic significance. Further clinical investigations are necessary to find out if definite differentiation between benign and malign findings is possible. PMID:7919882
2013-10-30
This video provides an overview of the Sandia National Laboratories developed 3-D World Model Building capability that provides users with an immersive, texture rich 3-D model of their environment in minutes using a laptop and color and depth camera.
None
2014-02-26
This video provides an overview of the Sandia National Laboratories developed 3-D World Model Building capability that provides users with an immersive, texture rich 3-D model of their environment in minutes using a laptop and color and depth camera.
Analysis of fluorescence decay curves by means of the Laplace transformation.
Gafni, A; Modlin, R L; Brand, L
1975-03-01
A computational procedure is described for the analysis of fluorescence decay data convolved with a lamp flash of finite width. The computer program calculates the ratio of the Laplace transforms of the decay and the lamp flash for different values of s to give the transforms of the impulse response for each value of s. These are set equal to the analytical Laplace transforms of the decay law involved. Solution of the nonlinear simultaneous equations yields the desired decay parameters. The method can be modified to analyze data that contains a component due to scattered light and can also provide essential information regarding transit time changes of the photomultiplier with changes in emission wavelength. The method was tested by the analysis of real and simulated data. The accuracy of the analysis depends on the degree of correlation among the parameters. PMID:1122338
NASA Astrophysics Data System (ADS)
Walsh, J. R.
2004-02-01
The Euro3D RTN is an EU funded Research Training Network to foster the exploitation of 3D spectroscopy in Europe. 3D spectroscopy is a general term for spectroscopy of an area of the sky and derives its name from its two spatial + one spectral dimensions. There are an increasing number of instruments which use integral field devices to achieve spectroscopy of an area of the sky, either using lens arrays, optical fibres or image slicers, to pack spectra of multiple pixels on the sky (``spaxels'') onto a 2D detector. On account of the large volume of data and the special methods required to reduce and analyse 3D data, there are only a few centres of expertise and these are mostly involved with instrument developments. There is a perceived lack of expertise in 3D spectroscopy spread though the astronomical community and its use in the armoury of the observational astronomer is viewed as being highly specialised. For precisely this reason the Euro3D RTN was proposed to train young researchers in this area and develop user tools to widen the experience with this particular type of data in Europe. The Euro3D RTN is coordinated by Martin M. Roth (Astrophysikalisches Institut Potsdam) and has been running since July 2002. The first Euro3D science conference was held in Cambridge, UK from 22 to 23 May 2003. The main emphasis of the conference was, in keeping with the RTN, to expose the work of the young post-docs who are funded by the RTN. In addition the team members from the eleven European institutes involved in Euro3D also presented instrumental and observational developments. The conference was organized by Andy Bunker and held at the Institute of Astronomy. There were over thirty participants and 26 talks covered the whole range of application of 3D techniques. The science ranged from Galactic planetary nebulae and globular clusters to kinematics of nearby galaxies out to objects at high redshift. Several talks were devoted to reporting recent observations with newly