Science.gov

Sample records for implicit integration scheme

  1. Performance and optimization of direct implicit time integration schemes for use in electrostatic particle simulation codes

    SciTech Connect

    Procassini, R.J.; Birdsall, C.K.; Morse, E.C.; Cohen, B.I.

    1988-01-01

    Implicit time integration schemes allow for the use of larger time steps than conventional explicit methods, thereby extending the applicability of kinetic particle simulation methods. This paper will describe a study of the performance and optimization of two such direct implicit schemes, which are used to follow the trajectories of charged particles in an electrostatic, particle-in-cell plasma simulation code. The direct implicit method that was used for this study is an alternative to the moment-equation implicit method. 10 refs., 7 figs., 4 tabs.

  2. Implicit and semi-implicit schemes: Algorithms

    NASA Astrophysics Data System (ADS)

    Keppens, R.; Tóth, G.; Botchev, M. A.; van der Ploeg, A.

    1999-06-01

    This study formulates general guidelines to extend an explicit code with a great variety of implicit and semi-implicit time integration schemes. The discussion is based on their specific implementation in the Versatile Advection Code, which is a general purpose software package for solving systems of non-linear hyperbolic (and/or parabolic) partial differential equations, using standard high resolution shock capturing schemes. For all combinations of explicit high resolution schemes with implicit and semi-implicit treatments, it is shown how second-order spatial and temporal accuracy for the smooth part of the solutions can be maintained. Strategies to obtain steady state and time accurate solutions implicitly are discussed. The implicit and semi-implicit schemes require the solution of large linear systems containing the Jacobian matrix. The Jacobian matrix itself is calculated numerically to ensure the generality of this implementation. Three options are discussed in terms of applicability, storage requirements and computational efficiency. One option is the easily implemented matrix-free approach, but the Jacobian matrix can also be calculated by using a general grid masking algorithm, or by an efficient implementation for a specific Lax-Friedrich-type total variation diminishing (TVD) spatial discretization. The choice of the linear solver depends on the dimensionality of the problem. In one dimension, a direct block tridiagonal solver can be applied, while in more than one spatial dimension, a conjugate gradient (CG)-type iterative solver is used. For advection-dominated problems, preconditioning is needed to accelerate the convergence of the iterative schemes. The modified block incomplete LU-preconditioner is implemented, which performs very well. Examples from two-dimensional hydrodynamic and magnetohydrodynamic computations are given. They model transonic stellar outflow and recover the complex magnetohydrodynamic bow shock flow in the switch-on regime

  3. IMEX-a : an adaptive, fifth order implicit-explicit integration scheme.

    SciTech Connect

    Brake, Matthew Robert

    2013-05-01

    This report presents an efficient and accurate method for integrating a system of ordinary differential equations, particularly those arising from a spatial discretization of partially differential equations. The algorithm developed, termed the IMEX a algorithm, belongs to a class of algorithms known as implicit-explicit (IMEX) methods. The explicit step is based on a fifth order Runge-Kutta explicit step known as the Dormand-Prince algorithm, which adaptively modifies the time step by calculating the error relative to a fourth order estimation. The implicit step, which follows the explicit step, is based on a backward Euler method, a special case of the generalized trapezoidal method. Reasons for choosing both of these methods, along with the algorithm development are presented. In applications that have less stringent accuracy requirements, several other methods are available through the IMEX a toolbox, each of which simplify the fifth order Dormand-Prince explicit step: the third order Bogacki-Shampine method, the second order Midpoint method, and the first order Euler method. The performance of the algorithm is evaluated on to examples. First, a two pawl system with contact is modeled. Results predicted by the IMEX a algorithm are compared to those predicted by six widely used integration schemes. The IMEX a algorithm is demonstrated to be significantly faster (by up to an order of magnitude) and at least as accurate as all of the other methods considered. A second example, an acoustic standing wave, is presented in order to assess the accuracy of the IMEX a algorithm. Finally, sample code is given in order to demonstrate the implementation of the proposed algorithm.

  4. A diagonally inverted LU implicit multigrid scheme

    NASA Technical Reports Server (NTRS)

    Yokota, Jeffrey W.; Caughey, David A.; Chima, Rodrick V.

    1988-01-01

    A new Diagonally Inverted LU Implicit scheme is developed within the framework of the multigrid method for the 3-D unsteady Euler equations. The matrix systems that are to be inverted in the LU scheme are treated by local diagonalizing transformations that decouple them into systems of scalar equations. Unlike the Diagonalized ADI method, the time accuracy of the LU scheme is not reduced since the diagonalization procedure does not destroy time conservation. Even more importantly, this diagonalization significantly reduces the computational effort required to solve the LU approximation and therefore transforms it into a more efficient method of numerically solving the 3-D Euler equations.

  5. Implicit schemes and parallel computing in unstructured grid CFD

    NASA Technical Reports Server (NTRS)

    Venkatakrishnam, V.

    1995-01-01

    The development of implicit schemes for obtaining steady state solutions to the Euler and Navier-Stokes equations on unstructured grids is outlined. Applications are presented that compare the convergence characteristics of various implicit methods. Next, the development of explicit and implicit schemes to compute unsteady flows on unstructured grids is discussed. Next, the issues involved in parallelizing finite volume schemes on unstructured meshes in an MIMD (multiple instruction/multiple data stream) fashion are outlined. Techniques for partitioning unstructured grids among processors and for extracting parallelism in explicit and implicit solvers are discussed. Finally, some dynamic load balancing ideas, which are useful in adaptive transient computations, are presented.

  6. Implicit TVD schemes for hyperbolic conservation laws in curvilinear coordinates

    NASA Technical Reports Server (NTRS)

    Yee, H. C.; Harten, A.

    1985-01-01

    The Harten (1983, 1984) total variation-diminishing (TVD) schemes, constituting a one-parameter explicit and implicit, second-order-accurate family, have the property of not generating spurious oscillations when applied to one-dimensional, nonlinear scalar hyperbolic conservation laws and constant coefficient hyperbolic systems. These methods are presently extended to the multidimensional hyperbolic conservation laws in curvilinear coordinates. Means by which to linearize the implicit operator and solution strategies, in order to improve the computation efficiency of the implicit algorithm, are discussed. Numerical experiments with steady state airfoil calculations indicate that the proposed linearized implicit TVD schemes are accurate and robust.

  7. Parallelization of implicit finite difference schemes in computational fluid dynamics

    NASA Technical Reports Server (NTRS)

    Decker, Naomi H.; Naik, Vijay K.; Nicoules, Michel

    1990-01-01

    Implicit finite difference schemes are often the preferred numerical schemes in computational fluid dynamics, requiring less stringent stability bounds than the explicit schemes. Each iteration in an implicit scheme involves global data dependencies in the form of second and higher order recurrences. Efficient parallel implementations of such iterative methods are considerably more difficult and non-intuitive. The parallelization of the implicit schemes that are used for solving the Euler and the thin layer Navier-Stokes equations and that require inversions of large linear systems in the form of block tri-diagonal and/or block penta-diagonal matrices is discussed. Three-dimensional cases are emphasized and schemes that minimize the total execution time are presented. Partitioning and scheduling schemes for alleviating the effects of the global data dependencies are described. An analysis of the communication and the computation aspects of these methods is presented. The effect of the boundary conditions on the parallel schemes is also discussed.

  8. Global Asymptotic Behavior of Iterative Implicit Schemes

    NASA Technical Reports Server (NTRS)

    Yee, H. C.; Sweby, P. K.

    1994-01-01

    The global asymptotic nonlinear behavior of some standard iterative procedures in solving nonlinear systems of algebraic equations arising from four implicit linear multistep methods (LMMs) in discretizing three models of 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed using the theory of dynamical systems. The iterative procedures include simple iteration and full and modified Newton iterations. The results are compared with standard Runge-Kutta explicit methods, a noniterative implicit procedure, and the Newton method of solving the steady part of the ODEs. Studies showed that aside from exhibiting spurious asymptotes, all of the four implicit LMMs can change the type and stability of the steady states of the differential equations (DEs). They also exhibit a drastic distortion but less shrinkage of the basin of attraction of the true solution than standard nonLMM explicit methods. The simple iteration procedure exhibits behavior which is similar to standard nonLMM explicit methods except that spurious steady-state numerical solutions cannot occur. The numerical basins of attraction of the noniterative implicit procedure mimic more closely the basins of attraction of the DEs and are more efficient than the three iterative implicit procedures for the four implicit LMMs. Contrary to popular belief, the initial data using the Newton method of solving the steady part of the DEs may not have to be close to the exact steady state for convergence. These results can be used as an explanation for possible causes and cures of slow convergence and nonconvergence of steady-state numerical solutions when using an implicit LMM time-dependent approach in computational fluid dynamics.

  9. Implicit Space-Time Conservation Element and Solution Element Schemes

    NASA Technical Reports Server (NTRS)

    Chang, Sin-Chung; Himansu, Ananda; Wang, Xiao-Yen

    1999-01-01

    Artificial numerical dissipation is in important issue in large Reynolds number computations. In such computations, the artificial dissipation inherent in traditional numerical schemes can overwhelm the physical dissipation and yield inaccurate results on meshes of practical size. In the present work, the space-time conservation element and solution element method is used to construct new and accurate implicit numerical schemes such that artificial numerical dissipation will not overwhelm physical dissipation. Specifically, these schemes have the property that numerical dissipation vanishes when the physical viscosity goes to zero. These new schemes therefore accurately model the physical dissipation even when it is extremely small. The new schemes presented are two highly accurate implicit solvers for a convection-diffusion equation. The two schemes become identical in the pure convection case, and in the pure diffusion case. The implicit schemes are applicable over the whole Reynolds number range, from purely diffusive equations to convection-dominated equations with very small viscosity. The stability and consistency of the schemes are analysed, and some numerical results are presented. It is shown that, in the inviscid case, the new schemes become explicit and their amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, their principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme.

  10. A factored implicit scheme for numerical weather prediction

    NASA Technical Reports Server (NTRS)

    Augenbaum, J. M.; Cohn, S. E.; Isaacson, E.; Dee, D. P.; Marchesin, D.

    1985-01-01

    An implicit method is proposed to factor the nonlinear partial differential equations governing fast and slow modes of dynamic motion in numerical weather prediction schemes. The method permits separate factorization of the slow and fast modes of the implicit operator. A simple two-dimensional version of the system of three-dimensional equations governing atmospheric dynamics over shallow water was analyzed to assess the accuracy of the proposed method. It is shown that the method has a small error which is comparable to other discretization errors in the overall scheme.

  11. An implicit-explicit Eulerian Godunov scheme for compressible flow

    SciTech Connect

    Collins, J.P.; Colella, P.; Glaz, H.M.

    1995-02-01

    A hybrid implicit-explicit scheme is developed for Eulerian hydrodynamics. The hybridization is a continuous switch and operates on each characteristic field separately. The explicit scheme is a version of the second-order Geodunov scheme; the implicit method is only first-order accurate in time but leads to a block tridiagonal matrix inversion for efficiency and is unconditionally stable for the case of linear advection. The methodology is described for the cases of linear advection, for nonlinear scalar problems, and for gas dynamics. An important element of our work is the use of a modified Engquist-Osher flux function in place of the Godunov flux. Several numerical results are presented to demonstrate the properties of the method, especially stable numerical shocks at very high CFL numbers and second-order accurate steady states. 24 refs., 6 figs., 2 tabs.

  12. The Implicit and Explicit alpha-mu Schemes

    NASA Technical Reports Server (NTRS)

    Chang, Sin-Chung; Himansu, Ananda

    1997-01-01

    Artificial numerical dissipation is an important issue in large Reynolds number computations. In such computations, the artificial dissipation inherent in traditional numerical schemes can overwhelm the physical dissipation and yield inaccurate results on meshes of practical size. In the present work, the space-time conservation element and solution element method is used to construct new and accurate numerical schemes such that artificial numerical dissipation will not overwhelm physical dissipation. Specifically, these schemes have the property that numerical dissipation vanishes when the physical viscosity goes to zero. These new schemes therefore accurately model the physical dissipation even when it is extremely small. The method of space-time conservation element and solution element, currently under development, is a nontraditional numerical method for solving conservation laws. The method is developed on the basis of local and global flux conservation in a space-time domain, in which space and time are treated in a unified manner. Explicit solvers for model and fluid dynamic conservation laws have previously been investigated. In this paper, we introduce a new concept in the design of implicit schemes, and use it to construct two highly accurate solvers for a convection-diffusion equation. The two schemes become identical in the pure convection case, and in the pure diffusion case. The implicit schemes are applicable over the whole Reynolds number range, from purely diffusive equations to purely inviscid (convective) equations. The stability and consistency of the schemes are analyzed, and some numerical results are presented. It is shown that, in the inviscid case, the new schemes become explicit and their amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, their principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme. We also construct an explicit solver

  13. Calculation of Supersonic Combustion Using Implicit Schemes

    NASA Technical Reports Server (NTRS)

    Yoon, Seokkwan; Kwak, Dochan (Technical Monitor)

    2003-01-01

    One of the technology goals of NASA for advanced space transportation is to develop highly efficient propulsion systems to reduce the cost of payload for space missions. Developments of rockets for the second generation Reusable Launch Vehicle (RLV) in the past several years have been focused on low-cost versions of conventional engines. However, recent changes in the Integrated Space Transportation Program to build a crew transportation vehicle to extend the life of the Space Shuttle fleet might suggest that air-breathing rockets could reemerge as a possible propulsion system for the third generation RLV to replace the Space Shuttle after 2015. The weight of the oxygen tank exceeds thirty percent of the total weight of the Space Shuttle at launch while the payload is only one percent of the total weight. The air-breathing rocket propulsion systems, which consume oxygen in the air, offer clear advantages by making vehicles lighter and more efficient. Experience in the National Aerospace Plane Program in the late 1980s indicates that scramjet engines can achieve high specific impulse for low hypersonic vehicle speeds. Whether taking a form of Rocket Based Combined Cycle (RBCC) or Turbine Based Combined Cycle (TBCC), the scramjet is an essential mode of operation for air-breathing rockets. It is well known that fuel-air mixing and rapid combustion are of crucial importance for the success of scramjet engines since the spreading rate of the supersonic mixing layer decreases as the Mach number increases. A factored form of the Gauss-Seidel relaxation method has been widely used in hypersonic flow research since its first application to non-equilibrium flows. However, difficulties in stability and convergence have been encountered when there is strong interaction between fluid motion and chemical reaction, such as multiple fuel injection problems. The present paper reports the results from investigation of the effect of modifications to the original algorithm on the

  14. Implicit lower-upper/approximate-factorization schemes for incompressible flows

    SciTech Connect

    Briley, W.R.; Neerarambam, S.S.; Whitfield, D.L.

    1996-10-01

    A lower-upper/approximate-factorization (LU/AF) scheme is developed for the incompressible Euler or Navier-Stokes equations. The LU/AF scheme contains an iteration parameter that can be adjusted to improve iterative convergence rate. The LU/AF scheme is to be used in conjunction with linearized implicit approximations and artificial compressibility to compute steady solutions, and within sub-iterations to compute unsteady solutions. Formulations based on time linearization with and without sub-iteration and on Newton linearization are developed using spatial difference operators. The spatial approximation used includes upwind differencing based on Roe`s approximate Riemann solver and van Leer`s MUSCL scheme, with numerically computed implicit flux linearizations. Simple one-dimensional diffusion and advection/diffusion problems are first studied analytically to provide insight for development of the Navier-Stokes algorithm. The optimal values of both time step and LU/AF parameter are determined for a test problem consisting of two-dimensional flow past a NACA 0012 airfoil, with a highly stretched grid. The optimal parameter provides a consistent improvement in convergence rate for four test cases having different grids and Reynolds numbers and, also, for an inviscid case. The scheme can be easily extended to three dimensions and adapted for compressible flows. 24 refs., 11 figs., 2 tabs.

  15. Implicit Block ACK Scheme for IEEE 802.11 WLANs

    PubMed Central

    Sthapit, Pranesh; Pyun, Jae-Young

    2016-01-01

    The throughput of IEEE 802.11 standard is significantly bounded by the associated Medium Access Control (MAC) overhead. Because of the overhead, an upper limit exists for throughput, which is bounded, including situations where data rates are extremely high. Therefore, an overhead reduction is necessary to achieve higher throughput. The IEEE 802.11e amendment introduced the block ACK mechanism, to reduce the number of control messages in MAC. Although the block ACK scheme greatly reduces overhead, further improvements are possible. In this letter, we propose an implicit block ACK method that further reduces the overhead associated with IEEE 802.11e’s block ACK scheme. The mathematical analysis results are presented for both the original protocol and the proposed scheme. A performance improvement of greater than 10% was achieved with the proposed implementation.

  16. An implicit, conservative, zonal-boundary scheme for Euler equation calculations

    NASA Technical Reports Server (NTRS)

    Rai, M. M.

    1985-01-01

    A zonal, or patched, grid approach is one in which the flow region of interest is divided into subregions which are then discretized independently, using existing grid generators. The equations of motion are integrated in each subregion in conjunction with zonal boundary schemes which allow proper information transfer across interfaces that separate subregions. The zonal approach greatly simplifies the treatment of complex geometries and also the addition of grid points to selected regions of the flow. A conservative, zonal boundary condition that could be used with explicit schemes was extended so that it can be used with existing second order accurate implicit integration schemes such as the Beam-Warming and Osher schemes. In the test case considered, the implicit schemes increased the rate of convergence considerably (by a factor of about 30 over that of the explicit scheme). Results demonstrating the time accuracy of the zonal scheme and the feasibility of performing calculations on zones that move relative to each other are also presented.

  17. Implicit preconditioned WENO scheme for steady viscous flow computation

    NASA Astrophysics Data System (ADS)

    Huang, Juan-Chen; Lin, Herng; Yang, Jaw-Yen

    2009-02-01

    A class of lower-upper symmetric Gauss-Seidel implicit weighted essentially nonoscillatory (WENO) schemes is developed for solving the preconditioned Navier-Stokes equations of primitive variables with Spalart-Allmaras one-equation turbulence model. The numerical flux of the present preconditioned WENO schemes consists of a first-order part and high-order part. For first-order part, we adopt the preconditioned Roe scheme and for the high-order part, we employ preconditioned WENO methods. For comparison purpose, a preconditioned TVD scheme is also given and tested. A time-derivative preconditioning algorithm is devised and a discriminant is devised for adjusting the preconditioning parameters at low Mach numbers and turning off the preconditioning at intermediate or high Mach numbers. The computations are performed for the two-dimensional lid driven cavity flow, low subsonic viscous flow over S809 airfoil, three-dimensional low speed viscous flow over 6:1 prolate spheroid, transonic flow over ONERA-M6 wing and hypersonic flow over HB-2 model. The solutions of the present algorithms are in good agreement with the experimental data. The application of the preconditioned WENO schemes to viscous flows at all speeds not only enhances the accuracy and robustness of resolving shock and discontinuities for supersonic flows, but also improves the accuracy of low Mach number flow with complicated smooth solution structures.

  18. Stability and accuracy of a semi-implicit Godunov scheme for mass transport

    NASA Astrophysics Data System (ADS)

    Bradford, Scott F.

    2004-06-01

    Semi-implicit, Godunov-type models are adapted for solving the two-dimensional, time-dependent, mass transport equation on a geophysical scale. The method uses Van Leer's MUSCL reconstruction in conjunction with an explicit, predictor-corrector method to discretize and integrate the advection and lateral diffusion portions of the governing equation to second-order spatial and temporal accuracy. Three classical schemes are investigated for computing advection: Lax-Wendroff, Warming-Beam, and Fromm. The proposed method uses second order, centred finite differences to spatially discretize the diffusion terms. In order to improve model stability and efficiency, vertical diffusion is implicitly integrated with the Crank-Nicolson method and implicit treatment of vertical diffusion in the predictor is also examined. Semi-discrete and Von Neumann analyses are utilized to compare the stability as well as the amplitude and phase accuracy of the proposed method with other explicit and semi-implicit schemes. Some linear, two-dimensional examples are solved and predictions are compared with the analytical solutions. Computational effort is also examined to illustrate the improved efficiency of the proposed model.

  19. Improvements to the RELAP5-3D Nearly-Implicit Numerical Scheme

    SciTech Connect

    Richard A. Riemke; Walter L. Weaver; RIchard R. Schultz

    2005-05-01

    The RELAP5-3D computer program has been improved with regard to its nearly-implicit numerical scheme for twophase flow and single-phase flow. Changes were made to the nearly-implicit numerical scheme finite difference momentum equations as follows: (1) added the velocity flip-flop mass/energy error mitigation logic, (2) added the modified Henry-Fauske choking model, (3) used the new time void fraction in the horizontal stratification force terms and gravity head, and (4) used an implicit form of the artificial viscosity. The code modifications allow the nearly-implicit numerical scheme to be more implicit and lead to enhanced numerical stability.

  20. Stable and fast semi-implicit integration of the stochastic Landau-Lifshitz equation.

    PubMed

    Mentink, J H; Tretyakov, M V; Fasolino, A; Katsnelson, M I; Rasing, Th

    2010-05-01

    We propose new semi-implicit numerical methods for the integration of the stochastic Landau-Lifshitz equation with built-in angular momentum conservation. The performance of the proposed integrators is tested on the 1D Heisenberg chain. For this system, our schemes show better stability properties and allow us to use considerably larger time steps than standard explicit methods. At the same time, these semi-implicit schemes are also of comparable accuracy to and computationally much cheaper than the standard midpoint implicit method. The results are of key importance for atomistic spin dynamics simulations and the study of spin dynamics beyond the macro spin approximation. PMID:21393676

  1. On the computational efficiency of particle dynamics simulations in the radiation belt: comparison between implicit and semi-implicit schemes

    NASA Astrophysics Data System (ADS)

    Camporeale, E.; Delzanno, G.; Zaharia, S. G.; Koller, J.

    2012-12-01

    The particle dynamics in the Earth's radiation belt is generally modeled by means of a two-dimensional diffusion equation for the particle distribution function in energy and pitch angle. In this work we survey and compare different numerical schemes for the solution of the diffusion equation, with the goal of outlining which is the optimal strategy from a numerical point of view. We focus on the general (and more computationally challenging) case where the mixed terms in the diffusion tensor are retained. We compare fully-implicit and semi-implicit schemes. For the former we have analyzed a direct solver based on a LU decomposition routine for sparse matrices, and an iterative ILU-preconditioned GMRES. For the semi-implicit scheme we have studied an Alternating Direction Implicit scheme. We present a convergence study for a realistic case that shows that the timestep and grid size are strongly constrained by the desired accuracy of the solution. We show that the fully-implicit scheme is to be preferred in most cases as the more computationally efficient.

  2. On the numerical simulation of particle dynamics in the radiation belt: 1. Implicit and semi-implicit schemes

    NASA Astrophysics Data System (ADS)

    Camporeale, E.; Delzanno, G. L.; Zaharia, S.; Koller, J.

    2013-06-01

    The particle dynamics in the Earth's radiation belt is generally modeled by means of a two-dimensional diffusion equation for the particle distribution function in energy and pitch angle. The goal of this paper is to survey and compare different numerical schemes for the solution of the diffusion equation, and to outline the optimal strategy from a numerical point of view. We focus on the general (and more computationally challenging) case where the mixed terms in the diffusion tensor are retained. In Part 1, we compare fully implicit and semi-implicit schemes. For the former, we have analyzed a direct solver based on a LU decomposition routine for sparse matrices, and an iterative incomplete LU preconditioned Generalized Minimal REsidual solver. For the semi-implicit scheme, we have studied an alternating direction implicit scheme. We present a convergence study for a realistic case that shows that the time step and grid size are strongly constrained by the desired accuracy of the solution. We show that the fully implicit scheme is to be preferred in most cases as the more computationally efficient.

  3. Implicit integration methods for dislocation dynamics

    SciTech Connect

    Gardner, D. J.; Woodward, C. S.; Reynolds, D. R.; Hommes, G.; Aubry, S.; Arsenlis, A.

    2015-01-20

    In dislocation dynamics simulations, strain hardening simulations require integrating stiff systems of ordinary differential equations in time with expensive force calculations, discontinuous topological events, and rapidly changing problem size. Current solvers in use often result in small time steps and long simulation times. Faster solvers may help dislocation dynamics simulations accumulate plastic strains at strain rates comparable to experimental observations. Here, this paper investigates the viability of high order implicit time integrators and robust nonlinear solvers to reduce simulation run times while maintaining the accuracy of the computed solution. In particular, implicit Runge-Kutta time integrators are explored as a way of providing greater accuracy over a larger time step than is typically done with the standard second-order trapezoidal method. In addition, both accelerated fixed point and Newton's method are investigated to provide fast and effective solves for the nonlinear systems that must be resolved within each time step. Results show that integrators of third order are the most effective, while accelerated fixed point and Newton's method both improve solver performance over the standard fixed point method used for the solution of the nonlinear systems.

  4. Implicit integration methods for dislocation dynamics

    DOE PAGESBeta

    Gardner, D. J.; Woodward, C. S.; Reynolds, D. R.; Hommes, G.; Aubry, S.; Arsenlis, A.

    2015-01-20

    In dislocation dynamics simulations, strain hardening simulations require integrating stiff systems of ordinary differential equations in time with expensive force calculations, discontinuous topological events, and rapidly changing problem size. Current solvers in use often result in small time steps and long simulation times. Faster solvers may help dislocation dynamics simulations accumulate plastic strains at strain rates comparable to experimental observations. Here, this paper investigates the viability of high order implicit time integrators and robust nonlinear solvers to reduce simulation run times while maintaining the accuracy of the computed solution. In particular, implicit Runge-Kutta time integrators are explored as a waymore » of providing greater accuracy over a larger time step than is typically done with the standard second-order trapezoidal method. In addition, both accelerated fixed point and Newton's method are investigated to provide fast and effective solves for the nonlinear systems that must be resolved within each time step. Results show that integrators of third order are the most effective, while accelerated fixed point and Newton's method both improve solver performance over the standard fixed point method used for the solution of the nonlinear systems.« less

  5. Implicit integration methods for dislocation dynamics

    NASA Astrophysics Data System (ADS)

    Gardner, D. J.; Woodward, C. S.; Reynolds, D. R.; Hommes, G.; Aubry, S.; Arsenlis, A.

    2015-03-01

    In dislocation dynamics simulations, strain hardening simulations require integrating stiff systems of ordinary differential equations in time with expensive force calculations, discontinuous topological events and rapidly changing problem size. Current solvers in use often result in small time steps and long simulation times. Faster solvers may help dislocation dynamics simulations accumulate plastic strains at strain rates comparable to experimental observations. This paper investigates the viability of high-order implicit time integrators and robust nonlinear solvers to reduce simulation run times while maintaining the accuracy of the computed solution. In particular, implicit Runge-Kutta time integrators are explored as a way of providing greater accuracy over a larger time step than is typically done with the standard second-order trapezoidal method. In addition, both accelerated fixed point and Newton's method are investigated to provide fast and effective solves for the nonlinear systems that must be resolved within each time step. Results show that integrators of third order are the most effective, while accelerated fixed point and Newton's method both improve solver performance over the standard fixed point method used for the solution of the nonlinear systems.

  6. A Second-Order Iterative Implicit Explicit Hybrid Scheme for Hyperbolic Systems of Conservation Laws

    NASA Astrophysics Data System (ADS)

    Dai, Wenlong; Woodward, Paul R.

    1996-10-01

    An iterative implicit-explicit hybrid scheme is proposed for hyperbolic systems of conservation laws. Each wave in a system may be implicitly, or explicitly, or partially implicitly and partially explicitly treated depending on its associated Courant number in each numerical cell, and the scheme is able to smoothly switch between implicit and explicit calculations. The scheme is of Godunov-type in both explicit and implicit regimes, is in a strict conservation form, and is accurate to second-order in both space and time for all Courant numbers. The computer code for the scheme is easy to vectorize. Multicolors proposed in this paper may reduce the number of iterations required to reach a converged solution by several orders for a large time step. The feature of the scheme is shown through numerical examples.

  7. On a fourth order accurate implicit finite difference scheme for hyperbolic conservation laws. II - Five-point schemes

    NASA Technical Reports Server (NTRS)

    Harten, A.; Tal-Ezer, H.

    1981-01-01

    This paper presents a family of two-level five-point implicit schemes for the solution of one-dimensional systems of hyperbolic conservation laws, which generalized the Crank-Nicholson scheme to fourth order accuracy (4-4) in both time and space. These 4-4 schemes are nondissipative and unconditionally stable. Special attention is given to the system of linear equations associated with these 4-4 implicit schemes. The regularity of this system is analyzed and efficiency of solution-algorithms is examined. A two-datum representation of these 4-4 implicit schemes brings about a compactification of the stencil to three mesh points at each time-level. This compact two-datum representation is particularly useful in deriving boundary treatments. Numerical results are presented to illustrate some properties of the proposed scheme.

  8. Generalized formulation of a class of explicit and implicit TVD schemes

    NASA Technical Reports Server (NTRS)

    Yee, H. C.

    1985-01-01

    A one parameter family of second order explicit and implicit total variation diminishing (TVD) schemes is reformulated so that a simpler and wider group of limiters is included. The resulting scheme can be viewed as a symmetrical algorithm with a variety of numerical dissipation terms that are designed for weak solutions of hyperbolic problems. This is a generalization of Roe and Davis's recent works to a wider class of symmetric schemes other than Lax-Wendroff. The main properties of the present class of schemes are that they can be implicit, and when steady state calculations are sought, the numerical solution is independent of the time step.

  9. A Low-Dispersion and Low-Dissipation Implicit Runge-Kutta Scheme

    PubMed Central

    Najafi-Yazdi, A.; Mongeau, L.

    2012-01-01

    A fourth-order, implicit, low-dispersion, and low-dissipation Runge-Kutta scheme is introduced. The scheme is optimized for minimal dissipation and dispersion errors. High order accuracy is achieved with fewer stages than standard explicit Runge-Kutta schemes. The scheme is designed to be As table for highly stiff problems. Possible applications include wall-bounded flows with solid boundaries in the computational domain, and sound generation by reacting flows. PMID:23243319

  10. A Low-Dispersion and Low-Dissipation Implicit Runge-Kutta Scheme.

    PubMed

    Najafi-Yazdi, A; Mongeau, L

    2013-01-15

    A fourth-order, implicit, low-dispersion, and low-dissipation Runge-Kutta scheme is introduced. The scheme is optimized for minimal dissipation and dispersion errors. High order accuracy is achieved with fewer stages than standard explicit Runge-Kutta schemes. The scheme is designed to be As table for highly stiff problems. Possible applications include wall-bounded flows with solid boundaries in the computational domain, and sound generation by reacting flows. PMID:23243319

  11. Hybrid, explicit-implicit, finite-volume schemes on unstructured grids for unsteady compressible flows

    NASA Astrophysics Data System (ADS)

    Timofeev, Evgeny; Norouzi, Farhang

    2016-06-01

    The motivation for using hybrid, explicit-implicit, schemes rather than fully implicit or explicit methods for some unsteady high-speed compressible flows with shocks is firstly discussed. A number of such schemes proposed in the past are briefly overviewed. A recently proposed hybridization approach is then introduced and used for the development of a hybrid, explicit-implicit, TVD (Total Variation Diminishing) scheme of the second order in space and time on smooth solutions in both, explicit and implicit, modes for the linear advection equation. Further generalizations of this finite-volume method for the Burgers, Euler and Navier-Stokes equations discretized on unstructured grids are mentioned in the concluding remarks.

  12. Implicit Predictor-Corrector finite difference scheme for the ideal MHD simulations

    NASA Astrophysics Data System (ADS)

    Tsai, T.; Yu, H.; Lai, S.

    2012-12-01

    A innovative simulation code for ideal magnetohydrodynamics (MHD) is developed. We present a multiple-dimensional MHD code based on high-order implicit predictor-corrector finite difference scheme (high-order IPCFD scheme). High-order IPCFD scheme adopts high-order predictor-corrector scheme for the time integration and high-order central difference method as the spatial derivative solver. We use Elimination-of-the-Runoff-Errors (ERE) technology to avoid the numerical oscillations and numerical instability in the simulation results. In one-dimensional MHD problem, our simulation results show good agreement with the Brio & Wu MHD shock tube problem. The divergent B constraint remains fully satisfied, that is the divergent B equals to zero throughout the simulation. When solving the two-dimensional (2D) linear wave in MHD plasma, we clearly obtain the group-velocity Friedrichs diagrams of the MHD waves. Here we demonstrate 2D simulation results of rotor problem, Orszag-Tang vortex system, vortex type K-H instability, and kink type K-H instability by using our IPCFD MHD code and discuss the advantage of our simulation code.

  13. Numerical solution of 3D Navier-Stokes equations with upwind implicit schemes

    NASA Technical Reports Server (NTRS)

    Marx, Yves P.

    1990-01-01

    An upwind MUSCL type implicit scheme for the three-dimensional Navier-Stokes equations is presented. Comparison between different approximate Riemann solvers (Roe and Osher) are performed and the influence of the reconstructions schemes on the accuracy of the solution as well as on the convergence of the method is studied. A new limiter is introduced in order to remove the problems usually associated with non-linear upwind schemes. The implementation of a diagonal upwind implicit operator for the three-dimensional Navier-Stokes equations is also discussed. Finally the turbulence modeling is assessed. Good prediction of separated flows are demonstrated if a non-equilibrium turbulence model is used.

  14. An implicit factored scheme for the compressible Navier-Stokes equations. II - The numerical ODE connection

    NASA Technical Reports Server (NTRS)

    Beam, R. M.; Warming, R. F.

    1979-01-01

    An attempt is made to establish a connection between linear multistep methods for applications to ordinary differential equations and their extension (by approximate factorization) to alternating direction implicit methods for partial differential equations. An earlier implicit factored scheme for the compressible Navier-Stokes equations is generalized by innovations that (1) increase the class of temporal difference schemes to include all linear multistep methods, (2) optimize the class of unconditionally stable factored schemes by a new choice of unknown variable, and (3) improve the computational efficiency by the introduction of quasi-one-leg methods.

  15. Application of locally one-dimensional semi-implicit scheme in phase-field equations

    NASA Astrophysics Data System (ADS)

    Cai, Dan; Zhang, Lijun; Du, Yong

    2015-07-01

    A locally one-dimensional (LOD) semi-implicit scheme is proposed for improving the numerical efficiency in the solving of parabolic partial differential equations in phase-field simulations. With LOD splitting, multi-dimensional parabolic problems can be numerically approximated by treating each of the spatial variables individually in single cycles. Additionally, each spatial variable can be treated in either real or Fourier space, allowing equations to be solved across a range of boundary conditions, including periodic, non-periodic, and even partial periodic. The proposed LOD semi-implicit scheme exhibits noticeable advantages over both explicit and implicit traditional schemes in terms of computational efficiency and accuracy, as demonstrated by two standard numerical tests. It is anticipated that future large-scale phase-field simulations will benefit greatly from the use of this LOD scheme.

  16. Development of the Semi-implicit Time Integration in KIM-SH

    NASA Astrophysics Data System (ADS)

    NAM, H.

    2015-12-01

    The Korea Institute of Atmospheric Prediction Systems (KIAPS) was founded in 2011 by the Korea Meteorological Administration (KMA) to develop Korea's own global Numerical Weather Prediction (NWP) system as nine year (2011-2019) project. The KIM-SH is a KIAPS integrated model-spectral element based in the HOMME. In KIM-SH, the explicit schemes are employed. We introduce the three- and two-time-level semi-implicit scheme in KIM-SH as the time integration. Explicit schemes however have a tendancy to be unstable and require very small timesteps while semi-implicit schemes are very stable and can have much larger timesteps.We define the linear and reference values, then by definition of semi-implicit scheme, we apply the linear solver as GMRES. The numerical results from experiments will be introduced with the current development status of the time integration in KIM-SH. Several numerical examples are shown to confirm the efficiency and reliability of the proposed schemes.

  17. Recent advances in the development of implicit schemes for the equations of fluid dynamics

    NASA Technical Reports Server (NTRS)

    Warming, R. F.; Beam, R. M.

    1981-01-01

    Innovations and extensions of implicit schemes for equations of fluid dynamics are presented. The notation and theory for linear multistep methods are reviewed, and extensions of work by Beam and Warming (1979) include the implementation of one-leg methods, ADI methods for equations with mixed derivatives, flux vector splitting, the P-dimensional wave equation, and boundary conditions. Numerical experiments indicate that implicit treatment of the boundary conditions is necessary for unconditional stability, and the improvement and implementation of the boundary condition theory should improve the implicit algorithms for gas dynamic equations.

  18. TTLEM - an implicit-explicit (IMEX) scheme for modelling landscape evolution in MATLAB

    NASA Astrophysics Data System (ADS)

    Campforts, Benjamin; Schwanghart, Wolfgang

    2016-04-01

    Landscape evolution models (LEM) are essential to unravel interdependent earth surface processes. They are proven very useful to bridge several temporal and spatial timescales and have been successfully used to integrate existing empirical datasets. There is a growing consensus that landscapes evolve at least as much in the horizontal as in the vertical direction urging for an efficient implementation of dynamic drainage networks. Here we present a spatially explicit LEM, which is based on the object-oriented function library TopoToolbox 2 (Schwanghart and Scherler, 2014). Similar to other LEMs, rivers are considered to be the main drivers for simulated landscape evolution as they transmit pulses of tectonic perturbations and set the base level of surrounding hillslopes. Highly performant graph algorithms facilitate efficient updates of the flow directions to account for planform changes in the river network and the calculation of flow-related terrain attributes. We implement the model using an implicit-explicit (IMEX) scheme, i.e. different integrators are used for different terms in the diffusion-incision equation. While linear diffusion is solved using an implicit scheme, we calculate incision explicitly. Contrary to previously published LEMS, however, river incision is solved using a total volume method which is total variation diminishing in order to prevent numerical diffusion when solving the stream power law (Campforts and Govers, 2015). We show that the use of this updated numerical scheme alters both landscape topography and catchment wide erosion rates at a geological time scale. Finally, the availability of a graphical user interface facilitates user interaction, making the tool very useful both for research and didactical purposes. References Campforts, B., Govers, G., 2015. Keeping the edge: A numerical method that avoids knickpoint smearing when solving the stream power law. J. Geophys. Res. Earth Surf. 120, 1189-1205. doi:10.1002/2014JF003376

  19. Integrating Implicit Induction Proofs into Certified Proof Environments

    NASA Astrophysics Data System (ADS)

    Stratulat, Sorin

    We give evidence of the direct integration and automated checking of implicit induction-based proofs inside certified reasoning environments, as that provided by the Coq proof assistant. This is the first step of a long term project focused on 1) mechanically certifying implicit induction proofs generated by automated provers like Spike, and 2) narrowing the gap between automated and interactive proof techniques inside proof assistants such that multiple induction steps can be executed completely automatically and mutual induction can be treated more conveniently. Contrary to the current approaches of reconstructing implicit induction proofs into scripts based on explicit induction tactics that integrate the usual proof assistants, our checking methodology is simpler and fits better for automation. The underlying implicit induction principles are separated and validated independently from the proof scripts that consist in a bunch of one-to-one translations of implicit induction proof steps. The translated steps can be checked independently, too, so the validation process fits well for parallelisation and for the management of large proof scripts. Moreover, our approach is more general; any kind of implicit induction proof can be considered because the limitations imposed by the proof reconstruction techniques no longer exist. An implementation that integrates automatic translators for generating fully checkable Coq scripts from Spike proofs is reported.

  20. Implicit unified gas-kinetic scheme for steady state solutions in all flow regimes

    NASA Astrophysics Data System (ADS)

    Zhu, Yajun; Zhong, Chengwen; Xu, Kun

    2016-06-01

    This paper presents an implicit unified gas-kinetic scheme (UGKS) for non-equilibrium steady state flow computation. The UGKS is a direct modeling method for flow simulation in all regimes with the updates of both macroscopic flow variables and microscopic gas distribution function. By solving the macroscopic equations implicitly, a predicted equilibrium state can be obtained first through iterations. With the newly predicted equilibrium state, the evolution equation of the gas distribution function and the corresponding collision term can be discretized in a fully implicit way for fast convergence through iterations as well. The lower-upper symmetric Gauss-Seidel (LU-SGS) factorization method is implemented to solve both macroscopic and microscopic equations, which improves the efficiency of the scheme. Since the UGKS is a direct modeling method and its physical solution depends on the mesh resolution and the local time step, a physical time step needs to be fixed before using an implicit iterative technique with a pseudo-time marching step. Therefore, the physical time step in the current implicit scheme is determined by the same way as that in the explicit UGKS for capturing the physical solution in all flow regimes, but the convergence to a steady state speeds up through the adoption of a numerical time step with large CFL number. Many numerical test cases in different flow regimes from low speed to hypersonic ones, such as the Couette flow, cavity flow, and the flow passing over a cylinder, are computed to validate the current implicit method. The overall efficiency of the implicit UGKS can be improved by one or two orders of magnitude in comparison with the explicit one.

  1. Massively Parallel and Scalable Implicit Time Integration Algorithms for Structural Dynamics

    NASA Technical Reports Server (NTRS)

    Farhat, Charbel

    1997-01-01

    Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because of the following additional facts: (a) explicit schemes are easier to parallelize than implicit ones, and (b) explicit schemes induce short range interprocessor communications that are relatively inexpensive, while the factorization methods used in most implicit schemes induce long range interprocessor communications that often ruin the sought-after speed-up. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet be offset by the speed of the currently available parallel hardware. Therefore, it is essential to develop efficient alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient when simulating the low-frequency dynamics of aerospace structures.

  2. Implicit integration in a case of integrative visual agnosia

    PubMed Central

    Aviezer, Hillel; Landau, Ayelet N.; Robertson, Lynn C.; Peterson, Mary A.; Soroker, Nachum; Sacher, Yaron; Bonneh, Yoram; Bentin, Shlomo

    2007-01-01

    We present a case (SE) with integrative visual agnosia following ischemic stroke affecting the right dorsal and the left ventral pathways of the visual system. Despite his inability to identify global hierarchical letters (Navon, 1977), and his dense object agnosia, SE showed normal global-to-local interference when responding to local letters in Navon hierarchical stimuli and significant picture-word identity priming in a semantic decision task for words. Since priming was absent if these features were scrambled, it stands to reason that these effects were not due to priming by distinctive features. The contrast between priming effects induced by coherent and scrambled stimuli is consistent with implicit but not explicit integration of features into a unified whole. We went on to show that possible/impossible object decisions were facilitated by words in a word-picture priming task, suggesting that prompts could activate perceptually integrated images in a backward fashion. We conclude that the absence of SE's ability to identify visual objects except through tedious serial construction reflects a deficit in accessing an integrated visual representation through bottom-up visual processing alone. However, top-down generated images can help activate these visual representations through semantic links. PMID:17339044

  3. Application of the implicit MacCormack scheme to the parabolized Navier-Stokes equations

    NASA Technical Reports Server (NTRS)

    Lawrence, J. L.; Tannehill, J. C.; Chaussee, D. S.

    1984-01-01

    MacCormack's implicit finite-difference scheme was used to solve the two-dimensional parabolized Navier-Stokes (PNS) equations. This method for solving the PNS equations does not require the inversion of block tridiagonal systems of algebraic equations and permits the original explicit MacCormack scheme to be employed in those regions where implicit treatment is not needed. The advantages and disadvantages of the present adaptation are discussed in relation to those of the conventional Beam-Warming scheme for a flat plate boundary layer test case. Comparisons are made for accuracy, stability, computer time, computer storage, and ease of implementation. The present method was also applied to a second test case of hypersonic laminar flow over a 15% compression corner. The computed results compare favorably with experiment and a numerical solution of the complete Navier-Stokes equations.

  4. A diagonal implicit scheme for computing flows with finite-rate chemistry

    NASA Technical Reports Server (NTRS)

    Eberhardt, Scott; Imlay, Scott

    1990-01-01

    A new algorithm for solving steady, finite-rate chemistry, flow problems is presented. The new scheme eliminates the expense of inverting large block matrices that arise when species conservation equations are introduced. The source Jacobian matrix is replaced by a diagonal matrix which is tailored to account for the fastest reactions in the chemical system. A point-implicit procedure is discussed and then the algorithm is included into the LU-SGS scheme. Solutions are presented for hypervelocity reentry and Hydrogen-Oxygen combustion. For the LU-SGS scheme a CFL number in excess of 10,000 has been achieved.

  5. Stable boundary approximations for a class of implicit schemes for the one-dimensional inviscid equations of gas dynamics

    NASA Technical Reports Server (NTRS)

    Yee, H. C.; Beam, R. M.; Warming, R. F.

    1981-01-01

    The applicability to practical calculations of recent theoretical developments in the stability analysis of difference approximations for initial-boundary-value problems of the hyperbolic type. For the numerical experiments, select the one-dimensional inviscid gas-dynamic equations in conservation-law form is selected. A class of implicit schemes based on linear multistep methods for ordinary differential equations is chosen and the use of space or space-time extrapolations as implicit or explicit boundary schemes is emphasized. Some numerical examples with various inflow-outflow conditions highlight the commonly discussed issues: explicit versus implicit boundary schemes, unconditionally stable schemes, and underspecification or overspecification of boundary conditions.

  6. Robust Integration Schemes for Generalized Viscoplasticity with Internal-State Variables

    NASA Technical Reports Server (NTRS)

    Saleeb, Atef F.; Li, W.; Wilt, Thomas E.

    1997-01-01

    The scope of the work in this presentation focuses on the development of algorithms for the integration of rate dependent constitutive equations. In view of their robustness; i.e., their superior stability and convergence properties for isotropic and anisotropic coupled viscoplastic-damage models, implicit integration schemes have been selected. This is the simplest in its class and is one of the most widely used implicit integrators at present.

  7. An implicit compact scheme solver with application to chemically reacting flows

    NASA Astrophysics Data System (ADS)

    Noskov, Mikhail; Smooke, Mitchell D.

    2005-03-01

    A novel, stable, implicit compact scheme solver that is higher order in space, suitable for modeling steady-state and time-dependent phenomena on nonuniform grids for one-dimensional configurations, is presented. Several properties of compact scheme discretizations are introduced to develop efficient algorithms for Jacobian matrix generation and Jacobian-vector multiplication using a new component form for Jacobian operations. Composite nonuniform grids are introduced that enable the implicit compact scheme solver to achieve sixth order accuracy. A robust Newton's method is employed with explicit generation of Jacobian matrices. Superior resolution characteristics of the implicit compact scheme solver are demonstrated with several steady-state and time-dependent problems for the Burgers equation. The example of the solution of stiff flame problem is given. An analysis of spectral properties of Jacobian matrices is presented, which shows that the condition number and the eigenvalue distributions behave similarly to those found in Jacobians associated with low-order discretizations. Two sparsification strategies are developed for the systematic approximation of a dense Jacobian aimed at the practical implementation of linear system preconditioning through partial Jacobians.

  8. Implicit flux-split Euler schemes for unsteady aerodynamic analysis involving unstructured dynamic meshes

    NASA Technical Reports Server (NTRS)

    Batina, John T.

    1990-01-01

    Improved algorithm for the solution of the time-dependent Euler equations are presented for unsteady aerodynamic analysis involving unstructured dynamic meshes. The improvements were developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves an implicit time-integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. For example, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady and unsteady flow results are presented for the NACA 0012 airfoil to demonstrate applications of the new Euler solvers. The unsteady results were obtained for the airfoil pitching harmonically about the quarter chord. The resulting instantaneous pressure distributions and lift and moment coefficients during a cycle of motion compare well with experimental data. A description of the Euler solvers is presented along with results and comparisons which assess the capability.

  9. A semi-implicit gas-kinetic scheme for smooth flows

    NASA Astrophysics Data System (ADS)

    Wang, Peng; Guo, Zhaoli

    2016-08-01

    In this paper, a semi-implicit gas-kinetic scheme (SIGKS) is derived for smooth flows based on the Bhatnagar-Gross-Krook (BGK) equation. As a finite-volume scheme, the evolution of the average flow variables in a control volume is under the Eulerian framework, whereas the construction of the numerical flux across the cell interface comes from the Lagrangian perspective. The adoption of the Lagrangian aspect makes the collision and the transport mechanisms intrinsically coupled together in the flux evaluation. As a result, the time step size is independent of the particle collision time and solely determined by the Courant-Friedrichs-Lewy (CFL) condition. An analysis of the reconstructed distribution function at the cell interface shows that the SIGKS can be viewed as a modified Lax-Wendroff type scheme with an additional term. Furthermore, the addition term coming from the implicitness in the reconstruction is expected to be able to enhance the numerical stability of the scheme. A number of numerical tests of smooth flows with low and moderate Mach numbers are performed to benchmark the SIGKS. The results show that the method has second-order spatial accuracy, and can give accurate numerical solutions in comparison with benchmark results. It is also demonstrated that the numerical stability of the proposed scheme is better than the original GKS for smooth flows.

  10. Implicit scheme for Maxwell equations solution in case of flat 3D domains

    NASA Astrophysics Data System (ADS)

    Boronina, Marina; Vshivkov, Vitaly

    2016-02-01

    We present a new finite-difference scheme for Maxwell's equations solution for three-dimensional domains with different scales in different directions. The stability condition of the standard leap-frog scheme requires decreasing of the time-step with decreasing of the minimal spatial step, which depends on the minimal domain size. We overcome the conditional stability by modifying the standard scheme adding implicitness in the direction of the smallest size. The new scheme satisfies the Gauss law for the electric and magnetic fields in the final- differences. The approximation order, the maintenance of the wave amplitude and propagation speed, the invariance of the wave propagation on angle with the coordinate axes are analyzed.

  11. A diagonally inverted LU implicit multigrid scheme for the 3-D Navier-Stokes equations and a two equation model of turbulence

    NASA Technical Reports Server (NTRS)

    Yokota, Jeffrey W.

    1988-01-01

    An LU implicit multigrid algorithm is developed to calculate 3-D compressible viscous flows. This scheme solves the full 3-D Reynolds-Averaged Navier-Stokes equation with a two-equation kappa-epsilon model of turbulence. The flow equations are integrated by an efficient, diagonally inverted, LU implicit multigrid scheme while the kappa-epsilon equations are solved, uncoupled from the flow equations, by a block LU implicit algorithm. The flow equations are solved within the framework of the multigrid method using a four-grid level W-cycle, while the kappa-epsilon equations are iterated only on the finest grid. This treatment of the Reynolds-Averaged Navier-Stokes equations proves to be an efficient method for calculating 3-D compressible viscous flows.

  12. A matrix free implicit scheme for solution of resistive magneto-hydrodynamics equations on unstructured grids

    NASA Astrophysics Data System (ADS)

    Sitaraman, H.; Raja, L. L.

    2013-10-01

    The resistive magneto-hydrodynamics (MHD) governing equations represent eight conservation equations for the evolution of density, momentum, energy and induced magnetic fields in an electrically conducting fluid, typically a plasma. A matrix free implicit method is developed to solve the conservation equations within the framework of an unstructured grid finite volume formulation. The analytic form of the convective flux Jacobian is derived on a general unstructured mesh and used in a Lower-Upper Symmetric Gauss Seidel (LU-SGS) technique developed as part of the implicit scheme. A grid coloring technique is also developed to create data parallelism in the algorithm. The computational efficiency of the matrix free method is compared with two common approaches: a global matrix solve technique that uses the GMRES (Generalized minimum residual) algorithm and an explicit method. The matrix-free method is observed to be overall computationally faster than the global matrix solve method and demonstrates excellent parallel scaling on multiple cores. The computational effort and memory requirements for the matrix free approach is comparable to the explicit approach which in turn is much lower than the global solve implicit approach. Both the matrix free and global solve implicit techniques exhibit superior steady state convergence compared to the explicit method.

  13. Stability of mixed time integration schemes for transient thermal analysis

    NASA Technical Reports Server (NTRS)

    Liu, W. K.; Lin, J. I.

    1982-01-01

    A current research topic in coupled-field problems is the development of effective transient algorithms that permit different time integration methods with different time steps to be used simultaneously in various regions of the problems. The implicit-explicit approach seems to be very successful in structural, fluid, and fluid-structure problems. This paper summarizes this research direction. A family of mixed time integration schemes, with the capabilities mentioned above, is also introduced for transient thermal analysis. A stability analysis and the computer implementation of this technique are also presented. In particular, it is shown that the mixed time implicit-explicit methods provide a natural framework for the further development of efficient, clean, modularized computer codes.

  14. An implicit midpoint difference scheme for the fractional Ginzburg-Landau equation

    NASA Astrophysics Data System (ADS)

    Wang, Pengde; Huang, Chengming

    2016-05-01

    This paper proposes and analyzes an efficient difference scheme for the nonlinear complex Ginzburg-Landau equation involving fractional Laplacian. The scheme is based on the implicit midpoint rule for the temporal discretization and a weighted and shifted Grünwald difference operator for the spatial fractional Laplacian. By virtue of a careful analysis of the difference operator, some useful inequalities with respect to suitable fractional Sobolev norms are established. Then the numerical solution is shown to be bounded, and convergent in the lh2 norm with the optimal order O (τ2 +h2) with time step τ and mesh size h. The a priori bound as well as the convergence order holds unconditionally, in the sense that no restriction on the time step τ in terms of the mesh size h needs to be assumed. Numerical tests are performed to validate the theoretical results and effectiveness of the scheme.

  15. Implicit approximate-factorization schemes for the low-frequency transonic equation

    NASA Technical Reports Server (NTRS)

    Ballhaus, W. F.; Steger, J. L.

    1975-01-01

    Two- and three-level implicit finite-difference algorithms for the low-frequency transonic small disturbance-equation are constructed using approximate factorization techniques. The schemes are unconditionally stable for the model linear problem. For nonlinear mixed flows, the schemes maintain stability by the use of conservatively switched difference operators for which stability is maintained only if shock propagation is restricted to be less than one spatial grid point per time step. The shock-capturing properties of the schemes were studied for various shock motions that might be encountered in problems of engineering interest. Computed results for a model airfoil problem that produces a flow field similar to that about a helicopter rotor in forward flight show the development of a shock wave and its subsequent propagation upstream off the front of the airfoil.

  16. An Initial Investigation of the Effects of Turbulence Models on the Convergence of the RK/Implicit Scheme

    NASA Technical Reports Server (NTRS)

    Swanson, R. C.; Rossow, C.-C.

    2008-01-01

    A three-stage Runge-Kutta (RK) scheme with multigrid and an implicit preconditioner has been shown to be an effective solver for the fluid dynamic equations. This scheme has been applied to both the compressible and essentially incompressible Reynolds-averaged Navier-Stokes (RANS) equations using the algebraic turbulence model of Baldwin and Lomax (BL). In this paper we focus on the convergence of the RK/implicit scheme when the effects of turbulence are represented by either the Spalart-Allmaras model or the Wilcox k-! model, which are frequently used models in practical fluid dynamic applications. Convergence behavior of the scheme with these turbulence models and the BL model are directly compared. For this initial investigation we solve the flow equations and the partial differential equations of the turbulence models indirectly coupled. With this approach we examine the convergence behavior of each system. Both point and line symmetric Gauss-Seidel are considered for approximating the inverse of the implicit operator of the flow solver. To solve the turbulence equations we use a diagonally dominant alternating direction implicit (DDADI) scheme. Computational results are presented for three airfoil flow cases and comparisons are made with experimental data. We demonstrate that the two-dimensional RANS equations and transport-type equations for turbulence modeling can be efficiently solved with an indirectly coupled algorithm that uses the RK/implicit scheme for the flow equations.

  17. Krylov single-step implicit integration factor WENO methods for advection-diffusion-reaction equations

    NASA Astrophysics Data System (ADS)

    Jiang, Tian; Zhang, Yong-Tao

    2016-04-01

    Implicit integration factor (IIF) methods were developed in the literature for solving time-dependent stiff partial differential equations (PDEs). Recently, IIF methods were combined with weighted essentially non-oscillatory (WENO) schemes in Jiang and Zhang (2013) [19] to efficiently solve stiff nonlinear advection-diffusion-reaction equations. The methods can be designed for arbitrary order of accuracy. The stiffness of the system is resolved well and the methods are stable by using time step sizes which are just determined by the non-stiff hyperbolic part of the system. To efficiently calculate large matrix exponentials, Krylov subspace approximation is directly applied to the implicit integration factor (IIF) methods. So far, the IIF methods developed in the literature are multistep methods. In this paper, we develop Krylov single-step IIF-WENO methods for solving stiff advection-diffusion-reaction equations. The methods are designed carefully to avoid generating positive exponentials in the matrix exponentials, which is necessary for the stability of the schemes. We analyze the stability and truncation errors of the single-step IIF schemes. Numerical examples of both scalar equations and systems are shown to demonstrate the accuracy, efficiency and robustness of the new methods.

  18. Implicit predictor-corrector central finite difference scheme for the equations of magnetohydrodynamic simulations

    NASA Astrophysics Data System (ADS)

    Tsai, T. C.; Yu, H.-S.; Hsieh, M.-S.; Lai, S. H.; Yang, Y.-H.

    2015-11-01

    Nowadays most of supercomputers are based on the frame of PC cluster; therefore, the efficiency of parallel computing is of importance especially with the increasing computing scale. This paper proposes a high-order implicit predictor-corrector central finite difference (iPCCFD) scheme and demonstrates its high efficiency in parallel computing. Of special interests are the large scale numerical studies such as the magnetohydrodynamic (MHD) simulations in the planetary magnetosphere. An iPCCFD scheme is developed based on fifth-order central finite difference method and fourth-order implicit predictor-corrector method in combination with elimination-of-the-round-off-errors (ERE) technique. We examine several numerical studies such as one-dimensional Brio-Wu shock tube problem, two-dimensional Orszag-Tang vortex system, vortex type K-H instability, kink type K-H instability, field loop advection, and blast wave. All the simulation results are consistent with many literatures. iPCCFD can minimize the numerical instabilities and noises along with the additional diffusion terms. All of our studies present relatively small numerical errors without employing any divergence-free reconstruction. In particular, we obtain fairly stable results in the two-dimensional Brio-Wu shock tube problem which well conserves ∇ ṡ B = 0 throughout the simulation. The ERE technique removes the accumulation of roundoff errors in the uniform or non-disturbed system. We have also shown that iPCCFD is characterized by the high order of accuracy and the low numerical dissipation in the circularly polarized Alfvén wave tests. The proposed iPCCFD scheme is a parallel-efficient and high precision numerical scheme for solving the MHD equations in hyperbolic conservation systems.

  19. Resonance in the dynamics of chemical systems simulated by the implicit midpoint scheme

    NASA Astrophysics Data System (ADS)

    Mandziuk, Margaret; Schlick, Tamar

    1995-05-01

    The numerical behavior of the symplectic, implicit midpoint method with a wide range of integration timesteps is examined through an application to a diatomic molecule governed by a Morse potential. Our oscillator with a 12.6 fs period exhibits notable, integrator induced, timestep- ( Δt) dependent resonances and we predict approximate values of Δt where they will occur. The particular case of a third-order resonance ( Δt ≈ 7 fs here) leads to instability, and higher-order resonances ( n = 4, 5) to large energetic fluctuations and/or corrupted phase diagrams. Significantly, for Δt > 10 fs the energy errors remain bound.

  20. Low-storage implicit/explicit Runge-Kutta schemes for the simulation of stiff high-dimensional ODE systems

    NASA Astrophysics Data System (ADS)

    Cavaglieri, Daniele; Bewley, Thomas

    2015-04-01

    Implicit/explicit (IMEX) Runge-Kutta (RK) schemes are effective for time-marching ODE systems with both stiff and nonstiff terms on the RHS; such schemes implement an (often A-stable or better) implicit RK scheme for the stiff part of the ODE, which is often linear, and, simultaneously, a (more convenient) explicit RK scheme for the nonstiff part of the ODE, which is often nonlinear. Low-storage RK schemes are especially effective for time-marching high-dimensional ODE discretizations of PDE systems on modern (cache-based) computational hardware, in which memory management is often the most significant computational bottleneck. In this paper, we develop and characterize eight new low-storage implicit/explicit RK schemes which have higher accuracy and better stability properties than the only low-storage implicit/explicit RK scheme available previously, the venerable second-order Crank-Nicolson/Runge-Kutta-Wray (CN/RKW3) algorithm that has dominated the DNS/LES literature for the last 25 years, while requiring similar storage (two, three, or four registers of length N) and comparable floating-point operations per timestep.

  1. Extension of the Nambu Jona-Lasinio model predictions at high densities and temperatures using an implicit regularization scheme

    NASA Astrophysics Data System (ADS)

    Farias, R. L. S.; Dallabona, G.; Krein, G.; Battistel, O. A.

    2008-06-01

    Traditional cutoff regularization schemes of the Nambu Jona-Lasinio model limit the applicability of the model to energy-momentum scales much below the value of the regularizing cutoff. In particular, the model cannot be used to study quark matter with Fermi momenta larger than the cutoff. In the present work, an extension of the model to high temperatures and densities recently proposed by Casalbuoni, Gatto, Nardulli, and Ruggieri is used in connection with an implicit regularization scheme. This is done by making use of scaling relations of the divergent one-loop integrals that relate these integrals at different energy-momentum scales. Fixing the pion decay constant at the chiral symmetry breaking scale in the vacuum, the scaling relations predict a running coupling constant that decreases as the regularization scale increases, implementing in a schematic way the property of asymptotic freedom of quantum chromodynamics. If the regularization scale is allowed to increase with density and temperature, the coupling will decrease with density and temperature, extending in this way the applicability of the model to high densities and temperatures. These results are obtained without specifying an explicit regularization. As an illustration of the formalism, numerical results are obtained for the finite density and finite temperature quark condensate and applied to the problem of color superconductivity at high quark densities and finite temperature.

  2. Implicit Affective Cues and Attentional Tuning: An Integrative Review

    PubMed Central

    Friedman, Ronald S.; Förster, Jens

    2010-01-01

    A large and growing number of studies support the notion that arousing positive emotional states expand, and that arousing negative states constrict, the scope of attention on both the perceptual and conceptual levels. However, these studies have predominantly involved the manipulation or measurement of conscious emotional experiences (e.g., subjective feelings of happiness or anxiety). This raises the question: Do cues that are merely associated with benign versus threatening situations, but that do not elicit conscious feelings of positive or negative emotional arousal, independently expand or contract attentional scope? Integrating theoretical advances in affective neuroscience, positive psychology, and social cognition, it is proposed that rudimentary intero- and exteroceptive stimuli may indeed become associated with the onset of arousing positive or negative emotional states and/or with appraisals that the environment is benign or threatening and thereby come to moderate the scope of attention in the absence of conscious emotional experience. Specifically, implicit “benign situation” cues are posited to broaden, and implicit “threatening situation” cues to narrow, the range of both perceptual as well as conceptual attentional selection. An extensive array of research findings involving a diverse set of such implicit affective cues (e.g., enactment of approach and avoidance behaviors, incidental exposure to colors signaling safety versus danger) is marshaled in support of this proposition. Potential alternative explanations for and moderators of these attentional tuning effects, as well as their higher-level neuropsychological underpinnings, are also discussed along with prospective extensions to a range of other situational cues and domains of social cognitive processing. PMID:20804240

  3. An efficient approximate factorization implicit scheme for the equations of gasdynamics

    NASA Technical Reports Server (NTRS)

    Barth, T. J.; Steger, J. L.

    1984-01-01

    An efficient implicit finite-difference algorithm for the gas dynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved while maintaining the same favorable stability characteristics and generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows about a NACA 0012 airfoil. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations. It is shown that the method offers advantages over the conventional Beam and Warming scheme and can retrofit existing Beam and Warming codes with minimal effort.

  4. Quasi-static response, implicit scheme and incremental problem in gradient plasticity

    NASA Astrophysics Data System (ADS)

    Nguyen, Quoc-Son

    2016-06-01

    This paper is devoted to the study of gradient plasticity at small strains. Some time-independent dissipative processes such as brittle damage can also be considered in the same framework. Our attention is focussed on the description of the constitutive equations, on the formulation of the governing equations in terms of the energy potential and the dissipation potential of the solid. A time-discretization by the implicit scheme of the evolution equation leads to the study of the incremental problem which is different from the rate problem. The increment of the response under an increment of the loads must satisfy a variational inequality and, if the energy potential is convex, an incremental minimum principle. In particular, a local minimum of the incremental minimum principle is a stable solution to the variational inequality.

  5. Convergence acceleration of implicit schemes in the presence of high aspect ratio grid cells

    NASA Technical Reports Server (NTRS)

    Buelow, B. E. O.; Venkateswaran, S.; Merkle, C. L.

    1993-01-01

    of the scheme. On the other hand, for implicit schemes, which are typically unconditionally stable, there appears to be room for improvement through careful tailoring of the time step definition based on results of linear stability analyses. In the present paper, we focus on the central-differenced alternating direction implicit (ADI) scheme. The understanding garnered from this analyses can then be applied to other implicit schemes. In order to systematically study the effects of aspect ratio and the methods of mitigating the associated problems, we use a two pronged approach. We use stability analyses as a tool for predicting numerical convergence behavior and numerical experiments on simple model problems to verify predicted trends. Based on these analyses, we determine that efficient convergence may be obtained at all aspect ratios by getting a combination of things right. Primary among these are the proper definition of the time step size, proper selection of viscous preconditioner and the precise treatment of boundary conditions. These algorithmic improvements are then applied to a variety of test cases to demonstrate uniform convergence at all aspect ratios.

  6. Convergence acceleration of implicit schemes in the presence of high aspect ratio grid cells

    NASA Astrophysics Data System (ADS)

    Buelow, B. E. O.; Venkateswaran, S.; Merkle, C. L.

    1993-07-01

    of the scheme. On the other hand, for implicit schemes, which are typically unconditionally stable, there appears to be room for improvement through careful tailoring of the time step definition based on results of linear stability analyses. In the present paper, we focus on the central-differenced alternating direction implicit (ADI) scheme. The understanding garnered from this analyses can then be applied to other implicit schemes. In order to systematically study the effects of aspect ratio and the methods of mitigating the associated problems, we use a two pronged approach. We use stability analyses as a tool for predicting numerical convergence behavior and numerical experiments on simple model problems to verify predicted trends. Based on these analyses, we determine that efficient convergence may be obtained at all aspect ratios by getting a combination of things right.

  7. An implicit finite volume nodal point scheme for the solution of two-dimensional compressible Navier-Stokes equations

    NASA Astrophysics Data System (ADS)

    Dutta, Vimala

    1993-07-01

    An implicit finite volume nodal point scheme has been developed for solving the two-dimensional compressible Navier-Stokes equations. The numerical scheme is evolved by efficiently combining the basic ideas of the implicit finite-difference scheme of Beam and Warming (1978) with those of nodal point schemes due to Hall (1985) and Ni (1982). The 2-D Navier-Stokes solver is implemented for steady, laminar/turbulent flows past airfoils by using C-type grids. Turbulence closure is achieved by employing the algebraic eddy-viscosity model of Baldwin and Lomax (1978). Results are presented for the NACA-0012 and RAE-2822 airfoil sections. Comparison of the aerodynamic coefficients with experimental results for the different test cases presented here establishes the validity and efficiency of the method.

  8. An implicit numerical scheme for the simulation of internal viscous flows on unstructured grids

    NASA Technical Reports Server (NTRS)

    Jorgenson, Philip C. E.; Pletcher, Richard H.

    1994-01-01

    The Navier-Stokes equations are solved numerically for two-dimensional steady viscous laminar flows. The grids are generated based on the method of Delaunay triangulation. A finite-volume approach is used to discretize the conservation law form of the compressible flow equations written in terms of primitive variables. A preconditioning matrix is added to the equations so that low Mach number flows can be solved economically. The equations are time marched using either an implicit Gauss-Seidel iterative procedure or a solver based on a conjugate gradient like method. A four color scheme is employed to vectorize the block Gauss-Seidel relaxation procedure. This increases the memory requirements minimally and decreases the computer time spent solving the resulting system of equations substantially. A factor of 7.6 speed up in the matrix solver is typical for the viscous equations. Numerical results are obtained for inviscid flow over a bump in a channel at subsonic and transonic conditions for validation with structured solvers. Viscous results are computed for developing flow in a channel, a symmetric sudden expansion, periodic tandem cylinders in a cross-flow, and a four-port valve. Comparisons are made with available results obtained by other investigators.

  9. Application of a lower-upper implicit scheme and an interactive grid generation for turbomachinery flow field simulations

    NASA Technical Reports Server (NTRS)

    Choo, Yung K.; Soh, Woo-Yung; Yoon, Seokkwan

    1989-01-01

    A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a turbine cascade. This approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to generate grids. This program uses the control point form of algebraic grid generation which uses a sparse collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct advantage of TURBO compared with other grid generation programs is that it allows the easy change of local mesh structure without affecting the grid outside the domain of dependence. Sample grids are generated by TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the LU implicit scheme on the grid generated by TURBO.

  10. Application of a lower-upper implicit scheme and an interactive grid generation for turbomachinery flow field simulations

    NASA Technical Reports Server (NTRS)

    Choo, Yung K.; Soh, Woo-Yung; Yoon, Seokkwan

    1989-01-01

    A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a tubine cascade. This approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to generate grids. This program uses the control point form of algebraic grid generation which uses a sparse collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct advantage of TURBO compared with other grid generation programs is that it allows the easy change of local mesh structure without affecting the grid outside the domain of independence. Sample grids are generated by TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the LU implicit scheme on the grid generated by TURBO.

  11. Associative and Propositional Processes in Evaluation: An Integrative Review of Implicit and Explicit Attitude Change

    ERIC Educational Resources Information Center

    Gawronski, Bertram; Bodenhausen, Galen V.

    2006-01-01

    A central theme in recent research on attitudes is the distinction between deliberate, "explicit" attitudes and automatic, "implicit" attitudes. The present article provides an integrative review of the available evidence on implicit and explicit attitude change that is guided by a distinction between associative and propositional processes.…

  12. Implicit Time Integration for Multiscale Molecular Dynamics Using Transcendental Padé Approximants.

    PubMed

    Abi Mansour, Andrew; Ortoleva, Peter J

    2016-04-12

    Molecular dynamics systems evolve through the interplay of collective and localized disturbances. As a practical consequence, there is a restriction on the time step imposed by the broad spectrum of time scales involved. To resolve this restriction, multiscale factorization was introduced for molecular dynamics as a method that exploits the separation of time scales by coevolving the coarse-grained and atom-resolved states via Trotter factorization. Developing a stable time-marching scheme for this coevolution, however, is challenging because the coarse-grained dynamical equations depend on the microstate; therefore, these equations cannot be expressed in closed form. The objective of this paper is to develop an implicit time integration scheme for multiscale simulation of large systems over long periods of time and with high accuracy. The scheme uses Padé approximants to account for both the stochastic and deterministic features of the coarse-grained dynamics. The method is demonstrated for a protein either undergoing a conformational change or migrating under the influence of an external force. The method shows promise in accelerating multiscale molecular dynamics without a loss of atomic precision or the need to conjecture the form of coarse-grained governing equations. PMID:26845510

  13. Implicit flux-split Euler schemes for unsteady aerodynamic analysis involving unstructured dynamic meshes

    NASA Technical Reports Server (NTRS)

    Batina, John T.

    1990-01-01

    Improved algorithms for the solution of the time-dependent Euler equations are presented for unsteady aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves an implicit time-integration shceme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. For example, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady and unsteady flow results are presented for the NACA 0012 airfoil to demonstrate applications of the new Euler solvers. The unsteady results were obtained for the airfoil pitching harmonically about the quarter chord. The resulting instantaneous pressure distributions and lift and moment coefficients during a cycle of motion compare well with experimental data. The paper presents a description of the Euler solvers along with results and comparisons which assess the capability.

  14. A point implicit time integration technique for slow transient flow problems

    SciTech Connect

    Kadioglu, Samet Y.; Berry, Ray A.; Martineau, Richard C.

    2015-05-01

    We introduce a point implicit time integration technique for slow transient flow problems. The method treats the solution variables of interest (that can be located at cell centers, cell edges, or cell nodes) implicitly and the rest of the information related to same or other variables are handled explicitly. The method does not require implicit iteration; instead it time advances the solutions in a similar spirit to explicit methods, except it involves a few additional function(s) evaluation steps. Moreover, the method is unconditionally stable, as a fully implicit method would be. This new approach exhibits the simplicity of implementation of explicit methods and the stability of implicit methods. It is specifically designed for slow transient flow problems of long duration wherein one would like to perform time integrations with very large time steps. Because the method can be time inaccurate for fast transient problems, particularly with larger time steps, an appropriate solution strategy for a problem that evolves from a fast to a slow transient would be to integrate the fast transient with an explicit or semi-implicit technique and then switch to this point implicit method as soon as the time variation slows sufficiently. We have solved several test problems that result from scalar or systems of flow equations. Our findings indicate the new method can integrate slow transient problems very efficiently; and its implementation is very robust.

  15. Implicit Personality Theory in the Classroom: An Integrative Approach.

    ERIC Educational Resources Information Center

    Embree, Marlowe C.

    1986-01-01

    Reports how 10 college students in a psychology of personality course increased their awareness of their own implicit theories of personality by completing an 18-item questionnaire dealing with various issues of personality. Instructor-supplied responses representing the perspectives of Freud, Skinner, Rogers, and Glasser were used to foster…

  16. Finite difference schemes for long-time integration

    NASA Technical Reports Server (NTRS)

    Haras, Zigo; Taasan, Shlomo

    1993-01-01

    Finite difference schemes for the evaluation of first and second derivatives are presented. These second order compact schemes were designed for long-time integration of evolution equations by solving a quadratic constrained minimization problem. The quadratic cost function measures the global truncation error while taking into account the initial data. The resulting schemes are applicable for integration times fourfold, or more, longer than similar previously studied schemes. A similar approach was used to obtain improved integration schemes.

  17. Junction-Generalized Riemann Problem for stiff hyperbolic balance laws in networks: An implicit solver and ADER schemes

    NASA Astrophysics Data System (ADS)

    Contarino, Christian; Toro, Eleuterio F.; Montecinos, Gino I.; Borsche, Raul; Kall, Jochen

    2016-06-01

    In this paper we design a new implicit solver for the Junction-Generalized Riemann Problem (J-GRP), which is based on a recently proposed implicit method for solving the Generalized Riemann Problem (GRP) for systems of hyperbolic balance laws. We use the new J-GRP solver to construct an ADER scheme that is globally explicit, locally implicit and with no theoretical accuracy barrier, in both space and time. The resulting ADER scheme is able to deal with stiff source terms and can be applied to non-linear systems of hyperbolic balance laws in domains consisting on networks of one-dimensional sub-domains. In this paper we specifically apply the numerical techniques to networks of blood vessels. We report on a test problem with exact solution for a simplified network of three vessels meeting at a single junction, which is then used to carry out a systematic convergence rate study of the proposed high-order numerical methods. Schemes up to fifth order of accuracy in space and time are implemented and tested. We then show the ability of the ADER scheme to deal with stiff sources through a numerical simulation in a network of vessels. An application to a physical test problem consisting of a network of 37 compliant silicon tubes (arteries) and 21 junctions, reveals that it is imperative to use high-order methods at junctions, in order to preserve the desired high order of accuracy in the full computational domain. For example, it is demonstrated that a second-order method throughout, gives comparable results to a method that is fourth order in the interior of the domain and first order at junctions.

  18. Implicit and explicit schemes for mass consistency preservation in hybrid particle/finite-volume algorithms for turbulent reactive flows

    SciTech Connect

    Popov, Pavel P. Pope, Stephen B.

    2014-01-15

    This work addresses the issue of particle mass consistency in Large Eddy Simulation/Probability Density Function (LES/PDF) methods for turbulent reactive flows. Numerical schemes for the implicit and explicit enforcement of particle mass consistency (PMC) are introduced, and their performance is examined in a representative LES/PDF application, namely the Sandia–Sydney Bluff-Body flame HM1. A new combination of interpolation schemes for velocity and scalar fields is found to better satisfy PMC than multilinear and fourth-order Lagrangian interpolation. A second-order accurate time-stepping scheme for stochastic differential equations (SDE) is found to improve PMC relative to Euler time stepping, which is the first time that a second-order scheme is found to be beneficial, when compared to a first-order scheme, in an LES/PDF application. An explicit corrective velocity scheme for PMC enforcement is introduced, and its parameters optimized to enforce a specified PMC criterion with minimal corrective velocity magnitudes.

  19. On a fourth order accurate implicit finite difference scheme for hyperbolic conservation laws. I - Nonstiff strongly dynamic problems

    NASA Technical Reports Server (NTRS)

    Harten, A.; Tal-Ezer, H.

    1981-01-01

    An implicit finite difference method of fourth order accuracy in space and time is introduced for the numerical solution of one-dimensional systems of hyperbolic conservation laws. The basic form of the method is a two-level scheme which is unconditionally stable and nondissipative. The scheme uses only three mesh points at level t and three mesh points at level t + delta t. The dissipative version of the basic method given is conditionally stable under the CFL (Courant-Friedrichs-Lewy) condition. This version is particularly useful for the numerical solution of problems with strong but nonstiff dynamic features, where the CFL restriction is reasonable on accuracy grounds. Numerical results are provided to illustrate properties of the proposed method.

  20. Implicit integrations for SPH in semi-Lagrangian approach: Application to the accretion disc modeling in a microquasar

    NASA Astrophysics Data System (ADS)

    Lanzafame, G.

    2013-03-01

    Current explicit integration techniques in fluid dynamics are deeply limited by the Courant-Friedrichs-Lewy condition of the time step progression, based on the adopted spatial resolution coupled with the maximum value between the kinetic velocity or the signal transmission speed in the computational domain. Eulerian implicit integration techniques, even though more time consuming, can allow us to perform stable computational fluid dynamics paying the price of a relatively larger inaccuracy in the calculations, without suffering such a strict temporal limitation. In this paper, we present a simple and effective scheme to perform free Lagrangian Smooth Particle Hydrodynamics (SPH) implicit integrations in the semi-Lagrangian approach without any Jacobian matrix inversion operations for viscous Navier-Stokes flows. Applications to SPH accretion disc simulation around a massive black hole (MBH) in a binary stellar system are shown, together with the comparison to the same results obtained according to the traditional explicit integration techniques. Some 1D and 2D critical tests are also discussed to check the validity of the technique.

  1. Implicit approximate-factorization schemes for the efficient solution of steady transonic flow problems

    NASA Technical Reports Server (NTRS)

    Ballhaus, W. F.; Jameson, A.; Albert, J.

    1977-01-01

    Implicit approximate-factorization algorithms (AF) are developed for the solution of steady-state transonic flow problems. The performance of the AF solution method is evaluated relative to that of the standard solution method for transonic flow problems, successive line over-relaxation (SLOR). Both methods are applied to the solution of the nonlinear, two-dimensional transonic small-disturbance equation. Results indicate that the AF method requires substantially less computer time than SLOR to solve the nonlinear finite-difference matrix equation for a transonic flow field. This increase in computational efficiency is achieved with no appreciable increase in computer storage or coding complexity.

  2. The Importance of Using Explicit and then Implicit Schemes in the Fast time-scale Rupturing at Oceanic-Continental Boundary

    NASA Astrophysics Data System (ADS)

    So, B.; Yuen, D. A.; Lee, S.

    2011-12-01

    Numerical modeling in geodynamics, such as subduction and lithospheric rupture, normally uses only one scheme (e.g., implicit or explicit). However every geodynamical phenomena has a multi-time-scale instability. So these problems cannot be solved completely by just one numerical scheme since implicit and explicit schemes have different characteristics and stabilize in different time step size. That's why modelers should select an appropriate scheme for their problems. In our fully coupled thermal-mechanical finite element modeling for asymmetric instability initiation induced by shear modulus contrast between oceanic and continental lithospheres, two simply attached lithospheres with different shear modulus but having same visco-plastic rheologies are compressed with constant velocity of few centimeters per year. We use both explicit and implicit schemes at stages of elastic energy release and strengthening of shear localization, respectively. Since elastic energy quickly propagates from top and bottom of lithosphere (>1500 km/Myr), the employment of explicit scheme is more suitable to understand thermal runaway effect than implicit scheme at the initial stage of stored elastic energy release (i.e., this stage needs small time step size). The calculation results with only-implicit scheme and explicit-implicit hybrid scheme are different, because the latter can calculate fast time-scale energy dissipation and temperature field better than the former. The small temperature difference between implicit and explicit schemes may cause the large difference at later stage due to thermal-mechanical feedback. To investigate the timing of the initiation of asymmetric instability crossing the interface, temperature and plastic energy distributions are calculated on fine grid composed of shear modulus contrast and activation energy for a period of 1 Myr. We found that asymmetric shear instabilities are induced by a elastic shear modulus contrast over an wide range of

  3. Convergence of a class of semi-implicit time-stepping schemes for nonsmooth rigid multibody dynamics.

    SciTech Connect

    Gavrea, B. I.; Anitescu, M.; Potra, F. A.; Mathematics and Computer Science; Univ. of Pennsylvania; Univ. of Maryland

    2008-01-01

    In this work we present a framework for the convergence analysis in a measure differential inclusion sense of a class of time-stepping schemes for multibody dynamics with contacts, joints, and friction. This class of methods solves one linear complementarity problem per step and contains the semi-implicit Euler method, as well as trapezoidal-like methods for which second-order convergence was recently proved under certain conditions. By using the concept of a reduced friction cone, the analysis includes, for the first time, a convergence result for the case that includes joints. An unexpected intermediary result is that we are able to define a discrete velocity function of bounded variation, although the natural discrete velocity function produced by our algorithm may have unbounded variation.

  4. An implicit difference scheme for the long-time evolution of localized solutions of a generalized Boussinesq system

    SciTech Connect

    Christov, C.I.; Maugin, G.A.

    1995-01-01

    We consider the nonlinear system of equations built up from a generalized Boussinesq equation coupled with a wave equation which is a model for the one-dimensional dynamics of phases in martensitic alloys. The strongly implicit scheme employing Newton`s quasilinearisation allows us to track the long time evolution of the localized solutions of the system. Two distinct classes of solutions are encountered for the pure Boussinesq equation. The first class consists of oscillatory pulses whose envelopes are localized waves. The second class consists of smoother solutions whose shapes are either heteroclinic (kinks) or homoclinic (bumps). The homoclinics decrease in amplitude with time while their support increases. An appropriate self-similar scaling is found analytically and confirmed by the direct numerical simulations to high accuracy. The rich phenomenology resulting from the coupling with the wave equation is also investigated. 11 refs., 12 figs., 2 tabs.

  5. Integration of the shallow water equations on the sphere using a vector semi-Lagrangian scheme with a multigrid solver

    NASA Technical Reports Server (NTRS)

    Bates, J. R.; Semazzi, F. H. M.; Higgins, R. W.; Barros, Saulo R. M.

    1990-01-01

    A vector semi-Lagrangian semi-implicit two-time-level finite-difference integration scheme for the shallow water equations on the sphere is presented. A C-grid is used for the spatial differencing. The trajectory-centered discretization of the momentum equation in vector form eliminates pole problems and, at comparable cost, gives greater accuracy than a previous semi-Lagrangian finite-difference scheme which used a rotated spherical coordinate system. In terms of the insensitivity of the results to increasing timestep, the new scheme is as successful as recent spectral semi-Lagrangian schemes. In addition, the use of a multigrid method for solving the elliptic equation for the geopotential allows efficient integration with an operation count which, at high resolution, is of lower order than in the case of the spectral models. The properties of the new scheme should allow finite-difference models to compete with spectral models more effectively than has previously been possible.

  6. Study of shock-induced combustion using an implicit TVD scheme

    NASA Technical Reports Server (NTRS)

    Yungster, Shayne

    1992-01-01

    The supersonic combustion flowfields associated with various hypersonic propulsion systems, such as the ram accelerator, the oblique detonation wave engine, and the scramjet, are being investigated using a new computational fluid dynamics (CFD) code. The code solves the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. It employs an iterative method and a second order differencing scheme to improve computational efficiency. The code is currently being applied to study shock wave/boundary layer interactions in premixed combustible gases, and to investigate the ram accelerator concept. Results obtained for a ram accelerator configuration indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outward and downstream. The combustion process creates a high pressure region over the back of the projectile resulting in a net positive thrust forward.

  7. On the accuracy and convergence of implicit numerical integration of finite element generated ordinary differential equations

    NASA Technical Reports Server (NTRS)

    Baker, A. J.; Soliman, M. O.

    1978-01-01

    A study of accuracy and convergence of linear functional finite element solution to linear parabolic and hyperbolic partial differential equations is presented. A variable-implicit integration procedure is employed for the resultant system of ordinary differential equations. Accuracy and convergence is compared for the consistent and two lumped assembly procedures for the identified initial-value matrix structure. Truncation error estimation is accomplished using Richardson extrapolation.

  8. Implicit integration of the time-dependent Ginzburg-Landau equations of superconductivity.

    SciTech Connect

    Gunter, D. O.; Kaper, H. G.; Leaf, G. K.; Mathematics and Computer Science

    2002-05-17

    This article is concerned with the integration of the time-dependent Ginzburg-Landau (TDGL) equations of superconductivity. Four algorithms, ranging from fully explicit to nonlinearly implicit, are presented and evaluated for stability, accuracy, and compute time. The benchmark problem for the evaluation is the equilibration of a vortex configuration in a superconductor that is embedded in a thin insulator and subject to an applied magnetic field.

  9. Higher Order Time Integration Schemes for the Unsteady Navier-Stokes Equations on Unstructured Meshes

    NASA Technical Reports Server (NTRS)

    Jothiprasad, Giridhar; Mavriplis, Dimitri J.; Caughey, David A.; Bushnell, Dennis M. (Technical Monitor)

    2002-01-01

    The efficiency gains obtained using higher-order implicit Runge-Kutta schemes as compared with the second-order accurate backward difference schemes for the unsteady Navier-Stokes equations are investigated. Three different algorithms for solving the nonlinear system of equations arising at each timestep are presented. The first algorithm (NMG) is a pseudo-time-stepping scheme which employs a non-linear full approximation storage (FAS) agglomeration multigrid method to accelerate convergence. The other two algorithms are based on Inexact Newton's methods. The linear system arising at each Newton step is solved using iterative/Krylov techniques and left preconditioning is used to accelerate convergence of the linear solvers. One of the methods (LMG) uses Richardson's iterative scheme for solving the linear system at each Newton step while the other (PGMRES) uses the Generalized Minimal Residual method. Results demonstrating the relative superiority of these Newton's methods based schemes are presented. Efficiency gains as high as 10 are obtained by combining the higher-order time integration schemes with the more efficient nonlinear solvers.

  10. A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

    NASA Technical Reports Server (NTRS)

    Lopez Ortega, Alejandro; Mikellides, Ioannis G.

    2014-01-01

    We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

  11. Higher Order Time Integration Schemes for the Unsteady Navier-Stokes Equations on Unstructured Meshes

    NASA Technical Reports Server (NTRS)

    Jothiprasad, Giridhar; Mavriplis, Dimitri J.; Caughey, David A.

    2002-01-01

    The rapid increase in available computational power over the last decade has enabled higher resolution flow simulations and more widespread use of unstructured grid methods for complex geometries. While much of this effort has been focused on steady-state calculations in the aerodynamics community, the need to accurately predict off-design conditions, which may involve substantial amounts of flow separation, points to the need to efficiently simulate unsteady flow fields. Accurate unsteady flow simulations can easily require several orders of magnitude more computational effort than a corresponding steady-state simulation. For this reason, techniques for improving the efficiency of unsteady flow simulations are required in order to make such calculations feasible in the foreseeable future. The purpose of this work is to investigate possible reductions in computer time due to the choice of an efficient time-integration scheme from a series of schemes differing in the order of time-accuracy, and by the use of more efficient techniques to solve the nonlinear equations which arise while using implicit time-integration schemes. This investigation is carried out in the context of a two-dimensional unstructured mesh laminar Navier-Stokes solver.

  12. Analysis of implicit second-order upwind-biased stencils

    NASA Technical Reports Server (NTRS)

    Roberts, Thomas W.; Warren, Gary P.

    1993-01-01

    Truncation error and stability properties of several implicit upwind schemes for the two-dimensional Euler equations are examined. The schemes use linear data reconstruction methods to achieve second-order flux integrations where the implicit Jacobian operators are first order. The stability properties of the schemes are examined by a Von Neumann analysis of the linearized, constant-coefficient Euler equations. The choice of the data reconstruction method used to evaluate the flux integral has a dramatic effect on the convergence properties of the implicit solution method. In particular, the typical one-dimensional data reconstruction methods used with structured grids exhibit poor convergence properties compared to the unstructured grid method considered. Of the schemes examined, the one with the superior convergence properties is well-suited for both unstructured and structured grids, which has important implications for the design of implicit methods.

  13. Semi-implicit integration factor methods on sparse grids for high-dimensional systems

    NASA Astrophysics Data System (ADS)

    Wang, Dongyong; Chen, Weitao; Nie, Qing

    2015-07-01

    Numerical methods for partial differential equations in high-dimensional spaces are often limited by the curse of dimensionality. Though the sparse grid technique, based on a one-dimensional hierarchical basis through tensor products, is popular for handling challenges such as those associated with spatial discretization, the stability conditions on time step size due to temporal discretization, such as those associated with high-order derivatives in space and stiff reactions, remain. Here, we incorporate the sparse grids with the implicit integration factor method (IIF) that is advantageous in terms of stability conditions for systems containing stiff reactions and diffusions. We combine IIF, in which the reaction is treated implicitly and the diffusion is treated explicitly and exactly, with various sparse grid techniques based on the finite element and finite difference methods and a multi-level combination approach. The overall method is found to be efficient in terms of both storage and computational time for solving a wide range of PDEs in high dimensions. In particular, the IIF with the sparse grid combination technique is flexible and effective in solving systems that may include cross-derivatives and non-constant diffusion coefficients. Extensive numerical simulations in both linear and nonlinear systems in high dimensions, along with applications of diffusive logistic equations and Fokker-Planck equations, demonstrate the accuracy, efficiency, and robustness of the new methods, indicating potential broad applications of the sparse grid-based integration factor method.

  14. Semi-implicit Integration Factor Methods on Sparse Grids for High-Dimensional Systems

    PubMed Central

    Wang, Dongyong; Chen, Weitao; Nie, Qing

    2015-01-01

    Numerical methods for partial differential equations in high-dimensional spaces are often limited by the curse of dimensionality. Though the sparse grid technique, based on a one-dimensional hierarchical basis through tensor products, is popular for handling challenges such as those associated with spatial discretization, the stability conditions on time step size due to temporal discretization, such as those associated with high-order derivatives in space and stiff reactions, remain. Here, we incorporate the sparse grids with the implicit integration factor method (IIF) that is advantageous in terms of stability conditions for systems containing stiff reactions and diffusions. We combine IIF, in which the reaction is treated implicitly and the diffusion is treated explicitly and exactly, with various sparse grid techniques based on the finite element and finite difference methods and a multi-level combination approach. The overall method is found to be efficient in terms of both storage and computational time for solving a wide range of PDEs in high dimensions. In particular, the IIF with the sparse grid combination technique is flexible and effective in solving systems that may include cross-derivatives and non-constant diffusion coefficients. Extensive numerical simulations in both linear and nonlinear systems in high dimensions, along with applications of diffusive logistic equations and Fokker-Planck equations, demonstrate the accuracy, efficiency, and robustness of the new methods, indicating potential broad applications of the sparse grid-based integration factor method. PMID:25897178

  15. Time Integration Schemes for the Unsteady Navier-stokes Equations

    NASA Technical Reports Server (NTRS)

    Bijl, Hester; Carpenter, Mark H.; Vatsa, Veer N.

    2001-01-01

    The efficiency and accuracy of several time integration schemes are investigated for the unsteady Navier-Stokes equations. This study focuses on the efficiency of higher-order Runge-Kutta schemes in comparison with the popular Backward Differencing Formulations. For this comparison an unsteady two-dimensional laminar flow problem is chosen, i.e., flow around a circular cylinder at Re = 1200. It is concluded that for realistic error tolerances (smaller than 10(exp -1)) fourth-and fifth-order Runge-Kutta schemes are the most efficient. For reasons of robustness and computer storage, the fourth-order Runge-Kutta method is recommended. The efficiency of the fourth-order Runge-Kutta scheme exceeds that of second-order Backward Difference Formula by a factor of 2.5 at engineering error tolerance levels (10(exp -1) to 10(exp -2)). Efficiency gains are more dramatic at smaller tolerances.

  16. Enabling an Integrated Rate-temporal Learning Scheme on Memristor

    NASA Astrophysics Data System (ADS)

    He, Wei; Huang, Kejie; Ning, Ning; Ramanathan, Kiruthika; Li, Guoqi; Jiang, Yu; Sze, Jiayin; Shi, Luping; Zhao, Rong; Pei, Jing

    2014-04-01

    Learning scheme is the key to the utilization of spike-based computation and the emulation of neural/synaptic behaviors toward realization of cognition. The biological observations reveal an integrated spike time- and spike rate-dependent plasticity as a function of presynaptic firing frequency. However, this integrated rate-temporal learning scheme has not been realized on any nano devices. In this paper, such scheme is successfully demonstrated on a memristor. Great robustness against the spiking rate fluctuation is achieved by waveform engineering with the aid of good analog properties exhibited by the iron oxide-based memristor. The spike-time-dependence plasticity (STDP) occurs at moderate presynaptic firing frequencies and spike-rate-dependence plasticity (SRDP) dominates other regions. This demonstration provides a novel approach in neural coding implementation, which facilitates the development of bio-inspired computing systems.

  17. Enabling an Integrated Rate-temporal Learning Scheme on Memristor

    PubMed Central

    He, Wei; Huang, Kejie; Ning, Ning; Ramanathan, Kiruthika; Li, Guoqi; Jiang, Yu; Sze, JiaYin; Shi, Luping; Zhao, Rong; Pei, Jing

    2014-01-01

    Learning scheme is the key to the utilization of spike-based computation and the emulation of neural/synaptic behaviors toward realization of cognition. The biological observations reveal an integrated spike time- and spike rate-dependent plasticity as a function of presynaptic firing frequency. However, this integrated rate-temporal learning scheme has not been realized on any nano devices. In this paper, such scheme is successfully demonstrated on a memristor. Great robustness against the spiking rate fluctuation is achieved by waveform engineering with the aid of good analog properties exhibited by the iron oxide-based memristor. The spike-time-dependence plasticity (STDP) occurs at moderate presynaptic firing frequencies and spike-rate-dependence plasticity (SRDP) dominates other regions. This demonstration provides a novel approach in neural coding implementation, which facilitates the development of bio-inspired computing systems. PMID:24755608

  18. Integrated optical 3D digital imaging based on DSP scheme

    NASA Astrophysics Data System (ADS)

    Wang, Xiaodong; Peng, Xiang; Gao, Bruce Z.

    2008-03-01

    We present a scheme of integrated optical 3-D digital imaging (IO3DI) based on digital signal processor (DSP), which can acquire range images independently without PC support. This scheme is based on a parallel hardware structure with aid of DSP and field programmable gate array (FPGA) to realize 3-D imaging. In this integrated scheme of 3-D imaging, the phase measurement profilometry is adopted. To realize the pipeline processing of the fringe projection, image acquisition and fringe pattern analysis, we present a multi-threads application program that is developed under the environment of DSP/BIOS RTOS (real-time operating system). Since RTOS provides a preemptive kernel and powerful configuration tool, with which we are able to achieve a real-time scheduling and synchronization. To accelerate automatic fringe analysis and phase unwrapping, we make use of the technique of software optimization. The proposed scheme can reach a performance of 39.5 f/s (frames per second), so it may well fit into real-time fringe-pattern analysis and can implement fast 3-D imaging. Experiment results are also presented to show the validity of proposed scheme.

  19. Path Integrals, BRST Identities, and Regularization Schemes in Nonstandard Gauges

    NASA Astrophysics Data System (ADS)

    Ren, Hai-cang

    2000-07-01

    The path integral of a gauge theory is studied in Coulomb-like gauges. The Christ-Lee terms of operator ordering are reproduced within the path integration framework. In the presence of fermions, a new operator term, in addition to that of Christ and Lee, is discovered. Such terms are found to be instrumental in restoring the invariance of the effective Lagrangian under a field-dependent gauge transformation, which underlies the BRST symmetry. A unitary regularization scheme which maintains manifest BRST symmetry and is free from energy divergences is proposed for a nonabelian gauge field.

  20. An Efficient, Semi-implicit Pressure-based Scheme Employing a High-resolution Finitie Element Method for Simulating Transient and Steady, Inviscid and Viscous, Compressible Flows on Unstructured Grids

    SciTech Connect

    Richard C. Martineau; Ray A. Berry

    2003-04-01

    A new semi-implicit pressure-based Computational Fluid Dynamics (CFD) scheme for simulating a wide range of transient and steady, inviscid and viscous compressible flow on unstructured finite elements is presented here. This new CFD scheme, termed the PCICEFEM (Pressure-Corrected ICE-Finite Element Method) scheme, is composed of three computational phases, an explicit predictor, an elliptic pressure Poisson solution, and a semiimplicit pressure-correction of the flow variables. The PCICE-FEM scheme is capable of second-order temporal accuracy by incorporating a combination of a time-weighted form of the two-step Taylor-Galerkin Finite Element Method scheme as an explicit predictor for the balance of momentum equations and the finite element form of a time-weighted trapezoid rule method for the semi-implicit form of the governing hydrodynamic equations. Second-order spatial accuracy is accomplished by linear unstructured finite element discretization. The PCICE-FEM scheme employs Flux-Corrected Transport as a high-resolution filter for shock capturing. The scheme is capable of simulating flows from the nearly incompressible to the high supersonic flow regimes. The PCICE-FEM scheme represents an advancement in mass-momentum coupled, pressurebased schemes. The governing hydrodynamic equations for this scheme are the conservative form of the balance of momentum equations (Navier-Stokes), mass conservation equation, and total energy equation. An operator splitting process is performed along explicit and implicit operators of the semi-implicit governing equations to render the PCICE-FEM scheme in the class of predictor-corrector schemes. The complete set of semi-implicit governing equations in the PCICE-FEM scheme are cast in this form, an explicit predictor phase and a semi-implicit pressure-correction phase with the elliptic pressure Poisson solution coupling the predictor-corrector phases. The result of this predictor-corrector formulation is that the pressure Poisson

  1. Fully Integrated EAS-Based Solid-Shell Finite Elements in Implicit Sheet Metal Forming Simulations

    SciTech Connect

    Valente, R.A. Fontes; Cardoso, R.P.R.; Alves de Sousa, R.J.; Parente, M.P.L.; Jorge, R.M. Natal

    2005-08-05

    In this communication sheet metal forming problems are analyzed with the Finite Element Method and a fully-integrated solid-shell element, based on the Enhanced Assumed Strain (EAS) method. Among the solid-shell element's distinguish features, it should be mentioned the solely use of the EAS approach in dealing with either transverse and volumetric-based locking pathologies, thus avoiding the inclusion of other mixed methods into the formulation. The adopted methodology is then able to successfully deal with small thickness shell problems within the incompressible range, aspects commonly appearing in sheet metal forming modeling with solid elements.Simulations of this type of forming processes are mainly solved resorting to membrane and shell-type finite elements, included in explicit commercial programs. Nevertheless, the presented solid-shell formulation, within a fully implicit approach, provides reliable solutions when compared to experimental results. It is also worth mentioning that the present solid-shell formulation encompasses a minimum set of enhancing strain variables, if compared to other fully integrated hexahedral finite elements in the literature.In order to assess the performance of the presented formulation, the S-Rail Forming problem of an aluminum alloy is described and analyzed, with the results being compared to experimental and numerical simulation data.

  2. A lightweight data integrity scheme for sensor networks.

    PubMed

    Kamel, Ibrahim; Juma, Hussam

    2011-01-01

    Limited energy is the most critical constraint that limits the capabilities of wireless sensor networks (WSNs). Most sensors operate on batteries with limited power. Battery recharging or replacement may be impossible. Security mechanisms that are based on public key cryptographic algorithms such as RSA and digital signatures are prohibitively expensive in terms of energy consumption and storage requirements, and thus unsuitable for WSN applications. This paper proposes a new fragile watermarking technique to detect unauthorized alterations in WSN data streams. We propose the FWC-D scheme, which uses group delimiters to keep the sender and receivers synchronized and help them to avoid ambiguity in the event of data insertion or deletion. The watermark, which is computed using a hash function, is stored in the previous group in a linked-list fashion to ensure data freshness and mitigate replay attacks, FWC-D generates a serial number SN that is attached to each group to help the receiver determines how many group insertions or deletions occurred. Detailed security analysis that compares the proposed FWC-D scheme with SGW, one of the latest integrity schemes for WSNs, shows that FWC-D is more robust than SGW. Simulation results further show that the proposed scheme is much faster than SGW. PMID:22163840

  3. A Lightweight Data Integrity Scheme for Sensor Networks

    PubMed Central

    Kamel, Ibrahim; Juma, Hussam

    2011-01-01

    Limited energy is the most critical constraint that limits the capabilities of wireless sensor networks (WSNs). Most sensors operate on batteries with limited power. Battery recharging or replacement may be impossible. Security mechanisms that are based on public key cryptographic algorithms such as RSA and digital signatures are prohibitively expensive in terms of energy consumption and storage requirements, and thus unsuitable for WSN applications. This paper proposes a new fragile watermarking technique to detect unauthorized alterations in WSN data streams. We propose the FWC-D scheme, which uses group delimiters to keep the sender and receivers synchronized and help them to avoid ambiguity in the event of data insertion or deletion. The watermark, which is computed using a hash function, is stored in the previous group in a linked-list fashion to ensure data freshness and mitigate replay attacks, FWC-D generates a serial number SN that is attached to each group to help the receiver determines how many group insertions or deletions occurred. Detailed security analysis that compares the proposed FWC-D scheme with SGW, one of the latest integrity schemes for WSNs, shows that FWC-D is more robust than SGW. Simulation results further show that the proposed scheme is much faster than SGW. PMID:22163840

  4. TRHD: Three-temperature radiation-hydrodynamics code with an implicit non-equilibrium radiation transport using a cell-centered monotonic finite volume scheme on unstructured-grids

    NASA Astrophysics Data System (ADS)

    Sijoy, C. D.; Chaturvedi, S.

    2015-05-01

    Three-temperature (3T), unstructured-mesh, non-equilibrium radiation hydrodynamics (RHD) code have been developed for the simulation of intense thermal radiation or high-power laser driven radiative shock hydrodynamics in two-dimensional (2D) axis-symmetric geometries. The governing hydrodynamics equations are solved using a compatible unstructured Lagrangian method based on a control volume differencing (CVD) scheme. A second-order predictor-corrector (PC) integration scheme is used for the temporal discretization of the hydrodynamics equations. For the radiation energy transport, frequency averaged gray model is used in which the flux-limited diffusion (FLD) approximation is used to recover the free-streaming limit of the radiation propagation in optically thin regions. The proposed RHD model allows to have different temperatures for the electrons and ions. In addition to this, the electron and thermal radiation temperatures are assumed to be in non-equilibrium. Therefore, the thermal relaxation between the electrons and ions and the coupling between the radiation and matter energies are required to be computed self-consistently. For this, the coupled flux limited electron heat conduction and the non-equilibrium radiation diffusion equations are solved simultaneously by using an implicit, axis-symmetric, cell-centered, monotonic, nonlinear finite volume (NLFV) scheme. In this paper, we have described the details of the 2D, 3T, non-equilibrium RHD code developed along with a suite of validation test problems to demonstrate the accuracy and performance of the algorithms. We have also conducted a performance analysis with different linearity preserving interpolation schemes that are used for the evaluation of the nodal values in the NLFV scheme. Finally, in order to demonstrate full capability of the code implementation, we have presented the simulation of laser driven thin Aluminum (Al) foil acceleration. The simulation results are found to be in good agreement

  5. Efficient implicit integration for finite-strain viscoplasticity with a nested multiplicative split

    NASA Astrophysics Data System (ADS)

    Shutov, A. V.

    2016-07-01

    An efficient and reliable stress computation algorithm is presented, which is based on implicit integration of the local evolution equations of multiplicative finite-strain plasticity/viscoplasticity. The algorithm is illustrated by an example involving a combined nonlinear isotropic/kinematic hardening; numerous backstress tensors are employed for a better description of the material behavior. The considered material model exhibits the so-called weak invariance under arbitrary isochoric changes of the reference configuration, and the presented algorithm retains this useful property. Even more: the weak invariance serves as a guide in constructing this algorithm. The constraint of inelastic incompressibility is exactly preserved as well. The proposed method is first-order accurate. Concerning the accuracy of the stress computation, the new algorithm is comparable to the Euler Backward method with a subsequent correction of incompressibility (EBMSC) and the classical exponential method (EM). Regarding the computational efficiency, the new algorithm is superior to the EBMSC and EM. Some accuracy tests are presented using parameters of the aluminum alloy 5754-O and the 42CrMo4 steel. FEM solutions of two boundary value problems using MSC.MARC are presented to show the correctness of the numerical implementation.

  6. How to Improve Adolescent Stress Responses: Insights From Integrating Implicit Theories of Personality and Biopsychosocial Models.

    PubMed

    Yeager, David S; Lee, Hae Yeon; Jamieson, Jeremy P

    2016-08-01

    This research integrated implicit theories of personality and the biopsychosocial model of challenge and threat, hypothesizing that adolescents would be more likely to conclude that they can meet the demands of an evaluative social situation when they were taught that people have the potential to change their socially relevant traits. In Study 1 (N = 60), high school students were assigned to an incremental-theory-of-personality or a control condition and then given a social-stress task. Relative to control participants, incremental-theory participants exhibited improved stress appraisals, more adaptive neuroendocrine and cardiovascular responses, and better performance outcomes. In Study 2 (N = 205), we used a daily-diary intervention to test high school students' stress reactivity outside the laboratory. Threat appraisals (Days 5-9 after intervention) and neuroendocrine responses (Days 8 and 9 after intervention only) were unrelated to the intensity of daily stressors when adolescents received the incremental-theory intervention. Students who received the intervention also had better grades over freshman year than those who did not. These findings offer new avenues for improving theories of adolescent stress and coping. PMID:27324267

  7. An implicit scheme for solving the anisotropic diffusion of heat and cosmic rays in the RAMSES code

    NASA Astrophysics Data System (ADS)

    Dubois, Yohan; Commerçon, Benoît

    2016-01-01

    Astrophysical plasmas are subject to a tight connection between magnetic fields and the diffusion of particles, which leads to an anisotropic transport of energy. Under the fluid assumption, this effect can be reduced to an advection-diffusion equation, thereby augmenting the equations of magnetohydrodynamics. We introduce a new method for solving the anisotropic diffusion equation using an implicit finite-volume method with adaptive mesh refinement and adaptive time-stepping in the ramses code. We apply this numerical solver to the diffusion of cosmic ray energy and diffusion of heat carried by electrons, which couple to the ion temperature. We test this new implementation against several numerical experiments and apply it to a simple supernova explosion with a uniform magnetic field.

  8. Implicit plasma simulation

    SciTech Connect

    Langdon, A.B.

    1985-03-03

    Implicit time integration methods have been used extensively in numerical modelling of slowly varying phenomena in systems that also support rapid variation. Examples include diffusion, hydrodynamics and reaction kinetics. This article discussed implementation of implicit time integration in plasma codes of the ''particle-in-cell'' family, and the benefits to be gained.

  9. An old integration scheme for compressible flows revisited, refurbished and put to work

    NASA Technical Reports Server (NTRS)

    Moretti, G.

    1978-01-01

    A scheme for integrating the Euler equations of compressible flow in any hyperbolic case is presented. The scheme relies on the concept of characteristics but is strictly a finite difference scheme. Improvements in accuracy and physical consistence due to the scheme are discussed and results of its application to complex flows are shown.

  10. A globally convergent matrix-free algorithm for implicit time-marching schemes arising in finite element analysis in fluids

    NASA Technical Reports Server (NTRS)

    Johan, Zdenek; Hughes, Thomas J. R.; Shakib, Farzin

    1991-01-01

    A solution procedure for solving nonlinear time-marching problems is presented. The nonsymmetric systems of equations arising from a Newton-type linearization of these time-marching problems are solved using an iterative strategy based on the generalized minimal residual (GMRES) algorithm. Matrix-free techniques leading to reduction in storage are presented. Incorporation of a linesearch algorithm in the Newton-GMRES scheme is discussed. An automatic time-increment control strategy is developed to increase the stability of the time-marching process. High-speed flow computations demonstrate the effectiveness of these algorithms.

  11. Evaluation of premeability-porosity relationships linked to mineral dissolution-precipitation using global implicit approach with a reduction scheme and operator splitting approach

    NASA Astrophysics Data System (ADS)

    Zolfaghari, R.; Shao, H.; Kolditz, O.

    2013-12-01

    Numerical simulation of reactive transport processes is essential in long term behavior assessment of hazardous materials. To simulate reactive transport processes global implicit approach (GIA) and operator splitting approach are commonly used. GIA has been getting more attentions due to advances in computational power and the lack of numerical accuracy and efficacy of operator splitting methods for simulating long term processes over the past few years. We have investigated the Efficiency and accuracy of these methods in handling slow reacting-processes in long term scenarios. GIA with reduction scheme proposed by Kräutel et al. (2010) and sequential non-iterative approach (SNIA) approach have been implemented into OpenGeoSys (OGS6) to solve reactive transport problems. The new reduction scheme in GIA uses a reformulation to reduce the number of coupled nonlinear partial differential equations by decoupling of equations and elimination of unknowns. The new reformulation divides components and species of the chemical system into decoupled linear reaction invariant components and coupled nonlinear reaction variant ones. A local chemical solver is used to handle the chemical problem in GIA and SNIA approaches. Equilibrium/ kinetic mineral reaction is treated as a complementarity problem in the local problem. In this context, a series of benchmarks have been adopted to assess the performance of GIA with reduction scheme and SNIA. The benchmarks objective is to simulate mineral dissolution-precipitation induced porosity changes and the resulting effects on the solute migration. The Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. The results produced by three codes of OGS6, OGS-PHREEQC and MIN3P have been compared and evaluated based on the benchmarks for the numerical accuracy and efficacy.

  12. Implicit Kalman filtering

    NASA Technical Reports Server (NTRS)

    Skliar, M.; Ramirez, W. F.

    1997-01-01

    For an implicitly defined discrete system, a new algorithm for Kalman filtering is developed and an efficient numerical implementation scheme is proposed. Unlike the traditional explicit approach, the implicit filter can be readily applied to ill-conditioned systems and allows for generalization to descriptor systems. The implementation of the implicit filter depends on the solution of the congruence matrix equation (A1)(Px)(AT1) = Py. We develop a general iterative method for the solution of this equation, and prove necessary and sufficient conditions for convergence. It is shown that when the system matrices of an implicit system are sparse, the implicit Kalman filter requires significantly less computer time and storage to implement as compared to the traditional explicit Kalman filter. Simulation results are presented to illustrate and substantiate the theoretical developments.

  13. Satellite attitude dynamics and estimation with the implicit midpoint method

    NASA Astrophysics Data System (ADS)

    Hellström, Christian; Mikkola, Seppo

    2009-07-01

    We describe the application of the implicit midpoint integrator to the problem of attitude dynamics for low-altitude satellites without the use of quaternions. Initially, we consider the satellite to rotate without external torques applied to it. We compare the numerical solution with the exact solution in terms of Jacobi's elliptic functions. Then, we include the gravity-gradient torque, where the implicit midpoint integrator proves to be a fast, simple and accurate method. Higher-order versions of the implicit midpoint scheme are compared to Gauss-Legendre Runge-Kutta methods in terms of accuracy and processing time. Finally, we investigate the performance of a parameter-adaptive Kalman filter based on the implicit midpoint integrator for the determination of the principal moments of inertia through observations.

  14. An implicit fast Fourier transform method for integration of the time dependent Schrodinger equation

    SciTech Connect

    Riley, M.E.; Ritchie, A.B.

    1997-12-31

    One finds that the conventional exponentiated split operator procedure is subject to difficulties when solving the time-dependent Schrodinger equation for Coulombic systems. By rearranging the kinetic and potential energy terms in the temporal propagator of the finite difference equations, one can find a propagation algorithm for three dimensions that looks much like the Crank-Nicholson and alternating direction implicit methods for one- and two-space-dimensional partial differential equations. The authors report investigations of this novel implicit split operator procedure. The results look promising for a purely numerical approach to certain electron quantum mechanical problems. A charge exchange calculation is presented as an example of the power of the method.

  15. Semi-implicit time-integrators for a scalable spectral element atmospheric model

    NASA Astrophysics Data System (ADS)

    Giraldo, Francis X.

    2005-07-01

    The Naval Research Laboratory's spectral element atmospheric model (NSEAM) for scalable computer architectures is presented. This new dynamical core is based on a high-order spectral element (SE) method in space and uses semi-implicit methods in time based on either the traditional second-order leapfrog (LF2) or second-order backward difference formulas (BDF2). The novelties of NSEAM are: it is geometrically flexible and thereby can accommodate any type of grid; LF2 or BDF2 are used to construct the semi-implicit method; and the horizontal operators are written, discretized, and solved in three-dimensional Cartesian space. The semi-implicit NSEAM is validated using: five baroclinic test cases; direct comparisons to the explicit version of NSEAM which has been extensively tested and the results previously reported in the literature; and comparisons with operational weather prediction and well-established climate models. A comparison with the US Navy's spectral transform global forecast model illustrates that NSEAM is 60% faster on an IBM SP4 using 96 processors for the current operational resolution of T239 L30. However, NSEAM can accommodate many more processors while continuing to scale efficiently even at higher grid resolutions. In fact, we show that at T498 L60, NSEAM scales linearly up to 384 processors.

  16. Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

    NASA Technical Reports Server (NTRS)

    Saleeb, Atef F.; Li, Wei

    1995-01-01

    This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the

  17. Marginalization: A Revisitation With Integration of Scholarship on Globalization, Intersectionality, Privilege, Microaggressions, and Implicit Biases.

    PubMed

    Hall, Joanne M; Carlson, Kelly

    2016-01-01

    In 1994, the concept of marginalization was explored in an article in Advances in Nursing Science. This is a revisitation of the concept incorporating new scholarship. This update is founded on feminism, postcolonialism, critical race theory, and discourse deconstruction, all viewpoints that have been explicated in nursing. The purpose of this analysis is to look at new scholarship and concepts useful to applying marginalization in nursing knowledge development from the standpoint of Bourdieu's macro, meso, and micro levels. New scholarship includes globalization, intersectionality, privilege, microaggressions, and implicit bias. Implications for decreasing health disparities through this new scholarship are discussed. PMID:27490876

  18. Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

    NASA Technical Reports Server (NTRS)

    Li, Wei; Saleeb, Atef F.

    1995-01-01

    This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of

  19. Implicit time-integration method for simultaneous solution of a coupled non-linear system

    NASA Astrophysics Data System (ADS)

    Watson, Justin Kyle

    Historically large physical problems have been divided into smaller problems based on the physics involved. This is no different in reactor safety analysis. The problem of analyzing a nuclear reactor for design basis accidents is performed by a handful of computer codes each solving a portion of the problem. The reactor thermal hydraulic response to an event is determined using a system code like TRAC RELAP Advanced Computational Engine (TRACE). The core power response to the same accident scenario is determined using a core physics code like Purdue Advanced Core Simulator (PARCS). Containment response to the reactor depressurization in a Loss Of Coolant Accident (LOCA) type event is calculated by a separate code. Sub-channel analysis is performed with yet another computer code. This is just a sample of the computer codes used to solve the overall problems of nuclear reactor design basis accidents. Traditionally each of these codes operates independently from each other using only the global results from one calculation as boundary conditions to another. Industry's drive to uprate power for reactors has motivated analysts to move from a conservative approach to design basis accident towards a best estimate method. To achieve a best estimate calculation efforts have been aimed at coupling the individual physics models to improve the accuracy of the analysis and reduce margins. The current coupling techniques are sequential in nature. During a calculation time-step data is passed between the two codes. The individual codes solve their portion of the calculation and converge to a solution before the calculation is allowed to proceed to the next time-step. This thesis presents a fully implicit method of simultaneous solving the neutron balance equations, heat conduction equations and the constitutive fluid dynamics equations. It discusses the problems involved in coupling different physics phenomena within multi-physics codes and presents a solution to these problems

  20. A more secure anonymous user authentication scheme for the integrated EPR information system.

    PubMed

    Wen, Fengtong

    2014-05-01

    Secure and efficient user mutual authentication is an essential task for integrated electronic patient record (EPR) information system. Recently, several authentication schemes have been proposed to meet this requirement. In a recent paper, Lee et al. proposed an efficient and secure password-based authentication scheme used smart cards for the integrated EPR information system. This scheme is believed to have many abilities to resist a range of network attacks. Especially, they claimed that their scheme could resist lost smart card attack. However, we reanalyze the security of Lee et al.'s scheme, and show that it fails to protect off-line password guessing attack if the secret information stored in the smart card is compromised. This also renders that their scheme is insecure against user impersonation attacks. Then, we propose a new user authentication scheme for integrated EPR information systems based on the quadratic residues. The new scheme not only resists a range of network attacks but also provides user anonymity. We show that our proposed scheme can provide stronger security. PMID:24760224

  1. Advanced receiver autonomous integrity monitoring (ARAIM) schemes with GNSS time offsets

    NASA Astrophysics Data System (ADS)

    Wu, Yun; Wang, Jinling; Jiang, Yiping

    2013-07-01

    Within the current Advanced Receiver Integrity Monitoring (ARAIM) scheme, the time offsets between different Global Navigation Satellite System (GNSS) constellations are estimated along with a position solution and the GNSS receiver clock error. This scheme is called the Time-offsets Estimated ARAIM, or the TOE ARAIM. In order to enhance the interoperability and compatibility between different constellations, the time offsets are expected to be broadcast to users in future multi-GNSS positioning and navigation applications. This paper describes two new ARAIM schemes to make use of the Broadcast Time Offsets (BTOs): Time Offset Observed ARAIM (TOO ARAIM) and Time Offset Synchronized ARAIM (TOS ARAIM). It has been shown that the VPL performances of these two new ARAIM schemes rely strongly on the accuracy of BTOs. By varying the error model of the BTOs, the simulation results also demonstrate that the proposed new TOO ARAIM scheme can outperform the existing TOE ARAIM scheme-even if the accuracy of BTOs for integrity is degraded to 4.5 m and the probability of a BTO fault is relaxed to 10-2 h-1. In addition, the new Time Offset Synchronized ARAIM scheme (TOS ARAIM) can also perform better than the existing TOE ARAIM scheme if the accuracy of BTOs for integrity can reach 0.75 m. As the TOO ARAIM also has a very relaxed requirement on BTOs and better VPL performance, the TOO ARAIM is regarded as a superior ARAIM scheme for multi-GNSS with BTOs available.

  2. Vehicle Scheduling Schemes for Commercial and Emergency Logistics Integration

    PubMed Central

    Li, Xiaohui; Tan, Qingmei

    2013-01-01

    In modern logistics operations, large-scale logistics companies, besides active participation in profit-seeking commercial business, also play an essential role during an emergency relief process by dispatching urgently-required materials to disaster-affected areas. Therefore, an issue has been widely addressed by logistics practitioners and caught researchers' more attention as to how the logistics companies achieve maximum commercial profit on condition that emergency tasks are effectively and performed satisfactorily. In this paper, two vehicle scheduling models are proposed to solve the problem. One is a prediction-related scheme, which predicts the amounts of disaster-relief materials and commercial business and then accepts the business that will generate maximum profits; the other is a priority-directed scheme, which, firstly groups commercial and emergency business according to priority grades and then schedules both types of business jointly and simultaneously by arriving at the maximum priority in total. Moreover, computer-based simulations are carried out to evaluate the performance of these two models by comparing them with two traditional disaster-relief tactics in China. The results testify the feasibility and effectiveness of the proposed models. PMID:24391724

  3. Vehicle scheduling schemes for commercial and emergency logistics integration.

    PubMed

    Li, Xiaohui; Tan, Qingmei

    2013-01-01

    In modern logistics operations, large-scale logistics companies, besides active participation in profit-seeking commercial business, also play an essential role during an emergency relief process by dispatching urgently-required materials to disaster-affected areas. Therefore, an issue has been widely addressed by logistics practitioners and caught researchers' more attention as to how the logistics companies achieve maximum commercial profit on condition that emergency tasks are effectively and performed satisfactorily. In this paper, two vehicle scheduling models are proposed to solve the problem. One is a prediction-related scheme, which predicts the amounts of disaster-relief materials and commercial business and then accepts the business that will generate maximum profits; the other is a priority-directed scheme, which, firstly groups commercial and emergency business according to priority grades and then schedules both types of business jointly and simultaneously by arriving at the maximum priority in total. Moreover, computer-based simulations are carried out to evaluate the performance of these two models by comparing them with two traditional disaster-relief tactics in China. The results testify the feasibility and effectiveness of the proposed models. PMID:24391724

  4. Multiscale integration schemes for jump-diffusion systems

    SciTech Connect

    Givon, D.; Kevrekidis, I.G.

    2008-12-09

    We study a two-time-scale system of jump-diffusion stochastic differential equations. We analyze a class of multiscale integration methods for these systems, which, in the spirit of [1], consist of a hybridization between a standard solver for the slow components and short runs for the fast dynamics, which are used to estimate the effect that the fast components have on the slow ones. We obtain explicit bounds for the discrepancy between the results of the multiscale integration method and the slow components of the original system.

  5. Implicit CAPTCHAs

    NASA Astrophysics Data System (ADS)

    Baird, Henry S.; Bentley, Jon L.

    2004-12-01

    We propose a design methodology for "implicit" CAPTCHAs to relieve drawbacks of present technology. CAPTCHAs are tests administered automatically over networks that can distinguish between people and machines and thus protect web services from abuse by programs masquerading as human users. All existing CAPTCHAs' challenges require a significant conscious effort by the person answering them -- e.g. reading and typing a nonsense word -- whereas implicit CAPTCHAs may require as little as a single click. Many CAPTCHAs distract and interrupt users, since the challenge is perceived as an irrelevant intrusion; implicit CAPTCHAs can be woven into the expected sequence of browsing using cues tailored to the site. Most existing CAPTCHAs are vulnerable to "farming-out" attacks in which challenges are passed to a networked community of human readers; by contrast, implicit CAPTCHAs are not "fungible" (in the sense of easily answerable in isolation) since they are meaningful only in the specific context of the website that is protected. Many existing CAPTCHAs irritate or threaten users since they are obviously tests of skill: implicit CAPTCHAs appear to be elementary and inevitable acts of browsing. It can often be difficult to detect when CAPTCHAs are under attack: implicit CAPTCHAs can be designed so that certain failure modes are correlated with failed bot attacks. We illustrate these design principles with examples.

  6. Implicit CAPTCHAs

    NASA Astrophysics Data System (ADS)

    Baird, Henry S.; Bentley, Jon L.

    2005-01-01

    We propose a design methodology for "implicit" CAPTCHAs to relieve drawbacks of present technology. CAPTCHAs are tests administered automatically over networks that can distinguish between people and machines and thus protect web services from abuse by programs masquerading as human users. All existing CAPTCHAs' challenges require a significant conscious effort by the person answering them -- e.g. reading and typing a nonsense word -- whereas implicit CAPTCHAs may require as little as a single click. Many CAPTCHAs distract and interrupt users, since the challenge is perceived as an irrelevant intrusion; implicit CAPTCHAs can be woven into the expected sequence of browsing using cues tailored to the site. Most existing CAPTCHAs are vulnerable to "farming-out" attacks in which challenges are passed to a networked community of human readers; by contrast, implicit CAPTCHAs are not "fungible" (in the sense of easily answerable in isolation) since they are meaningful only in the specific context of the website that is protected. Many existing CAPTCHAs irritate or threaten users since they are obviously tests of skill: implicit CAPTCHAs appear to be elementary and inevitable acts of browsing. It can often be difficult to detect when CAPTCHAs are under attack: implicit CAPTCHAs can be designed so that certain failure modes are correlated with failed bot attacks. We illustrate these design principles with examples.

  7. A Novel Multilayered RFID Tagged Cargo Integrity Assurance Scheme.

    PubMed

    Yang, Ming Hour; Luo, Jia Ning; Lu, Shao Yong

    2015-01-01

    To minimize cargo theft during transport, mobile radio frequency identification (RFID) grouping proof methods are generally employed to ensure the integrity of entire cargo loads. However, conventional grouping proofs cannot simultaneously generate grouping proofs for a specific group of RFID tags. The most serious problem of these methods is that nonexistent tags are included in the grouping proofs because of the considerable amount of time it takes to scan a high number of tags. Thus, applying grouping proof methods in the current logistics industry is difficult. To solve this problem, this paper proposes a method for generating multilayered offline grouping proofs. The proposed method provides tag anonymity; moreover, resolving disputes between recipients and transporters over the integrity of cargo deliveries can be expedited by generating grouping proofs and automatically authenticating the consistency between the receipt proof and pick proof. The proposed method can also protect against replay attacks, multi-session attacks, and concurrency attacks. Finally, experimental results verify that, compared with other methods for generating grouping proofs, the proposed method can efficiently generate offline grouping proofs involving several parties in a supply chain using mobile RFID. PMID:26512673

  8. A Novel Multilayered RFID Tagged Cargo Integrity Assurance Scheme

    PubMed Central

    Yang, Ming Hour; Luo, Jia Ning; Lu, Shao Yong

    2015-01-01

    To minimize cargo theft during transport, mobile radio frequency identification (RFID) grouping proof methods are generally employed to ensure the integrity of entire cargo loads. However, conventional grouping proofs cannot simultaneously generate grouping proofs for a specific group of RFID tags. The most serious problem of these methods is that nonexistent tags are included in the grouping proofs because of the considerable amount of time it takes to scan a high number of tags. Thus, applying grouping proof methods in the current logistics industry is difficult. To solve this problem, this paper proposes a method for generating multilayered offline grouping proofs. The proposed method provides tag anonymity; moreover, resolving disputes between recipients and transporters over the integrity of cargo deliveries can be expedited by generating grouping proofs and automatically authenticating the consistency between the receipt proof and pick proof. The proposed method can also protect against replay attacks, multi-session attacks, and concurrency attacks. Finally, experimental results verify that, compared with other methods for generating grouping proofs, the proposed method can efficiently generate offline grouping proofs involving several parties in a supply chain using mobile RFID. PMID:26512673

  9. An integrable high resolution all-optical analog-to-digital conversion scheme

    NASA Astrophysics Data System (ADS)

    Wei, Shile; Jian, Wu; Zhao, Lingjuan; Lu, Dan; Qiu, Jifang

    2014-05-01

    A novel 4 × 4 multimode interference couplers based phase-shifted photonic quantization scheme using multiwavelength mode locked pulse lasers as sampling source for all-optical analog-to-digital converter is proposed. Numerical analysis indicates that 8-bit quantization resolution operating at 40 GHz bandwidth could be achieved with an incident average optical power of 1.932 mW to each photodiode. The whole scheme can be integrated on a InP-based chip.

  10. Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme

    SciTech Connect

    Squire, J.; Tang, W. M.; Qin, H.

    2012-08-15

    A fully variational, unstructured, electromagnetic particle-in-cell integrator is developed for integration of the Vlasov-Maxwell equations. Using the formalism of discrete exterior calculus [Desbrun et al., e-print arXiv:math/0508341 (2005)], the field solver, interpolation scheme, and particle advance algorithm are derived through minimization of a single discrete field theory action. As a consequence of ensuring that the action is invariant under discrete electromagnetic gauge transformations, the integrator exactly conserves Gauss's law.

  11. Geometric Integration Of The Vlasov-Maxwell System With A Variational Particle-in-cell Scheme

    SciTech Connect

    J. Squire, H. Qin and W.M. Tang

    2012-03-27

    A fully variational, unstructured, electromagnetic particle-in-cell integrator is developed for integration of the Vlasov-Maxwell equations. Using the formalism of Discrete Exterior Calculus [1], the field solver, interpolation scheme and particle advance algorithm are derived through minimization of a single discrete field theory action. As a consequence of ensuring that the action is invariant under discrete electromagnetic gauge transformations, the integrator exactly conserves Gauss's law.

  12. Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme

    NASA Astrophysics Data System (ADS)

    Squire, J.; Qin, H.; Tang, W. M.

    2012-08-01

    A fully variational, unstructured, electromagnetic particle-in-cell integrator is developed for integration of the Vlasov-Maxwell equations. Using the formalism of discrete exterior calculus [Desbrun et al., e-print arXiv:math/0508341 (2005)], the field solver, interpolation scheme, and particle advance algorithm are derived through minimization of a single discrete field theory action. As a consequence of ensuring that the action is invariant under discrete electromagnetic gauge transformations, the integrator exactly conserves Gauss's law.

  13. New parallelizable schemes for integrating the Dissipative Particle Dynamics with Energy conservation.

    PubMed

    Homman, Ahmed-Amine; Maillet, Jean-Bernard; Roussel, Julien; Stoltz, Gabriel

    2016-01-14

    This work presents new parallelizable numerical schemes for the integration of dissipative particle dynamics with energy conservation. So far, no numerical scheme introduced in the literature is able to correctly preserve the energy over long times and give rise to small errors on average properties for moderately small time steps, while being straightforwardly parallelizable. We present in this article two new methods, both straightforwardly parallelizable, allowing to correctly preserve the total energy of the system. We illustrate the accuracy and performance of these new schemes both on equilibrium and nonequilibrium parallel simulations. PMID:26772559

  14. New parallelizable schemes for integrating the Dissipative Particle Dynamics with Energy conservation

    NASA Astrophysics Data System (ADS)

    Homman, Ahmed-Amine; Maillet, Jean-Bernard; Roussel, Julien; Stoltz, Gabriel

    2016-01-01

    This work presents new parallelizable numerical schemes for the integration of dissipative particle dynamics with energy conservation. So far, no numerical scheme introduced in the literature is able to correctly preserve the energy over long times and give rise to small errors on average properties for moderately small time steps, while being straightforwardly parallelizable. We present in this article two new methods, both straightforwardly parallelizable, allowing to correctly preserve the total energy of the system. We illustrate the accuracy and performance of these new schemes both on equilibrium and nonequilibrium parallel simulations.

  15. Comparison of Implicit Schemes to Solve Equations of Radiation Hydrodynamics with a Flux-limited Diffusion Approximation: Newton--Raphson, Operator Splitting, and Linearization

    NASA Astrophysics Data System (ADS)

    Tetsu, Hiroyuki; Nakamoto, Taishi

    2016-03-01

    Radiation is an important process of energy transport, a force, and a basis for synthetic observations, so radiation hydrodynamics (RHD) calculations have occupied an important place in astrophysics. However, although the progress in computational technology is remarkable, their high numerical cost is still a persistent problem. In this work, we compare the following schemes used to solve the nonlinear simultaneous equations of an RHD algorithm with the flux-limited diffusion approximation: the Newton-Raphson (NR) method, operator splitting, and linearization (LIN), from the perspective of the computational cost involved. For operator splitting, in addition to the traditional simple operator splitting (SOS) scheme, we examined the scheme developed by Douglas & Rachford (DROS). We solve three test problems (the thermal relaxation mode, the relaxation and the propagation of linear waves, and radiating shock) using these schemes and then compare their dependence on the time step size. As a result, we find the conditions of the time step size necessary for adopting each scheme. The LIN scheme is superior to other schemes if the ratio of radiation pressure to gas pressure is sufficiently low. On the other hand, DROS can be the most efficient scheme if the ratio is high. Although the NR scheme can be adopted independently of the regime, especially in a problem that involves optically thin regions, the convergence tends to be worse. In all cases, SOS is not practical.

  16. Aftermath: the implicit processes of integrating traumatic experience in the poetry of Siegfried Sassoon.

    PubMed

    Saks, Paul S

    2007-01-01

    The creation of narratives often allows individuals to bear witness to traumatic events. This study looked at connections between the processing of traumatic, affect laden experience and levels of symbolization and symmetry within the context of poetic expression. The sample for this pilot study is composed of selected works by Siegfried Sassoon (1886-1967), a British soldier-poet of the Great War. The language of the poems reflected the deepening trauma of the war experience by showing a progression toward paranoid (concrete)/symmetrical experiences. As the years passed and the poet was able to process the memory of the events, the poetry reflected a more balanced shift toward integration of depressive (symbolic)/asymmetrical experience. In terms of affect, the most significant changes were seen after Sassoon left the front and witnessed the flagrant dichotomy between civilian and military life. The results suggest a way in which traumatic events are processed. The routine horror and brutality of the Western Front initially lay outside of the realm of language and symbols and were thus highly concrete and unprocessed experiences. Time, place, and identity collapsed in on itself, leading to the increase of symmetrical experience, while the extreme "us versus them experience" of the trenches can be seen in the balance of asymmetrical experience. The study has implications for the treatment of war trauma, suggesting that writing provides a vehicle through which events can be processed and an internal sense of balance can be approached. PMID:18335612

  17. A Gauss-Kronrod-Trapezoidal integration scheme for modeling biological tissues with continuous fiber distributions.

    PubMed

    Hou, Chieh; Ateshian, Gerard A

    2016-01-01

    Fibrous biological tissues may be modeled using a continuous fiber distribution (CFD) to capture tension-compression nonlinearity, anisotropic fiber distributions, and load-induced anisotropy. The CFD framework requires spherical integration of weighted individual fiber responses, with fibers contributing to the stress response only when they are in tension. The common method for performing this integration employs the discretization of the unit sphere into a polyhedron with nearly uniform triangular faces (finite element integration or FEI scheme). Although FEI has proven to be more accurate and efficient than integration using spherical coordinates, it presents three major drawbacks: First, the number of elements on the unit sphere needed to achieve satisfactory accuracy becomes a significant computational cost in a finite element (FE) analysis. Second, fibers may not be in tension in some regions on the unit sphere, where the integration becomes a waste. Third, if tensed fiber bundles span a small region compared to the area of the elements on the sphere, a significant discretization error arises. This study presents an integration scheme specialized to the CFD framework, which significantly mitigates the first drawback of the FEI scheme, while eliminating the second and third completely. Here, integration is performed only over the regions of the unit sphere where fibers are in tension. Gauss-Kronrod quadrature is used across latitudes and the trapezoidal scheme across longitudes. Over a wide range of strain states, fiber material properties, and fiber angular distributions, results demonstrate that this new scheme always outperforms FEI, sometimes by orders of magnitude in the number of computational steps and relative accuracy of the stress calculation. PMID:26291492

  18. A comparison of two multi-variable integrator windup protection schemes

    NASA Technical Reports Server (NTRS)

    Mattern, Duane

    1993-01-01

    Two methods are examined for limit and integrator wind-up protection for multi-input, multi-output linear controllers subject to actuator constraints. The methods begin with an existing linear controller that satisfies the specifications for the nominal, small perturbation, linear model of the plant. The controllers are formulated to include an additional contribution to the state derivative calculations. The first method to be examined is the multi-variable version of the single-input, single-output, high gain, Conventional Anti-Windup (CAW) scheme. Except for the actuator limits, the CAW scheme is linear. The second scheme to be examined, denoted the Modified Anti-Windup (MAW) scheme, uses a scalar to modify the magnitude of the controller output vector while maintaining the vector direction. The calculation of the scalar modifier is a nonlinear function of the controller outputs and the actuator limits. In both cases the constrained actuator is tracked. These two integrator windup protection methods are demonstrated on a turbofan engine control system with five measurements, four control variables, and four actuators. The closed-loop responses of the two schemes are compared and contrasted during limit operation. The issue of maintaining the direction of the controller output vector using the Modified Anti-Windup scheme is discussed and the advantages and disadvantages of both of the IWP methods are presented.

  19. A unified scheme for the calculation of differentiated and undifferentiated molecular integrals over solid-harmonic Gaussians.

    PubMed

    Reine, Simen; Tellgren, Erik; Helgaker, Trygve

    2007-09-14

    Utilizing the fact that solid-harmonic combinations of Cartesian and Hermite Gaussian atomic orbitals are identical, a new scheme for the evaluation of molecular integrals over solid-harmonic atomic orbitals is presented, where the integration is carried out over Hermite rather than Cartesian atomic orbitals. Since Hermite Gaussians are defined as derivatives of spherical Gaussians, the corresponding molecular integrals become the derivatives of integrals over spherical Gaussians, whose transformation to the solid-harmonic basis is performed in the same manner as for integrals over Cartesian Gaussians, using the same expansion coefficients. The presented solid-harmonic Hermite scheme simplifies the evaluation of derivative molecular integrals, since differentiation by nuclear coordinates merely increments the Hermite quantum numbers, thereby providing a unified scheme for undifferentiated and differentiated four-center molecular integrals. For two- and three-center two-electron integrals, the solid-harmonic Hermite scheme is particularly efficient, significantly reducing the cost relative to the Cartesian scheme. PMID:17712455

  20. A secure and efficient password-based user authentication scheme using smart cards for the integrated EPR information system.

    PubMed

    Lee, Tian-Fu; Chang, I-Pin; Lin, Tsung-Hung; Wang, Ching-Cheng

    2013-06-01

    The integrated EPR information system supports convenient and rapid e-medicine services. A secure and efficient authentication scheme for the integrated EPR information system provides safeguarding patients' electronic patient records (EPRs) and helps health care workers and medical personnel to rapidly making correct clinical decisions. Recently, Wu et al. proposed an efficient password-based user authentication scheme using smart cards for the integrated EPR information system, and claimed that the proposed scheme could resist various malicious attacks. However, their scheme is still vulnerable to lost smart card and stolen verifier attacks. This investigation discusses these weaknesses and proposes a secure and efficient authentication scheme for the integrated EPR information system as alternative. Compared with related approaches, the proposed scheme not only retains a lower computational cost and does not require verifier tables for storing users' secrets, but also solves the security problems in previous schemes and withstands possible attacks. PMID:23553734

  1. Recursive multiport schemes for implementing quantum algorithms with photonic integrated circuits

    NASA Astrophysics Data System (ADS)

    Tabia, Gelo Noel M.

    2016-01-01

    We present recursive multiport schemes for implementing quantum Fourier transforms and the inversion step in Grover's algorithm on an integrated linear optics device. In particular, each scheme shows how to execute a quantum operation on 2 d modes using a pair of circuits for the same operation on d modes. The circuits operate on path-encoded qudits and realize d -dimensional unitary transformations on these states using linear optical networks with O (d2) optical elements. To evaluate the schemes against realistic errors, we ran simulations of proof-of-principle experiments using a simple fabrication model of silicon-based photonic integrated devices that employ directional couplers and thermo-optic modulators for beam splitters and phase shifters, respectively. We find that high-fidelity performance is achievable with our multiport circuits for 2-qubit and 3-qubit quantum Fourier transforms, and for quantum search on four-item and eight-item databases.

  2. Building a better leapfrog. [an algorithm for ensuring time symmetry in any integration scheme

    NASA Technical Reports Server (NTRS)

    Hut, Piet; Makino, Jun; Mcmillan, Steve

    1995-01-01

    In stellar dynamical computer simulations, as well as other types of simulations using particles, time step size is often held constant in order to guarantee a high degree of energy conservation. In many applications, allowing the time step size to change in time can offer a great saving in computational cost, but variable-size time steps usually imply a substantial degradation in energy conservation. We present a meta-algorithm' for choosing time steps in such a way as to guarantee time symmetry in any integration scheme, thus allowing vastly improved energy conservation for orbital calculations with variable time steps. We apply the algorithm to the familiar leapfrog scheme, and generalize to higher order integration schemes, showing how the stability properties of the fixed-step leapfrog scheme can be extended to higher order, variable-step integrators such as the Hermite method. We illustrate the remarkable properties of these time-symmetric integrators for the case of a highly eccentric elliptical Kepler orbit and discuss applications to more complex problems.

  3. Quantitative evaluation of numerical integration schemes for Lagrangian particle dispersion models

    NASA Astrophysics Data System (ADS)

    Ramli, Huda Mohd.; Esler, J. Gavin

    2016-07-01

    A rigorous methodology for the evaluation of integration schemes for Lagrangian particle dispersion models (LPDMs) is presented. A series of one-dimensional test problems are introduced, for which the Fokker-Planck equation is solved numerically using a finite-difference discretisation in physical space and a Hermite function expansion in velocity space. Numerical convergence errors in the Fokker-Planck equation solutions are shown to be much less than the statistical error associated with a practical-sized ensemble (N = 106) of LPDM solutions; hence, the former can be used to validate the latter. The test problems are then used to evaluate commonly used LPDM integration schemes. The results allow for optimal time-step selection for each scheme, given a required level of accuracy. The following recommendations are made for use in operational models. First, if computational constraints require the use of moderate to long time steps, it is more accurate to solve the random displacement model approximation to the LPDM rather than use existing schemes designed for long time steps. Second, useful gains in numerical accuracy can be obtained, at moderate additional computational cost, by using the relatively simple "small-noise" scheme of Honeycutt.

  4. Application of Intel Many Integrated Core (MIC) accelerators to the Pleim-Xiu land surface scheme

    NASA Astrophysics Data System (ADS)

    Huang, Melin; Huang, Bormin; Huang, Allen H.

    2015-10-01

    The land-surface model (LSM) is one physics process in the weather research and forecast (WRF) model. The LSM includes atmospheric information from the surface layer scheme, radiative forcing from the radiation scheme, and precipitation forcing from the microphysics and convective schemes, together with internal information on the land's state variables and land-surface properties. The LSM is to provide heat and moisture fluxes over land points and sea-ice points. The Pleim-Xiu (PX) scheme is one LSM. The PX LSM features three pathways for moisture fluxes: evapotranspiration, soil evaporation, and evaporation from wet canopies. To accelerate the computation process of this scheme, we employ Intel Xeon Phi Many Integrated Core (MIC) Architecture as it is a multiprocessor computer structure with merits of efficient parallelization and vectorization essentials. Our results show that the MIC-based optimization of this scheme running on Xeon Phi coprocessor 7120P improves the performance by 2.3x and 11.7x as compared to the original code respectively running on one CPU socket (eight cores) and on one CPU core with Intel Xeon E5-2670.

  5. Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme

    NASA Astrophysics Data System (ADS)

    Squire, Jonathan; Qin, Hong; Tang, William

    2012-10-01

    A fully variational, unstructured, electromagnetic particle-in-cell integrator is developed for integration of the Vlasov-Maxwell equations. Using the formalism of Discrete Exterior Calculus [1], the field solver, interpolation scheme and particle advance algorithm are derived through minimization of a single discrete field theory action. As a consequence of ensuring that the action is invariant under discrete electromagnetic gauge transformations, the integrator exactly conserves Gauss's law. This work was supported by USDOE Contract DE-AC02-09CH11466.[4pt] [1] M. Desbrun, A. N. Hirani, M. Leok, and J. E. Marsden, (2005), arXiv:math/0508341

  6. A Semi-Implicit, Three-Dimensional Model for Estuarine Circulation

    USGS Publications Warehouse

    Smith, Peter E.

    2006-01-01

    A semi-implicit, finite-difference method for the numerical solution of the three-dimensional equations for circulation in estuaries is presented and tested. The method uses a three-time-level, leapfrog-trapezoidal scheme that is essentially second-order accurate in the spatial and temporal numerical approximations. The three-time-level scheme is shown to be preferred over a two-time-level scheme, especially for problems with strong nonlinearities. The stability of the semi-implicit scheme is free from any time-step limitation related to the terms describing vertical diffusion and the propagation of the surface gravity waves. The scheme does not rely on any form of vertical/horizontal mode-splitting to treat the vertical diffusion implicitly. At each time step, the numerical method uses a double-sweep method to transform a large number of small tridiagonal equation systems and then uses the preconditioned conjugate-gradient method to solve a single, large, five-diagonal equation system for the water surface elevation. The governing equations for the multi-level scheme are prepared in a conservative form by integrating them over the height of each horizontal layer. The layer-integrated volumetric transports replace velocities as the dependent variables so that the depth-integrated continuity equation that is used in the solution for the water surface elevation is linear. Volumetric transports are computed explicitly from the momentum equations. The resulting method is mass conservative, efficient, and numerically accurate.

  7. A new computational scheme for the Dirac-Hartree-Fock method employing an efficient integral algorithm

    NASA Astrophysics Data System (ADS)

    Yanai, Takeshi; Nakajima, Takahito; Ishikawa, Yasuyuki; Hirao, Kimihiko

    2001-04-01

    A highly efficient computational scheme for four-component relativistic ab initio molecular orbital (MO) calculations over generally contracted spherical harmonic Gaussian-type spinors (GTSs) is presented. Benchmark calculations for the ground states of the group IB hydrides, MH, and dimers, M2 (M=Cu, Ag, and Au), by the Dirac-Hartree-Fock (DHF) method were performed with a new four-component relativistic ab initio MO program package oriented toward contracted GTSs. The relativistic electron repulsion integrals (ERIs), the major bottleneck in routine DHF calculations, are calculated efficiently employing the fast ERI routine SPHERICA, exploiting the general contraction scheme, and the accompanying coordinate expansion method developed by Ishida. Illustrative calculations clearly show the efficiency of our computational scheme.

  8. A secure and robust password-based remote user authentication scheme using smart cards for the integrated EPR information system.

    PubMed

    Das, Ashok Kumar

    2015-03-01

    An integrated EPR (Electronic Patient Record) information system of all the patients provides the medical institutions and the academia with most of the patients' information in details for them to make corrective decisions and clinical decisions in order to maintain and analyze patients' health. In such system, the illegal access must be restricted and the information from theft during transmission over the insecure Internet must be prevented. Lee et al. proposed an efficient password-based remote user authentication scheme using smart card for the integrated EPR information system. Their scheme is very efficient due to usage of one-way hash function and bitwise exclusive-or (XOR) operations. However, in this paper, we show that though their scheme is very efficient, their scheme has three security weaknesses such as (1) it has design flaws in password change phase, (2) it fails to protect privileged insider attack and (3) it lacks the formal security verification. We also find that another recently proposed Wen's scheme has the same security drawbacks as in Lee at al.'s scheme. In order to remedy these security weaknesses found in Lee et al.'s scheme and Wen's scheme, we propose a secure and efficient password-based remote user authentication scheme using smart cards for the integrated EPR information system. We show that our scheme is also efficient as compared to Lee et al.'s scheme and Wen's scheme as our scheme only uses one-way hash function and bitwise exclusive-or (XOR) operations. Through the security analysis, we show that our scheme is secure against possible known attacks. Furthermore, we simulate our scheme for the formal security verification using the widely-accepted AVISPA (Automated Validation of Internet Security Protocols and Applications) tool and show that our scheme is secure against passive and active attacks. PMID:25666924

  9. Computational Analysis of Muscular Dystrophy Sub-types Using A Novel Integrative Scheme

    PubMed Central

    Wang, Chen; Ha, Sook; Wang, Yue; Hoffman, Eric

    2012-01-01

    To construct biologically interpretable gene sets for muscular dystrophy (MD) sub-type classification, we propose a novel computational scheme to integrate protein-protein interaction (PPI) network, functional gene set information, and mRNA profiling data. The workflow of the proposed scheme includes the following three major steps: firstly, we apply an affinity propagation clustering (APC) approach to identify gene sub-networks associated with each MD sub-type, in which a new distance metric is proposed for APC to combine PPI network information and gene-gene co-expression relationship; secondly, we further incorporate functional gene set knowledge, which complements the physical PPI information, into our scheme for biomarker identification; finally, based on the constructed sub-networks and gene set features, we apply multi-class support vector machines (MSVMs) for MD sub-type classification, with which to highlight the biomarkers contributing to sub-type prediction. The experimental results show that our scheme can help identify sub-networks and gene sets that are more relevant to MD than those constructed by other conventional approaches. Moreover, our integrative strategy improves the prediction accuracy substantially, especially for those ’hard-to-classify’ sub-types. PMID:22773895

  10. High-order sampling schemes for path integrals and Gaussian chain simulations of polymers

    SciTech Connect

    Müser, Martin H.; Müller, Marcus

    2015-05-07

    In this work, we demonstrate that path-integral schemes, derived in the context of many-body quantum systems, benefit the simulation of Gaussian chains representing polymers. Specifically, we show how to decrease discretization corrections with little extra computation from the usual O(1/P{sup 2}) to O(1/P{sup 4}), where P is the number of beads representing the chains. As a consequence, high-order integrators necessitate much smaller P than those commonly used. Particular emphasis is placed on the questions of how to maintain this rate of convergence for open polymers and for polymers confined by a hard wall as well as how to ensure efficient sampling. The advantages of the high-order sampling schemes are illustrated by studying the surface tension of a polymer melt and the interface tension in a binary homopolymers blend.

  11. A Hash Based Remote User Authentication and Authenticated Key Agreement Scheme for the Integrated EPR Information System.

    PubMed

    Li, Chun-Ta; Weng, Chi-Yao; Lee, Cheng-Chi; Wang, Chun-Cheng

    2015-11-01

    To protect patient privacy and ensure authorized access to remote medical services, many remote user authentication schemes for the integrated electronic patient record (EPR) information system have been proposed in the literature. In a recent paper, Das proposed a hash based remote user authentication scheme using passwords and smart cards for the integrated EPR information system, and claimed that the proposed scheme could resist various passive and active attacks. However, in this paper, we found that Das's authentication scheme is still vulnerable to modification and user duplication attacks. Thereafter we propose a secure and efficient authentication scheme for the integrated EPR information system based on lightweight hash function and bitwise exclusive-or (XOR) operations. The security proof and performance analysis show our new scheme is well-suited to adoption in remote medical healthcare services. PMID:26349803

  12. Increasing the multiscale/multiphysics capability of CAM-SE using implicit time integration and GPU accelerators

    NASA Astrophysics Data System (ADS)

    Archibald, R.; Evans, K. J.; Worley, P.; Norman, M. R.; Lott, A.; Salinger, A.; Woodward, C. S.

    2014-12-01

    The recent focus on regional refinement in the Community Atmosphere Model (CAM5) has created a strong need to develop time-stepping methods capable of accelerating throughput on high performance computing for climate dynamics across multiple spatial and temporal scales. This research is focused on developing implicit methods that can be executed at scale on GPU based machines. Efforts to port the scalable spectral element dynamical core to incorporate these developments is presented, including both 2D and 3D benchmark test case results. The current implicit solver and preconditioner implementations utilize a Fortran interface package within the Trilinos project, third party software that allows fully tested, optimized, and robust code with a suite of parameter options to be included a priori. Merging this coding strategy with GPU libraries will be discussed along with beneficial optimization gains and data structure requirements to evaluate Trilinos binded residual calculations on GPU processors.

  13. Development of highly accurate approximate scheme for computing the charge transfer integral

    SciTech Connect

    Pershin, Anton; Szalay, Péter G.

    2015-08-21

    The charge transfer integral is a key parameter required by various theoretical models to describe charge transport properties, e.g., in organic semiconductors. The accuracy of this important property depends on several factors, which include the level of electronic structure theory and internal simplifications of the applied formalism. The goal of this paper is to identify the performance of various approximate approaches of the latter category, while using the high level equation-of-motion coupled cluster theory for the electronic structure. The calculations have been performed on the ethylene dimer as one of the simplest model systems. By studying different spatial perturbations, it was shown that while both energy split in dimer and fragment charge difference methods are equivalent with the exact formulation for symmetrical displacements, they are less efficient when describing transfer integral along the asymmetric alteration coordinate. Since the “exact” scheme was found computationally expensive, we examine the possibility to obtain the asymmetric fluctuation of the transfer integral by a Taylor expansion along the coordinate space. By exploring the efficiency of this novel approach, we show that the Taylor expansion scheme represents an attractive alternative to the “exact” calculations due to a substantial reduction of computational costs, when a considerably large region of the potential energy surface is of interest. Moreover, we show that the Taylor expansion scheme, irrespective of the dimer symmetry, is very accurate for the entire range of geometry fluctuations that cover the space the molecule accesses at room temperature.

  14. Development of highly accurate approximate scheme for computing the charge transfer integral.

    PubMed

    Pershin, Anton; Szalay, Péter G

    2015-08-21

    The charge transfer integral is a key parameter required by various theoretical models to describe charge transport properties, e.g., in organic semiconductors. The accuracy of this important property depends on several factors, which include the level of electronic structure theory and internal simplifications of the applied formalism. The goal of this paper is to identify the performance of various approximate approaches of the latter category, while using the high level equation-of-motion coupled cluster theory for the electronic structure. The calculations have been performed on the ethylene dimer as one of the simplest model systems. By studying different spatial perturbations, it was shown that while both energy split in dimer and fragment charge difference methods are equivalent with the exact formulation for symmetrical displacements, they are less efficient when describing transfer integral along the asymmetric alteration coordinate. Since the "exact" scheme was found computationally expensive, we examine the possibility to obtain the asymmetric fluctuation of the transfer integral by a Taylor expansion along the coordinate space. By exploring the efficiency of this novel approach, we show that the Taylor expansion scheme represents an attractive alternative to the "exact" calculations due to a substantial reduction of computational costs, when a considerably large region of the potential energy surface is of interest. Moreover, we show that the Taylor expansion scheme, irrespective of the dimer symmetry, is very accurate for the entire range of geometry fluctuations that cover the space the molecule accesses at room temperature. PMID:26298117

  15. Numerical integration of the stochastic Landau-Lifshitz-Gilbert equation in generic time-discretization schemes.

    PubMed

    Romá, Federico; Cugliandolo, Leticia F; Lozano, Gustavo S

    2014-08-01

    We introduce a numerical method to integrate the stochastic Landau-Lifshitz-Gilbert equation in spherical coordinates for generic discretization schemes. This method conserves the magnetization modulus and ensures the approach to equilibrium under the expected conditions. We test the algorithm on a benchmark problem: the dynamics of a uniformly magnetized ellipsoid. We investigate the influence of various parameters, and in particular, we analyze the efficiency of the numerical integration, in terms of the number of steps needed to reach a chosen long time with a given accuracy. PMID:25215839

  16. Validation Study of the Integral-Differential Scheme for Multi-Block Grids

    NASA Astrophysics Data System (ADS)

    Mrema, Honest Frank

    This MS Thesis seeks to validate the accuracy of the Integral-Differential Scheme (IDS). In the attempts to accomplish this task, research efforts were focused on the scheme's ability to capture the physics of known flow fields, as well as the scheme's ability to predict the features of flow field quantities that may be derived from experimental measurements. The IDS was developed with the goal of being computationally efficient, from a programming perspective, as well as being numerically accurate, stable, and robust, from a mathematical perspective. The IDS is designed to solve the full Navier-Stokes equations in their integral forms. Unlike traditional control volume schemes, the IDS is built upon two sets of cells: spatial and temporal cells. For 2-D flows, the IDS considers an elementary control volume as a collection of four spatial cells and a single temporal cell. Similar to other explicit CFD schemes, the IDS relies on the use of the Taylor series expansion and other traditional CFD criteria. It is of interest to note that there are previous IDS validation studies which were conducted at North Carolina A&T State University. These past studies mainly focused on the qualitative aspects of the flow field physics. Furthermore, in all cases, they focused on flow field problems that can be represented by single-block grids. In this analysis, the validation studies are focused on multi-block grids in which the physics of the flow field is made complicated due to the presence of shock waves and flow separation zones. Of interest to this MS Thesis are two supersonic flow field problems that are supported by experimental data; namely, the supersonic flow over a rearward-facing step problem and the supersonic flow over a cavity problem. The validation studies conducted herein demonstrated that the IDS was able to predict the experimental data in both cases.

  17. An Implicit Characteristic Based Method for Electromagnetics

    NASA Technical Reports Server (NTRS)

    Beggs, John H.; Briley, W. Roger

    2001-01-01

    An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

  18. A multi-dimensional nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell (PIC) algorithm

    NASA Astrophysics Data System (ADS)

    Chen, Guangye; Chacón, Luis; CoCoMans Team

    2014-10-01

    For decades, the Vlasov-Darwin model has been recognized to be attractive for PIC simulations (to avoid radiative noise issues) in non-radiative electromagnetic regimes. However, the Darwin model results in elliptic field equations that renders explicit time integration unconditionally unstable. Improving on linearly implicit schemes, fully implicit PIC algorithms for both electrostatic and electromagnetic regimes, with exact discrete energy and charge conservation properties, have been recently developed in 1D. This study builds on these recent algorithms to develop an implicit, orbit-averaged, time-space-centered finite difference scheme for the particle-field equations in multiple dimensions. The algorithm conserves energy, charge, and canonical-momentum exactly, even with grid packing. A simple fluid preconditioner allows efficient use of large timesteps, O (√{mi/me}c/veT) larger than the explicit CFL. We demonstrate the accuracy and efficiency properties of the of the algorithm with various numerical experiments in 2D3V.

  19. Advanced integration schemes for high-functionality/high-performance photonic integrated circuits

    NASA Astrophysics Data System (ADS)

    Raring, James W.; Sysak, Matthew N.; Tauke-Pedretti, Anna; Dummer, Matthew; Skogen, Erik J.; Barton, Jonathon S.; DenBaars, S. P.; Coldren, Larry A.

    2006-02-01

    The evolution of optical communication systems has facilitated the required bandwidth to meet the increasing data rate demands. However, as the peripheral technologies have progressed to meet the requirements of advanced systems, an abundance of viable solutions and products have emerged. The finite market for these products will inevitably force a paradigm shift upon the communications industry. Monolithic integration is a key technology that will facilitate this shift as it will provide solutions at low cost with reduced power dissipation and foot-print in the form of highly functional optical components based on photonic integrated circuits (PICs). In this manuscript, we discuss the advantages, potential applications, and challenges of photonic integration. After a brief overview of various integration techniques, we present our novel approaches to increase the performance of the individual components comprising highly functional PICs.

  20. Comparison of Aircraft Models and Integration Schemes for Interval Management in the TRACON

    NASA Technical Reports Server (NTRS)

    Neogi, Natasha; Hagen, George E.; Herencia-Zapana, Heber

    2012-01-01

    Reusable models of common elements for communication, computation, decision and control in air traffic management are necessary in order to enable simulation, analysis and assurance of emergent properties, such as safety and stability, for a given operational concept. Uncertainties due to faults, such as dropped messages, along with non-linearities and sensor noise are an integral part of these models, and impact emergent system behavior. Flight control algorithms designed using a linearized version of the flight mechanics will exhibit error due to model uncertainty, and may not be stable outside a neighborhood of the given point of linearization. Moreover, the communication mechanism by which the sensed state of an aircraft is fed back to a flight control system (such as an ADS-B message) impacts the overall system behavior; both due to sensor noise as well as dropped messages (vacant samples). Additionally simulation of the flight controller system can exhibit further numerical instability, due to selection of the integration scheme and approximations made in the flight dynamics. We examine the theoretical and numerical stability of a speed controller under the Euler and Runge-Kutta schemes of integration, for the Maintain phase for a Mid-Term (2035-2045) Interval Management (IM) Operational Concept for descent and landing operations. We model uncertainties in communication due to missed ADS-B messages by vacant samples in the integration schemes, and compare the emergent behavior of the system, in terms of stability, via the boundedness of the final system state. Any bound on the errors incurred by these uncertainties will play an essential part in a composable assurance argument required for real-time, flight-deck guidance and control systems,. Thus, we believe that the creation of reusable models, which possess property guarantees, such as safety and stability, is an innovative and essential requirement to assessing the emergent properties of novel airspace

  1. Impacts of biasing schemes in the one-loop integrated perturbation theory

    NASA Astrophysics Data System (ADS)

    Matsubara, Takahiko; Desjacques, Vincent

    2016-06-01

    The impact of biasing schemes on the clustering of tracers of the large-scale structure is analytically studied in the weakly nonlinear regime. For this purpose, we use the one-loop approximation of the integrated perturbation theory together with the renormalized bias functions of various, physically motivated Lagrangian bias schemes. These include the halo, peaks, and excursion set peaks model, for which we derive useful formulas for the evaluation of their renormalized bias functions. The shapes of the power spectra and correlation functions are affected by the different bias models at the level of a few percent on weakly nonlinear scales. These effects are studied quantitatively both in real and redshift space. The amplitude of the scale-dependent bias in the presence of primordial non-Gaussianity also depends on the details of the bias models. If left unaccounted for, these theoretical uncertainties could affect the robustness of the cosmological constraints extracted from galaxy clustering data.

  2. Evaluation of exercise-respiratory system modifications and integration schemes for physiological systems

    NASA Technical Reports Server (NTRS)

    Gallagher, R. R.

    1974-01-01

    Exercise subroutine modifications are implemented in an exercise-respiratory system model yielding improvement of system response to exercise forcings. A more physiologically desirable respiratory ventilation rate in addition to an improved regulation of arterial gas tensions and cerebral blood flow is observed. A respiratory frequency expression is proposed which would be appropriate as an interfacing element of the respiratory-pulsatile cardiovascular system. Presentation of a circulatory-respiratory system integration scheme along with its computer program listing is given. The integrated system responds to exercise stimulation for both nonstressed and stressed physiological states. Other integration possibilities are discussed with respect to the respiratory, pulsatile cardiovascular, thermoregulatory, and the long-term circulatory systems.

  3. Development of an integrated computerized scheme for metaphase chromosome image analysis: a robustness experiment

    NASA Astrophysics Data System (ADS)

    Wang, Xingwei; Zheng, Bin; Li, Shibo; Mulvihill, John J.; Wood, Marc C.; Yuan, Chaowei; Chen, Wei; Liu, Hong

    2008-02-01

    Our integrated computer-aided detection (CAD) scheme includes three basic modules. The first module detects whether a microscopic digital image depicts a metaphase chromosome cell. If a cell is detected, the scheme will justify whether it is analyzable with a decision tree. Once an analyzable cell is detected, the second module is applied to segment individual chromosomes and to compute two important features. Specifically, the scheme utilizes a modified thinning algorithm to identify the medial axis of a chromosome. By tracking perpendicular lines along the medial axis, the scheme computes four feature profiles, identifies centromeres, and assigns polarities of chromosomes based on a set of pre-optimized rules. The third module is followed to classify chromosomes into 24 types. In this module, each chromosome is initially represented by a vector of 31 features. A two-layer classifier with 8 artificial neural networks (ANN) is optimized by a genetic algorithm. A testing chromosome is first classified into one of the seven groups by the ANN in the first layer. Another ANN is then automatically selected from the seven ANNs in the second layer (one for each group) to further classify this chromosome into one of 24 types. To test the performance and robustness of this CAD scheme, we randomly selected and assembled an independent testing dataset. The dataset contains 100 microscopic digital images including 50 analyzable and 50 un-analyzable metphase cells identified by the experts. The centromere location, the corresponding polarity, and karyotype for each individual chromosome were recorded in the "truth" file. The performance of the CAD scheme applied to this image dataset is analyzed and compared with the results in the true file. The assessment accuracies are 93% for the first module, 90.8% for centromere identification and 93.2% for polarity assignment in the second module, over 96% for six chromosome groups and 81.8% for one group in the third module

  4. Integration of Full Particle Orbit in Toroidal Plasmas Using Boris Scheme

    NASA Astrophysics Data System (ADS)

    Wei, Xishuo; Xiao, Yong

    2014-10-01

    When studying particle dynamics in high frequency electromagnetic waves, such as low hybrid wave heating, it is important to integrate full particle orbit accurately to very long time in tokamaks. Here we derived a formulation under magnetic coordinate based on the Boris Scheme, which can be used effectively to push particles in long time scale. After several hundred gyro-periods, the banana orbit can be observed and the toroidal precession frequency can be measured. The toroidal precession frequency is found to match that from the guiding center simulation. This new method shows superior numeric properties than the traditional Runge-Kutta method in terms of conserving particle energy and magnetic moment.

  5. An evaluation of explicit time integration schemes for use with the generalized interpolation material point method

    NASA Astrophysics Data System (ADS)

    Wallstedt, P. C.; Guilkey, J. E.

    2008-11-01

    The stability and accuracy of the generalized interpolation material point (GIMP) Method is measured directly through carefully-formulated manufactured solutions over wide ranges of CFL numbers and mesh sizes. The manufactured solutions are described in detail. The accuracy and stability of several time integration schemes are compared via numerical experiments. The effect of various treatments of particle "size" are also considered. The hypothesis that GIMP is most accurate when particles remain contiguous and non-overlapping is confirmed by comparing manufactured solutions with and without this property.

  6. Integrating AlGaN/GaN high electron mobility transistor with Si: A comparative study of integration schemes

    SciTech Connect

    Mohan, Nagaboopathy; Raghavan, Srinivasan; Manikant,; Soman, Rohith

    2015-10-07

    AlGaN/GaN high electron mobility transistor stacks deposited on a single growth platform are used to compare the most common transition, AlN to GaN, schemes used for integrating GaN with Si. The efficiency of these transitions based on linearly graded, step graded, interlayer, and superlattice schemes on dislocation density reduction, stress management, surface roughness, and eventually mobility of the 2D-gas are evaluated. In a 500 nm GaN probe layer deposited, all of these transitions result in total transmission electron microscopy measured dislocations densities of 1 to 3 × 10{sup 9}/cm{sup 2} and <1 nm surface roughness. The 2-D electron gas channels formed at an AlGaN-1 nm AlN/GaN interface deposited on this GaN probe layer all have mobilities of 1600–1900 cm{sup 2}/V s at a carrier concentration of 0.7–0.9 × 10{sup 13}/cm{sup 2}. Compressive stress and changes in composition in GaN rich regions of the AlN-GaN transition are the most effective at reducing dislocation density. Amongst all the transitions studied the step graded transition is the one that helps to implement this feature of GaN integration in the simplest and most consistent manner.

  7. Integrating AlGaN/GaN high electron mobility transistor with Si: A comparative study of integration schemes

    NASA Astrophysics Data System (ADS)

    Mohan, Nagaboopathy; Singh, Manikant; Soman, Rohith; Raghavan, Srinivasan

    2015-10-01

    AlGaN/GaN high electron mobility transistor stacks deposited on a single growth platform are used to compare the most common transition, AlN to GaN, schemes used for integrating GaN with Si. The efficiency of these transitions based on linearly graded, step graded, interlayer, and superlattice schemes on dislocation density reduction, stress management, surface roughness, and eventually mobility of the 2D-gas are evaluated. In a 500 nm GaN probe layer deposited, all of these transitions result in total transmission electron microscopy measured dislocations densities of 1 to 3 × 109/cm2 and <1 nm surface roughness. The 2-D electron gas channels formed at an AlGaN-1 nm AlN/GaN interface deposited on this GaN probe layer all have mobilities of 1600-1900 cm2/V s at a carrier concentration of 0.7-0.9 × 1013/cm2. Compressive stress and changes in composition in GaN rich regions of the AlN-GaN transition are the most effective at reducing dislocation density. Amongst all the transitions studied the step graded transition is the one that helps to implement this feature of GaN integration in the simplest and most consistent manner.

  8. Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations. Chapter 1

    NASA Technical Reports Server (NTRS)

    Yee, Helen C.; Sjoegreen, B.; Mansour, Nagi N. (Technical Monitor)

    2001-01-01

    A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible; (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred; (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used; and (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for nonlinear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows. The minimization of employing very fine grids to overcome the production of spurious numerical solutions and/or instability due to under-resolved grids is also sought. The incremental studies to illustrate the performance of the approach are summarized. Extensive testing and full implementation of the approach is forthcoming. The results shown so far are very encouraging.

  9. Challenges and mitigation strategies for resist trim etch in resist-mandrel based SAQP integration scheme

    NASA Astrophysics Data System (ADS)

    Mohanty, Nihar; Franke, Elliott; Liu, Eric; Raley, Angelique; Smith, Jeffrey; Farrell, Richard; Wang, Mingmei; Ito, Kiyohito; Das, Sanjana; Ko, Akiteru; Kumar, Kaushik; Ranjan, Alok; O'Meara, David; Nawa, Kenjiro; Scheer, Steven; DeVillers, Anton; Biolsi, Peter

    2015-03-01

    Patterning the desired narrow pitch at 10nm technology node and beyond, necessitates employment of either extreme ultra violet (EUV) lithography or multi-patterning solutions based on 193nm-immersion lithography. With enormous challenges being faced in getting EUV lithography ready for production, multi-patterning solutions that leverage the already installed base of 193nm-immersion-lithography are poised to become the industry norm for 10 and 7nm technology nodes. For patterning sub-40nm pitch line/space features, self-aligned quadruple patterning (SAQP) with resist pattern as the first mandrel shows significant cost as well as design benefit, as compared to EUV lithography or other multi-patterning techniques. One of the most critical steps in this patterning scheme is the resist mandrel definition step which involves trimming / reformation of resist profile via plasma etch for achieving appropriate pitch after the final pattern. Being the first mandrel, the requirements for the Line Edge Roughness (LER) / Line Width Roughness (LWR); critical dimension uniformity (CDU); and profile in 3-dimensions for the resist trim / reformation etch is extremely aggressive. In this paper we highlight the unique challenges associated in developing resist trim / reformation plasma etch process for SAQP integration scheme and summarize our efforts in optimizing the trim etch chemistries, process steps and plasma etch parameters for meeting the mandrel definition targets. Finally, we have shown successful patterning of 30nm pitch patterns via the resist-mandrel SAQP scheme and its implementation for Si-fin formation at 7nm node.

  10. Numerical simulation of Stokes flow around particles via a hybrid Finite Difference-Boundary Integral scheme

    NASA Astrophysics Data System (ADS)

    Bhattacharya, Amitabh

    2013-11-01

    An efficient algorithm for simulating Stokes flow around particles is presented here, in which a second order Finite Difference method (FDM) is coupled to a Boundary Integral method (BIM). This method utilizes the strong points of FDM (i.e. localized stencil) and BIM (i.e. accurate representation of particle surface). Specifically, in each iteration, the flow field away from the particles is solved on a Cartesian FDM grid, while the traction on the particle surface (given the the velocity of the particle) is solved using BIM. The two schemes are coupled by matching the solution in an intermediate region between the particle and surrounding fluid. We validate this method by solving for flow around an array of cylinders, and find good agreement with Hasimoto's (J. Fluid Mech. 1959) analytical results.

  11. An Inverter Packaging Scheme for an Integrated Segmented Traction Drive System

    SciTech Connect

    Su, Gui-Jia; Tang, Lixin; Ayers, Curtis William; Wiles, Randy H

    2013-01-01

    The standard voltage source inverter (VSI), widely used in electric vehicle/hybrid electric vehicle (EV/HEV) traction drives, requires a bulky dc bus capacitor to absorb the large switching ripple currents and prevent them from shortening the battery s life. The dc bus capacitor presents a significant barrier to meeting inverter cost, volume, and weight requirements for mass production of affordable EVs/HEVs. The large ripple currents become even more problematic for the film capacitors (the capacitor technology of choice for EVs/HEVs) in high temperature environments as their ripple current handling capability decreases rapidly with rising temperatures. It is shown in previous work that segmenting the VSI based traction drive system can significantly decrease the ripple currents and thus the size of the dc bus capacitor. This paper presents an integrated packaging scheme to reduce the system cost of a segmented traction drive.

  12. A variational implementation of the implicit particle filter for the shallow water equations

    NASA Astrophysics Data System (ADS)

    Souopgui, I.; Morzfeld, M.; Hussaini, M.; Chorin, A. J.

    2013-12-01

    The estimation of initial conditions for shallow water equations is a well known test problem for operational data assimilation techniques. The state-of-the-art approach to this problem is the variational method (4D-Var), i.e. the computation of the mode of the posterior probability density function (pdf) via the adjoint technique. We add a sampling step to the variational method, thus turning a computation of the conditional mode (a biased estimator) into a computation of the conditional mean (the minimum least square error estimator). Our implementation relies on implicit sampling, which is a Monte Carlo (MC) sampling scheme. The idea in implicit sampling is to first search for the high-probability region of the posterior pdf and then to find samples in this region. Because the samples are concentrated in the high-probability region, fewer samples are required than with competing MC schemes and, thus, implicit sampling can be more efficient than other MC schemes. The search for the high-probability region can be done via a numerical minimization that is very similar to the minimization in 4D-Var. Here, we use existing 4D-Var code to implement the implicit sampling scheme. Once the minimization problem is solved, we obtain samples by solving algebraic equations with a random right-hand-side. These equations can be solved efficiently, so that the additional cost of our approach, compared to 4D-Var, is small. We present numerical experiments to demonstrate the applicability and efficiency of our approach. These numerical experiments mimic physical experiments done with the CORIOLIS turntable in Grenoble (France), which are used to study the drift of a vortex. In particular we consider shallow water equations on a square domain (2.5m x 2.5m) with open boundary conditions and discretize the equations with finite differences on a staggered grid of size 256 x 256 and a fourth order Runge-Kutta time integrator. Our goal is to estimate the initial state (velocities and

  13. Low-k/copper integration scheme suitable for ULSI manufacturing from 90nm to 45nm nodes

    NASA Astrophysics Data System (ADS)

    Nogami, T.; Lane, S.; Fukasawa, M.; Ida, K.; Angyal, M.; Chanda, K.; Chen, F.; Christiansen, C.; Cohen, S.; Cullinan, M.; Dziobkowski, C.; Fitzsimmons, J.; Flaitz, P.; Grill, A.; Gill, J.; Inoue, K.; Klymko, N.; Kumar, K.; Labelle, C.; Lane, M.; Li, B.; Liniger, E.; Madon, A.; Malone, K.; Martin, J.; McGahay, V.; McLaughlin, P.; Melville, I.; Minami, M.; Molis, S.; Nguyen, S.; Penny, C.; Restaino, D.; Sakamoto, A.; Sankar, M.; Sherwood, M.; Simonyi, E.; Shimooka, Y.; Tai, L.; Widodo, J.; Wildman, H.; Ono, M.; McHerron, D.; Nye, H.; Davis, C.; Sankaran, S.; Edelstein, D.; Ivers, T.

    2005-11-01

    This paper discusses low-k/copper integration schemes which has been in production in the 90 nm node, have been developed in the 65 nm node, and should be taken in the 45 nm node. While our baseline 65 nm BEOL process has been developed by extension and simple shrinkage of our PECVD SiCOH integration which has been in production in the 90 nm node with our SiCOH film having k=3.0, the 65 nm SiCOH integration has two other options to go to extend to lower capacitance. One is to add porosity to become ultra low-k (ULK). The other is to stay with low-k SiCOH, which is modified to have a "lower-k". The effective k- value attained with the lower-k (k=2.8) SiCOH processed in the "Direct CMP" scheme is very close to that with an ULK (k=2.5) SiCOH film built with the "Hard Mask Retention" scheme. This paper first describes consideration of these two damascene schemes, whose comparison leads to the conclusion that the lower-k SiCOH integration can have more advantages in terms of process simplicity and extendibility of our 90 nm scheme under certain assumptions. Then describing the k=2.8 SiCOH film development and its successful integration, damascene schemes for 45nm nodes are discussed based on our learning from development of the lower-k 65nm scheme. Capability of modern dry etchers to define the finer patterns, non-uniformity of CMP, and susceptibility to plasma and mechanical strength and adhesion of ULK are discussed as factors to hamper the applicability of ULK.

  14. High-order implicit time-marching methods for unsteady fluid flow simulation

    NASA Astrophysics Data System (ADS)

    Boom, Pieter David

    Unsteady computational fluid dynamics (CFD) is increasingly becoming a critical tool in the development of emerging technologies and modern aircraft. In spite of rapid mathematical and technological advancement, these simulations remain computationally intensive and time consuming. More efficient temporal integration will promote a wider use of unsteady analysis and extend its range of applicability. This thesis presents an investigation of efficient high-order implicit time-marching methods for application in unsteady compressible CFD. A generalisation of time-marching methods based on summation-by-parts (SBP) operators is described which reduces the number of stages required to obtain a prescribed order of accuracy, thus improving their efficiency. The classical accuracy and stability theory is formally extended for these generalised SBP (GSBP) methods, including superconvergence and nonlinear stability. Dual-consistent SBP and GSBP time-marching methods are shown to form a subclass of implicit Runge-Kutta methods, which enables extensions of nonlinear accuracy and stability results. A novel family of fully-implicit GSBP Runge-Kutta schemes based on Gauss quadrature are derived which are both algebraically stable and L-stable with order 2s - 1, where s is the number of stages. In addition, a numerical tool is developed for the construction and optimisation of general linear time-marching methods. The tool is applied to the development of several low-stage-order L-stable diagonally-implicit methods, including a diagonally-implicit GSBP Runge-Kutta scheme. The most notable and efficient method developed is a six-stage fifth-order L-stable stiffly-accurate explicit-first-stage singly-diagonally-implicit Runge-Kutta (ESDIRK5) method with stage order two. The theoretical results developed in this thesis are supported by numerical simulations, and the predicted relative efficiency of the schemes is realised.

  15. Development and evaluation of a building energy model integrated in the TEB scheme

    NASA Astrophysics Data System (ADS)

    Bueno, B.; Pigeon, G.; Norford, L. K.; Zibouche, K.; Marchadier, C.

    2012-03-01

    The use of air-conditioning systems is expected to increase as a consequence of global-scale and urban-scale climate warming. In order to represent future scenarios of urban climate and building energy consumption, the Town Energy Balance (TEB) scheme must be improved. This paper presents a new building energy model (BEM) that has been integrated in the TEB scheme. BEM-TEB makes it possible to represent the energy effects of buildings and building systems on the urban climate and to estimate the building energy consumption at city scale (~10 km) with a resolution of a neighbourhood (~100 m). The physical and geometric definition of buildings in BEM has been intentionally kept as simple as possible, while maintaining the required features of a comprehensive building energy model. The model considers a single thermal zone, where the thermal inertia of building materials associated with multiple levels is represented by a generic thermal mass. The model accounts for heat gains due to transmitted solar radiation, heat conduction through the enclosure, infiltration, ventilation, and internal heat gains. BEM allows for previously unavailable sophistication in the modelling of air-conditioning systems. It accounts for the dependence of the system capacity and efficiency on indoor and outdoor air temperatures and solves the dehumidification of the air passing through the system. Furthermore, BEM includes specific models for passive systems, such as window shadowing devices and natural ventilation. BEM has satisfactorily passed different evaluation processes, including testing its modelling assumptions, verifying that the chosen equations are solved correctly, and validating the model with field data.

  16. Development and evaluation of a building energy model integrated in the TEB scheme

    NASA Astrophysics Data System (ADS)

    Bueno, B.; Pigeon, G.; Norford, L. K.; Zibouche, K.

    2011-11-01

    The use of air-conditioning systems is expected to increase as a consequence of global-scale and urban-scale climate warming. In order to represent future scenarios of urban climate and building energy consumption, the Town Energy Budget (TEB) scheme must be improved. This paper presents a new building energy model (BEM) that has been integrated in the TEB scheme. BEM-TEB makes it possible to represent the energy effects of buildings and building systems on the urban climate and to estimate the building energy consumption at city scale (~10 km) with a resolution of a neighbourhood (~100 m). The physical and geometric definition of buildings in BEM has been intentionally kept as simple as possible, while maintaining the required features of a comprehensive building energy model. The model considers a single thermal zone, where the thermal inertia of building materials associated with multiple levels is represented by a generic thermal mass. The model accounts for heat gains due to transmitted solar radiation, heat conduction through the enclosure, infiltration, ventilation, and internal heat gains. As a difference with respect to other building parameterizations used in urban climate, BEM includes specific models for real air-conditioning systems. It accounts for the dependence of the system capacity and efficiency on indoor and outdoor air temperatures and solves the dehumidification of the air passing through the system. Furthermore, BEM includes specific models for passive systems, such as window shadowing devices and natural ventilation. BEM has satisfactorily passed different evaluation processes, including testing its modelling assumptions, verifying that the chosen equations are solved correctly, and validating the model with field data.

  17. An extension of A-stability to alternating direction implicit methods

    NASA Technical Reports Server (NTRS)

    Warming, R. F.; Beam, R. M.

    1978-01-01

    An alternating direction implicit (ADI) scheme was constructed by the method of approximate factorization. An A-stable linear multistep method (LMM) was used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation. Sufficient conditions for the A-stability of the LMM were determined by applying the theory of positive real functions to reduce the stability analysis of the partial differential equations to a simple algebraic test. A linear test equation for partial differential equations is defined and then used to analyze the stability of approximate factorization schemes. An ADI method for the three-dimensional heat equation is also presented.

  18. A whole-path importance-sampling scheme for Feynman path integral calculations of absolute partition functions and free energies.

    PubMed

    Mielke, Steven L; Truhlar, Donald G

    2016-01-21

    Using Feynman path integrals, a molecular partition function can be written as a double integral with the inner integral involving all closed paths centered at a given molecular configuration, and the outer integral involving all possible molecular configurations. In previous work employing Monte Carlo methods to evaluate such partition functions, we presented schemes for importance sampling and stratification in the molecular configurations that constitute the path centroids, but we relied on free-particle paths for sampling the path integrals. At low temperatures, the path sampling is expensive because the paths can travel far from the centroid configuration. We now present a scheme for importance sampling of whole Feynman paths based on harmonic information from an instantaneous normal mode calculation at the centroid configuration, which we refer to as harmonically guided whole-path importance sampling (WPIS). We obtain paths conforming to our chosen importance function by rejection sampling from a distribution of free-particle paths. Sample calculations on CH4 demonstrate that at a temperature of 200 K, about 99.9% of the free-particle paths can be rejected without integration, and at 300 K, about 98% can be rejected. We also show that it is typically possible to reduce the overhead associated with the WPIS scheme by sampling the paths using a significantly lower-order path discretization than that which is needed to converge the partition function. PMID:26801023

  19. A whole-path importance-sampling scheme for Feynman path integral calculations of absolute partition functions and free energies

    NASA Astrophysics Data System (ADS)

    Mielke, Steven L.; Truhlar, Donald G.

    2016-01-01

    Using Feynman path integrals, a molecular partition function can be written as a double integral with the inner integral involving all closed paths centered at a given molecular configuration, and the outer integral involving all possible molecular configurations. In previous work employing Monte Carlo methods to evaluate such partition functions, we presented schemes for importance sampling and stratification in the molecular configurations that constitute the path centroids, but we relied on free-particle paths for sampling the path integrals. At low temperatures, the path sampling is expensive because the paths can travel far from the centroid configuration. We now present a scheme for importance sampling of whole Feynman paths based on harmonic information from an instantaneous normal mode calculation at the centroid configuration, which we refer to as harmonically guided whole-path importance sampling (WPIS). We obtain paths conforming to our chosen importance function by rejection sampling from a distribution of free-particle paths. Sample calculations on CH4 demonstrate that at a temperature of 200 K, about 99.9% of the free-particle paths can be rejected without integration, and at 300 K, about 98% can be rejected. We also show that it is typically possible to reduce the overhead associated with the WPIS scheme by sampling the paths using a significantly lower-order path discretization than that which is needed to converge the partition function.

  20. Integrated subsurface water solutions for coastal environments through integrated pump&treat and aquifer storage and recovery (ASR) schemes

    NASA Astrophysics Data System (ADS)

    Perdikaki, Martha; Kallioras, Andreas; Christoforidis, Christophoros; Iossifidis, Dimitris; Zafeiropoulos, Anastasios; Dimitriadis, Klisthenis; Makropoulos, Christos; Raat, Klaasjan; van den Berg, Gerard

    2016-04-01

    Coastal wetlands in semi-arid regions, as in Circum-Mediterranean, are considered important ecosystems that provide valuable services to human population and the environment, such as: flood protection, erosion control, wildlife habitat, water quality, recreation and carbon sequestration. Un-managed surface and groundwater exploitation in these areas usually leads to deterioration of such sensitive ecosystems by means of water resources degradation and/or increased salinity. Groundwater usually plays a vital role for the sustainability of these hydrological systems, as the underlying aquifers operate as regulators for both quantity and quality of their waters. Multi-layer and multi-objective Managed Aquifer Recharge (MAR) systems can be proved effective groundwater engineered solutions for the restoration of deteriorated coastal wetlands in semi- and arid regions. The plain of Marathon is a typical Mediterranean environment that hosts a naturally occurring -and today degraded- coastal wetland with the characteristics of a distinct ecosystem linked to a typical coastal hydrogeological system of a semi-arid region; and therefore can serve as a model for similar systems world-wide. The geo-hydrological setting of the area involves a multi-layer aquifer system consisting of (i) an upper un-consolidated formation of depositional unit dominated mostly by fluvial sediments and (ii) the surrounding and underlying karstified marbles; both being linked to the investigated wetland and also subjected to seawater encroachment. A smart engineered MAR system via an optimised Pump & Treat system integrated with an Aquifer Storage and Recovery (ASR) scheme in this area would include the abstraction of brackish groundwater from the deeper karst aquifer at a location close to the shoreline and direct treatment with Reverse Osmosis (RO). for desalination. Two-fold re-use scheme of the purified effluent can then be engineered for (i) the restoration of the coastal wetland; and (ii

  1. Implicit solution of the material transport in Stokes flow simulation: Toward thermal convection simulation surrounded by free surface

    NASA Astrophysics Data System (ADS)

    Furuichi, Mikito; May, Dave A.

    2015-07-01

    We present implicit time integration schemes suitable for modeling free surface Stokes flow dynamics with marker in cell (MIC) based spatial discretization. Our target is for example thermal convection surrounded by deformable surface boundaries to simulate the long term planetary formation process. The numerical system becomes stiff when the dynamical balancing time scale for the increasing/decreasing load by surface deformation is very short compared with the time scale associated with thermal convection. Any explicit time integration scheme will require very small time steps; otherwise, serious numerical oscillation (spurious solutions) will occur. The implicit time integration scheme possesses a wider stability region than the explicit method; therefore, it is suitable for stiff problems. To investigate an efficient solution method for the stiff Stokes flow system, we apply first (backward Euler (BE)) and second order (trapezoidal method (TR) and trapezoidal rule-backward difference formula (TR-BDF2)) accurate implicit methods for the MIC solution scheme. The introduction of implicit time integration schemes results in nonlinear systems of equations. We utilize a Jacobian free Newton Krylov (JFNK) based Newton framework to solve the resulting nonlinear equations. In this work we also investigate two efficient implicit solution strategies to reduce the computational cost when solving stiff nonlinear systems. The two methods differ in how the advective term in the material transport evolution equation is treated. We refer to the method that employs Lagrangian update as "fully implicit" (Imp), whilst the method that employs Eulerian update is referred to as "semi-implicit" (SImp). Using a finite difference (FD) method, we have performed a series of numerical experiments which clarify the accuracy of solutions and trade-off between the computational cost associated with the nonlinear solver and time step size. In comparison with the general explicit Euler method

  2. Altered Implicit Category Learning in Anorexia Nervosa

    PubMed Central

    Shott, Megan E.; Filoteo, J. Vincent; Jappe, Leah M.; Pryor, Tamara; Maddox, W. Todd; Rollin, Michael D.H.; Hagman, Jennifer O.; Frank, Guido K.W.

    2012-01-01

    Objective Recent research has identified specific cognitive deficits in patients with anorexia nervosa (AN), including impairment in executive functioning and attention. Another such cognitive process, implicit category learning has been less studied in AN. This study examined whether implicit category learning is impaired in AN. Method Twenty-one women diagnosed with AN and 19 control women (CW) were administered an implicit category learning task in which they were asked to categorize simple perceptual stimuli (Gabor patches) into one of two categories. Category membership was based on a linear integration (i.e., an implicit task) of two stimulus dimensions (orientation and spatial frequency of the stimulus). Results AN individuals were less accurate on implicit category learning relative to age-matched CW. Model-based analyses indicated that, even when AN individuals used the appropriate (i.e., implicit) strategy they were still impaired relative to CW who also used the same strategy. In addition, task performance in AN patients was worse the higher they were in self-reported novelty seeking and the lower they were in sensitivity to punishment. Conclusions These results indicate that AN patients have implicit category learning deficits, and given this type of learning is thought to be mediated by striatal dopamine pathways, AN patients may have deficits in these neural systems. The finding of significant correlations with novelty seeking and sensitivity to punishment suggests that feedback sensitivity is related to implicit learning in AN. PMID:22201300

  3. Security Analysis and Improvement of 'a More Secure Anonymous User Authentication Scheme for the Integrated EPR Information System'.

    PubMed

    Islam, S K Hafizul; Khan, Muhammad Khurram; Li, Xiong

    2015-01-01

    Over the past few years, secure and privacy-preserving user authentication scheme has become an integral part of the applications of the healthcare systems. Recently, Wen has designed an improved user authentication system over the Lee et al.'s scheme for integrated electronic patient record (EPR) information system, which has been analyzed in this study. We have found that Wen's scheme still has the following inefficiencies: (1) the correctness of identity and password are not verified during the login and password change phases; (2) it is vulnerable to impersonation attack and privileged-insider attack; (3) it is designed without the revocation of lost/stolen smart card; (4) the explicit key confirmation and the no key control properties are absent, and (5) user cannot update his/her password without the help of server and secure channel. Then we aimed to propose an enhanced two-factor user authentication system based on the intractable assumption of the quadratic residue problem (QRP) in the multiplicative group. Our scheme bears more securities and functionalities than other schemes found in the literature. PMID:26263401

  4. Security Analysis and Improvement of ‘a More Secure Anonymous User Authentication Scheme for the Integrated EPR Information System’

    PubMed Central

    Islam, SK Hafizul; Khan, Muhammad Khurram; Li, Xiong

    2015-01-01

    Over the past few years, secure and privacy-preserving user authentication scheme has become an integral part of the applications of the healthcare systems. Recently, Wen has designed an improved user authentication system over the Lee et al.’s scheme for integrated electronic patient record (EPR) information system, which has been analyzed in this study. We have found that Wen’s scheme still has the following inefficiencies: (1) the correctness of identity and password are not verified during the login and password change phases; (2) it is vulnerable to impersonation attack and privileged-insider attack; (3) it is designed without the revocation of lost/stolen smart card; (4) the explicit key confirmation and the no key control properties are absent, and (5) user cannot update his/her password without the help of server and secure channel. Then we aimed to propose an enhanced two-factor user authentication system based on the intractable assumption of the quadratic residue problem (QRP) in the multiplicative group. Our scheme bears more securities and functionalities than other schemes found in the literature. PMID:26263401

  5. Community-based health insurance programmes and the national health insurance scheme of Nigeria: challenges to uptake and integration

    PubMed Central

    2014-01-01

    Background Nigeria has included a regulated community-based health insurance (CBHI) model within its National Health Insurance Scheme (NHIS). Uptake to date has been disappointing, however. The aim of this study is to review the present status of CBHI in SSA in general to highlight the issues that affect its successful integration within the NHIS of Nigeria and more widely in developing countries. Methods A literature survey using PubMed and EconLit was carried out to identify and review studies that report factors affecting implementation of CBHI in SSA with a focus on Nigeria. Results CBHI schemes with a variety of designs have been introduced across SSA but with generally disappointing results so far. Two exceptions are Ghana and Rwanda, both of which have introduced schemes with effective government control and support coupled with intensive implementation programmes. Poor support for CBHI is repeatedly linked elsewhere with failure to engage and account for the ‘real world’ needs of beneficiaries, lack of clear legislative and regulatory frameworks, inadequate financial support, and unrealistic enrolment requirements. Nigeria’s CBHI-type schemes for the informal sectors of its NHIS have been set up under an appropriate legislative framework, but work is needed to eliminate regressive financing, to involve scheme members in the setting up and management of programmes, to inform and educate more effectively, to eliminate lack of confidence in the schemes, and to address inequity in provision. Targeted subsidies should also be considered. Conclusions Disappointing uptake of CBHI-type NHIS elements in Nigeria can be addressed through closer integration of informal and formal programmes under the NHIS umbrella, with increasing involvement of beneficiaries in scheme design and management, improved communication and education, and targeted financial assistance. PMID:24559409

  6. Implicit and Multigrid Method for Ideal Multigrid Convergence: Direct Numerical Simulation of Separated Flow Around NACA 0012 Airfoil

    NASA Technical Reports Server (NTRS)

    Liu, Chao-Qun; Shan, H.; Jiang, L.

    1999-01-01

    Numerical investigation of flow separation over a NACA 0012 airfoil at large angles of attack has been carried out. The numerical calculation is performed by solving the full Navier-Stokes equations in generalized curvilinear coordinates. The second-order LU-SGS implicit scheme is applied for time integration. This scheme requires no tridiagonal inversion and is capable of being completely vectorized, provided the corresponding Jacobian matrices are properly selected. A fourth-order centered compact scheme is used for spatial derivatives. In order to reduce numerical oscillation, a sixth-order implicit filter is employed. Non-reflecting boundary conditions are imposed at the far-field and outlet boundaries to avoid possible non-physical wave reflection. Complex flow separation and vortex shedding phenomenon have been observed and discussed.

  7. Unconditionally stable time marching scheme for Reynolds stress models

    NASA Astrophysics Data System (ADS)

    Mor-Yossef, Y.

    2014-11-01

    Progress toward a stable and efficient numerical treatment for the compressible Favre-Reynolds-averaged Navier-Stokes equations with a Reynolds-stress model (RSM) is presented. The mean-flow and the Reynolds stress model equations are discretized using finite differences on a curvilinear coordinates mesh. The convective flux is approximated by a third-order upwind biased MUSCL scheme. The diffusive flux is approximated using second-order central differencing, based on a full-viscous stencil. The novel time-marching approach relies on decoupled, implicit time integration, that is, the five mean-flow equations are solved separately from the seven Reynolds-stress closure equations. The key idea is the use of the unconditionally positive-convergent implicit scheme (UPC), originally developed for two-equation turbulence models. The extension of the UPC scheme for RSM guarantees the positivity of the normal Reynolds-stress components and the turbulence (specific) dissipation rate for any time step. Thanks to the UPC matrix-free structure and the decoupled approach, the resulting computational scheme is very efficient. Special care is dedicated to maintain the implicit operator compact, involving only nearest neighbor grid points, while fully supporting the larger discretized residual stencil. Results obtained from two- and three-dimensional numerical simulations demonstrate the significant progress achieved in this work toward optimally convergent solution of Reynolds stress models. Furthermore, the scheme is shown to be unconditionally stable and positive.

  8. A robust data fusion scheme for integrated navigation systems employing fault detection methodology augmented with fuzzy adaptive filtering

    NASA Astrophysics Data System (ADS)

    Ushaq, Muhammad; Fang, Jiancheng

    2013-10-01

    Integrated navigation systems for various applications, generally employs the centralized Kalman filter (CKF) wherein all measured sensor data are communicated to a single central Kalman filter. The advantage of CKF is that there is a minimal loss of information and high precision under benign conditions. But CKF may suffer computational overloading, and poor fault tolerance. The alternative is the federated Kalman filter (FKF) wherein the local estimates can deliver optimal or suboptimal state estimate as per certain information fusion criterion. FKF has enhanced throughput and multiple level fault detection capability. The Standard CKF or FKF require that the system noise and the measurement noise are zero-mean and Gaussian. Moreover it is assumed that covariance of system and measurement noises remain constant. But if the theoretical and actual statistical features employed in Kalman filter are not compatible, the Kalman filter does not render satisfactory solutions and divergence problems also occur. To resolve such problems, in this paper, an adaptive Kalman filter scheme strengthened with fuzzy inference system (FIS) is employed to adapt the statistical features of contributing sensors, online, in the light of real system dynamics and varying measurement noises. The excessive faults are detected and isolated by employing Chi Square test method. As a case study, the presented scheme has been implemented on Strapdown Inertial Navigation System (SINS) integrated with the Celestial Navigation System (CNS), GPS and Doppler radar using FKF. Collectively the overall system can be termed as SINS/CNS/GPS/Doppler integrated navigation system. The simulation results have validated the effectiveness of the presented scheme with significantly enhanced precision, reliability and fault tolerance. Effectiveness of the scheme has been tested against simulated abnormal errors/noises during different time segments of flight. It is believed that the presented scheme can be

  9. Implementation of 5-layer thermal diffusion scheme in weather research and forecasting model with Intel Many Integrated Cores

    NASA Astrophysics Data System (ADS)

    Huang, Melin; Huang, Bormin; Huang, Allen H.

    2014-10-01

    For weather forecasting and research, the Weather Research and Forecasting (WRF) model has been developed, consisting of several components such as dynamic solvers and physical simulation modules. WRF includes several Land- Surface Models (LSMs). The LSMs use atmospheric information, the radiative and precipitation forcing from the surface layer scheme, the radiation scheme, and the microphysics/convective scheme all together with the land's state variables and land-surface properties, to provide heat and moisture fluxes over land and sea-ice points. The WRF 5-layer thermal diffusion simulation is an LSM based on the MM5 5-layer soil temperature model with an energy budget that includes radiation, sensible, and latent heat flux. The WRF LSMs are very suitable for massively parallel computation as there are no interactions among horizontal grid points. The features, efficient parallelization and vectorization essentials, of Intel Many Integrated Core (MIC) architecture allow us to optimize this WRF 5-layer thermal diffusion scheme. In this work, we present the results of the computing performance on this scheme with Intel MIC architecture. Our results show that the MIC-based optimization improved the performance of the first version of multi-threaded code on Xeon Phi 5110P by a factor of 2.1x. Accordingly, the same CPU-based optimizations improved the performance on Intel Xeon E5- 2603 by a factor of 1.6x as compared to the first version of multi-threaded code.

  10. Canopy reflectance, photosynthesis, and transpiration. III - A reanalysis using improved leaf models and a new canopy integration scheme

    NASA Technical Reports Server (NTRS)

    Sellers, P. J.; Berry, J. A.; Collatz, G. J.; Field, C. B.; Hall, F. G.

    1992-01-01

    The theoretical analyses of Sellers (1985, 1987), which linked canopy spectral reflectance properties to (unstressed) photosynthetic rates and conductances, are critically reviewed and significant shortcomings are identified. These are addressed in this article principally through the incorporation of a more sophisticated and realistic treatment of leaf physiological processes within a new canopy integration scheme. The results indicate that area-averaged spectral vegetation indices, as obtained from coarse resolution satellite sensors, may give good estimates of the area-integrals of photosynthesis and conductance even for spatially heterogenous (though physiologically uniform) vegetation covers.

  11. Time integration algorithms for the two-dimensional Euler equations on unstructured meshes

    NASA Technical Reports Server (NTRS)

    Slack, David C.; Whitaker, D. L.; Walters, Robert W.

    1994-01-01

    Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.

  12. Awareness of Implicit Attitudes

    PubMed Central

    Hahn, Adam; Judd, Charles M.; Hirsh, Holen K.; Blair, Irene V.

    2013-01-01

    Research on implicit attitudes has raised questions about how well people know their own attitudes. Most research on this question has focused on the correspondence between measures of implicit attitudes and measures of explicit attitudes, with low correspondence interpreted as showing that people have little awareness of their implicit attitudes. We took a different approach and directly asked participants to predict their results on upcoming IAT measures of implicit attitudes toward five different social groups. We found that participants were surprisingly accurate in their predictions. Across four studies, predictions were accurate regardless of whether implicit attitudes were described as true attitudes or culturally learned associations (Studies 1 and 2), regardless of whether predictions were made as specific response patterns (Study 1) or as conceptual responses (Studies 2–4), and regardless of how much experience or explanation participants received before making their predictions (Study 4). Study 3 further suggested that participants’ predictions reflected unique insight into their own implicit responses, beyond intuitions about how people in general might respond. Prediction accuracy occurred despite generally low correspondence between implicit and explicit measures of attitudes, as found in prior research. All together, the research findings cast doubt on the belief that attitudes or evaluations measured by the IAT necessarily reflect unconscious attitudes. PMID:24294868

  13. Assimilation of stream discharge for flood forecasting: Updating a semidistributed model with an integrated data assimilation scheme

    NASA Astrophysics Data System (ADS)

    Li, Yuan; Ryu, Dongryeol; Western, Andrew W.; Wang, Q. J.

    2015-05-01

    Real-time discharge observations can be assimilated into flood models to improve forecast accuracy; however, the presence of time lags in the routing process and a lack of methods to quantitatively represent different sources of uncertainties challenge the implementation of data assimilation techniques for operational flood forecasting. To address these issues, an integrated error parameter estimation and lag-aware data assimilation (IEELA) scheme was recently developed for a lumped model. The scheme combines an ensemble-based maximum a posteriori (MAP) error estimation approach with a lag-aware ensemble Kalman smoother (EnKS). In this study, the IEELA scheme is extended to a semidistributed model to provide for more general application in flood forecasting by including spatial and temporal correlations in model uncertainties between subcatchments. The result reveals that using a semidistributed model leads to more accurate forecasts than a lumped model in an open-loop scenario. The IEELA scheme improves the forecast accuracy significantly in both lumped and semidistributed models, and the superiority of the semidistributed model remains in the data assimilation scenario. However, the improvements resulting from IEELA are confined to the outlet of the catchment where the discharge observations are assimilated. Forecasts at "ungauged" internal locations are not improved, and in some instances, even become less accurate.

  14. Bidirectional fiber-wireless and fiber-IVLLC integrated system based on polarization-orthogonal modulation scheme.

    PubMed

    Lu, Hai-Han; Li, Chung-Yi; Chen, Hwan-Wei; Ho, Chun-Ming; Cheng, Ming-Te; Huang, Sheng-Jhe; Yang, Zih-Yi; Lin, Xin-Yao

    2016-07-25

    A bidirectional fiber-wireless and fiber-invisible laser light communication (IVLLC) integrated system that employs polarization-orthogonal modulation scheme for hybrid cable television (CATV)/microwave (MW)/millimeter-wave (MMW)/baseband (BB) signal transmission is proposed and demonstrated. To our knowledge, it is the first one that adopts a polarization-orthogonal modulation scheme in a bidirectional fiber-wireless and fiber-IVLLC integrated system with hybrid CATV/MW/MMW/BB signal. For downlink transmission, carrier-to-noise ratio (CNR), composite second-order (CSO), composite triple-beat (CTB), and bit error rate (BER) perform well over 40-km single-mode fiber (SMF) and 10-m RF/50-m optical wireless transport scenarios. For uplink transmission, good BER performance is obtained over 40-km SMF and 50-m optical wireless transport scenario. Such a bidirectional fiber-wireless and fiber-IVLLC integrated system for hybrid CATV/MW/MMW/BB signal transmission will be an attractive alternative for providing broadband integrated services, including CATV, Internet, and telecommunication services. It is shown to be a prominent one to present the advancements for the convergence of fiber backbone and RF/optical wireless feeder. PMID:27464174

  15. Fiber-optic integration and efficient detection schemes for optomechanical resonators

    NASA Astrophysics Data System (ADS)

    Cohen, Justin D.

    With the advent of the laser in the year 1960, the field of optics experienced a renaissance from what was considered to be a dull, solved subject to an active area of development, with applications and discoveries which are yet to be exhausted 55 years later. Light is now nearly ubiquitous not only in cutting-edge research in physics, chemistry, and biology, but also in modern technology and infrastructure. One quality of light, that of the imparted radiation pressure force upon reflection from an object, has attracted intense interest from researchers seeking to precisely monitor and control the motional degrees of freedom of an object using light. These optomechanical interactions have inspired myriad proposals, ranging from quantum memories and transducers in quantum information networks to precision metrology of classical forces. Alongside advances in micro- and nano-fabrication, the burgeoning field of optomechanics has yielded a class of highly engineered systems designed to produce strong interactions between light and motion. Optomechanical crystals are one such system in which the patterning of periodic holes in thin dielectric films traps both light and sound waves to a micro-scale volume. These devices feature strong radiation pressure coupling between high-quality optical cavity modes and internal nanomechanical resonances. Whether for applications in the quantum or classical domain, the utility of optomechanical crystals hinges on the degree to which light radiating from the device, having interacted with mechanical motion, can be collected and detected in an experimental apparatus consisting of conventional optical components such as lenses and optical fibers. While several efficient methods of optical coupling exist to meet this task, most are unsuitable for the cryogenic or vacuum integration required for many applications. The first portion of this dissertation will detail the development of robust and efficient methods of optically coupling

  16. A parallel algorithm for implicit depletant simulations

    NASA Astrophysics Data System (ADS)

    Glaser, Jens; Karas, Andrew S.; Glotzer, Sharon C.

    2015-11-01

    We present an algorithm to simulate the many-body depletion interaction between anisotropic colloids in an implicit way, integrating out the degrees of freedom of the depletants, which we treat as an ideal gas. Because the depletant particles are statistically independent and the depletion interaction is short-ranged, depletants are randomly inserted in parallel into the excluded volume surrounding a single translated and/or rotated colloid. A configurational bias scheme is used to enhance the acceptance rate. The method is validated and benchmarked both on multi-core processors and graphics processing units for the case of hard spheres, hemispheres, and discoids. With depletants, we report novel cluster phases in which hemispheres first assemble into spheres, which then form ordered hcp/fcc lattices. The method is significantly faster than any method without cluster moves and that tracks depletants explicitly, for systems of colloid packing fraction ϕc < 0.50, and additionally enables simulation of the fluid-solid transition.

  17. The versatile integrator of surface atmospheric processes - Part 2: evaluation of three topography-based runoff schemes

    NASA Astrophysics Data System (ADS)

    Niu, Guo-Yue; Yang, Zong-Liang

    2003-07-01

    Three different schemes of topography-based runoff production [versatile integrator of surface atmospheric processes (VISA)-TOP1, VISA-TOP2, and VISA-TOP3] are described for a land-surface model (LSM) developed for use with a general circulation model (GCM). The schemes' sensitivities to some key parameters are assessed for two catchments using data sets developed for the Project for Intercomparison of Land-Surface Parameterization Schemes (PILPS) Phase 2e. VISA-TOP1 differs from VISA-TOP2 only in how to treat oversaturated soil water from the soil layers. In VISA-TOP1, the oversaturated soil water is thrown out of the soil column; hence, it no longer plays a role in the ensuing soil water budgets. In VISA-TOP2, this oversaturated soil water is recharged back to the unsaturated soil layers above the water table; hence, it continues to involve in the water budgets. Unlike VISA-TOP1 and VISA-TOP2, VISA-TOP3 relaxes its dependence on the topographic parameters. The oversaturated soil water is treated the same in both VISA-TOP2 and VISA-TOP3. All three models reproduce the daily and seasonal cycles of streamflow provided that different values of the saturated hydraulic conductivity decay factor are used. The decay factor controls the timing and partitioning of subsurface runoff. In both VISA-TOP1 and VISA-TOP2, an anisotropic parameter explaining different hydraulic conductivities in the vertical and horizontal directions is critical for using the topographic index in the land-surface model. In the VISA-TOP2 scheme, the topography-controlled subsurface runoff is dominant because the oversaturated water is recharged to upper unsaturated soil layers to raise the water table. The water budgets in all these schemes show dramatically different responses to the decay factor, indicating that the calibrated parameters and the model formulations should not be separated.

  18. High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

    NASA Technical Reports Server (NTRS)

    Nielsen, Tanner B.; Carpenter, Mark H.; Fisher, Travis C.; Frankel, Steven H.

    2014-01-01

    This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.

  19. Parallel Implicit Kinetic Simulation with PARSEK

    NASA Astrophysics Data System (ADS)

    Stefano, Markidis; Giovanni, Lapenta

    2004-11-01

    Kinetic plasma simulation is the ultimate tool for plasma analysis. One of the prime tools for kinetic simulation is the particle in cell (PIC) method. The explicit or semi-implicit (i.e. implicit only on the fields) PIC method requires exceedingly small time steps and grid spacing, limited by the necessity to resolve the electron plasma frequency, the Debye length and the speed of light (for fully explicit schemes). A different approach is to consider fully implicit PIC methods where both particles and fields are discretized implicitly. This approach allows radically larger time steps and grid spacing, reducing the cost of a simulation by orders of magnitude while keeping the full kinetic treatment. In our previous work, simulations impossible for the explicit PIC method even on massively parallel computers have been made possible on a single processor machine using the implicit PIC code CELESTE3D [1]. We propose here another quantum leap: PARSEK, a parallel cousin of CELESTE3D, based on the same approach but sporting a radically redesigned software architecture (object oriented C++, where CELESTE3D was structured and written in FORTRAN77/90) and fully parallelized using MPI for both particle and grid communication. [1] G. Lapenta, J.U. Brackbill, W.S. Daughton, Phys. Plasmas, 10, 1577 (2003).

  20. Integration of electromagnetic induction sensor data in soil sampling scheme optimization using simulated annealing.

    PubMed

    Barca, E; Castrignanò, A; Buttafuoco, G; De Benedetto, D; Passarella, G

    2015-07-01

    Soil survey is generally time-consuming, labor-intensive, and costly. Optimization of sampling scheme allows one to reduce the number of sampling points without decreasing or even increasing the accuracy of investigated attribute. Maps of bulk soil electrical conductivity (EC a ) recorded with electromagnetic induction (EMI) sensors could be effectively used to direct soil sampling design for assessing spatial variability of soil moisture. A protocol, using a field-scale bulk EC a survey, has been applied in an agricultural field in Apulia region (southeastern Italy). Spatial simulated annealing was used as a method to optimize spatial soil sampling scheme taking into account sampling constraints, field boundaries, and preliminary observations. Three optimization criteria were used. the first criterion (minimization of mean of the shortest distances, MMSD) optimizes the spreading of the point observations over the entire field by minimizing the expectation of the distance between an arbitrarily chosen point and its nearest observation; the second criterion (minimization of weighted mean of the shortest distances, MWMSD) is a weighted version of the MMSD, which uses the digital gradient of the grid EC a data as weighting function; and the third criterion (mean of average ordinary kriging variance, MAOKV) minimizes mean kriging estimation variance of the target variable. The last criterion utilizes the variogram model of soil water content estimated in a previous trial. The procedures, or a combination of them, were tested and compared in a real case. Simulated annealing was implemented by the software MSANOS able to define or redesign any sampling scheme by increasing or decreasing the original sampling locations. The output consists of the computed sampling scheme, the convergence time, and the cooling law, which can be an invaluable support to the process of sampling design. The proposed approach has found the optimal solution in a reasonable computation time. The

  1. Electromagnetic direct implicit PIC simulation

    SciTech Connect

    Langdon, A.B.

    1983-03-29

    Interesting modelling of intense electron flow has been done with implicit particle-in-cell simulation codes. In this report, the direct implicit PIC simulation approach is applied to simulations that include full electromagnetic fields. The resulting algorithm offers advantages relative to moment implicit electromagnetic algorithms and may help in our quest for robust and simpler implicit codes.

  2. A family of time-staggered schemes for integrating hybrid DPD models for polymers: Algorithms and applications

    NASA Astrophysics Data System (ADS)

    Symeonidis, Vasileios; Karniadakis, George Em

    2006-10-01

    We propose new schemes for integrating the stochastic differential equations of dissipative particle dynamics (DPD) in simulations of dilute polymer solutions. The hybrid DPD models consist of hard potentials that describe the microscopic dynamics of polymers and soft potentials that describe the mesoscopic dynamics of the solvent. In particular, we develop extensions to the velocity-Verlet and Lowe's approaches - two representative DPD time-integrators - following a subcycling procedure whereby the solvent is advanced with a timestep much larger than the one employed in the polymer time-integration. The introduction of relaxation parameters allows optimization studies for accuracy while maintaining the low computational complexity of standard DPD algorithms. We demonstrate through equilibrium simulations that a 10-fold gain in efficiency can be obtained with the time-staggered algorithms without loss of accuracy compared to the non-staggered schemes. We then apply the new approach to investigate the scaling response of polymers in equilibrium as well as the dynamics of λ-phage DNA molecules subjected to shear.

  3. Orbit-averaged implicit particle codes

    NASA Astrophysics Data System (ADS)

    Cohen, B. I.; Freis, R. P.; Thomas, V.

    1982-03-01

    The merging of orbit-averaged particle code techniques with recently developed implicit methods to perform numerically stable and accurate particle simulations are reported. Implicitness and orbit averaging can extend the applicability of particle codes to the simulation of long time-scale plasma physics phenomena by relaxing time-step and statistical constraints. Difference equations for an electrostatic model are presented, and analyses of the numerical stability of each scheme are given. Simulation examples are presented for a one-dimensional electrostatic model. Schemes are constructed that are stable at large-time step, require fewer particles, and, hence, reduce input-output and memory requirements. Orbit averaging, however, in the unmagnetized electrostatic models tested so far is not as successful as in cases where there is a magnetic field. Methods are suggested in which orbit averaging should achieve more significant improvements in code efficiency.

  4. A transient FETI methodology for large-scale parallel implicit computations in structural mechanics, part 2

    NASA Technical Reports Server (NTRS)

    Farhat, Charbel; Crivelli, Luis

    1993-01-01

    Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallellize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet and perhaps will never be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient than explicit codes when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

  5. A transient FETI methodology for large-scale parallel implicit computations in structural mechanics

    NASA Technical Reports Server (NTRS)

    Farhat, Charbel; Crivelli, Luis; Roux, Francois-Xavier

    1992-01-01

    Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallelize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet -- and perhaps will never -- be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

  6. Implicit solvers for unstructured meshes

    NASA Technical Reports Server (NTRS)

    Venkatakrishnan, V.; Mavriplis, Dimitri J.

    1991-01-01

    Implicit methods were developed and tested for unstructured mesh computations. The approximate system which arises from the Newton linearization of the nonlinear evolution operator is solved by using the preconditioned GMRES (Generalized Minimum Residual) technique. Three different preconditioners were studied, namely, the incomplete LU factorization (ILU), block diagonal factorization, and the symmetric successive over relaxation (SSOR). The preconditioners were optimized to have good vectorization properties. SSOR and ILU were also studied as iterative schemes. The various methods are compared over a wide range of problems. Ordering of the unknowns, which affects the convergence of these sparse matrix iterative methods, is also studied. Results are presented for inviscid and turbulent viscous calculations on single and multielement airfoil configurations using globally and adaptively generated meshes.

  7. Implementation of the Kirchhoff integral for elastic waves in staggered-grid modeling schemes

    SciTech Connect

    Mittet, R. )

    1994-12-01

    Implementation of boundary conditions in finite-difference schemes is not straightforward for the elastic wave equation if a staggered grid formulation is used. Reverse time migration of VSP data requires a proper description of the recording surface so as not to excite false P- and S-waves. Such contributions may cause artifacts in the imaging procedure. The boundary conditions for the elastic stress tensor can be implemented numerically in a staggered coarse grid modeling scheme by using band-limited spatial delta-functions and band-limited first-order derivatives of these spatial delta-functions. A representation theorem for elastic waves is derived to test the implementation of the spatial part of the boundary condition. The implementation is tested in a 2-D numerical experiment for a closed, but curved, boundary S enclosing a volume V. The test condition is that within the volume V, the difference between the forward modeled field and the retropropagated field should be equal to zero. Both P- and S-waves are properly recovered in a 2-D reverse time modeling example. The numerical artifacts related to the proposed spatial approximation of the boundary condition are found to be negligible.

  8. On the Dynamics of Implicit Linear Multistep Methods

    NASA Technical Reports Server (NTRS)

    Yee, H. C.; Sweby, P. K.; Rai, Man Mohan (Technical Monitor)

    1995-01-01

    Some new guidelines on the usage of implicit linear multistep methods (LMMs) as time-dependent approaches for obtaining steady-state numerical solutions in computational fluid dynamics (CFD) are explored. The commonly used implicit LMMs in CFD belong to the class of superstable time discretizations. It can be shown that the nonlinear asymptotic behavior in terms of bifurcation diagrams and basins of attractions of these schemes can provide an improved range of initial data and time step over the linearized stability limit.

  9. Implicit Learning of Local Context in Autism Spectrum Disorder

    ERIC Educational Resources Information Center

    Kourkoulou, Anastasia; Leekam, Susan R.; Findlay, John M.

    2012-01-01

    Although previous research has reported impairments in implicit learning in individuals with ASD, research using one implicit learning paradigm, the contextual cueing task (Chun and Jiang in Cognitive Psychol 36:28-71, 1998), shows evidence of intact ability to integrate spatial contextual information. Using an adaptation of this paradigm, we…

  10. Case study of developing an integrated water and nitrogen scheme for agricultural systems on the North China Plain

    NASA Astrophysics Data System (ADS)

    Liu, Y.; Tao, F.; Luo, Y.; Ma, J.

    2013-12-01

    Appropriate irrigation and nitrogen fertilization, along with suitable crop management strategies, are essential prerequisites for optimum yields in agricultural systems. This research attempts to provide a scientific basis for sustainable agricultural production management for the North China Plain and other semi-arid regions. Based on a series of 72 treatments over 2003-2008, an optimized water and nitrogen scheme for winter wheat/summer maize cropping system was developed. Integrated systems incorporating 120 mm of water with 80 kg N ha-1 N fertilizer were used to simulate winter wheat yields in Hebei and 120 mm of water with 120 kg N ha-1 were used to simulate winter wheat yields in Shandong and Henan provinces in 2000-2007. Similarly, integrated treatments of 40 kg N ha-1 N fertilizer were used to simulate summer maize yields in Hebei, and 80 kg N ha-1 was used to simulate summer maize yields in Shandong and Henan provinces in 2000-2007. Under the optimized scheme, 341.74 107 mm ha-1 of water and 575.79 104 Mg of urea fertilizer could be saved per year under the wheat/maize rotation system. Despite slight drops in the yields of wheat and maize in some areas, water and fertilizer saving has tremendous long-term eco-environmental benefits.

  11. Parallel implicit unstructured grid Euler solvers

    NASA Technical Reports Server (NTRS)

    Venkatakrishnan, V.

    1994-01-01

    A mesh-vertex finite volume scheme for solving the Euler equations on triangular unstructured meshes is implemented on a multiple-instruction/multiple-data stream parallel computer. An explicit four-stage Runge-Kutta scheme is used to solve two-dimensional flow problems. A family of implicit schemes is also developed to solve these problems, where the linear system that arises at each time step is solved by a preconditioned GMRES algorithm. Two partitioning strategies are employed: one that partitions triangles and the other that partitions vertices. The choice of the preconditioner in a distributed memory setting is discussed. All of the methods are compared both in terms of elapsed times and convergence rates. It is shown that the implicit schemes offer adequate parallelism at the expense of minimal sequential overhead. The use of a global coarse grid to further minimize this overhead is also investigated. The schemes are implemented on a distributed memory parallel computer, the Intel iPSC/860.

  12. Parallel implicit unstructured grid Euler solvers

    NASA Technical Reports Server (NTRS)

    Venkatakrishnan, V.

    1994-01-01

    A mesh-vertex finite volume scheme for solving the Euler equations on triangular unstructured meshes is implemented on an MIMD (multiple instruction/multiple data stream) parallel computer. An explicit four-stage Runge-Kutta scheme is used to solve two-dimensional flow problems. A family of implicit schemes is also developed to solve these problems, where the linear system that arises at each time step is solved by a preconditioned GMRES algorithm. Two partitioning strategies are employed, one that partitions triangles and the other that partitions vertices. The choice of the preconditioner in a distributed memory setting is discussed. All the methods are compared both in terms of elapsed times and convergence rates. It is shown that the implicit schemes offer adequate parallelism at the expense of minimal sequential overhead. The use of a global coarse grid to further minimize this overhead is also investigated. The schemes are implemented on a distributed memory parallel computer, the iPSC/860.

  13. Minimally implicit Runge-Kutta methods for Resistive Relativistic MHD

    NASA Astrophysics Data System (ADS)

    Aloy, Miguel-Á.; Cordero-Carrión, Isabel

    2016-05-01

    The Relativistic Resistive Magnetohydrodynamic (RRMHD) equations are a hyperbolic system of partial differential equations used to describe the dynamics of relativistic magnetized fluids with a finite conductivity. Close to the ideal magnetohydrodynamic regime, the source term proportional to the conductivity becomes potentially stiff and cannot be handled with standard explicit time integration methods. We propose a new class of methods to deal with the stiffness fo the system, which we name Minimally Implicit Runge-Kutta methods. These methods avoid the development of numerical instabilities without increasing the computational costs in comparison with explicit methods, need no iterative extra loop in order to recover the primitive (physical) variables, the analytical inversion of the implicit operator is trivial and the several stages can actually be viewed as stages of explicit Runge-Kutta methods with an effective time-step. We test these methods with two different one-dimensional test beds in varied conductivity regimes, and show that our second-order schemes satisfy the theoretical expectations.

  14. Self-consistency based control scheme for magnetization dynamics

    SciTech Connect

    Albuquerque, G.; Miltat, J.; Thiaville, A.

    2001-06-01

    A numerical framework is presented for the solution of the Landau{endash}Lifshitz{endash}Gilbert equation of magnetization motion using a semi-implicit Crank{endash}Nicholson integration scheme. Along with the details of both space and time domain discretizations, we report on the development of a physically based self-consistency criterion that allows for a quantitative measurement of error in dynamic micromagnetic simulations. In essence, this criterion relies in recalculating from actual magnetization motion the imposed phenomenological damping constant. Test calculations were performed with special attention paid to the determination of suitable integration time steps. {copyright} 2001 American Institute of Physics.

  15. Implicit Attitudes in Prosopagnosia

    PubMed Central

    Knutson, Kristine M.; DeTucci, Karen A.; Grafman, Jordan

    2011-01-01

    We studied a male with acquired prosopagnosia using a battery of implicit association tests (IATs) to investigate whether observing faces varying by social category would activate the patient’s implicit social biases. We also asked him to categorize faces explicitly by race, gender, and political party. The patient, G.B., was marginally slower to categorize black compared to white faces. He showed congruency effects in the race and celebrity IATs, but not in the gender or political IATs. These results indicate that G.B. possesses an implicit social sensitivity to certain facial stimuli despite an inability to overtly recognize familiar faces. The results demonstrate that social biases can be retrieved based on facial stimuli via pathways bypassing the fusiform gyri. Thus the IAT effect can be added to the list of covert recognition effects found in prosopagnosia. PMID:21414330

  16. Implicit Spacecraft Gyro Calibration

    NASA Technical Reports Server (NTRS)

    Harman, Richard; Bar-Itzhack, Itzhack Y.

    2003-01-01

    This paper presents an implicit algorithm for spacecraft onboard instrument calibration, particularly to onboard gyro calibration. This work is an extension of previous work that was done where an explicit gyro calibration algorithm was applied to the AQUA spacecraft gyros. The algorithm presented in this paper was tested using simulated data and real data that were downloaded from the Microwave Anisotropy Probe (MAP) spacecraft. The calibration tests gave very good results. A comparison between the use of the implicit calibration algorithm used here with the explicit algorithm used for AQUA spacecraft indicates that both provide an excellent estimation of the gyro calibration parameters with similar accuracies.

  17. An advanced implicit solver for MHD

    NASA Astrophysics Data System (ADS)

    Udrea, Bogdan

    A new implicit algorithm has been developed for the solution of the time-dependent, viscous and resistive single fluid magnetohydrodynamic (MHD) equations. The algorithm is based on an approximate Riemann solver for the hyperbolic fluxes and central differencing applied on a staggered grid for the parabolic fluxes. The algorithm employs a locally aligned coordinate system that allows the solution to the Riemann problems to be solved in a natural direction, normal to cell interfaces. The result is an original scheme that is robust and reduces the complexity of the flux formulas. The evaluation of the parabolic fluxes is also implemented using a locally aligned coordinate system, this time on the staggered grid. The implicit formulation employed by WARP3 is a two level scheme that was applied for the first time to the single fluid MHD model. The flux Jacobians that appear in the implicit scheme are evaluated numerically. The linear system that results from the implicit discretization is solved using a robust symmetric Gauss-Seidel method. The code has an explicit mode capability so that implementation and test of new algorithms or new physics can be performed in this simpler mode. Last but not least the code was designed and written to run on parallel computers so that complex, high resolution runs can be per formed in hours rather than days. The code has been benchmarked against analytical and experimental gas dynamics and MHD results. The benchmarks consisted of one-dimensional Riemann problems and diffusion dominated problems, two-dimensional supersonic flow over a wedge, axisymmetric magnetoplasmadynamic (MPD) thruster simulation and three-dimensional supersonic flow over intersecting wedges and spheromak stability simulation. The code has been proven to be robust and the results of the simulations showed excellent agreement with analytical and experimental results. Parallel performance studies showed that the code performs as expected when run on parallel

  18. On the improvement of deconvolution with digitized data using a Runge-Kutta integration scheme

    NASA Technical Reports Server (NTRS)

    Houghton, J. R.; Townsend, M. A.; Packman, P. F.

    1977-01-01

    A relatively simple change in the treatment of the input function in numerical integration of high-order differential equations by Runge-Kutta methods provides substantial improvements in accuracy, particularly when the forcing function is in digitized form. The Runge-Kutta-Gill coefficients are modified to incorporate the changes; with pulse-type excitations, improvements on the order of 2 to 50 times greater accuracy are demonstrated.

  19. Implicit solvent methods for free energy estimation

    PubMed Central

    Decherchi, Sergio; Masetti, Matteo; Vyalov, Ivan; Rocchia, Walter

    2014-01-01

    Solvation is a fundamental contribution in many biological processes and especially in molecular binding. Its estimation can be performed by means of several computational approaches. The aim of this review is to give an overview of existing theories and methods to estimate solvent effects giving a specific focus on the category of implicit solvent models and their use in Molecular Dynamics. In many of these models, the solvent is considered as a continuum homogenous medium, while the solute can be represented at the atomic detail and at different levels of theory. Despite their degree of approximation, implicit methods are still widely employed due to their trade-off between accuracy and efficiency. Their derivation is rooted in the statistical mechanics and integral equations disciplines, some of the related details being provided here. Finally, methods that combine implicit solvent models and molecular dynamics simulation, are briefly described. PMID:25193298

  20. Emotion and Implicit Timing

    PubMed Central

    Droit-Volet, Sylvie

    2016-01-01

    This study examined the effects of emotion on implicit timing. In the implicit timing task used, the participants did not receive any temporal instructions. Instead they were simply asked and trained to press a key as quickly as possible after a stimulus (response stimulus) that was separated from a preceding stimulus by a given temporal interval (reference interval duration). However, in the testing phase, the interval duration was the reference interval duration or a shorter or longer interval duration. In addition, the participants attended two sessions: a first baseline session in which no stimulus was presented during the inter-stimulus intervals, and a second emotional session in which emotional facial expressions (angry, neutral and sad facial expressions) were presented during these intervals. Results showed faster RTs for interval durations close to the reference duration in both the baseline and the emotional conditions and yielded a U-shaped curve. This suggests that implicit processing of time persists in emotional contexts. In addition, the RT was faster for the facial expressions of anger than for those of neutrality and sadness. However, the U-shaped RT curve did not peak clearly at a shorter interval duration for the angry than for the other facial expressions. This lack of time distortion in an implicit timing task in response to arousing emotional stimuli questions the idea of an automatic speeding-up of the interval clock system involved in the representation of time. PMID:27380409

  1. Sexual Murderers' Implicit Theories

    ERIC Educational Resources Information Center

    Beech, Anthony; Fisher, Dawn; Ward, Tony

    2005-01-01

    Interviews with 28 sexual murderers were subjected to grounded theory analysis. Five implicit theories (ITs) were identified: dangerous world, male sex drive is uncontrollable, entitlement, women as sexual objects, and women as unknowable. These ITs were found to be identical to those identified in the literature as being present in rapists. The…

  2. Integration of polytetrafluoroethylene low-k dielectric material in a chemical vapor deposited aluminum metallization scheme

    NASA Astrophysics Data System (ADS)

    Wickland, Heidi Lee

    In order to increase speed and functionality, computer chip technology continues its evolution towards higher device density and reduced feature size. As interconnect signal delay has become prohibitively high for traditional Al/SiO2 architectures, there is a pressing need to pursue integration of low-k materials into back-end-of-the-line interconnect architectures. To this end, polytetrafluoroethylene (PTFE) was integrated with aluminum, implementing titanium nitride as a barrier layer. SPEEDFILM, a version of PTFE, has a low dielectric constant of 1.9--2.0. With the use of an adhesion promoter, SPEEDFILM can be easily spun-cast onto patterned aluminum substrates, allowing for implementation into existing interconnect process flows. CVD aluminum/PVD TiN binary stacks have been successfully grown on spin-cast films of PTFE. Elemental analysis revealed compositionally pure aluminum, showing virtually no fluorine contamination in as-deposited and annealed binary stacks, and resistivity as low as 3.2 muO-cm have been obtained. Thermal annealing studies indicate titanium nitride is an effective barrier against thermally-driven fluorine diffusion from PTFE. A fabrication flow for a two-level Al/PTFE test structure was identified in order to examine key integration issues and produce an electrically testable demonstration vehicle. Many integration milestones have been met. In particular, XPS studies of plasma-treated PTFE surfaces demonstrated that an NH3 plasma can defluorinate the PTFE surface, allowing sufficient adhesion of a subsequently-deposited SiO2 cap to survive chemical-mechanical polishing. First pass work also indicated promising etching and photoresist stripping results of SiO2/PTFE stacks, which are designed for use at the via level of this two-level structure. This is a particularly useful result since low-k dielectrics are often difficult to pattern due to their degradation upon exposure to traditional photoresist stripping chemistries. Finally, in

  3. Dynamics of spatial beams in quaternion description based on the Newmark integration scheme

    NASA Astrophysics Data System (ADS)

    Zupan, Eva; Saje, Miran; Zupan, Dejan

    2013-01-01

    The rotational quaternions represent a unique four dimensional parametrization of rotations in the three dimensional Euclidean space. In the present paper they are used as the basic rotational parameters in formulating the finite-element approach of geometrically exact beam-like structures. The classical concept of parameterizing the rotation matrix by the rotational vector is completely abandoned so that the only rotational parameters are the rotational quaternions representing both rotations and rotational strains in the beam. The space discretization based on the collocation method is used and the adjustment of the Newmark time-integration algorithm to the quaternion parameterizations of rotation is presented.

  4. Food-energy integrated development schemes in Brazil: FINEP's agro-energy communities program

    SciTech Connect

    La Rovere, E.L.; Baiardi, A.

    1983-12-01

    In Brazil social problems strongly affect the setting of development priorities. Effects of poverty and the need for imported petroleum products are discussed with reference to food and energy relationships. The object of the program is to design and test in real life conditions integrated food, energy and industrial raw materials production systems through the application of an appropriate technology specifically conceived in each case according to the participants of the ecosystem involved. The sophistication level for proposed solutions and the degree of participation of the local community are important considerations. Two projects have been started to date and two more are under consideration. 7 references.

  5. From Classical to Quantum and Back: A Hamiltonian Scheme for Adaptive Multiresolution Classical/Path-Integral Simulations.

    PubMed

    Kreis, Karsten; Tuckerman, Mark E; Donadio, Davide; Kremer, Kurt; Potestio, Raffaello

    2016-07-12

    Quantum delocalization of atomic nuclei affects the physical properties of many hydrogen-rich liquids and biological systems even at room temperature. In computer simulations, quantum nuclei can be modeled via the path-integral formulation of quantum statistical mechanics, which implies a substantial increase in computational overhead. By restricting the quantum description to a small spatial region, this cost can be significantly reduced. Herein, we derive a bottom-up, rigorous, Hamiltonian-based scheme that allows molecules to change from quantum to classical and vice versa on the fly as they diffuse through the system, both reducing overhead and making quantum grand-canonical simulations possible. The method is validated via simulations of low-temperature parahydrogen. Our adaptive resolution approach paves the way to efficient quantum simulations of biomolecules, membranes, and interfaces. PMID:27214610

  6. 40 Gb/s InGaAlAs EML module based on identical epitaxial layer integration scheme

    NASA Astrophysics Data System (ADS)

    Sun, Changzheng; Xiong, Bing; Wang, Jian; Cai, Pengfei; Xu, Jianming; Zhou, Qiwei; Yuan, He; Luo, Yi

    2008-11-01

    High-speed AlGaInAs multiple-quantum-well (MQW) electroabsorption modulated lasers (EMLs) based on identical epitaxial layer (IEL) integration scheme are developed for 40 Gb/s optical fiber communication systems. The electroabsorption modulator (EAM) section adopts a narrow high-mesa waveguide formed by inductively coupled plasma (ICP) dry etching technique, and a self-aligned planarization technique is employed to further reduce the device capacitance. Resonances are observed in the small signal modulation response of the packaged EML module, which are attributed to parallel-plate modes of the coplanar waveguide (CPW) transmission line used for modulation signal feeding and the residual reflection at the modulator facet, respectively. The influence of such resonances on the large signal eyediagram performance of the device is studied, and methods for their suppression are presented. Clear eye opening under 40 Gb/s non-return-to-zero (NRZ) modulation has been demonstrated for the optimized EML module.

  7. A membrane-integrated advanced scheme for treatment of industrial wastewater: dynamic modeling towards scale up.

    PubMed

    Kumar, Ramesh; Pal, Parimal

    2013-08-01

    Modeling and simulation was carried out for an advanced membrane-integrated hybrid treatment process that ensures reuse of water with conversion and recovery of ammoniacal nitrogen as value-added struvite fertilizer from coke wastewater. While toxic cyanide was largely removed in a pre-chemical treatment unit using Fenton's reagents under optimized conditions, more than 95% of NH4(+)-N could be recovered as a valuable by-product called struvite through addition of appropriate doses of magnesium and phosphate salts. Water could be turned reusable through a polishing treatment by nanofiltration membranes in a largely fouling free membrane module following a biodegradation step. Mathematical modeling of such an integrated process was done with Haldane-Andrew approach for the associated microbial degradation of phenol by Pseudomonas putida. Residual NH4(+) was degraded by nitrification and denitrification following the modified Monod kinetics. The model could successfully predict the plant performance as reflected in reasonably low relative error (0.03-0.18) and high Willmott d-index (>0.98). PMID:23735488

  8. A GPU-accelerated semi-implicit ADI method for incompressible and compressible Navier-Stokes equations

    NASA Astrophysics Data System (ADS)

    Ha, Sanghyun; You, Donghyun

    2015-11-01

    Utility of the computational power of Graphics Processing Units (GPUs) is elaborated for solutions of both incompressible and compressible Navier-Stokes equations. A semi-implicit ADI finite-volume method for integration of the incompressible and compressible Navier-Stokes equations, which are discretized on a structured arbitrary grid, is parallelized for GPU computations using CUDA (Compute Unified Device Architecture). In the semi-implicit ADI finite-volume method, the nonlinear convection terms and the linear diffusion terms are integrated in time using a combination of an explicit scheme and an ADI scheme. Inversion of multiple tri-diagonal matrices is found to be the major challenge in GPU computations of the present method. Some of the algorithms for solving tri-diagonal matrices on GPUs are evaluated and optimized for GPU-acceleration of the present semi-implicit ADI computations of incompressible and compressible Navier-Stokes equations. Supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning Grant NRF-2014R1A2A1A11049599.

  9. Implicit Learning as an Ability

    ERIC Educational Resources Information Center

    Kaufman, Scott Barry; DeYoung, Caroline G.; Gray, Jeremy R.; Jimenez, Luis; Brown, Jamie; Mackintosh, Nicholas

    2010-01-01

    The ability to automatically and implicitly detect complex and noisy regularities in the environment is a fundamental aspect of human cognition. Despite considerable interest in implicit processes, few researchers have conceptualized implicit learning as an ability with meaningful individual differences. Instead, various researchers (e.g., Reber,…

  10. The lambda-scheme. [for numerical integration of Euler equation of compressible gas flow

    NASA Technical Reports Server (NTRS)

    Moretti, G.

    1979-01-01

    A method for integrating the Euler equations of gas dynamics for compressible flows in any hyperbolic case is presented. This method is applied to the Mach number distribution over a stretch of an infinite duct having a variable cross section, and to the distribution in a channel opening into a vacuum with the Mach number equalling 1.04. An example of the ability of this method to handle two-dimensional unsteady flows is shown using the steady shock-and-isobars pattern reached asymptotically about an ablated blunt body with a free stream Mach number equalling 12. A final example is presented where the technique is applied to a three-dimensional steady supersonic flow, with a Mach number of 2 and an angle of attack of 5 deg.

  11. An integrated collision prediction and avoidance scheme for mobile robots in non-stationary environments

    NASA Technical Reports Server (NTRS)

    Kyriakopoulos, K. J.; Saridis, G. N.

    1993-01-01

    A formulation that makes possible the integration of collision prediction and avoidance stages for mobile robots moving in general terrains containing moving obstacles is presented. A dynamic model of the mobile robot and the dynamic constraints are derived. Collision avoidance is guaranteed if the distance between the robot and a moving obstacle is nonzero. A nominal trajectory is assumed to be known from off-line planning. The main idea is to change the velocity along the nominal trajectory so that collisions are avoided. A feedback control is developed and local asymptotic stability is proved if the velocity of the moving obstacle is bounded. Furthermore, a solution to the problem of inverse dynamics for the mobile robot is given. Simulation results verify the value of the proposed strategy.

  12. A MEMS-Based Flow Rate and Flow Direction Sensing Platform with Integrated Temperature Compensation Scheme.

    PubMed

    Ma, Rong-Hua; Wang, Dung-An; Hsueh, Tzu-Han; Lee, Chia-Yen

    2009-01-01

    This study develops a MEMS-based low-cost sensing platform for sensing gas flow rate and flow direction comprising four silicon nitride cantilever beams arranged in a cross-form configuration, a circular hot-wire flow meter suspended on a silicon nitride membrane, and an integrated resistive temperature detector (RTD). In the proposed device, the flow rate is inversely derived from the change in the resistance signal of the flow meter when exposed to the sensed air stream. To compensate for the effects of the ambient temperature on the accuracy of the flow rate measurements, the output signal from the flow meter is compensated using the resistance signal generated by the RTD. As air travels over the surface of the cross-form cantilever structure, the upstream cantilevers are deflected in the downward direction, while the downstream cantilevers are deflected in the upward direction. The deflection of the cantilever beams causes a corresponding change in the resistive signals of the piezoresistors patterned on their upper surfaces. The amount by which each beam deflects depends on both the flow rate and the orientation of the beam relative to the direction of the gas flow. Thus, following an appropriate compensation by the temperature-corrected flow rate, the gas flow direction can be determined through a suitable manipulation of the output signals of the four piezoresistors. The experimental results have confirmed that the resulting variation in the output signals of the integrated sensors can be used to determine not only the ambient temperature and the velocity of the air flow, but also its direction relative to the sensor with an accuracy of ± 7.5° error. PMID:22346708

  13. A MEMS-Based Flow Rate and Flow Direction Sensing Platform with Integrated Temperature Compensation Scheme

    PubMed Central

    Ma, Rong-Hua; Wang, Dung-An; Hsueh, Tzu-Han; Lee, Chia-Yen

    2009-01-01

    This study develops a MEMS-based low-cost sensing platform for sensing gas flow rate and flow direction comprising four silicon nitride cantilever beams arranged in a cross-form configuration, a circular hot-wire flow meter suspended on a silicon nitride membrane, and an integrated resistive temperature detector (RTD). In the proposed device, the flow rate is inversely derived from the change in the resistance signal of the flow meter when exposed to the sensed air stream. To compensate for the effects of the ambient temperature on the accuracy of the flow rate measurements, the output signal from the flow meter is compensated using the resistance signal generated by the RTD. As air travels over the surface of the cross-form cantilever structure, the upstream cantilevers are deflected in the downward direction, while the downstream cantilevers are deflected in the upward direction. The deflection of the cantilever beams causes a corresponding change in the resistive signals of the piezoresistors patterned on their upper surfaces. The amount by which each beam deflects depends on both the flow rate and the orientation of the beam relative to the direction of the gas flow. Thus, following an appropriate compensation by the temperature-corrected flow rate, the gas flow direction can be determined through a suitable manipulation of the output signals of the four piezoresistors. The experimental results have confirmed that the resulting variation in the output signals of the integrated sensors can be used to determine not only the ambient temperature and the velocity of the air flow, but also its direction relative to the sensor with an accuracy of ± 7.5° error. PMID:22346708

  14. Extrapolated implicit-explicit time stepping.

    SciTech Connect

    Constantinescu, E. M.; Sandu, A.; Mathematics and Computer Science; Virginia Polytechnic Inst. and State Univ.

    2010-01-01

    This paper constructs extrapolated implicit-explicit time stepping methods that allow one to efficiently solve problems with both stiff and nonstiff components. The proposed methods are based on Euler steps and can provide very high order discretizations of ODEs, index-1 DAEs, and PDEs in the method-of-lines framework. Implicit-explicit schemes based on extrapolation are simple to construct, easy to implement, and straightforward to parallelize. This work establishes the existence of perturbed asymptotic expansions of global errors, explains the convergence orders of these methods, and studies their linear stability properties. Numerical results with stiff ODE, DAE, and PDE test problems confirm the theoretical findings and illustrate the potential of these methods to solve multiphysics multiscale problems.

  15. Integrating Physical and Topographic Information Into a Fuzzy Scheme to Map Flooded Area by SAR

    PubMed Central

    Pierdicca, Nazzareno; Chini, Marco; Pulvirenti, Luca; Macina, Flavia

    2008-01-01

    A flood mapping procedure based on a fuzzy sets theory has been developed. The method is based on the integration of Synthetic Aperture Radar (SAR) measurements with additional data on the inundated area, such as a land cover map and a digital elevation model (DEM). The information on land cover has allowed us to account for both specular reflection, typical of open water, and double bounce backscattering, typical of forested and urban areas. DEM has been exploited to include simple hydraulic considerations on the dependence of inundation probability on surface characteristics. Contextual information has been taken into account too. The proposed algorithm has been tested on a flood occurred in Italy on November 1994. A pair of ERS-1 images, collected before and after (three days later) the flood, has been used. The results have been compared with the data provided by a ground survey carried out when the flood reached its maximum extension. Despite the temporal mismatch between the survey and the post-inundation SAR image, the comparison has yielded encouraging results, with the 87% of the pixels correctly classified as inundated.

  16. An integrated voice and data multiple-access scheme for a land-mobile satellite system

    NASA Technical Reports Server (NTRS)

    Li, V. O. K.; Yan, T.-Y.

    1984-01-01

    An analytical study is performed of the satellite requirements for a land mobile satellite system (LMSS). The spacecraft (MSAT-X) would be in GEO and would be compatible with multiple access by mobile radios and antennas and fixed stations. The FCC has received a petition from NASA to reserve the 821-825 and 866-870 MHz frequencies for the LMSS, while communications with fixed earth stations would be in the Ku band. MSAT-X transponders would alter the frequencies of signal and do no processing in the original configuration considered. Channel use would be governed by an integrated demand-assigned, multiple access protocol, which would divide channels into reservation and information channels, governed by a network management center. Further analyses will cover tradeoffs between data and voice users, probability of blocking, and the performance impacts of on-board switching and variable bandwidth assignment. Initial calculations indicate that a large traffic volume can be handled with acceptable delays and voice blocking probabilities.

  17. Optimizing the updated Goddard shortwave radiation Weather Research and Forecasting (WRF) scheme for Intel Many Integrated Core (MIC) architecture

    NASA Astrophysics Data System (ADS)

    Mielikainen, Jarno; Huang, Bormin; Huang, Allen H.-L.

    2015-05-01

    Intel Many Integrated Core (MIC) ushers in a new era of supercomputing speed, performance, and compatibility. It allows the developers to run code at trillions of calculations per second using the familiar programming model. In this paper, we present our results of optimizing the updated Goddard shortwave radiation Weather Research and Forecasting (WRF) scheme on Intel Many Integrated Core Architecture (MIC) hardware. The Intel Xeon Phi coprocessor is the first product based on Intel MIC architecture, and it consists of up to 61 cores connected by a high performance on-die bidirectional interconnect. The co-processor supports all important Intel development tools. Thus, the development environment is familiar one to a vast number of CPU developers. Although, getting a maximum performance out of Xeon Phi will require using some novel optimization techniques. Those optimization techniques are discusses in this paper. The results show that the optimizations improved performance of the original code on Xeon Phi 7120P by a factor of 1.3x.

  18. The time course of explicit and implicit categorization.

    PubMed

    Smith, J David; Zakrzewski, Alexandria C; Herberger, Eric R; Boomer, Joseph; Roeder, Jessica L; Ashby, F Gregory; Church, Barbara A

    2015-10-01

    Contemporary theory in cognitive neuroscience distinguishes, among the processes and utilities that serve categorization, explicit and implicit systems of category learning that learn, respectively, category rules by active hypothesis testing or adaptive behaviors by association and reinforcement. Little is known about the time course of categorization within these systems. Accordingly, the present experiments contrasted tasks that fostered explicit categorization (because they had a one-dimensional, rule-based solution) or implicit categorization (because they had a two-dimensional, information-integration solution). In Experiment 1, participants learned categories under unspeeded or speeded conditions. In Experiment 2, they applied previously trained category knowledge under unspeeded or speeded conditions. Speeded conditions selectively impaired implicit category learning and implicit mature categorization. These results illuminate the processing dynamics of explicit/implicit categorization. PMID:26025556

  19. Detector density and small field dosimetry: Integral versus point dose measurement schemes

    SciTech Connect

    Underwood, T. S. A. Hill, M. A.; Fenwick, J. D.

    2013-08-01

    Purpose: TheAlfonso et al. [Med. Phys.35, 5179–5186 (2008)] formalism for small field dosimetry proposes a set of correction factors (k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}}) which account for differences between the detector response in nonstandard (clinical) and machine-specific-reference fields. In this study, the Monte Carlo method was used to investigate the viability of such small field correction factors for four different detectors irradiated under a variety of conditions. Because k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} values for single detector position measurements are influenced by several factors, a new theoretical formalism for integrated-detector-position [dose area product (DAP)] measurements is also presented and was tested using Monte Carlo simulations. Methods: A BEAMnrc linac model was built and validated for a Varian Clinac iX accelerator. Using the egs++ geometry package, detailed virtual models were built for four different detectors: a PTW 60012 unshielded diode, a PTW 60003 Diamond detector, a PTW 31006 PinPoint (ionization chamber), and a PTW 31018 MicroLion (liquid-filled ionization chamber). The egs-chamber code was used to investigate the variation ofk{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} with detector type, detector construction, field size, off-axis position, and the azimuthal angle between the detector and beam axis. Simulations were also used to consider the DAP obtained by each detector: virtual detectors and water voxels were scanned through high resolution grids of positions extending far beyond the boundaries of the fields under consideration. Results: For each detector, the correction factor (k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s

  20. Asymptotic analysis of discrete schemes for non-equilibrium radiation diffusion

    NASA Astrophysics Data System (ADS)

    Cui, Xia; Yuan, Guang-wei; Shen, Zhi-jun

    2016-05-01

    Motivated by providing well-behaved fully discrete schemes in practice, this paper extends the asymptotic analysis on time integration methods for non-equilibrium radiation diffusion in [2] to space discretizations. Therein studies were carried out on a two-temperature model with Larsen's flux-limited diffusion operator, both the implicitly balanced (IB) and linearly implicit (LI) methods were shown asymptotic-preserving. In this paper, we focus on asymptotic analysis for space discrete schemes in dimensions one and two. First, in construction of the schemes, in contrast to traditional first-order approximations, asymmetric second-order accurate spatial approximations are devised for flux-limiters on boundary, and discrete schemes with second-order accuracy on global spatial domain are acquired consequently. Then by employing formal asymptotic analysis, the first-order asymptotic-preserving property for these schemes and furthermore for the fully discrete schemes is shown. Finally, with the help of manufactured solutions, numerical tests are performed, which demonstrate quantitatively the fully discrete schemes with IB time evolution indeed have the accuracy and asymptotic convergence as theory predicts, hence are well qualified for both non-equilibrium and equilibrium radiation diffusion.

  1. Hybrid wide-band, low-phase-noise scheme for Raman lasers in atom interferometry by integrating an acousto-optic modulator and a feedback loop.

    PubMed

    Wang, Kai; Yao, Zhanwei; Li, Runbing; Lu, Sibin; Chen, Xi; Wang, Jin; Zhan, Mingsheng

    2016-02-10

    We report a hybrid scheme for phase-coherent Raman lasers with low phase noise in a wide frequency range. In this scheme, a pair of Raman lasers with a frequency difference of 3.04 GHz is generated by the ±1-order diffracted lights of an acousto-optic modulator (1.52 GHz), where a feedback loop is simultaneously applied for suppressing the phase noise. The beat width of the Raman lasers is narrower than 3 Hz. In the low-frequency range, the phase noise of the Raman lasers is suppressed by 35 dB with the feedback. The phase noise is less than -109  dBc/Hz in the high-frequency range. The sensitivity of an atom gyroscope employing the hybrid Raman lasers can be implicitly improved 10 times. Due to the better high-frequency response, the sensitivity is not limited by the durations of Raman pulses. This work is important for improving the performance of atom-interferometer-based measurements. PMID:26906364

  2. Applications of high-resolution spatial discretization scheme and Jacobian-free Newton–Krylov method in two-phase flow problems

    SciTech Connect

    Zou, Ling; Zhao, Haihua; Zhang, Hongbin

    2015-09-01

    The majority of the existing reactor system analysis codes were developed using low-order numerical schemes in both space and time. In many nuclear thermal–hydraulics applications, it is desirable to use higher-order numerical schemes to reduce numerical errors. High-resolution spatial discretization schemes provide high order spatial accuracy in smooth regions and capture sharp spatial discontinuity without nonphysical spatial oscillations. In this work, we adapted an existing high-resolution spatial discretization scheme on staggered grids in two-phase flow applications. Fully implicit time integration schemes were also implemented to reduce numerical errors from operator-splitting types of time integration schemes. The resulting nonlinear system has been successfully solved using the Jacobian-free Newton–Krylov (JFNK) method. The high-resolution spatial discretization and high-order fully implicit time integration numerical schemes were tested and numerically verified for several two-phase test problems, including a two-phase advection problem, a two-phase advection with phase appearance/disappearance problem, and the water faucet problem. Numerical results clearly demonstrated the advantages of using such high-resolution spatial and high-order temporal numerical schemes to significantly reduce numerical diffusion and therefore improve accuracy. Our study also demonstrated that the JFNK method is stable and robust in solving two-phase flow problems, even when phase appearance/disappearance exists.

  3. 3D modelling of the active normal fault network in the Apulian Ridge (Eastern Mediterranean Sea): Integration of seismic and bathymetric data with implicit surface methods

    NASA Astrophysics Data System (ADS)

    Bistacchi, Andrea; Pellegrini, Caludio; Savini, Alessandra; Marchese, Fabio

    2016-04-01

    The Apulian ridge (North-eastern Ionian Sea, Mediterranean), interposed between the facing Apennines and Hellenides subduction zones (to the west and east respectively), is characterized by thick cretaceous carbonatic sequences and discontinuous tertiary deposits crosscut by a penetrative network of NNW-SSE normal faults. These are exposed onshore in Puglia, and are well represented offshore in a dataset composed of 2D seismics and wells collected by oil companies from the '60s to the '80s, more recent seismics collected during research projects in the '90s, recent very high resolution seismics (VHRS - Sparker and Chirp-sonar data), multibeam echosounder bathymetry, and sedimentological and geo-chronological analyses of sediment samples collected on the seabed. Faults are evident in 2D seismics at all scales, and their along-strike geometry and continuity can be characterized with multibeam bathymetric data, which show continuous fault scarps on the seabed (only partly reworked by currents and covered by landslides). Fault scarps also reveal the finite displacement accumulated in the Holocene-Pleistocene. We reconstructed a 3D model of the fault network and suitable geological boundaries (mainly unconformities due to the discontinuous distribution of quaternary and tertiary sediments) with implicit surface methods implemented in SKUA/GOCAD. This approach can be considered very effective and allowed reconstructing in details complex structures, like the frequent relay zones that are particularly well imaged by seafloor geomorphology. Mutual cross-cutting relationships have been recognized between fault scarps and submarine mass-wasting deposits (Holocene-Pleistocene), indicating that, at least in places, these features are coeval, hence the fault network should be considered active. At the regional scale, the 3D model allowed measuring the horizontal WSW-ENE stretching, which can be associated to the bending moment applied to the Apulian Plate by the combined effect

  4. Implicit measure for yoga research: Yoga implicit association test

    PubMed Central

    Ilavarasu, Judu V; Rajesh, Sasidharan K; Hankey, Alex

    2014-01-01

    Context: The implicit association test (IAT), a new tool for yoga research is presented. Implicit measures could be used in those situations where (1) The construct is difficult to self-report, (2) there is a threat of social desirability. Clinically, we can assess cognitive dissonance by evaluating incongruence between implicit and explicit measures. Explicit preferences are self-reported. Implicit preferences are what we inherently believe, often without our conscious awareness. Aims: The primary objective of this study is to provide a bird's eye view of the field, implicit cognition, with emphasis on the IAT and the secondary objective is to illustrate through an example of our study to develop an implicit tool to assess implicit preference toward yoga. Settings and Design: A total of 5 independent samples of total 69 students undergoing short and long-term yoga courses in a Yoga University were assessed for their implicit and explicit preferences towards yoga. Materials and Methods: The yoga-IAT (Y-IAT), explicit self-rating scale was administered through computers using the Inquisit program by Millisecond Software. Experimental and scoring materials are provided. Results: A moderate preference toward yoga was detected, with a lower implicit-explicit congruence, reflecting possible confound of social desirability in the self-report of preference toward yoga. Conclusions: Implicit measures may be used in the yoga field to assess constructs, which are difficult to self-report or may have social desirability threat. Y-IAT may be used to evaluate implicit preference toward yoga. PMID:25035621

  5. Explicit and Implicit Emotion Regulation: A Dual-Process Framework

    PubMed Central

    Gyurak, Anett; Gross, James J.; Etkin, Amit

    2012-01-01

    It is widely acknowledged that emotions can be regulated in an astonishing variety of ways. Most research to date has focused on explicit (effortful) forms of emotion regulation. However, there is growing research interest in implicit (automatic) forms of emotion regulation. To organize emerging findings, we present a dual-process framework that integrates explicit and implicit forms of emotion regulation, and argue that both forms of regulation are necessary for well-being. In the first section of this review, we provide a broad overview of the construct of emotion regulation, with an emphasis on explicit and implicit processes. In the second section, we focus on explicit emotion regulation, considering both neural mechanisms that are associated with these processes and their experiential and physiological consequences. In the third section, we turn to several forms of implicit emotion regulation, and integrate the burgeoning literature in this area. We conclude by outlining open questions and areas for future research. PMID:21432682

  6. Implicit Cognition and Addiction: A Tool for Explaining Paradoxical Behavior

    PubMed Central

    Stacy, Alan W.; Wiers, Reinout W.

    2012-01-01

    Research on implicit cognition and addiction has expanded greatly during the past decade. This research area provides new ways to understand why people engage in behaviors that they know are harmful or counterproductive in the long run. Implicit cognition takes a different view from traditional cognitive approaches to addiction by assuming that behavior is often not a result of a reflective decision that takes into account the pros and cons known by the individual. Instead of a cognitive algebra integrating many cognitions relevant to choice, implicit cognition assumes that the influential cognitions are the ones that are spontaneously activated during critical decision points. This selective review highlights many of the consistent findings supporting predictive effects of implicit cognition on substance use and abuse in adolescents and adults; reveals a recent integration with dual-process models; outlines the rapid evolution of different measurement tools; and introduces new routes for intervention. PMID:20192786

  7. Solving ODE Initial Value Problems With Implicit Taylor Series Methods

    NASA Technical Reports Server (NTRS)

    Scott, James R.

    2000-01-01

    In this paper we introduce a new class of numerical methods for integrating ODE initial value problems. Specifically, we propose an extension of the Taylor series method which significantly improves its accuracy and stability while also increasing its range of applicability. To advance the solution from t (sub n) to t (sub n+1), we expand a series about the intermediate point t (sub n+mu):=t (sub n) + mu h, where h is the stepsize and mu is an arbitrary parameter called an expansion coefficient. We show that, in general, a Taylor series of degree k has exactly k expansion coefficients which raise its order of accuracy. The accuracy is raised by one order if k is odd, and by two orders if k is even. In addition, if k is three or greater, local extrapolation can be used to raise the accuracy two additional orders. We also examine stability for the problem y'= lambda y, Re (lambda) less than 0, and identify several A-stable schemes. Numerical results are presented for both fixed and variable stepsizes. It is shown that implicit Taylor series methods provide an effective integration tool for most problems, including stiff systems and ODE's with a singular point.

  8. Dusty gas with smoothed particle hydrodynamics - II. Implicit timestepping and astrophysical drag regimes

    NASA Astrophysics Data System (ADS)

    Laibe, Guillaume; Price, Daniel J.

    2012-03-01

    In a companion paper, we have presented an algorithm for simulating two-fluid gas and dust mixtures in smoothed particle hydrodynamics (SPH). In this paper, we develop an implicit timestepping method that preserves the exact conservation of the both linear and angular momenta in the underlying SPH algorithm, but unlike previous schemes, allows the iterations to converge to arbitrary accuracy and is suited to the treatment of non-linear drag regimes. The algorithm presented in Paper I is also extended to deal with realistic astrophysical drag regimes, including both linear and non-linear Epstein and Stokes drag. The scheme is benchmarked against the test suite presented in Paper I, including (i) the analytic solutions of the DUSTYBOX problem and (ii) solutions of the DUSTYWAVE, DUSTYSHOCK, DUSTYSEDOV and DUSTYDISC obtained with explicit timestepping. We find that the implicit method is 1-10 times faster than the explicit temporal integration when the ratio r between the timestep and the drag stopping time is 1 ≲r≲ 1000.

  9. Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms

    NASA Technical Reports Server (NTRS)

    Yee, H. C.; Shinn, J. L.

    1986-01-01

    Some numerical aspects of finite-difference algorithms for nonlinear multidimensional hyperbolic conservation laws with stiff nonhomogenous (source) terms are discussed. If the stiffness is entirely dominated by the source term, a semi-implicit shock-capturing method is proposed provided that the Jacobian of the soruce terms possesses certain properties. The proposed semi-implicit method can be viewed as a variant of the Bussing and Murman point-implicit scheme with a more appropriate numerical dissipation for the computation of strong shock waves. However, if the stiffness is not solely dominated by the source terms, a fully implicit method would be a better choice. The situation is complicated by problems that are higher than one dimension, and the presence of stiff source terms further complicates the solution procedures for alternating direction implicit (ADI) methods. Several alternatives are discussed. The primary motivation for constructing these schemes was to address thermally and chemically nonequilibrium flows in the hypersonic regime. Due to the unique structure of the eigenvalues and eigenvectors for fluid flows of this type, the computation can be simplified, thus providing a more efficient solution procedure than one might have anticipated.

  10. Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms

    NASA Technical Reports Server (NTRS)

    Yee, H. C.; Shinn, Judy L.

    1987-01-01

    Some numerical aspects of finite-difference algorithms for nonlinear multidimensional hyperbolic conservation laws with stiff nonhomogeneous (source) terms are discussed. If the stiffness is entirely dominated by the source term, a semi-implicit shock-capturing method is proposed provided that the Jacobian of the source terms possesses certain properties. The proposed semi-implicit method can be viewed as a variant of the Bussing and Murman point-implicit scheme with a more appropriate numerical dissipation for the computation of strong shock waves. However, if the stiffness is not solely dominated by the source terms, a fully implicit method would be a better choice. The situation is complicated by problems that are higher than one dimension, and the presence of stiff source terms further complicates the solution procedures for alternating direction implicit (ADI) methods. Several alternatives are discussed. The primary motivation for constructing these schemes was to address thermally and chemically nonequilibrium flows in the hypersonic regime. Due to the unique structure of the eigenvalues and eigenvectors for fluid flows of this type, the computation can be simplified, thus providing a more efficient solution procedure than one might have anticipated.