Using 4th order Runge-Kutta method for solving a twisted Skyrme string equation
NASA Astrophysics Data System (ADS)
Hadi, Miftachul; Anderson, Malcolm; Husein, Andri
2016-03-01
We study numerical solution, especially using 4th order Runge-Kutta method, for solving a twisted Skyrme string equation. We find numerically that the value of minimum energy per unit length of vortex solution for a twisted Skyrmion string is 20.37 × 1060 eV/m.
A third order Runge-Kutta algorithm on a manifold
NASA Technical Reports Server (NTRS)
Crouch, P. E.; Grossman, R. G.; Yan, Y.
1992-01-01
A third order Runge-Kutta type algorithm is described with the property that it preserves certain geometric structures. In particular, if the algorithm is initialized on a Lie group, then the resulting iterates remain on the Lie group.
Generation and application of the equations of condition for high order Runge-Kutta methods
NASA Technical Reports Server (NTRS)
Haley, D. C.
1972-01-01
This thesis develops the equations of condition necessary for determining the coefficients for Runge-Kutta methods used in the solution of ordinary differential equations. The equations of condition are developed for Runge-Kutta methods of order four through order nine. Once developed, these equations are used in a comparison of the local truncation errors for several sets of Runge-Kutta coefficients for methods of order three up through methods of order eight.
Optimized fourth-order Runge-Kutta method for solving oscillatory problems
NASA Astrophysics Data System (ADS)
Hussain, Kasim; Ismail, Fudziah; Senu, Norazak; Rabiei, Faranak
2016-06-01
In this article, we develop a Runge-Kutta method with invalidation of phase lag, phase lag's derivatives and amplification error to solve second-order initial value problem (IVP) with oscillating solutions. The new method depends on the explicit Runge-Kutta method of algebraic order four. Numerical tests from its implementation to well-known oscillatory problems illustrate the robustness and competence of the new method as compared to the well-known Runge-Kutta methods in the scientific literature.
Diagonally Implicit Symplectic Runge-Kutta Methods with High Algebraic and Dispersion Order
Cong, Y. H.; Jiang, C. X.
2014-01-01
The numerical integration of Hamiltonian systems with oscillating solutions is considered in this paper. A diagonally implicit symplectic nine-stages Runge-Kutta method with algebraic order 6 and dispersion order 8 is presented. Numerical experiments with some Hamiltonian oscillatory problems are presented to show the proposed method is as competitive as the existing same type Runge-Kutta methods. PMID:24977178
Diagonally implicit symplectic Runge-Kutta methods with high algebraic and dispersion order.
Cong, Y H; Jiang, C X
2014-01-01
The numerical integration of Hamiltonian systems with oscillating solutions is considered in this paper. A diagonally implicit symplectic nine-stages Runge-Kutta method with algebraic order 6 and dispersion order 8 is presented. Numerical experiments with some Hamiltonian oscillatory problems are presented to show the proposed method is as competitive as the existing same type Runge-Kutta methods. PMID:24977178
Fourth-order 2N-storage Runge-Kutta schemes
NASA Technical Reports Server (NTRS)
Carpenter, Mark H.; Kennedy, Christopher A.
1994-01-01
A family of five-stage fourth-order Runge-Kutta schemes is derived; these schemes required only two storage locations. A particular scheme is identified that has desirable efficiency characteristics for hyperbolic and parabolic initial (boundary) value problems. This scheme is competitive with the classical fourth-order method (high-storage) and is considerably more efficient and accurate than existing third-order low-storage schemes.
Third-order 2N-storage Runge-Kutta schemes with error control
NASA Technical Reports Server (NTRS)
Carpenter, Mark H.; Kennedy, Christopher A.
1994-01-01
A family of four-stage third-order explicit Runge-Kutta schemes is derived that requires only two storage locations and has desirable stability characteristics. Error control is achieved by embedding a second-order scheme within the four-stage procedure. Certain schemes are identified that are as efficient and accurate as conventional embedded schemes of comparable order and require fewer storage locations.
Equations of condition for high order Runge-Kutta-Nystrom formulae
NASA Technical Reports Server (NTRS)
Bettis, D. G.
1974-01-01
Derivation of the equations of condition of order eight for a general system of second-order differential equations approximated by the basic Runge-Kutta-Nystrom algorithm. For this general case, the number of equations of condition is considerably larger than for the special case where the first derivative is not present. Specifically, it is shown that, for orders two through eight, the number of equations for each order is 1, 1, 1, 2, 3, 5, and 9 for the special case and is 1, 1, 2, 5, 13, 34, and 95 for the general case.
NASA Technical Reports Server (NTRS)
Fehlberg, E.
1973-01-01
New Runge-Kutta-Nystrom formulas of the eighth, seventh, sixth, and fifth order are derived for the special second-order (vector) differential equation x = f (t,x). In contrast to Runge-Kutta-Nystrom formulas of an earlier NASA report, these formulas provide a stepsize control procedure based on the leading term of the local truncation error in x. This new procedure is more accurate than the earlier Runge-Kutta-Nystrom procedure (with stepsize control based on the leading term of the local truncation error in x) when integrating close to singularities. Two central orbits are presented as examples. For these orbits, the accuracy and speed of the formulas of this report are compared with those of Runge-Kutta-Nystrom and Runge-Kutta formulas of earlier NASA reports.
Qualitatively stability of nonstandard 2-stage explicit Runge-Kutta methods of order two
NASA Astrophysics Data System (ADS)
Khalsaraei, M. M.; Khodadosti, F.
2016-02-01
When one solves differential equations, modeling physical phenomena, it is of great importance to take physical constraints into account. More precisely, numerical schemes have to be designed such that discrete solutions satisfy the same constraints as exact solutions. Nonstandard finite differences (NSFDs) schemes can improve the accuracy and reduce computational costs of traditional finite difference schemes. In addition NSFDs produce numerical solutions which also exhibit essential properties of solution. In this paper, a class of nonstandard 2-stage Runge-Kutta methods of order two (we call it nonstandard RK2) is considered. The preservation of some qualitative properties by this class of methods are discussed. In order to illustrate our results, we provide some numerical examples.
NASA Technical Reports Server (NTRS)
Fehlberg, E.
1974-01-01
Runge-Kutta-Nystrom formulas of the seventh, sixth, and fifth order were derived for the general second order (vector) differential equation written as the second derivative of x = f(t, x, the first derivative of x). The formulas include a stepsize control procedure, based on a complete coverage of the leading term of the local truncation error in x, and they require no more evaluations per step than the earlier Runge-Kutta formulas for the first derivative of x = f(t, x). The developed formulas are expected to be time saving in comparison to the Runge-Kutta formulas for first-order differential equations, since it is not necessary to convert the second-order differential equations into twice as many first-order differential equations. The examples shown saved from 25 percent to 60 percent more computer time than the earlier formulas for first-order differential equations, and are comparable in accuracy.
Evolutionary generation of 7th order Runge - Kutta - Nyström type methods for solving y(4) = f(x,y)
NASA Astrophysics Data System (ADS)
Papakostas, S. N.; Tsitmidelis, S.; Tsitouras, Ch.
2015-12-01
We present a 7th algebraic order Runge - Kutta - Nyström method for the solution of a special fourth order initial value problem. To achieve this, a set of non - linear equations is solved using differential evolution technique. Various numerical tests justify our efforts.
Evolutionary generation of high order Runge - Kutta - Nyström type pairs for solving y(4) = f (x,y)
NASA Astrophysics Data System (ADS)
Famelis, I. Th.; Tsitmidelis, S.; Tsitouras, Ch.
2016-06-01
We present a new Runge - Kutta - Nyström type pair of orders 8(6) for the solution of a special fourth order initial value problem. To achieve this, a set of non - linear equations is solved using differential evolution technique.
Starting methods for two-step Runge-Kutta methods of stage-order 3 and order 6
NASA Astrophysics Data System (ADS)
Verner, J. H.
2006-01-01
Jackiewicz and Tracogna [SIAM J. Numer. Anal. 32 (1995) 1390-1427] proposed a general formulation of two step Runge-Kutta (TSRK) methods. Using formulas for two-step pairs of TSRK methods constructed in [Japan JIAM 19 (2002) 227-248], Jackiewicz and Verner obtain results for order 8 pairs that fail to show this designated order. Hairer and Wanner [SIAM J. Numer. Anal. 34 (1997) 2087-2089] identify the problem by using B-series to formulate a complete set of order conditions for TSRK methods, and emphasize that special starting methods are necessary for the first step of implementation. They observe that for methods with stage order at least p-1, and design order p, starting methods of order at least p are sufficient. In this paper, the more general challenge to provide correct starting values for methods of low stage-order is met by showing how perturbed starting values should be selected for methods of order 6 and stage-order 3. The approach is sufficiently general that it may (and later will) be provided for such methods of higher orders. Evidence of the accompanying improvement in the implementation of TSRK methods illustrates that carefully designed starting methods are essential for efficient production codes based on methods of low stage-order.
Scaled Runge-Kutta algorithms for handling dense output
NASA Technical Reports Server (NTRS)
Horn, M. K.
1981-01-01
Low order Runge-Kutta algorithms are developed which determine the solution of a system of ordinary differential equations at any point within a given integration step, as well as at the end of each step. The scaled Runge-Kutta methods are designed to be used with existing Runge-Kutta formulas, using the derivative evaluations of these defining algorithms as the core of the system. For a slight increase in computing time, the solution may be generated within the integration step, improving the efficiency of the Runge-Kutta algorithms, since the step length need no longer be severely reduced to coincide with the desired output point. Scaled Runge-Kutta algorithms are presented for orders 3 through 5, along with accuracy comparisons between the defining algorithms and their scaled versions for a test problem.
A Runge-Kutta Nystrom algorithm.
NASA Technical Reports Server (NTRS)
Bettis, D. G.
1973-01-01
A Runge-Kutta algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative. The algorithm includes the Fehlberg step control procedure.
NASA Technical Reports Server (NTRS)
Fehlberg, E.
1972-01-01
The formulas include a stepsize control procedure, based on a complete coverage of the leading term of the truncation error in x. The formulas require fewer evaluations per stop than other Runge-Kutta-Nystrom formulas if the latter are operated by using the standard procedure for stepsize control. An example is presented. With results being of the same accuracy, Runge-Kutta-Nystrom formulas discussed save 50 percent or more computer time compared with other Runge-Kutta-Nystrom formulas.
Explicit Runge-Kutta method with trigonometrically-fitted for solving first order ODEs
NASA Astrophysics Data System (ADS)
Fawzi, Firas Adel; Senu, N.; Ismail, F.; Majid, Z. A.
2016-06-01
In this note, an explicit trigonometrically-fitted (RK) method is developed to determine the approximate solution of the first-order IVPs with oscillatory solution. The proposed method solves first order ODEs by first converting the second order ODEs to an equivalent first order; which is based on the RK method of order four. The numerical experiment performed shows the efficacy of our newly developed method.
NASA Astrophysics Data System (ADS)
Romeo, A.; Finocchio, G.; Carpentieri, M.; Torres, L.; Consolo, G.; Azzerboni, B.
2008-02-01
The Landau-Lifshitz-Gilbert (LLG) equation is the fundamental equation to describe magnetization dynamics in microscale and nanoscale magnetic systems. In this paper we present a brief overview of a time-domain numerical method related to the fifth order Runge-Kutta formula, which has been applied to the solution of the LLG equation successfully. We discuss advantages of the method, describing the results of a numerical experiment based on the standard problem #4. The results are in good agreement with the ones present in literature. By including thermal effects in our framework, our simulations show magnetization dynamics slightly dependent on the spatial discretization.
Composite group of explicit Runge-Kutta methods
NASA Astrophysics Data System (ADS)
Hamid, Fatin Nadiah Abd; Rabiei, Faranak; Ismail, Fudziah
2016-06-01
In this paper,the composite groups of Runge-Kutta (RK) method are proposed. The composite group of RK method of third and second order, RK3(2) and fourth and third order RK4(3) base on classical Runge-Kutta method are derived. The proposed methods are two-step in nature and have less number of function evaluations compared to the existing Runge-Kutta method. The order conditions up to order four are obtained using rooted trees and composite rule introduced by J. C Butcher. The stability regions of RK3(2) and RK4(3) methods are presented and initial value problems of first order ordinary differential equations are carried out. Numerical results are compared with existing Runge-Kutta method.
Galerkin/Runge-Kutta discretizations for semilinear parabolic equations
NASA Technical Reports Server (NTRS)
Keeling, Stephen L.
1987-01-01
A new class of fully discrete Galerkin/Runge-Kutta methods is constructed and analyzed for semilinear parabolic initial boundary value problems. Unlike any classical counterpart, this class offers arbitrarily high, optimal order convergence. In support of this claim, error estimates are proved, and computational results are presented. Furthermore, it is noted that special Runge-Kutta methods allow computations to be performed in parallel so that the final execution time can be reduced to that of a low order method.
Functional continuous Runge-Kutta methods for special systems
NASA Astrophysics Data System (ADS)
Eremin, A. S.; Olemskoy, I. V.
2016-06-01
We consider here numerical methods for systems of retarded functional differential equations of two equations in which the right-hand sides are cross-dependent of the unknown functions, i.e. the derivatives of unknowns don't depend on the same unknowns. It is shown that using the special structure of the system one can construct functional continuous methods of Runge-Kutta type with fewer stages than it is necessary in case of general Runge-Kutta functional continuous methods. Order conditions and example methods of orders three and four are presented. Test problems are solved, demonstrating the declared convergence order of the new methods.
Runge-Kutta Methods for Linear Ordinary Differential Equations
NASA Technical Reports Server (NTRS)
Zingg, David W.; Chisholm, Todd T.
1997-01-01
Three new Runge-Kutta methods are presented for numerical integration of systems of linear inhomogeneous ordinary differential equations (ODES) with constant coefficients. Such ODEs arise in the numerical solution of the partial differential equations governing linear wave phenomena. The restriction to linear ODEs with constant coefficients reduces the number of conditions which the coefficients of the Runge-Kutta method must satisfy. This freedom is used to develop methods which are more efficient than conventional Runge-Kutta methods. A fourth-order method is presented which uses only two memory locations per dependent variable, while the classical fourth-order Runge-Kutta method uses three. This method is an excellent choice for simulations of linear wave phenomena if memory is a primary concern. In addition, fifth- and sixth-order methods are presented which require five and six stages, respectively, one fewer than their conventional counterparts, and are therefore more efficient. These methods are an excellent option for use with high-order spatial discretizations.
Generalized disks of contractivity for explicit and implicit Runge-Kutta methods
NASA Technical Reports Server (NTRS)
Dahlquist, G.; Jeltsch, R.
1979-01-01
The A-contractivity of Runge-Kutta methods with respect to an inner product norm was investigated thoroughly by Butcher and Burrage (who used the term B-stability). Their theory is extended to contractivity in a region bounded by a circle through the origin. The largest possible circle is calculated for many known explicit Runge-Kutta methods. As a rule it is considerably smaller than the stability region, and in several cases it degenerates to a point. It is shown that an explicit Runge-Kutta method cannot be contractive in any circle of this class if it is more than fourth order accurate.
Two-derivative Runge-Kutta methods for differential equations
NASA Astrophysics Data System (ADS)
Chan, Robert P. K.; Wang, Shixiao; Tsai, Angela Y. J.
2012-09-01
Two-derivative Runge-Kutta (TDRK) methods are a special case of multi-derivative Runge-Kutta methods first studied by Kastlunger and Wanner [1, 2]. These methods incorporate derivatives of order higher than the first in their formulation but we consider only the first and second derivatives. In this paper we first present our study of both explicit [3] and implicit TDRK methods on stiff ODE problems. We then extend the applications of these TDRK methods to various partial differential equations [4]. In particular, we show how a 2-stage implicit TDRK method of order 4 and stage order 4 can be adapted to solve diffusion equations more efficiently than the popular Crank-Nicolson method.
A Low-Dispersion and Low-Dissipation Implicit Runge-Kutta Scheme
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
A Low-Dispersion and Low-Dissipation Implicit Runge-Kutta Scheme.
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
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.
Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations
NASA Technical Reports Server (NTRS)
Yu, S. T.; Tsai, Y.-L. P.; Hsieh, K. C.
1992-01-01
An investigation of the Runge-Kutta time-stepping, combined with compact difference schemes to solve the unsteady Euler equations, is presented. Initially, a generalized form of a N-step Runge-Kutta technique is derived. By comparing this generalized form with its Taylor's series counterpart, the criteria for the three-step and four-step schemes to be of third- and fourth-order accurate are obtained.
Low-dissipation and -dispersion Runge-Kutta schemes for computational acoustics
NASA Technical Reports Server (NTRS)
Hu, F. Q.; Hussaini, M. Y.; Manthey, J.
1994-01-01
In this paper, we investigate accurate and efficient time advancing methods for computational acoustics, where non-dissipative and non-dispersive properties are of critical importance. Our analysis pertains to the application of Runge-Kutta methods to high-order finite difference discretization. In many CFD applications multi-stage Runge-Kutta schemes have often been favored for their low storage requirements and relatively large stability limits. For computing acoustic waves, however, the stability consideration alone is not sufficient, since the Runge-Kutta schemes entail both dissipation and dispersion errors. The time step is now limited by the tolerable dissipation and dispersion errors in the computation. In the present paper, it is shown that if the traditional Runge-Kutta schemes are used for time advancing in acoustic problems, time steps greatly smaller than that allowed by the stability limit are necessary. Low-Dissipation and -Dispersion Runge-Kutta (LDDRE) schemes are proposed, based on an optimization that minimizes the dissipation and dispersion errors for wave propagation. Order optimizations of both single-step and two-step alternating schemes are considered. The proposed LDDRK schemes are remarkably more efficient than the classical Runge-Kutta schemes for acoustic computations. Moreover, low storage implementations of the optimized schemes are discussed. Special issues of implementing numerical boundary conditions in the LDDRK schemes are also addressed.
NASA Technical Reports Server (NTRS)
Hu, F. Q.; Hussaini, M. Y.; Manthey, J.
1995-01-01
We investigate accurate and efficient time advancing methods for computational aeroacoustics, where non-dissipative and non-dispersive properties are of critical importance. Our analysis pertains to the application of Runge-Kutta methods to high-order finite difference discretization. In many CFD applications, multi-stage Runge-Kutta schemes have often been favored for their low storage requirements and relatively large stability limits. For computing acoustic waves, however, the stability consideration alone is not sufficient, since the Runge-Kutta schemes entail both dissipation and dispersion errors. The time step is now limited by the tolerable dissipation and dispersion errors in the computation. In the present paper, it is shown that if the traditional Runge-Kutta schemes are used for time advancing in acoustic problems, time steps greatly smaller than that allowed by the stability limit are necessary. Low Dissipation and Dispersion Runge-Kutta (LDDRK) schemes are proposed, based on an optimization that minimizes the dissipation and dispersion errors for wave propagation. Optimizations of both single-step and two-step alternating schemes are considered. The proposed LDDRK schemes are remarkably more efficient than the classical Runge-Kutta schemes for acoustic computations. Numerical results of each Category of the Benchmark Problems are presented. Moreover, low storage implementations of the optimized schemes are discussed. Special issues of implementing numerical boundary conditions in the LDDRK schemes are also addressed.
Runge-Kutta based generalized convolution quadrature
NASA Astrophysics Data System (ADS)
Lopez-Fernandez, Maria; Sauter, Stefan
2016-06-01
We present the Runge-Kutta generalized convolution quadrature (gCQ) with variable time steps for the numerical solution of convolution equations for time and space-time problems. We present the main properties of the method and a convergence result.
Construction of IMEX methods with inherent Runge-Kutta stability
NASA Astrophysics Data System (ADS)
Braś, Michał; Izzo, Giuseppe; Jackiewicz, Zdzislaw
2016-06-01
We describe construction of implicit-explicit (IMEX) general linear methods (GLMs) with inherent Runge-Kutta stability (IRKS) for differential systems with non-stiff and stiff processes. We will use the extrapolation approach to remove implicitness in the non-stiff terms to compute unknown stage values in terms of stage derivatives at the previous step and those already computed in the current step. Highly stable IMEX GLMs of stage order equal to the order were derived up to the order four. These methods do not suffer from order reduction phenomenon which is confirmed by numerical experiments.
Accurate Monotonicity - Preserving Schemes With Runge-Kutta Time Stepping
NASA Technical Reports Server (NTRS)
Suresh, A.; Huynh, H. T.
1997-01-01
A new class of high-order monotonicity-preserving schemes for the numerical solution of conservation laws is presented. The interface value in these schemes is obtained by limiting a higher-order polynominal reconstruction. The limiting is designed to preserve accuracy near extrema and to work well with Runge-Kutta time stepping. Computational efficiency is enhanced by a simple test that determines whether the limiting procedure is needed. For linear advection in one dimension, these schemes are shown as well as the Euler equations also confirm their high accuracy, good shock resolution, and computational efficiency.
Scaled Runge-Kutta algorithms for treating the problem of dense output
NASA Technical Reports Server (NTRS)
Horn, M. K.
1982-01-01
A set of scaled Runge-Kutta algorithms for the third- through fifth-orders are developed to determine the solution at any point within the integration step at a relatively small increase in computing time. Each scaled algorithm is designed to be used with an existing Runge-Kutta formula, using the derivative evaluations of the defining algorithm along with an additional derivative evaluation (or two). Third-order, scaled algorithms are embedded within the existing formulas at no additional derivative expense. Such algorithms can easily be adopted to generate interpolating polynomials (or dependent variable stops) efficiently.
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.
Obtaining Runge-Kutta Solutions Between Time Steps
NASA Technical Reports Server (NTRS)
Horn, M. K.
1984-01-01
New interpolation method used with existing Runge-Kutta algorithms. Algorithm evaluates solution at intermediate point within integration step. Only few additional computations required to produce intermediate solution data. Runge-Kutta method provides accurate solution with larger time steps than allowable in other methods.
A Family of Trigonometrically-fitted Partitioned Runge-Kutta Symplectic Methods
Monovasilis, Th.; Kalogiratou, Z.; Simos, T. E.
2007-12-26
We are presenting a family of trigonometrically fitted partitioned Runge-Kutta symplectic methods of fourth order with six stages. The solution of the one dimensional time independent Schroedinger equation is considered by trigonometrically fitted symplectic integrators. The Schroedinger equation is first transformed into a Hamiltonian canonical equation. Numerical results are obtained for the one-dimensional harmonic oscillator and the exponential potential.
Amplification and Suppression of Round-Off Error in Runge-Kutta Methods
ERIC Educational Resources Information Center
Prentice, J. S. C.
2011-01-01
A simple nonstiff linear initial-value problem is used to demonstrate the amplification of round-off error in the course of using a second-order Runge-Kutta method. This amplification is understood in terms of an appropriate expression for the global error. An implicit method is then used to show how the roundoff error may actually be suppressed.…
On implicit Runge-Kutta methods for parallel computations
NASA Technical Reports Server (NTRS)
Keeling, Stephen L.
1987-01-01
Implicit Runge-Kutta methods which are well-suited for parallel computations are characterized. It is claimed that such methods are first of all, those for which the associated rational approximation to the exponential has distinct poles, and these are called multiply explicit (MIRK) methods. Also, because of the so-called order reduction phenomenon, there is reason to require that these poles be real. Then, it is proved that a necessary condition for a q-stage, real MIRK to be A sub 0-stable with maximal order q + 1 is that q = 1, 2, 3, or 5. Nevertheless, it is shown that for every positive integer q, there exists a q-stage, real MIRK which is I-stable with order q. Finally, some useful examples of algebraically stable MIRKs are given.
Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations
NASA Technical Reports Server (NTRS)
Kennedy, Christopher A.; Carpenter, Mark H.
2002-01-01
Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one- dimensional convection-diffusion-reaction (CDR) equations. Accuracy, stability, conservation, and dense-output are first considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, (N = 2), additive Runge-Kutta (ARK(sub 2)) methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms of the partitioned method are of equal order to those of the elemental methods. Derived ARK(sub 2) methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z['] yields -infinity, and retain high stability efficiency in the absence of stiffness, z['] yield 0. Extrapolation-type stage- value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK(sub 2) error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK(sub 2) methods.
Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations
NASA Technical Reports Server (NTRS)
Kennedy, Christopher A.; Carpenter, Mark H.
2001-01-01
Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one-dimensional convection-diffusion-reaction (CDR) equations. First, accuracy, stability, conservation, and dense output are considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, N = 2, additive Runge-Kutta ARK2 methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms are of equal order to those of the elemental methods. Derived ARK2 methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z(exp [I]) goes to infinity, and retain high stability efficiency in the absence of stiffness, z(exp [I]) goes to zero. Extrapolation-type stage-value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK2 error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK2 methods while offering error control and stage-value predictors.
NASA Astrophysics Data System (ADS)
Kalogiratou, Z.; Monovasilis, Th.; Psihoyios, G.; Simos, T. E.
2014-03-01
In this work we review single step methods of the Runge-Kutta type with special properties. Among them are methods specially tuned to integrate problems that exhibit a pronounced oscillatory character and such problems arise often in celestial mechanics and quantum mechanics. Symplectic methods, exponentially and trigonometrically fitted methods, minimum phase-lag and phase-fitted methods are presented. These are Runge-Kutta, Runge-Kutta-Nyström and Partitioned Runge-Kutta methods. The theory of constructing such methods is given as well as several specific methods. In order to present the performance of the methods we have tested 58 methods from all categories. We consider the two dimensional harmonic oscillator, the two body problem, the pendulum problem and the orbital problem studied by Stiefel and Bettis. Also we have tested the methods on the computation of the eigenvalues of the one dimensional time independent Schrödinger equation with the harmonic oscillator, the doubly anharmonic oscillator and the exponential potentials.
Multirate Runge-Kutta schemes for advection equations
NASA Astrophysics Data System (ADS)
Schlegel, Martin; Knoth, Oswald; Arnold, Martin; Wolke, Ralf
2009-04-01
Explicit time integration methods can be employed to simulate a broad spectrum of physical phenomena. The wide range of scales encountered lead to the problem that the fastest cell of the simulation dictates the global time step. Multirate time integration methods can be employed to alter the time step locally so that slower components take longer and fewer time steps, resulting in a moderate to substantial reduction of the computational cost, depending on the scenario to simulate [S. Osher, R. Sanders, Numerical approximations to nonlinear conservation laws with locally varying time and space grids, Math. Comput. 41 (1983) 321-336; H. Tang, G. Warnecke, A class of high resolution schemes for hyperbolic conservation laws and convection-diffusion equations with varying time and pace grids, SIAM J. Sci. Comput. 26 (4) (2005) 1415-1431; E. Constantinescu, A. Sandu, Multirate timestepping methods for hyperbolic conservation laws, SIAM J. Sci. Comput. 33 (3) (2007) 239-278]. In air pollution modeling the advection part is usually integrated explicitly in time, where the time step is constrained by a locally varying Courant-Friedrichs-Lewy (CFL) number. Multirate schemes are a useful tool to decouple different physical regions so that this constraint becomes a local instead of a global restriction. Therefore it is of major interest to apply multirate schemes to the advection equation. We introduce a generic recursive multirate Runge-Kutta scheme that can be easily adapted to an arbitrary number of refinement levels. It preserves the linear invariants of the system and is of third order accuracy when applied to certain explicit Runge-Kutta methods as base method.
NASA Technical Reports Server (NTRS)
Lear, W. M.
1974-01-01
The integration is discussed of the vector differential equation X = F(x, t) from time t sub i to t sub (i = 1) where only the values of x sub i are available for the the integration. No previous values of x or x prime are used. Using an orbit integration problem, comparisons are made between Taylor series integrators and various types and orders of Runge-Kutta integrators. A fourth order Runge-Kutta type integrator for orbital work is presented, and approximate (there may be no exact) fifth order Runge-Kutta integrators are discussed. Also discussed and compared is a self starting integrator ising delta f/delta x. A numerical method for controlling the accuracy of integration is given, and the special equations for accurately integrating accelerometer data are shown.
Galerkin/Runge-Kutta discretizations for parabolic equations with time-dependent coefficients
NASA Technical Reports Server (NTRS)
Keeling, Stephen L.
1989-01-01
A new class of fully discrete Galerkin/Runge-Kutta methods is constructed and analyzed for linear parabolic initial boundary value problems with time dependent coefficients. Unlike any classical counterpart, this class offers arbitrarily high order convergence while significantly avoiding what has been called order reduction. In support of this claim, error estimates are proved, and computational results are presented. Additionally, since the time stepping equations involve coefficient matrices changing at each time step, a preconditioned iterative technique is used to solve the linear systems only approximately. Nevertheless, the resulting algorithm is shown to preserve the original convergence rate while using only the order of work required by the base scheme applied to a linear parabolic problem with time independent coefficients. Furthermore, it is noted that special Runge-Kutta methods allow computations to be performed in parallel so that the final execution time can be reduced to that of a low order method.
Galerkin/Runge-Kutta discretizations for parabolic equations with time dependent coefficients
NASA Technical Reports Server (NTRS)
Keeling, Stephen L.
1987-01-01
A new class of fully discrete Galerkin/Runge-Kutta methods is constructed and analyzed for linear parabolic initial boundary value problems with time dependent coefficients. Unlike any classical counterpart, this class offers arbitrarily high order convergence while significantly avoiding what has been called order reduction. In support of this claim, error estimates are proved, and computational results are presented. Additionally, since the time stepping equations involve coefficient matrices changing at each time step, a preconditioned iterative technique is used to solve the linear systems only approximately. Nevertheless, the resulting algorithm is shown to preserve the original convergence rate while using only the order of work required by the base scheme applied to a linear parabolic problem with time independent coefficients. Furthermore, it is noted that special Runge-Kutta methods allow computations to be performed in parallel so that the final execution time can be reduced to that of a low order method.
A Runge-Kutta discontinuous finite element method for high speed flows
NASA Technical Reports Server (NTRS)
Bey, Kim S.; Oden, J. T.
1991-01-01
A Runge-Kutta discontinuous finite element method is developed for hyperbolic systems of conservation laws in two space variables. The discontinuous Galerkin spatial approximation to the conservation laws results in a system of ordinary differential equations which are marched in time using Runge-Kutta methods. Numerical results for the two-dimensional Burger's equation show that the method is (p+1)-order accurate in time and space, where p is the degree of the polynomial approximation of the solution within an element and is capable of capturing shocks over a single element without oscillations. Results for this problem also show that the accuracy of the solution in smooth regions is unaffected by the local projection and that the accuracy in smooth regions increases as p increases. Numerical results for the Euler equations show that the method captures shocks without oscillations and with higher resolution than a first-order scheme.
Tremblay, Jean Christophe; Carrington, Tucker Jr.
2004-12-15
If the Hamiltonian is time dependent it is common to solve the time-dependent Schroedinger equation by dividing the propagation interval into slices and using an (e.g., split operator, Chebyshev, Lanczos) approximate matrix exponential within each slice. We show that a preconditioned adaptive step size Runge-Kutta method can be much more efficient. For a chirped laser pulse designed to favor the dissociation of HF the preconditioned adaptive step size Runge-Kutta method is about an order of magnitude more efficient than the time sliced method.
A stiffly-stable implicit Runge-Kutta algorithm for CFD applications
NASA Technical Reports Server (NTRS)
Baker, A. J.; Iannelli, G. S.
1988-01-01
A stiffly-stable implicit Runge-Kutta integration algorithm is derived for CFD applications spanning the range of semidiscrete theories. The algorithm family contains the one-step 'theta' algorithms, including backwards Euler and the trapezoidal rule, and provides a versatile framework to identify expressions governing algorithm stability characteristics. Parameters of a Runge-Kutta optimal implicit algorithm, second-order accurate in time and stiffly-stable, are established. This algorithm is implemented within a weak statement finite element semidiscrete formulation for one- and two-dimensional conservation law systems. Numerical results are compared to theta-algorithm solutions, for unsteady quasi-one-dimensional Euler predictions with shocks, and for a specially derived two-dimensional conservation law system modeling the Euler equations.
NASA Technical Reports Server (NTRS)
Jameson, A.; Schmidt, Wolfgang; Turkel, Eli
1981-01-01
A new combination of a finite volume discretization in conjunction with carefully designed dissipative terms of third order, and a Runge Kutta time stepping scheme, is shown to yield an effective method for solving the Euler equations in arbitrary geometric domains. The method has been used to determine the steady transonic flow past an airfoil using an O mesh. Convergence to a steady state is accelerated by the use of a variable time step determined by the local Courant member, and the introduction of a forcing term proportional to the difference between the local total enthalpy and its free stream value.
An unconditionally stable Runge-Kutta method for unsteady flows
NASA Technical Reports Server (NTRS)
Jorgenson, Philip C. E.; Chima, Rodrick V.
1989-01-01
A quasi-three-dimensional analysis was developed for unsteady rotor-stator interaction in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream surface thickness. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. The equations are integrated in time using a four-stage Runge-Kutta scheme with a constant time step. Implicit residual smoothing was employed to accelerate the solution of the time accurate computations. The scheme is described and accuracy analyses are given. Results are shown for a supersonic through-flow fan designed for NASA Lewis. The rotor:stator blade ratio was taken as 1:1. Results are also shown for the first stage of the Space Shuttle Main Engine high pressure fuel turbopump. Here the blade ratio is 2:3. Implicit residual smoothing was used to increase the time step limit of the unsmoothed scheme by a factor of six with negligible differences in the unsteady results. It is felt that the implicitly smoothed Runge-Kutta scheme is easily competitive with implicit schemes for unsteady flows while retaining the simplicity of an explicit scheme.
An unconditionally stable Runge-Kutta method for unsteady flows
NASA Technical Reports Server (NTRS)
Jorgenson, Philip C. E.; Chima, Rodrick V.
1988-01-01
A quasi-three dimensional analysis was developed for unsteady rotor-stator interaction in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body fitted coordinate system. It accounts for the effects of rotation, radius change, and stream surface thickness. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. The equations are integrated in time using a four stage Runge-Kutta scheme with a constant time step. Implicit residual smoothing was employed to accelerate the solution of the time accurate computations. The scheme is described and accuracy analyses are given. Results are shown for a supersonic through-flow fan designed for NASA Lewis. The rotor:stator blade ratio was taken as 1:1. Results are also shown for the first stage of the Space Shuttle Main Engine high pressure fuel turbopump. Here the blade ratio is 2:3. Implicit residual smoothing was used to increase the time step limit of the unsmoothed scheme by a factor of six with negligible differences in the unsteady results. It is felt that the implicitly smoothed Runge-Kutta scheme is easily competitive with implicit schemes for unsteady flows while retaining the simplicity of an explicit scheme.
A computer program for determining truncation error coefficients for Runge-Kutta methods
NASA Technical Reports Server (NTRS)
Horn, M. K.
1980-01-01
The basic structure of a program to generate the truncation error coefficients for Runge-Kutta (RK) methods is reformulated to reduce storage requirements significantly and to accommodate variable dimensioning. This FORTRAN program, SUBROUTINE RKEQ, determines truncation error coefficients for RK algorithms for orders 1 through 10 and extends the order of coefficients through 12 with the 11th- and 12th-order terms determined following the patterns used to establish the lower order coefficients. Both subroutines (the original and RKEQ) are also written to treat RK m-fold methods which utilize m known derivatives of f to increase the order of the algorithm. Setting m = 0 gives the classical RK algorithm.
NASA Technical Reports Server (NTRS)
Subramanian, S. V.; Bozzola, R.
1987-01-01
Numerical solutions of the unsteady Euler equations are obtained using the classical fourth order Runge Kutta time marching scheme. This method is fully explicit and is applied to the governing equations in the finite volume, conservation law form. In order to determine the efficiency of this scheme for solving turbomachinery flows, steady blade-to-blade solutions are obtained for compressor and turbine cascades under subsonic and transonic flow conditions. Computed results are compared with other numerical methods and wind tunnel measurements. The study also focuses on other important numerical aspects influencing the performance of the algorithm and the solution accuracy such as grid types, boundary conditions and artificial viscosity. For this purpose, H, O, and C type computational grids as well as characteristic and extrapolation type boundary conditions are included in solution procedures.
NASA Technical Reports Server (NTRS)
Subramanian, S. V.; Bozzola, R.
1985-01-01
Numerical solutions of the unsteady Euler equations are obtained using the classical fourth order Runge Kutta time marching scheme. This method is fully explicit and is applied to the governing equations in the finite volume, conservation law form. In order to determine the efficiency of this scheme for solving turbomachinery flows, steady blade-to-blade solutions are obtained for compressor and turbine cascades under subsonic and transonic flow conditions. Computed results are compared with other numerical methods and wind tunnel measurements. The present study also focuses on other important numerical aspects influencing the performance of the algorithm and the solution accuracy such as grid types, boundary conditions, and artificial viscosity. For this purpose, H, O, and C type computational grids as well as characteristic and extrapolation type boundary conditions are included in the solution procedure.
An explicit Runge-Kutta method for turbulent reacting flow calculations
NASA Technical Reports Server (NTRS)
Boretti, A. A.
1989-01-01
The paper presents a numerical method for the solution of the conservation equations governing steady, reacting, turbulent viscous flow in two-dimensional geometries, in both Cartesian and axisymmetric coordinates. These equations are written in Favre-averaged form and closed with a first order model. A two-equation K-epsilon model, where low Reynolds number and compressibility effects are included, and a modified eddy-break up model are used to simulate fluid mechanics turbulence, chemistry and turbulence-combustion interaction. The solution is obtained by using a pseudo-unsteady method with improved perturbation propagation properties. The equations are discretized in space by using a finite volume formulation. An explicit multi-stage dissipative Runge-Kutta algorithm is then used to advance the flow equations in the pseudo-time. The method is applied to the computation of both diffusion and premixed turbulent reacting flows. The computed temperature distributions compare favorably with experimental data.
GPU acceleration of Runge Kutta-Fehlberg and its comparison with Dormand-Prince method
NASA Astrophysics Data System (ADS)
Seen, Wo Mei; Gobithaasan, R. U.; Miura, Kenjiro T.
2014-07-01
There is a significant reduction of processing time and speedup of performance in computer graphics with the emergence of Graphic Processing Units (GPUs). GPUs have been developed to surpass Central Processing Unit (CPU) in terms of performance and processing speed. This evolution has opened up a new area in computing and researches where highly parallel GPU has been used for non-graphical algorithms. Physical or phenomenal simulations and modelling can be accelerated through General Purpose Graphic Processing Units (GPGPU) and Compute Unified Device Architecture (CUDA) implementations. These phenomena can be represented with mathematical models in the form of Ordinary Differential Equations (ODEs) which encompasses the gist of change rate between independent and dependent variables. ODEs are numerically integrated over time in order to simulate these behaviours. The classical Runge-Kutta (RK) scheme is the common method used to numerically solve ODEs. The Runge Kutta Fehlberg (RKF) scheme has been specially developed to provide an estimate of the principal local truncation error at each step, known as embedding estimate technique. This paper delves into the implementation of RKF scheme for GPU devices and compares its result with Dorman Prince method. A pseudo code is developed to show the implementation in detail. Hence, practitioners will be able to understand the data allocation in GPU, formation of RKF kernels and the flow of data to/from GPU-CPU upon RKF kernel evaluation. The pseudo code is then written in C Language and two ODE models are executed to show the achievable speedup as compared to CPU implementation. The accuracy and efficiency of the proposed implementation method is discussed in the final section of this paper.
Flux-vector splitting and Runge-Kutta methods for the Euler equations
NASA Technical Reports Server (NTRS)
Turkel, E.; Vanleer, B.
1984-01-01
Runge-Kutta schemes have been used as a method of solving the Euler equations exterior to an airfoil. In the past this has been coupled with central differences and an artificial vesocity in space. In this study the Runge-Kutta time-stepping scheme is coupled with an upwinded space approximation based on flux-vector splitting. Several acceleration techniques are also considered including a local time step, residual smoothing and multigrid.
Rendering log aesthetic curves via Runge-Kutta method
NASA Astrophysics Data System (ADS)
Gobithaasan, R. U.; Meng, T. Y.; Piah, A. R. M.; Miura, K. T.
2014-07-01
Log Aesthetic Curves (LAC) are visually pleasing curves which has been developed using monotonic curvature profile. Hence, it can be easily implemented in product design environment, e.g, Rhino 3D CAD systems. LAC is generally represented in an integral form of its turning angle. Traditionally, Gaussian-Kronrod method has been used to render this curve which consumes less than one second for a given interval. Recently, Incomplete Gamma Function was proposed to represent LAC analytically which decreases the computation time up to 13 times. However, only certain value of shape parameters (denoted as α) which dictates the types of curves generated for LAC, can be used to compute LAC. In this paper, the classical Runge-Kutta (RK4) method is proposed to evaluate LAC numerically to reduce the LAC computation time for arbitrary, α. The preliminary result looks promising where the evaluation time is decreased tremendously. This paper also demonstrates the accuracy control of LAC by reducing the stepsize of RK4. The computation time and the accuracy for various α, are also illustrated in the last section of this paper.
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.
NASA Astrophysics Data System (ADS)
Kulikov, G. Yu.
2015-06-01
A technique for constructing nested implicit Runge-Kutta methods in the class of mono-implicit formulas of this type is studied. These formulas are highly efficient in practice, since the dimension of the original system of differential equations is preserved, which is not possible in the case of implicit multistage Runge-Kutta formulas of the general from. On the other hand, nested implicit Runge-Kutta methods inherit all major properties of general formulas of this form, such as A-stability, symmetry, and symplecticity in a certain sense. Moreover, they can have sufficiently high stage and classical orders and, without requiring high extra costs, can ensure dense output of integration results of the same accuracy as the order of the underlying method. Thus, nested methods are efficient when applied to the numerical integration of differential equations of various sorts, including stiff and nonstiff problems, Hamiltonian systems, and invertible equations. In this paper, previously proposed nested methods based on the Gauss quadrature formulas are generalized to Lobatto-type methods. Additionally, a unified technique for constructing all such methods is proposed. Its performance is demonstrated as applied to embedded examples of nested implicit formulas of various orders. All the methods constructed are supplied with tools for local error estimation and automatic variable-stepsize mesh generation based on an optimal stepsize selection. These numerical methods are verified by solving test problems with known solutions. Additionally, a comparative analysis of these methods with Matlab built-in solvers is presented.
NASA Astrophysics Data System (ADS)
Diele, F.; Marangi, C.; Ragni, S.
2009-08-01
Direct numerical approximation of a continuous-time infinite horizon control problem, requires to recast the model as a discrete-time, finite-horizon control model. The quality of the optimization results can be heavily degraded if the discretization process does not take into account features of the original model to be preserved. Restricting their attention to optimal growh problems with a steady state, Mercenier and Michel in [1] and [2], studied the conditions to be imposed for ensuring that discrete first-order approximation models have the same steady states as the infinite-horizon continuous-times counterpart. Here we show that Mercenier and Michel scheme is a first order partitioned Runge-Kutta method applied to the state-costate differential system which arises from the Pontryagin maximum principle. The main consequence is that it is possible to consider high order schemes which generalize that algorithm by preserving the steady-growth invariance of the solutions with respect to the discretization process. Numerical examples show the efficiency and accuracy of the proposed methods when applied to the classical Ramsey growth model.
Exponential Runge-Kutta integrators for modelling Predator-Prey interactions
NASA Astrophysics Data System (ADS)
Diele, F.; Marangi, C.; Ragni, S.
2012-09-01
Spatially explicit models consisting of reaction-diffusion partial differential equations are considered in order to model prey-predator interactions, since it is known that the role of spatial processes reveals of great interest in the study of the effects of habitat fragmentation on biodiversity. As almost all of the realistic models in biology, these models are nonlinear and their solution is not known in closed form. Our aim is approximating the solution itself by means of exponential Runge-Kutta integrators. Moreover, we apply the shift-and-invert Krylov approach in order to evaluate the entire functions needed for implementing the exponential method. This numerical procedure reveals to be very eff cient in avoiding numerical instability during the simulation, since it allows us to adopt high order in the accuracy. This work has received funding from the European Union's Seventh Framework Programme FP7/2007-2013, SPA.2010.1.1-04: "Stimulating the development of GMES services in specif c are", under grant agreement 263435, project title: Biodiversity Multi-Source Monitoring System:from Space To Species (BIOSOS) coordinated by CNR-ISSIA, Bari-Italy (http://www.biosos.eu).
Diagonally Implicit Runge-Kutta Methods for Ordinary Differential Equations. A Review
NASA Technical Reports Server (NTRS)
Kennedy, Christopher A.; Carpenter, Mark H.
2016-01-01
A review of diagonally implicit Runge-Kutta (DIRK) methods applied to rst-order ordinary di erential equations (ODEs) is undertaken. The goal of this review is to summarize the characteristics, assess the potential, and then design several nearly optimal, general purpose, DIRK-type methods. Over 20 important aspects of DIRKtype methods are reviewed. A design study is then conducted on DIRK-type methods having from two to seven implicit stages. From this, 15 schemes are selected for general purpose application. Testing of the 15 chosen methods is done on three singular perturbation problems. Based on the review of method characteristics, these methods focus on having a stage order of two, sti accuracy, L-stability, high quality embedded and dense-output methods, small magnitudes of the algebraic stability matrix eigenvalues, small values of aii, and small or vanishing values of the internal stability function for large eigenvalues of the Jacobian. Among the 15 new methods, ESDIRK4(3)6L[2]SA is recommended as a good default method for solving sti problems at moderate error tolerances.
Low-Storage, Explicit Runge-Kutta Schemes for the Compressible Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Kennedy, Chistopher A.; Carpenter, Mark H.; Lewis, R. Michael
1999-01-01
The derivation of storage explicit Runge-Kutta (ERK) schemes has been performed in the context of integrating the compressible Navier-Stokes equations via direct numerical simulation. Optimization of ERK methods is done across the broad range of properties, such as stability and accuracy efficiency, linear and nonlinear stability, error control reliability, step change stability, and dissipation/dispersion accuracy, subject to varying degrees of memory economization. Following van der Houwen and Wray, 16 ERK pairs are presented using from two to five registers of memory per equation, per grid point and having accuracies from third- to fifth-order. Methods have been assessed using the differential equation testing code DETEST, and with the 1D wave equation. Two of the methods have been applied to the DNS of a compressible jet as well as methane-air and hydrogen-air flames. Derived 3(2) and 4(3) pairs are competitive with existing full-storage methods. Although a substantial efficiency penalty accompanies use of two- and three-register, fifth-order methods, the best contemporary full-storage methods can be pearl), matched while still saving two to three registers of memory.
Modified Runge-Kutta methods for solving ODES. M.S. Thesis
NASA Technical Reports Server (NTRS)
Vanvu, T.
1981-01-01
A class of Runge-Kutta formulas is examined which permit the calculation of an accurate solution anywhere in the interval of integration. This is used in a code which seldom has to reject a step; rather it takes a reduced step if the estimated error is too large. The absolute stability implications of this are examined.
NASA Astrophysics Data System (ADS)
Langer, Stefan
2014-11-01
For unstructured finite volume methods an agglomeration multigrid with an implicit multistage Runge-Kutta method as a smoother is developed for solving the compressible Reynolds averaged Navier-Stokes (RANS) equations. The implicit Runge-Kutta method is interpreted as a preconditioned explicit Runge-Kutta method. The construction of the preconditioner is based on an approximate derivative. The linear systems are solved approximately with a symmetric Gauss-Seidel method. To significantly improve this solution method grid anisotropy is treated within the Gauss-Seidel iteration in such a way that the strong couplings in the linear system are resolved by tridiagonal systems constructed along these directions of strong coupling. The agglomeration strategy is adapted to this procedure by taking into account exactly these anisotropies in such a way that a directional coarsening is applied along these directions of strong coupling. Turbulence effects are included by a Spalart-Allmaras model, and the additional transport-type equation is approximately solved in a loosely coupled manner with the same method. For two-dimensional and three-dimensional numerical examples and a variety of differently generated meshes we show the wide range of applicability of the solution method. Finally, we exploit the GMRES method to determine approximate spectral information of the linearized RANS equations. This approximate spectral information is used to discuss and compare characteristics of multistage Runge-Kutta methods.
A Runge-Kutta discontinuous Galerkin approach to solve reactive flows: The hyperbolic operator
Billet, G.; Ryan, J.
2011-02-20
A Runge-Kutta discontinuous Galerkin method to solve the hyperbolic part of reactive Navier-Stokes equations written in conservation form is presented. Complex thermodynamics laws are taken into account. Particular care has been taken to solve the stiff gaseous interfaces correctly with no restrictive hypothesis. 1D and 2D test cases are presented.
NASA Technical Reports Server (NTRS)
Cockrell, C. R.
1989-01-01
Numerical solutions of the differential equation which describe the electric field within an inhomogeneous layer of permittivity, upon which a perpendicularly-polarized plane wave is incident, are considered. Richmond's method and the Runge-Kutta method are compared for linear and exponential profiles of permittivities. These two approximate solutions are also compared with the exact solutions.
NASA Astrophysics Data System (ADS)
Wang, Xiaoqiang; Ju, Lili; Du, Qiang
2016-07-01
The Willmore flow formulated by phase field dynamics based on the elastic bending energy model has been widely used to describe the shape transformation of biological lipid vesicles. In this paper, we develop and investigate some efficient and stable numerical methods for simulating the unconstrained phase field Willmore dynamics and the phase field Willmore dynamics with fixed volume and surface area constraints. The proposed methods can be high-order accurate and are completely explicit in nature, by combining exponential time differencing Runge-Kutta approximations for time integration with spectral discretizations for spatial operators on regular meshes. We also incorporate novel linear operator splitting techniques into the numerical schemes to improve the discrete energy stability. In order to avoid extra numerical instability brought by use of large penalty parameters in solving the constrained phase field Willmore dynamics problem, a modified augmented Lagrange multiplier approach is proposed and adopted. Various numerical experiments are performed to demonstrate accuracy and stability of the proposed methods.
NASA Technical Reports Server (NTRS)
Shih, C. C.
1973-01-01
A theoretical investigation of gas flow inside a multilayer insulation system has been made for the case of the broadside pumping process. A set of simultaneous first-order differential equations for the temperature and pressure of the gas mixture was obtained by considering the diffusion mechanism of the gas molecules through the perforations on the insulation layers. A modified Runge-Kutta method was used for numerical experiment. The numerical stability problem was investigated. It has been shown that when the relaxation time is small compared with the time period over which the gas properties change appreciably, the set of differential equations can be replaced by a set of algebraic equations for solution. Numerical examples were given, and comparisons with experimental data were made.
Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations
NASA Technical Reports Server (NTRS)
Yu, Sheng-Tao
1992-01-01
Recent development using compact difference schemes to solve the Navier-Stokes equations show spectral-like accuracy. A study was made of the numerical characteristics of various combinations of the Runge-Kutta (RK) methods and compact difference schemes to calculate the unsteady Euler equations. The accuracy of finite difference schemes is assessed based on the evaluations of dissipative error. The objectives are reducing the numerical damping and, at the same time, preserving numerical stability. While this approach has tremendous success solving steady flows, numerical characteristics of unsteady calculations remain largely unclear. For unsteady flows, in addition to the dissipative errors, phase velocity and harmonic content of the numerical results are of concern. As a result of the discretization procedure, the simulated unsteady flow motions actually propagate in a dispersive numerical medium. Consequently, the dispersion characteristics of the numerical schemes which relate the phase velocity and wave number may greatly impact the numerical accuracy. The aim is to assess the numerical accuracy of the simulated results. To this end, the Fourier analysis is to provide the dispersive correlations of various numerical schemes. First, a detailed investigation of the existing RK methods is carried out. A generalized form of an N-step RK method is derived. With this generalized form, the criteria are derived for the three and four-step RK methods to be third and fourth-order time accurate for the non-linear equations, e.g., flow equations. These criteria are then applied to commonly used RK methods such as Jameson's 3-step and 4-step schemes and Wray's algorithm to identify the accuracy of the methods. For the spatial discretization, compact difference schemes are presented. The schemes are formulated in the operator-type to render themselves suitable for the Fourier analyses. The performance of the numerical methods is shown by numerical examples. These examples
NASA Astrophysics Data System (ADS)
Huang, Yong; Shi, Guo-Dong; Zhu, Ke-Yong
2016-06-01
This paper adopts the Runge-Kutta ray tracing method to obtain the ray-trajectory numerical solution in a two-dimensional gradient index medium. The emitting, absorbing and scattering processes are simulated by the Monte Carlo method. The temperature field and ray trajectory in the medium are obtained by the three methods, the Runge-Kutta ray tracing method, the ray tracing method with the cell model and the discrete curved ray tracing method with the linear refractive index cell model. Comparing the results of the three methods, it is found that the results by the Monte Carlo Runge-Kutta ray tracing method are of the highest accuracy. To improve the computational speed, the variable step-size Runge-Kutta ray tracing method is proposed, and the maximum relative error between the temperature field in the nonscattering medium by this method and the benchmark solution is less than 0.5%. The results also suggest that the Runge-Kutta ray tracing method would make the radiative transfer solution in the three-dimensional graded index media much easier.
An unconditionally stable Runge-Kutta method for unsteady rotor-stator interaction
NASA Technical Reports Server (NTRS)
Chima, Rodrick V.; Jorgenson, Philip C. E.
1989-01-01
A quasi-three-dimensional analysis has been developed for unsteady rotor-stator interaction in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stress-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations are integrated in time using an explicit four-stage Runge-Kutta scheme with a constant time step. Implicit residual smoothing is used to increase the stability limit of the time-accurate computations. The scheme is described, and stability and accuracy analyses are given.
Application of Runge-Kutta scheme for high-speed inviscid internal flows
NASA Technical Reports Server (NTRS)
Moitra, A.; Turkel, E.; Kumar, A.
1986-01-01
A multi-stage Runge-Kutta method is analyzed for solving the two-dimensional Euler equations for external and internal flow problems. Subsonic, supersonic and, highly supersonic flows are studied. Various techniques for accelerating the convergence to a steady state are described and analyzed. Effects of the grid aspect ratio on the rate of convergence are evaluated. An enthalpy damping technique applicable to supersonic flows is described in detail. Numerical results for supersonic flows containing both oblique and normal shocks are presented confirming the efficiency of the method.
Application of a Runge-Kutta scheme for high-speed inviscid internal flows
NASA Technical Reports Server (NTRS)
Moitra, A.; Turkel, E.; Kumar, A.
1986-01-01
A multi-stage Runge-Kutta method is analyzed for solving the two-dimensional Euler equations for external and internal flow problems. Subsonic, supersonic and, highly supersonic flows are studied. Various techniques for accelerating the convergence to a steady state are described and analyzed. Effects of the grid aspect ratio on the rate of convergence are evaluated. An enthalpy damping technique applicable to supersonic flows is described in detail. Numerical results for supersonic flows containing both oblique and normal shocks are presented confirming the efficiency of the method.
Convergence Acceleration of Runge-Kutta Schemes for Solving the Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Swanson, Roy C., Jr.; Turkel, Eli; Rossow, C.-C.
2007-01-01
The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 can be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. This RK/implicit scheme is used as a smoother for multigrid. Fourier analysis is applied to determine damping properties. Numerical dissipation operators based on the Roe scheme, a matrix dissipation, and the CUSP scheme are considered in evaluating the RK/implicit scheme. In addition, the effect of the number of RK stages is examined. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. Turbulent flows over an airfoil and wing at subsonic and transonic conditions are computed. The effects of the cell aspect ratio on convergence are investigated for Reynolds numbers between 5:7 x 10(exp 6) and 100 x 10(exp 6). It is demonstrated that the implicit preconditioner can reduce the computational time of a well-tuned standard RK scheme by a factor between four and ten.
Elkina, N V; Fedotov, A M; Herzing, C; Ruhl, H
2014-05-01
The Landau-Lifshitz equation provides an efficient way to account for the effects of radiation reaction without acquiring the nonphysical solutions typical for the Lorentz-Abraham-Dirac equation. We solve the Landau-Lifshitz equation in its covariant four-vector form in order to control both the energy and momentum of radiating particles. Our study reveals that implicit time-symmetric collocation methods of the Runge-Kutta-Nyström type are superior in accuracy and better at maintaining the mass-shell condition than their explicit counterparts. We carry out an extensive study of numerical accuracy by comparing the analytical and numerical solutions of the Landau-Lifshitz equation. Finally, we present the results of the simulation of particle scattering by a focused laser pulse. Due to radiation reaction, particles are less capable of penetrating into the focal region compared to the case where radiation reaction is neglected. Our results are important for designing forthcoming experiments with high intensity laser fields. PMID:25353922
An explicit Runge-Kutta method for unsteady rotor/stator interaction
NASA Technical Reports Server (NTRS)
Jorgenson, Philip C. E.; Chima, Rodrick V.
1988-01-01
A quasi-three-dimensional rotor/stator analysis has been developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations are integrated in time using a four-stage Runge-Kutta scheme with a constant timestep. Results are shown for the first stage of the Space Shuttle Main Engine high pressure fuel turbopump. Euler and Navier-Stokes results are compared on the scaled single- and multi-passage machine. The method is relatively fast and the quasi-three-dimensional formulation is applicable to a wide range of turbomachinery geometries.
Steady and Unsteady Numerical Solution of Generalized Newtonian Fluids Flow by Runge-Kutta method
NASA Astrophysics Data System (ADS)
Keslerová, R.; Kozel, K.; Prokop, V.
2010-09-01
In this paper the laminar viscous incompressible flow for generalized Newtonian (Newtonian and non-Newtonian) fluids is considered. The governing system of equations is the system of Navier-Stokes equations and the continuity equation. The steady and unsteady numerical solution for this system is computed by finite volume method combined with an artificial compressibility method. For time discretization the explicit multistage Runge-Kutta numerical scheme is considered. Steady state solution is achieved for t→∞ using steady boundary conditions and followed by steady residual behavior. The dual time-stepping method is considered for unsteady computation. The high artificial compressibility coefficient is used in the artificial compressibility method applied in the dual time τ. The steady and unsteady numerical results of Newtonian and non-Newtonian (shear thickening and shear thinning) fluids flow in the branching channel are presented.
Recent advances in Runge-Kutta schemes for solving 3-D Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Vatsa, Veer N.; Wedan, Bruce W.; Abid, Ridha
1989-01-01
A thin-layer Navier-Stokes has been developed for solving high Reynolds number, turbulent flows past aircraft components under transonic flow conditions. The computer code has been validated through data comparisons for flow past isolated wings, wing-body configurations, prolate spheroids and wings mounted inside wind-tunnels. The basic code employs an explicit Runge-Kutta time-stepping scheme to obtain steady state solution to the unsteady governing equations. Significant gain in the efficiency of the code has been obtained by implementing a multigrid acceleration technique to achieve steady-state solutions. The improved efficiency of the code has made it feasible to conduct grid-refinement and turbulence model studies in a reasonable amount of computer time. The non-equilibrium turbulence model of Johnson and King has been extended to three-dimensional flows and excellent agreement with pressure data has been obtained for transonic separated flow over a transport type of wing.
On spurious steady-state solutions of explicit Runge-Kutta schemes
NASA Technical Reports Server (NTRS)
Sweby, P. K.; Yee, H. C.; Griffiths, D. F.
1990-01-01
The bifurcation diagram associated with the logistic equation v sup n+1 = av sup n (1-v sup n) is by now well known, as is its equivalence to solving the ordinary differential equation u prime = alpha u (1-u) by the explicit Euler difference scheme. It has also been noted by Iserles that other popular difference schemes may not only exhibit period doubling and chaotic phenomena but also possess spurious fixed points. Runge-Kutta schemes applied to both the equation u prime = alpha u (1-u) and the cubic equation u prime = alpha u (1-u)(b-u) were studied computationally and analytically and their behavior was contrasted with the explicit Euler scheme. Their spurious fixed points and periodic orbits were noted. In particular, it was observed that these may appear below the linearized stability limits of the scheme and, consequently, computation may lead to erroneous results.
Runge-Kutta model-based nonlinear observer for synchronization and control of chaotic systems.
Beyhan, Selami
2013-07-01
This paper proposes a novel nonlinear gradient-based observer for synchronization and observer-based control of chaotic systems. The model is based on a Runge-Kutta model of the chaotic system where the evolution of the states or parameters is derived based on the error-square minimization. The stability and convergence conditions of observer and control methods are analyzed using a Lyapunov stability approach. In numerical simulations, the proposed observer and well-known sliding-mode observer are compared for the synchronization of a Lü chaotic system and observer-based stabilization of a Chen chaotic system. The noisy case for synchronization and parameter uncertainty case for stabilization are also considered for both observer-based methods. PMID:23672740
An explicit Runge-Kutta method for 3D turbulent incompressible flows
NASA Technical Reports Server (NTRS)
Sung, Chao-Ho; Lin, Cheng-Wen; Hung, C. M.
1988-01-01
A computer code has been developed to solve for the steady-state solution of the 3D incompressible Reynolds-averaged Navier-Stokes equations. The approach is based on the cell-center, central-difference, finite-volume formulation and an explicit one-step, multistage Runge-Kutta time-stepping scheme. The Baldwin-Lomax turbulence model is used. Techniques to accelerate the rate of convergence to a steady-state solution include the preconditioned method, the local time stepping, and the implicit residual smoothing. Improvements in computational efficiency have been demonstrated in several areas. This numerical procedure has been used to simulate the turbulent horseshoe vortex flow around an airfoil/flat-plate juncture.
NASA Technical Reports Server (NTRS)
Kanevsky, Alex
2004-01-01
My goal is to develop and implement efficient, accurate, and robust Implicit-Explicit Runge-Kutta (IMEX RK) methods [9] for overcoming geometry-induced stiffness with applications to computational electromagnetics (CEM), computational fluid dynamics (CFD) and computational aeroacoustics (CAA). IMEX algorithms solve the non-stiff portions of the domain using explicit methods, and isolate and solve the more expensive stiff portions using implicit methods. Current algorithms in CEM can only simulate purely harmonic (up to lOGHz plane wave) EM scattering by fighter aircraft, which are assumed to be pure metallic shells, and cannot handle the inclusion of coatings, penetration into and radiation out of the aircraft. Efficient MEX RK methods could potentially increase current CEM capabilities by 1-2 orders of magnitude, allowing scientists and engineers to attack more challenging and realistic problems.
NASA Astrophysics Data System (ADS)
Igumnov, Leonid; Ipatov, Aleksandr; Belov, Aleksandr; Petrov, Andrey
2015-09-01
The report presents the development of the time-boundary element methodology and a description of the related software based on a stepped method of numerical inversion of the integral Laplace transform in combination with a family of Runge-Kutta methods for analyzing 3-D mixed initial boundary-value problems of the dynamics of inhomogeneous elastic and poro-elastic bodies. The results of the numerical investigation are presented. The investigation methodology is based on direct-approach boundary integral equations of 3-D isotropic linear theories of elasticity and poroelasticity in Laplace transforms. Poroelastic media are described using Biot models with four and five base functions. With the help of the boundary-element method, solutions in time are obtained, using the stepped method of numerically inverting Laplace transform on the nodes of Runge-Kutta methods. The boundary-element method is used in combination with the collocation method, local element-by-element approximation based on the matched interpolation model. The results of analyzing wave problems of the effect of a non-stationary force on elastic and poroelastic finite bodies, a poroelastic half-space (also with a fictitious boundary) and a layered half-space weakened by a cavity, and a half-space with a trench are presented. Excitation of a slow wave in a poroelastic medium is studied, using the stepped BEM-scheme on the nodes of Runge-Kutta methods.
Runge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured meshes
NASA Astrophysics Data System (ADS)
Zhu, Jun; Zhong, Xinghui; Shu, Chi-Wang; Qiu, Jianxian
2013-09-01
In this paper we generalize a new type of limiters based on the weighted essentially non-oscillatory (WENO) finite volume methodology for the Runge-Kutta discontinuous Galerkin (RKDG) methods solving nonlinear hyperbolic conservation laws, which were recently developed in [32] for structured meshes, to two-dimensional unstructured triangular meshes. The key idea of such limiters is to use the entire polynomials of the DG solutions from the troubled cell and its immediate neighboring cells, and then apply the classical WENO procedure to form a convex combination of these polynomials based on smoothness indicators and nonlinear weights, with suitable adjustments to guarantee conservation. The main advantage of this new limiter is its simplicity in implementation, especially for the unstructured meshes considered in this paper, as only information from immediate neighbors is needed and the usage of complicated geometric information of the meshes is largely avoided. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good performance of this procedure.
Navier-Stokes calculations for DFVLR F5-wing in wind tunnel using Runge-Kutta time-stepping scheme
NASA Technical Reports Server (NTRS)
Vatsa, V. N.; Wedan, B. W.
1988-01-01
A three-dimensional Navier-Stokes code using an explicit multistage Runge-Kutta type of time-stepping scheme is used for solving the transonic flow past a finite wing mounted inside a wind tunnel. Flow past the same wing in free air was also computed to assess the effect of wind-tunnel walls on such flows. Numerical efficiency is enhanced through vectorization of the computer code. A Cyber 205 computer with 32 million words of internal memory was used for these computations.
Cui, Hengfei; Wang, Desheng; Wan, Min; Zhang, Jun-Mei; Zhao, Xiaodan; Tan, Ru San; Huang, Weimin; Xiong, Wei; Duan, Yuping; Zhou, Jiayin; Luo, Tong; Kassab, Ghassan S; Zhong, Liang
2016-06-01
The CT angiography (CTA) is a clinically indicated test for the assessment of coronary luminal stenosis that requires centerline extractions. There is currently no centerline extraction algorithm that is automatic, real-time and very accurate. Therefore, we sought to (i) develop a hybrid approach by incorporating fast marching and Runge-Kutta based methods for the extraction of coronary artery centerlines from CTA; (ii) evaluate the accuracy of the present method compared to Van's method by using ground truth centerline as a reference; (iii) evaluate the coronary lumen area of our centerline method in comparison with the intravascular ultrasound (IVUS) as the standard of reference. The proposed method was found to be more computationally efficient, and performed better than the Van's method in terms of overlap measures (i.e., OV: [Formula: see text] vs. [Formula: see text]; OF: [Formula: see text] vs. [Formula: see text]; and OT: [Formula: see text] vs. [Formula: see text], all [Formula: see text]). In comparison with IVUS derived coronary lumen area, the proposed approach was more accurate than the Van's method. This hybrid approach by incorporating fast marching and Runge-Kutta based methods could offer fast and accurate extraction of centerline as well as the lumen area. This method may garner wider clinical potential as a real-time coronary stenosis assessment tool. PMID:27140197
NASA Technical Reports Server (NTRS)
Carpenter, Mark H.; Gottlieb, David; Abarbanel, Saul; Don, Wai-Sun
1993-01-01
The conventional method of imposing time dependent boundary conditions for Runge-Kutta (RK) time advancement reduces the formal accuracy of the space-time method to first order locally, and second order globally, independently of the spatial operator. This counter intuitive result is analyzed in this paper. Two methods of eliminating this problem are proposed for the linear constant coefficient case: (1) impose the exact boundary condition only at the end of the complete RK cycle, (2) impose consistent intermediate boundary conditions derived from the physical boundary condition and its derivatives. The first method, while retaining the RK accuracy in all cases, results in a scheme with much reduced CFL condition, rendering the RK scheme less attractive. The second method retains the same allowable time step as the periodic problem. However it is a general remedy only for the linear case. For non-linear hyperbolic equations the second method is effective only for for RK schemes of third order accuracy or less. Numerical studies are presented to verify the efficacy of each approach.
NASA Astrophysics Data System (ADS)
Korneev, B. A.; Levchenko, V. D.
2016-03-01
In this paper we present the Runge-Kutta discontinuous Galerkin method (RKDG method) for the numerical solution of the Euler equations of gas dynamics. The method is being tested on a series of Riemann problems in the one-dimensional case. For the implementation of the method in the three-dimensional case, a DiamondTorre algorithm is proposed. It belongs to the class of the locally recursive non-locally asynchronous algorithms (LRnLA). With the help of this algorithm a significant increase of speed of calculations is achieved. As an example of the three-dimensional computing, a problem of the interaction of a bubble with a shock wave is considered.
Analysis of numerical stability and amplification matrices: Fourth-order Runge-Kutta methods
NASA Technical Reports Server (NTRS)
Kennedy, E. W.
1979-01-01
Amplification matrices, numerical kernels, stable, and exponentially stable numerical solutions are examined. The various techniques involved in these concepts are applied to certain systems that have Jordan forms, which are nondiagonal, with particular interest in the case of imaginary or zero eigenvalues.
Xu, Zhiliang; Chen, Xu-Yan; Liu, Yingjie
2014-12-01
We present a new formulation of the Runge-Kutta discontinuous Galerkin (RKDG) method [9, 8, 7, 6] for solving conservation Laws with increased CFL numbers. The new formulation requires the computed RKDG solution in a cell to satisfy additional conservation constraint in adjacent cells and does not increase the complexity or change the compactness of the RKDG method. Numerical computations for solving one-dimensional and two-dimensional scalar and systems of nonlinear hyperbolic conservation laws are performed with approximate solutions represented by piecewise quadratic and cubic polynomials, respectively. The hierarchical reconstruction [17, 33] is applied as a limiter to eliminate spurious oscillations in discontinuous solutions. From both numerical experiments and the analytic estimate of the CFL number of the newly formulated method, we find that: 1) this new formulation improves the CFL number over the original RKDG formulation by at least three times or more and thus reduces the overall computational cost; and 2) the new formulation essentially does not compromise the resolution of the numerical solutions of shock wave problems compared with ones computed by the RKDG method. PMID:25414520
Xu, Zhiliang; Chen, Xu-Yan; Liu, Yingjie
2014-01-01
We present a new formulation of the Runge-Kutta discontinuous Galerkin (RKDG) method [9, 8, 7, 6] for solving conservation Laws with increased CFL numbers. The new formulation requires the computed RKDG solution in a cell to satisfy additional conservation constraint in adjacent cells and does not increase the complexity or change the compactness of the RKDG method. Numerical computations for solving one-dimensional and two-dimensional scalar and systems of nonlinear hyperbolic conservation laws are performed with approximate solutions represented by piecewise quadratic and cubic polynomials, respectively. The hierarchical reconstruction [17, 33] is applied as a limiter to eliminate spurious oscillations in discontinuous solutions. From both numerical experiments and the analytic estimate of the CFL number of the newly formulated method, we find that: 1) this new formulation improves the CFL number over the original RKDG formulation by at least three times or more and thus reduces the overall computational cost; and 2) the new formulation essentially does not compromise the resolution of the numerical solutions of shock wave problems compared with ones computed by the RKDG method. PMID:25414520
Chen, Zhaoxia; Li, Juan; Zhang, Ruqiang; You, Xiong
2015-01-01
Oscillation is one of the most important phenomena in the chemical reaction systems in living cells. The general purpose simulation algorithms fail to take into account this special character and produce unsatisfying results. In order to enhance the accuracy of the integrator, the second-order derivative is incorporated in the scheme. The oscillatory feature of the solution is captured by the integrators with an exponential fitting property. Three practical exponentially fitted TDRK (EFTDRK) methods are derived. To test the effectiveness of the new EFTDRK methods, the two-gene system with cross-regulation and the circadian oscillation of the period protein in Drosophila are simulated. Each EFTDRK method has the best fitting frequency which minimizes the global error. The numerical results show that the new EFTDRK methods are more accurate and more efficient than their prototype TDRK methods or RK methods of the same order and the traditional exponentially fitted RK method in the literature. PMID:26633991
Chen, Zhaoxia; Li, Juan; Zhang, Ruqiang; You, Xiong
2015-01-01
Oscillation is one of the most important phenomena in the chemical reaction systems in living cells. The general purpose simulation algorithms fail to take into account this special character and produce unsatisfying results. In order to enhance the accuracy of the integrator, the second-order derivative is incorporated in the scheme. The oscillatory feature of the solution is captured by the integrators with an exponential fitting property. Three practical exponentially fitted TDRK (EFTDRK) methods are derived. To test the effectiveness of the new EFTDRK methods, the two-gene system with cross-regulation and the circadian oscillation of the period protein in Drosophila are simulated. Each EFTDRK method has the best fitting frequency which minimizes the global error. The numerical results show that the new EFTDRK methods are more accurate and more efficient than their prototype TDRK methods or RK methods of the same order and the traditional exponentially fitted RK method in the literature. PMID:26633991
NASA Technical Reports Server (NTRS)
Abarbanel, Saul; Gottlieb, David; Carpenter, Mark H.
1994-01-01
It has been previously shown that the temporal integration of hyperbolic partial differential equations (PDE's) may, because of boundary conditions, lead to deterioration of accuracy of the solution. A procedure for removal of this error in the linear case has been established previously. In the present paper we consider hyperbolic (PDE's) (linear and non-linear) whose boundary treatment is done via the SAT-procedure. A methodology is present for recovery of the full order of accuracy, and has been applied to the case of a 4th order explicit finite difference scheme.
NASA Technical Reports Server (NTRS)
Cannizzaro, Frank E.; Ash, Robert L.
1992-01-01
A state-of-the-art computer code has been developed that incorporates a modified Runge-Kutta time integration scheme, upwind numerical techniques, multigrid acceleration, and multi-block capabilities (RUMM). A three-dimensional thin-layer formulation of the Navier-Stokes equations is employed. For turbulent flow cases, the Baldwin-Lomax algebraic turbulence model is used. Two different upwind techniques are available: van Leer's flux-vector splitting and Roe's flux-difference splitting. Full approximation multi-grid plus implicit residual and corrector smoothing were implemented to enhance the rate of convergence. Multi-block capabilities were developed to provide geometric flexibility. This feature allows the developed computer code to accommodate any grid topology or grid configuration with multiple topologies. The results shown in this dissertation were chosen to validate the computer code and display its geometric flexibility, which is provided by the multi-block structure.
Multi-Dimensional Asymptotically Stable 4th Order Accurate Schemes for the Diffusion Equation
NASA Technical Reports Server (NTRS)
Abarbanel, Saul; Ditkowski, Adi
1996-01-01
An algorithm is presented which solves the multi-dimensional diffusion equation on co mplex shapes to 4th-order accuracy and is asymptotically stable in time. This bounded-error result is achieved by constructing, on a rectangular grid, a differentiation matrix whose symmetric part is negative definite. The differentiation matrix accounts for the Dirichlet boundary condition by imposing penalty like terms. Numerical examples in 2-D show that the method is effective even where standard schemes, stable by traditional definitions fail.
Computational aspects of the nonlinear normal mode initialization of the GLAS 4th order GCM
NASA Technical Reports Server (NTRS)
Navon, I. M.; Bloom, S. C.; Takacs, L.
1984-01-01
Using the normal modes of the GLAS 4th Order Model, a Machenhauer nonlinear normal mode initialization (NLNMI) was carried out for the external vertical mode using the GLAS 4th Order shallow water equations model for an equivalent depth corresponding to that associated with the external vertical mode. A simple procedure was devised which was directed at identifying computational modes by following the rate of increase of BAL sub M, the partial (with respect to the zonal wavenumber m) sum of squares of the time change of the normal mode coefficients (for fixed vertical mode index) varying over the latitude index L of symmetric or antisymmetric gravity waves. A working algorithm is presented which speeds up the convergence of the iterative Machenhauer NLNMI. A 24 h integration using the NLNMI state was carried out using both Matsuno and leap-frog time-integration schemes; these runs were then compared to a 24 h integration starting from a non-initialized state. The maximal impact of the nonlinear normal mode initialization was found to occur 6-10 hours after the initial time.
FAST DISPLACEMENT PROBABILITY PROFILE APPROXIMATION FROM HARDI USING 4TH-ORDER TENSORS.
Barmpoutis, Angelos; Vemuri, Baba C; Forder, John R
2008-05-14
Cartesian tensor basis have been widely used to approximate spherical functions. In Medical Imaging, tensors of various orders have been used to model the diffusivity function in Diffusion-weighted MRI data sets. However, it is known that the peaks of the diffusivity do not correspond to orientations of the underlying fibers and hence the displacement probability profiles should be employed instead. In this paper, we present a novel representation of the probability profile by a 4(th) order tensor, which is a smooth spherical function that can approximate single-fibers as well as multiple-fiber structures. We also present a method for efficiently estimating the unknown tensor coefficients of the probability profile directly from a given high-angular resolution diffusion-weighted (HARDI) data set. The accuracy of our model is validated by experiments on synthetic and real HARDI datasets from a fixed rat spinal cord. PMID:20046536
Continuum Kinetic Plasma Modeling Using a Conservative 4th-Order Method with AMR
NASA Astrophysics Data System (ADS)
Vogman, Genia; Colella, Phillip
2012-10-01
When the number of particles in a Debye sphere is large, a plasma can be accurately represented by a distribution function, which can be treated as a continuous incompressible fluid in phase space. In the most general case the evolution of such a distribution function is described by the 6D Boltzmann-Maxwell partial differential equation system. To address the challenges associated with solving a 6D hyperbolic governing equation, a simpler 3D Vlasov-Poisson system is considered. A 4th-order accurate Vlasov-Poisson model has been developed in one spatial and two velocity dimensions. The governing equation is cast in conservation law form and is solved with a finite volume representation. Adaptive mesh refinement (AMR) is used to allow for efficient use of computational resources while maintaining desired levels of resolution. The model employs a flux limiter to remedy non-physical effects such as numerical dispersion. The model is tested on the two-stream, beam-plasma, and Dory-Guest-Harris instabilities. All results are compared with linear theory.
Numerical integration of second order differential equations
NASA Technical Reports Server (NTRS)
Shanks, E. B.
1971-01-01
Performance characteristics of higher order approximations of Runge-Kutta type are analyzed, and performance predictors for time required on machine and for error size are developed. Technique is useful in evaluating system performance, analyzing material characteristics, and designing inertial guidance and nuclear instrumentation and materials.
NASA Technical Reports Server (NTRS)
Zingg, D. W.; Lomax, H.
1993-01-01
A six-stage low-storage Runge-Kutta time-marching method is presented and shown to be an efficient method for use with high-accuracy spatial difference operators for wave propagation problems. The accuracy of the method for inhomogeneous ordinary differential equations is demonstrated through numerical solutions of the linear convection equation with forced boundary conditions. Numerical experiments are presented simulating a sine wave and a Gaussian pulse propagating into and through the domain. For practical levels of mesh refinement corresponding to roughly ten points per wavelength, the six-stage Runge-Kutta method is more accurate than the popular fourth-order Runge-Kutta method. Further numerical experiments are presented which show that the numerical boundary scheme at an inflow boundary can be a significant source of error when high-accuracy spatial discretizations are used.
High-Order Residual-Distribution Hyperbolic Advection-Diffusion Schemes: 3rd-, 4th-, and 6th-Order
NASA Technical Reports Server (NTRS)
Mazaheri, Alireza R.; Nishikawa, Hiroaki
2014-01-01
In this paper, spatially high-order Residual-Distribution (RD) schemes using the first-order hyperbolic system method are proposed for general time-dependent advection-diffusion problems. The corresponding second-order time-dependent hyperbolic advection- diffusion scheme was first introduced in [NASA/TM-2014-218175, 2014], where rapid convergences over each physical time step, with typically less than five Newton iterations, were shown. In that method, the time-dependent hyperbolic advection-diffusion system (linear and nonlinear) was discretized by the second-order upwind RD scheme in a unified manner, and the system of implicit-residual-equations was solved efficiently by Newton's method over every physical time step. In this paper, two techniques for the source term discretization are proposed; 1) reformulation of the source terms with their divergence forms, and 2) correction to the trapezoidal rule for the source term discretization. Third-, fourth, and sixth-order RD schemes are then proposed with the above techniques that, relative to the second-order RD scheme, only cost the evaluation of either the first derivative or both the first and the second derivatives of the source terms. A special fourth-order RD scheme is also proposed that is even less computationally expensive than the third-order RD schemes. The second-order Jacobian formulation was used for all the proposed high-order schemes. The numerical results are then presented for both steady and time-dependent linear and nonlinear advection-diffusion problems. It is shown that these newly developed high-order RD schemes are remarkably efficient and capable of producing the solutions and the gradients to the same order of accuracy of the proposed RD schemes with rapid convergence over each physical time step, typically less than ten Newton iterations.
A Very High Order, Adaptable MESA Implementation for Aeroacoustic Computations
NASA Technical Reports Server (NTRS)
Dydson, Roger W.; Goodrich, John W.
2000-01-01
Since computational efficiency and wave resolution scale with accuracy, the ideal would be infinitely high accuracy for problems with widely varying wavelength scales. Currently, many of the computational aeroacoustics methods are limited to 4th order accurate Runge-Kutta methods in time which limits their resolution and efficiency. However, a new procedure for implementing the Modified Expansion Solution Approximation (MESA) schemes, based upon Hermitian divided differences, is presented which extends the effective accuracy of the MESA schemes to 57th order in space and time when using 128 bit floating point precision. This new approach has the advantages of reducing round-off error, being easy to program. and is more computationally efficient when compared to previous approaches. Its accuracy is limited only by the floating point hardware. The advantages of this new approach are demonstrated by solving the linearized Euler equations in an open bi-periodic domain. A 500th order MESA scheme can now be created in seconds, making these schemes ideally suited for the next generation of high performance 256-bit (double quadruple) or higher precision computers. This ease of creation makes it possible to adapt the algorithm to the mesh in time instead of its converse: this is ideal for resolving varying wavelength scales which occur in noise generation simulations. And finally, the sources of round-off error which effect the very high order methods are examined and remedies provided that effectively increase the accuracy of the MESA schemes while using current computer technology.
Comparison of rhodomine-WT and sodium chloride tracer transport in a 4th order arctic river
NASA Astrophysics Data System (ADS)
Smull, E. M.; Wlostowski, A. N.; Gooseff, M. N.; Bowden, W. B.; Wollheim, W. M.
2012-12-01
Conservative tracers are useful for tracking a parcel of water through a river reach and understanding tracer transport phenomena (i.e. advection, dispersion, and transient storage). Rhodomine- WT (RWT) and sodium chloride (NaCl) are two popular stream tracers. NaCl is considered to be conservative and relatively inexpensive, yet it cannot be detected at very low concentrations. On the other hand, RWT can be detected at very low concentrations (<0.1 ppb), but it is known to photo-degrade and sorb to organic materials. Previous work has compared these tracers with small-scale laboratory analyses and field experiments on small headwater streams. The limitations and advantages to each of these tracers, as applied to large river slug injections, are not clearly understood. This work seeks to answer the following questions: 1) Does RWT improve the tracer window of detection (time of tracer arrival to time of tracer non-detection), compared to NaCl? 2) Are there differences in the late-time tailing behavior of each tracer? More specifically, can we compare RWT and NaCl breakthrough curve tail shapes to understand processes contributing to late time solute transport (transient storage or sorption-desorption)? During the summer of 2012, combined slug additions of RWT and NaCl were injected into a 1.5-kilometer reach on the Kuparuk River, a 4th order tundra river underlain by continuous permafrost located on Alaska's North Slope. Fluorescence and electrical conductivity were continuously logged at the upstream and downstream ends of the reach. Preliminary results show that the window of detection is expanded when using RWT under both high and low flow conditions by 0.2 times the advective transport timescale. Tail shapes are more similar under higher discharge conditions and dissimilar under lower discharge conditions. For example, using an exponential regression model (c(t) = eat) to quantify tail shapes, at Q = 500 l/s the exponential coefficient ratio, aRWT:aNaCl, is 0
Higher-Order Compact Schemes for Numerical Simulation of Incompressible Flows
NASA Technical Reports Server (NTRS)
Wilson, Robert V.; Demuren, Ayodeji O.; Carpenter, Mark
1998-01-01
A higher order accurate numerical procedure has been developed for solving incompressible Navier-Stokes equations for 2D or 3D fluid flow problems. It is based on low-storage Runge-Kutta schemes for temporal discretization and fourth and sixth order compact finite-difference schemes for spatial discretization. The particular difficulty of satisfying the divergence-free velocity field required in incompressible fluid flow is resolved by solving a Poisson equation for pressure. It is demonstrated that for consistent global accuracy, it is necessary to employ the same order of accuracy in the discretization of the Poisson equation. Special care is also required to achieve the formal temporal accuracy of the Runge-Kutta schemes. The accuracy of the present procedure is demonstrated by application to several pertinent benchmark problems.
A third-order multistep time discretization for a Chebyshev tau spectral method
NASA Astrophysics Data System (ADS)
Vreman, A. W.; Kuerten, J. G. M.
2016-01-01
A time discretization scheme based on the third-order backward difference formula has been embedded into a Chebyshev tau spectral method for the Navier-Stokes equations. The time discretization is a variant of the second-order backward scheme proposed by Krasnov et al. (2008) [3]. High-resolution direct numerical simulations of turbulent incompressible channel flow have been performed to compare the backward scheme to the Runge-Kutta scheme proposed by Spalart et al. (1991) [2]. It is shown that the Runge-Kutta scheme leads to a poor convergence of some third-order spatial derivatives in the direct vicinity of the wall, derivatives that represent the diffusion of wall-tangential vorticity. The convergence at the wall is shown to be significantly improved if the backward scheme is applied.
On the Total Variation of High-Order Semi-Discrete Central Schemes for Conservation Laws
NASA Technical Reports Server (NTRS)
Bryson, Steve; Levy, Doron
2004-01-01
We discuss a new fifth-order, semi-discrete, central-upwind scheme for solving one-dimensional systems of conservation laws. This scheme combines a fifth-order WENO reconstruction, a semi-discrete central-upwind numerical flux, and a strong stability preserving Runge-Kutta method. We test our method with various examples, and give particular attention to the evolution of the total variation of the approximations.
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.
DNS and LES of Turbulent Backward-Facing Step Flow Using 2ND-and 4TH-Order Discretization
NASA Astrophysics Data System (ADS)
Meri, Adnan; Wengle, Hans
Results are presented from a Direct Numerical Simulation (DNS) and Large-Eddy Simulations (LES) of turbulent flow over a backward-facing step (Reh=3300) with a fully developed channel flow (Rcτ=180) utilized asatime-dependent inflow condition. Numerical solutions using a fourth-order compact (Hermitian) scheme, which was formulated directly for anon-equidistant and staggered grid in [1] are compared with numerical solutions using the classical second-order central scheme. There sults from LES (using the dynamic subgrid scale model) are evaluated against a corresponding DNS reference data set (fourth-order solution).
NASA Astrophysics Data System (ADS)
Zhao, Ye; Gu, Zhuquan; Liu, Yafeng
2012-07-01
In this paper, the Neumann system for the 4th-order eigenvalue problem Ly = (∂4+ q∂2+∂2 q+ ip∂+ i∂ p+ y = Λy) has been given. By means of the Neumann constraint condition, the perfect constraint set Γ and the relations between the potentials { q, p, r} and the eigenvector y are obtained. Then, based on the Euler-Lagrange function and Legendre transformations, a reasonable Jacobi-Ostrogradsky coordinate system has been found, which can be equal to the real Hamiltonian canonical coordinate system in R 8 N . Using Cao's method and Moser's constraint manifold, the Lax pairs of the evolution equation hierarchy with the 4th-order eigenvalue problems are nonlinearized. So a new finite-dimensional integrable Hamilton system on the constraint submanifold R 8 N-4 is generated. Moreover, the solutions of the evolution equations for the infinite-dimensional soliton systems are obtained by the involutive flow of the finite-dimensional completely integrable systems.
ERIC Educational Resources Information Center
Baydo-Reed, Katie
2010-01-01
Following the bombing of Pearl Harbor on Dec. 7, 1941, U.S. officials issued a series of proclamations that violated the civil and human rights of the vast majority of Japanese Americans in the United States--ostensibly to protect the nation from further Japanese aggression. The proclamations culminated in Executive Order 9066, which gave the…
Linear 3 and 5-step methods using Taylor series expansion for solving special 3rd order ODEs
NASA Astrophysics Data System (ADS)
Rajabi, Marzieh; Ismail, Fudziah; Senu, Norazak
2016-06-01
Some new linear 3 and 5-step methods for solving special third order ordinary differential equations directly are constructed using Taylor's series expansion. A set of test problems are solved using the new method and the results are compared when the problem is reduced to a system of first order ordinary differential equations and then using the existing Runge-Kutta method. The numerical results have clearly shown the advantage and competency of the new methods.
NASA Technical Reports Server (NTRS)
Navon, I. M.; Bloom, S.; Takacs, L. L.
1985-01-01
An attempt was made to use the GLAS global 4th order shallow water equations to perform a Machenhauer nonlinear normal mode initialization (NLNMI) for the external vertical mode. A new algorithm was defined for identifying and filtering out computational modes which affect the convergence of the Machenhauer iterative procedure. The computational modes and zonal waves were linearly initialized and gravitational modes were nonlinearly initialized. The Machenhauer NLNMI was insensitive to the absence of high zonal wave numbers. The effects of the Machenhauer scheme were evaluated by performing 24 hr integrations with nondissipative and dissipative explicit time integration models. The NLNMI was found to be inferior to the Rasch (1984) pseudo-secant technique for obtaining convergence when the time scales of nonlinear forcing were much smaller than the time scales expected from the natural frequency of the mode.
NASA Astrophysics Data System (ADS)
O'Daniel, S. J.; Amerson, B. E.; Lambert, M. B.
2014-12-01
Persistent societal interest in improving water quality and recovering imperiled, native, aquatic species has expanded the scope of stream restoration to include the hyporheic zone as a focus. Despite the lack of detailed studies, hyporheic restoration is often invoked as a means to achieve multiple objectives including moderation of water temperature, delay of seasonal flows and increasing the localized volume of floodplain water. We present interim results from an ongoing case study that monitors the changes as a result of stream restoration of the hyporheic zone of a 4th order, alluvial floodplain in northeast Oregon, USA, Meacham Creek. Active and passive restoration of 2.5 km of Meacham Creek has altered the creek from a single-threaded, incised and bedrock-dominated channel to a perched, alluvial channel that seasonally exchanges overbank flows with the surrounding floodplain. Our results suggest that the stream restoration effort on Meacham Creek has increased the volume of annual hyporheic storage and created a more diverse distribution of flowpath lengths within the restoration site. Furthermore, our monitoring indicates that hyporheic process response to stream restoration, analogous to other geomorphic processes, conforms to a systematic hierarchy where nested flow paths range in length and residence time from meters and hours at the habitat scale to tens of meters and months at the floodplain scale. We assert that scale-explicit and measurement-focused restoration planning has a greater likelihood of meeting the stated objectives and result in improved water quality and encourage recovery of many native aquatic species.
NASA Astrophysics Data System (ADS)
Mongeon, Michael C.
1996-03-01
This paper investigates the development of printer device profiles used in color document printing system environments when devices with intrinsically different gamut capabilities communicate with one another in a common (CIELAB) color space. While the main thrust of this activity focuses on the output printer, namely the Xerox 5760 printer, and its rendition of some device independent image description, characterizations are provided which investigate relative areas of photographic, monitor, and printer gamuts using a visual hue leaf comparison between devices. The printer is modeled using 4th-order polynomial regression which maps the device independent CIELAB image representation into device dependent printer CMYK. This technique results in 1.89 AEEavg over the training data set. Some key properties of the proposed calibration method are as follows: (1) Linearized CMYK tone reproduction curves with respect to AEEpaper to improve the distribution of calibration data in color space. (2) Application of GCR strategy and linearization to the calibration target prior to the regression on the measured CIELAB and original CMY values. Each strategy employs a K addition/No CMY removal method which maximizes printer gamut and relies on the regression to determine the appropriate CMY removal. The following GCR strategies are explored: CMY only (0% K addition), 50% K addition, 100% K addition, and non-linear K addition. A library of image processing algorithms is included, using LabView object oriented programming, which provides a modular approach for key color processing tasks. In the user interface, an image is selected with appropriate GCR strategy, and the program operates on the image. In general, the pictorial image quality is excellent for each GCR strategy with subtle differences between GCR approaches. Quantitative analysis of Q60 color matching performance is included.
NASA Astrophysics Data System (ADS)
Scarponi, D.; Kaufman, D.; Bright, J.; Kowalewski, M.
2009-04-01
Single fossiliferous beds contain biotic remnants that commonly vary in age over a time span of hundreds to thousands of years. Multiple recent studies suggest that such temporal mixing is a widespread phenomenon in marine depositional systems. This research focuses on quantitative estimates of temporal mixing obtained by direct dating of individual corbulid bivalve shells (Lentidium mediterraneum and Corbula gibba) from Po plain marine units of the Holocene 4th-order depositional sequence, including Transgressive Systems Tract [TST] and Highstand Systems Tract [HST]. These units displays a distinctive succession of facies consisting of brackish to marginal marine retrogradational deposits, (early TST), overlain by fully marine fine to coarse gray sands (late TST), and capped with progradational deltaic clays and sands (HST). More than 300 corbulid specimens, representing 19 shell-rich horizons evenly distributed along the depositional sequence and sampled from 9 cores, have been dated by means of aspartic acid racemization calibrated using 23 AMS-radiocarbon dates (14 dates for Lentidium mediterraneum and 9 dates for Corbula gibba, respectively). The results indicate that the scale of time-averaging is comparable when similar depositional environments from the same systems tract are compared across cores. However, time averaging is notably different when similar depositional environments from TST and HST segments of the sequence are compared. Specifically, late HST horizons (n=8) display relatively low levels of time-averaging: the mean within-horizon range of shell ages is 537 years and standard deviation averages 165 years. In contrast, late TST horizons (n=7) are dramatically more time-averaged: mean range of 5104 years and mean standard deviations of 1420 years. Thus, late TST horizons experience a 1 order of magnitude higher time-averaging than environmentally comparable late HST horizons. In conclusion the HST and TST systems tracts of the Po Plain display
Contractivity-preserving explicit Hermite-Obrechkoff ODE solver of order 13
NASA Astrophysics Data System (ADS)
Nguyen-Ba, Truong; Desjardins, Steven J.; Sharp, Philip W.; Vaillancourt, Rémi
2013-12-01
A new optimal, explicit, Hermite-Obrechkoff method of order 13, denoted by HO(13), that is contractivity-preserving (CP) and has nonnegative coefficients is constructed for solving nonstiff first-order initial value problems. Based on the CP conditions, the new 9-derivative HO(13) has maximum order 13. The new method usually requires significantly fewer function evaluations and significantly less CPU time than the Taylor method of order 13 and the Runge-Kutta method DP(8,7)13M to achieve the same global error when solving standard -body problems.
Liang, Xiao; Khaliq, Abdul Q. M.; Xing, Yulong
2015-01-23
In this paper, we study a local discontinuous Galerkin method combined with fourth order exponential time differencing Runge-Kutta time discretization and a fourth order conservative method for solving the nonlinear Schrödinger equations. Based on different choices of numerical fluxes, we propose both energy-conserving and energy-dissipative local discontinuous Galerkin methods, and have proven the error estimates for the semi-discrete methods applied to linear Schrödinger equation. The numerical methods are proven to be highly efficient and stable for long-range soliton computations. Finally, extensive numerical examples are provided to illustrate the accuracy, efficiency and reliability of the proposed methods.
A second-order accurate immersed boundary-lattice Boltzmann method for particle-laden flows
NASA Astrophysics Data System (ADS)
Zhou, Qiang; Fan, Liang-Shih
2014-07-01
A new immersed boundary-lattice Boltzmann method (IB-LBM) is presented for fully resolved simulations of incompressible viscous flows laden with rigid particles. The immersed boundary method (IBM) recently developed by Breugem (2012) [19] is adopted in the present method, development including the retraction technique, the multi-direct forcing method and the direct account of the inertia of the fluid contained within the particles. The present IB-LBM is, however, formulated with further improvement with the implementation of the high-order Runge-Kutta schemes in the coupled fluid-particle interaction. The major challenge to implement high-order Runge-Kutta schemes in the LBM is that the flow information such as density and velocity cannot be directly obtained at a fractional time step from the LBM since the LBM only provides the flow information at an integer time step. This challenge can be, however, overcome as given in the present IB-LBM by extrapolating the flow field around particles from the known flow field at the previous integer time step. The newly calculated fluid-particle interactions from the previous fractional time steps of the current integer time step are also accounted for in the extrapolation. The IB-LBM with high-order Runge-Kutta schemes developed in this study is validated by several benchmark applications. It is demonstrated, for the first time, that the IB-LBM has the capacity to resolve the translational and rotational motion of particles with the second-order accuracy. The optimal retraction distances for spheres and tubes that help the method achieve the second-order accuracy are found to be around 0.30 and -0.47 times of the lattice spacing, respectively. Simulations of the Stokes flow through a simple cubic lattice of rotational spheres indicate that the lift force produced by the Magnus effect can be very significant in view of the magnitude of the drag force when the practical rotating speed of the spheres is encountered. This finding
McCorquodale, Peter; Ullrich, Paul; Johansen, Hans; Colella, Phillip
2015-09-04
We present a high-order finite-volume approach for solving the shallow-water equations on the sphere, using multiblock grids on the cubed-sphere. This approach combines a Runge--Kutta time discretization with a fourth-order accurate spatial discretization, and includes adaptive mesh refinement and refinement in time. Results of tests show fourth-order convergence for the shallow-water equations as well as for advection in a highly deformational flow. Hierarchical adaptive mesh refinement allows solution error to be achieved that is comparable to that obtained with uniform resolution of the most refined level of the hierarchy, but with many fewer operations.
Liang, Xiao; Khaliq, Abdul Q.M.; Xing, Yulong
2015-01-23
In this paper, we study a local discontinuous Galerkin method combined with fourth order exponential time differencing Runge-Kutta time discretization and a fourth order conservative method for solving the nonlinear Schrödinger equations. Based on different choices of numerical fluxes, we propose both energy-conserving and energy-dissipative local discontinuous Galerkin methods, and have proven the error estimates for the semi-discrete methods applied to linear Schrödinger equation. The numerical methods are proven to be highly efficient and stable for long-range soliton computations. Finally, extensive numerical examples are provided to illustrate the accuracy, efficiency and reliability of the proposed methods.
Guzik, S; McCorquodale, P; Colella, P
2011-12-16
A fourth-order accurate finite-volume method is presented for solving time-dependent hyperbolic systems of conservation laws on mapped grids that are adaptively refined in space and time. Novel considerations for formulating the semi-discrete system of equations in computational space combined with detailed mechanisms for accommodating the adapting grids ensure that conservation is maintained and that the divergence of a constant vector field is always zero (freestream-preservation property). Advancement in time is achieved with a fourth-order Runge-Kutta method.
NASA Astrophysics Data System (ADS)
Olemskoy, I. V.; Eremin, A. S.
2016-06-01
We construct here an embedded Dormand-Prince pair of explicit methods of orders 6 and 4 for systems of ordinary differential equations with special structure, namely with two parts, in which the right-hand sides are dependent only on the unknown functions from the other group. The number of stages is six, which is fewer than for general explicit Runge-Kutta methods. The comparison to Dormand-Prince method of the same computation cost is made showing the higher accuracy of the suggested method.
The 4th Thermodynamic Principle?
Montero Garcia, Jose de la Luz; Novoa Blanco, Jesus Francisco
2007-04-28
It should be emphasized that the 4th Principle above formulated is a thermodynamic principle and, at the same time, is mechanical-quantum and relativist, as it should inevitably be and its absence has been one of main the theoretical limitations of the physical theory until today.We show that the theoretical discovery of Dimensional Primitive Octet of Matter, the 4th Thermodynamic Principle, the Quantum Hexet of Matter, the Global Hexagonal Subsystem of Fundamental Constants of Energy and the Measurement or Connected Global Scale or Universal Existential Interval of the Matter is that it is possible to be arrived at a global formulation of the four 'forces' or fundamental interactions of nature. The Einstein's golden dream is possible.
Seismic Waves, 4th order accurate
Energy Science and Technology Software Center (ESTSC)
2013-08-16
SW4 is a program for simulating seismic wave propagation on parallel computers. SW4 colves the seismic wave equations in Cartesian corrdinates. It is therefore appropriate for regional simulations, where the curvature of the earth can be neglected. SW4 implements a free surface boundary condition on a realistic topography, absorbing super-grid conditions on the far-field boundaries, and a kinematic source model consisting of point force and/or point moment tensor source terms. SW4 supports a fully 3-Dmore » heterogeneous material model that can be specified in several formats. SW4 can output synthetic seismograms in an ASCII test format, or in the SAC finary format. It can also present simulation information as GMT scripts, whixh can be used to create annotated maps. Furthermore, SW4 can output the solution as well as the material model along 2-D grid planes.« less
Seismic Waves, 4th order accurate
2013-08-16
SW4 is a program for simulating seismic wave propagation on parallel computers. SW4 colves the seismic wave equations in Cartesian corrdinates. It is therefore appropriate for regional simulations, where the curvature of the earth can be neglected. SW4 implements a free surface boundary condition on a realistic topography, absorbing super-grid conditions on the far-field boundaries, and a kinematic source model consisting of point force and/or point moment tensor source terms. SW4 supports a fully 3-D heterogeneous material model that can be specified in several formats. SW4 can output synthetic seismograms in an ASCII test format, or in the SAC finary format. It can also present simulation information as GMT scripts, whixh can be used to create annotated maps. Furthermore, SW4 can output the solution as well as the material model along 2-D grid planes.
A high-order finite-volume method for hyperbolic conservation laws on locally-refined grids
McCorquodale, Peter; Colella, Phillip
2011-01-28
We present a fourth-order accurate finite-volume method for solving time-dependent hyperbolic systems of conservation laws on Cartesian grids with multiple levels of refinement. The underlying method is a generalization of that in [5] to nonlinear systems, and is based on using fourth-order accurate quadratures for computing fluxes on faces, combined with fourth-order accurate Runge?Kutta discretization in time. To interpolate boundary conditions at refinement boundaries, we interpolate in time in a manner consistent with the individual stages of the Runge-Kutta method, and interpolate in space by solving a least-squares problem over a neighborhood of each target cell for the coefficients of a cubic polynomial. The method also uses a variation on the extremum-preserving limiter in [8], as well as slope flattening and a fourth-order accurate artificial viscosity for strong shocks. We show that the resulting method is fourth-order accurate for smooth solutions, and is robust in the presence of complex combinations of shocks and smooth flows.
Sanders, Ross H; Gonjo, Tomohiro; McCabe, Carla B
2015-06-01
The purpose of this study was to explore the reliability of estimating three-dimensional (3D) linear kinematics and kinetics of a swimmer derived from digitized video and to assess the effect of framing rate and smoothing window size. A stroke cycle of two high-level front crawl swimmers and one high level backstroke swimmer was recorded by four underwater and two above water video cameras. One of the front crawl swimmers was recorded and digitized at 50 Hz with a window for smoothing by 4(th) order Butterworth digital filter extending 10 frames beyond the start and finish of the stroke cycle, while the other front crawl and backstroke swimmer were recorded and digitized at 25 Hz with the window extending five frames beyond the start and finish of the stroke cycle. Each camera view of the stroke cycle was digitized five times yielding five independent 3D data sets from which whole body centre of mass (CM) component velocities and accelerations were derived together with wrist and ankle linear velocities. Coefficients of reliability ranging from r = 0.942 to r = 0.999 indicated that both methods are sufficiently reliable to identify real differences in net force production during the pulls of the right and left hands. Reliability of digitizing was better for front crawl when digitizing at 50Hz with 10 frames extension than at 25 Hz with 5 frames extension (p < 0.01) and better for backstroke than front crawl (p < 0.01). However, despite the extension and reflection of data, errors were larger in the first 15% of the stroke cycle than the period between 15 and 85% of the stroke cycle for CM velocity and acceleration and for foot speed (p < 0.01). Key pointsAn inverse dynamics based on 3D position data digitized from multiple camera views above and below the water surface is sufficiently reliable to yield insights regarding force production in swimming additional to those of other approaches.The ability to link the force profiles to swimming actions and technique is
NASA Astrophysics Data System (ADS)
Ying, Teh Yuan; Yaacob, Nazeeruddin
2013-04-01
In this paper, a new implicit Runge-Kutta method which based on a 4-point Gauss-Kronrod-Radau II quadrature formula is developed. The resulting implicit method is a 4-stage sixth order Gauss-Kronrod-Radau IIA method, or in brief as GKRM(4,6)-IIA. GKRM(4,6)-IIA requires four function of evaluations at each integration step and it gives accuracy of order six. In addition, GKRM(4,6)-IIA has stage order four and being L-stable. Numerical experiments compare the accuracy between GKRM(4,6)-IIA and the classical 3-stage sixth order Gauss-Legendre method in solving some test problems. Numerical results reveal that GKRM(4,6)-IIA is more accurate than the 3-stage sixth order Gauss-Legendre method because GKRM(4,6)-IIA has higher stage order.
NASA Astrophysics Data System (ADS)
Ying, Teh Yuan; Yaacob, Nazeeruddin
2013-04-01
In this paper, a new implicit Runge-Kutta method which based on a 7-point Gauss-Kronrod-Lobatto quadrature formula is developed. The resulting implicit method is a 7-stage tenth order Gauss-Kronrod-Lobatto IIIA method, or in brief as GKLM(7,10)-IIIA. GKLM(7,10)-IIIA requires seven function of evaluations at each integration step and it gives accuracy of order ten. In addition, GKLM(7,10)-IIIA has stage order seven and being A-stable. Numerical experiments compare the accuracy between GKLM(7,10)-IIIA and the classical 5-stage tenth order Gauss-Legendre method in solving some test problems. Numerical results reveal that GKLM(7,10)-IIIA is more accurate than the 5-stage tenth order Gauss-Legendre method because GKLM(7,10)-IIIA has higher stage order.
17. 4th floor roof, view south, 4th and 5th floor ...
17. 4th floor roof, view south, 4th and 5th floor setback to left and atrium structure to right - Sheffield Farms Milk Plant, 1075 Webster Avenue (southwest corner of 166th Street), Bronx, Bronx County, NY
Finite-volume application of high-order ENO schemes to two-dimensional boundary-value problems
NASA Technical Reports Server (NTRS)
Casper, Jay
1991-01-01
Finite-volume applications of high-order accurate ENO schemes to two-dimensional boundary-value problems are studied. These schemes achieve high-order spatial accuracy, in smooth regions, by a piecewise polynomial approximation of the solution from cell averages. In addition, this spatial operation involves an adaptive stencil algorithm in order to avoid the oscillatory behavior that is associated with interpolation across steep gradients. High-order TVD Runge-Kutta methods are employed for time integration, thus making these schemes best suited for unsteady problems. Fifth- and sixth-order accurate applications are validated through a grid refinement study involving the solutions of scalar hyperbolic equations. A previously proposed extension for the Euler equations of gas dynamics is tested, including its application to solutions of boundary-value problems involving solid walls and curvilinear coordinates.
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.
Direct calculations of waves in fluid flows using a high-order compact difference scheme
NASA Technical Reports Server (NTRS)
Yu, Sheng-Tao; Hultgren, Lennart S.; Liu, Nan-Suey
1993-01-01
The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expressions for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations, and the numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the CFL limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about 6 to 8 grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 to 300 time steps. Three numerical examples of waves in compressible flow are included.
NASA Technical Reports Server (NTRS)
Fatemi, Emad; Jerome, Joseph; Osher, Stanley
1989-01-01
A micron n+ - n - n+ silicon diode is simulated via the hydrodynamic model for carrier transport. The numerical algorithms employed are for the non-steady case, and a limiting process is used to reach steady state. The simulation employs shock capturing algorithms, and indeed shocks, or very rapid transition regimes, are observed in the transient case for the coupled system, consisting of the potential equation and the conservation equations describing charge, momentum, and energy transfer for the electron carriers. These algorithms, termed essentially non-oscillatory, were successfully applied in other contexts to model the flow in gas dynamics, magnetohydrodynamics, and other physical situations involving the conservation laws in fluid mechanics. The method here is first order in time, but the use of small time steps allows for good accuracy. Runge-Kutta methods allow one to achieve higher accuracy in time if desired. The spatial accuracy is of high order in regions of smoothness.
Cobb, J.W.
1995-02-01
There is an increasing need for more accurate numerical methods for large-scale nonlinear magneto-fluid turbulence calculations. These methods should not only increase the current state of the art in terms of accuracy, but should also continue to optimize other desired properties such as simplicity, minimized computation, minimized memory requirements, and robust stability. This includes the ability to stably solve stiff problems with long time-steps. This work discusses a general methodology for deriving higher-order numerical methods. It also discusses how the selection of various choices can affect the desired properties. The explicit discussion focuses on third-order Runge-Kutta methods, including general solutions and five examples. The study investigates the linear numerical analysis of these methods, including their accuracy, general stability, and stiff stability. Additional appendices discuss linear multistep methods, discuss directions for further work, and exhibit numerical analysis results for some other commonly used lower-order methods.
NASA Astrophysics Data System (ADS)
Campoamor-Stursberg, R.; Rodríguez, M. A.; Winternitz, P.
2016-01-01
Ordinary differential equations (ODEs) and ordinary difference systems (OΔSs) invariant under the actions of the Lie groups {{SL}}x(2),{{SL}}y(2) and {{SL}}x(2)× {{SL}}y(2) of projective transformations of the independent variables x and dependent variables y are constructed. The ODEs are continuous limits of the OΔSs, or conversely, the OΔSs are invariant discretizations of the ODEs. The invariant OΔSs are used to calculate numerical solutions of the invariant ODEs of order up to five. The solutions of the invariant numerical schemes are compared to numerical solutions obtained by standard Runge-Kutta methods and to exact solutions, when available. The invariant method performs at least as well as standard ones and much better in the vicinity of singularities of solutions.
NASA Astrophysics Data System (ADS)
Balan, Aravind; May, Georg; Schöberl, Joachim
2012-03-01
Numerical schemes using piecewise polynomial approximation are very popular for high order discretization of conservation laws. While the most widely used numerical scheme under this paradigm appears to be the Discontinuous Galerkin method, the Spectral Difference scheme has often been found attractive as well, because of its simplicity of formulation and implementation. However, recently it has been shown that the scheme is not linearly stable on triangles. In this paper we present an alternate formulation of the scheme, featuring a new flux interpolation technique using Raviart-Thomas spaces, which proves stable under a similar linear analysis in which the standard scheme failed. We demonstrate viability of the concept by showing linear stability both in the semi-discrete sense and for time stepping schemes of the SSP Runge-Kutta type. Furthermore, we present convergence studies, as well as case studies in compressible flow simulation using the Euler equations.
NASA Astrophysics Data System (ADS)
JavanNezhad, R.; Meshkatee, A. H.; Ghader, S.; Ahmadi-Givi, F.
2016-09-01
This study is devoted to application of the fourth-order compact MacCormack scheme to spatial differencing of the conservative form of two-dimensional and non-hydrostatic equation of a dry atmosphere. To advance the solution in time a four-stage Runge-Kutta method is used. To perform the simulations, two test cases including evolution of a warm bubble and a cold bubble in a neutral atmosphere with open and rigid boundaries are employed. In addition, the second-order MacCormack and the standard fourth-order compact MacCormack schemes are used to perform the simulations. Qualitative and quantitative assessment of the numerical results for different test cases exhibit the superiority of the fourth-order compact MacCormack scheme on the second-order method.
NASA Astrophysics Data System (ADS)
Li, G. Q.; Zhu, Z. H.
2015-12-01
Dynamic modeling of tethered spacecraft with the consideration of elasticity of tether is prone to the numerical instability and error accumulation over long-term numerical integration. This paper addresses the challenges by proposing a globally stable numerical approach with the nodal position finite element method (NPFEM) and the implicit, symplectic, 2-stage and 4th order Gaussian-Legendre Runge-Kutta time integration. The NPFEM eliminates the numerical error accumulation by using the position instead of displacement of tether as the state variable, while the symplectic integration enforces the energy and momentum conservation of the discretized finite element model to ensure the global stability of numerical solution. The effectiveness and robustness of the proposed approach is assessed by an elastic pendulum problem, whose dynamic response resembles that of tethered spacecraft, in comparison with the commonly used time integrators such as the classical 4th order Runge-Kutta schemes and other families of non-symplectic Runge-Kutta schemes. Numerical results show that the proposed approach is accurate and the energy of the corresponding numerical model is conservative over the long-term numerical integration. Finally, the proposed approach is applied to the dynamic modeling of deorbiting process of tethered spacecraft over a long period.
NASA Technical Reports Server (NTRS)
Shu, Chi-Wang
2004-01-01
This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier
Fourth order difference methods for hyperbolic IBVP's
NASA Technical Reports Server (NTRS)
Gustafsson, Bertil; Olsson, Pelle
1994-01-01
Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.
IMEX Runge-Kutta Schemes and Hyperbolic Systems of Conservation Laws with Stiff Diffusive Relaxation
NASA Astrophysics Data System (ADS)
Boscarino, S.; Pareschi, L.; Russo, G.
2009-09-01
Hyperbolic system of conservation laws often have relaxation terms that, under a suitable scaling, lead to a reduced system of parabolic or hyperbolic type. The development of numerical methods to solve systems of this form his an active area of research. These systems in addition to the stiff relaxation term have the convection term stiff too. In this paper we will mainly concentrate on the study of the stiff regime. In fact in this stiff regime most of the popular methods for the solution of these system fail to capture the correct behavior of the relaxation limit unless the small relaxation rate is numericaly resolved. We will show how to overcome this difficulties and how to construct numerical schemes with the correct asymnptotic limit, i.e., the correct zero-relaxation limit should be preserved at a discrete level.
A high-order gas-kinetic Navier-Stokes flow solver
Li Qibing; Xu Kun; Fu Song
2010-09-20
The foundation for the development of modern compressible flow solver is based on the Riemann solution of the inviscid Euler equations. The high-order schemes are basically related to high-order spatial interpolation or reconstruction. In order to overcome the low-order wave interaction mechanism due to the Riemann solution, the temporal accuracy of the scheme can be improved through the Runge-Kutta method, where the dynamic deficiencies in the first-order Riemann solution is alleviated through the sub-step spatial reconstruction in the Runge-Kutta process. The close coupling between the spatial and temporal evolution in the original nonlinear governing equations seems weakened due to its spatial and temporal decoupling. Many recently developed high-order methods require a Navier-Stokes flux function under piece-wise discontinuous high-order initial reconstruction. However, the piece-wise discontinuous initial data and the hyperbolic-parabolic nature of the Navier-Stokes equations seem inconsistent mathematically, such as the divergence of the viscous and heat conducting terms due to initial discontinuity. In this paper, based on the Boltzmann equation, we are going to present a time-dependent flux function from a high-order discontinuous reconstruction. The theoretical basis for such an approach is due to the fact that the Boltzmann equation has no specific requirement on the smoothness of the initial data and the kinetic equation has the mechanism to construct a dissipative wave structure starting from an initially discontinuous flow condition on a time scale being larger than the particle collision time. The current high-order flux evaluation method is an extension of the second-order gas-kinetic BGK scheme for the Navier-Stokes equations (BGK-NS). The novelty for the easy extension from a second-order to a higher order is due to the simple particle transport and collision mechanism on the microscopic level. This paper will present a hierarchy to construct such
NASA Astrophysics Data System (ADS)
Balac, Stéphane; Fernandez, Arnaud
2016-02-01
The computer program SPIP is aimed at solving the Generalized Non-Linear Schrödinger equation (GNLSE), involved in optics e.g. in the modelling of light-wave propagation in an optical fibre, by the Interaction Picture method, a new efficient alternative method to the Symmetric Split-Step method. In the SPIP program a dedicated costless adaptive step-size control based on the use of a 4th order embedded Runge-Kutta method is implemented in order to speed up the resolution.
A neuro approach to solve fuzzy Riccati differential equations
NASA Astrophysics Data System (ADS)
Shahrir, Mohammad Shazri; Kumaresan, N.; Kamali, M. Z. M.; Ratnavelu, Kurunathan
2015-10-01
There are many applications of optimal control theory especially in the area of control systems in engineering. In this paper, fuzzy quadratic Riccati differential equation is estimated using neural networks (NN). Previous works have shown reliable results using Runge-Kutta 4th order (RK4). The solution can be achieved by solving the 1st Order Non-linear Differential Equation (ODE) that is found commonly in Riccati differential equation. Research has shown improved results relatively to the RK4 method. It can be said that NN approach shows promising results with the advantage of continuous estimation and improved accuracy that can be produced over RK4.
A neuro approach to solve fuzzy Riccati differential equations
Shahrir, Mohammad Shazri; Kumaresan, N. Kamali, M. Z. M.; Ratnavelu, Kurunathan
2015-10-22
There are many applications of optimal control theory especially in the area of control systems in engineering. In this paper, fuzzy quadratic Riccati differential equation is estimated using neural networks (NN). Previous works have shown reliable results using Runge-Kutta 4th order (RK4). The solution can be achieved by solving the 1st Order Non-linear Differential Equation (ODE) that is found commonly in Riccati differential equation. Research has shown improved results relatively to the RK4 method. It can be said that NN approach shows promising results with the advantage of continuous estimation and improved accuracy that can be produced over RK4.
Compact high order schemes for the Euler equations
NASA Technical Reports Server (NTRS)
Abarbanel, Saul; Kumar, Ajay
1988-01-01
An implicit approximate factorization (AF) algorithm is constructed which has the following characteistics. In 2-D: The scheme is unconditionally stable, has a 3 x 3 stencil and at steady state has a fourth order spatial accuracy. The temporal evolution is time accurate either to first or second order through choice of parameter. In 3-D: The scheme has almost the same properties as in 2-D except that it is now only conditionally stable, with the stability condition (the CFL number) being dependent on the cell aspect ratios, delta y/delta x and delta z/delta x. The stencil is still compact and fourth order accuracy at steady state is maintained. Numerical experiments on a 2-D shock-reflection problem show the expected improvement over lower order schemes, not only in accuracy (measured by the L sub 2 error) but also in the dispersion. It is also shown how the same technique is immediately extendable to Runge-Kutta type schemes resulting in improved stability in addition to the enhanced accuracy.
Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows
NASA Technical Reports Server (NTRS)
Sjoegreen, Bjoern; Yee, H. C.; Tang, Harry (Technical Monitor)
2002-01-01
Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think
NASA Astrophysics Data System (ADS)
Li Chun Fong, Lena C. M.; Lach, Grzegorz; Le Roy, Robert J.; Dattani, Nikesh S.
2015-06-01
The 13.81(8)s half-life of the halo nucleonic atom 11Be is orders of magnitude longer than those for any other halo nucleonic atom known, and makes Be-based diatomics the most promising candidates for the formation of the first halo nucleonic molecules. However, the 4e^- species LiH and BeH^+ are some of the first molecules for which the highest accuracy ab initio methods are not accessible, so empirical potential energy functions will be important for making predictions and for benchmarking how ab initio calculations break down at this transition from 3e^- to 4e^-. BeH^+ is also very light, and has one of the most extensive data sets involving a tritium isotopologue, making it a very useful benchmark for studying Born-Oppenheimer breakdown. We therefore seek to determine an empirical analytic potential energy function for BeH^+ that has as much precision as possible. To this end, all available spectroscopic data for all stable isotopologues of BeH^+ are analyzed in a standard direct-potential-fit procedure that uses least-squares fits to optimize the parameters defining an analytic potential. The ``Morse/Long-range'' (MLR) model used for the potential energy function incorporates the inverse-power long-range tail required by theory, and the calculation of the leading long-range coefficients C_4, C_6, C_7, and C_8 include non-adiabatic terms, and up to 4th order QED corrections. As a by-product, we have calculated some fundamental properties of 1e^- systems with unprecedented precision, such as the dipole, quadrupole, octupole, non-adiabatic, and mixed higher order polarizabilities of hydrogen, deuterium, and tritium. We provide good first estimates for the transition energies for the halo nucleonic species 11BeH^+ and 14BeH^+.
Mathematical modeling of intrinsic Josephson junctions with capacitive and inductive couplings
NASA Astrophysics Data System (ADS)
Rahmonov, I. R.; Shukrinov, Yu M.; Zemlyanaya, E. V.; Sarhadov, I.; Andreeva, O.
2012-11-01
We investigate the current voltage characteristics (CVC) of intrinsic Josephson junctions (IJJ) with two types of couplings between junctions: capacitive and inductive. The IJJ model is described by a system of coupled sine-Gordon equations which is solved numerically by the 4th order Runge-Kutta method. The method of numerical simulation and numerical results are presented. The magnetic field distribution is calculated as the function of coordinate and time at different values of the bias current. The influence of model parameters on the CVC is studied. The behavior of the IJJ in dependence on coupling parameters is discussed.
Improving Social Interaction among 4th Grade Students through Social Skills Instruction.
ERIC Educational Resources Information Center
Dunleavy, Shannon; Karwowski, Sandra; Shudes-Eitel, Jennifer
This action research project implemented a program for improving social skills in order to establish positive interaction among 4th grade students at a northern Chicago suburban school. Social skills deficiency was documented through behavior checklists and referrals, teacher observations and student reflection. Teachers reported that low incomes,…
An almost symmetric Strang splitting scheme for the construction of high order composition methods.
Einkemmer, Lukas; Ostermann, Alexander
2014-12-01
In this paper we consider splitting methods for nonlinear ordinary differential equations in which one of the (partial) flows that results from the splitting procedure cannot be computed exactly. Instead, we insert a well-chosen state [Formula: see text] into the corresponding nonlinearity [Formula: see text], which results in a linear term [Formula: see text] whose exact flow can be determined efficiently. Therefore, in the spirit of splitting methods, it is still possible for the numerical simulation to satisfy certain properties of the exact flow. However, Strang splitting is no longer symmetric (even though it is still a second order method) and thus high order composition methods are not easily attainable. We will show that an iterated Strang splitting scheme can be constructed which yields a method that is symmetric up to a given order. This method can then be used to attain high order composition schemes. We will illustrate our theoretical results, up to order six, by conducting numerical experiments for a charged particle in an inhomogeneous electric field, a post-Newtonian computation in celestial mechanics, and a nonlinear population model and show that the methods constructed yield superior efficiency as compared to Strang splitting. For the first example we also perform a comparison with the standard fourth order Runge-Kutta methods and find significant gains in efficiency as well better conservation properties. PMID:25473146
Toward a consistent framework for high order mesh refinement schemes in numerical relativity
NASA Astrophysics Data System (ADS)
Mongwane, Bishop
2015-05-01
It has now become customary in the field of numerical relativity to couple high order finite difference schemes to mesh refinement algorithms. To this end, different modifications to the standard Berger-Oliger adaptive mesh refinement algorithm have been proposed. In this work we present a fourth order stable mesh refinement scheme with sub-cycling in time for numerical relativity. We do not use buffer zones to deal with refinement boundaries but explicitly specify boundary data for refined grids. We argue that the incompatibility of the standard mesh refinement algorithm with higher order Runge Kutta methods is a manifestation of order reduction phenomena, caused by inconsistent application of boundary data in the refined grids. Our scheme also addresses the problem of spurious reflections that are generated when propagating waves cross mesh refinement boundaries. We introduce a transition zone on refined levels within which the phase velocity of propagating modes is allowed to decelerate in order to smoothly match the phase velocity of coarser grids. We apply the method to test problems involving propagating waves and show a significant reduction in spurious reflections.
Fourth-order compact schemes for the numerical simulation of coupled Burgers' equation
NASA Astrophysics Data System (ADS)
Bhatt, H. P.; Khaliq, A. Q. M.
2016-03-01
This paper introduces two new modified fourth-order exponential time differencing Runge-Kutta (ETDRK) schemes in combination with a global fourth-order compact finite difference scheme (in space) for direct integration of nonlinear coupled viscous Burgers' equations in their original form without using any transformations or linearization techniques. One scheme is a modification of the Cox and Matthews ETDRK4 scheme based on (1 , 3) -Padé approximation and other is a modification of Krogstad's ETDRK4-B scheme based on (2 , 2) -Padé approximation. Efficient versions of the proposed schemes are obtained by using a partial fraction splitting technique of rational functions. The stability properties of the proposed schemes are studied by plotting the stability regions, which provide an explanation of their behavior for dispersive and dissipative problems. The order of convergence of the schemes is examined empirically and found that the modification of ETDRK4 converges with the expected rate even if the initial data are nonsmooth. On the other hand, modification of ETDRK4-B suffers with order reduction if the initial data are nonsmooth. Several numerical experiments are carried out in order to demonstrate the performance and adaptability of the proposed schemes. The numerical results indicate that the proposed schemes provide better accuracy than other schemes available in the literature. Moreover, the results show that the modification of ETDRK4 is reliable and yields more accurate results than modification of ETDRK4-B, while solving problems with nonsmooth data or with high Reynolds number.
A Fourth Order Difference Scheme for the Maxwell Equations on Yee Grid
Fathy, Aly E; Wilson, Joshua L
2008-09-01
The Maxwell equations are solved by a long-stencil fourth order finite difference method over a Yee grid, in which different physical variables are located at staggered mesh points. A careful treatment of the numerical values near the boundary is introduced, which in turn leads to a 'symmetric image' formula at the 'ghost' grid points. Such a symmetric formula assures the stability of the boundary extrapolation. In turn, the fourth order discrete curl operator for the electric and magnetic vectors gives a complete set of eigenvalues in the purely imaginary axis. To advance the dynamic equations, the four-stage Runge-Kutta method is utilized, which results in a full fourth order accuracy in both time and space. A stability constraint for the time step is formulated at both the theoretical and numerical levels, using an argument of stability domain. An accuracy check is presented to verify the fourth order precision, using a comparison between exact solution and numerical solutions at a fixed final time. In addition, some numerical simulations of a loss-less rectangular cavity are also carried out and the frequency is measured precisely.
Effects of spatial order of accuracy on the computation of vortical flowfields
NASA Technical Reports Server (NTRS)
Ekaterinaris, J. A.
1993-01-01
The effect of the order-of-accuracy, used for the spatial discretization, on the resolution of the leading edge vortices over sharp-edged delta wings is investigated. The flowfield is computed using a viscous/inviscid zonal approach. The viscous flow in the vicinity of the wing is computed using the conservative formulation of the compressible, thin-layer Navier-Stokes equations. The leeward-side vortical flowfield and the other flow regions away from the surface are computed as inviscid. The time integration is performed with both an explicit fourth-order Runge-Kutta scheme and an implicit, factorized, iterative scheme. High-order-accurate inviscid fluxes are computed using both a conservative and a non-conservative (primitive variable) formulation. The nonlinear, inviscid terms of the primitive variable form of the governing equations are evaluated with a finite-difference numerical scheme based on the sign of the eigenvalues. High-order, upwind-biased, finite difference formulas are used to evaluate the derivatives of the nonlinear convective terms. Computed results are compared with available experimental data, and comparisons of the flowfield in the vicinity of the vortex cores are presented.
A second order residual based predictor-corrector approach for time dependent pollutant transport
NASA Astrophysics Data System (ADS)
Pavan, S.; Hervouet, J.-M.; Ricchiuto, M.; Ata, R.
2016-08-01
We present a second order residual distribution scheme for scalar transport problems in shallow water flows. The scheme, suitable for the unsteady cases, is obtained adapting to the shallow water context the explicit Runge-Kutta schemes for scalar equations [1]. The resulting scheme is decoupled from the hydrodynamics yet the continuity equation has to be considered in order to respect some important numerical properties at discrete level. Beyond the classical characteristics of the residual formulation presented in [1,2], we introduce the possibility to iterate the corrector step in order to improve the accuracy of the scheme. Another novelty is that the scheme is based on a precise monotonicity condition which guarantees the respect of the maximum principle. We thus end up with a scheme which is mass conservative, second order accurate and monotone. These properties are checked in the numerical tests, where the proposed approach is also compared to some finite volume schemes on unstructured grids. The results obtained show the interest in adopting the predictor-corrector scheme for pollutant transport applications, where conservation of the mass, monotonicity and accuracy are the most relevant concerns.
NASA Astrophysics Data System (ADS)
Núñez-de la Rosa, Jonatan; Munz, Claus-Dieter
2016-02-01
The present work describes the building blocks of a new code for computational magnetohydrodynamics based on very high order finite volume methods on Cartesian meshes. Spatial high-order accuracy is obtained with a weighted essentially non-oscillatory (WENO) reconstruction operator up to seventh order, while the time discretization is performed with a fourth-order strong-stability preserving Runge-Kutta method. Based on a shock-detection approach, the reconstruction operator employs a very high order WENO scheme in smooth flow regions and a third-order WENO scheme in those parts of the flow with discontinuities or shocks. The generalized Lagrange multiplier method is employed to enforce the solenoidal constraint on the magnetic field. Extensive numerical computations in one and two space dimensions are reported. Convergence rates for smooth flows verify the high-order accuracy of the scheme, and tests with strong shocks, including the Orszag-Tang vortex, the cylindrical blast wave problem, the rotor problem, and the Kelvin-Helmholtz instability, confirm the robustness and stability of the approach.
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.
A high-order element-based Galerkin Method for the global shallow water equations.
Nair, Ramachandran D.; Tufo, Henry M.; Levy, Michael Nathan
2010-08-01
The shallow water equations are used as a test for many atmospheric models because the solution mimics the horizontal aspects of atmospheric dynamics while the simplicity of the equations make them useful for numerical experiments. This study describes a high-order element-based Galerkin method for the global shallow water equations using absolute vorticity, divergence, and fluid depth (atmospheric thickness) as the prognostic variables, while the wind field is a diagnostic variable that can be calculated from the stream function and velocity potential (the Laplacians of which are the vorticity and divergence, respectively). The numerical method employed to solve the shallow water system is based on the discontinuous Galerkin and spectral element methods. The discontinuous Galerkin method, which is inherently conservative, is used to solve the equations governing two conservative variables - absolute vorticity and atmospheric thickness (mass). The spectral element method is used to solve the divergence equation and the Poisson equations for the velocity potential and the stream function. Time integration is done with an explicit strong stability-preserving second-order Runge-Kutta scheme and the wind field is updated directly from the vorticity and divergence at each stage, and the computational domain is the cubed sphere. A stable steady-state test is run and convergence results are provided, showing that the method is high-order accurate. Additionally, two tests without analytic solutions are run with comparable results to previous high-resolution runs found in the literature.
High-order central Hermite WENO schemes: Dimension-by-dimension moment-based reconstructions
NASA Astrophysics Data System (ADS)
Tao, Zhanjing; Li, Fengyan; Qiu, Jianxian
2016-08-01
In this paper, a class of high-order central finite volume schemes is proposed for solving one- and two-dimensional hyperbolic conservation laws. Formulated on staggered meshes, the methods involve Hermite WENO (HWENO) spatial reconstructions, and Lax-Wendroff type discretizations or the natural continuous extension of Runge-Kutta methods in time. Differently from the central Hermite WENO methods we developed previously in Tao et al. (2015) [34], the spatial reconstructions, a core ingredient of the methods, are based on the zeroth-order and the first-order moments of the solution, and are implemented through a dimension-by-dimension strategy when the spatial dimension is higher than one. This leads to much simpler implementation of the methods in higher dimension and better cost efficiency. Meanwhile, the proposed methods have the attractive features of the general central Hermite WENO methods such as being compact in reconstruction and requiring neither flux splitting nor numerical fluxes, while being accurate and essentially non-oscillatory. A collection of one- and two-dimensional numerical examples is presented to demonstrate high resolution and robustness of the methods in capturing smooth and non-smooth solutions.
Validation of a High-Order Prefactored Compact Scheme on Nonlinear Flows with Complex Geometries
NASA Technical Reports Server (NTRS)
Hixon, Ray; Mankbadi, Reda R.; Povinelli, L. A. (Technical Monitor)
2000-01-01
Three benchmark problems are solved using a sixth-order prefactored compact scheme employing an explicit 10th-order filter with optimized fourth-order Runge-Kutta time stepping. The problems solved are the following: (1) propagation of sound waves through a transonic nozzle; (2) shock-sound interaction; and (3) single airfoil gust response. In the first two problems, the spatial accuracy of the scheme is tested on a stretched grid, and the effectiveness of boundary conditions is shown. The solution stability and accuracy near a shock discontinuity is shown as well. Also, 1-D nonlinear characteristic boundary conditions will be evaluated. In the third problem, a nonlinear Euler solver will be used that solves the equations in generalized curvilinear coordinates using the chain rule transformation. This work, continuing earlier work on flat-plate cascades and Joukowski airfoils, will focus mainly on the effect of the grid and boundary conditions on the accuracy of the solution. The grids were generated using a commercially available grid generator, GridPro/az3000.
NASA Astrophysics Data System (ADS)
Jang, Juhi; Li, Fengyan; Qiu, Jing-Mei; Xiong, Tao
2015-01-01
In this paper, we develop a family of high order asymptotic preserving schemes for some discrete-velocity kinetic equations under a diffusive scaling, that in the asymptotic limit lead to macroscopic models such as the heat equation, the porous media equation, the advection-diffusion equation, and the viscous Burgers' equation. Our approach is based on the micro-macro reformulation of the kinetic equation which involves a natural decomposition of the equation to the equilibrium and non-equilibrium parts. To achieve high order accuracy and uniform stability as well as to capture the correct asymptotic limit, two new ingredients are employed in the proposed methods: discontinuous Galerkin (DG) spatial discretization of arbitrary order of accuracy with suitable numerical fluxes; high order globally stiffly accurate implicit-explicit (IMEX) Runge-Kutta scheme in time equipped with a properly chosen implicit-explicit strategy. Formal asymptotic analysis shows that the proposed scheme in the limit of ε → 0 is a consistent high order discretization for the limiting equation. Numerical results are presented to demonstrate the stability and high order accuracy of the proposed schemes together with their performance in the limit. Our methods are also tested for the continuous-velocity one-group transport equation in slab geometry and for several examples with spatially varying parameters.
NASA Astrophysics Data System (ADS)
Pan, Liang; Xu, Kun
2016-08-01
In this paper, for the first time a third-order compact gas-kinetic scheme is proposed on unstructured meshes for the compressible viscous flow computations. The possibility to design such a third-order compact scheme is due to the high-order gas evolution model, where a time-dependent gas distribution function at cell interface not only provides the fluxes across a cell interface, but also presents a time accurate solution for flow variables at cell interface. As a result, both cell averaged and cell interface flow variables can be used for the initial data reconstruction at the beginning of next time step. A weighted least-square procedure has been used for the initial reconstruction. Therefore, a compact third-order gas-kinetic scheme with the involvement of neighboring cells only can be developed on unstructured meshes. In comparison with other conventional high-order schemes, the current method avoids the Gaussian point integration for numerical fluxes along a cell interface and the multi-stage Runge-Kutta method for temporal accuracy. The third-order compact scheme is numerically stable under CFL condition CFL ≈ 0.5. Due to its multidimensional gas-kinetic formulation and the coupling of inviscid and viscous terms, even with unstructured meshes, the boundary layer solution and vortex structure can be accurately captured by the current scheme. At the same time, the compact scheme can capture strong shocks as well.
European Code against Cancer, 4th Edition: Cancer screening.
Armaroli, Paola; Villain, Patricia; Suonio, Eero; Almonte, Maribel; Anttila, Ahti; Atkin, Wendy S; Dean, Peter B; de Koning, Harry J; Dillner, Lena; Herrero, Rolando; Kuipers, Ernst J; Lansdorp-Vogelaar, Iris; Minozzi, Silvia; Paci, Eugenio; Regula, Jaroslaw; Törnberg, Sven; Segnan, Nereo
2015-12-01
In order to update the previous version of the European Code against Cancer and formulate evidence-based recommendations, a systematic search of the literature was performed according to the methodology agreed by the Code Working Groups. Based on the review, the 4th edition of the European Code against Cancer recommends: "Take part in organized cancer screening programmes for: Bowel cancer (men and women); Breast cancer (women); Cervical cancer (women)." Organized screening programs are preferable because they provide better conditions to ensure that the Guidelines for Quality Assurance in Screening are followed in order to achieve the greatest benefit with the least harm. Screening is recommended only for those cancers where a demonstrated life-saving effect substantially outweighs the potential harm of examining very large numbers of people who may otherwise never have, or suffer from, these cancers, and when an adequate quality of the screening is achieved. EU citizens are recommended to participate in cancer screening each time an invitation from the national or regional screening program is received and after having read the information materials provided and carefully considered the potential benefits and harms of screening. Screening programs in the European Union vary with respect to the age groups invited and to the interval between invitations, depending on each country's cancer burden, local resources, and the type of screening test used For colorectal cancer, most programs in the EU invite men and women starting at the age of 50-60 years, and from then on every 2 years if the screening test is the guaiac-based fecal occult blood test or fecal immunochemical test, or every 10 years or more if the screening test is flexible sigmoidoscopy or total colonoscopy. Most programs continue sending invitations to screening up to the age of 70-75 years. For breast cancer, most programs in the EU invite women starting at the age of 50 years, and not before the age
Computations of Flow Over a Hump Model Using Higher Order Method With Turbulence Modeling
NASA Technical Reports Server (NTRS)
Balakumar, Ponnampalam
2004-01-01
Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k-omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as they are observed in the experiments. The computed separation regions are much longer than that are observed. However, the percentage reduction in the separation region in the steady suction case is closer to that was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than that are measured and it points towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.
Computations of Flow over a Hump Model Using Higher Order Method with Turbulence Modeling
NASA Technical Reports Server (NTRS)
Balakumar, P.
2005-01-01
Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k - omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially. nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as those observed in the experiments. The computed separation regions are much longer than those observed experimentally. However, the percentage reduction in the separation region in the steady suction case is closer to what was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than the measured values pointing towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.
Direct calculations of waves in fluid flows using high-order compact difference scheme
NASA Technical Reports Server (NTRS)
Yu, Sheng-Tao; Hultgren, Lennart S.; Liu, Nan-Suey
1994-01-01
The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expression for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations. The numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the Courant-Friedrichs-Lewy (CFL) limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about six - eight grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 - 300 time steps. Three numerical examples of waves in compressible flow are included: (1) sound propagation in a duct with linear shear, (2) linear wave growth in a compressible free shear layer, and (3) vortex pairing in a compressible free shear layer perturbed at two frequencies.
Parallel Adjective High-Order CFD Simulations Characterizing SOFIA Cavity Acoustics
NASA Technical Reports Server (NTRS)
Barad, Michael F.; Brehm, Christoph; Kiris, Cetin C.; Biswas, Rupak
2016-01-01
This paper presents large-scale MPI-parallel computational uid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft fuselage of a Boeing 747SP. These simulations focus on how the unsteady ow eld inside and over the cavity interferes with the optical path and mounting structure of the telescope. A temporally fourth-order accurate Runge-Kutta, and spatially fth-order accurate WENO- 5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh re nement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32k CPU cores and 4 billion compu- tational cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregular numerical cost associated with blocks con- taining boundaries. Limits to scaling beyond 32k cores are identi ed, and targeted code optimizations are discussed.
NASA Astrophysics Data System (ADS)
Teyssier, Romain; Fromang, Sébastien; Dormy, Emmanuel
2006-10-01
We propose to extend the well-known MUSCL-Hancock scheme for Euler equations to the induction equation modeling the magnetic field evolution in kinematic dynamo problems. The scheme is based on an integral form of the underlying conservation law which, in our formulation, results in a “finite-surface” scheme for the induction equation. This naturally leads to the well-known “constrained transport” method, with additional continuity requirement on the magnetic field representation. The second ingredient in the MUSCL scheme is the predictor step that ensures second order accuracy both in space and time. We explore specific constraints that the mathematical properties of the induction equations place on this predictor step, showing that three possible variants can be considered. We show that the most aggressive formulations (referred to as C-MUSCL and U-MUSCL) reach the same level of accuracy as the other one (referred to as Runge Kutta), at a lower computational cost. More interestingly, these two schemes are compatible with the adaptive mesh refinement (AMR) framework. It has been implemented in the AMR code RAMSES. It offers a novel and efficient implementation of a second order scheme for the induction equation. We have tested it by solving two kinematic dynamo problems in the low diffusion limit. The construction of this scheme for the induction equation constitutes a step towards solving the full MHD set of equations using an extension of our current methodology.
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.
4th International Plant Biomechanics Conference Proceedings (Abstracts)
Frank W. Telewski; Lothar H. Koehler; Frank W. Ewers
2003-07-20
The 4th International Plant Biomechanics Conference facilitated an interdisciplinary exchange between scientists, engineers, and educators addressing the major questions encountered in the field of Plant Biomechanics. Subjects covered by the conference include: Evolution; Ecology; Mechanoreception; Cell Walls; Genetic Modification; Applied Biomechanics of Whole Plants, Plant Products, Fibers & Composites; Fluid Dynamics; Wood & Trees; Fracture Mechanics; Xylem Pressure & Water Transport; Modeling; and Introducing Plant Biomechanics in Secondary School Education.
NASA Astrophysics Data System (ADS)
Greene, Patrick T.; Eldredge, Jeff D.; Zhong, Xiaolin; Kim, John
2016-07-01
In this paper, we present a method for performing uniformly high-order direct numerical simulations of high-speed flows over arbitrary geometries. The method was developed with the goal of simulating and studying the effects of complex isolated roughness elements on the stability of hypersonic boundary layers. The simulations are carried out on Cartesian grids with the geometries imposed by a third-order cut-stencil method. A fifth-order hybrid weighted essentially non-oscillatory scheme was implemented to capture any steep gradients in the flow created by the geometries and a third-order Runge-Kutta method is used for time advancement. A multi-zone refinement method was also utilized to provide extra resolution at locations with expected complex physics. The combination results in a globally fourth-order scheme in space and third order in time. Results confirming the method's high order of convergence are shown. Two-dimensional and three-dimensional test cases are presented and show good agreement with previous results. A simulation of Mach 3 flow over the logo of the Ubuntu Linux distribution is shown to demonstrate the method's capabilities for handling complex geometries. Results for Mach 6 wall-bounded flow over a three-dimensional cylindrical roughness element are also presented. The results demonstrate that the method is a promising tool for the study of hypersonic roughness-induced transition.
a Numerical Comparison of Langrange and Kane's Methods of AN Arm Segment
NASA Astrophysics Data System (ADS)
Rambely, Azmin Sham; Halim, Norhafiza Ab.; Ahmad, Rokiah Rozita
A 2-D model of a two-link kinematic chain is developed using two dynamics equations of motion, namely Kane's and Lagrange Methods. The dynamics equations are reduced to first order differential equation and solved using modified Euler and fourth order Runge Kutta to approximate the shoulder and elbow joint angles during a smash performance in badminton. Results showed that Runge-Kutta produced a better and exact approximation than that of modified Euler and both dynamic equations produced better absolute errors.
NASA Astrophysics Data System (ADS)
Abdul Hakeem, A. K.; Vishnu Ganesh, N.; Ganga, B.
2015-05-01
The magnetic field effect on a steady two dimensional laminar radiative flow of an incompressible viscous water based nanofluid over a stretching/shrinking sheet with second order slip boundary condition is investigated both analytically and numerically. The governing partial differential equations are reduced to nonlinear ordinary differential equations by means of Lie symmetry group transformations. The dimensionless governing equations for this investigation are solved analytically using hyper-geometric function and numerically by the fourth order Runge-Kutta method with the shooting technique. A unique exact solution exists for momentum equation in stretching sheet case and dual solutions are obtained for shrinking sheet case which has upper and lower branches. It is found that the lower branch solution vanishes in the presence of higher magnetic field. The velocity and temperature profiles, the local skin friction coefficient and the reduced Nusselt number are examined and discussed for different spherical nanoparticles such as Au, Ag, Cu, Al, Al2 O3 and TiO2. A comparative study between the previously published results and the present analytical and numerical results for a special case is found to be in good agreement.
NASA Astrophysics Data System (ADS)
Yu, Cong
2011-03-01
The force-free (or low inertia) limit of magnetohydrodynamics (MHD) can be applied to many astrophysical objects, including black holes, neutron stars and accretion discs, where the electromagnetic field is so strong that the inertia and pressure of the plasma can be ignored. This is difficult to achieve with the standard MHD numerical methods because they still have to deal with plasma inertial terms even when these terms are much smaller than the electromagnetic terms. Under the force-free approximation, the plasma dynamics is entirely determined by the magnetic field. The plasma provides the currents and charge densities required by the dynamics of electromagnetic fields, but these currents carry no inertia. We present a high-order Godunov scheme to study such force-free electrodynamics. We have implemented weighted essentially non-oscillatory (WENO) spatial interpolations in our scheme. An exact Riemann solver is implemented, which requires spectral decomposition into characteristic waves. We advance the magnetic field with the constrained transport (CT) scheme to preserve the divergence-free condition to machine round-off error. We apply the third-order total variation diminishing (TVD) Runge-Kutta scheme for the temporal integration. The mapping from face-centred variables to volume-centred variables is carefully considered. Extensive testing are performed to demonstrate the ability of our scheme to address force-free electrodynamics correctly. We finally apply the scheme to study relativistic magnetically dominated tearing instabilities and neutron star magnetospheres.
A Review of High-Order and Optimized Finite-Difference Methods for Simulating Linear Wave Phenomena
NASA Technical Reports Server (NTRS)
Zingg, David W.
1996-01-01
This paper presents a review of high-order and optimized finite-difference methods for numerically simulating the propagation and scattering of linear waves, such as electromagnetic, acoustic, or elastic waves. The spatial operators reviewed include compact schemes, non-compact schemes, schemes on staggered grids, and schemes which are optimized to produce specific characteristics. The time-marching methods discussed include Runge-Kutta methods, Adams-Bashforth methods, and the leapfrog method. In addition, the following fourth-order fully-discrete finite-difference methods are considered: a one-step implicit scheme with a three-point spatial stencil, a one-step explicit scheme with a five-point spatial stencil, and a two-step explicit scheme with a five-point spatial stencil. For each method studied, the number of grid points per wavelength required for accurate simulation of wave propagation over large distances is presented. Recommendations are made with respect to the suitability of the methods for specific problems and practical aspects of their use, such as appropriate Courant numbers and grid densities. Avenues for future research are suggested.
Bifrost: A 4th Generation Launch Architecture Concept
NASA Astrophysics Data System (ADS)
Rohrschneider, R. R.; Young, D.; St.Germain, B.; Brown, N.; Crowley, J.; Maatsch, J.; Olds, J. R.
2002-01-01
A 4th generation launch architecture is studied for the purpose of drastically reducing launch costs and hence enabling new large mass missions such as space solar power and human exploration of other planets. The architecture consists of a magnetic levitation launch tube placed on the equator with the exit end elevated to approximately 20 km. Several modules exist for sending manned and unmanned payloads into Earth orbit. Analysis of the launch tube operations, launch trajectories, module aerodynamics, propulsion modules, and system costs are presented. Using the hybrid logistics module, it is possible to place payloads into low Earth orbit for just over 100 per lb.
High Order Difference Method for Low Mach Number Aeroacoustics
NASA Technical Reports Server (NTRS)
Mueller, B.; Yee, H. C.; Mansour, Nagi (Technical Monitor)
2001-01-01
A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved
Li, Wentao; Zhang, Dong H; Sun, Zhigang
2014-10-23
An efficient fourth-order split operator (named 4A6a in the main text), which was presented in the work by Blanes and Moan and was a partitioned Runge-Kutta method ( J. Comput. Appl. Math. 2002 , 142 , 313 ), is recommended for general usage in a reactive scattering calculation by the time-dependent quantum wavepacket method. This 4A6a propagator is constructed in a TVT form, that is, splitting in kinetic-potential-kinetic form, which is an optimal one among a series of higher-order split operators in examining with several typical triatomic reactive scattering processes, H + H2, H + H2(+), H + NH, H + O2, and F + HD reactions. A detailed comparison between the performances of higher-order split operators in the VTV form, that is, splitting in a potential-kinetic-potential, which was reported by Sun et al. ( Phys. Chem. Chem. Phys. 2012 , 14 , 1827 ), and in the TVT form reported in the current work suggests that the recommended 4A6a operator in the TVT form always has good numerical efficiency. This fact may suggest that this fourth propagator in the TVT form can be safely chosen without any further examination, at least among all of the higher-order split operators tested in this work, to apply in an efficient time-dependent wavepacket numerical calculation for describing a triatomic reactive scattering process. PMID:25268464
NASA Astrophysics Data System (ADS)
Motheau, E.; Abraham, J.
2016-05-01
A novel and efficient algorithm is presented in this paper to deal with DNS of turbulent reacting flows under the low-Mach-number assumption, with detailed chemistry and a quasi-spectral accuracy. The temporal integration of the equations relies on an operating-split strategy, where chemical reactions are solved implicitly with a stiff solver and the convection-diffusion operators are solved with a Runge-Kutta-Chebyshev method. The spatial discretisation is performed with high-order compact schemes, and a FFT based constant-coefficient spectral solver is employed to solve a variable-coefficient Poisson equation. The numerical implementation takes advantage of the 2DECOMP&FFT libraries developed by [1], which are based on a pencil decomposition method of the domain and are proven to be computationally very efficient. An enhanced pressure-correction method is proposed to speed up the achievement of machine precision accuracy. It is demonstrated that a second-order accuracy is reached in time, while the spatial accuracy ranges from fourth-order to sixth-order depending on the set of imposed boundary conditions. The software developed to implement the present algorithm is called HOLOMAC, and its numerical efficiency opens the way to deal with DNS of reacting flows to understand complex turbulent and chemical phenomena in flames.
A 95 GHz, 4th harmonic gyro-oscillator
Hargreaves, T.A.; Scheitrum, G.P.; Bemis, T.; Higgins, L.
1994-12-31
There is currently an interest in medium power ({approximately}100 kW), compact 95 GHz amplifiers for future radar applications. Size, weight, and efficiency are critical for airborne applications. Litton has been investigating a 4th harmonic, 4-cavity gyro-amplifier. The key to success of the amplifier is the axis-encircling electron beam from a new type of electron gun, the advanced center post (ACP) gun. Gun simulations incorporating the actual magnetic field and thermal velocity spread in the emitted electrons show that axial velocity spreads of less than 2% are attainable, which is significantly better than other gun concepts. The amplifier utilizes coaxial-magnetron-type cavities operating in the {pi} mode. In this cavity, vanes extend nearly down to the electron beam`s outside diameter. The majority of the RF stored energy in the system is in the coaxial cavity, so that the resonant frequency and quality factor of each coaxial magnetron cavity may be adjusted by varying only the coaxial cavity. Several components are being tested individually. To test the cavity design, a 4th harmonic oscillator based on a coaxial magnetron cavity has been designed. Results of the oscillator testing will be presented.
Some numerical methods for integrating systems of first-order ordinary differential equations
NASA Technical Reports Server (NTRS)
Clark, N. W.
1969-01-01
Report on numerical methods of integration includes the extrapolation methods of Bulirsch-Stoer and Neville. A comparison is made nith the Runge-Kutta and Adams-Moulton methods, and circumstances are discussed under which the extrapolation method may be preferred.
Federal Register 2010, 2011, 2012, 2013, 2014
2011-06-29
... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zones; July 4th Weekend Fireworks Displays Within... under Executive Order 13132, Federalism, if it has a substantial direct effect on State or local governments and would either preempt State law or impose a substantial direct cost of compliance on them....
ERIC Educational Resources Information Center
Nam, Younkyeong; Karahan, Engin; Roehrig, Gillian
2016-01-01
Geologic time scale is a very important concept for understanding long-term earth system events such as climate change. This study examines forty-three 4th-8th grade Native American--particularly Ojibwe tribe--students' understanding of relative ordering and absolute time of Earth's significant geological and biological events. This study also…
Kinetic and reactor models for HDT of middle distillates
Cotta, R.M.; Filho, R.M.
1996-12-31
Hydrodesulfurization (HDS) and hydrodenitrogenation (HDN) of middle distillates over a commercial Ni-Mo/y-Al{sub 2}O{sub 3} has been studied under wide operating conditions just as 340 to 380{degrees}C and 38 to 98 atm. A Power Law model was presented to each one of those reactions. The parameters of kinetic equations were estimated solving the ordinary differential equations by the 4 order Runge-Kutta-Gill algorithm and Marquardt method for searching of set of kinetic parameters (kinetic constants as well as the orders of reactions). An adiabatic diesel hydrotreating trickle-bed reactor packed with the same catalyst was simulated numerically in order to check up the behavior of this specific reaction system. One dimensional pseudo-homogeneous model was used in this work. For each feed, the mass and energy balance equations were integrated along the length of the catalytic bed using the 4th Runge-Kutta-Gill method. The performance of two industrial reactors was checked. 5 refs., 2 tabs.
Special Issue: 4th International Workshop on Space Radiation (IWSRR)
NASA Technical Reports Server (NTRS)
Cucinotta, Francis A.
2007-01-01
This special issue of the journal "Radiation and Environmental Biophysics" contains 20 peer-reviewed papers contributed by leading space radiation researcher's world-wide attending the 4th IWSRR. Manuscripts cover a broad range of topics ranging from radiation environments and transport in shielding and planetary surfaces to new results in understanding the biological effects of protons and high-charge and energy (HZE) nuclei on the risk of cancer, and degenerative diseases such as central nervous system effects, heart disease, and cataracts. The issue provides a snapshot of the state-of-the-art of the research in this field, demonstrating both the important results gathered in the past few years with experiments at accelerators, and the need for more research to quantify the risk and develop countermeasures.
[Time--the 4th dimision in medicine and psychotherapy].
Bergmann, Günther
2003-01-01
Time is presented as well in his historical meaning and as 4th dimension in its medical and psychotherapeutic context. In this medical and psychotherapeutic process it has an important function and is a variable of a process procedure. The difference between "kairos" = (the right point of time) and "chronos" = (the period of time) is historically meanful. The subjective experienced time is as well emphasized by the development of time in the relation to the development of the "self" as in the subjective experience of time in medical and psychotherapeutic situations. There are also changed conceptions and understandings of time running parallel to the development of nature sciences. The importance of time is explained for the medical practice and the meeting with the patient--especially for chronic diseases. The connection of confidence and time is particularly emphasized in the systemic approach. PMID:12764877
European Code against Cancer 4th Edition: Breastfeeding and cancer.
Scoccianti, Chiara; Key, Timothy J; Anderson, Annie S; Armaroli, Paola; Berrino, Franco; Cecchini, Michele; Boutron-Ruault, Marie-Christine; Leitzmann, Michael; Norat, Teresa; Powers, Hilary; Schüz, Joachim; Wiseman, Martin; Romieu, Isabelle
2015-12-01
Breast cancer is the most frequent cancer in women, and incidence rates have been rising in European Union (EU) countries over recent decades due in part to a sharp decline in breastfeeding practices. Evidence for a protective association between breastfeeding and the risk of breast cancer at all ages is convincing, and modest protective relationships between breastfeeding and the risk of endometrial and ovarian cancers have been suggested. The reduction in breast cancer risk is estimated at 2% for an increase of 5 months of lifetime breastfeeding. The longer women breastfeed, the more they are protected against breast cancer. In addition, breastfeeding is associated with several health benefits for both the mother and the breastfed child. Taking all this evidence into account, the 4th edition of the European Code against Cancer recommends: "Breastfeeding reduces the mother's cancer risk. If you can, breastfeed your baby". PMID:26116994
The Epilepsy Foundation's 4th Biennial Epilepsy Pipeline Update Conference.
French, Jacqueline A; Schachter, Steven C; Sirven, Joseph; Porter, Roger
2015-05-01
On June 5 and 6, 2014, the Epilepsy Foundation held its 4th Biennial Epilepsy Pipeline Update Conference, an initiative of the Epilepsy Therapy Project, which showcased the most promising epilepsy innovations from health-care companies and academic laboratories dedicated to pioneering and advancing drugs, biologics, technologies, devices, and diagnostics for epilepsy. Speakers and attendees included emerging biotech and medical technology companies, major pharmaceutical and device companies, as well as investigators and innovators at the cutting-edge of epilepsy. The program included panel discussions on collaboration between small and large companies, how to get products in need of funding to the marketplace, who is currently funding epilepsy and CNS innovation, and how the NIH facilitates early-stage drug development. Finally, the conference featured the third annual "Shark Tank" competition. The presentations are summarized in this paper, which is followed by a compilation of the meeting poster abstracts. PMID:25922152
A second-order parareal algorithm for fractional PDEs
NASA Astrophysics Data System (ADS)
Wu, Shu-Lin
2016-02-01
We are concerned with using the parareal (parallel-in-time) algorithm for large scale ODEs system U‧ (t) + AU (t) + dAα U (t) = F (t) arising frequently in semi-discretizing time-dependent PDEs with spatial fractional operators, where d > 0 is a constant, α ∈ (0 , 1) and A is a spare and symmetric positive definite (SPD) matrix. The parareal algorithm is iterative and is characterized by two propagators F and G, which are respectively associated with small temporal mesh size Δt and large temporal mesh size ΔT. The two mesh sizes satisfy ΔT = JΔt with J ≥ 2 being an integer, which is called mesh ratio. Let Tunitf and Tunitg be respectively the computational cost of the two propagators for moving forward one time step. Then, it is well understood that the speedup of the parareal algorithm, namely E, satisfies E = O (clog (1 / ρ)) , where c : = Tunitf / Tunitg and ρ is the convergence factor. A larger E corresponds a more efficient parareal solver. For G = Backward-Euler and some choices of F, previous studies show that ρ can be a satisfactory quantity. Particularly, for F = 2nd-order DIRK (diagonally implicit Runge-Kutta), it holds ρ ≈1/3 for any choice of the mesh ratio J. In this paper, we continue to consider F = 2nd-order DIRK, but with a new choice for G, the IMEX (implicit-explicit) Euler method, where the 'implicit' and 'explicit' computation is respectively associated with A and dAα. Compared to the widely used Backward-Euler method, this choice on the one hand increases c (this point is apparent), and interestingly on the other hand it can also make the convergence factor ρ smaller: ρ ≈1/5! Numerical results are provided to support our conclusions.
High-Order Space-Time Methods for Conservation Laws
NASA Technical Reports Server (NTRS)
Huynh, H. T.
2013-01-01
Current high-order methods such as discontinuous Galerkin and/or flux reconstruction can provide effective discretization for the spatial derivatives. Together with a time discretization, such methods result in either too small a time step size in the case of an explicit scheme or a very large system in the case of an implicit one. To tackle these problems, two new high-order space-time schemes for conservation laws are introduced: the first is explicit and the second, implicit. The explicit method here, also called the moment scheme, achieves a Courant-Friedrichs-Lewy (CFL) condition of 1 for the case of one-spatial dimension regardless of the degree of the polynomial approximation. (For standard explicit methods, if the spatial approximation is of degree p, then the time step sizes are typically proportional to 1/p(exp 2)). Fourier analyses for the one and two-dimensional cases are carried out. The property of super accuracy (or super convergence) is discussed. The implicit method is a simplified but optimal version of the discontinuous Galerkin scheme applied to time. It reduces to a collocation implicit Runge-Kutta (RK) method for ordinary differential equations (ODE) called Radau IIA. The explicit and implicit schemes are closely related since they employ the same intermediate time levels, and the former can serve as a key building block in an iterative procedure for the latter. A limiting technique for the piecewise linear scheme is also discussed. The technique can suppress oscillations near a discontinuity while preserving accuracy near extrema. Preliminary numerical results are shown
Comparative study of numerical schemes of TVD3, UNO3-ACM and optimized compact scheme
NASA Technical Reports Server (NTRS)
Lee, Duck-Joo; Hwang, Chang-Jeon; Ko, Duck-Kon; Kim, Jae-Wook
1995-01-01
Three different schemes are employed to solve the benchmark problem. The first one is a conventional TVD-MUSCL (Monotone Upwind Schemes for Conservation Laws) scheme. The second scheme is a UNO3-ACM (Uniformly Non-Oscillatory Artificial Compression Method) scheme. The third scheme is an optimized compact finite difference scheme modified by us: the 4th order Runge Kutta time stepping, the 4th order pentadiagonal compact spatial discretization with the maximum resolution characteristics. The problems of category 1 are solved by using the second (UNO3-ACM) and third (Optimized Compact) schemes. The problems of category 2 are solved by using the first (TVD3) and second (UNO3-ACM) schemes. The problem of category 5 is solved by using the first (TVD3) scheme. It can be concluded from the present calculations that the Optimized Compact scheme and the UN03-ACM show good resolutions for category 1 and category 2 respectively.
76 FR 37649 - Safety Zone; Northern California Annual Fireworks Events, July 4th Fireworks Display
Federal Register 2010, 2011, 2012, 2013, 2014
2011-06-28
... SECURITY Coast Guard 33 CFR Part 165 Safety Zone; Northern California Annual Fireworks Events, July 4th Fireworks Display AGENCY: Coast Guard, DHS. ACTION: Notice of enforcement of regulation. SUMMARY: The Coast Guard will enforce safety zone for the annual July 4th Fireworks Display (Tahoe City 4th of...
77 FR 39172 - Safety Zone: Skagway Harbor, Skagway, AK for 4th of July Fireworks
Federal Register 2010, 2011, 2012, 2013, 2014
2012-07-02
... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zone: Skagway Harbor, Skagway, AK for 4th of July... 4th of July Fireworks display. Persons desiring to transit within this safety zone must contact the... Inlet for the City of Skagway, Alaska, sponsored 4th of July fireworks display. A no-action...
PREFACE: 4th International Symposium on Functional Materials (ISFM2011)
NASA Astrophysics Data System (ADS)
Yin, Shu; Sekino, Tohru; Tanaka, Shun-ichiro; Sato, Tsugio; Lu, Li; Xue, Dongfeng
2012-01-01
The 4th International Symposium on Functional Materials (ISFM2011) was held in Sendai, Japan, on 2-6 August 2011. This Special Issue of Journal of Physics: Conference Series (JPCS) consists of partial manuscripts which were presented at ISFM2011. Advanced materials have experienced a dramatic increase in demand for research, development and applications. The aim of the International Symposium on Functional Materials (ISFM) was to provide an overview of the present status with historical background and to foresee future trends in the field of functional materials. The 4th symposium, ISFM 2011, covered a wide variety of topics within state-of-the-art advanced materials science and technology, and focused especially on four major categories including: Environmental Materials, Electronic Materials, Energy Materials and Biomedical Materials. As you know, a massive earthquake and the Tsunami that followed occurred near the Tohoku region on 11 March 2011. After the earthquake, although there were many difficulties in continuing to organize the symposium, we received warm encouragement from many researchers and societies, especially from the members of the International Advisory Committee and Organizing Committee, so that ISFM2011 could be held on schedule. We are honored that ISFM2011 was the first formal international academic conference held in the Tohoku area of Japan after the 11 March earthquake. About 140 participants from 14 countries took part in the ISFM2011 symposium, which included five plenary talks by world-leading scientists, 32 invited talks, and many oral and poster presentations. We are delighted to see that many researchers are interested in the synthesis and the properties as well as the applications of functional materials. Many fruitful and exciting research achievements were presented in the symposium. We believe that this symposium provided a good chance for scientists to communicate and exchange opinions with each other. We would also like to
IMEX-a : an adaptive, fifth order implicit-explicit integration scheme.
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.
NASA Astrophysics Data System (ADS)
Belmonte, Juan Antonio
2015-08-01
The pyramids of Egypt, notably those of the 4th Dinasty as Giza, have always be considered an unmistikable part of human world heritage as the only surviving wonders of the Ancient World. Their majesty, technical hability and innovative character have always beeen considered as representative of ancient Egyptian ingenuity. However, past and present fringe theories about the pyramids and astronomy have always polluted the role of our discipline in the design, construction and symbolism of these impressive monuments. This is indeed unfear. Fortunately, things have started to change in the last couple of decades and now astronomy is interpreted as a neccessary tool for the correct interpretation of the astral eschatology present in the 5th and 6th Dynasty Texts of the Pyramids. Although the pyramid complexes of the 4th Dynasty are mute, there is however recent research showing that a strong astral symbolism could be hidden in many aspects of the complex architecture and in the design of these exceptional monuments. This idea comes from several hints obtained not only from planning and construction, but also from epigraphy and the analysis of celestial and local landscapes. Chronology also plays a most relevant role on this. The pyramid complexes of the 4th Dynasty at Meidum, Dahshur, Giza and Abu Rowash -- all of which enjoy UNESCO World Heritage recognition -- willl be scrutinized. As a consequence, we will show how astronomy can certainly enhance the face value of these extraordinary monuments as a definitive proof of the ancient Egyptian quest for Ma'at, i.e. their perennial obsesion for Cosmic Order.
Managing haemophilia for life: 4th Haemophilia Global Summit.
Astermark, J; Dolan, G; Hilberg, T; Jiménez-Yuste, V; Laffan, M; Lassila, R; Lobet, S; Martinoli, C; Perno, C-F
2014-07-01
The 4th Haemophilia Global Summit was held in Potsdam, Germany, in September 2013 and brought together an international faculty of haemophilia experts and delegates from multidisciplinary backgrounds. The programme was designed by an independent Scientific Steering Committee of haemophilia experts and explored global perspectives in haemophilia care, discussing practical approaches to the optimal management of haemophilia now and in the future. The topics outlined in this supplement were selected by the Scientific Steering Committee for their relevance and potential to influence haemophilia care globally. In this supplement from the meeting, Jan Astermark reviews current understanding of risk factors for the development of inhibitory antibodies and discusses whether this risk can be modulated and minimized. Factors key to the improvement of joint health in people with haemophilia are explored, with Carlo Martinoli and Víctor Jiménez-Yuste discussing the utility of ultrasound for the early detection of haemophilic arthropathy. Other aspects of care necessary for the prevention and management of joint disease in people with haemophilia are outlined by Thomas Hilberg and Sébastian Lobet, who highlight the therapeutic benefits of physiotherapy and sports therapy. Riitta Lassila and Carlo-Federico Perno describe current knowledge surrounding the risk of transmission of infectious agents via clotting factor concentrates. Finally, different types of extended half-life technology are evaluated by Mike Laffan, with a focus on the practicalities and challenges associated with these products. PMID:24924596
European Code against Cancer 4th Edition: Diet and cancer.
Norat, Teresa; Scoccianti, Chiara; Boutron-Ruault, Marie-Christine; Anderson, Annie; Berrino, Franco; Cecchini, Michele; Espina, Carolina; Key, Tim; Leitzmann, Michael; Powers, Hilary; Wiseman, Martin; Romieu, Isabelle
2015-12-01
Lifestyle factors, including diet, have long been recognised as potentially important determinants of cancer risk. In addition to the significant role diet plays in affecting body fatness, a risk factor for several cancers, experimental studies have indicated that diet may influence the cancer process in several ways. Prospective studies have shown that dietary patterns characterised by higher intakes of fruits, vegetables, and whole-grain foods, and lower intakes of red and processed meats and salt, are related to reduced risks of death and cancer, and that a healthy diet can improve overall survival after diagnosis of breast and colorectal cancers. There is evidence that high intakes of fruit and vegetables may reduce the risk of cancers of the aerodigestive tract, and the evidence that dietary fibre protects against colorectal cancer is convincing. Red and processed meats increase the risk of colorectal cancer. Diets rich in high-calorie foods, such as fatty and sugary foods, may lead to increased calorie intake, thereby promoting obesity and leading to an increased risk of cancer. There is some evidence that sugary drinks are related to an increased risk of pancreatic cancer. Taking this evidence into account, the 4th edition of the European Code against Cancer recommends that people have a healthy diet to reduce their risk of cancer: they should eat plenty of whole grains, pulses, vegetables and fruits; limit high-calorie foods (foods high in sugar or fat); avoid sugary drinks and processed meat; and limit red meat and foods high in salt. PMID:26164653
Adaptive Numerical Dissipative Control in High Order Schemes for Multi-D Non-Ideal MHD
NASA Technical Reports Server (NTRS)
Yee, H. C.; Sjoegreen, B.
2004-01-01
derivatives and a fourth-order Runge-Kutta method are denoted.
miR-155 Inhibition Sensitizes CD4+ Th Cells for TREG Mediated Suppression
Rust, Werner; Labhart, Paul; Alexiadis, Vassili; Becker, Christian; Hafner, Mathias; Weith, Andreas; Lenter, Martin C.; Jonuleit, Helmut; Schmitt, Edgar; Mennerich, Detlev
2009-01-01
Background In humans and mice naturally occurring CD4+CD25+ regulatory T cells (nTregs) are a thymus-derived subset of T cells, crucial for the maintenance of peripheral tolerance by controlling not only potentially autoreactive T cells but virtually all cells of the adaptive and innate immune system. Recent work using Dicer-deficient mice irrevocably demonstrated the importance of miRNAs for nTreg cell-mediated tolerance. Principal Findings DNA-Microarray analyses of human as well as murine conventional CD4+ Th cells and nTregs revealed a strong up-regulation of mature miR-155 (microRNA-155) upon activation in both populations. Studying miR-155 expression in FoxP3-deficient scurfy mice and performing FoxP3 ChIP-Seq experiments using activated human T lymphocytes, we show that the expression and maturation of miR-155 seem to be not necessarily regulated by FoxP3. In order to address the functional relevance of elevated miR-155 levels, we transfected miR-155 inhibitors or mature miR-155 RNAs into freshly-isolated human and mouse primary CD4+ Th cells and nTregs and investigated the resulting phenotype in nTreg suppression assays. Whereas miR-155 inhibition in conventional CD4+ Th cells strengthened nTreg cell-mediated suppression, overexpression of mature miR-155 rendered these cells unresponsive to nTreg cell-mediated suppression. Conclusion Investigation of FoxP3 downstream targets, certainly of bound and regulated miRNAs revealed the associated function between the master regulator FoxP3 and miRNAs as regulators itself. miR-155 is shown to be crucially involved in nTreg cell mediated tolerance by regulating the susceptibility of conventional human as well as murine CD4+ Th cells to nTreg cell-mediated suppression. PMID:19777054
European Code against Cancer, 4th Edition: Tobacco and cancer.
Leon, Maria E; Peruga, Armando; McNeill, Ann; Kralikova, Eva; Guha, Neela; Minozzi, Silvia; Espina, Carolina; Schüz, Joachim
2015-12-01
Tobacco use, and in particular cigarette smoking, is the single largest preventable cause of cancer in the European Union (EU). All tobacco products contain a wide range of carcinogens. The main cancer-causing agents in tobacco smoke are polycyclic aromatic hydrocarbons, tobacco-specific N-nitrosamines, aromatic amines, aldehydes, and certain volatile organic compounds. Tobacco consumers are also exposed to nicotine, leading to tobacco addiction in many users. Cigarette smoking causes cancer in multiple organs and is the main cause of lung cancer, responsible for approximately 82% of cases. In 2012, about 313,000 new cases of lung cancer and 268,000 lung cancer deaths were reported in the EU; 28% of adults in the EU smoked tobacco, and the overall prevalence of current use of smokeless tobacco products was almost 2%. Smokeless tobacco products, a heterogeneous category, are also carcinogenic but cause a lower burden of cancer deaths than tobacco smoking. One low-nitrosamine product, snus, is associated with much lower cancer risk than other smokeless tobacco products. Smoking generates second-hand smoke (SHS), an established cause of lung cancer, and inhalation of SHS by non-smokers is still common in indoor workplaces as well as indoor public places, and more so in the homes of smokers. Several interventions have proved effective for stopping smoking; the most effective intervention is the use of a combination of pharmacotherapy and behavioural support. Scientific evidence leads to the following two recommendations for individual action on tobacco in the 4th edition of the European Code Against Cancer: (1) "Do not smoke. Do not use any form of tobacco"; (2) "Make your home smoke-free. Support smoke-free policies in your workplace". PMID:26272517
PREFACE: 4th International Hadron Physics Conference (TROIA'14)
NASA Astrophysics Data System (ADS)
Dağ, Hüseyin; Erkol, Güray; Küçükarslan, Ayşe; Özpineci, Altuğ
2014-11-01
The 4th International Conference on Hadron Physics, TROIA'14, was held at Canakkale, Turkey on 1-5 July 2014. Ozyegin University, Middle East Technical University, Canakkale Onsekiz Mart University, Turkish Atomic Energy Authority and HadronPhysics2 Consortium sponsored the conference. It aimed at bringing together the experts and the young scientists working on experimental and theoretical hadron physics. About 50 participants from 10 countries attended the conference. The topics covered included: . Chiral Perturbation Theory . QCD Sum Rules . Effective Field Theory . Exotic Hadrons . Hadron Properties from Lattice QCD . Experimental Results and Future Perspectives . Hadronic Distribution Amplitudes The conference presentations were organized such that the morning sessions contained invited talks and afternoon sessions were devoted to contributed talks. The speakers of the invited talks were: C. Alexandrou, A. Gal, L. Tolos, J.R. Pelaez and M. Schindler. We had also guest speakers D. A. Demir and T. Senger. The conference venue was a resort hotel around Canakkale. As a social program, a guided full-day excursion to the excavation site of the ancient Troia town and Assos was organized. We believe that this conference provided a medium for young scientists and experts in the field to effectively communicate and share ideas. We would like to express our sincere thanks to supporting agencies and to all participants for their contributions and stimulating discussions. We are also grateful to the Scientific Secretary, Bora Işıldak, and all other members of the Organizing Committee for their patience and efforts. 30.10.2014 The Editors
High order asymptotic preserving nodal discontinuous Galerkin IMEX schemes for the BGK equation
NASA Astrophysics Data System (ADS)
Xiong, Tao; Jang, Juhi; Li, Fengyan; Qiu, Jing-Mei
2015-03-01
In this paper, we develop high-order asymptotic preserving (AP) schemes for the BGK equation in a hyperbolic scaling, which leads to the macroscopic models such as the Euler and compressible Navier-Stokes equations in the asymptotic limit. Our approaches are based on the so-called micro-macro formulation of the kinetic equation which involves a natural decomposition of the problem to the equilibrium and the non-equilibrium parts. The proposed methods are formulated for the BGK equation with constant or spatially variant Knudsen number. The new ingredients for the proposed methods to achieve high order accuracy are the following: we introduce discontinuous Galerkin (DG) discretization of arbitrary order of accuracy with nodal Lagrangian basis functions in space; we employ a high order globally stiffly accurate implicit-explicit (IMEX) Runge-Kutta (RK) scheme as time discretization. Two versions of the schemes are proposed: Scheme I is a direct formulation based on the micro-macro decomposition of the BGK equation, while Scheme II, motivated by the asymptotic analysis for the continuous problem, utilizes certain properties of the projection operator. Compared with Scheme I, Scheme II not only has better computational efficiency (the computational cost is reduced by half roughly), but also allows the establishment of a formal asymptotic analysis. Specifically, it is demonstrated that when 0 < ε ≪ 1, Scheme II, up to O (ε2), becomes a local DG discretization with an explicit RK method for the macroscopic compressible Navier-Stokes equations, a method in a similar spirit to the ones in Bassi and Rebay (1997) [3], Cockburn and Shu (1998) [16]. Numerical results are presented for a wide range of Knudsen number to illustrate the effectiveness and high order accuracy of the methods.
NASA Astrophysics Data System (ADS)
Schaal, Kevin; Bauer, Andreas; Chandrashekar, Praveen; Pakmor, Rüdiger; Klingenberg, Christian; Springel, Volker
2015-11-01
Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in current astrophysical codes, especially also in light of evolving computer technology, which favours certain computational approaches over others. Here we introduce the new adaptive mesh refinement (AMR) code TENET, which employs a high-order discontinuous Galerkin (DG) scheme for hydrodynamics. The Euler equations in this method are solved in a weak formulation with a polynomial basis by means of explicit Runge-Kutta time integration and Gauss-Legendre quadrature. This approach offers significant advantages over commonly employed second-order finite-volume (FV) solvers. In particular, the higher order capability renders it computationally more efficient, in the sense that the same precision can be obtained at significantly less computational cost. Also, the DG scheme inherently conserves angular momentum in regions where no limiting takes place, and it typically produces much smaller numerical diffusion and advection errors than an FV approach. A further advantage lies in a more natural handling of AMR refinement boundaries, where a fall-back to first order can be avoided. Finally, DG requires no wide stencils at high order, and offers an improved data locality and a focus on local computations, which is favourable for current and upcoming highly parallel supercomputers. We describe the formulation and implementation details of our new code, and demonstrate its performance and accuracy with a set of two- and three-dimensional test problems. The results confirm that DG schemes have a high potential for astrophysical applications.
Overview of the NASA Glenn Flux Reconstruction Based High-Order Unstructured Grid Code
NASA Technical Reports Server (NTRS)
Spiegel, Seth C.; DeBonis, James R.; Huynh, H. T.
2016-01-01
A computational fluid dynamics code based on the flux reconstruction (FR) method is currently being developed at NASA Glenn Research Center to ultimately provide a large- eddy simulation capability that is both accurate and efficient for complex aeropropulsion flows. The FR approach offers a simple and efficient method that is easy to implement and accurate to an arbitrary order on common grid cell geometries. The governing compressible Navier-Stokes equations are discretized in time using various explicit Runge-Kutta schemes, with the default being the 3-stage/3rd-order strong stability preserving scheme. The code is written in modern Fortran (i.e., Fortran 2008) and parallelization is attained through MPI for execution on distributed-memory high-performance computing systems. An h- refinement study of the isentropic Euler vortex problem is able to empirically demonstrate the capability of the FR method to achieve super-accuracy for inviscid flows. Additionally, the code is applied to the Taylor-Green vortex problem, performing numerous implicit large-eddy simulations across a range of grid resolutions and solution orders. The solution found by a pseudo-spectral code is commonly used as a reference solution to this problem, and the FR code is able to reproduce this solution using approximately the same grid resolution. Finally, an examination of the code's performance demonstrates good parallel scaling, as well as an implementation of the FR method with a computational cost/degree- of-freedom/time-step that is essentially independent of the solution order of accuracy for structured geometries.
A laboratory model of post-Newtonian gravity with high power lasers and 4th generation light sources
NASA Astrophysics Data System (ADS)
Gregori, G.; Levy, M. C.; Wadud, M. A.; Crowley, B. J. B.; Bingham, R.
2016-04-01
Using the post-Newtonian formalism of gravity, we attempt to calculate the x-ray Thomson scattering cross section of electrons that are accelerated in the field of a high intensity optical laser. We show that our results are consistent with previous calculations, suggesting that the combination of high power laser and 4th generation light sources may become a powerful platform to test models exploring high order corrections to the Newtonian gravity.
General Chemistry Collection for Students (CD-ROM), Abstract of Special Issue 16, 4th Edition
NASA Astrophysics Data System (ADS)
2000-07-01
bookstore. The cost per CD can be quite low when large numbers are ordered (as little as $3 each), making this a cost-effective method of allowing students access to the software they need whenever and wherever they desire. Other JCE Software CDs can also be adopted. Network licenses to distribute the software to your students via your local campus network can also be arranged. Contact us for details on purchasing multiple user licenses. Price and Ordering An order form is inserted in this issue that provides prices and other ordering information. If this card is not available or if you need additional information, contact: JCE Software, University of Wisconsin-Madison, 1101 University Avenue, Madison, WI 53706-1396; phone; 608/262-5153 or 800/991-5534; fax: 608/265-8094; email: jcesoft@chem.wisc.edu. Table 1. Contents of the General Chemistry Collection, 4th Edition
Quark masses and mixings in the RS1 model with a condensing 4th generation
NASA Astrophysics Data System (ADS)
Hernández, A. E. Cárcamo; Dib, Claudio O.; Neill, Nicolás A.; Zerwekh, Alfonso R.
2012-02-01
We study the hierarchy of quark masses and mixings in a model based on a 5-dimensional spacetime with constant curvature of Randall-Sundrum type with two branes, where the Electroweak Symmetry Breaking is caused dynamically by the condensation of a 4th generation of quarks, due to underlying physics from the 5D bulk and the first KK gluons. We first study the hierarchy of quark masses and mixings that can be obtained from purely adjusting the profile localizations, finding that realistic masses are not reproduced unless non trivial hierarchies of underlying 4-fermion interactions from the bulk are included. Then we study global U(1) symmetries that can be imposed in order to obtain non-symmetric modified Fritzsch-like textures in the mass matrices that reproduce reasonably well quark masses and CKM mixings.
Urban Infrasound Observations - Examples from July 4th 2012
NASA Astrophysics Data System (ADS)
McComas, S.; Hayward, C.; Golden, P.; McKenna, M.; Simpson, C.
2012-12-01
, the Heroy Building Rooftop Array, is a two-element 30m line on a single rooftop. Large-scale fireworks displays in Dallas on 4 July 2012 provided an opportunity to identify and characterize known signals in an urban setting. The identified events were associated with one of these fireworks displays about 2 km from the arrays. Signals from these sources were used to tune processing parameters for an automatic coherent detection process, Progressive Multichannel Correlation Method (PMCC). PMCC was then used to scan the data for all possible firework sources in the urban environment and determine temporal, back azimuth, apparent velocity, and frequency information about the sources. The signal frequencies seen were 10-80 Hz and documented the details of the nearly 30 minute firework show. The resulting PMCC analysis showed potential to effectively identify other, lower frequency sources in the urban environment. These data were also is used to characterize the noise environment. Significant roof-to-roof noise differences may be related to the building configurations and mechanical equipment, as well as the interactions of the winds with the structures. During the evening of July 4th , additional ground deployed infrasound gauges provided a comparison of free surface and rooftop measurements. Permission to publish was granted by Director, Geotechnical and Structures Laboratory.
75 FR 35649 - Safety Zone; Northern California Annual Fireworks Events, July 4th Fireworks Display
Federal Register 2010, 2011, 2012, 2013, 2014
2010-06-23
... SECURITY Coast Guard 33 CFR Part 165 Safety Zone; Northern California Annual Fireworks Events, July 4th Fireworks Display AGENCY: Coast Guard, DHS. ACTION: Notice of enforcement of regulation. SUMMARY: The Coast Guard will enforce the Tahoe City 4th of July Fireworks Display safety zone, from 9 a.m. through 10...
The Effects of Cooperative Learning Strategies on Vocabulary Skills of 4th Grade Students
ERIC Educational Resources Information Center
Bilen, Didem; Tavil, Zekiye Müge
2015-01-01
This study was carried out to investigate the effects of cooperative learning strategies on the vocabulary skills of 4th grade students. The study was also designed to ascertain the attitudes of the students in the experimental group towards cooperative learning. Out of 96 4th grade students enrolled in the private school where the study took…
The school nutrition program's role in weight management of 4th grade elementary students
Technology Transfer Automated Retrieval System (TEKTRAN)
We are attempting to uncover the school nutrition program's role in weight management of 4th grade elementary students. Data was collected within a time frame for the food frequency questionnaire (FFQ) set at two months at the WT Cheney Elementary School and South Wood Elementary for 4th grade stud...
75 FR 34636 - Safety Zone; Jameson Beach 4th of July Fireworks Display
Federal Register 2010, 2011, 2012, 2013, 2014
2010-06-18
... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zone; Jameson Beach 4th of July Fireworks Display... temporary safety zone in the navigable waters of Lake Tahoe, for the Jameson Beach 4th of July Fireworks... has a substantial direct effect on State or local governments and would either preempt State law...
75 FR 34639 - Safety Zone; Reedville July 4th Celebration, Cockrell's Creek, Reedville, VA
Federal Register 2010, 2011, 2012, 2013, 2014
2010-06-18
... Celebration, Cockrell's Creek, Reedville, VA in the Federal Register (75 FR 26157). We received no comments on... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zone; Reedville July 4th Celebration, Cockrell's... the Reedville July 4th Celebration event. This action is intended to restrict vessel traffic...
76 FR 37650 - Safety Zone; 4th of July Festival Berkeley Marina Fireworks Display Berkeley, CA
Federal Register 2010, 2011, 2012, 2013, 2014
2011-06-28
... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zone; 4th of July Festival Berkeley Marina... Berkeley Pier, Berkeley, CA in support of the 4th of July Festival Berkeley Marina Fireworks Display... used in the fireworks display. Background and Purpose The City of Berkeley Marina will sponsor the...
75 FR 26157 - Safety Zone; Reedville July 4th Celebration, Cockrell's Creek, Reedville, VA
Federal Register 2010, 2011, 2012, 2013, 2014
2010-05-11
... Federal Register (73 FR 3316). Public Meeting We do not now plan to hold a public meeting. But you may... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zone; Reedville July 4th Celebration, Cockrell's..., VA in support of the Reedville July 4th Celebration event. This action is intended to restrict...
Science Content Courses: Workshop in Food Chemistry for 4th Grade School Teachers
ERIC Educational Resources Information Center
Chaiyapechara, S.; Dong, F. M.
2004-01-01
A science content course in food chemistry was offered as a 4-day summer workshop from 1999 to 2001 to 4th grade school teachers in the Seattle School District. The objectives of the workshop were to increase the teachers' knowledge of food science, to perform simple experiments that could be used in the 4th grade classroom, and to help the…
Control of the new 4th-order hyper-chaotic system with one input
NASA Astrophysics Data System (ADS)
Loría, Antonio
2010-06-01
We solve the problem of chaos suppression of Lü's hyper-chaotic system via feedback control. We use only one control input and moreover the controller is a simple proportional feedback and uses the measurement of only one variable. We show that this simple control law suffices to stabilize the hyper-chaotic system to the zero equilibrium globally and asymptotically. We present stability proofs based on Lyapunov's direct method and integration of solutions. As a corollary of our main result we draw the conclusion that the system is globally stabilizable by simply varying one parameter, when possible. Simulation experiments that show the effectiveness of our method are also presented.
NASA Astrophysics Data System (ADS)
Fang, Ming-chung; Lee, Zi-yi
2013-08-01
This paper develops a nonlinear mathematical model to simulate the dynamic motion behavior of the barge equipped with the portable outboard Dynamic Positioning (DP) system in short-crested waves. The self-tuning Proportional-Derivative (PD) controller based on the neural network algorithm is applied to control the thrusters for optimal adjustment of the barge position in waves. In addition to the wave, the current, the wind and the nonlinear drift force are also considered in the calculations. The time domain simulations for the six-degree-of-freedom motions of the barge with the DP system are solved by the 4th order Runge-Kutta method which can compromise the efficiency and the accuracy of the simulations. The technique of the portable alternative DP system developed here can serve as a practical tool to assist those ships without being equipped with the DP facility while the dynamic positioning missions are needed.
NASA Astrophysics Data System (ADS)
Tukaram Aghav, Sandip; Achyut Gangal, Shashikala
2014-06-01
In this paper, the main work is focused on designing and simplifying the orbit determination algorithm which will be used for Low Earth Orbit (LEO) navigation. The various data processing algorithms, state estimation algorithms and modeling forces were studied in detail, and simplified algorithm is selected to reduce hardware burden and computational cost. This is done by using raw navigation solution provided by GPS Navigation sensor. A fixed step-size Runge-Kutta 4th order numerical integration method is selected for orbit propagation. Both, the least square and Extended Kalman Filter (EKF) orbit estimation algorithms are developed and the results of the same are compared with each other. EKF algorithm converges faster than least square algorithm. EKF algorithm satisfies the criterions of low computation burden which is required for autonomous orbit determination. Simple static force models also feasible to reduce the hardware burden and computational cost.
Charged particle tracking at Titan, and further applications
NASA Astrophysics Data System (ADS)
Bebesi, Zsofia; Erdos, Geza; Szego, Karoly
2016-04-01
We use the CAPS ion data of Cassini to investigate the dynamics and origin of Titan's atmospheric ions. We developed a 4th order Runge-Kutta method to calculate particle trajectories in a time reversed scenario. The test particle magnetic field environment imitates the curved magnetic environment in the vicinity of Titan. The minimum variance directions along the S/C trajectory have been calculated for all available Titan flybys, and we assumed a homogeneous field that is perpendicular to the minimum variance direction. Using this method the magnetic field lines have been calculated along the flyby orbits so we could select those observational intervals when Cassini and the upper atmosphere of Titan were magnetically connected. We have also taken the Kronian magnetodisc into consideration, and used different upstream magnetic field approximations depending on whether Titan was located inside of the magnetodisc current sheet, or in the lobe regions. We also discuss the code's applicability to comets.
NASA Astrophysics Data System (ADS)
Núñez-de la Rosa, Jonatan; Munz, Claus-Dieter
2016-04-01
In this work we discuss the extension of the XTROEM-FV code to relativistic hydrodynamics and magnetohydrodynamics. XTROEM-FV is a simulation package for computational astrophysics based on very high order finite-volume methods on Cartesian coordinates. Arbitrary spatial high order of accuracy is achieved with a WENO reconstruction operator, and the time evolution is carried out with a strong-stability preserving Runge-Kutta scheme. In XTROEM-FV has been implemented a cheap, robust, and accurate shock capturing strategy for handling complex shock waves problems, typical in an astrophysical environment. The divergence constraint of the magnetic field is tackled with the generalized Lagrange multiplier divergence cleaning approach. Numerical computations of smooth flows for the relativistic hydrodynamics and magnetohydrodynamics equations are performed and confirm the high order accuracy of the main reconstruction algorithm for such kind of flows. XTROEM-FV has been subject to a comprehensive numerical benchmark, especially for complex flows configurations within an astrophysical context. Computations of problems with shocks with very high order reconstruction operators up to seventh order are reported. For instance, one-dimensional shock tubes problems for relativistic hydrodynamics and magnetohydrodynamics, as well as two-dimensional flows like the relativistic double Mach reflection problem, the interaction of a shock wave with a bubble, the relativistic Orszag-Tang vortex, the cylindrical blast wave problem, the rotor problem, the Kelvin-Helmholtz instability, and an astrophysical slab jet. XTROEM-FV represents a new attempt to simulate astrophysical flow phenomena with very high order numerical methods.
NASA Astrophysics Data System (ADS)
Núñez-de la Rosa, Jonatan; Munz, Claus-Dieter
2016-07-01
In this work, we discuss the extension of the XTROEM-FV code to relativistic hydrodynamics and magnetohydrodynamics. XTROEM-FV is a simulation package for computational astrophysics based on very high order finite-volume methods on Cartesian coordinates. Arbitrary spatial high order of accuracy is achieved with a weighted essentially non-oscillatory (WENO) reconstruction operator, and the time evolution is carried out with a strong stability preserving Runge-Kutta scheme. In XTROEM-FV has been implemented a cheap, robust, and accurate shock-capturing strategy for handling complex shock waves problems, typical in an astrophysical environment. The divergence constraint of the magnetic field is tackled with the generalized Lagrange multiplier divergence cleaning approach. Numerical computations of smooth flows for the relativistic hydrodynamics and magnetohydrodynamics equations are performed and confirm the high-order accuracy of the main reconstruction algorithm for such kind of flows. XTROEM-FV has been subject to a comprehensive numerical benchmark, especially for complex flows configurations within an astrophysical context. Computations of problems with shocks with very high order reconstruction operators up to seventh order are reported. For instance, one-dimensional shock tubes problems for relativistic hydrodynamics and magnetohydrodynamics, as well as two-dimensional flows like the relativistic double Mach reflection problem, the interaction of a shock wave with a bubble, the relativistic Orszag-Tang vortex, the cylindrical blast wave problem, the rotor problem, the Kelvin-Helmholtz instability, and an astrophysical slab jet. XTROEM-FV represents a new attempt to simulate astrophysical flow phenomena with very high order numerical methods.
29. VIEW OF 4TH FLOOR'S TELEPHONE RACKS WITH CABLE TRAYS ...
29. VIEW OF 4TH FLOOR'S TELEPHONE RACKS WITH CABLE TRAYS ABOVE. THESE ARE NEWER APPARATUS AND NOT ORIGINAL. - Pacific Telephone & Telegraph Company Building, 1519 Franklin Street, Oakland, Alameda County, CA
16. 4th floor roof, view west, north side of setback ...
16. 4th floor roof, view west, north side of setback to left and delivery stair bulkhead to right - Sheffield Farms Milk Plant, 1075 Webster Avenue (southwest corner of 166th Street), Bronx, Bronx County, NY
TID Test Results for 4th Generation iPad(TradeMark)
NASA Technical Reports Server (NTRS)
Guertin, S. M.; Allen, G. R.; McClure, S. S.; LaBel, K. A.
2013-01-01
TID testing of 4th generation iPads is reported. Of iPad subsystems, results indicate that the charging circuitry and display drivers fail at lowest TID levels. Details of construction are investigated for additional testing of components.
18. DETAILED OFFSHORE VIEW OF 4TH TEE, LOOKING NORTHWEST, SHOWING ...
18. DETAILED OFFSHORE VIEW OF 4TH TEE, LOOKING NORTHWEST, SHOWING TRANSITION FROM WOOD BENTS TO CONCRETE BENTS - Huntington Beach Municipal Pier, Pacific Coast Highway at Main Street, Huntington Beach, Orange County, CA
Physics Computing '92: Proceedings of the 4th International Conference
NASA Astrophysics Data System (ADS)
de Groot, Robert A.; Nadrchal, Jaroslav
1993-04-01
* Ordered Particle Simulations for Serial and MIMD Parallel Computers * "NOLP" -- Program Package for Laser Plasma Nonlinear Optics * Algorithms to Solve Nonlinear Least Square Problems * Distribution of Hydrogen Atoms in Pd-H Computed by Molecular Dynamics * A Ray Tracing of Optical System for Protein Crystallography Beamline at Storage Ring-SIBERIA-2 * Vibrational Properties of a Pseudobinary Linear Chain with Correlated Substitutional Disorder * Application of the Software Package Mathematica in Generalized Master Equation Method * Linelist: An Interactive Program for Analysing Beam-foil Spectra * GROMACS: A Parallel Computer for Molecular Dynamics Simulations * GROMACS Method of Virial Calculation Using a Single Sum * The Interactive Program for the Solution of the Laplace Equation with the Elimination of Singularities for Boundary Functions * Random-Number Generators: Testing Procedures and Comparison of RNG Algorithms * Micro-TOPIC: A Tokamak Plasma Impurities Code * Rotational Molecular Scattering Calculations * Orthonormal Polynomial Method for Calibrating of Cryogenic Temperature Sensors * Frame-based System Representing Basis of Physics * The Role of Massively Data-parallel Computers in Large Scale Molecular Dynamics Simulations * Short-range Molecular Dynamics on a Network of Processors and Workstations * An Algorithm for Higher-order Perturbation Theory in Radiative Transfer Computations * Hydrostochastics: The Master Equation Formulation of Fluid Dynamics * HPP Lattice Gas on Transputers and Networked Workstations * Study on the Hysteresis Cycle Simulation Using Modeling with Different Functions on Intervals * Refined Pruning Techniques for Feed-forward Neural Networks * Random Walk Simulation of the Motion of Transient Charges in Photoconductors * The Optical Hysteresis in Hydrogenated Amorphous Silicon * Diffusion Monte Carlo Analysis of Modern Interatomic Potentials for He * A Parallel Strategy for Molecular Dynamics Simulations of Polar Liquids on
European Code against Cancer 4th Edition: Infections and Cancer.
Villain, Patricia; Gonzalez, Paula; Almonte, Maribel; Franceschi, Silvia; Dillner, Joakim; Anttila, Ahti; Park, Jin Young; De Vuyst, Hugo; Herrero, Rolando
2015-12-01
Of the 2,635,000 new cancer cases (excluding non-melanoma skin cancers) occurring in the European Union (EU) in 2012, it is estimated that approximately 185,000 are related to infection with human papillomaviruses (HPVs), hepatitis B and C viruses (HBV and HCV), and Helicobacter pylori (H. pylori). Chronic infection with these agents can lead to cancers of the cervix uteri, liver, and stomach, respectively. Chronic infection with HCV can also lead to B-cell non-Hodgkin lymphoma. Human immunodeficiency virus (HIV) infection continues to be of major public health importance in several EU countries and increases cancer risk via HIV-induced immunosuppression. The fourth edition of the European Code Against Cancer presents recommendations on effective and safe preventive interventions in order to reduce the risk of infection-related cancers in EU citizens. Based on current available evidence, the fourth edition recommends that parents ensure the participation of their children in vaccination programs against HBV (for newborns) and HPV (for girls). In the 'Questions and Answers' (Q&As) section about vaccination and infections in the website for the European Code Against Cancer, individuals who are at risk of chronic HBV or HCV are advised to seek medical advice about testing and obtaining treatment when appropriate. Individuals most at risk of HIV are advised to consult their doctor or healthcare provider to access counselling and, if needed, testing and treatment without delay. Information about H. pylori testing and treatment is also provided as testing might currently be offered in some high-risk areas in Europe. The rationale and supporting evidence for the recommendations on vaccination in the European Code Against Cancer, and for the main recommendations on vaccination and infection in the Q&As, are explained in the present review. PMID:26589774
Physics Computing '92: Proceedings of the 4th International Conference
NASA Astrophysics Data System (ADS)
de Groot, Robert A.; Nadrchal, Jaroslav
1993-04-01
* Ordered Particle Simulations for Serial and MIMD Parallel Computers * "NOLP" -- Program Package for Laser Plasma Nonlinear Optics * Algorithms to Solve Nonlinear Least Square Problems * Distribution of Hydrogen Atoms in Pd-H Computed by Molecular Dynamics * A Ray Tracing of Optical System for Protein Crystallography Beamline at Storage Ring-SIBERIA-2 * Vibrational Properties of a Pseudobinary Linear Chain with Correlated Substitutional Disorder * Application of the Software Package Mathematica in Generalized Master Equation Method * Linelist: An Interactive Program for Analysing Beam-foil Spectra * GROMACS: A Parallel Computer for Molecular Dynamics Simulations * GROMACS Method of Virial Calculation Using a Single Sum * The Interactive Program for the Solution of the Laplace Equation with the Elimination of Singularities for Boundary Functions * Random-Number Generators: Testing Procedures and Comparison of RNG Algorithms * Micro-TOPIC: A Tokamak Plasma Impurities Code * Rotational Molecular Scattering Calculations * Orthonormal Polynomial Method for Calibrating of Cryogenic Temperature Sensors * Frame-based System Representing Basis of Physics * The Role of Massively Data-parallel Computers in Large Scale Molecular Dynamics Simulations * Short-range Molecular Dynamics on a Network of Processors and Workstations * An Algorithm for Higher-order Perturbation Theory in Radiative Transfer Computations * Hydrostochastics: The Master Equation Formulation of Fluid Dynamics * HPP Lattice Gas on Transputers and Networked Workstations * Study on the Hysteresis Cycle Simulation Using Modeling with Different Functions on Intervals * Refined Pruning Techniques for Feed-forward Neural Networks * Random Walk Simulation of the Motion of Transient Charges in Photoconductors * The Optical Hysteresis in Hydrogenated Amorphous Silicon * Diffusion Monte Carlo Analysis of Modern Interatomic Potentials for He * A Parallel Strategy for Molecular Dynamics Simulations of Polar Liquids on
Cutting orientations for non-complex parts in 4th axis machining
NASA Astrophysics Data System (ADS)
Osman Zahid, M. N.; Case, K.; Watts, D. M.
2016-02-01
The application of Computer Numerically Controlled (CNC) machining for Rapid Manufacturing processes (CNC-RM) exploits the innate potential of 4th axis machining. The use of an indexer allows the workpiece to be rotated to various orientations which directly increased the region accessible to the cutting tool. However, in order to avoid thin webs and preserve tool life, cutting must be executed with a minimum of three orientations even for geometrically simple parts. Recent findings have suggested the separation of cutting orientations into roughing and finishing operations. Thus, the selection of orientations in finishing processes becomes more flexible and independent. This study was conducted to identify the effects of using a minimum of two cutting orientations in finishing operations for CNC-RM applications. This method is only applicable for non-complex parts where all the features can be machined from two directions. The results of the study illustrate the positive effects of minimizing the number of orientations. Despite improvement in machining operations, the complexity in defining the cutting orientations was also reduced.
European Code against Cancer 4th Edition: Ionising and non-ionising radiation and cancer.
McColl, Neil; Auvinen, Anssi; Kesminiene, Ausrele; Espina, Carolina; Erdmann, Friederike; de Vries, Esther; Greinert, Rüdiger; Harrison, John; Schüz, Joachim
2015-12-01
Ionising radiation can transfer sufficient energy to ionise molecules, and this can lead to chemical changes, including DNA damage in cells. Key evidence for the carcinogenicity of ionising radiation comes from: follow-up studies of the survivors of the atomic bombings in Japan; other epidemiological studies of groups that have been exposed to radiation from medical, occupational or environmental sources; experimental animal studies; and studies of cellular responses to radiation. Considering exposure to environmental ionising radiation, inhalation of naturally occurring radon is the major source of radiation in the population - in doses orders of magnitude higher than those from nuclear power production or nuclear fallout. Indoor exposure to radon and its decay products is an important cause of lung cancer; radon may cause approximately one in ten lung cancers in Europe. Exposures to radon in buildings can be reduced via a three-step process of identifying those with potentially elevated radon levels, measuring radon levels, and reducing exposure by installation of remediation systems. In the 4th Edition of the European Code against Cancer it is therefore recommended to: "Find out if you are exposed to radiation from naturally high radon levels in your home. Take action to reduce high radon levels". Non-ionising types of radiation (those with insufficient energy to ionise molecules) - including extremely low-frequency electric and magnetic fields as well as radiofrequency electromagnetic fields - are not an established cause of cancer and are therefore not addressed in the recommendations to reduce cancer risk. PMID:26126928
Seligmann, Hervé
2016-01-01
In mitochondria, secondary structures punctuate post-transcriptional RNA processing. Recently described transcripts match the human mitogenome after systematic deletions of every 4th, respectively every 4th and 5th nucleotides, called delRNAs. Here I explore predicted stem-loop hairpin formation by delRNAs, and their associations with delRNA transcription and detected peptides matching their translation. Despite missing 25, respectively 40% of the nucleotides in the original sequence, del-transformed sequences form significantly more secondary structures than corresponding randomly shuffled sequences, indicating biological function, independently of, and in combination with, previously detected delRNA and thereof translated peptides. Self-hybridization decreases delRNA abundances, indicating downregulation. Systematic deletions of the human mitogenome reveal new, unsuspected coding and structural informations. PMID:27018206
75 FR 38721 - Safety Zone; Munising 4th of July Fireworks, South Bay, Lake Superior, Munising, MI
Federal Register 2010, 2011, 2012, 2013, 2014
2010-07-06
... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zone; Munising 4th of July Fireworks, South Bay... is intended to restrict vessels from a portion of South Bay during the Munising 4th of July Fireworks... from hazards associated with the Munising 4th of July Fireworks display. Based on the explosive...
Federal Register 2010, 2011, 2012, 2013, 2014
2010-06-17
... Chicago's July 4th Celebration Fireworks, Chicago, Illinois in the Federal Register (75 FR 22330). We... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zones; City of Chicago's July 4th Celebration... associated with the City of Chicago's July 4th Celebration Fireworks. The Captain of the Port, Sector...
75 FR 34379 - Safety Zone; Mackinac Island 4th of July Fireworks, Lake Huron, Mackinac Island, MI
Federal Register 2010, 2011, 2012, 2013, 2014
2010-06-17
... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zone; Mackinac Island 4th of July Fireworks, Lake... intended to restrict vessels from a portion of Lake Huron during the Mackinac Island 4th of July Fireworks... with the Mackinac Island 4th of July fireworks display. The fireworks display will occur between 9:45...
Federal Register 2010, 2011, 2012, 2013, 2014
2010-07-06
... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zone; Sault Sainte Marie 4th of July Fireworks... the Sault Sainte Marie 4th of July Fireworks display, July ] 4, 2010. This temporary safety zone is... with the Sault Sainte Marie 4th of July Fireworks display. The fireworks display is planned to...
Hauck, Cory D; Alldredge, Graham; Tits, Andre
2012-01-01
We present a numerical algorithm to implement entropy-based (M{sub N}) moment models in the context of a simple, linear kinetic equation for particles moving through a material slab. The closure for these models - as is the case for all entropy-based models - is derived through the solution of constrained, convex optimization problem. The algorithm has two components. The first component is a discretization of the moment equations which preserves the set of realizable moments, thereby ensuring that the optimization problem has a solution (in exact arithmetic). The discretization is a second-order kinetic scheme which uses MUSCL-type limiting in space and a strong-stability-preserving, Runge-Kutta time integrator. The second component of the algorithm is a Newton-based solver for the dual optimization problem, which uses an adaptive quadrature to evaluate integrals in the dual objective and its derivatives. The accuracy of the numerical solution to the dual problem plays a key role in the time step restriction for the kinetic scheme. We study in detail the difficulties in the dual problem that arise near the boundary of realizable moments, where quadrature formulas are less reliable and the Hessian of the dual objection function is highly ill-conditioned. Extensive numerical experiments are performed to illustrate these difficulties. In cases where the dual problem becomes 'too difficult' to solve numerically, we propose a regularization technique to artificially move moments away from the realizable boundary in a way that still preserves local particle concentrations. We present results of numerical simulations for two challenging test problems in order to quantify the characteristics of the optimization solver and to investigate when and how frequently the regularization is needed.
Development of a variable time-step transient NEW code: SPANDEX
Aviles, B.N. )
1993-01-01
This paper describes a three-dimensional, variable time-step transient multigroup diffusion theory code, SPANDEX (space-time nodal expansion method). SPANDEX is based on the static nodal expansion method (NEM) code, NODEX (Ref. 1), and employs a nonlinear algorithm and a fifth-order expansion of the transverse-integrated fluxes. The time integration scheme in SPANDEX is a fourth-order implicit generalized Runge-Kutta method (GRK) with on-line error control and variable time-step selection. This Runge-Kutta method has been applied previously to point kinetics and one-dimensional finite difference transient analysis. This paper describes the application of the Runge-Kutta method to three-dimensional reactor transient analysis in a multigroup NEM code.
European Code against Cancer 4th Edition: Environment, occupation and cancer.
Espina, Carolina; Straif, Kurt; Friis, Søren; Kogevinas, Manolis; Saracci, Rodolfo; Vainio, Harri; Schüz, Joachim
2015-12-01
People are exposed throughout life to a wide range of environmental and occupational pollutants from different sources at home, in the workplace or in the general environment - exposures that normally cannot be directly controlled by the individual. Several chemicals, metals, dusts, fibres, and occupations have been established to be causally associated with an increased risk of specific cancers, such as cancers of the lung, skin and urinary bladder, and mesothelioma. Significant amounts of air pollutants - mainly from road transport and industry - continue to be emitted in the European Union (EU); an increased occurrence of lung cancer has been attributed to air pollution even in areas below the EU limits for daily air pollution. Additionally, a wide range of pesticides as well as industrial and household chemicals may lead to widespread human exposure, mainly through food and water. For most environmental pollutants, the most effective measures are regulations and community actions aimed at reducing and eliminating the exposures. Thus, it is imperative to raise awareness about environmental and occupational carcinogens in order to motivate individuals to be proactive in advocating protection and supporting initiatives aimed at reducing pollution. Regulations are not homogeneous across EU countries, and protective measures in the workplace are not used consistently by all workers all the time; compliance with regulations needs to be continuously monitored and enforced. Therefore, the recommendation on Environment and Occupation of the 4th edition of the European Code against Cancer, focusing on what individuals can do to reduce their cancer risk, reads: "In the workplace, protect yourself against cancer-causing substances by following health and safety instructions." PMID:26164655
20. TYPICAL VIEW OF FRONT WINDOWS FROM 4TH TO 9TH ...
20. TYPICAL VIEW OF FRONT WINDOWS FROM 4TH TO 9TH FLOOR WITH WHITE GLAZED TERRA COTTA SILL AND HEADERS. MULLIONS ARE ORANGE BROWN BRICKS LIKE THE WALLS. BRICKS ARE IN FLEMISH BOND PATTERN. - Pacific Telephone & Telegraph Company Building, 1519 Franklin Street, Oakland, Alameda County, CA
77 FR 39408 - Safety Zone; Buffalo July 4th Fireworks, Lake Erie, Buffalo, NY
Federal Register 2010, 2011, 2012, 2013, 2014
2012-07-03
... DHS Department of Homeland Security FR Federal Register NPRM Notice of Proposed Rulemaking A... SECURITY Coast Guard 33 CFR Part 165 RIN 1625-AA00 Safety Zone; Buffalo July 4th Fireworks, Lake Erie, Buffalo, NY AGENCY: Coast Guard, DHS. ACTION: Temporary final rule. SUMMARY: The Coast Guard...