A numerical comparison of discrete Kalman filtering algorithms: An orbit determination case study
NASA Technical Reports Server (NTRS)
Thornton, C. L.; Bierman, G. J.
1976-01-01
The numerical stability and accuracy of various Kalman filter algorithms are thoroughly studied. Numerical results and conclusions are based on a realistic planetary approach orbit determination study. The case study results of this report highlight the numerical instability of the conventional and stabilized Kalman algorithms. Numerical errors associated with these algorithms can be so large as to obscure important mismodeling effects and thus give misleading estimates of filter accuracy. The positive result of this study is that the Bierman-Thornton U-D covariance factorization algorithm is computationally efficient, with CPU costs that differ negligibly from the conventional Kalman costs. In addition, accuracy of the U-D filter using single-precision arithmetic consistently matches the double-precision reference results. Numerical stability of the U-D filter is further demonstrated by its insensitivity of variations in the a priori statistics.
Comparative Study of Algorithms for the Numerical Simulation of Lattice QCD
Luz, Fernando H. P.; Mendes, Tereza
2010-11-12
Large-scale numerical simulations are the prime method for a nonperturbative study of QCD from first principles. Although the lattice simulation of the pure-gauge (or quenched-QCD) case may be performed very efficiently on parallel machines, there are several additional difficulties in the simulation of the full-QCD case, i.e. when dynamical quark effects are taken into account. We discuss the main aspects of full-QCD simulations, describing the most common algorithms. We present a comparative analysis of performance for two versions of the hybrid Monte Carlo method (the so-called R and RHMC algorithms), as provided in the MILC software package. We consider two degenerate flavors of light quarks in the staggered formulation, having in mind the case of finite-temperature QCD.
A numerical comparison of discrete Kalman filtering algorithms - An orbit determination case study
NASA Technical Reports Server (NTRS)
Thornton, C. L.; Bierman, G. J.
1976-01-01
An improved Kalman filter algorithm based on a modified Givens matrix triangularization technique is proposed for solving a nonstationary discrete-time linear filtering problem. The proposed U-D covariance factorization filter uses orthogonal transformation technique; measurement and time updating of the U-D factors involve separate application of Gentleman's fast square-root-free Givens rotations. Numerical stability and accuracy of the algorithm are compared with those of the conventional and stabilized Kalman filters and the Potter-Schmidt square-root filter, by applying these techniques to a realistic planetary navigation problem (orbit determination for the Saturn approach phase of the Mariner Jupiter-Saturn Mission, 1977). The new algorithm is shown to combine the numerical precision of square root filtering with the efficiency of the original Kalman algorithm.
Biphasic indentation of articular cartilage--II. A numerical algorithm and an experimental study.
Mow, V C; Gibbs, M C; Lai, W M; Zhu, W B; Athanasiou, K A
1989-01-01
Part I (Mak et al., 1987, J. Biomechanics 20, 703-714) presented the theoretical solutions for the biphasic indentation of articular cartilage under creep and stress-relaxation conditions. In this study, using the creep solution, we developed an efficient numerical algorithm to compute all three material coefficients of cartilage in situ on the joint surface from the indentation creep experiment. With this method we determined the average values of the aggregate modulus. Poisson's ratio and permeability for young bovine femoral condylar cartilage in situ to be HA = 0.90 MPa, vs = 0.39 and k = 0.44 x 10(-15) m4/Ns respectively, and those for patellar groove cartilage to be HA = 0.47 MPa, vs = 0.24, k = 1.42 x 10(-15) m4/Ns. One surprising finding from this study is that the in situ Poisson's ratio of cartilage (0.13-0.45) may be much less than those determined from measurements performed on excised osteochondral plugs (0.40-0.49) reported in the literature. We also found the permeability of patellar groove cartilage to be several times higher than femoral condyle cartilage. These findings may have important implications on understanding the functional behavior of cartilage in situ and on methods used to determine the elastic moduli of cartilage using the indentation experiments. PMID:2613721
Numerical Study of Three-Dimensional Flows Using Unfactored Upwind-Relaxation Sweeping Algorithm
NASA Astrophysics Data System (ADS)
Zha, G.-C.; Bilgen, E.
1996-05-01
The linear stability analysis of the unfactored upwind relaxation-sweeping (URS) algorithm for 3D flow field calculations has been carried out and it is shown that the URS algorithm is unconditionally stable. The algorithm is independent of the global sweeping direction selection. However, choosing the direction with relatively low variable gradient as the global sweeping direction results in a higher degree of stability. Three-dimensional compressible Euler equations are solved by using the implicit URS algorithm to study internal flows of a non-axisymmetric nozzle with a circular-to-rectangular transition duct and complex shock wave structures for a 3D channel flow. The efficiency and robustness of the URS algorithm has been demonstrated.
NASA Astrophysics Data System (ADS)
Miao, J. C.; Zhu, P.; Shi, G. L.; Chen, G. L.
2008-01-01
Numerical stability is an important issue for any integral procedure. Since sub-cycling algorithm has been presented by Belytschko et al. (Comput Methods Appl Mech Eng 17/18: 259-275, 1979), various kinds of these integral procedures were developed in later 20 years and their stability were widely studied. However, on how to apply the sub-cycling to flexible multi-body dynamics (FMD) is still a lack of investigation up to now. A particular sub-cycling algorithm for the FMD based on the central difference method was introduced in detail in part I (Miao et al. in Comp Mech doi: 10.1007/s00466-007-0183-9) of this paper. Adopting an integral approximation operator method, stability of the presented algorithm is transformed to a generalized eigenvalue problem in the paper and is discussed by solving the problem later. Numerical examples are performed to verify the availability and efficiency of the algorithm further.
Trees, bialgebras and intrinsic numerical algorithms
NASA Technical Reports Server (NTRS)
Crouch, Peter; Grossman, Robert; Larson, Richard
1990-01-01
Preliminary work about intrinsic numerical integrators evolving on groups is described. Fix a finite dimensional Lie group G; let g denote its Lie algebra, and let Y(sub 1),...,Y(sub N) denote a basis of g. A class of numerical algorithms is presented that approximate solutions to differential equations evolving on G of the form: dot-x(t) = F(x(t)), x(0) = p is an element of G. The algorithms depend upon constants c(sub i) and c(sub ij), for i = 1,...,k and j is less than i. The algorithms have the property that if the algorithm starts on the group, then it remains on the group. In addition, they also have the property that if G is the abelian group R(N), then the algorithm becomes the classical Runge-Kutta algorithm. The Cayley algebra generated by labeled, ordered trees is used to generate the equations that the coefficients c(sub i) and c(sub ij) must satisfy in order for the algorithm to yield an rth order numerical integrator and to analyze the resulting algorithms.
Numerical Algorithms Based on Biorthogonal Wavelets
NASA Technical Reports Server (NTRS)
Ponenti, Pj.; Liandrat, J.
1996-01-01
Wavelet bases are used to generate spaces of approximation for the resolution of bidimensional elliptic and parabolic problems. Under some specific hypotheses relating the properties of the wavelets to the order of the involved operators, it is shown that an approximate solution can be built. This approximation is then stable and converges towards the exact solution. It is designed such that fast algorithms involving biorthogonal multi resolution analyses can be used to resolve the corresponding numerical problems. Detailed algorithms are provided as well as the results of numerical tests on partial differential equations defined on the bidimensional torus.
Numerical Algorithm for Delta of Asian Option
Zhang, Boxiang; Yu, Yang; Wang, Weiguo
2015-01-01
We study the numerical solution of the Greeks of Asian options. In particular, we derive a close form solution of Δ of Asian geometric option and use this analytical form as a control to numerically calculate Δ of Asian arithmetic option, which is known to have no explicit close form solution. We implement our proposed numerical method and compare the standard error with other classical variance reduction methods. Our method provides an efficient solution to the hedging strategy with Asian options. PMID:26266271
Numerical Algorithm for Delta of Asian Option.
Zhang, Boxiang; Yu, Yang; Wang, Weiguo
2015-01-01
We study the numerical solution of the Greeks of Asian options. In particular, we derive a close form solution of Δ of Asian geometric option and use this analytical form as a control to numerically calculate Δ of Asian arithmetic option, which is known to have no explicit close form solution. We implement our proposed numerical method and compare the standard error with other classical variance reduction methods. Our method provides an efficient solution to the hedging strategy with Asian options. PMID:26266271
Stochastic Formal Correctness of Numerical Algorithms
NASA Technical Reports Server (NTRS)
Daumas, Marc; Lester, David; Martin-Dorel, Erik; Truffert, Annick
2009-01-01
We provide a framework to bound the probability that accumulated errors were never above a given threshold on numerical algorithms. Such algorithms are used for example in aircraft and nuclear power plants. This report contains simple formulas based on Levy's and Markov's inequalities and it presents a formal theory of random variables with a special focus on producing concrete results. We selected four very common applications that fit in our framework and cover the common practices of systems that evolve for a long time. We compute the number of bits that remain continuously significant in the first two applications with a probability of failure around one out of a billion, where worst case analysis considers that no significant bit remains. We are using PVS as such formal tools force explicit statement of all hypotheses and prevent incorrect uses of theorems.
Experimentally constructing finite difference algorithms in numerical relativity
NASA Astrophysics Data System (ADS)
Anderson, Matthew; Neilsen, David; Matzner, Richard
2002-04-01
Computational studies of gravitational waves require numerical algorithms with long-term stability (necessary for convergence). However, constructing stable finite difference algorithms (FDA) for the ADM formulation of the Einstein equations, especially in multiple dimensions, has proven difficult. Most FDA's are constructed using rules of thumb gained from experience with simple model equations. To search for FDA's with improved stability, we adopt a brute-force approach, where we systematically test thousands of numerical schemes. We sort the spatial derivatives of the Einstein equations into groups, and parameterize each group by finite difference type (centered or upwind) and order. Furthermore, terms proportional to the constraints are added to the evolution equations with additional parameters. A spherically symmetric, excised Schwarzschild black hole (one dimension) and linearized waves in multiple dimensions are used as model systems to evaluate the different numerical schemes.
Adaptive Numerical Algorithms in Space Weather Modeling
NASA Technical Reports Server (NTRS)
Toth, Gabor; vanderHolst, Bart; Sokolov, Igor V.; DeZeeuw, Darren; Gombosi, Tamas I.; Fang, Fang; Manchester, Ward B.; Meng, Xing; Nakib, Dalal; Powell, Kenneth G.; Stout, Quentin F.; Glocer, Alex; Ma, Ying-Juan; Opher, Merav
2010-01-01
Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical
Adaptive numerical algorithms in space weather modeling
NASA Astrophysics Data System (ADS)
Tóth, Gábor; van der Holst, Bart; Sokolov, Igor V.; De Zeeuw, Darren L.; Gombosi, Tamas I.; Fang, Fang; Manchester, Ward B.; Meng, Xing; Najib, Dalal; Powell, Kenneth G.; Stout, Quentin F.; Glocer, Alex; Ma, Ying-Juan; Opher, Merav
2012-02-01
Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different relevant physics in different domains. A multi-physics system can be modeled by a software framework comprising several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamic (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit
Research on numerical algorithms for large space structures
NASA Technical Reports Server (NTRS)
Denman, E. D.
1982-01-01
Numerical algorithms for large space structures were investigated with particular emphasis on decoupling method for analysis and design. Numerous aspects of the analysis of large systems ranging from the algebraic theory to lambda matrices to identification algorithms were considered. A general treatment of the algebraic theory of lambda matrices is presented and the theory is applied to second order lambda matrices.
A Polynomial Time, Numerically Stable Integer Relation Algorithm
NASA Technical Reports Server (NTRS)
Ferguson, Helaman R. P.; Bailey, Daivd H.; Kutler, Paul (Technical Monitor)
1998-01-01
Let x = (x1, x2...,xn be a vector of real numbers. X is said to possess an integer relation if there exist integers a(sub i) not all zero such that a1x1 + a2x2 + ... a(sub n)Xn = 0. Beginning in 1977 several algorithms (with proofs) have been discovered to recover the a(sub i) given x. The most efficient of these existing integer relation algorithms (in terms of run time and the precision required of the input) has the drawback of being very unstable numerically. It often requires a numeric precision level in the thousands of digits to reliably recover relations in modest-sized test problems. We present here a new algorithm for finding integer relations, which we have named the "PSLQ" algorithm. It is proved in this paper that the PSLQ algorithm terminates with a relation in a number of iterations that is bounded by a polynomial in it. Because this algorithm employs a numerically stable matrix reduction procedure, it is free from the numerical difficulties, that plague other integer relation algorithms. Furthermore, its stability admits an efficient implementation with lower run times oil average than other algorithms currently in Use. Finally, this stability can be used to prove that relation bounds obtained from computer runs using this algorithm are numerically accurate.
Linsen, Sarah; Torbeyns, Joke; Verschaffel, Lieven; Reynvoet, Bert; De Smedt, Bert
2016-03-01
There are two well-known computation methods for solving multi-digit subtraction items, namely mental and algorithmic computation. It has been contended that mental and algorithmic computation differentially rely on numerical magnitude processing, an assumption that has already been examined in children, but not yet in adults. Therefore, in this study, we examined how numerical magnitude processing was associated with mental and algorithmic computation, and whether this association with numerical magnitude processing was different for mental versus algorithmic computation. We also investigated whether the association between numerical magnitude processing and mental and algorithmic computation differed for measures of symbolic versus nonsymbolic numerical magnitude processing. Results showed that symbolic, and not nonsymbolic, numerical magnitude processing was associated with mental computation, but not with algorithmic computation. Additional analyses showed, however, that the size of this association with symbolic numerical magnitude processing was not significantly different for mental and algorithmic computation. We also tried to further clarify the association between numerical magnitude processing and complex calculation by also including relevant arithmetical subskills, i.e. arithmetic facts, needed for complex calculation that are also known to be dependent on numerical magnitude processing. Results showed that the associations between symbolic numerical magnitude processing and mental and algorithmic computation were fully explained by individual differences in elementary arithmetic fact knowledge. PMID:26914586
A Numerical Instability in an ADI Algorithm for Gyrokinetics
E.A. Belli; G.W. Hammett
2004-12-17
We explore the implementation of an Alternating Direction Implicit (ADI) algorithm for a gyrokinetic plasma problem and its resulting numerical stability properties. This algorithm, which uses a standard ADI scheme to divide the field solve from the particle distribution function advance, has previously been found to work well for certain plasma kinetic problems involving one spatial and two velocity dimensions, including collisions and an electric field. However, for the gyrokinetic problem we find a severe stability restriction on the time step. Furthermore, we find that this numerical instability limitation also affects some other algorithms, such as a partially implicit Adams-Bashforth algorithm, where the parallel motion operator v{sub {parallel}} {partial_derivative}/{partial_derivative}z is treated implicitly and the field terms are treated with an Adams-Bashforth explicit scheme. Fully explicit algorithms applied to all terms can be better at long wavelengths than these ADI or partially implicit algorithms.
Research on numerical algorithms for large space structures
NASA Technical Reports Server (NTRS)
Denman, E. D.
1981-01-01
Numerical algorithms for analysis and design of large space structures are investigated. The sign algorithm and its application to decoupling of differential equations are presented. The generalized sign algorithm is given and its application to several problems discussed. The Laplace transforms of matrix functions and the diagonalization procedure for a finite element equation are discussed. The diagonalization of matrix polynomials is considered. The quadrature method and Laplace transforms is discussed and the identification of linear systems by the quadrature method investigated.
Technical Report: Scalable Parallel Algorithms for High Dimensional Numerical Integration
Masalma, Yahya; Jiao, Yu
2010-10-01
We implemented a scalable parallel quasi-Monte Carlo numerical high-dimensional integration for tera-scale data points. The implemented algorithm uses the Sobol s quasi-sequences to generate random samples. Sobol s sequence was used to avoid clustering effects in the generated random samples and to produce low-discrepancy random samples which cover the entire integration domain. The performance of the algorithm was tested. Obtained results prove the scalability and accuracy of the implemented algorithms. The implemented algorithm could be used in different applications where a huge data volume is generated and numerical integration is required. We suggest using the hyprid MPI and OpenMP programming model to improve the performance of the algorithms. If the mixed model is used, attention should be paid to the scalability and accuracy.
A hybrid artificial bee colony algorithm for numerical function optimization
NASA Astrophysics Data System (ADS)
Alqattan, Zakaria N.; Abdullah, Rosni
2015-02-01
Artificial Bee Colony (ABC) algorithm is one of the swarm intelligence algorithms; it has been introduced by Karaboga in 2005. It is a meta-heuristic optimization search algorithm inspired from the intelligent foraging behavior of the honey bees in nature. Its unique search process made it as one of the most competitive algorithm with some other search algorithms in the area of optimization, such as Genetic algorithm (GA) and Particle Swarm Optimization (PSO). However, the ABC performance of the local search process and the bee movement or the solution improvement equation still has some weaknesses. The ABC is good in avoiding trapping at the local optimum but it spends its time searching around unpromising random selected solutions. Inspired by the PSO, we propose a Hybrid Particle-movement ABC algorithm called HPABC, which adapts the particle movement process to improve the exploration of the original ABC algorithm. Numerical benchmark functions were used in order to experimentally test the HPABC algorithm. The results illustrate that the HPABC algorithm can outperform the ABC algorithm in most of the experiments (75% better in accuracy and over 3 times faster).
Understanding disordered systems through numerical simulation and algorithm development
NASA Astrophysics Data System (ADS)
Sweeney, Sean Michael
Disordered systems arise in many physical contexts. Not all matter is uniform, and impurities or heterogeneities can be modeled by fixed random disorder. Numerous complex networks also possess fixed disorder, leading to applications in transportation systems, telecommunications, social networks, and epidemic modeling, to name a few. Due to their random nature and power law critical behavior, disordered systems are difficult to study analytically. Numerical simulation can help overcome this hurdle by allowing for the rapid computation of system states. In order to get precise statistics and extrapolate to the thermodynamic limit, large systems must be studied over many realizations. Thus, innovative algorithm development is essential in order reduce memory or running time requirements of simulations. This thesis presents a review of disordered systems, as well as a thorough study of two particular systems through numerical simulation, algorithm development and optimization, and careful statistical analysis of scaling properties. Chapter 1 provides a thorough overview of disordered systems, the history of their study in the physics community, and the development of techniques used to study them. Topics of quenched disorder, phase transitions, the renormalization group, criticality, and scale invariance are discussed. Several prominent models of disordered systems are also explained. Lastly, analysis techniques used in studying disordered systems are covered. In Chapter 2, minimal spanning trees on critical percolation clusters are studied, motivated in part by an analytic perturbation expansion by Jackson and Read that I check against numerical calculations. This system has a direct mapping to the ground state of the strongly disordered spin glass. We compute the path length fractal dimension of these trees in dimensions d = {2, 3, 4, 5} and find our results to be compatible with the analytic results suggested by Jackson and Read. In Chapter 3, the random bond Ising
An efficient cuckoo search algorithm for numerical function optimization
NASA Astrophysics Data System (ADS)
Ong, Pauline; Zainuddin, Zarita
2013-04-01
Cuckoo search algorithm which reproduces the breeding strategy of the best known brood parasitic bird, the cuckoos has demonstrated its superiority in obtaining the global solution for numerical optimization problems. However, the involvement of fixed step approach in its exploration and exploitation behavior might slow down the search process considerably. In this regards, an improved cuckoo search algorithm with adaptive step size adjustment is introduced and its feasibility on a variety of benchmarks is validated. The obtained results show that the proposed scheme outperforms the standard cuckoo search algorithm in terms of convergence characteristic while preserving the fascinating features of the original method.
Multiresolution representation and numerical algorithms: A brief review
NASA Technical Reports Server (NTRS)
Harten, Amiram
1994-01-01
In this paper we review recent developments in techniques to represent data in terms of its local scale components. These techniques enable us to obtain data compression by eliminating scale-coefficients which are sufficiently small. This capability for data compression can be used to reduce the cost of many numerical solution algorithms by either applying it to the numerical solution operator in order to get an approximate sparse representation, or by applying it to the numerical solution itself in order to reduce the number of quantities that need to be computed.
Fast Quantum Algorithms for Numerical Integrals and Stochastic Processes
NASA Technical Reports Server (NTRS)
Abrams, D.; Williams, C.
1999-01-01
We discuss quantum algorithms that calculate numerical integrals and descriptive statistics of stochastic processes. With either of two distinct approaches, one obtains an exponential speed increase in comparison to the fastest known classical deterministic algotithms and a quadratic speed increase incomparison to classical Monte Carlo methods.
A novel bee swarm optimization algorithm for numerical function optimization
NASA Astrophysics Data System (ADS)
Akbari, Reza; Mohammadi, Alireza; Ziarati, Koorush
2010-10-01
The optimization algorithms which are inspired from intelligent behavior of honey bees are among the most recently introduced population based techniques. In this paper, a novel algorithm called bee swarm optimization, or BSO, and its two extensions for improving its performance are presented. The BSO is a population based optimization technique which is inspired from foraging behavior of honey bees. The proposed approach provides different patterns which are used by the bees to adjust their flying trajectories. As the first extension, the BSO algorithm introduces different approaches such as repulsion factor and penalizing fitness (RP) to mitigate the stagnation problem. Second, to maintain efficiently the balance between exploration and exploitation, time-varying weights (TVW) are introduced into the BSO algorithm. The proposed algorithm (BSO) and its two extensions (BSO-RP and BSO-RPTVW) are compared with existing algorithms which are based on intelligent behavior of honey bees, on a set of well known numerical test functions. The experimental results show that the BSO algorithms are effective and robust; produce excellent results, and outperform other algorithms investigated in this consideration.
Determining the Numerical Stability of Quantum Chemistry Algorithms.
Knizia, Gerald; Li, Wenbin; Simon, Sven; Werner, Hans-Joachim
2011-08-01
We present a simple, broadly applicable method for determining the numerical properties of quantum chemistry algorithms. The method deliberately introduces random numerical noise into computations, which is of the same order of magnitude as the floating point precision. Accordingly, repeated runs of an algorithm give slightly different results, which can be analyzed statistically to obtain precise estimates of its numerical stability. This noise is produced by automatic code injection into regular compiler output, so that no substantial programming effort is required, only a recompilation of the affected program sections. The method is applied to investigate: (i) the numerical stability of the three-center Obara-Saika integral evaluation scheme for high angular momenta, (ii) if coupled cluster perturbative triples can be evaluated with single precision arithmetic, (iii) how to implement the density fitting approximation in Møller-Plesset perturbation theory (MP2) most accurately, and (iv) which parts of density fitted MP2 can be safely evaluated with single precision arithmetic. In the integral case, we find a numerical instability in an equation that is used in almost all integral programs. Due to the results of (ii) and (iv), we conjecture that single precision arithmetic can be applied whenever a calculation is done in an orthogonal basis set and excessively long linear sums are avoided. PMID:26606614
An algorithm for the numerical solution of linear differential games
Polovinkin, E S; Ivanov, G E; Balashov, M V; Konstantinov, R V; Khorev, A V
2001-10-31
A numerical algorithm for the construction of stable Krasovskii bridges, Pontryagin alternating sets, and also of piecewise program strategies solving two-person linear differential (pursuit or evasion) games on a fixed time interval is developed on the basis of a general theory. The aim of the first player (the pursuer) is to hit a prescribed target (terminal) set by the phase vector of the control system at the prescribed time. The aim of the second player (the evader) is the opposite. A description of numerical algorithms used in the solution of differential games of the type under consideration is presented and estimates of the errors resulting from the approximation of the game sets by polyhedra are presented.
A numerical algorithm for magnetohydrodynamics of ablated materials.
Lu, Tianshi; Du, Jian; Samulyak, Roman
2008-07-01
A numerical algorithm for the simulation of magnetohydrodynamics in partially ionized ablated material is described. For the hydro part, the hyperbolic conservation laws with electromagnetic terms is solved using techniques developed for free surface flows; for the electromagnetic part, the electrostatic approximation is applied and an elliptic equation for electric potential is solved. The algorithm has been implemented in the frame of front tracking, which explicitly tracks geometrically complex evolving interfaces. An elliptic solver based on the embedded boundary method were implemented for both two- and three-dimensional simulations. A surface model on the interface between the solid target and the ablated vapor has also been developed as well as a numerical model for the equation of state which accounts for atomic processes in the ablated material. The code has been applied to simulations of the pellet ablation in a magnetically confined plasma and the laser-ablated plasma plume expansion in magnetic fields. PMID:19051925
Algorithms for the Fractional Calculus: A Selection of Numerical Methods
NASA Technical Reports Server (NTRS)
Diethelm, K.; Ford, N. J.; Freed, A. D.; Luchko, Yu.
2003-01-01
Many recently developed models in areas like viscoelasticity, electrochemistry, diffusion processes, etc. are formulated in terms of derivatives (and integrals) of fractional (non-integer) order. In this paper we present a collection of numerical algorithms for the solution of the various problems arising in this context. We believe that this will give the engineer the necessary tools required to work with fractional models in an efficient way.
Copps, Kevin D.; Carnes, Brian R.
2008-04-01
We examine algorithms for the finite element approximation of thermal contact models. We focus on the implementation of thermal contact algorithms in SIERRA Mechanics. Following the mathematical formulation of models for tied contact and resistance contact, we present three numerical algorithms: (1) the multi-point constraint (MPC) algorithm, (2) a resistance algorithm, and (3) a new generalized algorithm. We compare and contrast both the correctness and performance of the algorithms in three test problems. We tabulate the convergence rates of global norms of the temperature solution on sequentially refined meshes. We present the results of a parameter study of the effect of contact search tolerances. We outline best practices in using the software for predictive simulations, and suggest future improvements to the implementation.
Predictive Lateral Logic for Numerical Entry Guidance Algorithms
NASA Technical Reports Server (NTRS)
Smith, Kelly M.
2016-01-01
Recent entry guidance algorithm development123 has tended to focus on numerical integration of trajectories onboard in order to evaluate candidate bank profiles. Such methods enjoy benefits such as flexibility to varying mission profiles and improved robustness to large dispersions. A common element across many of these modern entry guidance algorithms is a reliance upon the concept of Apollo heritage lateral error (or azimuth error) deadbands in which the number of bank reversals to be performed is non-deterministic. This paper presents a closed-loop bank reversal method that operates with a fixed number of bank reversals defined prior to flight. However, this number of bank reversals can be modified at any point, including in flight, based on contingencies such as fuel leaks where propellant usage must be minimized.
Algorithm-Based Fault Tolerance for Numerical Subroutines
NASA Technical Reports Server (NTRS)
Tumon, Michael; Granat, Robert; Lou, John
2007-01-01
A software library implements a new methodology of detecting faults in numerical subroutines, thus enabling application programs that contain the subroutines to recover transparently from single-event upsets. The software library in question is fault-detecting middleware that is wrapped around the numericalsubroutines. Conventional serial versions (based on LAPACK and FFTW) and a parallel version (based on ScaLAPACK) exist. The source code of the application program that contains the numerical subroutines is not modified, and the middleware is transparent to the user. The methodology used is a type of algorithm- based fault tolerance (ABFT). In ABFT, a checksum is computed before a computation and compared with the checksum of the computational result; an error is declared if the difference between the checksums exceeds some threshold. Novel normalization methods are used in the checksum comparison to ensure correct fault detections independent of algorithm inputs. In tests of this software reported in the peer-reviewed literature, this library was shown to enable detection of 99.9 percent of significant faults while generating no false alarms.
Numerical Studies of Topological phases
NASA Astrophysics Data System (ADS)
Geraedts, Scott
years ago, it can still produce novel phenomena. Of much recent interest is the existence of non-Abelian anyons in FQHE systems. Though it is possible to construct wave functions that realize such particles, whether these wavefunctions are the ground state is a difficult quantitative question that must be answered numerically. In this thesis I describe progress using a density-matrix renormalization group algorithm to study a bilayer system thought to host non-Abelian anyons. We find phase diagrams in terms of experimentally relevant parameters, and also find evidence for a non-Abelian phase known as the 'interlayer Pfaffian'.
Razali, Azhani Mohd Abdullah, Jaafar
2015-04-29
Single Photon Emission Computed Tomography (SPECT) is a well-known imaging technique used in medical application, and it is part of medical imaging modalities that made the diagnosis and treatment of disease possible. However, SPECT technique is not only limited to the medical sector. Many works are carried out to adapt the same concept by using high-energy photon emission to diagnose process malfunctions in critical industrial systems such as in chemical reaction engineering research laboratories, as well as in oil and gas, petrochemical and petrochemical refining industries. Motivated by vast applications of SPECT technique, this work attempts to study the application of SPECT on a Pebble Bed Reactor (PBR) using numerical phantom of pebbles inside the PBR core. From the cross-sectional images obtained from SPECT, the behavior of pebbles inside the core can be analyzed for further improvement of the PBR design. As the quality of the reconstructed image is largely dependent on the algorithm used, this work aims to compare two image reconstruction algorithms for SPECT, namely the Expectation Maximization Algorithm and the Exact Inversion Formula. The results obtained from the Exact Inversion Formula showed better image contrast and sharpness, and shorter computational time compared to the Expectation Maximization Algorithm.
NASA Astrophysics Data System (ADS)
Razali, Azhani Mohd; Abdullah, Jaafar
2015-04-01
Single Photon Emission Computed Tomography (SPECT) is a well-known imaging technique used in medical application, and it is part of medical imaging modalities that made the diagnosis and treatment of disease possible. However, SPECT technique is not only limited to the medical sector. Many works are carried out to adapt the same concept by using high-energy photon emission to diagnose process malfunctions in critical industrial systems such as in chemical reaction engineering research laboratories, as well as in oil and gas, petrochemical and petrochemical refining industries. Motivated by vast applications of SPECT technique, this work attempts to study the application of SPECT on a Pebble Bed Reactor (PBR) using numerical phantom of pebbles inside the PBR core. From the cross-sectional images obtained from SPECT, the behavior of pebbles inside the core can be analyzed for further improvement of the PBR design. As the quality of the reconstructed image is largely dependent on the algorithm used, this work aims to compare two image reconstruction algorithms for SPECT, namely the Expectation Maximization Algorithm and the Exact Inversion Formula. The results obtained from the Exact Inversion Formula showed better image contrast and sharpness, and shorter computational time compared to the Expectation Maximization Algorithm.
Numerical Studies of Collisionless Current Layers
NASA Technical Reports Server (NTRS)
Quest, Kevin B.
1993-01-01
The purpose of this proposal was to investigate collisionless current layers using a variety of analytic and numerical tools. The first year of the contract was dedicated to analytical studies, to the porting and adaption of codes being used in this study, and to the numerical simulation of collisionless current layers. The second year entailed the development of multi-dimensional hybrid algorithms as well as the re-examination of the problem of integro-differential equations that occur in the linear stage of plasma instabilities.
A Parallel Compact Multi-Dimensional Numerical Algorithm with Aeroacoustics Applications
NASA Technical Reports Server (NTRS)
Povitsky, Alex; Morris, Philip J.
1999-01-01
In this study we propose a novel method to parallelize high-order compact numerical algorithms for the solution of three-dimensional PDEs (Partial Differential Equations) in a space-time domain. For this numerical integration most of the computer time is spent in computation of spatial derivatives at each stage of the Runge-Kutta temporal update. The most efficient direct method to compute spatial derivatives on a serial computer is a version of Gaussian elimination for narrow linear banded systems known as the Thomas algorithm. In a straightforward pipelined implementation of the Thomas algorithm processors are idle due to the forward and backward recurrences of the Thomas algorithm. To utilize processors during this time, we propose to use them for either non-local data independent computations, solving lines in the next spatial direction, or local data-dependent computations by the Runge-Kutta method. To achieve this goal, control of processor communication and computations by a static schedule is adopted. Thus, our parallel code is driven by a communication and computation schedule instead of the usual "creative, programming" approach. The obtained parallelization speed-up of the novel algorithm is about twice as much as that for the standard pipelined algorithm and close to that for the explicit DRP algorithm.
Ozone Uncertainties Study Algorithm (OUSA)
NASA Technical Reports Server (NTRS)
Bahethi, O. P.
1982-01-01
An algorithm to carry out sensitivities, uncertainties and overall imprecision studies to a set of input parameters for a one dimensional steady ozone photochemistry model is described. This algorithm can be used to evaluate steady state perturbations due to point source or distributed ejection of H2O, CLX, and NOx, besides, varying the incident solar flux. This algorithm is operational on IBM OS/360-91 computer at NASA/Goddard Space Flight Center's Science and Applications Computer Center (SACC).
Ozone Uncertainties Study Algorithm (OUSA)
NASA Astrophysics Data System (ADS)
Bahethi, O. P.
An algorithm to carry out sensitivities, uncertainties and overall imprecision studies to a set of input parameters for a one dimensional steady ozone photochemistry model is described. This algorithm can be used to evaluate steady state perturbations due to point source or distributed ejection of H2O, CLX, and NOx, besides, varying the incident solar flux. This algorithm is operational on IBM OS/360-91 computer at NASA/Goddard Space Flight Center's Science and Applications Computer Center (SACC).
Case study of isosurface extraction algorithm performance
Sutton, P M; Hansen, C D; Shen, H; Schikore, D
1999-12-14
Isosurface extraction is an important and useful visualization method. Over the past ten years, the field has seen numerous isosurface techniques published leaving the user in a quandary about which one should be used. Some papers have published complexity analysis of the techniques yet empirical evidence comparing different methods is lacking. This case study presents a comparative study of several representative isosurface extraction algorithms. It reports and analyzes empirical measurements of execution times and memory behavior for each algorithm. The results show that asymptotically optimal techniques may not be the best choice when implemented on modern computer architectures.
A fast algorithm for numerical solutions to Fortet's equation
NASA Astrophysics Data System (ADS)
Brumen, Gorazd
2008-10-01
A fast algorithm for computation of default times of multiple firms in a structural model is presented. The algorithm uses a multivariate extension of the Fortet's equation and the structure of Toeplitz matrices to significantly improve the computation time. In a financial market consisting of M[not double greater-than sign]1 firms and N discretization points in every dimension the algorithm uses O(nlogn·M·M!·NM(M-1)/2) operations, where n is the number of discretization points in the time domain. The algorithm is applied to firm survival probability computation and zero coupon bond pricing.
François, Marianne M.
2015-05-28
A review of recent advances made in numerical methods and algorithms within the volume tracking framework is presented. The volume tracking method, also known as the volume-of-fluid method has become an established numerical approach to model and simulate interfacial flows. Its advantage is its strict mass conservation. However, because the interface is not explicitly tracked but captured via the material volume fraction on a fixed mesh, accurate estimation of the interface position, its geometric properties and modeling of interfacial physics in the volume tracking framework remain difficult. Several improvements have been made over the last decade to address these challenges. In this study, the multimaterial interface reconstruction method via power diagram, curvature estimation via heights and mean values and the balanced-force algorithm for surface tension are highlighted.
François, Marianne M.
2015-05-28
A review of recent advances made in numerical methods and algorithms within the volume tracking framework is presented. The volume tracking method, also known as the volume-of-fluid method has become an established numerical approach to model and simulate interfacial flows. Its advantage is its strict mass conservation. However, because the interface is not explicitly tracked but captured via the material volume fraction on a fixed mesh, accurate estimation of the interface position, its geometric properties and modeling of interfacial physics in the volume tracking framework remain difficult. Several improvements have been made over the last decade to address these challenges.more » In this study, the multimaterial interface reconstruction method via power diagram, curvature estimation via heights and mean values and the balanced-force algorithm for surface tension are highlighted.« less
Stochastic algorithms for the analysis of numerical flame simulations
Bell, John B.; Day, Marcus S.; Grcar, Joseph F.; Lijewski, Michael J.
2004-04-26
Recent progress in simulation methodologies and high-performance parallel computers have made it is possible to perform detailed simulations of multidimensional reacting flow phenomena using comprehensive kinetics mechanisms. As simulations become larger and more complex, it becomes increasingly difficult to extract useful information from the numerical solution, particularly regarding the interactions of the chemical reaction and diffusion processes. In this paper we present a new diagnostic tool for analysis of numerical simulations of reacting flow. Our approach is based on recasting an Eulerian flow solution in a Lagrangian frame. Unlike a conventional Lagrangian view point that follows the evolution of a volume of the fluid, we instead follow specific chemical elements, e.g., carbon, nitrogen, etc., as they move through the system . From this perspective an ''atom'' is part of some molecule of a species that is transported through the domain by advection and diffusion. Reactions cause the atom to shift from one chemical host species to another and the subsequent transport of the atom is given by the movement of the new species. We represent these processes using a stochastic particle formulation that treats advection deterministically and models diffusion and chemistry as stochastic processes. In this paper, we discuss the numerical issues in detail and demonstrate that an ensemble of stochastic trajectories can accurately capture key features of the continuum solution. The capabilities of this diagnostic are then demonstrated by applications to study the modulation of carbon chemistry during a vortex-flame interaction, and the role of cyano chemistry in rm NO{sub x} production for a steady diffusion flame.
Numerical analysis of EPR spectra. 7. The simplex algorithm
NASA Astrophysics Data System (ADS)
Beckwith, Athelstan L. J.; Brumby, Steven
The Simplex algorithm is well suited to the least-squares analysis of highly complex EPR spectra. The application of the algorithm to the analysis of the spectra of benzo[ a]pyrenyl-6-oxy, chloro(methoxycarbonyl)methyl, and cyano(methoxy)methyl free radicals is described.
A comparative study of AGN feedback algorithms
NASA Astrophysics Data System (ADS)
Wurster, J.; Thacker, R. J.
2013-05-01
Modelling active galactic nuclei (AGN) feedback in numerical simulations is both technically and theoretically challenging, with numerous approaches having been published in the literature. We present a study of five distinct approaches to modelling AGN feedback within gravitohydrodynamic simulations of major mergers of Milky Way-sized galaxies. To constrain differences to only be between AGN feedback models, all simulations start from the same initial conditions and use the same star formation algorithm. Most AGN feedback algorithms have five key aspects: the black hole accretion rate, energy feedback rate and method, particle accretion algorithm, black hole advection algorithm and black hole merger algorithm. All models follow different accretion histories, and in some cases, accretion rates differ by up to three orders of magnitude at any given time. We consider models with either thermal or kinetic feedback, with the associated energy deposited locally around the black hole. Each feedback algorithm modifies the region around the black hole to different extents, yielding gas densities and temperatures within r ˜ 200 pc that differ by up to six orders of magnitude at any given time. The particle accretion algorithms usually maintain good agreement between the total mass accreted by dot{M} dt and the total mass of gas particles removed from the simulation, although not all algorithms guarantee this to be true. The black hole advection algorithms dampen inappropriate dragging of the black holes by two-body interactions. Advecting the black hole a limited distance based upon local mass distributions has many desirably properties, such as avoiding large artificial jumps and allowing the possibility of the black hole remaining in a gas void. Lastly, two black holes instantly merge when given criteria are met, and we find a range of merger times for different criteria. This is important since the AGN feedback rate changes across the merger in a way that is dependent
Particle-In-Cell Multi-Algorithm Numerical Test-Bed
NASA Astrophysics Data System (ADS)
Meyers, M. D.; Yu, P.; Tableman, A.; Decyk, V. K.; Mori, W. B.
2015-11-01
We describe a numerical test-bed that allows for the direct comparison of different numerical simulation schemes using only a single code. It is built from the UPIC Framework, which is a set of codes and modules for constructing parallel PIC codes. In this test-bed code, Maxwell's equations are solved in Fourier space in two dimensions. One can readily examine the numerical properties of a real space finite difference scheme by including its operators' Fourier space representations in the Maxwell solver. The fields can be defined at the same location in a simulation cell or can be offset appropriately by half-cells, as in the Yee finite difference time domain scheme. This allows for the accurate comparison of numerical properties (dispersion relations, numerical stability, etc.) across finite difference schemes, or against the original spectral scheme. We have also included different options for the charge and current deposits, including a strict charge conserving current deposit. The test-bed also includes options for studying the analytic time domain scheme, which eliminates numerical dispersion errors in vacuum. We will show examples from the test-bed that illustrate how the properties of some numerical instabilities vary between different PIC algorithms. Work supported by the NSF grant ACI 1339893 and DOE grant DE-SC0008491.
NASA Astrophysics Data System (ADS)
Alfonso, Lester; Zamora, Jose; Cruz, Pedro
2015-04-01
The stochastic approach to coagulation considers the coalescence process going in a system of a finite number of particles enclosed in a finite volume. Within this approach, the full description of the system can be obtained from the solution of the multivariate master equation, which models the evolution of the probability distribution of the state vector for the number of particles of a given mass. Unfortunately, due to its complexity, only limited results were obtained for certain type of kernels and monodisperse initial conditions. In this work, a novel numerical algorithm for the solution of the multivariate master equation for stochastic coalescence that works for any type of kernels and initial conditions is introduced. The performance of the method was checked by comparing the numerically calculated particle mass spectrum with analytical solutions obtained for the constant and sum kernels, with an excellent correspondence between the analytical and numerical solutions. In order to increase the speedup of the algorithm, software parallelization techniques with OpenMP standard were used, along with an implementation in order to take advantage of new accelerator technologies. Simulations results show an important speedup of the parallelized algorithms. This study was funded by a grant from Consejo Nacional de Ciencia y Tecnologia de Mexico SEP-CONACYT CB-131879. The authors also thanks LUFAC® Computacion SA de CV for CPU time and all the support provided.
Numerical algorithms for steady and unsteady incompressible Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Hafez, Mohammed; Dacles, Jennifer
1989-01-01
The numerical analysis of the incompressible Navier-Stokes equations are becoming important tools in the understanding of some fluid flow problems which are encountered in research as well as in industry. With the advent of the supercomputers, more realistic problems can be studied with a wider choice of numerical algorithms. An alternative formulation is presented for viscous incompressible flows. The incompressible Navier-Stokes equations are cast in a velocity/vorticity formulation. This formulation consists of solving the Poisson equations for the velocity components and the vorticity transport equation. Two numerical algorithms for the steady two-dimensional laminar flows are presented. The first method is based on the actual partial differential equations. This uses a finite-difference approximation of the governing equations on a staggered grid. The second method uses a finite element discretization with the vorticity transport equation approximated using a Galerkin approximation and the Poisson equations are obtained using a least squares method. The equations are solved efficiently using Newton's method and a banded direct matrix solver (LINPACK). The method is extended to steady three-dimensional laminar flows and applied to a cubic driven cavity using finite difference schemes and a staggered grid arrangement on a Cartesian mesh. The equations are solved iteratively using a plane zebra relaxation scheme. Currently, a two-dimensional, unsteady algorithm is being developed using a generalized coordinate system. The equations are discretized using a finite-volume approach. This work will then be extended to three-dimensional flows.
Numerical Optimization Algorithms and Software for Systems Biology
Saunders, Michael
2013-02-02
The basic aims of this work are: to develop reliable algorithms for solving optimization problems involving large stoi- chiometric matrices; to investigate cyclic dependency between metabolic and macromolecular biosynthetic networks; and to quantify the significance of thermodynamic constraints on prokaryotic metabolism.
Numerical study of rock blasting
NASA Astrophysics Data System (ADS)
Stefanov, Yu. P.; Bakeev, R. A.; Yudin, A. S.; Kuznetsova, N. S.
2015-10-01
The paper presents numerical simulation results on fracture of a concrete block due to dynamic explosive loads applied to the walls of a blast hole. Considered in the study is the influence of the pulse shape and rock properties on the pattern of irreversible deformation and cracking. It is found that a fractured zone bounded by a plastically deformed contour always arises around the explosion site. Comparison of elastoplastic deformation and fracture induced in the concrete block by explosion pulses of different durations and amplitudes shows that shorter pulses with higher amplitudes and steeper rise times provide a higher blasting efficiency.
A Numerical Study of Boson Star Binaries
NASA Astrophysics Data System (ADS)
Mundim, Bruno C.
2010-02-01
This thesis describes a numerical study of binary boson stars within the context of an approximation to general relativity. The approximation we adopt places certain restrictions on the dynamical variables of general relativity (conformal flatness of the 3-metric), and on the time-slicing of the spacetime (maximal slicing). The resulting modeling problem requires the solution of a coupled nonlinear system of 4 hyperbolic, and 5 elliptic partial differential equations (PDEs) in three space dimensions and time. We approximately solve this system as an initial-boundary value problem, using finite difference techniques and well known, computationally efficient numerical algorithms such as the multigrid method in the case of the elliptic equations. Careful attention is paid to the issue of code validation, and a key part of the thesis is the demonstration that, as the basic scale of finite difference discretization is reduced, our numerical code generates results that converge to a solution of the continuum system of PDEs as desired. The thesis concludes with a discussion of results from some initial explorations of the orbital dynamics of boson star binaries. In particular, we describe calculations in which motion of such a binary is followed for more than two orbital periods, which is a significant advance over previous studies. We also present results from computations in which the boson stars merge, and where there is evidence for black hole formation.
Computational Fluid Dynamics. [numerical methods and algorithm development
NASA Technical Reports Server (NTRS)
1992-01-01
This collection of papers was presented at the Computational Fluid Dynamics (CFD) Conference held at Ames Research Center in California on March 12 through 14, 1991. It is an overview of CFD activities at NASA Lewis Research Center. The main thrust of computational work at Lewis is aimed at propulsion systems. Specific issues related to propulsion CFD and associated modeling will also be presented. Examples of results obtained with the most recent algorithm development will also be presented.
A bibliography on parallel and vector numerical algorithms
NASA Technical Reports Server (NTRS)
Ortega, James M.; Voigt, Robert G.; Romine, Charles H.
1988-01-01
This is a bibliography on numerical methods. It also includes a number of other references on machine architecture, programming language, and other topics of interest to scientific computing. Certain conference proceedings and anthologies which have been published in book form are also listed.
A bibliography on parallel and vector numerical algorithms
NASA Technical Reports Server (NTRS)
Ortega, J. M.; Voigt, R. G.
1987-01-01
This is a bibliography of numerical methods. It also includes a number of other references on machine architecture, programming language, and other topics of interest to scientific computing. Certain conference proceedings and anthologies which have been published in book form are listed also.
Numerical Laplace Transform Inversion Employing the Gaver-Stehfest Algorithm.
ERIC Educational Resources Information Center
Jacquot, Raymond G.; And Others
1985-01-01
Presents a technique for the numerical inversion of Laplace Transforms and several examples employing this technique. Limitations of the method in terms of available computer word length and the effects of these limitations on approximate inverse functions are also discussed. (JN)
Numerical stability analysis of the pseudo-spectral analytical time-domain PIC algorithm
Godfrey, Brendan B.; Vay, Jean-Luc; Haber, Irving
2014-02-01
The pseudo-spectral analytical time-domain (PSATD) particle-in-cell (PIC) algorithm solves the vacuum Maxwell's equations exactly, has no Courant time-step limit (as conventionally defined), and offers substantial flexibility in plasma and particle beam simulations. It is, however, not free of the usual numerical instabilities, including the numerical Cherenkov instability, when applied to relativistic beam simulations. This paper derives and solves the numerical dispersion relation for the PSATD algorithm and compares the results with corresponding behavior of the more conventional pseudo-spectral time-domain (PSTD) and finite difference time-domain (FDTD) algorithms. In general, PSATD offers superior stability properties over a reasonable range of time steps. More importantly, one version of the PSATD algorithm, when combined with digital filtering, is almost completely free of the numerical Cherenkov instability for time steps (scaled to the speed of light) comparable to or smaller than the axial cell size.
Numerical studies of frontal dynamics
NASA Technical Reports Server (NTRS)
Keyser, Daniel
1986-01-01
Efforts concentrated on the development of a two dimensional primitive equation (PE) model of frontogenesis that simultaneously incorporates the frontagenetical mechanisms of confluence and horizontal shear. Applying this model to study the effects of upper level frontogenesis, it appeared to be dominated by tilting effects associated with cross front variation of vertical motion, in which subsidence is maximized within and to the warm side of the frontal zone. Results suggest that aspects characteristic of three-dimensional baroclinic waves may be abstracted to a significant extent in a two dimensional framework. They also show that upper-level frontogenesis and tropopause folding can occur in the absence of three-dimensional curvature effects, commonly believed to be necessary for realistic upper-level frontogenesis. An implication of the dominant tilting effects is that they may have to be adequately resolved by numerical weather prediction models, thus requiring better horizontal and vertical resolution.
Shock focusing flow field simulated by a high-resolution numerical algorithm
NASA Astrophysics Data System (ADS)
Jung, Y. G.; Chang, K. S.
2012-11-01
Shock-focusing concave reflector is a very simple and effective tool to obtain a high-pressure pulse wave near the physical focal point. In the past, many optical images were obtained through experimental studies. However, measurement of field variables is not easy because the phenomenon is of short duration and the magnitude of shock waves is varied from pulse to pulse due to poor reproducibility. Using a wave propagation algorithm and the Cartesian embedded boundary method, we have successfully obtained numerical schlieren images that resemble the experimental results. By the numerical results, various field variables, such as pressure, density and vorticity, become available for the better understanding and design of shock focusing devices.
Thrombosis modeling in intracranial aneurysms: a lattice Boltzmann numerical algorithm
NASA Astrophysics Data System (ADS)
Ouared, R.; Chopard, B.; Stahl, B.; Rüfenacht, D. A.; Yilmaz, H.; Courbebaisse, G.
2008-07-01
The lattice Boltzmann numerical method is applied to model blood flow (plasma and platelets) and clotting in intracranial aneurysms at a mesoscopic level. The dynamics of blood clotting (thrombosis) is governed by mechanical variations of shear stress near wall that influence platelets-wall interactions. Thrombosis starts and grows below a shear rate threshold, and stops above it. Within this assumption, it is possible to account qualitatively well for partial, full or no occlusion of the aneurysm, and to explain why spontaneous thrombosis is more likely to occur in giant aneurysms than in small or medium sized aneurysms.
A Numerical Algorithm for the Solution of a Phase-Field Model of Polycrystalline Materials
Dorr, M R; Fattebert, J; Wickett, M E; Belak, J F; Turchi, P A
2008-12-04
We describe an algorithm for the numerical solution of a phase-field model (PFM) of microstructure evolution in polycrystalline materials. The PFM system of equations includes a local order parameter, a quaternion representation of local orientation and a species composition parameter. The algorithm is based on the implicit integration of a semidiscretization of the PFM system using a backward difference formula (BDF) temporal discretization combined with a Newton-Krylov algorithm to solve the nonlinear system at each time step. The BDF algorithm is combined with a coordinate projection method to maintain quaternion unit length, which is related to an important solution invariant. A key element of the Newton-Krylov algorithm is the selection of a preconditioner to accelerate the convergence of the Generalized Minimum Residual algorithm used to solve the Jacobian linear system in each Newton step. Results are presented for the application of the algorithm to 2D and 3D examples.
Numerical Algorithms for Precise and Efficient Orbit Propagation and Positioning
NASA Astrophysics Data System (ADS)
Bradley, Ben K.
Motivated by the growing space catalog and the demands for precise orbit determination with shorter latency for science and reconnaissance missions, this research improves the computational performance of orbit propagation through more efficient and precise numerical integration and frame transformation implementations. Propagation of satellite orbits is required for astrodynamics applications including mission design, orbit determination in support of operations and payload data analysis, and conjunction assessment. Each of these applications has somewhat different requirements in terms of accuracy, precision, latency, and computational load. This dissertation develops procedures to achieve various levels of accuracy while minimizing computational cost for diverse orbit determination applications. This is done by addressing two aspects of orbit determination: (1) numerical integration used for orbit propagation and (2) precise frame transformations necessary for force model evaluation and station coordinate rotations. This dissertation describes a recently developed method for numerical integration, dubbed Bandlimited Collocation Implicit Runge-Kutta (BLC-IRK), and compare its efficiency in propagating orbits to existing techniques commonly used in astrodynamics. The BLC-IRK scheme uses generalized Gaussian quadratures for bandlimited functions. It requires significantly fewer force function evaluations than explicit Runge-Kutta schemes and approaches the efficiency of the 8th-order Gauss-Jackson multistep method. Converting between the Geocentric Celestial Reference System (GCRS) and International Terrestrial Reference System (ITRS) is necessary for many applications in astrodynamics, such as orbit propagation, orbit determination, and analyzing geoscience data from satellite missions. This dissertation provides simplifications to the Celestial Intermediate Origin (CIO) transformation scheme and Earth orientation parameter (EOP) storage for use in positioning and
NASA Astrophysics Data System (ADS)
Dong, Suchuan
2015-11-01
This talk focuses on simulating the motion of a mixture of N (N>=2) immiscible incompressible fluids with given densities, dynamic viscosities and pairwise surface tensions. We present an N-phase formulation within the phase field framework that is thermodynamically consistent, in the sense that the formulation satisfies the conservations of mass/momentum, the second law of thermodynamics and Galilean invariance. We also present an efficient algorithm for numerically simulating the N-phase system. The algorithm has overcome the issues caused by the variable coefficient matrices associated with the variable mixture density/viscosity and the couplings among the (N-1) phase field variables and the flow variables. We compare simulation results with the Langmuir-de Gennes theory to demonstrate that the presented method produces physically accurate results for multiple fluid phases. Numerical experiments will be presented for several problems involving multiple fluid phases, large density contrasts and large viscosity contrasts to demonstrate the capabilities of the method for studying the interactions among multiple types of fluid interfaces. Support from NSF and ONR is gratefully acknowledged.
A stable and efficient numerical algorithm for unconfined aquifer analysis.
Keating, Elizabeth; Zyvoloski, George
2009-01-01
The nonlinearity of equations governing flow in unconfined aquifers poses challenges for numerical models, particularly in field-scale applications. Existing methods are often unstable, do not converge, or require extremely fine grids and small time steps. Standard modeling procedures such as automated model calibration and Monte Carlo uncertainty analysis typically require thousands of model runs. Stable and efficient model performance is essential to these analyses. We propose a new method that offers improvements in stability and efficiency and is relatively tolerant of coarse grids. It applies a strategy similar to that in the MODFLOW code to the solution of Richard's equation with a grid-dependent pressure/saturation relationship. The method imposes a contrast between horizontal and vertical permeability in gridblocks containing the water table, does not require "dry" cells to convert to inactive cells, and allows recharge to flow through relatively dry cells to the water table. We establish the accuracy of the method by comparison to an analytical solution for radial flow to a well in an unconfined aquifer with delayed yield. Using a suite of test problems, we demonstrate the efficiencies gained in speed and accuracy over two-phase simulations, and improved stability when compared to MODFLOW. The advantages for applications to transient unconfined aquifer analysis are clearly demonstrated by our examples. We also demonstrate applicability to mixed vadose zone/saturated zone applications, including transport, and find that the method shows great promise for these types of problem as well. PMID:19341374
A stable and efficient numerical algorithm for unconfined aquifer analysis
Keating, Elizabeth; Zyvoloski, George
2008-01-01
The non-linearity of equations governing flow in unconfined aquifers poses challenges for numerical models, particularly in field-scale applications. Existing methods are often unstable, do not converge, or require extremely fine grids and small time steps. Standard modeling procedures such as automated model calibration and Monte Carlo uncertainty analysis typically require thousands of forward model runs. Stable and efficient model performance is essential to these analyses. We propose a new method that offers improvements in stability and efficiency, and is relatively tolerant of coarse grids. It applies a strategy similar to that in the MODFLOW code to solution of Richard's Equation with a grid-dependent pressure/saturation relationship. The method imposes a contrast between horizontal and vertical permeability in gridblocks containing the water table. We establish the accuracy of the method by comparison to an analytical solution for radial flow to a well in an unconfined aquifer with delayed yield. Using a suite of test problems, we demonstrate the efficiencies gained in speed and accuracy over two-phase simulations, and improved stability when compared to MODFLOW. The advantages for applications to transient unconfined aquifer analysis are clearly demonstrated by our examples. We also demonstrate applicability to mixed vadose zone/saturated zone applications, including transport, and find that the method shows great promise for these types of problem, as well.
Comparison of Fully Numerical Predictor-Corrector and Apollo Skip Entry Guidance Algorithms
NASA Astrophysics Data System (ADS)
Brunner, Christopher W.; Lu, Ping
2012-09-01
The dramatic increase in computational power since the Apollo program has enabled the development of numerical predictor-corrector (NPC) entry guidance algorithms that allow on-board accurate determination of a vehicle's trajectory. These algorithms are sufficiently mature to be flown. They are highly adaptive, especially in the face of extreme dispersion and off-nominal situations compared with reference-trajectory following algorithms. The performance and reliability of entry guidance are critical to mission success. This paper compares the performance of a recently developed fully numerical predictor-corrector entry guidance (FNPEG) algorithm with that of the Apollo skip entry guidance. Through extensive dispersion testing, it is clearly demonstrated that the Apollo skip entry guidance algorithm would be inadequate in meeting the landing precision requirement for missions with medium (4000-7000 km) and long (>7000 km) downrange capability requirements under moderate dispersions chiefly due to poor modeling of atmospheric drag. In the presence of large dispersions, a significant number of failures occur even for short-range missions due to the deviation from planned reference trajectories. The FNPEG algorithm, on the other hand, is able to ensure high landing precision in all cases tested. All factors considered, a strong case is made for adopting fully numerical algorithms for future skip entry missions.
Variationally consistent discretization schemes and numerical algorithms for contact problems
NASA Astrophysics Data System (ADS)
Wohlmuth, Barbara
We consider variationally consistent discretization schemes for mechanical contact problems. Most of the results can also be applied to other variational inequalities, such as those for phase transition problems in porous media, for plasticity or for option pricing applications from finance. The starting point is to weakly incorporate the constraint into the setting and to reformulate the inequality in the displacement in terms of a saddle-point problem. Here, the Lagrange multiplier represents the surface forces, and the constraints are restricted to the boundary of the simulation domain. Having a uniform inf-sup bound, one can then establish optimal low-order a priori convergence rates for the discretization error in the primal and dual variables. In addition to the abstract framework of linear saddle-point theory, complementarity terms have to be taken into account. The resulting inequality system is solved by rewriting it equivalently by means of the non-linear complementarity function as a system of equations. Although it is not differentiable in the classical sense, semi-smooth Newton methods, yielding super-linear convergence rates, can be applied and easily implemented in terms of a primal-dual active set strategy. Quite often the solution of contact problems has a low regularity, and the efficiency of the approach can be improved by using adaptive refinement techniques. Different standard types, such as residual- and equilibrated-based a posteriori error estimators, can be designed based on the interpretation of the dual variable as Neumann boundary condition. For the fully dynamic setting it is of interest to apply energy-preserving time-integration schemes. However, the differential algebraic character of the system can result in high oscillations if standard methods are applied. A possible remedy is to modify the fully discretized system by a local redistribution of the mass. Numerical results in two and three dimensions illustrate the wide range of
Numerical optimization algorithm for rotationally invariant multi-orbital slave-boson method
NASA Astrophysics Data System (ADS)
Quan, Ya-Min; Wang, Qing-wei; Liu, Da-Yong; Yu, Xiang-Long; Zou, Liang-Jian
2015-06-01
We develop a generalized numerical optimization algorithm for the rotationally invariant multi-orbital slave boson approach, which is applicable for arbitrary boundary constraints of high-dimensional objective function by combining several classical optimization techniques. After constructing the calculation architecture of rotationally invariant multi-orbital slave boson model, we apply this optimization algorithm to find the stable ground state and magnetic configuration of two-orbital Hubbard models. The numerical results are consistent with available solutions, confirming the correctness and accuracy of our present algorithm. Furthermore, we utilize it to explore the effects of the transverse Hund's coupling terms on metal-insulator transition, orbital selective Mott phase and magnetism. These results show the quick convergency and robust stable character of our algorithm in searching the optimized solution of strongly correlated electron systems.
NASA Astrophysics Data System (ADS)
Balashov, V. A.; Savenkov, E. B.
2015-10-01
The applicability of numerical algorithms based on a quasi-hydrodynamic system of equations for computing viscous heat-conducting compressible gas flows at Mach numbers M = 10-2-10-1 is studied numerically. The numerical algorithm is briefly described, and the results obtained for a number of two- and three-dimensional test problems are presented and compared with earlier numerical data.
Extremal polynomials and methods of optimization of numerical algorithms
Lebedev, V I
2004-10-31
Chebyshev-Markov-Bernstein-Szegoe polynomials C{sub n}(x) extremal on [-1,1] with weight functions w(x)=(1+x){sup {alpha}}(1- x){sup {beta}}/{radical}(S{sub l}(x)) where {alpha},{beta}=0,1/2 and S{sub l}(x)={pi}{sub k=1}{sup m}(1-c{sub k}T{sub l{sub k}}(x))>0 are considered. A universal formula for their representation in trigonometric form is presented. Optimal distributions of the nodes of the weighted interpolation and explicit quadrature formulae of Gauss, Markov, Lobatto, and Rado types are obtained for integrals with weight p(x)=w{sup 2}(x)(1-x{sup 2}){sup -1/2}. The parameters of optimal Chebyshev iterative methods reducing the error optimally by comparison with the initial error defined in another norm are determined. For each stage of the Fedorenko-Bakhvalov method iteration parameters are determined which take account of the results of the previous calculations. Chebyshev filters with weight are constructed. Iterative methods of the solution of equations containing compact operators are studied.
Extremal polynomials and methods of optimization of numerical algorithms
NASA Astrophysics Data System (ADS)
Lebedev, V. I.
2004-10-01
Chebyshëv-Markov-Bernstein-Szegö polynomials C_n(x) extremal on \\lbrack -1,1 \\rbrack with weight functions w(x)=(1+x)^\\alpha(1- x)^\\beta/\\sqrt{S_l(x)} where \\alpha,\\beta=0,\\frac12 and S_l(x)=\\prod_{k=1}^m(1-c_kT_{l_k}(x))>0 are considered. A universal formula for their representation in trigonometric form is presented. Optimal distributions of the nodes of the weighted interpolation and explicit quadrature formulae of Gauss, Markov, Lobatto, and Rado types are obtained for integrals with weight p(x)=w^2(x)(1-x^2)^{-1/2}. The parameters of optimal Chebyshëv iterative methods reducing the error optimally by comparison with the initial error defined in another norm are determined. For each stage of the Fedorenko-Bakhvalov method iteration parameters are determined which take account of the results of the previous calculations. Chebyshëv filters with weight are constructed. Iterative methods of the solution of equations containing compact operators are studied.
Wang, Peng; Zhu, Zhouquan; Huang, Shuai
2013-01-01
This paper presents a novel biologically inspired metaheuristic algorithm called seven-spot ladybird optimization (SLO). The SLO is inspired by recent discoveries on the foraging behavior of a seven-spot ladybird. In this paper, the performance of the SLO is compared with that of the genetic algorithm, particle swarm optimization, and artificial bee colony algorithms by using five numerical benchmark functions with multimodality. The results show that SLO has the ability to find the best solution with a comparatively small population size and is suitable for solving optimization problems with lower dimensions. PMID:24385879
On the impact of communication complexity in the design of parallel numerical algorithms
NASA Technical Reports Server (NTRS)
Gannon, D.; Vanrosendale, J.
1984-01-01
This paper describes two models of the cost of data movement in parallel numerical algorithms. One model is a generalization of an approach due to Hockney, and is suitable for shared memory multiprocessors where each processor has vector capabilities. The other model is applicable to highly parallel nonshared memory MIMD systems. In the second model, algorithm performance is characterized in terms of the communication network design. Techniques used in VLSI complexity theory are also brought in, and algorithm independent upper bounds on system performance are derived for several problems that are important to scientific computation.
Zhu, Zhouquan
2013-01-01
This paper presents a novel biologically inspired metaheuristic algorithm called seven-spot ladybird optimization (SLO). The SLO is inspired by recent discoveries on the foraging behavior of a seven-spot ladybird. In this paper, the performance of the SLO is compared with that of the genetic algorithm, particle swarm optimization, and artificial bee colony algorithms by using five numerical benchmark functions with multimodality. The results show that SLO has the ability to find the best solution with a comparatively small population size and is suitable for solving optimization problems with lower dimensions. PMID:24385879
NASA Astrophysics Data System (ADS)
Simpson, Matthew J.; Landman, Kerry A.; Newgreen, Donald F.
2006-08-01
A numerical algorithm to simulate chemotactic and/or diffusive migration on a one-dimensional growing domain is developed. The domain growth can be spatially nonuniform and the growth-derived advection term must be discretised. The hyperbolic terms in the conservation equations associated with chemotactic migration and domain growth are accurately discretised using an explicit central scheme. Generality of the algorithm is maintained using an operator split technique to simulate diffusive migration implicitly. The resulting algorithm is applicable for any combination of diffusive and/or chemotactic migration on a growing domain with a general growth-induced velocity field. The accuracy of the algorithm is demonstrated by testing the results against some simple analytical solutions and in an inter-code comparison. The new algorithm demonstrates that the form of nonuniform growth plays a critical role in determining whether a population of migratory cells is able to overcome the domain growth and fully colonise the domain.
Burtsev, S.; Camassa, R.; Timofeyev, I.
1998-11-20
The authors implement two different algorithms for computing numerically the direct Zakharov-Shabat eigenvalue problem on the infinite line. The first algorithm replaces the potential in the eigenvalue problem by a piecewise-constant approximation, which allows one to solve analytically the corresponding ordinary differential equation. The resulting algorithm is of second order in the step size. The second algorithm uses the fourth-order Runge-Kutta method. They test and compare the performance of these two algorithms on three exactly solvable potentials. They find that even though the Runge-Kutta method is of higher order, this extra accuracy can be lost because of the additional dependence of its numerical error on the eigenvalue. this limits the usefulness of the Runge-Kutta algorithm to a region inside the unit circle around the origin in the complex plane of the eigenvalues. For the computation of the continuous spectrum density, this limitation is particularly severe, as revealed by the spectral decomposition of the L{sup 2}-norm of a solution to the nonlinear Schroedinger equation. They show that no such limitations exist for the piecewise-constant algorithm. In particular, this scheme converges uniformly for both continuous and discrete spectrum components.
Gyrotactic trapping: A numerical study
NASA Astrophysics Data System (ADS)
Ghorai, S.
2016-04-01
Gyrotactic trapping is a mechanism proposed by Durham et al. ["Disruption of vertical motility by shear triggers formation of thin Phytoplankton layers," Science 323, 1067-1070 (2009)] to explain the formation of thin phytoplankton layer just below the ocean surface. This mechanism is examined numerically using a rational model based on the generalized Taylor dispersion theory. The crucial role of sedimentation speed in the thin layer formation is demonstrated. The effects of variation in different parameters on the thin layer formation are also investigated.
A new free-surface stabilization algorithm for geodynamical modelling: Theory and numerical tests
NASA Astrophysics Data System (ADS)
Andrés-Martínez, Miguel; Morgan, Jason P.; Pérez-Gussinyé, Marta; Rüpke, Lars
2015-09-01
The surface of the solid Earth is effectively stress free in its subaerial portions, and hydrostatic beneath the oceans. Unfortunately, this type of boundary condition is difficult to treat computationally, and for computational convenience, numerical models have often used simpler approximations that do not involve a normal stress-loaded, shear-stress free top surface that is free to move. Viscous flow models with a computational free surface typically confront stability problems when the time step is bigger than the viscous relaxation time. The small time step required for stability (< 2 Kyr) makes this type of model computationally intensive, so there remains a need to develop strategies that mitigate the stability problem by making larger (at least ∼10 Kyr) time steps stable and accurate. Here we present a new free-surface stabilization algorithm for finite element codes which solves the stability problem by adding to the Stokes formulation an intrinsic penalization term equivalent to a portion of the future load at the surface nodes. Our algorithm is straightforward to implement and can be used with both Eulerian or Lagrangian grids. It includes α and β parameters to respectively control both the vertical and the horizontal slope-dependent penalization terms, and uses Uzawa-like iterations to solve the resulting system at a cost comparable to a non-stress free surface formulation. Four tests were carried out in order to study the accuracy and the stability of the algorithm: (1) a decaying first-order sinusoidal topography test, (2) a decaying high-order sinusoidal topography test, (3) a Rayleigh-Taylor instability test, and (4) a steep-slope test. For these tests, we investigate which α and β parameters give the best results in terms of both accuracy and stability. We also compare the accuracy and the stability of our algorithm with a similar implicit approach recently developed by Kaus et al. (2010). We find that our algorithm is slightly more accurate
NASA Technical Reports Server (NTRS)
Yee, H. C.; Sweby, P. K.; Griffiths, D. F.
1991-01-01
Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.
NASA Technical Reports Server (NTRS)
Yee, H. C.; Sweby, P. K.; Griffiths, D. F.
1990-01-01
Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.
PolyPole-1: An accurate numerical algorithm for intra-granular fission gas release
NASA Astrophysics Data System (ADS)
Pizzocri, D.; Rabiti, C.; Luzzi, L.; Barani, T.; Van Uffelen, P.; Pastore, G.
2016-09-01
The transport of fission gas from within the fuel grains to the grain boundaries (intra-granular fission gas release) is a fundamental controlling mechanism of fission gas release and gaseous swelling in nuclear fuel. Hence, accurate numerical solution of the corresponding mathematical problem needs to be included in fission gas behaviour models used in fuel performance codes. Under the assumption of equilibrium between trapping and resolution, the process can be described mathematically by a single diffusion equation for the gas atom concentration in a grain. In this paper, we propose a new numerical algorithm (PolyPole-1) to efficiently solve the fission gas diffusion equation in time-varying conditions. The PolyPole-1 algorithm is based on the analytic modal solution of the diffusion equation for constant conditions, combined with polynomial corrective terms that embody the information on the deviation from constant conditions. The new algorithm is verified by comparing the results to a finite difference solution over a large number of randomly generated operation histories. Furthermore, comparison to state-of-the-art algorithms used in fuel performance codes demonstrates that the accuracy of PolyPole-1 is superior to other algorithms, with similar computational effort. Finally, the concept of PolyPole-1 may be extended to the solution of the general problem of intra-granular fission gas diffusion during non-equilibrium trapping and resolution, which will be the subject of future work.
Numerical study of combustion processes in afterburners
NASA Technical Reports Server (NTRS)
Zhou, Xiaoqing; Zhang, Xiaochun
1986-01-01
Mathematical models and numerical methods are presented for computer modeling of aeroengine afterburners. A computer code GEMCHIP is described briefly. The algorithms SIMPLER, for gas flow predictions, and DROPLET, for droplet flow calculations, are incorporated in this code. The block correction technique is adopted to facilitate convergence. The method of handling irregular shapes of combustors and flameholders is described. The predicted results for a low-bypass-ratio turbofan afterburner in the cases of gaseous combustion and multiphase spray combustion are provided and analyzed, and engineering guides for afterburner optimization are presented.
A numerical algorithm for the explicit calculation of SU(N) and SL(N,C) Clebsch-Gordan coefficients
Alex, Arne; Delft, Jan von; Kalus, Matthias; Huckleberry, Alan
2011-02-15
We present an algorithm for the explicit numerical calculation of SU(N) and SL(N,C) Clebsch-Gordan coefficients, based on the Gelfand-Tsetlin pattern calculus. Our algorithm is well suited for numerical implementation; we include a computer code in an appendix. Our exposition presumes only familiarity with the representation theory of SU(2).
NASA Technical Reports Server (NTRS)
Nacozy, P. E.
1984-01-01
The equations of motion are developed for a perfectly flexible, inelastic tether with a satellite at its extremity. The tether is attached to a space vehicle in orbit. The tether is allowed to possess electrical conductivity. A numerical solution algorithm to provide the motion of the tether and satellite system is presented. The resulting differential equations can be solved by various existing standard numerical integration computer programs. The resulting differential equations allow the introduction of approximations that can lead to analytical, approximate general solutions. The differential equations allow more dynamical insight of the motion.
Analysis of V-cycle multigrid algorithms for forms defined by numerical quadrature
Bramble, J.H. . Dept. of Mathematics); Goldstein, C.I.; Pasciak, J.E. . Applied Mathematics Dept.)
1994-05-01
The authors describe and analyze certain V-cycle multigrid algorithms with forms defined by numerical quadrature applied to the approximation of symmetric second-order elliptic boundary value problems. This approach can be used for the efficient solution of finite element systems resulting from numerical quadrature as well as systems arising from finite difference discretizations. The results are based on a regularity free theory and hence apply to meshes with local grid refinement as well as the quasi-uniform case. It is shown that uniform (independent of the number of levels) convergence rates often hold for appropriately defined V-cycle algorithms with as few as one smoothing per grid. These results hold even on applications without full elliptic regularity, e.g., a domain in R[sup 2] with a crack.
NASA Astrophysics Data System (ADS)
Kim, J.; Sonnenthal, E. L.; Rutqvist, J.
2011-12-01
Rigorous modeling of coupling between fluid, heat, and geomechanics (thermo-poro-mechanics), in fractured porous media is one of the important and difficult topics in geothermal reservoir simulation, because the physics are highly nonlinear and strongly coupled. Coupled fluid/heat flow and geomechanics are investigated using the multiple interacting continua (MINC) method as applied to naturally fractured media. In this study, we generalize constitutive relations for the isothermal elastic dual porosity model proposed by Berryman (2002) to those for the non-isothermal elastic/elastoplastic multiple porosity model, and derive the coupling coefficients of coupled fluid/heat flow and geomechanics and constraints of the coefficients. When the off-diagonal terms of the total compressibility matrix for the flow problem are zero, the upscaled drained bulk modulus for geomechanics becomes the harmonic average of drained bulk moduli of the multiple continua. In this case, the drained elastic/elastoplastic moduli for mechanics are determined by a combination of the drained moduli and volume fractions in multiple porosity materials. We also determine a relation between local strains of all multiple porosity materials in a gridblock and the global strain of the gridblock, from which we can track local and global elastic/plastic variables. For elastoplasticity, the return mapping is performed for all multiple porosity materials in the gridblock. For numerical implementation, we employ and extend the fixed-stress sequential method of the single porosity model to coupled fluid/heat flow and geomechanics in multiple porosity systems, because it provides numerical stability and high accuracy. This sequential scheme can be easily implemented by using a porosity function and its corresponding porosity correction, making use of the existing robust flow and geomechanics simulators. We implemented the proposed modeling and numerical algorithm to the reaction transport simulator
Numerical advection algorithms and their role in atmospheric transport and chemistry models
NASA Technical Reports Server (NTRS)
Rood, Richard B.
1987-01-01
During the last 35 years, well over 100 algorithms for modeling advection processes have been described and tested. This review summarizes the development and improvements that have taken place. The nature of the errors caused by numerical approximation to the advection equation are highlighted. Then the particular devices that have been proposed to remedy these errors are discussed. The extensive literature comparing transport algorithms is reviewed. Although there is no clear cut 'best' algorithm, several conclusions can be made. Spectral and pseudospectral techniques consistently provide the highest degree of accuracy, but expense and difficulties assuring positive mixing ratios are serious drawbacks. Schemes which consider fluid slabs bounded by grid points (volume schemes), rather than the simple specification of constituent values at the grid points, provide accurate positive definite results.
The double pendulum: a numerical study
NASA Astrophysics Data System (ADS)
Calvão, A. M.; Penna, T. J. P.
2015-07-01
Analysis and characterization of dynamical systems is a common task in computational physics. It frequently demands new algorithms for finding solutions and new techniques for analysing the results. Here we review some of these algorithms and techniques in the study of the double pendulum, which, despite being a very simple mechanical system, can display complex behaviour. Even though it has been studied before (Yu and Bi 1998 J. Sound Vib. 217 691736; Stachowiak and Okada 2006 Chaos, Solitons & Fractals 29 417422; Rafat, Wheatland and Bedding 2009 Am. J. Phys. 77 216-23; Levien and Tan 1993 Am. J. Phys. 61 103844), here we present a deeper discussion of the several methods and algorithms that are used in typical studies of dynamical systems. In addition, we present new results obtained through techniques commonly used in the analysis of complex systems.
Hu, Shaoxing; Xu, Shike; Wang, Duhu; Zhang, Aiwu
2015-01-01
Aiming at addressing the problem of high computational cost of the traditional Kalman filter in SINS/GPS, a practical optimization algorithm with offline-derivation and parallel processing methods based on the numerical characteristics of the system is presented in this paper. The algorithm exploits the sparseness and/or symmetry of matrices to simplify the computational procedure. Thus plenty of invalid operations can be avoided by offline derivation using a block matrix technique. For enhanced efficiency, a new parallel computational mechanism is established by subdividing and restructuring calculation processes after analyzing the extracted “useful” data. As a result, the algorithm saves about 90% of the CPU processing time and 66% of the memory usage needed in a classical Kalman filter. Meanwhile, the method as a numerical approach needs no precise-loss transformation/approximation of system modules and the accuracy suffers little in comparison with the filter before computational optimization. Furthermore, since no complicated matrix theories are needed, the algorithm can be easily transplanted into other modified filters as a secondary optimization method to achieve further efficiency. PMID:26569247
Hu, Shaoxing; Xu, Shike; Wang, Duhu; Zhang, Aiwu
2015-01-01
Aiming at addressing the problem of high computational cost of the traditional Kalman filter in SINS/GPS, a practical optimization algorithm with offline-derivation and parallel processing methods based on the numerical characteristics of the system is presented in this paper. The algorithm exploits the sparseness and/or symmetry of matrices to simplify the computational procedure. Thus plenty of invalid operations can be avoided by offline derivation using a block matrix technique. For enhanced efficiency, a new parallel computational mechanism is established by subdividing and restructuring calculation processes after analyzing the extracted "useful" data. As a result, the algorithm saves about 90% of the CPU processing time and 66% of the memory usage needed in a classical Kalman filter. Meanwhile, the method as a numerical approach needs no precise-loss transformation/approximation of system modules and the accuracy suffers little in comparison with the filter before computational optimization. Furthermore, since no complicated matrix theories are needed, the algorithm can be easily transplanted into other modified filters as a secondary optimization method to achieve further efficiency. PMID:26569247
Algorithms for verbal autopsies: a validation study in Kenyan children.
Quigley, M. A.; Armstrong Schellenberg, J. R.; Snow, R. W.
1996-01-01
The verbal autopsy (VA) questionnaire is a widely used method for collecting information on cause-specific mortality where the medical certification of deaths in childhood is incomplete. This paper discusses review by physicians and expert algorithms as approaches to ascribing cause of deaths from the VA questionnaire and proposes an alternative, data-derived approach. In this validation study, the relatives of 295 children who had died in hospital were interviewed using a VA questionnaire. The children were assigned causes of death using data-derived algorithms obtained under logistic regression and using expert algorithms. For most causes of death, the data-derived algorithms and expert algorithms yielded similar levels of diagnostic accuracy. However, a data-derived algorithm for malaria gave a sensitivity of 71% (95% Cl: 58-84%), which was significantly higher than the sensitivity of 47% obtained under an expert algorithm. The need for exploring this and other ways in which the VA technique can be improved are discussed. The implications of less-than-perfect sensitivity and specificity are explored using numerical examples. Misclassification bias should be taken into consideration when planning and evaluating epidemiological studies. PMID:8706229
Numerical studies of solar chromospheric jets
NASA Astrophysics Data System (ADS)
Iijima, Haruhisa
2016-03-01
The solar chromospheric jet is one of the most characteristic structures near the solar surface. The quantitative understanding of chromospheric jets is of substantial importance for not only the partially ionized phenomena in the chromosphere but also the energy input and dissipation processes in the corona. In this dissertation, the formation and dynamics of chromospheric jets are investigated using the radiation magnetohydrodynamic simulations. We newly develop a numerical code for the radiation magnetohydrodynamic simulations of the comprehensive modeling of solar atmosphere. Because the solar chromosphere is highly nonlinear, magnetic pressure dominated, and turbulent, a robust and high-resolution numerical scheme is required. In Chapter 2, we propose a new algorithm for the simulation of magnetohydrodynamics. Through the test problems and accuracy analyses, the proposed scheme is proved to satisfy the requirements. In Chapter 3, the effect of the non-local radiation energy transport, Spitzer-type thermal conduction, latent heat of partial ionization and molecule formation, and gravity are implemented to the magnetohydrodynamic code. The numerical schemes for the radiation transport and thermal conduction is carefully chosen in a view of the efficiency and compatibility with the parallel computation. Based on the developed radiation magnetohydrodynamic code, the formation and dynamics of chromospheric jets are investigated. In Chapter 4, we investigate the dependence of chromospheric jets on the coronal temperature in the two-dimensional simulations. Various scale of chromospheric jets with the parabolic trajectory are found with the maximum height of 2–8 Mm, lifetime of 2–7 min, maximum upward velocity of 10– 50 km/s, and deceleration of 100–350 m/s2. We find that chromospheric jets are more elongated under the cool corona and shorter under the hot corona. We also find that the pressure gradient force caused by the periodic shock waves accelerates
NASA Technical Reports Server (NTRS)
Daso, E. O.
1986-01-01
An implicit approximate factorization algorithm is employed to quantify the parametric effects of Courant number and artificial smoothing on numerical solutions of the unsteady 3-D Euler equations for a windmilling propeller (low speed) flow field. The results show that propeller global or performance chracteristics vary strongly with Courant number and artificial dissipation parameters, though the variation is such less severe at high Courant numbers. Candidate sets of Courant number and dissipation parameters could result in parameter-dependent solutions. Parameter-independent numerical solutions can be obtained if low values of the dissipation parameter-time step ratio are used in the computations. Furthermore, it is realized that too much artificial damping can degrade numerical stability. Finally, it is demonstrated that highly resolved meshes may, in some cases, delay convergence, thereby suggesting some optimum cell size for a given flow solution. It is suspected that improper boundary treatment may account for the cell size constraint.
NASA Astrophysics Data System (ADS)
Degtyarev, Alexander; Khramushin, Vasily
2016-02-01
The paper deals with the computer implementation of direct computational experiments in fluid mechanics, constructed on the basis of the approach developed by the authors. The proposed approach allows the use of explicit numerical scheme, which is an important condition for increasing the effciency of the algorithms developed by numerical procedures with natural parallelism. The paper examines the main objects and operations that let you manage computational experiments and monitor the status of the computation process. Special attention is given to a) realization of tensor representations of numerical schemes for direct simulation; b) realization of representation of large particles of a continuous medium motion in two coordinate systems (global and mobile); c) computing operations in the projections of coordinate systems, direct and inverse transformation in these systems. Particular attention is paid to the use of hardware and software of modern computer systems.
Evolutionary algorithms, simulated annealing, and Tabu search: a comparative study
NASA Astrophysics Data System (ADS)
Youssef, Habib; Sait, Sadiq M.; Adiche, Hakim
1998-10-01
Evolutionary algorithms, simulated annealing (SA), and Tabu Search (TS) are general iterative algorithms for combinatorial optimization. The term evolutionary algorithm is used to refer to any probabilistic algorithm whose design is inspired by evolutionary mechanisms found in biological species. Most widely known algorithms of this category are Genetic Algorithms (GA). GA, SA, and TS have been found to be very effective and robust in solving numerous problems from a wide range of application domains.Furthermore, they are even suitable for ill-posed problems where some of the parameters are not known before hand. These properties are lacking in all traditional optimization techniques. In this paper we perform a comparative study among GA, SA, and TS. These algorithms have many similarities, but they also possess distinctive features, mainly in their strategies for searching the solution state space. the three heuristics are applied on the same optimization problem and compared with respect to (1) quality of the best solution identified by each heuristic, (2) progress of the search from initial solution(s) until stopping criteria are met, (3) the progress of the cost of the best solution as a function of time, and (4) the number of solutions found at successive intervals of the cost function. The benchmark problem was is the floorplanning of very large scale integrated circuits. This is a hard multi-criteria optimization problem. Fuzzy logic is used to combine all objective criteria into a single fuzzy evaluation function, which is then used to rate competing solutions.
Numerical algorithms for computations of feedback laws arising in control of flexible systems
NASA Technical Reports Server (NTRS)
Lasiecka, Irena
1989-01-01
Several continuous models will be examined, which describe flexible structures with boundary or point control/observation. Issues related to the computation of feedback laws are examined (particularly stabilizing feedbacks) with sensors and actuators located either on the boundary or at specific point locations of the structure. One of the main difficulties is due to the great sensitivity of the system (hyperbolic systems with unbounded control actions), with respect to perturbations caused either by uncertainty of the model or by the errors introduced in implementing numerical algorithms. Thus, special care must be taken in the choice of the appropriate numerical schemes which eventually lead to implementable finite dimensional solutions. Finite dimensional algorithms are constructed on a basis of a priority analysis of the properties of the original, continuous (infinite diversional) systems with the following criteria in mind: (1) convergence and stability of the algorithms and (2) robustness (reasonable insensitivity with respect to the unknown parameters of the systems). Examples with mixed finite element methods and spectral methods are provided.
Chen, Deng-kai; Gu, Rong; Gu, Yu-feng; Yu, Sui-huai
2016-01-01
Consumers' Kansei needs reflect their perception about a product and always consist of a large number of adjectives. Reducing the dimension complexity of these needs to extract primary words not only enables the target product to be explicitly positioned, but also provides a convenient design basis for designers engaging in design work. Accordingly, this study employs a numerical design structure matrix (NDSM) by parameterizing a conventional DSM and integrating genetic algorithms to find optimum Kansei clusters. A four-point scale method is applied to assign link weights of every two Kansei adjectives as values of cells when constructing an NDSM. Genetic algorithms are used to cluster the Kansei NDSM and find optimum clusters. Furthermore, the process of the proposed method is presented. The details of the proposed approach are illustrated using an example of electronic scooter for Kansei needs clustering. The case study reveals that the proposed method is promising for clustering Kansei needs adjectives in product emotional design.
Godfrey, Brendan B.; Vay, Jean-Luc
2013-09-01
Rapidly growing numerical instabilities routinely occur in multidimensional particle-in-cell computer simulations of plasma-based particle accelerators, astrophysical phenomena, and relativistic charged particle beams. Reducing instability growth to acceptable levels has necessitated higher resolution grids, high-order field solvers, current filtering, etc. except for certain ratios of the time step to the axial cell size, for which numerical growth rates and saturation levels are reduced substantially. This paper derives and solves the cold beam dispersion relation for numerical instabilities in multidimensional, relativistic, electromagnetic particle-in-cell programs employing either the standard or the Cole–Karkkainnen finite difference field solver on a staggered mesh and the common Esirkepov current-gathering algorithm. Good overall agreement is achieved with previously reported results of the WARP code. In particular, the existence of select time steps for which instabilities are minimized is explained. Additionally, an alternative field interpolation algorithm is proposed for which instabilities are almost completely eliminated for a particular time step in ultra-relativistic simulations.
A Novel Quantum-Behaved Bat Algorithm with Mean Best Position Directed for Numerical Optimization
Zhu, Wenyong; Liu, Zijuan; Duan, Qingyan; Cao, Long
2016-01-01
This paper proposes a novel quantum-behaved bat algorithm with the direction of mean best position (QMBA). In QMBA, the position of each bat is mainly updated by the current optimal solution in the early stage of searching and in the late search it also depends on the mean best position which can enhance the convergence speed of the algorithm. During the process of searching, quantum behavior of bats is introduced which is beneficial to jump out of local optimal solution and make the quantum-behaved bats not easily fall into local optimal solution, and it has better ability to adapt complex environment. Meanwhile, QMBA makes good use of statistical information of best position which bats had experienced to generate better quality solutions. This approach not only inherits the characteristic of quick convergence, simplicity, and easy implementation of original bat algorithm, but also increases the diversity of population and improves the accuracy of solution. Twenty-four benchmark test functions are tested and compared with other variant bat algorithms for numerical optimization the simulation results show that this approach is simple and efficient and can achieve a more accurate solution. PMID:27293424
A Novel Quantum-Behaved Bat Algorithm with Mean Best Position Directed for Numerical Optimization.
Zhu, Binglian; Zhu, Wenyong; Liu, Zijuan; Duan, Qingyan; Cao, Long
2016-01-01
This paper proposes a novel quantum-behaved bat algorithm with the direction of mean best position (QMBA). In QMBA, the position of each bat is mainly updated by the current optimal solution in the early stage of searching and in the late search it also depends on the mean best position which can enhance the convergence speed of the algorithm. During the process of searching, quantum behavior of bats is introduced which is beneficial to jump out of local optimal solution and make the quantum-behaved bats not easily fall into local optimal solution, and it has better ability to adapt complex environment. Meanwhile, QMBA makes good use of statistical information of best position which bats had experienced to generate better quality solutions. This approach not only inherits the characteristic of quick convergence, simplicity, and easy implementation of original bat algorithm, but also increases the diversity of population and improves the accuracy of solution. Twenty-four benchmark test functions are tested and compared with other variant bat algorithms for numerical optimization the simulation results show that this approach is simple and efficient and can achieve a more accurate solution. PMID:27293424
Comparative Study of Two Automatic Registration Algorithms
NASA Astrophysics Data System (ADS)
Grant, D.; Bethel, J.; Crawford, M.
2013-10-01
The Iterative Closest Point (ICP) algorithm is prevalent for the automatic fine registration of overlapping pairs of terrestrial laser scanning (TLS) data. This method along with its vast number of variants, obtains the least squares parameters that are necessary to align the TLS data by minimizing some distance metric between the scans. The ICP algorithm uses a "model-data" concept in which the scans obtain differential treatment in the registration process depending on whether they were assigned to be the "model" or "data". For each of the "data" points, corresponding points from the "model" are sought. Another concept of "symmetric correspondence" was proposed in the Point-to-Plane (P2P) algorithm, where both scans are treated equally in the registration process. The P2P method establishes correspondences on both scans and minimizes the point-to-plane distances between the scans by simultaneously considering the stochastic properties of both scans. This paper studies both the ICP and P2P algorithms in terms of their consistency in registration parameters for pairs of TLS data. The question being investigated in this paper is, should scan A be registered to scan B, will the parameters be the same if scan B were registered to scan A? Experiments were conducted with eight pairs of real TLS data which were registered by the two algorithms in the forward (scan A to scan B) and backward (scan B to scan A) modes and the results were compared. The P2P algorithm was found to be more consistent than the ICP algorithm. The differences in registration accuracy between the forward and backward modes were negligible when using the P2P algorithm (mean difference of 0.03 mm). However, the ICP had a mean difference of 4.26 mm. Each scan was also transformed by the forward and backward parameters of the two algorithms and the misclosure computed. The mean misclosure for the P2P algorithm was 0.80 mm while that for the ICP algorithm was 5.39 mm. The conclusion from this study is
A semi-numerical algorithm for instability of compressible multilayered structures
NASA Astrophysics Data System (ADS)
Tang, Shan; Yang, Yang; Peng, Xiang He; Liu, Wing Kam; Huang, Xiao Xu; Elkhodary, Khalil
2015-07-01
A computational method is proposed for the analysis and prediction of instability (wrinkling or necking) of multilayered compressible plates and sheets made by metals or polymers under plane strain conditions. In previous works, a basic assumption (or a physical argument) that has been frequently made is that materials are incompressible to simplify mathematical derivations. To account for the compressibility of metals and polymers (the lower Poisson's ratio leads to the more compressible material), we propose a combined semi-numerical algorithm and finite element method for instability analysis. Our proposed algorithm is herein verified by comparing its predictions with published results in literature for thin films with polymer/metal substrates and for polymer/metal systems. The new combined method is then used to predict the effects of compressibility on instability behaviors. Results suggest potential utility for compressibility in the design of multilayered structures.
Numerical studies of solar chromospheric jets
NASA Astrophysics Data System (ADS)
Iijima, Haruhisa
2016-03-01
The solar chromospheric jet is one of the most characteristic structures near the solar surface. The quantitative understanding of chromospheric jets is of substantial importance for not only the partially ionized phenomena in the chromosphere but also the energy input and dissipation processes in the corona. In this dissertation, the formation and dynamics of chromospheric jets are investigated using the radiation magnetohydrodynamic simulations. We newly develop a numerical code for the radiation magnetohydrodynamic simulations of the comprehensive modeling of solar atmosphere. Because the solar chromosphere is highly nonlinear, magnetic pressure dominated, and turbulent, a robust and high-resolution numerical scheme is required. In Chapter 2, we propose a new algorithm for the simulation of magnetohydrodynamics. Through the test problems and accuracy analyses, the proposed scheme is proved to satisfy the requirements. In Chapter 3, the effect of the non-local radiation energy transport, Spitzer-type thermal conduction, latent heat of partial ionization and molecule formation, and gravity are implemented to the magnetohydrodynamic code. The numerical schemes for the radiation transport and thermal conduction is carefully chosen in a view of the efficiency and compatibility with the parallel computation. Based on the developed radiation magnetohydrodynamic code, the formation and dynamics of chromospheric jets are investigated. In Chapter 4, we investigate the dependence of chromospheric jets on the coronal temperature in the two-dimensional simulations. Various scale of chromospheric jets with the parabolic trajectory are found with the maximum height of 2-8 Mm, lifetime of 2-7 min, maximum upward velocity of 10-50 km/s, and deceleration of 100-350 m/s2. We find that chromospheric jets are more elongated under the cool corona and shorter under the hot corona. We also find that the pressure gradient force caused by the periodic shock waves accelerates some of the
NASA Astrophysics Data System (ADS)
Bor, E.; Turduev, M.; Kurt, H.
2016-08-01
Photonic structure designs based on optimization algorithms provide superior properties compared to those using intuition-based approaches. In the present study, we numerically and experimentally demonstrate subwavelength focusing of light using wavelength scale absorption-free dielectric scattering objects embedded in an air background. An optimization algorithm based on differential evolution integrated into the finite-difference time-domain method was applied to determine the locations of each circular dielectric object with a constant radius and refractive index. The multiobjective cost function defined inside the algorithm ensures strong focusing of light with low intensity side lobes. The temporal and spectral responses of the designed compact photonic structure provided a beam spot size in air with a full width at half maximum value of 0.19λ, where λ is the wavelength of light. The experiments were carried out in the microwave region to verify numerical findings, and very good agreement between the two approaches was found. The subwavelength light focusing is associated with a strong interference effect due to nonuniformly arranged scatterers and an irregular index gradient. Improving the focusing capability of optical elements by surpassing the diffraction limit of light is of paramount importance in optical imaging, lithography, data storage, and strong light-matter interaction.
NASA Astrophysics Data System (ADS)
Wilkie, George J.; Dorland, William
2016-05-01
The δf particle-in-cell algorithm has been a useful tool in studying the physics of plasmas, particularly turbulent magnetized plasmas in the context of gyrokinetics. The reduction in noise due to not having to resolve the full distribution function indicates an efficiency advantage over the standard ("full-f") particle-in-cell. Despite its successes, the algorithm behaves strangely in some circumstances. In this work, we document a fully resolved numerical instability that occurs in the simplest of multiple-species test cases: the electrostatic ΩH mode. There is also a poorly understood numerical instability that occurs when one is under-resolved in particle number, which may require a prohibitively large number of particles to stabilize. Both of these are independent of the time-stepping scheme, and we conclude that they exist if the time advancement were exact. The exact analytic form of the algorithm is presented, and several schemes for mitigating these instabilities are also presented.
Bor, E; Turduev, M; Kurt, H
2016-01-01
Photonic structure designs based on optimization algorithms provide superior properties compared to those using intuition-based approaches. In the present study, we numerically and experimentally demonstrate subwavelength focusing of light using wavelength scale absorption-free dielectric scattering objects embedded in an air background. An optimization algorithm based on differential evolution integrated into the finite-difference time-domain method was applied to determine the locations of each circular dielectric object with a constant radius and refractive index. The multiobjective cost function defined inside the algorithm ensures strong focusing of light with low intensity side lobes. The temporal and spectral responses of the designed compact photonic structure provided a beam spot size in air with a full width at half maximum value of 0.19λ, where λ is the wavelength of light. The experiments were carried out in the microwave region to verify numerical findings, and very good agreement between the two approaches was found. The subwavelength light focusing is associated with a strong interference effect due to nonuniformly arranged scatterers and an irregular index gradient. Improving the focusing capability of optical elements by surpassing the diffraction limit of light is of paramount importance in optical imaging, lithography, data storage, and strong light-matter interaction. PMID:27477060
Bor, E.; Turduev, M.; Kurt, H.
2016-01-01
Photonic structure designs based on optimization algorithms provide superior properties compared to those using intuition-based approaches. In the present study, we numerically and experimentally demonstrate subwavelength focusing of light using wavelength scale absorption-free dielectric scattering objects embedded in an air background. An optimization algorithm based on differential evolution integrated into the finite-difference time-domain method was applied to determine the locations of each circular dielectric object with a constant radius and refractive index. The multiobjective cost function defined inside the algorithm ensures strong focusing of light with low intensity side lobes. The temporal and spectral responses of the designed compact photonic structure provided a beam spot size in air with a full width at half maximum value of 0.19λ, where λ is the wavelength of light. The experiments were carried out in the microwave region to verify numerical findings, and very good agreement between the two approaches was found. The subwavelength light focusing is associated with a strong interference effect due to nonuniformly arranged scatterers and an irregular index gradient. Improving the focusing capability of optical elements by surpassing the diffraction limit of light is of paramount importance in optical imaging, lithography, data storage, and strong light-matter interaction. PMID:27477060
Bu, Sunyoung; Huang, Jingfang; Boyer, Treavor H.; Miller, Cass T.
2010-01-01
The focus of this work is on the modeling of an ion exchange process that occurs in drinking water treatment applications. The model formulation consists of a two-scale model in which a set of microscale diffusion equations representing ion exchange resin particles that vary in size and age are coupled through a boundary condition with a macroscopic ordinary differential equation (ODE), which represents the concentration of a species in a well-mixed reactor. We introduce a new age-averaged model (AAM) that averages all ion exchange particle ages for a given size particle to avoid the expensive Monte-Carlo simulation associated with previous modeling applications. We discuss two different numerical schemes to approximate both the original Monte Carlo algorithm and the new AAM for this two-scale problem. The first scheme is based on the finite element formulation in space coupled with an existing backward-difference-formula-based ODE solver in time. The second scheme uses an integral equation based Krylov deferred correction (KDC) method and a fast elliptic solver (FES) for the resulting elliptic equations. Numerical results are presented to validate the new AAM algorithm, which is also shown to be more computationally efficient than the original Monte Carlo algorithm. We also demonstrate that the higher order KDC scheme is more efficient than the traditional finite element solution approach and this advantage becomes increasingly important as the desired accuracy of the solution increases. We also discuss issues of smoothness, which affect the efficiency of the KDC-FES approach, and outline additional algorithmic changes that would further improve the efficiency of these developing methods for a wide range of applications. PMID:20577570
NASA Astrophysics Data System (ADS)
Bu, Sunyoung; Huang, Jingfang; Boyer, Treavor H.; Miller, Cass T.
2010-07-01
The focus of this work is on the modeling of an ion exchange process that occurs in drinking water treatment applications. The model formulation consists of a two-scale model in which a set of microscale diffusion equations representing ion exchange resin particles that vary in size and age are coupled through a boundary condition with a macroscopic ordinary differential equation (ODE), which represents the concentration of a species in a well-mixed reactor. We introduce a new age-averaged model (AAM) that averages all ion exchange particle ages for a given size particle to avoid the expensive Monte-Carlo simulation associated with previous modeling applications. We discuss two different numerical schemes to approximate both the original Monte-Carlo algorithm and the new AAM for this two-scale problem. The first scheme is based on the finite element formulation in space coupled with an existing backward difference formula-based ODE solver in time. The second scheme uses an integral equation based Krylov deferred correction (KDC) method and a fast elliptic solver (FES) for the resulting elliptic equations. Numerical results are presented to validate the new AAM algorithm, which is also shown to be more computationally efficient than the original Monte-Carlo algorithm. We also demonstrate that the higher order KDC scheme is more efficient than the traditional finite element solution approach and this advantage becomes increasingly important as the desired accuracy of the solution increases. We also discuss issues of smoothness, which affect the efficiency of the KDC-FES approach, and outline additional algorithmic changes that would further improve the efficiency of these developing methods for a wide range of applications.
NASA Astrophysics Data System (ADS)
Imada, Masatoshi; Kashima, Tsuyoshi
2000-09-01
A numerical algorithm for studying strongly correlated electron systems is proposed. The groundstate wavefunction is projected out after a numerical renormalization procedure in the path integral formalism. The wavefunction is expressed from the optimized linear combination of retained states in the truncated Hilbert space with a numerically chosen basis. This algorithm does not suffer from the negative sign problem and can be applied to any type of Hamiltonian in any dimension. The efficiency is tested in examples of the Hubbard model where the basis of Slater determinants is numerically optimized. We show results on fast convergence and accuracy achieved with a small number of retained states.
New Concepts in Breast Cancer Emerge from Analyzing Clinical Data Using Numerical Algorithms
Retsky, Michael
2009-01-01
A small international group has recently challenged fundamental concepts in breast cancer. As a guiding principle in therapy, it has long been assumed that breast cancer growth is continuous. However, this group suggests tumor growth commonly includes extended periods of quasi-stable dormancy. Furthermore, surgery to remove the primary tumor often awakens distant dormant micrometastases. Accordingly, over half of all relapses in breast cancer are accelerated in this manner. This paper describes how a numerical algorithm was used to come to these conclusions. Based on these findings, a dormancy preservation therapy is proposed. PMID:19440287
NASA Technical Reports Server (NTRS)
Whyte, W. A.; Heyward, A. O.; Ponchak, D. S.; Spence, R. L.; Zuzek, J. E.
1988-01-01
A detailed description of a Numerical Arc Segmentation Algorithm for a Radio Conference (NASARC) software package for communication satellite systems planning is presented. This software provides a method of generating predetermined arc segments for use in the development of an allotment planning procedure to be carried out at the 1988 World Administrative Radio Conference (WARC - 88) on the use of the GEO and the planning of space services utilizing GEO. The features of the NASARC software package are described, and detailed information is given about the function of each of the four NASARC program modules. The results of a sample world scenario are presented and discussed.
NASA Technical Reports Server (NTRS)
Shia, Run-Lie; Ha, Yuk Lung; Wen, Jun-Shan; Yung, Yuk L.
1990-01-01
Extensive testing of the advective scheme proposed by Prather (1986) has been carried out in support of the California Institute of Technology-Jet Propulsion Laboratory two-dimensional model of the middle atmosphere. The original scheme is generalized to include higher-order moments. In addition, it is shown how well the scheme works in the presence of chemistry as well as eddy diffusion. Six types of numerical experiments including simple clock motion and pure advection in two dimensions have been investigated in detail. By comparison with analytic solutions, it is shown that the new algorithm can faithfully preserve concentration profiles, has essentially no numerical diffusion, and is superior to a typical fourth-order finite difference scheme.
A Numerical Algorithm for Complex Biological Flow in Irregular Microdevice Geometries
Nonaka, A; Miller, G H; Marshall, T; Liepmann, D; Gulati, S; Trebotich, D; Colella, P
2003-12-15
We present a numerical algorithm to simulate non-Newtonian flow in complex microdevice components. The model consists of continuum viscoelastic incompressible flow in irregular microscale geometries. Our numerical approach is the projection method of Bell, Colella and Glaz (BCG) to impose the incompressibility constraint coupled with the polymeric stress splitting discretization of Trebotich, Colella and Miller (TCM). In this approach we exploit the hyperbolic structure of the equations of motion to achieve higher resolution in the presence of strong gradients and to gain an order of magnitude in the timestep. We also extend BCG and TCM to an embedded boundary method to treat irregular domain geometries which exist in microdevices. Our method allows for particle representation in a continuum fluid. We present preliminary results for incompressible viscous flow with comparison to flow of DNA and simulants in microchannels and other components used in chem/bio microdevices.
NASA Astrophysics Data System (ADS)
Acebrón, Juan A.; Rodríguez-Rozas, Ángel
2013-10-01
An efficient numerical method based on a probabilistic representation for the Vlasov-Poisson system of equations in the Fourier space has been derived. This has been done theoretically for arbitrary dimensional problems, and particularized to unidimensional problems for numerical purposes. Such a representation has been validated theoretically in the linear regime comparing the solution obtained with the classical results of the linear Landau damping. The numerical strategy followed requires generating suitable random trees combined with a Padé approximant for approximating accurately a given divergent series. Such series are obtained by summing the partial contributions to the solution coming from trees with arbitrary number of branches. These contributions, coming in general from multi-dimensional definite integrals, are efficiently computed by a quasi-Monte Carlo method. It is shown how the accuracy of the method can be effectively increased by considering more terms of the series. The new representation was used successfully to develop a Probabilistic Domain Decomposition method suited for massively parallel computers, which improves the scalability found in classical methods. Finally, a few numerical examples based on classical phenomena such as the non-linear Landau damping, and the two streaming instability are given, illustrating the remarkable performance of the algorithm, when compared the results with those obtained using a classical method.
Numerical and experimental study of rotating jet flows
NASA Astrophysics Data System (ADS)
Shin, Seungwon; Che, Zhizhao; Kahouadji, Lyes; Matar, Omar; Chergui, Jalel; Juric, Damir
2015-11-01
Rotating jets are investigated through experimental measurements and numerical simulations. The experiments are performed on a rotating jet rig and the effects of a range of parameters controlling the liquid jet are investigated, e.g. jet flow rate, rotation speed, jet diameter, etc. Different regimes of the jet morphology are identified, and the dependence on several dimensionless numbers is studied, e.g. Reynolds number, Weber number, etc. The breakup process of droplets is visualized through high speed imaging. Full three-dimensional direct numerical simulations are performed using BLUE, a massively parallel two-phase flow code. The novel interface algorithms in BLUE track the gas-liquid interface through a wide dynamic range including ligament formation, break up and rupture. EPSRC Programme Grant, MEMPHIS, EP/K0039761/1.
Fast numerical algorithms for fitting multiresolution hybrid shape models to brain MRI.
Vemuri, B C; Guo, Y; Lai, S H; Leonard, C M
1997-09-01
In this paper, we present new and fast numerical algorithms for shape recovery from brain MRI using multiresolution hybrid shape models. In this modeling framework, shapes are represented by a core rigid shape characterized by a superquadric function and a superimposed displacement function which is characterized by a membrane spline discretized using the finite-element method. Fitting the model to brain MRI data is cast as an energy minimization problem which is solved numerically. We present three new computational methods for model fitting to data. These methods involve novel mathematical derivations that lead to efficient numerical solutions of the model fitting problem. The first method involves using the nonlinear conjugate gradient technique with a diagonal Hessian preconditioner. The second method involves the nonlinear conjugate gradient in the outer loop for solving global parameters of the model and a preconditioned conjugate gradient scheme for solving the local parameters of the model. The third method involves the nonlinear conjugate gradient in the outer loop for solving the global parameters and a combination of the Schur complement formula and the alternating direction-implicit method for solving the local parameters of the model. We demonstrate the efficiency of our model fitting methods via experiments on several MR brain scans. PMID:9873915
NASA Astrophysics Data System (ADS)
Bouallegue, Kais; Chaari, Abdessattar
In this study, one propose to study a numeric type strategy permitting the generation of any shape of path in view of the scheduling of the trajectories for a car-like mobile robot where the planned motions considered are continuous sequences in the space of the robot. These paths are programmed in order to have some types of closed or open trajectories. One is interested in the motion control of the robot from an initial position to a final position while optimizing the consumed energy in its alternated circular motion on both sides of the segment joining these two points. In this study, one presents a new method based on a numeric approach conceived from the kinematics equations of the robot. This new technique of numeric, adaptive and dynamic control of the robot is implemented on DSP21065L of the SHARC family. This algorithm assures the robot control of an initial position of departure to a final position of arrival without the existence of obstacles.
Numerical and Experimental Study of Levee Breach
NASA Astrophysics Data System (ADS)
Elalfy, E. Y.; LaRocque, L.; Riahi-Nezhad, C. K.; Chaudhry, H.
2014-12-01
Levees are constructed along rivers and channels for flood protection. Failure of these levees can cause loss of life and property damage. A better understanding of the flow field from a levee breach allows the decision maker to assess risks and to prepare emergency plans. For this purpose, a two-dimensional numerical model is developed to simulate the levee breach. The model solves the shallow-water equations using the MacCormack explicit, finite- difference two-step, predictor-corrector scheme. The scheme is second-order accurate in time and space. The artificial viscosity technique is used to smooth the high-frequency oscillations in the computed results. The numerical results compare satisfactorily with the experimental results. A parametric study is carried-out to investigate the effect of main channel width, breach width on the computed flow field.
Numerical study of flow turning phenomenon
NASA Astrophysics Data System (ADS)
Baum, J. D.; Levine, J. N.
1986-01-01
A research project is currently being conducted that is to provide an understanding of the physical mechanisms by which energy is exchanged between the mean and acoustic flowfields in resonant combustion chambers, giving particular attention to solid rocket motors. The present paper is concerned with progress which has been made toward the understanding of the 'flow turning' phenomenon. This term is used to describe the loss of acoustic energy by the acoustic field in the combustor resulting from the inflow of combustion products through the lateral boundary of a combustion chamber containing longitudinal acoustic waves. Attention is given to the modeling of flow turning, acoustic refraction, the numerical solution, numerical results, acoustic wave propagation with no mean flow, and a flow turning study. The discussed research verifies the existence of the flow turning loss phenomenon.
NASA Astrophysics Data System (ADS)
Li, Yiming
2007-12-01
This symposium is an open forum for discussion on the current trends and future directions of physical modeling, mathematical theory, and numerical algorithm in electrical and electronic engineering. The goal is for computational scientists and engineers, computer scientists, applied mathematicians, physicists, and researchers to present their recent advances and exchange experience. We welcome contributions from researchers of academia and industry. All papers to be presented in this symposium have carefully been reviewed and selected. They include semiconductor devices, circuit theory, statistical signal processing, design optimization, network design, intelligent transportation system, and wireless communication. Welcome to this interdisciplinary symposium in International Conference of Computational Methods in Sciences and Engineering (ICCMSE 2007). Look forward to seeing you in Corfu, Greece!
Unquenched Studies Using the Truncated Determinant Algorithm
A. Duncan, E. Eichten and H. Thacker
2001-11-29
A truncated determinant algorithm is used to study the physical effects of the quark eigenmodes associated with eigenvalues below 420 MeV. This initial high statistics study focuses on coarse (6{sup 4}) lattices (with O(a{sup 2}) improved gauge action), light internal quark masses and large physical volumes. Three features of full QCD are examined: topological charge distributions, string breaking as observed in the static energy and the eta prime mass.
Dong, S.
2015-02-15
We present a family of physical formulations, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of N (N⩾2) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The N-phase formulations stem from a phase field model we developed in a recent work based on the conservations of mass/momentum, and the second law of thermodynamics. The introduction of general order parameters leads to an extremely strongly-coupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. We show that the increased complexity in the form of the phase field equations associated with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the (N−1) strongly-coupled phase field equations for general order parameters into 2(N−1) Helmholtz-type equations that are completely de-coupled from one another. This leads to a computational complexity comparable to that for the simplified phase field equations associated with certain special choice of the order parameters. We demonstrate the capabilities of the method developed herein using several test problems involving multiple fluid phases and large contrasts in densities and viscosities among the multitude of fluids. In particular, by comparing simulation results with the Langmuir–de Gennes theory of floating liquid lenses we show that the method using general order parameters produces physically accurate results for multiple fluid phases.
A Numerical Study of Feathering Instability
NASA Astrophysics Data System (ADS)
Lee, Wing-Kit; Wang, Hsiang-Hsu
2016-06-01
The stability of a spiral shock of self-gravitating, magnetized interstellar medium is studied by performing two-dimensional numerical simulations of a local patch of tight-winding spiral arm. As previously suggested by the linear studies, two types of instabilities are identified, namely, wiggle instability and feathering instability. The former instability occurs in the hydrodynamics limit and results in short wavelength perturbations. On the other hand, the feathering instability requires both self-gravitating and magnetic fields and results in wider structures.
Numerical Arc Segmentation Algorithm for a Radio Conference (NASARC), version 4.0: User's manual
NASA Technical Reports Server (NTRS)
Whyte, Wayne A., Jr.; Heyward, Ann O.; Ponchak, Denise S.; Spence, Rodney L.; Zuzek, John E.
1988-01-01
The information in the NASARC (Version 4.0) Technical Manual (NASA-TM-101453) and NASARC (Version 4.0) User's Manual (NASA-TM-101454) relates to the state of Numerical Arc Segmentation Algorithm for a Radio Conference (NASARC) software development through November 1, 1988. The Technical Manual describes the NASARC concept and the algorithms used to implement the concept. The User's Manual provides information on computer system considerations, installation instructions, description of input files, and program operation instructions. Significant revisions were incorporated in the Version 4.0 software over prior versions. These revisions have further enhanced the modeling capabilities of the NASARC procedure and provide improved arrangements of predetermined arcs within the geostationary orbit. Array dimensions within the software were structured to fit within the currently available 12-megabyte memory capacity of the International Frequency Registration Board (IFRB) computer facility. A piecewise approach to predetermined arc generation in NASARC (Version 4.) allows worldwide planning problem scenarios to be accommodated within computer run time and memory constraints with enhanced likelihood and ease of solution.
Numerical Arc Segmentation Algorithm for a Radio Conference-NASARC (version 4.0) technical manual
NASA Technical Reports Server (NTRS)
Whyte, Wayne A., Jr.; Heyward, Ann O.; Ponchak, Denise S.; Spence, Rodney L.; Zuzek, John E.
1988-01-01
The information contained in the NASARC (Version 4.0) Technical Manual and NASARC (Version 4.0) User's Manual relates to the Numerical Arc Segmentation Algorithm for a Radio Conference (NASARC) software development through November 1, 1988. The Technical Manual describes the NASARC concept and the algorithms used to implement the concept. The User's Manual provides information on computer system considerations, installation instructions, description of input files, and program operation instructions. Significant revisions were incorporated in the Version 4.0 software over prior versions. These revisions have further enhanced the modeling capabilities of the NASARC procedure and provide improved arrangements of predetermined arcs within the geostationary orbits. Array dimensions within the software were structured to fit within the currently available 12 megabyte memory capacity of the International Frequency Registration Board (IFRB) computer facility. A piecewise approach to predetermined arc generation in NASARC (Version 4.0) allows worldwide planning problem scenarios to be accommodated within computer run time and memory constraints with enhanced likelihood and ease of solution.
Numerical Arc Segmentation Algorithm for a Radio Conference-NASARC, Version 2.0: User's Manual
NASA Technical Reports Server (NTRS)
Whyte, Wayne A., Jr.; Heyward, Ann O.; Ponchak, Denise S.; Spence, Rodney L.; Zuzek, John E.
1987-01-01
The information contained in the NASARC (Version 2.0) Technical Manual (NASA TM-100160) and the NASARC (Version 2.0) User's Manual (NASA TM-100161) relates to the state of the Numerical Arc Segmentation Algorithm for a Radio Conference (NASARC) software development through October 16, 1987. The technical manual describes the NASARC concept and the algorithms which are used to implement it. The User's Manual provides information on computer system considerations, installation instructions, description of input files, and program operation instructions. Significant revisions have been incorporated in the Version 2.0 software over prior versions. These revisions have enhanced the modeling capabilities of the NASARC procedure while greatly reducing the computer run time and memory requirements. Array dimensions within the software have been structured to fit into the currently available 6-megabyte memory capacity of the International Frequency Registration Board (IFRB) computer facility. A piecewise approach to predetermined arc generation in NASARC (Version 2.0) allows worldwide scenarios to be accommodated within these memory constraints while at the same time reducing computer run time.
Numerical arc segmentation algorithm for a radio conference-NASARC (version 2.0) technical manual
NASA Technical Reports Server (NTRS)
Whyte, Wayne A., Jr.; Heyward, Ann O.; Ponchak, Denise S.; Spence, Rodney L.; Zuzek, John E.
1987-01-01
The information contained in the NASARC (Version 2.0) Technical Manual (NASA TM-100160) and NASARC (Version 2.0) User's Manual (NASA TM-100161) relates to the state of NASARC software development through October 16, 1987. The Technical Manual describes the Numerical Arc Segmentation Algorithm for a Radio Conference (NASARC) concept and the algorithms used to implement the concept. The User's Manual provides information on computer system considerations, installation instructions, description of input files, and program operating instructions. Significant revisions have been incorporated in the Version 2.0 software. These revisions have enhanced the modeling capabilities of the NASARC procedure while greatly reducing the computer run time and memory requirements. Array dimensions within the software have been structured to fit within the currently available 6-megabyte memory capacity of the International Frequency Registration Board (IFRB) computer facility. A piecewise approach to predetermined arc generation in NASARC (Version 2.0) allows worldwide scenarios to be accommodated within these memory constraints while at the same time effecting an overall reduction in computer run time.
Experimental and Numerical Studies of Oceanic Overflow
NASA Astrophysics Data System (ADS)
Gibson, Thomas; Hohman, Fred; Morrison, Theresa; Reckinger, Shanon; Reckinger, Scott
2014-11-01
Oceanic overflows occur when dense water flows down a continental slope into less dense ambient water. The resulting density driven plumes occur naturally in various regions of the global ocean and affect the large-scale circulation. General circulation models currently rely on parameterizations for representing dense overflows due to resolution restrictions. The work presented here involves a direct qualitative and quantitative comparison between physical laboratory experiments and lab-scale numerical simulations. Laboratory experiments are conducted using a rotating square tank customized for idealized overflow and a high-resolution camera mounted on the table in the rotating reference frame for data collection. Corresponding numerical simulations are performed using the MIT general circulation model (MITgcm) run in the non-hydrostatic configuration. Resolution and numerical parameter studies are presented to ensure accuracy of the simulation. Laboratory and computational experiments are compared across a wide range of physical parameters, including Coriolis parameter, inflow density anomaly, and dense inflow volumetric flow rate. The results are analyzed using various calculated metrics, such as the plume velocity. Funding for this project is provided by the National Science Foundation.
Belief network algorithms: A study of performance
Jitnah, N.
1996-12-31
This abstract gives an overview of the work. We present a survey of Belief Network algorithms and propose a domain characterization system to be used as a basis for algorithm comparison and for predicting algorithm performance.
NASA Astrophysics Data System (ADS)
Kitaura, F. S.; Enßlin, T. A.
2008-09-01
We address the inverse problem of cosmic large-scale structure reconstruction from a Bayesian perspective. For a linear data model, a number of known and novel reconstruction schemes, which differ in terms of the underlying signal prior, data likelihood and numerical inverse extraregularization schemes are derived and classified. The Bayesian methodology presented in this paper tries to unify and extend the following methods: Wiener filtering, Tikhonov regularization, ridge regression, maximum entropy and inverse regularization techniques. The inverse techniques considered here are the asymptotic regularization, the Jacobi, Steepest Descent, Newton-Raphson, Landweber-Fridman and both linear and non-linear Krylov methods based on Fletcher-Reeves, Polak-Ribière and Hestenes-Stiefel conjugate gradients. The structures of the up-to-date highest performing algorithms are presented, based on an operator scheme, which permits one to exploit the power of fast Fourier transforms. Using such an implementation of the generalized Wiener filter in the novel ARGO software package, the different numerical schemes are benchmarked with one-, two- and three-dimensional problems including structured white and Poissonian noise, data windowing and blurring effects. A novel numerical Krylov scheme is shown to be superior in terms of performance and fidelity. These fast inverse methods ultimately will enable the application of sampling techniques to explore complex joint posterior distributions. We outline how the space of the dark matter density field, the peculiar velocity field and the power spectrum can jointly be investigated by a Gibbs-sampling process. Such a method can be applied for the redshift distortions correction of the observed galaxies and for time-reversal reconstructions of the initial density field.
Numerical Study of Tip Vortex Flows
NASA Technical Reports Server (NTRS)
Dacles-Mariani, Jennifer; Hafez, Mohamed
1998-01-01
This paper presents an overview and summary of the many different research work related to tip vortex flows and wake/trailing vortices as applied to practical engineering problems. As a literature survey paper, it outlines relevant analytical, theoretical, experimental and computational study found in literature. It also discusses in brief some of the fundamental aspects of the physics and its complexities. An appendix is also included. The topics included in this paper are: 1) Analytical Vortices; 2) Experimental Studies; 3) Computational Studies; 4) Wake Vortex Control and Management; 5) Wake Modeling; 6) High-Lift Systems; 7) Issues in Numerical Studies; 8) Instabilities; 9) Related Topics; 10) Visualization Tools for Vertical Flows; 11) Further Work Needed; 12) Acknowledgements; 13) References; and 14) Appendix.
Numerical Studies of Doped Iron Pnictides
NASA Astrophysics Data System (ADS)
Bishop, Christopher; Liang, Shuhua; Moreo, Adriana; Dagotto, Elbio
The phase diagram of electron-doped pnictides is studied varying the temperature, electronic density, and isotropic disorder strength and dilution via numerical studies of a three-orbital spin-fermion model with lattice degrees of freedom. Doping introduces disorder but in theoretical studies the effect of the randomly located dopants is difficult to address. Numerically the effects of electronic doping, regulated by a chemical potential, and impurity disorder at randomly selected sites can be independently controlled. It was found that the reduction with doping of the Neel and the structural transition temperatures, and the stabilization of a nematic state, is mainly controlled by the magnetic dilution due to the disorder. Fermi surface changes due to doping affect only slightly both critical temperatures. Our findings are compatible with neutron scattering and STM results, unveiling a patchy network of locally magnetically ordered anisotropic clusters, despite the isotropic disorder. The fragile tendency to nematicity intrinsic of translational invariant electronic systems needs to be supplemented by disorder and dilution to stabilize the robust nematic phase experimentally found in electron-doped 122 pnictides. National Science Foundation Grant No. DMR-1404375.
A Numerical Climate Observing Network Design Study
NASA Technical Reports Server (NTRS)
Stammer, Detlef
2003-01-01
This project was concerned with three related questions of an optimal design of a climate observing system: 1. The spatial sampling characteristics required from an ARGO system. 2. The degree to which surface observations from ARGO can be used to calibrate and test satellite remote sensing observations of sea surface salinity (SSS) as it is anticipated now. 3. The more general design of an climate observing system as it is required in the near future for CLIVAR in the Atlantic. An important question in implementing an observing system is that of the sampling density required to observe climate-related variations in the ocean. For that purpose this project was concerned with the sampling requirements for the ARGO float system, but investigated also other elements of a climate observing system. As part of this project we studied the horizontal and vertical sampling characteristics of a global ARGO system which is required to make it fully complementary to altimeter data with the goal to capture climate related variations on large spatial scales (less thanAttachment: 1000 km). We addressed this question in the framework of a numerical model study in the North Atlantic with an 1/6 horizontal resolution. The advantage of a numerical design study is the knowledge of the full model state. Sampled by a synthetic float array, model results will therefore allow to test and improve existing deployment strategies with the goal to make the system as optimal and cost-efficient as possible. Attachment: "Optimal observations for variational data assimilation".
NASA Astrophysics Data System (ADS)
Ersoy, Ozlem; Dag, Idris
2015-12-01
The solutions of the reaction-diffusion system are given by method of collocation based on the exponential B-splines. Thus the reaction-diffusion systemturns into an iterative banded algebraic matrix equation. Solution of the matrix equation is carried out byway of Thomas algorithm. The present methods test on both linear and nonlinear problems. The results are documented to compare with some earlier studies by use of L∞ and relative error norm for problems respectively.
NASA Astrophysics Data System (ADS)
García, Hermes A.; Guerrero-Bolaño, Francisco J.; Obregón-Neira, Nelson
2010-05-01
Due to both mathematical tractability and efficiency on computational resources, it is very common to find in the realm of numerical modeling in hydro-engineering that regular linearization techniques have been applied to nonlinear partial differential equations properly obtained in environmental flow studies. Sometimes this simplification is also made along with omission of nonlinear terms involved in such equations which in turn diminishes the performance of any implemented approach. This is the case for example, for contaminant transport modeling in streams. Nowadays, a traditional and one of the most common used water quality model such as QUAL2k, preserves its original algorithm, which omits nonlinear terms through linearization techniques, in spite of the continuous algorithmic development and computer power enhancement. For that reason, the main objective of this research was to generate a flexible tool for non-linear water quality modeling. The solution implemented here was based on two genetic algorithms, used in a nested way in order to find two different types of solutions sets: the first set is composed by the concentrations of the physical-chemical variables used in the modeling approach (16 variables), which satisfies the non-linear equation system. The second set, is the typical solution of the inverse problem, the parameters and constants values for the model when it is applied to a particular stream. From a total of sixteen (16) variables, thirteen (13) was modeled by using non-linear coupled equation systems and three (3) was modeled in an independent way. The model used here had a requirement of fifty (50) parameters. The nested genetic algorithm used for the numerical solution of a non-linear equation system proved to serve as a flexible tool to handle with the intrinsic non-linearity that emerges from the interactions occurring between multiple variables involved in water quality studies. However because there is a strong data limitation in
Numerical study of a magnesium hydride tank
NASA Astrophysics Data System (ADS)
Delhomme, Baptiste; de Rango, Patricia; Marty, Philippe
2012-11-01
Hydrogen storage in metal hydride tanks (MHT) is a very promising solution. Several experimental tanks, studied by different teams, have already proved the feasibility and the interesting performances of this solution. However, in much cases, an optimization of tank geometry is still needed in order to perform fast hydrogen loading. The development of efficient numerical tools is a key issue for MHT design and optimization. We propose a simple model representing a metal hydride tank exchanging its heat of reaction with a thermal fluid flow. In this model, the radial and axial discretisations have been decoupled by using Matlab® one-dimensional tools. Calculations are compared to experimental results obtained in a previous study. A good agreement is found for the loading case. The discharging case shows some discrepancies, which are discussed in this paper.
Numerical aerodynamic simulation facility feasibility study
NASA Technical Reports Server (NTRS)
1979-01-01
There were three major issues examined in the feasibility study. First, the ability of the proposed system architecture to support the anticipated workload was evaluated. Second, the throughput of the computational engine (the flow model processor) was studied using real application programs. Third, the availability reliability, and maintainability of the system were modeled. The evaluations were based on the baseline systems. The results show that the implementation of the Numerical Aerodynamic Simulation Facility, in the form considered, would indeed be a feasible project with an acceptable level of risk. The technology required (both hardware and software) either already exists or, in the case of a few parts, is expected to be announced this year. Facets of the work described include the hardware configuration, software, user language, and fault tolerance.
NASA Technical Reports Server (NTRS)
Whyte, W. A.; Heyward, A. O.; Ponchak, D. S.; Spence, R. L.; Zuzek, J. E.
1988-01-01
The Numerical Arc Segmentation Algorithm for a Radio Conference (NASARC) provides a method of generating predetermined arc segments for use in the development of an allotment planning procedure to be carried out at the 1988 World Administrative Radio Conference (WARC) on the Use of the Geostationary Satellite Orbit and the Planning of Space Services Utilizing It. Through careful selection of the predetermined arc (PDA) for each administration, flexibility can be increased in terms of choice of system technical characteristics and specific orbit location while reducing the need for coordination among administrations. The NASARC software determines pairwise compatibility between all possible service areas at discrete arc locations. NASARC then exhaustively enumerates groups of administrations whose satellites can be closely located in orbit, and finds the arc segment over which each such compatible group exists. From the set of all possible compatible groupings, groups and their associated arc segments are selected using a heuristic procedure such that a PDA is identified for each administration. Various aspects of the NASARC concept and how the software accomplishes specific features of allotment planning are discussed.
ERIC Educational Resources Information Center
Gonzalez-Vega, Laureano
1999-01-01
Using a Computer Algebra System (CAS) to help with the teaching of an elementary course in linear algebra can be one way to introduce computer algebra, numerical analysis, data structures, and algorithms. Highlights the advantages and disadvantages of this approach to the teaching of linear algebra. (Author/MM)
A numerical study of thin flame representations
Rotman, D.A.; Pindera, M.Z.
1989-08-11
In studies of reacting flows, the flame may be viewed as a moving discontinuity endowed with certain properties; notably, it acts as a source of velocity and vorticity. Asymptotic analysis shows this to be justified provided that the flame curvature is small compared to the flame thickness. Such an approach is useful when one is interested in the hydrodynamic effects of the flame on the surrounding flowfield. In numerical models of this kind it is customary to treat the discontinuity as a collection of discrete velocity blobs. In this study, we show that the velocities associated with such a representation can be very non-smooth, particularly very near to the flame surface. As an alternative, we propose the use of a finite line source as the basic flame element. Comparisons of the two flame representations are made for several simple test cases as well as for a flame propagating through an enclosure forming the tulip shape. The results show that the use of line sources eliminates spurious fluctuations in nearfield velocities thus allowing for a more accurate calculation of flame propagation and flame-flowfield interactions. 7 refs., 15 figs.
NASA Astrophysics Data System (ADS)
Korneev, Boris; Levchenko, Vadim
2016-02-01
Interaction between a shock wave and an inhomogeneity in fluid has complicated behavior, including vortex and turbulence generating, mixing, shock wave scattering and reflection. In the present paper we deal with the numerical simulation of the considered process. The Euler equations of unsteady inviscid compressible three-dimensional flow are used into the four-equation model of multicomponent flow. These equations are discretized using the RKDG numerical method. It is implemented with the help of the DiamondTorre algorithm, so the effective GPGPU solver is obtained having outstanding computing properties. With its use we carry out several sets of numerical experiments of shock-bubble interaction problem. The bubble deformation and mixture formation is observed.
Empirical study of parallel LRU simulation algorithms
NASA Technical Reports Server (NTRS)
Carr, Eric; Nicol, David M.
1994-01-01
This paper reports on the performance of five parallel algorithms for simulating a fully associative cache operating under the LRU (Least-Recently-Used) replacement policy. Three of the algorithms are SIMD, and are implemented on the MasPar MP-2 architecture. Two other algorithms are parallelizations of an efficient serial algorithm on the Intel Paragon. One SIMD algorithm is quite simple, but its cost is linear in the cache size. The two other SIMD algorithm are more complex, but have costs that are independent on the cache size. Both the second and third SIMD algorithms compute all stack distances; the second SIMD algorithm is completely general, whereas the third SIMD algorithm presumes and takes advantage of bounds on the range of reference tags. Both MIMD algorithm implemented on the Paragon are general and compute all stack distances; they differ in one step that may affect their respective scalability. We assess the strengths and weaknesses of these algorithms as a function of problem size and characteristics, and compare their performance on traces derived from execution of three SPEC benchmark programs.
NASA Technical Reports Server (NTRS)
Spratlin, Kenneth Milton
1987-01-01
An adaptive numeric predictor-corrector guidance is developed for atmospheric entry vehicles which utilize lift to achieve maximum footprint capability. Applicability of the guidance design to vehicles with a wide range of performance capabilities is desired so as to reduce the need for algorithm redesign with each new vehicle. Adaptability is desired to minimize mission-specific analysis and planning. The guidance algorithm motivation and design are presented. Performance is assessed for application of the algorithm to the NASA Entry Research Vehicle (ERV). The dispersions the guidance must be designed to handle are presented. The achievable operational footprint for expected worst-case dispersions is presented. The algorithm performs excellently for the expected dispersions and captures most of the achievable footprint.
The numerical study of first order wetting transition with two defect lines
NASA Astrophysics Data System (ADS)
Wu, X. T.
2016-09-01
The first order wetting transition with two defect lines, one near a wall and another at a distance N1, in the d = 2 Ising model is studied by the bond propagation algorithm. The numerical calculations are carried out on very large lattices with size up to 1602 × 160. The finite size effects of the first order transition in that model are discussed. The magnetization profile is also calculated. The numerical results agree with the exact results very well.
Comparative study of heart sound localization algorithms
NASA Astrophysics Data System (ADS)
Moukadem, A.; Dieterlen, A.; Hueber, N.; Brandt, C.; Raymond, P.
2011-05-01
The purpose of this document is to present a comparative study of five algorithms of heart sound localization, one of which, is a method based on radial basis function networks applied in a novel approach. The advantages and disadvantages of each method are evaluated according to a data base of 50 subjects in which there are 25 healthy subjects selected from the University Hospital of Strasbourg (HUS) and from theMARS500 project (Moscow) and 25 subjects with cardiac pathologies selected from the HUS. This study is made under the control of an experienced cardiologist. The performance of each method is evaluated by calculating the area under a receiver operating curve (AUC) and the robustness is shown against different levels of additive white Gaussian noise.
Supersymmetric Q-balls: A numerical study
Campanelli, L.; Ruggieri, M.
2008-02-15
We study numerically a class of nontopological solitons, the Q-balls, arising in a supersymmetric extension of the standard model with low-energy, gauge-mediated symmetry breaking. Taking into account the exact form of the supersymmetric potential giving rise to Q-balls, we find that there is a lower limit on the value of the charge Q in order to make them classically stable: Q > or approx. 5x10{sup 2}Q{sub cr}, where Q{sub cr} is constant depending on the parameters defining the potential and can be in the range 1 < or approx. Q{sub cr} < or approx. 10{sup 8} {sup divide} {sup 16}. If Q is the baryon number, stability with respect to the decay into protons requires Q > or approx. 10{sup 17}Q{sub cr}, while if the gravitino mass is greater then m{sub 3/2} > or approx. 61 MeV, no stable gauge-mediation supersymmetric Q-balls exist. Finally, we find that energy and radius of Q-balls can be parametrized as E{approx}{xi}{sub E}Q{sup 3/4} and R{approx}{xi}{sub R}Q{sup 1/4}, where {xi}{sub E} and {xi}{sub R} are slowly varying functions of the charge.
Numerical Study of a Convective Turbulence Encounter
NASA Technical Reports Server (NTRS)
Proctor, Fred H.; Hamilton, David W.; Bowles, Roland L.
2002-01-01
A numerical simulation of a convective turbulence event is investigated and compared with observational data. The specific case was encountered during one of NASA's flight tests and was characterized by severe turbulence. The event was associated with overshooting convective turrets that contained low to moderate radar reflectivity. Model comparisons with observations are quite favorable. Turbulence hazard metrics are proposed and applied to the numerical data set. Issues such as adequate grid size are examined.
Angus, Simon D.; Piotrowska, Monika Joanna
2014-01-01
Multi-dose radiotherapy protocols (fraction dose and timing) currently used in the clinic are the product of human selection based on habit, received wisdom, physician experience and intra-day patient timetabling. However, due to combinatorial considerations, the potential treatment protocol space for a given total dose or treatment length is enormous, even for relatively coarse search; well beyond the capacity of traditional in-vitro methods. In constrast, high fidelity numerical simulation of tumor development is well suited to the challenge. Building on our previous single-dose numerical simulation model of EMT6/Ro spheroids, a multi-dose irradiation response module is added and calibrated to the effective dose arising from 18 independent multi-dose treatment programs available in the experimental literature. With the developed model a constrained, non-linear, search for better performing cadidate protocols is conducted within the vicinity of two benchmarks by genetic algorithm (GA) techniques. After evaluating less than 0.01% of the potential benchmark protocol space, candidate protocols were identified by the GA which conferred an average of 9.4% (max benefit 16.5%) and 7.1% (13.3%) improvement (reduction) on tumour cell count compared to the two benchmarks, respectively. Noticing that a convergent phenomenon of the top performing protocols was their temporal synchronicity, a further series of numerical experiments was conducted with periodic time-gap protocols (10 h to 23 h), leading to the discovery that the performance of the GA search candidates could be replicated by 17–18 h periodic candidates. Further dynamic irradiation-response cell-phase analysis revealed that such periodicity cohered with latent EMT6/Ro cell-phase temporal patterning. Taken together, this study provides powerful evidence towards the hypothesis that even simple inter-fraction timing variations for a given fractional dose program may present a facile, and highly cost
Numerical study of a microscopic artificial swimmer
NASA Astrophysics Data System (ADS)
Gauger, Erik; Stark, Holger
2006-08-01
We present a detailed numerical study of a microscopic artificial swimmer realized recently by Dreyfus in experiments [Dreyfus , Nature 437, 862 (2005)]. It consists of an elastic filament composed of superparamagnetic particles that are linked together by DNA strands. Attached to a load particle, the resulting swimmer is actuated by an oscillating external magnetic field so that it performs a nonreciprocal motion in order to move forward. We model the superparamagnetic filament by a bead-spring configuration that resists bending like a rigid rod and whose beads experience friction with the surrounding fluid and hydrodynamic interactions with each other. We show that, aside from finite-size effects, its dynamics is governed by the dimensionless sperm number, the magnitude of the magnetic field, and the angular amplitude of the field’s oscillating direction. Then we study the mean velocity and the efficiency of the swimmer as a function of these parameters and the size of the load particle. In particular, we clarify that the real velocity of the swimmer is influenced by two main factors, namely the shape of the beating filament (determined by the sperm number and the magnetic-field strength) and the oscillation frequency. Furthermore, the load size influences the performance of the swimmer and has to be chosen as a compromise between the largest swimming velocity and the best efficiency. Finally, we demonstrate that the direction of the swimming velocity changes in a symmetry-breaking transition when the angular amplitude of the field’s oscillating direction is increased, in agreement with experiments.
Numerical Study on Mixed-mode Fracture in Reinforced Concrete
NASA Astrophysics Data System (ADS)
Yu, Rena C.; Saucedo, Luis; Ruiz, Gonzalo
2010-05-01
The object of this work is to model the propagation of fracture in mixed-mode in lightly reinforced concrete beams. When a notched beam does not have enough shear reinforcement, fracture can initiate and propagate unstably and lead to failure through diagonal tension. In order to study this phenomenon numerically, a model capable of dealing with both static and dynamic crack propagation as well as the natural transition of those two regimes is necessary. We adopt a cohesive model for concrete fracture and an interface model for the deterioration between concrete and steel re-bar, both combined with an insertion algorithm. The static process is solved by dynamic relaxation (DR) method together with a modified technique [1] to enhance convergence rate. The same DR method is used to detect a dynamic process and switch to a dynamic calculation. The numerically obtained load-displacement curves, load-CMOD curves and crack patterns fit reasonably well with their experimental counterparts, having in mind that we fed the calculations only with parameters measured experimentally.
Numerical study of Taylor bubbles with adaptive unstructured meshes
NASA Astrophysics Data System (ADS)
Xie, Zhihua; Pavlidis, Dimitrios; Percival, James; Pain, Chris; Matar, Omar; Hasan, Abbas; Azzopardi, Barry
2014-11-01
The Taylor bubble is a single long bubble which nearly fills the entire cross section of a liquid-filled circular tube. This type of bubble flow regime often occurs in gas-liquid slug flows in many industrial applications, including oil-and-gas production, chemical and nuclear reactors, and heat exchangers. The objective of this study is to investigate the fluid dynamics of Taylor bubbles rising in a vertical pipe filled with oils of extremely high viscosity (mimicking the ``heavy oils'' found in the oil-and-gas industry). A modelling and simulation framework is presented here which can modify and adapt anisotropic unstructured meshes to better represent the underlying physics of bubble rise and reduce the computational effort without sacrificing accuracy. The numerical framework consists of a mixed control-volume and finite-element formulation, a ``volume of fluid''-type method for the interface capturing based on a compressive control volume advection method, and a force-balanced algorithm for the surface tension implementation. Numerical examples of some benchmark tests and the dynamics of Taylor bubbles are presented to show the capability of this method. EPSRC Programme Grant, MEMPHIS, EP/K0039761/1.
Mathematical modeling of electrocardiograms: a numerical study.
Boulakia, Muriel; Cazeau, Serge; Fernández, Miguel A; Gerbeau, Jean-Frédéric; Zemzemi, Nejib
2010-03-01
This paper deals with the numerical simulation of electrocardiograms (ECG). Our aim is to devise a mathematical model, based on partial differential equations, which is able to provide realistic 12-lead ECGs. The main ingredients of this model are classical: the bidomain equations coupled to a phenomenological ionic model in the heart, and a generalized Laplace equation in the torso. The obtention of realistic ECGs relies on other important features--including heart-torso transmission conditions, anisotropy, cell heterogeneity and His bundle modeling--that are discussed in detail. The numerical implementation is based on state-of-the-art numerical methods: domain decomposition techniques and second order semi-implicit time marching schemes, offering a good compromise between accuracy, stability and efficiency. The numerical ECGs obtained with this approach show correct amplitudes, shapes and polarities, in all the 12 standard leads. The relevance of every modeling choice is carefully discussed and the numerical ECG sensitivity to the model parameters investigated. PMID:20033779
An efficient algorithm for numerical computations of continuous densities of states
NASA Astrophysics Data System (ADS)
Langfeld, K.; Lucini, B.; Pellegrini, R.; Rago, A.
2016-06-01
In Wang-Landau type algorithms, Monte-Carlo updates are performed with respect to the density of states, which is iteratively refined during simulations. The partition function and thermodynamic observables are then obtained by standard integration. In this work, our recently introduced method in this class (the LLR approach) is analysed and further developed. Our approach is a histogram free method particularly suited for systems with continuous degrees of freedom giving rise to a continuum density of states, as it is commonly found in lattice gauge theories and in some statistical mechanics systems. We show that the method possesses an exponential error suppression that allows us to estimate the density of states over several orders of magnitude with nearly constant relative precision. We explain how ergodicity issues can be avoided and how expectation values of arbitrary observables can be obtained within this framework. We then demonstrate the method using compact U(1) lattice gauge theory as a show case. A thorough study of the algorithm parameter dependence of the results is performed and compared with the analytically expected behaviour. We obtain high precision values for the critical coupling for the phase transition and for the peak value of the specific heat for lattice sizes ranging from 8^4 to 20^4. Our results perfectly agree with the reference values reported in the literature, which covers lattice sizes up to 18^4. Robust results for the 20^4 volume are obtained for the first time. This latter investigation, which, due to strong metastabilities developed at the pseudo-critical coupling of the system, so far has been out of reach even on supercomputers with importance sampling approaches, has been performed to high accuracy with modest computational resources. This shows the potential of the method for studies of first order phase transitions. Other situations where the method is expected to be superior to importance sampling techniques are pointed
Numerical studies of supersonic/hypersonic combustion
Yoon, W.S.; Chung, T.J. )
1992-01-01
This paper is concerned with the development of direct numerical simulations of turbulence interacting with shock waves and chemical reactions using unstructured adaptive finite element h-p methods. Reliable methods for resolving the complicated time and length scales involved in turbulence interacting with shock waves and chemical reactions are not yet available. Direct numerical simulations are here developed via Taylor-Galerkin finite element implicit scheme, with mesh refinements and spectral orders optimized such that errors are reduced where gradients of variables are large. 31 refs.
Verification Studies for Multi-Fluid Plasma Algorithms with Applications to Fast MHD Physics
NASA Astrophysics Data System (ADS)
Becker, Joe; Hakim, Ammar; Loverich, John; Stoltz, Peter
2011-10-01
In this paper we present a series of verification studies for finite volume algorithms in Nautilus, a numerical solver for fluid plasmas. Results include a set of typical Euler, Maxwell, MHD and Two-fluid benchmarks. In addition results and algorithms for a set of hyperbolic gauge cleaning schemes that can be applied to the MHD and Two-fluid systems using finite volume type methods will be presented. Finally we move onto applications in field reversed configuration (FRC) plasmas.
Numerical study of mixing in hypervelocity flows
Krishnamurthy, R.; Woods, D.M.; Chandra, S.
1996-12-31
Results are reported from an analysis of data obtained with a combustor model tested in the new pulse facility at California Institute of Technology. Data analysis was performed using the numerical code GASP. Comparisons were drawn between the predictions of the algebraic Baldwin-Lomax and the two equation k-{epsilon} turbulence models. It is concluded that the use of the simpler Baldwin-Lomax model is sufficient as it saves computational resources and yields adequate predictions, for the conditions considered here.
A numerical study of confined turbulent jets
NASA Technical Reports Server (NTRS)
Zhu, J.; Shih, T.-H.
1993-01-01
A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.
BLUM,T.
1999-09-14
The RIKEN BNL Research Center hosted its 19th workshop April 27th through May 1, 1999. The topic was Numerical Algorithms at Non-Zero Chemical Potential. QCD at a non-zero chemical potential (non-zero density) poses a long-standing unsolved challenge for lattice gauge theory. Indeed, it is the primary unresolved issue in the fundamental formulation of lattice gauge theory. The chemical potential renders conventional lattice actions complex, practically excluding the usual Monte Carlo techniques which rely on a positive definite measure for the partition function. This ''sign'' problem appears in a wide range of physical systems, ranging from strongly coupled electronic systems to QCD. The lack of a viable numerical technique at non-zero density is particularly acute since new exotic ''color superconducting'' phases of quark matter have recently been predicted in model calculations. A first principles confirmation of the phase diagram is desirable since experimental verification is not expected soon. At the workshop several proposals for new algorithms were made: cluster algorithms, direct simulation of Grassman variables, and a bosonization of the fermion determinant. All generated considerable discussion and seem worthy of continued investigation. Several interesting results using conventional algorithms were also presented: condensates in four fermion models, SU(2) gauge theory in fundamental and adjoint representations, and lessons learned from strong; coupling, non-zero temperature and heavy quarks applied to non-zero density simulations.
Zhu, Xinjun; Chen, Zhanqing; Tang, Chen; Mi, Qinghua; Yan, Xiusheng
2013-03-20
In this paper, we are concerned with denoising in experimentally obtained electronic speckle pattern interferometry (ESPI) speckle fringe patterns with poor quality. We extend the application of two existing oriented partial differential equation (PDE) filters, including the second-order single oriented PDE filter and the double oriented PDE filter, to two experimentally obtained ESPI speckle fringe patterns with very poor quality, and compare them with other efficient filtering methods, including the adaptive weighted filter, the improved nonlinear complex diffusion PDE, and the windowed Fourier transform method. All of the five filters have been illustrated to be efficient denoising methods through previous comparative analyses in published papers. The experimental results have demonstrated that the two oriented PDE models are applicable to low-quality ESPI speckle fringe patterns. Then for solving the main shortcoming of the two oriented PDE models, we develop the numerically fast algorithms based on Gauss-Seidel strategy for the two oriented PDE models. The proposed numerical algorithms are capable of accelerating the convergence greatly, and perform significantly better in terms of computational efficiency. Our numerically fast algorithms are extended automatically to some other PDE filtering models. PMID:23518722
Numerical Study of Orbits around Europa
NASA Astrophysics Data System (ADS)
Mourao, Decio; Carvalho, Jean Paulo; Vilhena de Moraes, Rodolpho; Cardoso dos Santos, Josué; Campos de Carvalho Costa, Luis Fernando
NASA's Galileo spacecraft probe recently discovered what appears to be a body of liquid water locked inside the icy shell of Jupiter’s moon Europa. The improved likelihood of life on Europa motivated new mission proposals. In this work we used numerical simulations to compare several possible orbits of satellites near the surface of Europa. We spread a set of particles around the satellite with different initially conditions, from 50 to 500km altitude and inclinations higher than 35 degrees, and we monitored the evolution of the test particles during the numerical integrations. We consider the effect of the oblateness of Europa by considering the C22,J2 and J4 parameters and Jupiter gravitational perturbation. These perturbations were first accounted for separately in order to better understand the importance of each effect, and then considered jointly. All particles collide with the Europa surface in a few days. The oblateness of Jupiter alone causes particles with high inclination to collide with the surface of Europa, while the oblateness of Europa affects low orbits decreasing the lifetime of most of the particles. We identified a stable region of orbits with initial altitudes around 300 km of altitude and 90 degrees of inclination. Particles in this region survived more than 200 days. In most of the simulations pericenter initial values near 90 or 270 degrees favor a higher lifetime for the particles, even when considering Europa oblateness.
Numerical Studies of Impurities in Fusion Plasmas
DOE R&D Accomplishments Database
Hulse, R. A.
1982-09-01
The coupled partial differential equations used to describe the behavior of impurity ions in magnetically confined controlled fusion plasmas require numerical solution for cases of practical interest. Computer codes developed for impurity modeling at the Princeton Plasma Physics Laboratory are used as examples of the types of codes employed for this purpose. These codes solve for the impurity ionization state densities and associated radiation rates using atomic physics appropriate for these low-density, high-temperature plasmas. The simpler codes solve local equations in zero spatial dimensions while more complex cases require codes which explicitly include transport of the impurity ions simultaneously with the atomic processes of ionization and recombination. Typical applications are discussed and computational results are presented for selected cases of interest.
Numerical Simulation of Turbulent MHD Flows Using an Iterative PNS Algorithm
NASA Technical Reports Server (NTRS)
Kato, Hiromasa; Tannehill, John C.; Mehta, Unmeel B.
2003-01-01
A new parabolized Navier-Stokes (PNS) algorithm has been developed to efficiently compute magnetohydrodynamic (MHD) flows in the low magnetic Reynolds number regime. In this regime, the electrical conductivity is low and the induced magnetic field is negligible compared to the applied magnetic field. The MHD effects are modeled by introducing source terms into the PNS equation which can then be solved in a very efficient manner. To account for upstream (elliptic) effects, the flowfields are computed using multiple streamwise sweeps with an iterated PNS algorithm. Turbulence has been included by modifying the Baldwin-Lomax turbulence model to account for MHD effects. The new algorithm has been used to compute both laminar and turbulent, supersonic, MHD flows over flat plates and supersonic viscous flows in a rectangular MHD accelerator. The present results are in excellent agreement with previous complete Navier-Stokes calculations.
Structure of the Gabor matrix and efficient numerical algorithms for discrete Gabor expansions
NASA Astrophysics Data System (ADS)
Qiu, Sigang; Feichtinger, Hans G.
1994-09-01
The standard way to obtain suitable coefficients for the (non-orthogonal) Gabor expansion of a general signal for a given Gabor atom g and a pair of lattice constants in the (discrete) time/frequency plane, requires to compute the dual Gabor window function g- first. In this paper, we present an explicit description of the sparsity, the block and banded structure of the Gabor frame matrix G. On this basis efficient algorithms are developed for computing g- by solving the linear equation g- * G equals g with the conjugate- gradients method. Using the dual Gabor wavelet, a fast Gabor reconstruction algorithm with very low computational complexity is proposed.
A numerical algorithm suggested by problems of transport in periodic media - The matrix case.
NASA Technical Reports Server (NTRS)
Allen, R. C., Jr.; Burgmeier, J. W.; Mundorff, P.; Wing, G. M.
1972-01-01
Extension of Allen and Wing's (1970) previous work on problems of transport in periodic media to the matrix case. A method in the form of a complete set of equations is presented that may be used without any further analytical work by investigators interested in computing solutions to problems of the type the method is designed to handle. All the formulas have been checked out numerically, and their effectiveness is demonstrated by numerical examples.
Numerical Studies of Boundary-Layer Receptivity
NASA Technical Reports Server (NTRS)
Reed, Helen L.
1995-01-01
Direct numerical simulations (DNS) of the acoustic receptivity process on a semi-infinite flat plate with a modified-super-elliptic (MSE) leading edge are performed. The incompressible Navier-Stokes equations are solved in stream-function/vorticity form in a general curvilinear coordinate system. The steady basic-state solution is found by solving the governing equations using an alternating direction implicit (ADI) procedure which takes advantage of the parallelism present in line-splitting techniques. Time-harmonic oscillations of the farfield velocity are applied as unsteady boundary conditions to the unsteady disturbance equations. An efficient time-harmonic scheme is used to produce the disturbance solutions. Buffer-zone techniques have been applied to eliminate wave reflection from the outflow boundary. The spatial evolution of Tollmien-Schlichting (T-S) waves is analyzed and compared with experiment and theory. The effects of nose-radius, frequency, Reynolds number, angle of attack, and amplitude of the acoustic wave are investigated. This work is being performed in conjunction with the experiments at the Arizona State University Unsteady Wind Tunnel under the direction of Professor William Saric. The simulations are of the same configuration and parameters used in the wind-tunnel experiments.
Externally fed star formation: a numerical study
NASA Astrophysics Data System (ADS)
Mohammadpour, Motahareh; Stahler, Steven W.
2013-08-01
We investigate, through a series of numerical calculations, the evolution of dense cores that are accreting external gas up to and beyond the point of star formation. Our model clouds are spherical, unmagnetized configurations with fixed outer boundaries, across which gas enters subsonically. When we start with any near-equilibrium state, we find that the cloud's internal velocity also remains subsonic for an extended period, in agreement with observations. However, the velocity becomes supersonic shortly before the star forms. Consequently, the accretion rate building up the protostar is much greater than the benchmark value c_s^3/G, where cs is the sound speed in the dense core. This accretion spike would generate a higher luminosity than those seen in even the most embedded young stars. Moreover, we find that the region of supersonic infall surrounding the protostar races out to engulf much of the cloud, again in violation of the observations, which show infall to be spatially confined. Similar problematic results have been obtained by all other hydrodynamic simulations to date, regardless of the specific infall geometry or boundary conditions adopted. Low-mass star formation is evidently a quasi-static process, in which cloud gas moves inward subsonically until the birth of the star itself. We speculate that magnetic tension in the cloud's deep interior helps restrain the infall prior to this event.
Numerical study of homogeneous nanodroplet growth.
Quang, Tran Si Bui; Leong, Fong Yew; Mirsaidov, Utkur M
2015-01-15
We investigate the axisymmetric homogeneous growth of 10-100 nm water nanodroplets on a substrate surface. The main mechanism of droplet growth is attributed to the accumulation of laterally diffusing water monomers, formed by the absorption of water vapour in the environment onto the substrate. Under assumptions of quasi-steady thermodynamic equilibrium, the nanodroplet evolves according to the augmented Young-Laplace equation. Using continuum theory, we model the dynamics of nanodroplet growth including the coupled effects of disjoining pressure, contact angle and monomer diffusion. Our numerical results show that the initial droplet growth is dominated by monomer diffusion, and the steady late growth rate of droplet radius follows a power law of 1/3, which is unaffected by the substrate disjoining pressure. Instead, the disjoining pressure modifies the growth rate of the droplet height, which then follows a power law of 1/4. We demonstrate how spatial depletion of monomers could lead to a growth arrest of the nanodroplet, as observed experimentally. This work has further implications on the growth kinetics, transport and phase transition of liquids at the nanoscale. PMID:25454424
Numerical studies of 2-dimensional flows
NASA Technical Reports Server (NTRS)
Moretti, G.
1985-01-01
A formulation of the lambda scheme for the analysis of two dimensional inviscid, compressible, unsteady transonic flows is presented. The scheme uses generalized Riemann variables to determine the appropriate two point, one sided finite difference approximation for each derivative in the unsteady Euler equations. These finite differences are applied at the predictor and corrector levels with shock updating at each level. The weaker oblique shocks are captured, but strong near normal shocks are fitted into the flow using the Rankine-Hugoniot relations. This code is demonstrated with a numerical example of a duct flow problem with developing normal and oblique shock waves. The technique is implemented in a code which has been made efficient by streamlining to a minimal number of operations and by eliminating branch statements. The scheme is shown to provide an accurate analysis of the flow, including formation, motions, and interactions of shocks; the results obtained on a relatively coarse mesh are comparable to those obtained by other methods on much finer meshes.
An Approach for Assessing RNA-seq Quantification Algorithms in Replication Studies
Wu, Po-Yen; Phan, John H.; Wang, May D.
2016-01-01
One way to gain a more comprehensive picture of the complex function of a cell is to study the transcriptome. A promising technology for studying the transcriptome is RNA sequencing, an application of which is to quantify elements in the transcriptome and to link quantitative observations to biology. Although numerous quantification algorithms are publicly available, no method of systematically assessing these algorithms has been developed. To meet the need for such an assessment, we present an approach that includes (1) simulated and real datasets, (2) three alignment strategies, and (3) six quantification algorithms. Examining the normalized root-mean-square error, the percentage error of the coefficient of variation, and the distribution of the coefficient of variation, we found that quantification algorithms with the input of sequence alignment reported in the transcriptomic coordinate usually performed better in terms of the multiple metrics proposed in this study.
Numerical studies of constraints and gravitational wave extraction in general relativity
NASA Astrophysics Data System (ADS)
Fiske, David Robert
Within classical physics, general relativity is the theory of gravity. Its equations are non-linear partial differential equations for which relatively few closed form solutions are known. Because of the growing observational need for solutions representing gravitational waves from astrophysically plausible sources, a subfield of general relativity; numerical relativity, has a emerged with the goal of generating numerical solutions to the Einstein equations. This dissertation focuses on two fundamental problems in modern numerical relativity: (1)Creating a theoretical treatment of the constraints in the presence of constraint-violating numerical errors, and (2)Designing and implementing an algorithm to compute the spherical harmonic decomposition of radiation quantities for comparison with observation. On the issue of the constraints, I present a novel and generic procedure for incorporating the constraints into the equations of motion of the theory in a way designed to make the constraint hypersurface an attractor of the evolution. In principle, the prescription generates non- linear corrections for the Einstein equations. The dissertation presents numerical evidence that the correction terms do work in the case of two formulations of the Maxwell equations and two formulations of the linearized Einstein equations. On the issue of radiation extraction, I provide the first in-depth analysis of a novel algorithm, due originally to Misner, for computing spherical harmonic components on a cubic grid. I compute explicitly how the truncation error in the algorithm depends on its various parameters, and I also provide a detailed analysis showing how to implement the method on grids in which explicit symmetries are enforced via boundary conditions. Finally, I verify these error estimates and symmetry arguments with a numerical study using a solution of the linearized Einstein equations known as a Teukolsky wave. The algorithm performs well and the estimates prove true both
Numerical study of fluid motion in bioreactor with two mixers
NASA Astrophysics Data System (ADS)
Zheleva, I.; Lecheva, A.
2015-10-01
Numerical study of hydrodynamic laminar behavior of a viscous fluid in bioreactor with multiple mixers is provided in the present paper. The reactor is equipped with two disk impellers. The fluid motion is studied in stream function-vorticity formulation. The calculations are made by a computer program, written in MATLAB. The fluid structure is described and numerical results are graphically presented and commented.
Numerical study of fluid motion in bioreactor with two mixers
Zheleva, I.; Lecheva, A.
2015-10-28
Numerical study of hydrodynamic laminar behavior of a viscous fluid in bioreactor with multiple mixers is provided in the present paper. The reactor is equipped with two disk impellers. The fluid motion is studied in stream function-vorticity formulation. The calculations are made by a computer program, written in MATLAB. The fluid structure is described and numerical results are graphically presented and commented.
Numerical studies of variable-range hopping in one-dimensional systems
NASA Astrophysics Data System (ADS)
Rodin, A. S.; Fogler, M. M.
2010-03-01
We report on our recent numerical study [1] of hopping transport in disordered one-dimensional systems. A fast new algorithm, based on Dijkstra shortest-path algorithm, is devised to find the lowest-resistance path through the hopping network at arbitrary electric field. Probability distribution functions of individual resistances on the path and the net resistance are calculated and fitted to compact analytic formulas. Qualitative differences between statistics of resistance fluctuations in Ohmic and non-Ohmic regimes are elucidated. The results are compared with prior theoretical and experimental work on the subject.[6pt] [1] A. S. Rodin and M. M. Fogler, Phys. Rev. B 80, 155435 (2009).
Comparative study of skew detection algorithms
NASA Astrophysics Data System (ADS)
Amin, Adnan; Fischer, Stephen; Parkinson, Anthony F.; Shiu, Ricky
1996-10-01
Document image processing has become an increasingly important technology in the automation of office documentation tasks. Automatic document scanners such as text readers and optical character recognition systems are an essential component of systems capable of those tasks. One of the problems in this field is that the document to be read is not always placed correctly on a flat-bed scanner. This means that the document may be skewed on the scanner bed, resulting in a skewed image. This skew has a detrimental effect on document analysis, document understanding, and character segmentation and recognition. Consequently, detecting the skew of a document image and correcting it are important issues in realizing a practical document reader. We describe a new algorithm for skew detection. We then compare the performance and results of this skew detection algorithm to other published methods from O'Gorman, Hinds, Le, Baird, Postl, and Akiyama. Finally, we discuss the theory of skew detection and the different approaches taken to solve the problem of skew in documents. The skew correction algorithm we propose has been shown to be extremely fast, with run times averaging under 0.25 CPU seconds to calculate the angle on a DEC 5000/20 workstation.
Kisselev, V B; Roberti, L; Perona, G
1995-12-20
The recently developed finite-element method for solution of the radiative transfer equation has been extended to compute the full azimuthal dependence of the radiance in a vertically inhomogeneous plane-parallel medium. The physical processes that are included in the algorithm are multiple scattering and bottom boundary bidirectional reflectivity. The incident radiation is a parallel flux on the top boundary that is characteristic for illumination of the atmosphere by the Sun in the UV, visible, and near-infrared regions of the electromagnetic spectrum. The theoretical basis is presented together with a number of applications to realistic atmospheres. The method is shown to be accurate even with a low number of grid points for most of the considered situations. The FORTRAN code for this algorithm is developed and is available for applications. PMID:21068966
NASA Astrophysics Data System (ADS)
Kisselev, Viatcheslav B.; Roberti, Laura; Perona, Giovanni
1995-12-01
The recently developed finite-element method for solution of the radiative transfer equation has been extended to compute the full azimuthal dependence of the radiance in a vertically inhomogeneous plane-parallel medium. The physical processes that are included in the algorithm are multiple scattering and bottom boundary bidirectional reflectivity. The incident radiation is a parallel flux on the top boundary that is characteristic for illumination of the atmosphere by the Sun in the UV, visible, and near-infrared regions of the electromagnetic spectrum. The theoretical basis is presented together with a number of applications to realistic atmospheres. The method is shown to be accurate even with a low number of grid points for most of the considered situations. The fortran code for this algorithm is developed and is available for applications.
Parallel technology for numerical modeling of fluid dynamics problems by high-accuracy algorithms
NASA Astrophysics Data System (ADS)
Gorobets, A. V.
2015-04-01
A parallel computation technology for modeling fluid dynamics problems by finite-volume and finite-difference methods of high accuracy is presented. The development of an algorithm, the design of a software implementation, and the creation of parallel programs for computations on large-scale computing systems are considered. The presented parallel technology is based on a multilevel parallel model combining various types of parallelism: with shared and distributed memory and with multiple and single instruction streams to multiple data flows.
Scanning of wind turbine upwind conditions: numerical algorithm and first applications
NASA Astrophysics Data System (ADS)
Calaf, Marc; Cortina, Gerard; Sharma, Varun; Parlange, Marc B.
2014-11-01
Wind turbines still obtain in-situ meteorological information by means of traditional wind vane and cup anemometers installed at the turbine's nacelle, right behind the blades. This has two important drawbacks: 1-turbine misalignment with the mean wind direction is common and energy losses are experienced; 2-the near-blade monitoring does not provide any time to readjust the profile of the wind turbine to incoming turbulence gusts. A solution is to install wind Lidar devices on the turbine's nacelle. This technique is currently under development as an alternative to traditional in-situ wind anemometry because it can measure the wind vector at substantial distances upwind. However, at what upwind distance should they interrogate the atmosphere? A new flexible wind turbine algorithm for large eddy simulations of wind farms that allows answering this question, will be presented. The new wind turbine algorithm timely corrects the turbines' yaw misalignment with the changing wind. The upwind scanning flexibility of the algorithm also allows to track the wind vector and turbulent kinetic energy as they approach the wind turbine's rotor blades. Results will illustrate the spatiotemporal evolution of the wind vector and the turbulent kinetic energy as the incoming flow approaches the wind turbine under different atmospheric stability conditions. Results will also show that the available atmospheric wind power is larger during daytime periods at the cost of an increased variance.
AN ACCURATE AND EFFICIENT ALGORITHM FOR NUMERICAL SIMULATION OF CONDUCTION-TYPE PROBLEMS. (R824801)
A modification of the finite analytic numerical method for conduction-type (diffusion) problems is presented. The finite analytic discretization scheme is derived by means of the Fourier series expansion for the most general case of nonuniform grid and variabl...
Systematic analytical and numerical studies of highly correlated electron systems
NASA Astrophysics Data System (ADS)
Tsai, Shan-Wen
Strong electron correlations in condensed matter systems give rise to a wide range of striking physical properties, producing phenomena as varied as high temperature superconductivity, metal-insulator transitions and the integer and fractional quantum Hall effects. Quantum critical systems also exhibit strong correlations between a large number of degrees of freedom. In this thesis we study these complicated systems using a combination of analytical and numerical approaches. We perform systematic investigations, which adds to the robustness of our results. We develop a new method, based on the density-matrix renormalization-group (DMRG) algorithm combined with finite-size scaling analysis, to study critical behavior in quantum spin chains and extract critical exponents. Accurate results are obtained for spin-1/2 antiferromagnetic chains and the spin-1 chain at the critical point separating the Haldane and the dimerized phases. Disorder in a system can change its properties drastically. Plateau transitions in the integer quantum Hall effect provide the clearest example of quantum critical behavior in a disordered system. We provide analytical proof that the Chalker-Coddington model, which is used to describe the plateau transitions, is quantum critical. Starting from a field theory based on this model, equivalent to a non-Hermitian supersymmetric spin chain, we prove quantum criticality by a Lieb-Schultz-Mattis type theorem. This approach was motivated by numerical results obtained using the DMRG/finite-size scaling method. Our generalized LSM theorem also applies to the spin quantum Hall effect, which can appear in disordered d-wave superconductors with broken time-reversal symmetry. The last part of the thesis is a renormalization-group study of two dimensional interacting electron systems. We obtain results relevant to high-temperature superconductors and also to the family of kappa - (BEDT - TTF)2X organic superconductors. At half filling, the fully nested
NASA Astrophysics Data System (ADS)
Wang, Jiong; Steinmann, Paul
2016-05-01
This is part II of this series of papers. The aim of the current paper was to solve the governing PDE system derived in part I numerically, such that the procedure of variant reorientation in a magnetic shape memory alloy (MSMA) sample can be simulated. The sample to be considered in this paper has a 3D cuboid shape and is subject to typical magnetic and mechanical loading conditions. To investigate the demagnetization effect on the sample's response, the surrounding space of the sample is taken into account. By considering the different properties of the independent variables, an iterative numerical algorithm is proposed to solve the governing system. The related mathematical formulas and some techniques facilitating the numerical calculations are introduced. Based on the results of numerical simulations, the distributions of some important physical quantities (e.g., magnetization, demagnetization field, and mechanical stress) in the sample can be determined. Furthermore, the properties of configurational force on the twin interfaces are investigated. By virtue of the twin interface movement criteria derived in part I, the whole procedure of magnetic field- or stress-induced variant reorientations in the MSMA sample can be properly simulated.
Numerical study of time domain analogy applied to noise prediction from rotating blades
NASA Astrophysics Data System (ADS)
Fedala, D.; Kouidri, S.; Rey, R.
2009-04-01
Aeroacoustic formulations in time domain are frequently used to model the aerodynamic sound of airfoils, the time data being more accessible. The formulation 1A developed by Farassat, an integral solution of the Ffowcs Williams and Hawkings equation, holds great interest because of its ability to handle surfaces in arbitrary motion. The aim of this work is to study the numerical sensitivity of this model to specified parameters used in the calculation. The numerical algorithms, spatial and time discretizations, and approximations used for far-field acoustic simulation are presented. An approach of quantifying of the numerical errors resulting from implementation of formulation 1A is carried out based on Isom's and Tam's test cases. A helicopter blade airfoil, as defined by Farassat to investigate Isom's case, is used in this work. According to Isom, the acoustic response of a dipole source with a constant aerodynamic load, ρ0c02, is equal to the thickness noise contribution. Discrepancies are observed when the two contributions are computed numerically. In this work, variations of these errors, which depend on the temporal resolution, Mach number, source-observer distance, and interpolation algorithm type, are investigated. The results show that the spline interpolating algorithm gives the minimum error. The analysis is then extended to Tam's test case. Tam's test case has the advantage of providing an analytical solution for the first harmonic of the noise produced by a specific force distribution.
A Comparative Study of Interferometric Regridding Algorithms
NASA Technical Reports Server (NTRS)
Hensley, Scott; Safaeinili, Ali
1999-01-01
THe paper discusses regridding options: (1) The problem of interpolating data that is not sampled on a uniform grid, that is noisy, and contains gaps is a difficult problem. (2) Several interpolation algorithms have been implemented: (a) Nearest neighbor - Fast and easy but shows some artifacts in shaded relief images. (b) Simplical interpolator - uses plane going through three points containing point where interpolation is required. Reasonably fast and accurate. (c) Convolutional - uses a windowed Gaussian approximating the optimal prolate spheroidal weighting function for a specified bandwidth. (d) First or second order surface fitting - Uses the height data centered in a box about a given point and does a weighted least squares surface fit.
NASA Technical Reports Server (NTRS)
Carter, Richard G.
1989-01-01
For optimization problems associated with engineering design, parameter estimation, image reconstruction, and other optimization/simulation applications, low accuracy function and gradient values are frequently much less expensive to obtain than high accuracy values. Here, researchers investigate the computational performance of trust region methods for nonlinear optimization when high accuracy evaluations are unavailable or prohibitively expensive, and confirm earlier theoretical predictions when the algorithm is convergent even with relative gradient errors of 0.5 or more. The proper choice of the amount of accuracy to use in function and gradient evaluations can result in orders-of-magnitude savings in computational cost.
An Effective Hybrid Firefly Algorithm with Harmony Search for Global Numerical Optimization
Guo, Lihong; Wang, Gai-Ge; Wang, Heqi; Wang, Dinan
2013-01-01
A hybrid metaheuristic approach by hybridizing harmony search (HS) and firefly algorithm (FA), namely, HS/FA, is proposed to solve function optimization. In HS/FA, the exploration of HS and the exploitation of FA are fully exerted, so HS/FA has a faster convergence speed than HS and FA. Also, top fireflies scheme is introduced to reduce running time, and HS is utilized to mutate between fireflies when updating fireflies. The HS/FA method is verified by various benchmarks. From the experiments, the implementation of HS/FA is better than the standard FA and other eight optimization methods. PMID:24348137
NASA Technical Reports Server (NTRS)
Gunzburger, M. D.; Nicolaides, R. A.
1986-01-01
Substructuring methods are in common use in mechanics problems where typically the associated linear systems of algebraic equations are positive definite. Here these methods are extended to problems which lead to nonpositive definite, nonsymmetric matrices. The extension is based on an algorithm which carries out the block Gauss elimination procedure without the need for interchanges even when a pivot matrix is singular. Examples are provided wherein the method is used in connection with finite element solutions of the stationary Stokes equations and the Helmholtz equation, and dual methods for second-order elliptic equations.
Numerical Study of the Structure of Metastable Configurations for the Thomson Problem
NASA Astrophysics Data System (ADS)
Bondarenko, A. N.; Bugueva, T. V.; Kozinkin, L. A.
2016-05-01
A numerical method is proposed for solving the Thomson problem - finding stable positions for a system of N point charges distributed on a sphere that minimize the potential energy of the system. The behavior of this system is essentially nonlinear, and the number of metastable structures grows exponentially with N. This makes the problem of finding all stable configurations extremely difficult. The results of testing of the developed algorithm and of numerical study of the properties of the local potential energy minima for a system of point charges are presented.
NASA Astrophysics Data System (ADS)
Hu, P.; Shi, D. Y.; Ying, L.; Shen, G. Z.; Chang, Y.; Liu, W. Q.
2013-05-01
Thermal-mechanical-transformation coupled theoretical model for hot stamping and rheological behavior of high strength steel at elevated temperatures were obtained through non-isothermal and isothermal tensile tests respectively in this work. The static explicit finite element equations for hot stamping were proposed based on thermal-mechanical-transformation coupled constitutive laws and nonlinear, large deformation analysis. According to these equations, the hot stamping module of KMAS (King Mesh Analysis System) was developed for the numerical simulation of sheet metal forming at elevated temperatures. Afterwards, the hot stamping simulation of a typical B-pillar conducted by the KMAS software was compared to the experiment. The comparison consists of the following sides: temperature distribution, thickness distribution and martensite fraction. The good agreement between numerical simulation and the experiment confirms that the multi-field coupled constitutive laws and the KMAS software can predict hot stamping process accurately.
A Numerical Study of Nonlinear Wave Interactions
NASA Astrophysics Data System (ADS)
de Bakker, A.; Tissier, M.; Ruessink, G.
2014-12-01
Nonlinear triad interactions redistribute energy among a wave field, which transforms the shape of the incident short waves (f = 0.05 - 2 Hz) and generates energy at infragravity frequencies (f = 0.005-0.05 Hz). Recently, it has been suggested that infragravity energy may dissipate by energy transfers from infragravity frequencies to either the (former) short-wave spectral peak, or through infragravity-infragravity self-interactions that cause the infragravity waves to steepen and to eventually break. To investigate these infragravity dissipation mechanisms, we use the non-hydrostatic SWASH model. In this study, we first validate the model with the high-resolution GLOBEX laboratory data set and then explore the dependence of the energy transfers, with a focus on infragravity frequencies, on beach slope. Consistent with previous studies we find that SWASH is able to reproduce the transformation and corresponding nonlinear energy transfers of shoreward propagating waves to great detail. Bispectral analysis is used to study the coupling between wave frequencies; nonlinear energy transfers are then quantified using the Boussinesq coupling coefficient. To obtain more detailed insight we divide the nonlinear interactions in four categories based on triads including 1) infragravity frequencies only, 2) two infragravity frequencies and one short-wave frequency, 3) one infragravity frequency and two short-wave frequencies and 4) short-wave frequencies only. Preliminary results suggest that interactions are rather weak on gently beach slopes (1:80) and, in the innermost part of the surf zone, are dominated by infragravity-infragravity interactions. On steeper slopes (1:20), interactions are stronger, but entirely dominated by those involving short-wave frequencies only. The dependence of the transfers on offshore wave conditions and beach shape will be explored too. Funded by NWO.
Stamnes, K; Tsay, S C; Wiscombe, W; Jayaweera, K
1988-06-15
We summarize an advanced, thoroughly documented, and quite general purpose discrete ordinate algorithm for time-independent transfer calculations in vertically inhomogeneous, nonisothermal, plane-parallel media. Atmospheric applications ranging from the UV to the radar region of the electromagnetic spectrum are possible. The physical processes included are thermal emission, scattering, absorption, and bidirectional reflection and emission at the lower boundary. The medium may be forced at the top boundary by parallel or diffuse radiation and by internal and boundary thermal sources as well. We provide a brief account of the theoretical basis as well as a discussion of the numerical implementation of the theory. The recent advances made by ourselves and our collaborators-advances in both formulation and numerical solution-are all incorporated in the algorithm. Prominent among these advances are the complete conquest of two illconditioning problems which afflicted all previous discrete ordinate implementations: (1) the computation of eigenvalues and eigenvectors and (2) the inversion of the matrix determining the constants of integration. Copies of the FORTRAN program on microcomputer diskettes are available for interested users. PMID:20531783
NASA Astrophysics Data System (ADS)
Zhang, Lisha
We present fast and robust numerical algorithms for 3-D scattering from perfectly electrical conducting (PEC) and dielectric random rough surfaces in microwave remote sensing. The Coifman wavelets or Coiflets are employed to implement Galerkin's procedure in the method of moments (MoM). Due to the high-precision one-point quadrature, the Coiflets yield fast evaluations of the most off-diagonal entries, reducing the matrix fill effort from O(N2) to O( N). The orthogonality and Riesz basis of the Coiflets generate well conditioned impedance matrix, with rapid convergence for the conjugate gradient solver. The resulting impedance matrix is further sparsified by the matrix-formed standard fast wavelet transform (SFWT). By properly selecting multiresolution levels of the total transformation matrix, the solution precision can be enhanced while matrix sparsity and memory consumption have not been noticeably sacrificed. The unified fast scattering algorithm for dielectric random rough surfaces can asymptotically reduce to the PEC case when the loss tangent grows extremely large. Numerical results demonstrate that the reduced PEC model does not suffer from ill-posed problems. Compared with previous publications and laboratory measurements, good agreement is observed.
A Numerical Study of Superconducting Cavity Components
B.C. Yunn; J.J. Bisognano
1990-09-10
Computer programs which solve Maxwell's equations in three dimensions are becoming an invaluable tool in the design of RF structures for particle accelerators. In particular, the lack of cylindrical symmetry of superconducting cavities with waveguide couplers demands a 3-D analysis for a reasonable description of a number of important phenomena. A set of codes, collectively known as MAFIA, developed by Weiland and his collaborators, has been used at CEBAF to study its five-cell superconducting accelerating cavities. The magnitude of RF crosstalk between cavities is found to depend critically on the breaking of cylindrical symmetry by the fundamental power couplers. A model of the higher order mode coupler exhibits an unexpected mode which is in good agreement with measurement.
NASA Astrophysics Data System (ADS)
Harries, Tim J.
2015-04-01
We present a set of new numerical methods that are relevant to calculating radiation pressure terms in hydrodynamics calculations, with a particular focus on massive star formation. The radiation force is determined from a Monte Carlo estimator and enables a complete treatment of the detailed microphysics, including polychromatic radiation and anisotropic scattering, in both the free-streaming and optically thick limits. Since the new method is computationally demanding we have developed two new methods that speed up the algorithm. The first is a photon packet splitting algorithm that enables efficient treatment of the Monte Carlo process in very optically thick regions. The second is a parallelization method that distributes the Monte Carlo workload over many instances of the hydrodynamic domain, resulting in excellent scaling of the radiation step. We also describe the implementation of a sink particle method that enables us to follow the accretion on to, and the growth of, the protostars. We detail the results of extensive testing and benchmarking of the new algorithms.
Numerical study of insect free hovering flight
NASA Astrophysics Data System (ADS)
Wu, Di; Yeo, Khoon Seng; Lim, Tee Tai; Fluid lab, Mechanical Engineering, National University of Singapore Team
2012-11-01
In this paper we present the computational fluid dynamics study of three-dimensional flow field around a free hovering fruit fly integrated with unsteady FSI analysis and the adaptive flight control system for the first time. The FSI model being specified for fruitfly hovering is achieved by coupling a structural problem based on Newton's second law with a rigorous CFD solver concerning generalized finite difference method. In contrast to the previous hovering flight research, the wing motion employed here is not acquired from experimental data but governed by our proposed control systems. Two types of hovering control strategies i.e. stroke plane adjustment mode and paddling mode are explored, capable of generating the fixed body position and orientation characteristic of hovering flight. Hovering flight associated with multiple wing kinematics and body orientations are shown as well, indicating the means by which fruitfly actually maintains hovering may have considerable freedom and therefore might be influenced by many other factors beyond the physical and aerodynamic requirements. Additionally, both the near- and far-field flow and vortex structure agree well with the results from other researchers, demonstrating the reliability of our current model.
A numerical study of aircraft empennage buffet
NASA Astrophysics Data System (ADS)
Findlay, David Bruce
1999-10-01
A method to predict tightly-coupled dynamic aeroelastic vertical tail buffet was presented. Analysis of high angle of attack vertical tail buffet was performed. A Navier-Stokes fluid dynamics method was coupled with a modal structural dynamics method. The approach was to improve upon existing methods to evaluate complex geometric arrangements with general multi-zone interfacing. The method was demonstrated through a step- wise approach beginning with a simple configuration and building up to a complete aircraft at high angle of attack with flexible tail surfaces. Results compared well with in-flight and Full-scale wind tunnel measured trends and frequency content. Comparisons with measured absolute values of buffet loads showed the computations to be under-predicting the test data. This was primarily attributed to insufficient grid resolution, in particular in the vicinity of the main vortex flow. The demanding computational requirements of full-configuration tail buffet prediction limited the fidelity. The primary contribution of the present study was the extension and demonstration of a tightly-coupled aeroelastic computational fluid dynamics/structural dynamics based analysis method for analysis of aircraft empennage buffet. The focus was on improving the development process associated with characterizing empennage buffet loads and the resulting structural response. The intent was to establish a computationally based alternative approach to the experimentally based process currently employed. The computational method was employed to provide far greater insight into the flow physics phenomena associated with specific configurations and conditions of interest.
Numerical Study of Explosive Dispersal of Particles
NASA Astrophysics Data System (ADS)
Rollin, Bertrand; Annamalai, Subramanian; Neal, Christopher; Jackson, Thomas; Balachandar, S.
2014-11-01
Recent experiments have shown that when a layer of solid particles is explosively dispersed, a multiphase instability front occurs, which leads to the formation of aerodynamically stable jet-like particle structures. We aim at replicating these experimental observations using highly resolved large-scale simulations, to improve our understanding of particulate front instabilities and jetting phenomenon. We consider a cylindrical core of high pressure and density gas generated from energetic material. Throughout the length of the cylinder, an annular region of micron-sized inert spherical particles surrounds the charge. The particles are treated as point particles, the gas is treated as a continuum, and a rigorous two-way coupled compressible multiphase formulation is used. The jets are believed to have their origin during the early phase of rapid acceleration of the bed of particles. Therefore, this work focuses on capturing the early-time behavior and growth of the instabilities caused by the presence of particles. The accuracy of our predictive simulations will be studied by comparing the shock radius, particle front location, and other relevant metrics against the data extracted from experimental results. This work is supported by the U.S. DoE, National Nuclear Security Administration, Advanced Simulation and Computing Program, as a Cooperative Agreement under the Predictive Science Academic Alliance Program, under Contract No. DE-NA0002378.
Numerical Studies of Properties of Confined Helium
NASA Technical Reports Server (NTRS)
Manousakis, Efstratios
2003-01-01
We carry out state of the art simulations of properties of confined liquid helium near the superfluid transition to a degree of accuracy which allows to make predictions for the outcome of fundamental physics experiments in microgravity. First we report our results for the finite-size scaling behavior of heat capacity of superfluids for cubic and parallel-plate geometry. This allows us to study the crossover from zero and two dimensions to three dimensions. Our calculated scaling functions are in good agreement with recently measured specific heat scaling functions for the above mentioned geometries. We also present our results of a quantum simulation of submonolayer of molecular hydrogen deposited on an ideal graphite substrate using path-integral quantum Monte Carlo simulation. We find that the monolayer phase diagram is rich and very similar to that of helium monolayer. We are able to uncover the main features of the complex monolayer phase diagram, such as the commensurate solid phases and the commensurate to incommensurate transition, in agreement with the experiments and to find some features which are missing from the experimental analysis.
Study of image matching algorithm and sub-pixel fitting algorithm in target tracking
NASA Astrophysics Data System (ADS)
Yang, Ming-dong; Jia, Jianjun; Qiang, Jia; Wang, Jian-yu
2015-03-01
Image correlation matching is a tracking method that searched a region most approximate to the target template based on the correlation measure between two images. Because there is no need to segment the image, and the computation of this method is little. Image correlation matching is a basic method of target tracking. This paper mainly studies the image matching algorithm of gray scale image, which precision is at sub-pixel level. The matching algorithm used in this paper is SAD (Sum of Absolute Difference) method. This method excels in real-time systems because of its low computation complexity. The SAD method is introduced firstly and the most frequently used sub-pixel fitting algorithms are introduced at the meantime. These fitting algorithms can't be used in real-time systems because they are too complex. However, target tracking often requires high real-time performance, we put forward a fitting algorithm named paraboloidal fitting algorithm based on the consideration above, this algorithm is simple and realized easily in real-time system. The result of this algorithm is compared with that of surface fitting algorithm through image matching simulation. By comparison, the precision difference between these two algorithms is little, it's less than 0.01pixel. In order to research the influence of target rotation on precision of image matching, the experiment of camera rotation was carried on. The detector used in the camera is a CMOS detector. It is fixed to an arc pendulum table, take pictures when the camera rotated different angles. Choose a subarea in the original picture as the template, and search the best matching spot using image matching algorithm mentioned above. The result shows that the matching error is bigger when the target rotation angle is larger. It's an approximate linear relation. Finally, the influence of noise on matching precision was researched. Gaussian noise and pepper and salt noise were added in the image respectively, and the image
Dynamical Approach Study of Spurious Numerics in Nonlinear Computations
NASA Technical Reports Server (NTRS)
Yee, H. C.; Mansour, Nagi (Technical Monitor)
2002-01-01
The last two decades have been an era when computation is ahead of analysis and when very large scale practical computations are increasingly used in poorly understood multiscale complex nonlinear physical problems and non-traditional fields. Ensuring a higher level of confidence in the predictability and reliability (PAR) of these numerical simulations could play a major role in furthering the design, understanding, affordability and safety of our next generation air and space transportation systems, and systems for planetary and atmospheric sciences, and in understanding the evolution and origin of life. The need to guarantee PAR becomes acute when computations offer the ONLY way of solving these types of data limited problems. Employing theory from nonlinear dynamical systems, some building blocks to ensure a higher level of confidence in PAR of numerical simulations have been revealed by the author and world expert collaborators in relevant fields. Five building blocks with supporting numerical examples were discussed. The next step is to utilize knowledge gained by including nonlinear dynamics, bifurcation and chaos theories as an integral part of the numerical process. The third step is to design integrated criteria for reliable and accurate algorithms that cater to the different multiscale nonlinear physics. This includes but is not limited to the construction of appropriate adaptive spatial and temporal discretizations that are suitable for the underlying governing equations. In addition, a multiresolution wavelets approach for adaptive numerical dissipation/filter controls for high speed turbulence, acoustics and combustion simulations will be sought. These steps are corner stones for guarding against spurious numerical solutions that are solutions of the discretized counterparts but are not solutions of the underlying governing equations.
NASA Astrophysics Data System (ADS)
Lakshminarayana, B.; Ho, Y.; Basson, A.
1993-07-01
The objective of this research is to simulate steady and unsteady viscous flows, including rotor/stator interaction and tip clearance effects in turbomachinery. The numerical formulation for steady flow developed here includes an efficient grid generation scheme, particularly suited to computational grids for the analysis of turbulent turbomachinery flows and tip clearance flows, and a semi-implicit, pressure-based computational fluid dynamics scheme that directly includes artificial dissipation, and is applicable to both viscous and inviscid flows. The values of these artificial dissipation is optimized to achieve accuracy and convergency in the solution. The numerical model is used to investigate the structure of tip clearance flows in a turbine nozzle. The structure of leakage flow is captured accurately, including blade-to-blade variation of all three velocity components, pitch and yaw angles, losses and blade static pressures in the tip clearance region. The simulation also includes evaluation of such quantities of leakage mass flow, vortex strength, losses, dominant leakage flow regions and the spanwise extent affected by the leakage flow. It is demonstrated, through optimization of grid size and artificial dissipation, that the tip clearance flow field can be captured accurately. The above numerical formulation was modified to incorporate time accurate solutions. An inner loop iteration scheme is used at each time step to account for the non-linear effects. The computation of unsteady flow through a flat plate cascade subjected to a transverse gust reveals that the choice of grid spacing and the amount of artificial dissipation is critical for accurate prediction of unsteady phenomena. The rotor-stator interaction problem is simulated by starting the computation upstream of the stator, and the upstream rotor wake is specified from the experimental data. The results show that the stator potential effects have appreciable influence on the upstream rotor wake
NASA Technical Reports Server (NTRS)
Lakshminarayana, B.; Ho, Y.; Basson, A.
1993-01-01
The objective of this research is to simulate steady and unsteady viscous flows, including rotor/stator interaction and tip clearance effects in turbomachinery. The numerical formulation for steady flow developed here includes an efficient grid generation scheme, particularly suited to computational grids for the analysis of turbulent turbomachinery flows and tip clearance flows, and a semi-implicit, pressure-based computational fluid dynamics scheme that directly includes artificial dissipation, and is applicable to both viscous and inviscid flows. The values of these artificial dissipation is optimized to achieve accuracy and convergency in the solution. The numerical model is used to investigate the structure of tip clearance flows in a turbine nozzle. The structure of leakage flow is captured accurately, including blade-to-blade variation of all three velocity components, pitch and yaw angles, losses and blade static pressures in the tip clearance region. The simulation also includes evaluation of such quantities of leakage mass flow, vortex strength, losses, dominant leakage flow regions and the spanwise extent affected by the leakage flow. It is demonstrated, through optimization of grid size and artificial dissipation, that the tip clearance flow field can be captured accurately. The above numerical formulation was modified to incorporate time accurate solutions. An inner loop iteration scheme is used at each time step to account for the non-linear effects. The computation of unsteady flow through a flat plate cascade subjected to a transverse gust reveals that the choice of grid spacing and the amount of artificial dissipation is critical for accurate prediction of unsteady phenomena. The rotor-stator interaction problem is simulated by starting the computation upstream of the stator, and the upstream rotor wake is specified from the experimental data. The results show that the stator potential effects have appreciable influence on the upstream rotor wake
NASA Astrophysics Data System (ADS)
Kuraz, Michal
2016-06-01
This paper presents pseudo-deterministic catchment runoff model based on the Richards equation model [1] - the governing equation for the subsurface flow. The subsurface flow in a catchment is described here by two-dimensional variably saturated flow (unsaturated and saturated). The governing equation is the Richards equation with a slight modification of the time derivative term as considered e.g. by Neuman [2]. The nonlinear nature of this problem appears in unsaturated zone only, however the delineation of the saturated zone boundary is a nonlinear computationally expensive issue. The simple one-dimensional Boussinesq equation was used here as a rough estimator of the saturated zone boundary. With this estimate the dd-adaptivity algorithm (see Kuraz et al. [4, 5, 6]) could always start with an optimal subdomain split, so it is now possible to avoid solutions of huge systems of linear equations in the initial iteration level of our Richards equation based runoff model.
Gas-kinetic numerical studies of three-dimensional complex flows on spacecraft re-entry
NASA Astrophysics Data System (ADS)
Li, Zhi-Hui; Zhang, Han-Xin
2009-03-01
The gas-kinetic numerical algorithm solving the Boltzmann model equation is extended and developed to study the three-dimensional hypersonic flows of spacecraft re-entry into the atmosphere in perfect gas. In this study, the simplified velocity distribution function equation for various flow regimes is presented on the basis of the kinetic Boltzmann-Shakhov model. The discrete velocity ordinate technique and numerical quadrature methods, such as the Gauss quadrature formulas with the weight function 2/ π1/2exp(- V2) and the Gauss-Legendre numerical quadrature rule, are studied to resolve the barrier in simulating complex flows from low Mach numbers to hypersonic problems. Specially, the gas-kinetic finite-difference scheme is constructed for the computation of three-dimensional flow problems, which directly captures the time evolution of the molecular velocity distribution function. The gas-kinetic boundary conditions and numerical procedures are studied and implemented by directly acting on the velocity distribution function. The HPF (high performance fortran) parallel implementation technique for the gas-kinetic numerical method is developed and applied to study the hypersonic flows around three-dimensional complex bodies. The main purpose of the current research is to provide a way to extend the gas-kinetic numerical algorithm to the flow computation of three-dimensional complex hypersonic problems with high Mach numbers. To verify the current method and simulate gas transport phenomena covering various flow regimes, the three-dimensional hypersonic flows around sphere and spacecraft shape with different Knudsen numbers and Mach numbers are studied by HPF parallel computing. Excellent results have been obtained for all examples computed.
NASA Astrophysics Data System (ADS)
Angeli, D.; Stalio, E.; Corticelli, M. A.; Barozzi, G. S.
2015-11-01
A parallel algorithm is presented for the Direct Numerical Simulation of buoyancy- induced flows in open or partially confined periodic domains, containing immersed cylindrical bodies of arbitrary cross-section. The governing equations are discretized by means of the Finite Volume method on Cartesian grids. A semi-implicit scheme is employed for the diffusive terms, which are treated implicitly on the periodic plane and explicitly along the homogeneous direction, while all convective terms are explicit, via the second-order Adams-Bashfort scheme. The contemporary solution of velocity and pressure fields is achieved by means of a projection method. The numerical resolution of the set of linear equations resulting from discretization is carried out by means of efficient and highly parallel direct solvers. Verification and validation of the numerical procedure is reported in the paper, for the case of flow around an array of heated cylindrical rods arranged in a square lattice. Grid independence is assessed in laminar flow conditions, and DNS results in turbulent conditions are presented for two different grids and compared to available literature data, thus confirming the favorable qualities of the method.
A critical evaluation of numerical algorithms and flow physics in complex supersonic flows
NASA Astrophysics Data System (ADS)
Aradag, Selin
In this research, two different complex supersonic flows are selected to apply CFD to Navier-Stokes simulations. First test case is "Supersonic Flow over an Open Rectangular Cavity". Open cavity flow fields are remarkably complicated with internal and external regions that are coupled via self-sustained shear layer oscillations. Supersonic flow past a cavity has numerous applications in store carriage and release. Internal carriage of stores, which can be modeled using a cavity configuration, is used for supersonic aircraft in order to reduce radar cross section, aerodynamic drag and aerodynamic heating. Supersonic, turbulent, three-dimensional unsteady flow past an open rectangular cavity is simulated, to understand the physics and three-dimensional nature of the cavity flow oscillations. Influences of numerical parameters such as numerical flux scheme, computation time and flux limiter on the computed flow are determined. Two dimensional simulations are also performed for comparison purposes. The next test case is "The Computational Design of Boeing/AFOSR Mach 6 Wind Tunnel". Due to huge differences between geometrical scales, this problem is both challenging and computationally intensive. It is believed that most of the experimental data obtained from conventional ground testing facilities are not reliable due to high levels of noise associated with the acoustic fluctuations from the turbulent boundary layers on the wind tunnel walls. Therefore, it is very important to have quiet testing facilities for hypersonic flow research. The Boeing/AFOSR Mach 6 Wind tunnel in Purdue University has been designed as a quiet tunnel for which the noise level is an order of magnitude lower than that in conventional wind tunnels. However, quiet flow is achieved in the Purdue Mach 6 tunnel for only low Reynolds numbers. Early transition of the nozzle wall boundary layer has been identified as the cause of the test section noise. Separation bubbles on the bleed lip and associated
Numerical study of free surface flow around large obstacles
NASA Astrophysics Data System (ADS)
Zhang, Yanming
In this thesis a numerical model was developed to study three-dimensional turbulent flows around large obstacles in an open channel. With this numerical model, a series of numerical tests was carried out, and the properties of turbulent flows around a single obstacle or a cluster of obstacles were investigated. The origin of this study was to study the flow properties around fish habitat structures. Actually, the numerical model can be applied to the study of general turbulent flows under free surfaces. In the numerical model the three-dimensional Reynolds-averaged Navier-Stokes equations in conjunction with k-epsilon turbulence model were solved in a free surface fitted coordinate system. First, different forms of governing equations for turbulent flow were investigated, and a concise form of fully transformed governing equations in a general curvilinear coordinate system was derived. In the numerical solution the FAVOR (Fractional Area/Volume Obstacle Representation) technique was extended into the free surface fitted coordinate system. With this feature the problem of complex turbulent flow with a free surface and general shaped obstacles could be solved efficiently. To locate the free surface, a method based on integrating the momentum equation in the vertical direction was developed. After study and tests of several popular difference schemes, a QUICK scheme with UMIST limiter was adopted in this numerical model. Several test cases were presented to demonstrate the present numerical model. The first test case was to simulate a submerged hydraulic jump. The calculated velocity, free surface profile and turbulence properties of the flow showed a close match with the experimental data. The second test was a submerged hydraulic jump with a baffle sill. The comparison between numerical and experimental data indicated that the current numerical model could catch the general flow structures of the submerged hydraulic jumps. The last two test cases were flows around a
NASA Astrophysics Data System (ADS)
Li, Cong; Lei, Jianshe
2014-10-01
In this paper, we focus on the influences of various parameters in the niching genetic algorithm inversion procedure on the results, such as various objective functions, the number of the models in each subpopulation, and the critical separation radius. The frequency-waveform integration (F-K) method is applied to synthesize three-component waveform data with noise in various epicentral distances and azimuths. Our results show that if we use a zero-th-lag cross-correlation function, then we will obtain the model with a faster convergence and a higher precision than other objective functions. The number of models in each subpopulation has a great influence on the rate of convergence and computation time, suggesting that it should be obtained through tests in practical problems. The critical separation radius should be determined carefully because it directly affects the multi-extreme values in the inversion. We also compare the inverted results from full-band waveform data and surface-wave frequency-band (0.02-0.1 Hz) data, and find that the latter is relatively poorer but still has a higher precision, suggesting that surface-wave frequency-band data can also be used to invert for the crustal structure.
NASA Astrophysics Data System (ADS)
Li, Cong; Lei, Jianshe
2014-09-01
In this paper, we focus on the influences of various parameters in the niching genetic algorithm inversion procedure on the results, such as various objective functions, the number of the models in each subpopulation, and the critical separation radius. The frequency-waveform integration (F-K) method is applied to synthesize three-component waveform data with noise in various epicentral distances and azimuths. Our results show that if we use a zero-th-lag cross-correlation function, then we will obtain the model with a faster convergence and a higher precision than other objective functions. The number of models in each subpopulation has a great influence on the rate of convergence and computation time, suggesting that it should be obtained through tests in practical problems. The critical separation radius should be determined carefully because it directly affects the multi-extreme values in the inversion. We also compare the inverted results from full-band waveform data and surface-wave frequency-band (0.02-0.1 Hz) data, and find that the latter is relatively poorer but still has a higher precision, suggesting that surface-wave frequency-band data can also be used to invert for the crustal structure.
Comparative study of texture detection and classification algorithms
NASA Astrophysics Data System (ADS)
Koltsov, P. P.
2011-08-01
A description and results of application of the computer system PETRA (performance evaluation of texture recognition algorithms) are given. This system is designed for the comparative study of texture analysis algorithms; it includes a database of textured images and a collection of software implementations of texture analysis algorithms. The functional capabilities of the system are illustrated using texture classification examples. Test examples are taken from the Brodatz album, MeasTech database, and a set of aerospace images. Results of a comparative evaluation of five well-known texture analysis methods are described—Gabor filters, Laws masks, ring/wedge filters, gray-level cooccurrence matrices (GLCMs), and autoregression image model.
Numerical study of dynamical mass generation in QED3
NASA Astrophysics Data System (ADS)
Bashir, A.; Huet, A.; Raya, A.
2006-05-01
We carry out a numerical study of dynamical generation of fermion masses by solving the Schwinger-Dyson equation for the fermion propagator in three-dimensional quenched Quantum Electrodynamics (QED3) in various gauges. We employ an ansatz for the three-point vertex which satisfies the Ward-Green-Takahashi identity, namely, the Ball-Chiu Vertex. We discuss the advantages of our numerical method over some earlier ones.
Real-space, mean-field algorithm to numerically calculate long-range interactions
NASA Astrophysics Data System (ADS)
Cadilhe, A.; Costa, B. V.
2016-02-01
Long-range interactions are known to be of difficult treatment in statistical mechanics models. There are some approaches that introduce a cutoff in the interactions or make use of reaction field approaches. However, those treatments suffer the illness of being of limited use, in particular close to phase transitions. The use of open boundary conditions allows the sum of the long-range interactions over the entire system to be done, however, this approach demands a sum over all degrees of freedom in the system, which makes a numerical treatment prohibitive. Techniques like the Ewald summation or fast multipole expansion account for the exact interactions but are still limited to a few thousands of particles. In this paper we introduce a novel mean-field approach to treat long-range interactions. The method is based in the division of the system in cells. In the inner cell, that contains the particle in sight, the 'local' interactions are computed exactly, the 'far' contributions are then computed as the average over the particles inside a given cell with the particle in sight for each of the remaining cells. Using this approach, the large and small cells limits are exact. At a fixed cell size, the method also becomes exact in the limit of large lattices. We have applied the procedure to the two-dimensional anisotropic dipolar Heisenberg model. A detailed comparison between our method, the exact calculation and the cutoff radius approximation were done. Our results show that the cutoff-cell approach outperforms any cutoff radius approach as it maintains the long-range memory present in these interactions, contrary to the cutoff radius approximation. Besides that, we calculated the critical temperature and the critical behavior of the specific heat of the anisotropic Heisenberg model using our method. The results are in excellent agreement with extensive Monte Carlo simulations using Ewald summation.
Study and application on accelerated algorithm of ray-casting
NASA Astrophysics Data System (ADS)
Sun, Xiaoping; Wu, Jian; Cui, Zhiming; Ma, Jianlin
2007-12-01
Medical image 3D reconstruct is an important application filed for volume rendering, for it special using, it required fast interactive speed and high image quality. The ray casting algorithm (RCA) is a widely used basic volume rendering algorithm. It can get high quality image but the rendering speed is very slowly for powerful computing capacity. Due to these shortcomings and deficiencies, the accelerated ray casting algorithm is presented in this paper to improve its rendering speed and apply it to medical image 3D reconstruct. Firstly, accelerate algorithms for ray casting are fully studied and compared. Secondly, improved tri-linear interpolation technology has been selected and extended to continuous ray casting in order to reduce matrix computation by matrix transformation characteristics of re-sampling points. Then ray interval casting technology is used to reduce the number of rays. Utilizing volume data sets cropping technology that improving boundary box technique avoids the sampling in empty voxel. Finally, the synthesized accelerate algorithm has been proposed. The result shown that compare with standard ray casting algorithm, the accelerate algorithm not only improve the rendering speed but also produce the required quality images.
Experimental and numerical study on condensation in transonic steam flow
NASA Astrophysics Data System (ADS)
Majkut, Mirosław; Dykas, Sławomir; Strozik, Michał; Smołka, Krystian
2015-09-01
The present paper describes an experimental and numerical study of steam condensing flow in a linear cascade of turbine stator blades. The experimental research was performed on the facility of a small scale steam power plant located at Silesian University of Technology in Gliwice, Poland. The test rig of the facility allows us to perform the tests of steam transonic flows for the conditions corresponding to these which prevail in the low-pressure (LP) condensing steam turbine stages. The experimental data of steam condensing flow through the blade-to- blade stator channel were compared with numerical results obtained using the in-house CFD numerical code TraCoFlow. Obtained results confirmed a good quality of the performed experiment and numerical calculations.
Study on the numerical schemes for hypersonic flow simulation
NASA Astrophysics Data System (ADS)
Nagdewe, S. P.; Shevare, G. R.; Kim, Heuy-Dong
2009-10-01
Hypersonic flow is full of complex physical and chemical processes, hence its investigation needs careful analysis of existing schemes and choosing a suitable scheme or designing a brand new scheme. The present study deals with two numerical schemes Harten, Lax, and van Leer with Contact (HLLC) and advection upstream splitting method (AUSM) to effectively simulate hypersonic flow fields, and accurately predict shock waves with minimal diffusion. In present computations, hypersonic flows have been modeled as a system of hyperbolic equations with one additional equation for non-equilibrium energy and relaxing source terms. Real gas effects, which appear typically in hypersonic flows, have been simulated through energy relaxation method. HLLC and AUSM methods are modified to incorporate the conservation laws for non-equilibrium energy. Numerical implementation have shown that non-equilibrium energy convect with mass, and hence has no bearing on the basic numerical scheme. The numerical simulation carried out shows good comparison with experimental data available in literature. Both numerical schemes have shown identical results at equilibrium. Present study has demonstrated that real gas effects in hypersonic flows can be modeled through energy relaxation method along with either AUSM or HLLC numerical scheme.
NASA Astrophysics Data System (ADS)
Tang, Yu-Hang; Karniadakis, George Em
2014-11-01
We present a scalable dissipative particle dynamics simulation code, fully implemented on the Graphics Processing Units (GPUs) using a hybrid CUDA/MPI programming model, which achieves 10-30 times speedup on a single GPU over 16 CPU cores and almost linear weak scaling across a thousand nodes. A unified framework is developed within which the efficient generation of the neighbor list and maintaining particle data locality are addressed. Our algorithm generates strictly ordered neighbor lists in parallel, while the construction is deterministic and makes no use of atomic operations or sorting. Such neighbor list leads to optimal data loading efficiency when combined with a two-level particle reordering scheme. A faster in situ generation scheme for Gaussian random numbers is proposed using precomputed binary signatures. We designed custom transcendental functions that are fast and accurate for evaluating the pairwise interaction. The correctness and accuracy of the code is verified through a set of test cases simulating Poiseuille flow and spontaneous vesicle formation. Computer benchmarks demonstrate the speedup of our implementation over the CPU implementation as well as strong and weak scalability. A large-scale simulation of spontaneous vesicle formation consisting of 128 million particles was conducted to further illustrate the practicality of our code in real-world applications. Catalogue identifier: AETN_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AETN_v1_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 3 No. of lines in distributed program, including test data, etc.: 1 602 716 No. of bytes in distributed program, including test data, etc.: 26 489 166 Distribution format: tar.gz Programming language: C/C++, CUDA C/C++, MPI. Computer: Any computers having nVidia GPGPUs with compute capability 3.0. Operating system: Linux. Has the code been
A computational study of routing algorithms for realistic transportation networks
Jacob, R.; Marathe, M.V.; Nagel, K.
1998-12-01
The authors carry out an experimental analysis of a number of shortest path (routing) algorithms investigated in the context of the TRANSIMS (Transportation Analysis and Simulation System) project. The main focus of the paper is to study how various heuristic and exact solutions, associated data structures affected the computational performance of the software developed especially for realistic transportation networks. For this purpose the authors have used Dallas Fort-Worth road network with very high degree of resolution. The following general results are obtained: (1) they discuss and experimentally analyze various one-one shortest path algorithms, which include classical exact algorithms studied in the literature as well as heuristic solutions that are designed to take into account the geometric structure of the input instances; (2) they describe a number of extensions to the basic shortest path algorithm. These extensions were primarily motivated by practical problems arising in TRANSIMS and ITS (Intelligent Transportation Systems) related technologies. Extensions discussed include--(i) time dependent networks, (ii) multi-modal networks, (iii) networks with public transportation and associated schedules. Computational results are provided to empirically compare the efficiency of various algorithms. The studies indicate that a modified Dijkstra`s algorithm is computationally fast and an excellent candidate for use in various transportation planning applications as well as ITS related technologies.
Numerical study of subcritical flow with fluid injection
NASA Technical Reports Server (NTRS)
Balasubramanian, R.
1990-01-01
It is suggested that the study of synthetic flows, where controlled experiments can be performed, is useful in understanding turbulent flow structures. The early states of formation of hairpin structures in shear flows and the subsequent evolution of these structures is studied in shear flows and the subsequent evolution of these structures is studied through numerical simulations, by developing full-time dependent three-dimensional flow solution of an initially laminar (subcritical) flow in which injection of fluid through a narrow streamwise slot from the bottom wall of a plate is carried out. Details of the numerical approach and significance of the present findings are reported in this work.
Experimental and Numerical Study of Free-Field Blast Mitigation
NASA Astrophysics Data System (ADS)
Allen, R. M.; Kirkpatrick, D. J.; Longbottom, A. W.; Milne, A. M.; Bourne, N. K.
2004-07-01
The development of a fundamental understanding of the mechanisms governing the attenuation of explosives effects by a surrounding mitigant material or system would benefit many civilian and military applications. Current approaches rely almost exclusively on empirical data, few if any truly predictive models exist. Dstl has recently pursued an experimental programme investigating the mitigation of effects from detonating explosives in support of general requirements to attenuate blast and fragmentation. The physical properties of a range of mitigant materials have been studied at a more fundamental level, both experimentally and numerically. A preliminary numerical parameter study has been undertaken by FGE using two-phase numerical simulations to complement the experimental studies. Initial work used idealised equations of state for generic mitigants but more recently material characterisation experiments have been undertaken at RMCS. Results confirm that porosity and particle density are dominant factors affecting the efficiency of the mitigant in reducing free-field blast.
A Numerical Study of Hypersonic Forebody/Inlet Integration Problem
NASA Technical Reports Server (NTRS)
Kumar, Ajay
1991-01-01
A numerical study of hypersonic forebody/inlet integration problem is presented in the form of the view-graphs. The following topics are covered: physical/chemical modeling; solution procedure; flow conditions; mass flow rate at inlet face; heating and skin friction loads; 3-D forebogy/inlet integration model; and sensitivity studies.
Numerical studies of laminar and turbulent drag reduction
NASA Technical Reports Server (NTRS)
Balasubramanian, R.; Orszag, S. A.
1981-01-01
Two-dimensional incompressible flow over wavy surfaces is studied numerically by spectral methods. Turbulence effects are modeled. Results for symmetric and asymmetric wave forms are presented. Effect of propagating surface waves on drag reduction is studied. Comparisons between computer simulations and experimental results are made.
Numerical Study of Magnetic Damping During Unidirectional Solidification
NASA Technical Reports Server (NTRS)
Li, Ben Q.
1997-01-01
A fully 3-D numerical model is developed to represent magnetic damping of complex fluid flow, heat transfer and electromagnetic field distributions in a melt cavity. The model is developed based on our in-house finite element code for the fluid flow, heat transfer and electromagnetic field calculations. The computer code has been tested against benchmark test problems that are solved by other commercial codes as well as analytical solutions whenever available. The numerical model is tested against numerical and experimental results for water reported in literature. With the model so tested, various numerical simulations are carried out for the Sn-35.5% Pb melt convection and temperature distribution in a cylindrical cavity with and without the presence of a transverse magnetic field. Numerical results show that magnetic damping can be effectively applied to reduce turbulence and flow levels in the melt undergoing solidification and over a certain threshold value a higher magnetic field resulted in a higher velocity reduction. It is found also that for a fully 3-D representation of the magnetic damping effects, the electric field induced in the melt by the applied DC magnetic field does not vanish, as some researchers suggested, and must be included even for molten metal and semiconductors. Also, for the study of the melt flow instability, a long enough time has to be applied to ensure the final fluid flow recirculation pattern. Moreover, our numerical results suggested that there seems to exist a threshold value of applied magnetic field, above which magnetic damping becomes possible and below which the convection in the melt is actually enhanced. Because of the limited financial resource allocated for the project, we are unable to carry out extensive study on this effect, which should warrant further theoretical and experimental study. In that endeavor, the developed numerical model should be very useful; and the model should serve as a useful tool for exploring
Numerical aerodynamic simulation facility preliminary study: Executive study
NASA Technical Reports Server (NTRS)
1977-01-01
A computing system was designed with the capability of providing an effective throughput of one billion floating point operations per second for three dimensional Navier-Stokes codes. The methodology used in defining the baseline design, and the major elements of the numerical aerodynamic simulation facility are described.
Study on gas kinetic unified algorithm for flows from rarefied transition to continuum
NASA Astrophysics Data System (ADS)
Li, Zhi-Hui; Zhang, Han-Xin
2004-01-01
The modified BGK equation adapted to various flow regimes can be presented by the aid of the basic characteristics on molecular movement and collision approaching to equilibrium. The discrete velocity ordinate method is developed and applied to the velocity distribution function to remove its continuous dependency on the velocity space, and then the velocity distribution function equation is cast into hyperbolic conservation law form with nonlinear source terms. Based on the unsteady time-splitting method and the non-oscillatory, containing no free parameters, and dissipative (NND) scheme, the gas kinetic finite difference second-order scheme is constructed for the computation of the discrete velocity distribution functions. The mathematical model on the interaction of molecules with solid surface is studied and used in the numerical method. Four types of numerical quadrature rules, such as the modified Gauss-Hermite formula, the composite Newton-Cotes integration method, the Gauss-Legendre numerical quadrature rule, and the Golden Section number-theoretic integral method, are developed and applied to the discretized velocity space to evaluate the macroscopic flow parameters at each point in the physical space. As a result, a unified simplified gas kinetic algorithm is established for the flows from rarefied transition to continuum regime. Based on analyzing the inner parallel degree of the unified algorithm, the parallel strategy adapted to the gas kinetic numerical algorithm is studied, and then the HPF parallel processing software for the unified algorithm is developed. To test the present method, the one-dimensional shock-tube problems, the flows past two-dimensional circular cylinder, and the flows around three-dimensional sphere and spacecraft shape with various Knudsen numbers are simulated. The computational results are found in high resolution of the flow fields and good agreement with the theoretical, DSMC, N-S, and experimental results.
Numerical Algorithm Based on Haar-Sinc Collocation Method for Solving the Hyperbolic PDEs
Javadi, H. H. S.; Navidi, H. R.
2014-01-01
The present study investigates the Haar-Sinc collocation method for the solution of the hyperbolic partial telegraph equations. The advantages of this technique are that not only is the convergence rate of Sinc approximation exponential but the computational speed also is high due to the use of the Haar operational matrices. This technique is used to convert the problem to the solution of linear algebraic equations via expanding the required approximation based on the elements of Sinc functions in space and Haar functions in time with unknown coefficients. To analyze the efficiency, precision, and performance of the proposed method, we presented four examples through which our claim was confirmed. PMID:25485295
Numerical study of damage growth in particulate composites
Kwon, Y.W.; Liu, C.T.
1999-10-01
A numerical study was conducted to simulate and predict damage initiation and growth around the crack tip in particulate composite specimens made of hard particles embedded in a soft rubber-like matrix material. Therefore, damage evolution in the matrix material around crack tips was investigated. The progressive damage was modeled using a micro/macro-approach which combined two levels of analyses like the micro-level and the macro-level analyses. Damage description was undertaken at the microlevel using a simplified three-dimensional unit-cell model and an isotropic continuum damage theory. The numerical study examined both thin and thick specimens with a short or long edge crack to understand the effects of specimen thickness and crack size on the damage initiation, growth, and saturation. Numerical results were compared with experimental data.
Numerical study of axisymmetric collapses of submarine granular > columns
NASA Astrophysics Data System (ADS)
Monsorno, Davide; Varsakelis, Christos
2014-11-01
In this talk, we report on the results of a numerical study of the axisymmetric collapse of subaqueous granular columns. Our study is based on a 2-pressure, 2-velocity continuum flow model for fluid-saturated granular materials. This model is integrated via a multi-phase projection method that incorporates a regularization method for the treatment of material interfaces. In our simulations, a dense column of a granular material immersed in water is placed on a horizontal plane and is allowed to collapse and spread due to its weight. Emphasis is placed on the run-out distance and the termination height and their correlation with the aspect ratio, the volume fraction and the diameter of the grains. Comparisons against experimental measurements and previous numerical predictions are also performed. Finally, in order to examine and quantify the role of the interstitial fluid, we compare our numerical predictions against experimental results from column collapses of dry granular materials.
A numerical method to study the dynamics of capillary fluid systems
NASA Astrophysics Data System (ADS)
Herrada, M. A.; Montanero, J. M.
2016-02-01
We propose a numerical approach to study both the nonlinear dynamics and linear stability of capillary fluid systems. In the nonlinear analysis, the time-dependent fluid region is mapped onto a fixed numerical domain through a coordinate transformation. The hydrodynamic equations are spatially discretized with the Chebyshev spectral collocation technique, while an implicit time advancement is performed using second-order backward finite differences. The resulting algebraic equations are solved with the iterative Newton-Raphson technique. The most novel aspect of the method is the fact that the elements of the Jacobian of the discretized system of equations are symbolic functions calculated before running the simulation. These functions are evaluated numerically in the Newton-Raphson iterations to find the solution at each time step, which reduces considerably the computing time. Besides, this numerical procedure can be easily adapted to solve the eigenvalue problem which determines the linear global modes of the capillary system. Therefore, both the nonlinear dynamics and the linear stability analysis can be conducted with essentially the same algorithm. We validate this numerical approach by studying the dynamics of a liquid bridge close to its minimum volume stability limit. The results are virtually the same as those obtained with other methods. The proposed approach proves to be much more computationally efficient than those other methods. Finally, we show the versatility of the method by calculating the linear global modes of a gravitational jet.
A Numerical/Experimental Study of Nitinol Actuator Springs
NASA Astrophysics Data System (ADS)
Auricchio, Ferdinando; Scalet, Giulia; Urbano, Marco
2014-07-01
This study deals with the numerical modeling, simulation and experimental analysis of shape-memory alloy (SMA) helicoidal springs. An experimental campaign is conducted on both SMA straight wires and helicoidal springs that experienced the same annealing process. Then, we use such experimental results to investigate three phenomenological constitutive models able to represent SMA macroscopic behavior. In particular, after the identification of all the material parameters from experimental results on SMA wires, we inspect the thermo-mechanical behavior of SMA helicoidal springs by comparing numerical predictions to experimental data. Finally, we discuss models capabilities and some aspects characterizing SMA material behavior.
Numerical study of the small scale structures in Boussinesq convection
NASA Technical Reports Server (NTRS)
Weinan, E.; Shu, Chi-Wang
1992-01-01
Two-dimensional Boussinesq convection is studied numerically using two different methods: a filtered pseudospectral method and a high order accurate Essentially Nonoscillatory (ENO) scheme. The issue whether finite time singularity occurs for initially smooth flows is investigated. The numerical results suggest that the collapse of the bubble cap is unlikely to occur in resolved calculations. The strain rate corresponding to the intensification of the density gradient across the front saturates at the bubble cap. We also found that the cascade of energy to small scales is dominated by the formulation of thin and sharp fronts across which density jumps.
Experimental and numerical study of pulsating transversal jets
NASA Astrophysics Data System (ADS)
Goldfeld, M. A.; Fedorova, N. N.; Fedorchenko, I. A.; Pozdnyakov, G. A.; Timofeev, K. Yu.; Zhakharova, Yu. V.
2015-06-01
Paper presents results of joint experimental and numerical investigation of pulsating jet penetration into still air and supersonic flow. Goal of the study is to investigate two-dimensional (2D) Hartmann generator (HG) properties and clear up its possibilities in providing better mixing between air and secondary (injected) gases.
Numerical study of the Kerr solution in rotating coordinates
NASA Astrophysics Data System (ADS)
Bai, S.; Izquierdo, G.; Klein, C.
2016-06-01
The Kerr solution in coordinates corotating with the horizon is studied as a testbed for a spacetime with a helical Killing vector in the Ernst picture. The solution is numerically constructed by solving the Ernst equation with a spectral method and a Newton iteration. We discuss convergence of the iteration for several initial iterates and different values of the Kerr parameters.
Numerical studies of gravity effects in two-phase reservoirs
Bodvarsson, G.S.; Cox, B.L.
1986-06-01
Numerical studies are performed to investigate the effects of localized feed zones on the pressure transients in two-phase reservoirs. It is shown that gravity effects can significantly affect the pressure transients, because of the large difference in the density of liquid water and vapor. Pressure transients for shallow and deep feed zones and the resulting fluid flow patterns are discussed.
Numerical study of porosity in titanium dental castings.
Wu, M; Sahm, P R; Augthun, M; Spiekermann, H; Schädlich-Stubenrauch, J
1999-09-01
A commercial software package, MAGMASOFT (MAGMA Giessereitechnologie GmbH, Aachen, Germany), was used to study shrinkage and gas porosity in titanium dental castings. A geometrical model for two simplified tooth crowns connected by a connector bar was created. Both mold filling and solidification of this casting model were numerically simulated. Shrinkage porosity was quantitatively predicted by means of a built-in feeding criterion. The risk of gas pore formation was investigated using the numerical filling and solidification results. The results of the numerical simulations were compared with experiments, which were carried out on a centrifugal casting machine with an investment block mold. The block mold was made of SiO2 based slurry with a 1 mm thick Zr2 face coat to reduce metal-mold reactions. Both melting and casting were carried out under protective argon (40 kPa). The finished castings were sectioned and the shrinkage porosity determined. The experimentally determined shrinkage porosity coincided with the predicted numerical simulation results. No apparent gas porosity was found in these model castings. Several running and gating systems for the above model casting were numerically simulated. An optimized running and gating system design was then experimentally cast, which resulted in porosity-free castings. PMID:15348102
Study of genetic direct search algorithms for function optimization
NASA Technical Reports Server (NTRS)
Zeigler, B. P.
1974-01-01
The results are presented of a study to determine the performance of genetic direct search algorithms in solving function optimization problems arising in the optimal and adaptive control areas. The findings indicate that: (1) genetic algorithms can outperform standard algorithms in multimodal and/or noisy optimization situations, but suffer from lack of gradient exploitation facilities when gradient information can be utilized to guide the search. (2) For large populations, or low dimensional function spaces, mutation is a sufficient operator. However for small populations or high dimensional functions, crossover applied in about equal frequency with mutation is an optimum combination. (3) Complexity, in terms of storage space and running time, is significantly increased when population size is increased or the inversion operator, or the second level adaptation routine is added to the basic structure.
Numerical Model Studies of the Martian Mesoscale Circulations
NASA Technical Reports Server (NTRS)
Segal, Moti; Arritt, Raymond W.
1997-01-01
The study objectives were to evaluate by numerical modeling various possible mesoscale circulation on Mars and related atmospheric boundary layer processes. The study was in collaboration with J. Tillman of the University of Washington (who supported the study observationally). Interaction has been made with J. Prusa of Iowa State University in numerical modeling investigation of dynamical effects of topographically-influenced flow. Modeling simulations included evaluations of surface physical characteristics on: (i) the Martian atmospheric boundary layer and (ii) their impact on thermally and dynamically forced mesoscale flows. Special model evaluations were made in support of selection of the Pathfinder landing sites. J. Tillman's finding of VL-2 inter-annual temperature difference was followed by model simulations attempting to point out the forcing for this feature. Publication of the results in the reviewed literature in pending upon completion of the manuscripts in preparation as indicated later.
Numerical studies of HF Doppler variations caused by ionospheric disturbances
NASA Astrophysics Data System (ADS)
Takefu, M.; Hiroshige, N.
HF Doppler variations caused by ionospheric disturbances are studied using an ionosphere model containing sinusoidal traveling electron density fluctuations. The present study uses a more realistic ionosphere model and a more accurate numerical method than previous works using corrugated specular reflector models. The study gives a clue to estimate the TID-associated fluctuations of ionospheric electron density by means of HF Doppler measurements. It is shown that some kinds of characteristic HF Doppler traces result depending on the wavelength of the disturbance and its traveling direction. Numerical results suggest that more or less 5 percent of the background electron density can explain most of the quasi-periodic variations on the observed HF Doppler records.
Numerical studies of the nonlinear properties of composites
NASA Astrophysics Data System (ADS)
Zhang, X.; Stroud, D.
1994-01-01
Using both numerical and analytical techniques, we investigate various ways to enhance the cubic nonlinear susceptibility χe of a composite material. We start from the exact relation χe =tsumipiχi<(E.E)2>i,lin/ E40, where χi and pi are the cubic nonlinear susceptibility and volume fraction of the ith component, E0 is the applied electric field, and
Analytical and Numerical Studies of Coherent Field Configurations
NASA Astrophysics Data System (ADS)
Muller, Hans-Reinhard
Motivated by the goal of investigating the dynamics of the electroweak phase transition in the early Universe, a study of dynamical aspects of nonlinear field-theoretical systems is performed. Symmetric and asymmetric double-well potentials in the φ4 theory are used as a model for these systems. In the first part, it is shown that in this model, a class of localized, time-dependent, spherically-symmetric objects dubbed oscillons exists. The most distinctive feature of these objects is that they are extremely long-lived. The properties of oscillons are studied by numerical and analytical means. In the second part of the study, the matching between (1+1)-dimensional nonlinear field theories coupled to an external stochastic environment and their lattice simulations is investigated. In particular, a method is developed to obtain numerical results which are lattice-spacing independent, and to extract the correct effective potential which emerges from the simulations. As an application, the thermal production of kinks is studied, obtaining the lattice-spacing independent number density of kinks and the effective barrier for kink production, i.e., the effective kink mass. Within its range of validity, the approach can be used to match numerical simulations to continuum studies of the emergence of coherent field structures in cosmology and condensed matter physics.
3D-radiation hydro simulations of disk-planet interactions. I. Numerical algorithm and test cases
NASA Astrophysics Data System (ADS)
Klahr, H.; Kley, W.
2006-01-01
We study the evolution of an embedded protoplanet in a circumstellar disk using the 3D-Radiation Hydro code TRAMP, and treat the thermodynamics of the gas properly in three dimensions. The primary interest of this work lies in the demonstration and testing of the numerical method. We show how far numerical parameters can influence the simulations of gap opening. We study a standard reference model under various numerical approximations. Then we compare the commonly used locally isothermal approximation to the radiation hydro simulation using an equation for the internal energy. Models with different treatments of the mass accretion process are compared. Often mass accumulates in the Roche lobe of the planet creating a hydrostatic atmosphere around the planet. The gravitational torques induced by the spiral pattern of the disk onto the planet are not strongly affected in the average magnitude, but the short time scale fluctuations are stronger in the radiation hydro models. An interesting result of this work lies in the analysis of the temperature structure around the planet. The most striking effect of treating the thermodynamics properly is the formation of a hot pressure-supported bubble around the planet with a pressure scale height of H/R ≈ 0.5 rather than a thin Keplerian circumplanetary accretion disk.
Numerical study on 3D composite morphing actuators
NASA Astrophysics Data System (ADS)
Oishi, Kazuma; Saito, Makoto; Anandan, Nishita; Kadooka, Kevin; Taya, Minoru
2015-04-01
There are a number of actuators using the deformation of electroactive polymer (EAP), where fewer papers seem to have focused on the performance of 3D morphing actuators based on the analytical approach, due mainly to their complexity. The present paper introduces a numerical analysis approach on the large scale deformation and motion of a 3D half dome shaped actuator composed of thin soft membrane (passive material) and EAP strip actuators (EAP active coupon with electrodes on both surfaces), where the locations of the active EAP strips is a key parameter. Simulia/Abaqus Static and Implicit analysis code, whose main feature is the high precision contact analysis capability among structures, are used focusing on the whole process of the membrane to touch and wrap around the object. The unidirectional properties of the EAP coupon actuator are used as input data set for the material properties for the simulation and the verification of our numerical model, where the verification is made as compared to the existing 2D solution. The numerical results can demonstrate the whole deformation process of the membrane to wrap around not only smooth shaped objects like a sphere or an egg, but also irregularly shaped objects. A parametric study reveals the proper placement of the EAP coupon actuators, with the modification of the dome shape to induce the relevant large scale deformation. The numerical simulation for the 3D soft actuators shown in this paper could be applied to a wider range of soft 3D morphing actuators.
Numerical Study on Cryogenic Coflowing Jets under Transcritical Conditions
NASA Astrophysics Data System (ADS)
Tani, Hiroumi; Teramoto, Susumu; Okamoto, Koji; Yamanishi, Nobuhiro
2012-11-01
A numerical and experimental study is presented on cryogenic coflowing jets under transcritical conditions for a better understanding of the propellant mixing in supercritical-pressure rocket engines. The major concerns are dominant flow structures in the mixing of cryogenic coflowing jets under transcritical conditions. Experimentally, in advance of detailed numerical simulations, cryogenic nitrogen/gaseous nitrogen coaxial jets were visualized by the backlighting photography technique. It was observed that a dense nitrogen core has a shear-layer instability near the injector exit and eventually breaks up into large lumps which dissolve and fade away downstream. In numerical simulations, LES technique was employed for more detailed discussion on the flow structures. LES of a cryogenic nitrogen/gaseous nitrogen coflowing plane jet was conducted with the same density and velocity ratios of inner/outer jets as the experiments. As observed in the experiments, the shear-layer instability in the inner mixing layers is predominant near the injector exit. After roll-up and paring, the shear-layer instability waves become large-scale vortices. They cause coherent vortex structures which become dominant in the downstream and break the dense core into lumps. Strouhal numbers of the shear-layer instability and the dense lump shedding in the numerical simulations were comparable to those measured in the experiments, respectively.
NASA Astrophysics Data System (ADS)
Saito, Kyosuke; Tanabe, Tadao; Oyama, Yutaka
2016-04-01
We have presented a numerical analysis to describe the behavior of a second harmonic generation (SHG) in THz regime by taking into account for both linear and nonlinear optical susceptibility. We employed a nonlinear finite-difference-time-domain (nonlinear FDTD) method to simulate SHG output characteristics in THz photonic crystal waveguide based on semi insulating gallium phosphide crystal. Unique phase matching conditions originated from photonic band dispersions with low group velocity are appeared, resulting in SHG output characteristics. This numerical study provides spectral information of SHG output in THz PC waveguide. THz PC waveguides is one of the active nonlinear optical devices in THz regime, and nonlinear FDTD method is a powerful tool to design photonic nonlinear THz devices.
A Numerical Study of the Effects of Curvature and Convergence on Dilution Jet Mixing
NASA Technical Reports Server (NTRS)
Holdeman, J. D.; Reynolds, R.; White, C.
1987-01-01
An analytical program was conducted to assemble and assess a three-dimensional turbulent viscous flow computer code capable of analyzing the flow field in the transition liners of small gas turbine engines. This code is of the TEACH type with hybrid numerics, and uses the power law and SIMPLER algorithms, an orthogonal curvilinear coordinate system, and an algebraic Reynolds stress turbulence model. The assessments performed in this study, consistent with results in the literature, showed that in its present form this code is capable of predicting trends and qualitative results. The assembled code was used to perform a numerical experiment to investigate the effects of curvature and convergence in the transition liner on the mixing of single and opposed rows of cool dilution jets injected into a hot mainstream flow.
A numerical study of the effects of curvature and convergence on dilution jet mixing
NASA Technical Reports Server (NTRS)
Holdeman, J. D.; Reynolds, R.; White, C.
1987-01-01
An analytical program was conducted to assemble and assess a three-dimensional turbulent viscous flow computer code capable of analyzing the flow field in the transition liners of small gas turbine engines. This code is of the TEACH type with hybrid numerics, and uses the power law and SIMPLER algorithms, an orthogonal curvilinear coordinate system, and an algebraic Reynolds stress turbulence model. The assessments performed in this study, consistent with results in the literature, showed that in its present form this code is capable of predicting trends and qualitative results. The assembled code was used to perform a numerical experiment to investigate the effects of curvature and convergence in the transition liner on the mixing of single and opposed rows of cool dilution jets injected into a hot mainstream flow.
On the Numerical Study of Heavy Rainfall in Taiwan
NASA Technical Reports Server (NTRS)
Tao, Wei-Kuo; Chen, Ching-Sen; Chen, Yi-Leng; Jou, Ben Jong-Dao; Lin, Pay-Liam; Starr, David OC. (Technical Monitor)
2001-01-01
Heavy rainfall events are frequently observed over the western side of the CMR (central mountain range), which runs through Taiwan in a north-south orientation, in a southwesterly flow regime and over the northeastern side of the CMR in a northeasterly flow regime. Previous studies have revealed the mechanisms by which the heavy rainfall events are formed. Some of them have examined characteristics of the heavy rainfall via numerical simulations. In this paper, some of the previous numerical studies on heavy rainfall events around Taiwan during the Mei-Yu season (May and June), summer (non-typhoon cases) and autumn will be reviewed. Associated mechanisms proposed from observational studies will be reviewed first, and then characteristics of numerically simulated heavy rainfall events will be presented. The formation mechanisms of heavy rainfall from simulated results and from observational analysis are then compared and discussed. Based on these previous modeling studies, we will also discuss what are the major observations and modeling processes which will be needed for understanding the heavy precipitation in the future.
Numerical Studies of a Fluidic Diverter for Flow Control
NASA Technical Reports Server (NTRS)
Gokoglu, Suleyman A.; Kuczmarski, Maria A.; Culley, Dennis E.; Raghu, Surya
2009-01-01
The internal flow structure in a specific fluidic diverter is studied over a range from low subsonic to sonic inlet conditions by a time-dependent numerical analysis. The understanding will aid in the development of fluidic diverters with minimum pressure losses and advanced designs of flow control actuators. The velocity, temperature and pressure fields are calculated for subsonic conditions and the self-induced oscillatory behavior of the flow is successfully predicted. The results of our numerical studies have excellent agreement with our experimental measurements of oscillation frequencies. The acoustic speed in the gaseous medium is determined to be a key factor for up to sonic conditions in governing the mechanism of initiating the oscillations as well as determining its frequency. The feasibility of employing plasma actuation with a minimal perturbation level is demonstrated in steady-state calculations to also produce oscillation frequencies of our own choosing instead of being dependent on the fixed-geometry fluidic device.
Numerical study of fractional nonlinear Schrödinger equations
Klein, Christian; Sparber, Christof; Markowich, Peter
2014-01-01
Using a Fourier spectral method, we provide a detailed numerical investigation of dispersive Schrödinger-type equations involving a fractional Laplacian in an one-dimensional case. By an appropriate choice of the dispersive exponent, both mass and energy sub- and supercritical regimes can be identified. This allows us to study the possibility of finite time blow-up versus global existence, the nature of the blow-up, the stability and instability of nonlinear ground states and the long-time dynamics of solutions. The latter is also studied in a semiclassical setting. Moreover, we numerically construct ground state solutions of the fractional nonlinear Schrödinger equation. PMID:25484604
Numerical study of unsteady processes in a Faraday MHD generator
NASA Astrophysics Data System (ADS)
Vinogradova, G. N.; Panchenko, V. P.
1981-07-01
A numerical study is presented on the unsteady processes occurring in a Faraday MHD generator with a high power-conversion efficiency. A supersonic MHD generator operating with an equilibrium plasma and designed to convert energy in a system using a thermonuclear reactor is considered, and the steady operating modes are established for cases when an ohmic load is connected, disconnected, or reduced. A magnetic field is assumed to be generated by a suitable profiling of the external magnetic field, and the working medium is modeled by an ideal gas. Partial differential equations are solved numerically by using a central difference predictor-corrector scheme. The study can be applied to problems (e.g., transient times, nominal parameter maximal values and rates of change, methods of regulating the generator and switching it on and off) arising during the design of MHD generators.