Time-symmetric integration in astrophysics

Science.gov (United States)

Hernandez, David M.; Bertschinger, Edmund

2018-04-01

Calculating the long-term solution of ordinary differential equations, such as those of the N-body problem, is central to understanding a wide range of dynamics in astrophysics, from galaxy formation to planetary chaos. Because generally no analytic solution exists to these equations, researchers rely on numerical methods that are prone to various errors. In an effort to mitigate these errors, powerful symplectic integrators have been employed. But symplectic integrators can be severely limited because they are not compatible with adaptive stepping and thus they have difficulty in accommodating changing time and length scales. A promising alternative is time-reversible integration, which can handle adaptive time-stepping, but the errors due to time-reversible integration in astrophysics are less understood. The goal of this work is to study analytically and numerically the errors caused by time-reversible integration, with and without adaptive stepping. We derive the modified differential equations of these integrators to perform the error analysis. As an example, we consider the trapezoidal rule, a reversible non-symplectic integrator, and show that it gives secular energy error increase for a pendulum problem and for a Hénon-Heiles orbit. We conclude that using reversible integration does not guarantee good energy conservation and that, when possible, use of symplectic integrators is favoured. We also show that time-symmetry and time-reversibility are properties that are distinct for an integrator.

Symplectic integrators with adaptive time steps

Science.gov (United States)

Richardson, A. S.; Finn, J. M.

2012-01-01

In recent decades, there have been many attempts to construct symplectic integrators with variable time steps, with rather disappointing results. In this paper, we identify the causes for this lack of performance, and find that they fall into two categories. In the first, the time step is considered a function of time alone, Δ = Δ(t). In this case, backward error analysis shows that while the algorithms remain symplectic, parametric instabilities may arise because of resonance between oscillations of Δ(t) and the orbital motion. In the second category the time step is a function of phase space variables Δ = Δ(q, p). In this case, the system of equations to be solved is analyzed by introducing a new time variable τ with dt = Δ(q, p) dτ. The transformed equations are no longer in Hamiltonian form, and thus do not benefit from integration methods which would be symplectic for Hamiltonian systems. We analyze two methods for integrating the transformed equations which do, however, preserve the structure of the original equations. The first is an extended phase space method, which has been successfully used in previous studies of adaptive time step symplectic integrators. The second, novel, method is based on a non-canonical mixed-variable generating function. Numerical trials for both of these methods show good results, without parametric instabilities or spurious growth or damping. It is then shown how to adapt the time step to an error estimate found by backward error analysis, in order to optimize the time-stepping scheme. Numerical results are obtained using this formulation and compared with other time-stepping schemes for the extended phase space symplectic method.

Exponential integrators in time-dependent density-functional calculations

Science.gov (United States)

Kidd, Daniel; Covington, Cody; Varga, Kálmán

2017-12-01

The integrating factor and exponential time differencing methods are implemented and tested for solving the time-dependent Kohn-Sham equations. Popular time propagation methods used in physics, as well as other robust numerical approaches, are compared to these exponential integrator methods in order to judge the relative merit of the computational schemes. We determine an improvement in accuracy of multiple orders of magnitude when describing dynamics driven primarily by a nonlinear potential. For cases of dynamics driven by a time-dependent external potential, the accuracy of the exponential integrator methods are less enhanced but still match or outperform the best of the conventional methods tested.

EVALUATION OF THE POUNDING FORCES DURING EARTHQUAKE USING EXPLICIT DYNAMIC TIME INTEGRATION METHOD

Directory of Open Access Journals (Sweden)

Nica George Bogdan

2017-09-01

Full Text Available Pounding effects during earthquake is a subject of high significance for structural engineers performing in the urban areas. In this paper, two ways to account for structural pounding are used in a MATLAB code, namely classical stereomechanics approach and nonlinear viscoelastic impact element. The numerical study is performed on SDOF structures acted by ELCentro recording. While most of the studies available in the literature are related to Newmark implicit time integration method, in this study the equations of motion are numerical integrated using central finite difference method, an explicit method, having the main advantage that in the displacement at the ith+1 step is calculated based on the loads from the ith step. Thus, the collision is checked and the pounding forces are taken into account into the equation of motion in an easier manner than in an implicit integration method. First, a comparison is done using available data in the literature. Both linear and nonlinear behavior of the structures during earthquake is further investigated. Several layout scenarios are also investigated, in which one or more weak buildings are adjacent to a stiffer building. One of the main findings in this paper is related to the behavior of a weak structure located between two stiff structures.

A hybrid method combining the FDTD and a time domain boundary-integral equation marching-on-in-time algorithm

Directory of Open Access Journals (Sweden)

A. Becker

2003-01-01

Full Text Available In this paper a hybrid method combining the FDTD/FIT with a Time Domain Boundary-Integral Marching-on-in-Time Algorithm (TD-BIM is presented. Inhomogeneous regions are modelled with the FIT-method, an alternative formulation of the FDTD. Homogeneous regions (which is in the presented numerical example the open space are modelled using a TD-BIM with equivalent electric and magnetic currents flowing on the boundary between the inhomogeneous and the homogeneous regions. The regions are coupled by the tangential magnetic fields just outside the inhomogeneous regions. These fields are calculated by making use of a Mixed Potential Integral Formulation for the magnetic field. The latter consists of equivalent electric and magnetic currents on the boundary plane between the homogeneous and the inhomogeneous region. The magnetic currents result directly from the electric fields of the Yee lattice. Electric currents in the same plane are calculated by making use of the TD-BIM and using the electric field of the Yee lattice as boundary condition. The presented hybrid method only needs the interpolations inherent in FIT and no additional interpolation. A numerical result is compared to a calculation that models both regions with FDTD.

A discontinuous galerkin time domain-boundary integral method for analyzing transient electromagnetic scattering

KAUST Repository

Li, Ping

2014-07-01

This paper presents an algorithm hybridizing discontinuous Galerkin time domain (DGTD) method and time domain boundary integral (BI) algorithm for 3-D open region electromagnetic scattering analysis. The computational domain of DGTD is rigorously truncated by analytically evaluating the incoming numerical flux from the outside of the truncation boundary through BI method based on the Huygens\\' principle. The advantages of the proposed method are that it allows the truncation boundary to be conformal to arbitrary (convex/ concave) scattering objects, well-separated scatters can be truncated by their local meshes without losing the physics (such as coupling/multiple scattering) of the problem, thus reducing the total mesh elements. Furthermore, low frequency waves can be efficiently absorbed, and the field outside the truncation domain can be conveniently calculated using the same BI formulation. Numerical examples are benchmarked to demonstrate the accuracy and versatility of the proposed method.

A higher order space-time Galerkin scheme for time domain integral equations

KAUST Repository

Pray, Andrew J.

2014-12-01

Stability of time domain integral equation (TDIE) solvers has remained an elusive goal formany years. Advancement of this research has largely progressed on four fronts: 1) Exact integration, 2) Lubich quadrature, 3) smooth temporal basis functions, and 4) space-time separation of convolutions with the retarded potential. The latter method\\'s efficacy in stabilizing solutions to the time domain electric field integral equation (TD-EFIE) was previously reported for first-order surface descriptions (flat elements) and zeroth-order functions as the temporal basis. In this work, we develop the methodology necessary to extend the scheme to higher order surface descriptions as well as to enable its use with higher order basis functions in both space and time. These basis functions are then used in a space-time Galerkin framework. A number of results are presented that demonstrate convergence in time. The viability of the space-time separation method in producing stable results is demonstrated experimentally for these examples.

Numerical Time Integration Methods for a Point Absorber Wave Energy Converter

DEFF Research Database (Denmark)

Zurkinden, Andrew Stephen; Kramer, Morten

2012-01-01

on a discretization of the convolution integral. The calculation of the convolution integral is performed at each time step regardless of the chosen numerical scheme. In the second model the convolution integral is replaced by a system of linear ordinary differential equations. The formulation of the state...

Integrable Time-Dependent Quantum Hamiltonians

Science.gov (United States)

Sinitsyn, Nikolai A.; Yuzbashyan, Emil A.; Chernyak, Vladimir Y.; Patra, Aniket; Sun, Chen

2018-05-01

We formulate a set of conditions under which the nonstationary Schrödinger equation with a time-dependent Hamiltonian is exactly solvable analytically. The main requirement is the existence of a non-Abelian gauge field with zero curvature in the space of system parameters. Known solvable multistate Landau-Zener models satisfy these conditions. Our method provides a strategy to incorporate time dependence into various quantum integrable models while maintaining their integrability. We also validate some prior conjectures, including the solution of the driven generalized Tavis-Cummings model.

Expansion methods for solving integral equations with multiple time lags using Bernstein polynomial of the second kind

Directory of Open Access Journals (Sweden)

Mahmoud Paripour

2014-08-01

Full Text Available In this paper, the Bernstein polynomials are used to approximatethe solutions of linear integral equations with multiple time lags (IEMTL through expansion methods (collocation method, partition method, Galerkin method. The method is discussed in detail and illustrated by solving some numerical examples. Comparison between the exact and approximated results obtained from these methods is carried out

A hybrid time-domain discontinuous galerkin-boundary integral method for electromagnetic scattering analysis

KAUST Repository

Li, Ping; Shi, Yifei; Jiang, Lijun; Bagci, Hakan

2014-01-01

A scheme hybridizing discontinuous Galerkin time-domain (DGTD) and time-domain boundary integral (TDBI) methods for accurately analyzing transient electromagnetic scattering is proposed. Radiation condition is enforced using the numerical flux on the truncation boundary. The fields required by the flux are computed using the TDBI from equivalent currents introduced on a Huygens' surface enclosing the scatterer. The hybrid DGTDBI ensures that the radiation condition is mathematically exact and the resulting computation domain is as small as possible since the truncation boundary conforms to scatterer's shape and is located very close to its surface. Locally truncated domains can also be defined around each disconnected scatterer additionally reducing the size of the overall computation domain. Numerical examples demonstrating the accuracy and versatility of the proposed method are presented. © 2014 IEEE.

A hybrid time-domain discontinuous galerkin-boundary integral method for electromagnetic scattering analysis

KAUST Repository

Li, Ping

2014-05-01

A scheme hybridizing discontinuous Galerkin time-domain (DGTD) and time-domain boundary integral (TDBI) methods for accurately analyzing transient electromagnetic scattering is proposed. Radiation condition is enforced using the numerical flux on the truncation boundary. The fields required by the flux are computed using the TDBI from equivalent currents introduced on a Huygens\\' surface enclosing the scatterer. The hybrid DGTDBI ensures that the radiation condition is mathematically exact and the resulting computation domain is as small as possible since the truncation boundary conforms to scatterer\\'s shape and is located very close to its surface. Locally truncated domains can also be defined around each disconnected scatterer additionally reducing the size of the overall computation domain. Numerical examples demonstrating the accuracy and versatility of the proposed method are presented. © 2014 IEEE.

A higher order space-time Galerkin scheme for time domain integral equations

KAUST Repository

Pray, Andrew J.; Beghein, Yves; Nair, Naveen V.; Cools, Kristof; Bagci, Hakan; Shanker, Balasubramaniam

2014-01-01

Stability of time domain integral equation (TDIE) solvers has remained an elusive goal formany years. Advancement of this research has largely progressed on four fronts: 1) Exact integration, 2) Lubich quadrature, 3) smooth temporal basis functions, and 4) space-time separation of convolutions with the retarded potential. The latter method's efficacy in stabilizing solutions to the time domain electric field integral equation (TD-EFIE) was previously reported for first-order surface descriptions (flat elements) and zeroth-order functions as the temporal basis. In this work, we develop the methodology necessary to extend the scheme to higher order surface descriptions as well as to enable its use with higher order basis functions in both space and time. These basis functions are then used in a space-time Galerkin framework. A number of results are presented that demonstrate convergence in time. The viability of the space-time separation method in producing stable results is demonstrated experimentally for these examples.

Doubling time measurement by method of digital integration of pulses from a CFU7 detector

International Nuclear Information System (INIS)

Gauthier, Guy

1968-01-01

The author reports an experimental study which aimed at measuring the doubling time with a CFU7 fission chamber by using an integration method on the Siloette pile with a new core and with a spent core, and at comparing results with those obtained with a specific instrument which receives information from an ionisation chamber. The interest of this method relies on the fact that the fission chamber is insensitive to gamma radiations

Time-splitting combined with exponential wave integrator fourier pseudospectral method for Schrödinger-Boussinesq system

Science.gov (United States)

Liao, Feng; Zhang, Luming; Wang, Shanshan

2018-02-01

In this article, we formulate an efficient and accurate numerical method for approximations of the coupled Schrödinger-Boussinesq (SBq) system. The main features of our method are based on: (i) the applications of a time-splitting Fourier spectral method for Schrödinger-like equation in SBq system, (ii) the utilizations of exponential wave integrator Fourier pseudospectral for spatial derivatives in the Boussinesq-like equation. The scheme is fully explicit and efficient due to fast Fourier transform. The numerical examples are presented to show the efficiency and accuracy of our method.

Practical method of calculating time-integrated concentrations at medium and large distances

International Nuclear Information System (INIS)

Cagnetti, P.; Ferrara, V.

1980-01-01

Previous reports have covered the possibility of calculating time-integrated concentrations (TICs) for a prolonged release, based on concentration estimates for a brief release. This study proposes a simple method of evaluating concentrations in the air at medium and large distances, for a brief release. It is known that the stability of the atmospheric layers close to ground level influence diffusion only over short distances. Beyond some tens of kilometers, as the pollutant cloud progressively reaches higher layers, diffusion is affected by factors other than the stability at ground level, such as wind shear for intermediate distances and the divergence and rotational motion of air masses towards the upper limit of the mesoscale and on the synoptic scale. Using the data available in the literature, expressions for sigmasub(y) and sigmasub(z) are proposed for transfer times corresponding to those for up to distances of several thousand kilometres, for two initial diffusion situations (up to distances of 10 - 20 km), those characterized by stable and neutral conditions respectively. Using this method simple hand calculations can be made for any problem relating to the diffusion of radioactive pollutants over long distances

Efficient orbit integration by manifold correction methods.

Science.gov (United States)

Fukushima, Toshio

2005-12-01

Triggered by a desire to investigate, numerically, the planetary precession through a long-term numerical integration of the solar system, we developed a new formulation of numerical integration of orbital motion named manifold correct on methods. The main trick is to rigorously retain the consistency of physical relations, such as the orbital energy, the orbital angular momentum, or the Laplace integral, of a binary subsystem. This maintenance is done by applying a correction to the integrated variables at each integration step. Typical methods of correction are certain geometric transformations, such as spatial scaling and spatial rotation, which are commonly used in the comparison of reference frames, or mathematically reasonable operations, such as modularization of angle variables into the standard domain [-pi, pi). The form of the manifold correction methods finally evolved are the orbital longitude methods, which enable us to conduct an extremely precise integration of orbital motions. In unperturbed orbits, the integration errors are suppressed at the machine epsilon level for an indefinitely long period. In perturbed cases, on the other hand, the errors initially grow in proportion to the square root of time and then increase more rapidly, the onset of which depends on the type and magnitude of the perturbations. This feature is also realized for highly eccentric orbits by applying the same idea as used in KS-regularization. In particular, the introduction of time elements greatly enhances the performance of numerical integration of KS-regularized orbits, whether the scaling is applied or not.

A variable timestep generalized Runge-Kutta method for the numerical integration of the space-time diffusion equations

International Nuclear Information System (INIS)

Aviles, B.N.; Sutton, T.M.; Kelly, D.J. III.

1991-09-01

A generalized Runge-Kutta method has been employed in the numerical integration of the stiff space-time diffusion equations. The method is fourth-order accurate, using an embedded third-order solution to arrive at an estimate of the truncation error for automatic timestep control. The efficiency of the Runge-Kutta method is enhanced by a block-factorization technique that exploits the sparse structure of the matrix system resulting from the space and energy discretized form of the time-dependent neutron diffusion equations. Preliminary numerical evaluation using a one-dimensional finite difference code shows the sparse matrix implementation of the generalized Runge-Kutta method to be highly accurate and efficient when compared to an optimized iterative theta method. 12 refs., 5 figs., 4 tabs

Numerov iteration method for second order integral-differential equation

International Nuclear Information System (INIS)

Zeng Fanan; Zhang Jiaju; Zhao Xuan

1987-01-01

In this paper, Numerov iterative method for second order integral-differential equation and system of equations are constructed. Numerical examples show that this method is better than direct method (Gauss elimination method) in CPU time and memoy requireing. Therefore, this method is an efficient method for solving integral-differential equation in nuclear physics

A point implicit time integration technique for slow transient flow problems

Energy Technology Data Exchange (ETDEWEB)

Kadioglu, Samet Y., E-mail: kadioglu@yildiz.edu.tr [Department of Mathematical Engineering, Yildiz Technical University, 34210 Davutpasa-Esenler, Istanbul (Turkey); Berry, Ray A., E-mail: ray.berry@inl.gov [Idaho National Laboratory, P.O. Box 1625, MS 3840, Idaho Falls, ID 83415 (United States); Martineau, Richard C. [Idaho National Laboratory, P.O. Box 1625, MS 3840, Idaho Falls, ID 83415 (United States)

2015-05-15

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

A point implicit time integration technique for slow transient flow problems

International Nuclear Information System (INIS)

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

2015-01-01

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

The Use of Conditional Probability Integral Transformation Method for Testing Accelerated Failure Time Models

Directory of Open Access Journals (Sweden)

Abdalla Ahmed Abdel-Ghaly

2016-06-01

Full Text Available This paper suggests the use of the conditional probability integral transformation (CPIT method as a goodness of fit (GOF technique in the field of accelerated life testing (ALT, specifically for validating the underlying distributional assumption in accelerated failure time (AFT model. The method is based on transforming the data into independent and identically distributed (i.i.d Uniform (0, 1 random variables and then applying the modified Watson statistic to test the uniformity of the transformed random variables. This technique is used to validate each of the exponential, Weibull and lognormal distributions' assumptions in AFT model under constant stress and complete sampling. The performance of the CPIT method is investigated via a simulation study. It is concluded that this method performs well in case of exponential and lognormal distributions. Finally, a real life example is provided to illustrate the application of the proposed procedure.

Time integration of tensor trains

OpenAIRE

Lubich, Christian; Oseledets, Ivan; Vandereycken, Bart

2014-01-01

A robust and efficient time integrator for dynamical tensor approximation in the tensor train or matrix product state format is presented. The method is based on splitting the projector onto the tangent space of the tensor manifold. The algorithm can be used for updating time-dependent tensors in the given data-sparse tensor train / matrix product state format and for computing an approximate solution to high-dimensional tensor differential equations within this data-sparse format. The formul...

An analytical method for PID controller tuning with specified gain and phase margins for integral plus time delay processes.

Science.gov (United States)

Hu, Wuhua; Xiao, Gaoxi; Li, Xiumin

2011-04-01

In this paper, an analytical method is proposed for proportional-integral/proportional-derivative/proportional-integral-derivative (PI/PD/PID) controller tuning with specified gain and phase margins (GPMs) for integral plus time delay (IPTD) processes. Explicit formulas are also obtained for estimating the GPMs resulting from given PI/PD/PID controllers. The proposed method indicates a general form of the PID parameters and unifies a large number of existing rules as PI/PD/PID controller tuning with various GPM specifications. The GPMs realized by existing PID tuning rules are computed and documented as a reference for control engineers to tune the PID controllers. Copyright © 2011 ISA. Published by Elsevier Ltd. All rights reserved.

Path integral solution for some time-dependent potential

International Nuclear Information System (INIS)

Storchak, S.N.

1989-12-01

The quantum-mechanical problem with a time-dependent potential is solved by the path integral method. The solution is obtained by the application of the previously derived general formula for rheonomic homogeneous point transformation and reparametrization in the path integral. (author). 4 refs

Rigorous time slicing approach to Feynman path integrals

CERN Document Server

Fujiwara, Daisuke

2017-01-01

This book proves that Feynman's original definition of the path integral actually converges to the fundamental solution of the Schrödinger equation at least in the short term if the potential is differentiable sufficiently many times and its derivatives of order equal to or higher than two are bounded. The semi-classical asymptotic formula up to the second term of the fundamental solution is also proved by a method different from that of Birkhoff. A bound of the remainder term is also proved. The Feynman path integral is a method of quantization using the Lagrangian function, whereas Schrödinger's quantization uses the Hamiltonian function. These two methods are believed to be equivalent. But equivalence is not fully proved mathematically, because, compared with Schrödinger's method, there is still much to be done concerning rigorous mathematical treatment of Feynman's method. Feynman himself defined a path integral as the limit of a sequence of integrals over finite-dimensional spaces which is obtained by...

Multiple time scale methods in tokamak magnetohydrodynamics

International Nuclear Information System (INIS)

Jardin, S.C.

1984-01-01

Several methods are discussed for integrating the magnetohydrodynamic (MHD) equations in tokamak systems on other than the fastest time scale. The dynamical grid method for simulating ideal MHD instabilities utilizes a natural nonorthogonal time-dependent coordinate transformation based on the magnetic field lines. The coordinate transformation is chosen to be free of the fast time scale motion itself, and to yield a relatively simple scalar equation for the total pressure, P = p + B 2 /2μ 0 , which can be integrated implicitly to average over the fast time scale oscillations. Two methods are described for the resistive time scale. The zero-mass method uses a reduced set of two-fluid transport equations obtained by expanding in the inverse magnetic Reynolds number, and in the small ratio of perpendicular to parallel mobilities and thermal conductivities. The momentum equation becomes a constraint equation that forces the pressure and magnetic fields and currents to remain in force balance equilibrium as they evolve. The large mass method artificially scales up the ion mass and viscosity, thereby reducing the severe time scale disparity between wavelike and diffusionlike phenomena, but not changing the resistive time scale behavior. Other methods addressing the intermediate time scales are discussed

A New time Integration Scheme for Cahn-hilliard Equations

KAUST Repository

Schaefer, R.

2015-06-01

In this paper we present a new integration scheme that can be applied to solving difficult non-stationary non-linear problems. It is obtained by a successive linearization of the Crank- Nicolson scheme, that is unconditionally stable, but requires solving non-linear equation at each time step. We applied our linearized scheme for the time integration of the challenging Cahn-Hilliard equation, modeling the phase separation in fluids. At each time step the resulting variational equation is solved using higher-order isogeometric finite element method, with B- spline basis functions. The method was implemented in the PETIGA framework interfaced via the PETSc toolkit. The GMRES iterative solver was utilized for the solution of a resulting linear system at every time step. We also apply a simple adaptivity rule, which increases the time step size when the number of GMRES iterations is lower than 30. We compared our method with a non-linear, two stage predictor-multicorrector scheme, utilizing a sophisticated step length adaptivity. We controlled the stability of our simulations by monitoring the Ginzburg-Landau free energy functional. The proposed integration scheme outperforms the two-stage competitor in terms of the execution time, at the same time having a similar evolution of the free energy functional.

A New time Integration Scheme for Cahn-hilliard Equations

KAUST Repository

Schaefer, R.; Smol-ka, M.; Dalcin, L; Paszyn'ski, M.

2015-01-01

In this paper we present a new integration scheme that can be applied to solving difficult non-stationary non-linear problems. It is obtained by a successive linearization of the Crank- Nicolson scheme, that is unconditionally stable, but requires solving non-linear equation at each time step. We applied our linearized scheme for the time integration of the challenging Cahn-Hilliard equation, modeling the phase separation in fluids. At each time step the resulting variational equation is solved using higher-order isogeometric finite element method, with B- spline basis functions. The method was implemented in the PETIGA framework interfaced via the PETSc toolkit. The GMRES iterative solver was utilized for the solution of a resulting linear system at every time step. We also apply a simple adaptivity rule, which increases the time step size when the number of GMRES iterations is lower than 30. We compared our method with a non-linear, two stage predictor-multicorrector scheme, utilizing a sophisticated step length adaptivity. We controlled the stability of our simulations by monitoring the Ginzburg-Landau free energy functional. The proposed integration scheme outperforms the two-stage competitor in terms of the execution time, at the same time having a similar evolution of the free energy functional.

Mixed time slicing in path integral simulations

International Nuclear Information System (INIS)

Steele, Ryan P.; Zwickl, Jill; Shushkov, Philip; Tully, John C.

2011-01-01

A simple and efficient scheme is presented for using different time slices for different degrees of freedom in path integral calculations. This method bridges the gap between full quantization and the standard mixed quantum-classical (MQC) scheme and, therefore, still provides quantum mechanical effects in the less-quantized variables. Underlying the algorithm is the notion that time slices (beads) may be 'collapsed' in a manner that preserves quantization in the less quantum mechanical degrees of freedom. The method is shown to be analogous to multiple-time step integration techniques in classical molecular dynamics. The algorithm and its associated error are demonstrated on model systems containing coupled high- and low-frequency modes; results indicate that convergence of quantum mechanical observables can be achieved with disparate bead numbers in the different modes. Cost estimates indicate that this procedure, much like the MQC method, is most efficient for only a relatively few quantum mechanical degrees of freedom, such as proton transfer. In this regime, however, the cost of a fully quantum mechanical simulation is determined by the quantization of the least quantum mechanical degrees of freedom.

Feasibility of real-time calculation of correlation integral derived statistics applied to EGG time series

NARCIS (Netherlands)

van den Broek, PLC; van Egmond, J; van Rijn, CM; Takens, F; Coenen, AML; Booij, LHDJ

2005-01-01

Background: This study assessed the feasibility of online calculation of the correlation integral (C(r)) aiming to apply C(r)derived statistics. For real-time application it is important to reduce calculation time. It is shown how our method works for EEG time series. Methods: To achieve online

Feasibility of real-time calculation of correlation integral derived statistics applied to EEG time series

NARCIS (Netherlands)

Broek, P.L.C. van den; Egmond, J. van; Rijn, C.M. van; Takens, F.; Coenen, A.M.L.; Booij, L.H.D.J.

2005-01-01

This study assessed the feasibility of online calculation of the correlation integral (C(r)) aiming to apply C(r)-derived statistics. For real-time application it is important to reduce calculation time. It is shown how our method works for EEG time series. Methods: To achieve online calculation of

Conservative multi-implicit integral deferred correction methods with adaptive mesh refinement

International Nuclear Information System (INIS)

Layton, A.T.

2004-01-01

In most models of reacting gas dynamics, the characteristic time scales of chemical reactions are much shorter than the hydrodynamic and diffusive time scales, rendering the reaction part of the model equations stiff. Moreover, nonlinear forcings may introduce into the solutions sharp gradients or shocks, the robust behavior and correct propagation of which require the use of specialized spatial discretization procedures. This study presents high-order conservative methods for the temporal integration of model equations of reacting flows. By means of a method of lines discretization on the flux difference form of the equations, these methods compute approximations to the cell-averaged or finite-volume solution. The temporal discretization is based on a multi-implicit generalization of integral deferred correction methods. The advection term is integrated explicitly, and the diffusion and reaction terms are treated implicitly but independently, with the splitting errors present in traditional operator splitting methods reduced via the integral deferred correction procedure. To reduce computational cost, time steps used to integrate processes with widely-differing time scales may differ in size. (author)

Improving Music Genre Classification by Short-Time Feature Integration

DEFF Research Database (Denmark)

Meng, Anders; Ahrendt, Peter; Larsen, Jan

2005-01-01

Many different short-time features, using time windows in the size of 10-30 ms, have been proposed for music segmentation, retrieval and genre classification. However, often the available time frame of the music to make the actual decision or comparison (the decision time horizon) is in the range...... of seconds instead of milliseconds. The problem of making new features on the larger time scale from the short-time features (feature integration) has only received little attention. This paper investigates different methods for feature integration and late information fusion for music genre classification...

Comparison between time-step-integration and probabilistic methods in seismic analysis of a linear structure

International Nuclear Information System (INIS)

Schneeberger, B.; Breuleux, R.

1977-01-01

Assuming that earthquake ground motion is a stationary time function, the seismic analysis of a linear structure can be done by probailistic methods using the 'power spectral density function' (PSD), instead of applying the more traditional time-step-integration using earthquake time histories (TH). A given structure was analysed both by PSD and TH methods computing and comparing 'floor response spectra'. The analysis using TH was performed for two different TH and different frequency intervals for the 'floor-response-spectra'. The analysis using PSD first produced PSD functions of the responses of the floors and these were then converted into 'foor-response-spectra'. Plots of the resulting 'floor-response-spectra' show: (1) The agreement of TH and PSD results is quite close. (2) The curves produced by PSD are much smoother than those produced by TH and mostly form an enelope of the latter. (3) The curves produced by TH are quite jagged with the location and magnitude of the peaks depending on the choice of frequencies at which the 'floor-response-spectra' were evaluated and on the choice of TH. (Auth.)

Pure radiation in space-time models that admit integration of the eikonal equation by the separation of variables method

Science.gov (United States)

Osetrin, Evgeny; Osetrin, Konstantin

2017-11-01

We consider space-time models with pure radiation, which admit integration of the eikonal equation by the method of separation of variables. For all types of these models, the equations of the energy-momentum conservation law are integrated. The resulting form of metric, energy density, and wave vectors of radiation as functions of metric for all types of spaces under consideration is presented. The solutions obtained can be used for any metric theories of gravitation.

A modified precise integration method based on Magnus expansion for transient response analysis of time varying dynamical structure

International Nuclear Information System (INIS)

Yue, Cong; Ren, Xingmin; Yang, Yongfeng; Deng, Wangqun

2016-01-01

This paper provides a precise and efficacious methodology for manifesting forced vibration response with respect to the time-variant linear rotational structure subjected to unbalanced excitation. A modified algorithm based on time step precise integration method and Magnus expansion is developed for instantaneous dynamic problems. The iterative solution is achieved by the ideology of transition and dimensional increment matrix. Numerical examples on a typical accelerating rotation system considering gyroscopic moment and mass unbalance force comparatively demonstrate the validity, effectiveness and accuracy with Newmark-β method. It is shown that the proposed algorithm has high accuracy without loss efficiency.

A multi-domain spectral method for time-fractional differential equations

Science.gov (United States)

Chen, Feng; Xu, Qinwu; Hesthaven, Jan S.

2015-07-01

This paper proposes an approach for high-order time integration within a multi-domain setting for time-fractional differential equations. Since the kernel is singular or nearly singular, two main difficulties arise after the domain decomposition: how to properly account for the history/memory part and how to perform the integration accurately. To address these issues, we propose a novel hybrid approach for the numerical integration based on the combination of three-term-recurrence relations of Jacobi polynomials and high-order Gauss quadrature. The different approximations used in the hybrid approach are justified theoretically and through numerical examples. Based on this, we propose a new multi-domain spectral method for high-order accurate time integrations and study its stability properties by identifying the method as a generalized linear method. Numerical experiments confirm hp-convergence for both time-fractional differential equations and time-fractional partial differential equations.

Formulation of an explicit-multiple-time-step time integration method for use in a global primitive equation grid model

Science.gov (United States)

Chao, W. C.

1982-01-01

With appropriate modifications, a recently proposed explicit-multiple-time-step scheme (EMTSS) is incorporated into the UCLA model. In this scheme, the linearized terms in the governing equations that generate the gravity waves are split into different vertical modes. Each mode is integrated with an optimal time step, and at periodic intervals these modes are recombined. The other terms are integrated with a time step dictated by the CFL condition for low-frequency waves. This large time step requires a special modification of the advective terms in the polar region to maintain stability. Test runs for 72 h show that EMTSS is a stable, efficient and accurate scheme.

Optimal distribution of integration time for intensity measurements in Stokes polarimetry.

Science.gov (United States)

Li, Xiaobo; Liu, Tiegen; Huang, Bingjing; Song, Zhanjie; Hu, Haofeng

2015-10-19

We consider the typical Stokes polarimetry system, which performs four intensity measurements to estimate a Stokes vector. We show that if the total integration time of intensity measurements is fixed, the variance of the Stokes vector estimator depends on the distribution of the integration time at four intensity measurements. Therefore, by optimizing the distribution of integration time, the variance of the Stokes vector estimator can be decreased. In this paper, we obtain the closed-form solution of the optimal distribution of integration time by employing Lagrange multiplier method. According to the theoretical analysis and real-world experiment, it is shown that the total variance of the Stokes vector estimator can be significantly decreased about 40% in the case discussed in this paper. The method proposed in this paper can effectively decrease the measurement variance and thus statistically improves the measurement accuracy of the polarimetric system.

Precise digital integration in wide time range: theory and realization

International Nuclear Information System (INIS)

Batrakov, A.M.; Pavlenko, A.V.

2017-01-01

The digital integration method based on using high-speed precision analog-to-digital converters (ADC) has become widely used over the recent years. The paper analyzes the limitations of this method that are caused by the signal properties, ADC sampling rate and noise spectral density of the ADC signal path. This analysis allowed creating digital integrators with accurate synchronization and achieving an integration error of less than 10 −5 in the time range from microseconds to tens of seconds. The structure of the integrator is described and its basic parameters are presented. The possibilities of different ADC chips in terms of their applicability to digital integrators are discussed. A comparison with other integrating devices is presented.

DRK methods for time-domain oscillator simulation

NARCIS (Netherlands)

Sevat, M.F.; Houben, S.H.M.J.; Maten, ter E.J.W.; Di Bucchianico, A.; Mattheij, R.M.M.; Peletier, M.A.

2006-01-01

This paper presents a new Runge-Kutta type integration method that is well-suited for time-domain simulation of oscillators. A unique property of the new method is that its damping characteristics can be controlled by a continuous parameter.

On Tuning PI Controllers for Integrating Plus Time Delay Systems

Directory of Open Access Journals (Sweden)

David Di Ruscio

2010-10-01

Full Text Available Some analytical results concerning PI controller tuning based on integrator plus time delay models are worked out and presented. A method for obtaining PI controller parameters, Kp=alpha/(k*tau, and, Ti=beta*tau, which ensures a given prescribed maximum time delay error, dtau_max, to time delay, tau, ratio parameter delta=dau_max/tau, is presented. The corner stone in this method, is a method product parameter, c=alpha*beta. Analytical relations between the PI controller parameters, Ti, and, Kp, and the time delay error parameter, delta, is presented, and we propose the setting, beta=c/a*(delta+1, and, alpha=a/(delta+1, which gives, Ti=c/a*(delta+1*tau, and Kp=a/((delta+1*k*tau, where the parameter, a, is constant in the method product parameter, c=alpha*beta. It also turns out that the integral time, Ti, is linear in, delta, and the proportional gain, Kp, inversely proportional to, delta+1. For the original Ziegler Nichols (ZN method this parameter is approximately, c=2.38, and the presented method may e.g., be used to obtain new modified ZN parameters with increased robustness margins, also documented in the paper.

Elucidating dynamic metabolic physiology through network integration of quantitative time-course metabolomics

DEFF Research Database (Denmark)

Bordbar, Aarash; Yurkovich, James T.; Paglia, Giuseppe

2017-01-01

The increasing availability of metabolomics data necessitates novel methods for deeper data analysis and interpretation. We present a flux balance analysis method that allows for the computation of dynamic intracellular metabolic changes at the cellular scale through integration of time-course ab......The increasing availability of metabolomics data necessitates novel methods for deeper data analysis and interpretation. We present a flux balance analysis method that allows for the computation of dynamic intracellular metabolic changes at the cellular scale through integration of time...

Non-integrability of time-dependent spherically symmetric Yang-Mills equations

Energy Technology Data Exchange (ETDEWEB)

Matinyan, S G; Prokhorenko, E B; Savvidy, G K

1988-03-07

The integrability of time-dependent spherically symmetric Yang-Mills equations is studied using the Fermi-Pasta-Ulam method. It is shown that the motion of this system is ergodic, while the system itself is non-integrable, i.e. manifests dynamical chaos.

analysis of large electromagnetic pulse simulators using the electric field integral equation method in time domain

International Nuclear Information System (INIS)

Jamali, J.; Aghajafari, R.; Moini, R.; Sadeghi, H.

2002-01-01

A time-domain approach is presented to calculate electromagnetic fields inside a large Electromagnetic Pulse (EMP) simulator. This type of EMP simulator is used for studying the effect of electromagnetic pulses on electrical apparatus in various structures such as vehicles, a reoplanes, etc. The simulator consists of three planar transmission lines. To solve the problem, we first model the metallic structure of the simulator as a grid of conducting wires. The numerical solution of the governing electric field integral equation is then obtained using the method of moments in time domain. To demonstrate the accuracy of the model, we consider a typical EMP simulator. The comparison of our results with those obtained experimentally in the literature validates the model introduced in this paper

Global Format for Conservative Time Integration in Nonlinear Dynamics

DEFF Research Database (Denmark)

Krenk, Steen

2014-01-01

The widely used classic collocation-based time integration procedures like Newmark, Generalized-alpha etc. generally work well within a framework of linear problems, but typically may encounter problems, when used in connection with essentially nonlinear structures. These problems are overcome....... In the present paper a conservative time integration algorithm is developed in a format using only the internal forces and the associated tangent stiffness at the specific time integration points. Thus, the procedure is computationally very similar to a collocation method, consisting of a series of nonlinear...... equivalent static load steps, easily implemented in existing computer codes. The paper considers two aspects: representation of nonlinear internal forces in a form that implies energy conservation, and the option of an algorithmic damping with the purpose of extracting energy from undesirable high...

An integrating factor matrix method to find first integrals

International Nuclear Information System (INIS)

Saputra, K V I; Quispel, G R W; Van Veen, L

2010-01-01

In this paper we develop an integrating factor matrix method to derive conditions for the existence of first integrals. We use this novel method to obtain first integrals, along with the conditions for their existence, for two- and three-dimensional Lotka-Volterra systems with constant terms. The results are compared to previous results obtained by other methods.

An energy-stable time-integrator for phase-field models

KAUST Repository

Vignal, Philippe

2016-12-27

We introduce a provably energy-stable time-integration method for general classes of phase-field models with polynomial potentials. We demonstrate how Taylor series expansions of the nonlinear terms present in the partial differential equations of these models can lead to expressions that guarantee energy-stability implicitly, which are second-order accurate in time. The spatial discretization relies on a mixed finite element formulation and isogeometric analysis. We also propose an adaptive time-stepping discretization that relies on a first-order backward approximation to give an error-estimator. This error estimator is accurate, robust, and does not require the computation of extra solutions to estimate the error. This methodology can be applied to any second-order accurate time-integration scheme. We present numerical examples in two and three spatial dimensions, which confirm the stability and robustness of the method. The implementation of the numerical schemes is done in PetIGA, a high-performance isogeometric analysis framework.

An energy-stable time-integrator for phase-field models

KAUST Repository

Vignal, Philippe; Collier, N.; Dalcin, Lisandro; Brown, D.L.; Calo, V.M.

2016-01-01

We introduce a provably energy-stable time-integration method for general classes of phase-field models with polynomial potentials. We demonstrate how Taylor series expansions of the nonlinear terms present in the partial differential equations of these models can lead to expressions that guarantee energy-stability implicitly, which are second-order accurate in time. The spatial discretization relies on a mixed finite element formulation and isogeometric analysis. We also propose an adaptive time-stepping discretization that relies on a first-order backward approximation to give an error-estimator. This error estimator is accurate, robust, and does not require the computation of extra solutions to estimate the error. This methodology can be applied to any second-order accurate time-integration scheme. We present numerical examples in two and three spatial dimensions, which confirm the stability and robustness of the method. The implementation of the numerical schemes is done in PetIGA, a high-performance isogeometric analysis framework.

Statistical Methods in Integrative Genomics

Science.gov (United States)

Richardson, Sylvia; Tseng, George C.; Sun, Wei

2016-01-01

Statistical methods in integrative genomics aim to answer important biology questions by jointly analyzing multiple types of genomic data (vertical integration) or aggregating the same type of data across multiple studies (horizontal integration). In this article, we introduce different types of genomic data and data resources, and then review statistical methods of integrative genomics, with emphasis on the motivation and rationale of these methods. We conclude with some summary points and future research directions. PMID:27482531

Integral methods in low-frequency electromagnetics

CERN Document Server

Solin, Pavel; Karban, Pavel; Ulrych, Bohus

2009-01-01

A modern presentation of integral methods in low-frequency electromagnetics This book provides state-of-the-art knowledge on integral methods in low-frequency electromagnetics. Blending theory with numerous examples, it introduces key aspects of the integral methods used in engineering as a powerful alternative to PDE-based models. Readers will get complete coverage of: The electromagnetic field and its basic characteristics An overview of solution methods Solutions of electromagnetic fields by integral expressions Integral and integrodifferential methods

Indoor integrated navigation and synchronous data acquisition method for Android smartphone

Science.gov (United States)

Hu, Chunsheng; Wei, Wenjian; Qin, Shiqiao; Wang, Xingshu; Habib, Ayman; Wang, Ruisheng

2015-08-01

Smartphones are widely used at present. Most smartphones have cameras and kinds of sensors, such as gyroscope, accelerometer and magnet meter. Indoor navigation based on smartphone is very important and valuable. According to the features of the smartphone and indoor navigation, a new indoor integrated navigation method is proposed, which uses MEMS (Micro-Electro-Mechanical Systems) IMU (Inertial Measurement Unit), camera and magnet meter of smartphone. The proposed navigation method mainly involves data acquisition, camera calibration, image measurement, IMU calibration, initial alignment, strapdown integral, zero velocity update and integrated navigation. Synchronous data acquisition of the sensors (gyroscope, accelerometer and magnet meter) and the camera is the base of the indoor navigation on the smartphone. A camera data acquisition method is introduced, which uses the camera class of Android to record images and time of smartphone camera. Two kinds of sensor data acquisition methods are introduced and compared. The first method records sensor data and time with the SensorManager of Android. The second method realizes open, close, data receiving and saving functions in C language, and calls the sensor functions in Java language with JNI interface. A data acquisition software is developed with JDK (Java Development Kit), Android ADT (Android Development Tools) and NDK (Native Development Kit). The software can record camera data, sensor data and time at the same time. Data acquisition experiments have been done with the developed software and Sumsang Note 2 smartphone. The experimental results show that the first method of sensor data acquisition is convenient but lost the sensor data sometimes, the second method is much better in real-time performance and much less in data losing. A checkerboard image is recorded, and the corner points of the checkerboard are detected with the Harris method. The sensor data of gyroscope, accelerometer and magnet meter have

Applications of hybrid time-frequency methods in nonlinear structural dynamics

International Nuclear Information System (INIS)

Politopoulos, I.; Piteau, Ph.; Borsoi, L.; Antunes, J.

2014-01-01

This paper presents a study on methods which may be used to compute the nonlinear response of systems whose linear properties are determined in the frequency or Laplace domain. Typically, this kind of situation may arise in soil-structure and fluid-structure interaction problems. In particular three methods are investigated: (a) the hybrid time-frequency method, (b) the computation of the convolution integral which requires an inverse Fourier or Laplace transform of the system's transfer function, and (c) the identification of an equivalent system defined in the time domain which may be solved with classical time integration methods. These methods are illustrated by their application to some simple, one degree of freedom, non-linear systems and their advantages and drawbacks are highlighted. (authors)

Detection of anatomical changes in lung cancer patients with 2D time-integrated, 2D time-resolved and 3D time-integrated portal dosimetry: a simulation study

Science.gov (United States)

Wolfs, Cecile J. A.; Brás, Mariana G.; Schyns, Lotte E. J. R.; Nijsten, Sebastiaan M. J. J. G.; van Elmpt, Wouter; Scheib, Stefan G.; Baltes, Christof; Podesta, Mark; Verhaegen, Frank

2017-08-01

The aim of this work is to assess the performance of 2D time-integrated (2D-TI), 2D time-resolved (2D-TR) and 3D time-integrated (3D-TI) portal dosimetry in detecting dose discrepancies between the planned and (simulated) delivered dose caused by simulated changes in the anatomy of lung cancer patients. For six lung cancer patients, tumor shift, tumor regression and pleural effusion are simulated by modifying their CT images. Based on the modified CT images, time-integrated (TI) and time-resolved (TR) portal dose images (PDIs) are simulated and 3D-TI doses are calculated. The modified and original PDIs and 3D doses are compared by a gamma analysis with various gamma criteria. Furthermore, the difference in the D 95% (ΔD 95%) of the GTV is calculated and used as a gold standard. The correlation between the gamma fail rate and the ΔD 95% is investigated, as well the sensitivity and specificity of all combinations of portal dosimetry method, gamma criteria and gamma fail rate threshold. On the individual patient level, there is a correlation between the gamma fail rate and the ΔD 95%, which cannot be found at the group level. The sensitivity and specificity analysis showed that there is not one combination of portal dosimetry method, gamma criteria and gamma fail rate threshold that can detect all simulated anatomical changes. This work shows that it will be more beneficial to relate portal dosimetry and DVH analysis on the patient level, rather than trying to quantify a relationship for a group of patients. With regards to optimizing sensitivity and specificity, different combinations of portal dosimetry method, gamma criteria and gamma fail rate should be used to optimally detect certain types of anatomical changes.

Detection of anatomical changes in lung cancer patients with 2D time-integrated, 2D time-resolved and 3D time-integrated portal dosimetry: a simulation study.

Science.gov (United States)

Wolfs, Cecile J A; Brás, Mariana G; Schyns, Lotte E J R; Nijsten, Sebastiaan M J J G; van Elmpt, Wouter; Scheib, Stefan G; Baltes, Christof; Podesta, Mark; Verhaegen, Frank

2017-07-12

The aim of this work is to assess the performance of 2D time-integrated (2D-TI), 2D time-resolved (2D-TR) and 3D time-integrated (3D-TI) portal dosimetry in detecting dose discrepancies between the planned and (simulated) delivered dose caused by simulated changes in the anatomy of lung cancer patients. For six lung cancer patients, tumor shift, tumor regression and pleural effusion are simulated by modifying their CT images. Based on the modified CT images, time-integrated (TI) and time-resolved (TR) portal dose images (PDIs) are simulated and 3D-TI doses are calculated. The modified and original PDIs and 3D doses are compared by a gamma analysis with various gamma criteria. Furthermore, the difference in the D 95% (ΔD 95% ) of the GTV is calculated and used as a gold standard. The correlation between the gamma fail rate and the ΔD 95% is investigated, as well the sensitivity and specificity of all combinations of portal dosimetry method, gamma criteria and gamma fail rate threshold. On the individual patient level, there is a correlation between the gamma fail rate and the ΔD 95% , which cannot be found at the group level. The sensitivity and specificity analysis showed that there is not one combination of portal dosimetry method, gamma criteria and gamma fail rate threshold that can detect all simulated anatomical changes. This work shows that it will be more beneficial to relate portal dosimetry and DVH analysis on the patient level, rather than trying to quantify a relationship for a group of patients. With regards to optimizing sensitivity and specificity, different combinations of portal dosimetry method, gamma criteria and gamma fail rate should be used to optimally detect certain types of anatomical changes.

Combined Monte Carlo and path-integral method for simulated library of time-resolved reflectance curves from layered tissue models

Science.gov (United States)

Wilson, Robert H.; Vishwanath, Karthik; Mycek, Mary-Ann

2009-02-01

Monte Carlo (MC) simulations are considered the "gold standard" for mathematical description of photon transport in tissue, but they can require large computation times. Therefore, it is important to develop simple and efficient methods for accelerating MC simulations, especially when a large "library" of related simulations is needed. A semi-analytical method involving MC simulations and a path-integral (PI) based scaling technique generated time-resolved reflectance curves from layered tissue models. First, a zero-absorption MC simulation was run for a tissue model with fixed scattering properties in each layer. Then, a closed-form expression for the average classical path of a photon in tissue was used to determine the percentage of time that the photon spent in each layer, to create a weighted Beer-Lambert factor to scale the time-resolved reflectance of the simulated zero-absorption tissue model. This method is a unique alternative to other scaling techniques in that it does not require the path length or number of collisions of each photon to be stored during the initial simulation. Effects of various layer thicknesses and absorption and scattering coefficients on the accuracy of the method will be discussed.

General new time formalism in the path integral

International Nuclear Information System (INIS)

Pak, N.K.; Sokmen, I.

1983-08-01

We describe a general method of applying point canonical transformations to the path integral followed by the corresponding new time transformations aimed at reducing an arbitrary one-dimensional problem into an exactly solvable form. Our result is independent of operator ordering ambiguities by construction. (author)

An integral nodal variational method for multigroup criticality calculations

International Nuclear Information System (INIS)

Lewis, E.E.; Tsoulfanidis, N.

2003-01-01

An integral formulation of the variational nodal method is presented and applied to a series of benchmark critically problems. The method combines an integral transport treatment of the even-parity flux within the spatial node with an odd-parity spherical harmonics expansion of the Lagrange multipliers at the node interfaces. The response matrices that result from this formulation are compatible with those in the VARIANT code at Argonne National Laboratory. Either homogeneous or heterogeneous nodes may be employed. In general, for calculations requiring higher-order angular approximations, the integral method yields solutions with comparable accuracy while requiring substantially less CPU time and memory than the standard spherical harmonics expansion using the same spatial approximations. (author)

A method to evaluate process performance by integrating time and resources

Science.gov (United States)

Wang, Yu; Wei, Qingjie; Jin, Shuang

2017-06-01

The purpose of process mining is to improve the existing process of the enterprise, so how to measure the performance of the process is particularly important. However, the current research on the performance evaluation method is still insufficient. The main methods of evaluation are mainly using time or resource. These basic statistics cannot evaluate process performance very well. In this paper, a method of evaluating the performance of the process based on time dimension and resource dimension is proposed. This method can be used to measure the utilization and redundancy of resources in the process. This paper will introduce the design principle and formula of the evaluation algorithm. Then, the design and the implementation of the evaluation method will be introduced. Finally, we will use the evaluating method to analyse the event log from a telephone maintenance process and propose an optimization plan.

The imaginary-time path integral and non-time-reversal-invariant saddle points of the Euclidean action

International Nuclear Information System (INIS)

Dasgupta, I.

1998-01-01

We discuss new bounce-like (but non-time-reversal-invariant) solutions to Euclidean equations of motion, which we dub boomerons. In the Euclidean path integral approach to quantum theories, boomerons make an imaginary contribution to the vacuum energy. The fake vacuum instability can be removed by cancelling boomeron contributions against contributions from time reversed boomerons (anti-boomerons). The cancellation rests on a sign choice whose significance is not completely understood in the path integral method. (orig.)

Methods optimization for the first time core critical

International Nuclear Information System (INIS)

Yan Liang

2014-01-01

The PWR reactor core commissioning programs the content of the first critical reactor physics experiment, and describes thc physical test method. However, all the methods arc not exactly the same but efficient. This article aims to enhance the reactor for the first time in the process of critical safety, shorten the overall time of critical physical test for the first time, and improve the integrity of critical physical test data for the first time and accuracy, eventually to improve the operation of the plant economic benefit adopting sectional dilution, power feedback for Doppler point improvement of physical test methods, and so on. (author)

Optimal distribution of integration time for intensity measurements in degree of linear polarization polarimetry.

Science.gov (United States)

Li, Xiaobo; Hu, Haofeng; Liu, Tiegen; Huang, Bingjing; Song, Zhanjie

2016-04-04

We consider the degree of linear polarization (DOLP) polarimetry system, which performs two intensity measurements at orthogonal polarization states to estimate DOLP. We show that if the total integration time of intensity measurements is fixed, the variance of the DOLP estimator depends on the distribution of integration time for two intensity measurements. Therefore, by optimizing the distribution of integration time, the variance of the DOLP estimator can be decreased. In this paper, we obtain the closed-form solution of the optimal distribution of integration time in an approximate way by employing Delta method and Lagrange multiplier method. According to the theoretical analyses and real-world experiments, it is shown that the variance of the DOLP estimator can be decreased for any value of DOLP. The method proposed in this paper can effectively decrease the measurement variance and thus statistically improve the measurement accuracy of the polarimetry system.

The Crank Nicolson Time Integrator for EMPHASIS.

Energy Technology Data Exchange (ETDEWEB)

McGregor, Duncan Alisdair Odum; Love, Edward; Kramer, Richard Michael Jack

2018-03-01

We investigate the use of implicit time integrators for finite element time domain approxi- mations of Maxwell's equations in vacuum. We discretize Maxwell's equations in time using Crank-Nicolson and in 3D space using compatible finite elements. We solve the system by taking a single step of Newton's method and inverting the Eddy-Current Schur complement allowing for the use of standard preconditioning techniques. This approach also generalizes to more complex material models that can include the Unsplit PML. We present verification results and demonstrate performance at CFL numbers up to 1000.

An analytical nodal method for time-dependent one-dimensional discrete ordinates problems

International Nuclear Information System (INIS)

Barros, R.C. de

1992-01-01

In recent years, relatively little work has been done in developing time-dependent discrete ordinates (S N ) computer codes. Therefore, the topic of time integration methods certainly deserves further attention. In this paper, we describe a new coarse-mesh method for time-dependent monoenergetic S N transport problesm in slab geometry. This numerical method preserves the analytic solution of the transverse-integrated S N nodal equations by constants, so we call our method the analytical constant nodal (ACN) method. For time-independent S N problems in finite slab geometry and for time-dependent infinite-medium S N problems, the ACN method generates numerical solutions that are completely free of truncation errors. Bsed on this positive feature, we expect the ACN method to be more accurate than conventional numerical methods for S N transport calculations on coarse space-time grids

Energy conservation in Newmark based time integration algorithms

DEFF Research Database (Denmark)

Krenk, Steen

2006-01-01

Energy balance equations are established for the Newmark time integration algorithm, and for the derived algorithms with algorithmic damping introduced via averaging, the so-called a-methods. The energy balance equations form a sequence applicable to: Newmark integration of the undamped equations...... of motion, an extended form including structural damping, and finally the generalized form including structural as well as algorithmic damping. In all three cases the expression for energy, appearing in the balance equation, is the mechanical energy plus some additional terms generated by the discretization...

Path integration on space times with symmetry

International Nuclear Information System (INIS)

Low, S.G.

1985-01-01

Path integration on space times with symmetry is investigated using a definition of path integration of Gaussian integrators. Gaussian integrators, systematically developed using the theory of projective distributions, may be defined in terms of a Jacobi operator Green function. This definition of the path integral yields a semiclassical expansion of the propagator which is valid on caustics. The semiclassical approximation to the free particle propagator on symmetric and reductive homogeneous spaces is computed in terms of the complete solution of the Jacobi equation. The results are used to test the validity of using the Schwinger-DeWitt transform to compute an approximation to the coincidence limit of a field theory Green function from a WKB propagator. The method is found not to be valid except for certain special cases. These cases include manifolds constructed from the direct product of flat space and group manifolds, on which the free particle WKB approximation is exact and two sphere. The multiple geodesic contribution to 2 > on Schwarzschild in the neighborhood of rho = 3M is computed using the transform

A Scalable, Timing-Safe, Network-on-Chip Architecture with an Integrated Clock Distribution Method

DEFF Research Database (Denmark)

Bjerregaard, Tobias; Stensgaard, Mikkel Bystrup; Sparsø, Jens

2007-01-01

Growing system sizes together with increasing performance variability are making globally synchronous operation hard to realize. Mesochronous clocking constitutes a possible solution to the problems faced. The most fundamental of problems faced when communicating between mesochronously clocked re...... is based purely on local observations. It is demonstrated with a 90 nm CMOS standard cell network-on-chip design which implements completely timing-safe, global communication in a modular system......Growing system sizes together with increasing performance variability are making globally synchronous operation hard to realize. Mesochronous clocking constitutes a possible solution to the problems faced. The most fundamental of problems faced when communicating between mesochronously clocked...... regions concerns the possibility of data corruption caused by metastability. This paper presents an integrated communication and mesochronous clocking strategy, which avoids timing related errors while maintaining a globally synchronous system perspective. The architecture is scalable as timing integrity...

[A new measurement method of time-resolved spectrum].

Science.gov (United States)

Shi, Zhi-gang; Huang, Shi-hua; Liang, Chun-jun; Lei, Quan-sheng

2007-02-01

A new method for measuring time-resolved spectrum (TRS) is brought forward. Programming with assemble language controlled the micro-control-processor (AT89C51), and a kind of peripheral circuit constituted the drive circuit, which drived the stepping motor to run the monochromator. So the light of different kinds of expected wavelength could be obtained. The optical signal was transformed to electrical signal by optical-to-electrical transform with the help of photomultiplier tube (Hamamatsu 1P28). The electrical signal of spectrum data was transmitted to the oscillograph. Connecting the two serial interfaces of RS232 between the oscillograph and computer, the electrical signal of spectrum data could be transmitted to computer for programming to draw the attenuation curve and time-resolved spectrum (TRS) of the swatch. The method for measuring time-resolved spectrum (TRS) features parallel measurement in time scale but serial measurement in wavelength scale. Time-resolved spectrum (TRS) and integrated emission spectrum of Tb3+ in swatch Tb(o-BBA)3 phen were measured using this method. Compared with the real time-resolved spectrum (TRS). It was validated to be feasible, credible and convenient. The 3D spectra of fluorescence intensity-wavelength-time, and the integrated spectrum of the swatch Tb(o-BBA)3 phen are given.

Non-integrability of time-dependent spherically symmetric Yang-Mills equations

International Nuclear Information System (INIS)

Matinyan, S.G.; Prokhorenko, E.V.; Savvidy, G.K.

1986-01-01

The integrability of time-dependent spherically symmetric Yang-Mills equations is studied using the Fermi-Pasta-Ulam method. The phase space of this system is shown to have no quasi-periodic motion specific for integrable systems. In particular, the well-known Wu-Yang static solution is unstable, so its vicinity in phase is the stochasticity region

Indirect methods for wake potential integration

International Nuclear Information System (INIS)

Zagorodnov, I.

2006-05-01

The development of the modern accelerator and free-electron laser projects requires to consider wake fields of very short bunches in arbitrary three dimensional structures. To obtain the wake numerically by direct integration is difficult, since it takes a long time for the scattered fields to catch up to the bunch. On the other hand no general algorithm for indirect wake field integration is available in the literature so far. In this paper we review the know indirect methods to compute wake potentials in rotationally symmetric and cavity-like three dimensional structures. For arbitrary three dimensional geometries we introduce several new techniques and test them numerically. (Orig.)

Alternative containment integrity test methods, an overview of possible techniques

International Nuclear Information System (INIS)

Spletzer, B.L.

1986-01-01

A study is being conducted to develop and analyze alternative methods for testing of containment integrity. The study is focused on techniques for continuously monitoring containment integrity to provide rapid detection of existing leaks, thus providing greater certainty of the integrity of the containment at any time. The study is also intended to develop techniques applicable to the currently required Type A integrated leakage rate tests. A brief discussion of the range of alternative methods currently being considered is presented. The methods include applicability to all major containment types, operating and shutdown plant conditions, and quantitative and qualitative leakage measurements. The techniques are analyzed in accordance with the current state of knowledge of each method. The bulk of the techniques discussed are in the conceptual stage, have not been tested in actual plant conditions, and are presented here as a possible future direction for evaluating containment integrity. Of the methods considered, no single method provides optimum performance for all containment types. Several methods are limited in the types of containment for which they are applicable. The results of the study to date indicate that techniques for continuous monitoring of containment integrity exist for many plants and may be implemented at modest cost

On the initial condition problem of the time domain PMCHWT surface integral equation

KAUST Repository

Uysal, Ismail Enes

2017-05-13

Non-physical, linearly increasing and constant current components are induced in marching on-in-time solution of time domain surface integral equations when initial conditions on time derivatives of (unknown) equivalent currents are not enforced properly. This problem can be remedied by solving the time integral of the surface integral for auxiliary currents that are defined to be the time derivatives of the equivalent currents. Then the equivalent currents are obtained by numerically differentiating the auxiliary ones. In this work, this approach is applied to the marching on-in-time solution of the time domain Poggio-Miller-Chan-Harrington-Wu-Tsai surface integral equation enforced on dispersive/plasmonic scatterers. Accuracy of the proposed method is demonstrated by a numerical example.

A new integral method for solving the point reactor neutron kinetics equations

International Nuclear Information System (INIS)

Li Haofeng; Chen Wenzhen; Luo Lei; Zhu Qian

2009-01-01

A numerical integral method that efficiently provides the solution of the point kinetics equations by using the better basis function (BBF) for the approximation of the neutron density in one time step integrations is described and investigated. The approach is based on an exact analytic integration of the neutron density equation, where the stiffness of the equations is overcome by the fully implicit formulation. The procedure is tested by using a variety of reactivity functions, including step reactivity insertion, ramp input and oscillatory reactivity changes. The solution of the better basis function method is compared to other analytical and numerical solutions of the point reactor kinetics equations. The results show that selecting a better basis function can improve the efficiency and accuracy of this integral method. The better basis function method can be used in real time forecasting for power reactors in order to prevent reactivity accidents.

Two new solutions to the third-order symplectic integration method

International Nuclear Information System (INIS)

Iwatsu, Reima

2009-01-01

Two new solutions are obtained for the symplecticity conditions of explicit third-order partitioned Runge-Kutta time integration method. One of them has larger stability limit and better dispersion property than the Ruth's method.

An iterated Radau method for time-dependent PDE's

NARCIS (Netherlands)

S. Pérez-Rodríguez; S. González-Pinto; B.P. Sommeijer (Ben)

2008-01-01

htmlabstractThis paper is concerned with the time integration of semi-discretized, multi-dimensional PDEs of advection-diffusion-reaction type. To cope with the stiffness of these ODEs, an implicit method has been selected, viz., the two-stage, third-order Radau IIA method. The main topic of this

A block-iterative nodal integral method for forced convection problems

International Nuclear Information System (INIS)

Decker, W.J.; Dorning, J.J.

1992-01-01

A new efficient iterative nodal integral method for the time-dependent two- and three-dimensional incompressible Navier-Stokes equations has been developed. Using the approach introduced by Azmy and Droning to develop nodal mehtods with high accuracy on coarse spatial grids for two-dimensional steady-state problems and extended to coarse two-dimensional space-time grids by Wilson et al. for thermal convection problems, we have developed a new iterative nodal integral method for the time-dependent Navier-Stokes equations for mechanically forced convection. A new, extremely efficient block iterative scheme is employed to invert the Jacobian within each of the Newton-Raphson iterations used to solve the final nonlinear discrete-variable equations. By taking advantage of the special structure of the Jacobian, this scheme greatly reduces memory requirements. The accuracy of the overall method is illustrated by appliying it to the time-dependent version of the classic two-dimensional driven cavity problem of computational fluid dynamics

A time-domain method to generate artificial time history from a given reference response spectrum

Energy Technology Data Exchange (ETDEWEB)

Shin, Gang Sik [Korea Institute of Nuclear Safety, Daejeon (Korea, Republic of); Song, Oh Seop [Dept. of Mechanical Engineering, Chungnam National University, Daejeon (Korea, Republic of)

2016-06-15

Seismic qualification by test is widely used as a way to show the integrity and functionality of equipment that is related to the overall safety of nuclear power plants. Another means of seismic qualification is by direct integration analysis. Both approaches require a series of time histories as an input. However, in most cases, the possibility of using real earthquake data is limited. Thus, artificial time histories are widely used instead. In many cases, however, response spectra are given. Thus, most of the artificial time histories are generated from the given response spectra. Obtaining the response spectrum from a given time history is straightforward. However, the procedure for generating artificial time histories from a given response spectrum is difficult and complex to understand. Thus, this paper presents a simple time-domain method for generating a time history from a given response spectrum; the method was shown to satisfy conditions derived from nuclear regulatory guidance.

A time-domain method to generate artificial time history from a given reference response spectrum

International Nuclear Information System (INIS)

Shin, Gang Sik; Song, Oh Seop

2016-01-01

Seismic qualification by test is widely used as a way to show the integrity and functionality of equipment that is related to the overall safety of nuclear power plants. Another means of seismic qualification is by direct integration analysis. Both approaches require a series of time histories as an input. However, in most cases, the possibility of using real earthquake data is limited. Thus, artificial time histories are widely used instead. In many cases, however, response spectra are given. Thus, most of the artificial time histories are generated from the given response spectra. Obtaining the response spectrum from a given time history is straightforward. However, the procedure for generating artificial time histories from a given response spectrum is difficult and complex to understand. Thus, this paper presents a simple time-domain method for generating a time history from a given response spectrum; the method was shown to satisfy conditions derived from nuclear regulatory guidance

An open-chain imaginary-time path-integral sampling approach to the calculation of approximate symmetrized quantum time correlation functions

Science.gov (United States)

Cendagorta, Joseph R.; Bačić, Zlatko; Tuckerman, Mark E.

2018-03-01

We introduce a scheme for approximating quantum time correlation functions numerically within the Feynman path integral formulation. Starting with the symmetrized version of the correlation function expressed as a discretized path integral, we introduce a change of integration variables often used in the derivation of trajectory-based semiclassical methods. In particular, we transform to sum and difference variables between forward and backward complex-time propagation paths. Once the transformation is performed, the potential energy is expanded in powers of the difference variables, which allows us to perform the integrals over these variables analytically. The manner in which this procedure is carried out results in an open-chain path integral (in the remaining sum variables) with a modified potential that is evaluated using imaginary-time path-integral sampling rather than requiring the generation of a large ensemble of trajectories. Consequently, any number of path integral sampling schemes can be employed to compute the remaining path integral, including Monte Carlo, path-integral molecular dynamics, or enhanced path-integral molecular dynamics. We believe that this approach constitutes a different perspective in semiclassical-type approximations to quantum time correlation functions. Importantly, we argue that our approximation can be systematically improved within a cumulant expansion formalism. We test this approximation on a set of one-dimensional problems that are commonly used to benchmark approximate quantum dynamical schemes. We show that the method is at least as accurate as the popular ring-polymer molecular dynamics technique and linearized semiclassical initial value representation for correlation functions of linear operators in most of these examples and improves the accuracy of correlation functions of nonlinear operators.

An open-chain imaginary-time path-integral sampling approach to the calculation of approximate symmetrized quantum time correlation functions.

Science.gov (United States)

Cendagorta, Joseph R; Bačić, Zlatko; Tuckerman, Mark E

2018-03-14

We introduce a scheme for approximating quantum time correlation functions numerically within the Feynman path integral formulation. Starting with the symmetrized version of the correlation function expressed as a discretized path integral, we introduce a change of integration variables often used in the derivation of trajectory-based semiclassical methods. In particular, we transform to sum and difference variables between forward and backward complex-time propagation paths. Once the transformation is performed, the potential energy is expanded in powers of the difference variables, which allows us to perform the integrals over these variables analytically. The manner in which this procedure is carried out results in an open-chain path integral (in the remaining sum variables) with a modified potential that is evaluated using imaginary-time path-integral sampling rather than requiring the generation of a large ensemble of trajectories. Consequently, any number of path integral sampling schemes can be employed to compute the remaining path integral, including Monte Carlo, path-integral molecular dynamics, or enhanced path-integral molecular dynamics. We believe that this approach constitutes a different perspective in semiclassical-type approximations to quantum time correlation functions. Importantly, we argue that our approximation can be systematically improved within a cumulant expansion formalism. We test this approximation on a set of one-dimensional problems that are commonly used to benchmark approximate quantum dynamical schemes. We show that the method is at least as accurate as the popular ring-polymer molecular dynamics technique and linearized semiclassical initial value representation for correlation functions of linear operators in most of these examples and improves the accuracy of correlation functions of nonlinear operators.

Terminal current interpolation for multirate time integration of hierarchical IC models

NARCIS (Netherlands)

Verhoeven, A.; Maten, ter E.J.W.; Dohmen, J.J.; Tasic, B.; Mattheij, R.M.M.; Fitt, A.D.; Norbury, J.; Ockendon, H.; Wilson, E.

2010-01-01

Multirate time-integration methods [3–5] appear to be attractive for initial value problems for DAEs with latency or multirate behaviour. Latency means that parts of the circuit are constant or slowly time-varying during a certain time interval, while multirate behaviour means that some variables

Integral equation methods for electromagnetics

CERN Document Server

Volakis, John

2012-01-01

This text/reference is a detailed look at the development and use of integral equation methods for electromagnetic analysis, specifically for antennas and radar scattering. Developers and practitioners will appreciate the broad-based approach to understanding and utilizing integral equation methods and the unique coverage of historical developments that led to the current state-of-the-art. In contrast to existing books, Integral Equation Methods for Electromagnetics lays the groundwork in the initial chapters so students and basic users can solve simple problems and work their way up to the mo

Classification Method in Integrated Information Network Using Vector Image Comparison

Directory of Open Access Journals (Sweden)

Zhou Yuan

2014-05-01

Full Text Available Wireless Integrated Information Network (WMN consists of integrated information that can get data from its surrounding, such as image, voice. To transmit information, large resource is required which decreases the service time of the network. In this paper we present a Classification Approach based on Vector Image Comparison (VIC for WMN that improve the service time of the network. The available methods for sub-region selection and conversion are also proposed.

Steffensen's Integral Inequality on Time Scales

Directory of Open Access Journals (Sweden)

Ozkan Umut Mutlu

2007-01-01

Full Text Available We establish generalizations of Steffensen's integral inequality on time scales via the diamond- dynamic integral, which is defined as a linear combination of the delta and nabla integrals.

Parareal algorithms with local time-integrators for time fractional differential equations

Science.gov (United States)

Wu, Shu-Lin; Zhou, Tao

2018-04-01

It is challenge work to design parareal algorithms for time-fractional differential equations due to the historical effect of the fractional operator. A direct extension of the classical parareal method to such equations will lead to unbalance computational time in each process. In this work, we present an efficient parareal iteration scheme to overcome this issue, by adopting two recently developed local time-integrators for time fractional operators. In both approaches, one introduces auxiliary variables to localized the fractional operator. To this end, we propose a new strategy to perform the coarse grid correction so that the auxiliary variables and the solution variable are corrected separately in a mixed pattern. It is shown that the proposed parareal algorithm admits robust rate of convergence. Numerical examples are presented to support our conclusions.

Diverse methods for integrable models

NARCIS (Netherlands)

Fehér, G.

2017-01-01

This thesis is centered around three topics, sharing integrability as a common theme. This thesis explores different methods in the field of integrable models. The first two chapters are about integrable lattice models in statistical physics. The last chapter describes an integrable quantum chain.

Iterative algorithm for the volume integral method for magnetostatics problems

International Nuclear Information System (INIS)

Pasciak, J.E.

1980-11-01

Volume integral methods for solving nonlinear magnetostatics problems are considered in this paper. The integral method is discretized by a Galerkin technique. Estimates are given which show that the linearized problems are well conditioned and hence easily solved using iterative techniques. Comparisons of iterative algorithms with the elimination method of GFUN3D shows that the iterative method gives an order of magnitude improvement in computational time as well as memory requirements for large problems. Computational experiments for a test problem as well as a double layer dipole magnet are given. Error estimates for the linearized problem are also derived

Connection between Feynman integrals having different values of the space-time dimension

International Nuclear Information System (INIS)

Tarasov, O.V.

1996-05-01

A systematic algorithm for obtaining recurrence relations for dimensionally regularized Feynman integrals w.r.t. the space-time dimension d is proposed. The relation between d and d-2 dimensional integrals is given in terms of a differential operator for which an explicit formula can be obtained for each Feynman diagram. We show how the method works for one-, two- and three-loop integrals. The new recurrence relations w.r.t. d are complementary to the recurrence relations which derive from the method of integration by parts. We find that the problem of the irreducible numerators in Feynman integrals can be naturally solved in the framework of the proposed generalized recurrence relations. (orig.)

Development of a precise long-time digital integrator for magnetic measurements in a tokamak

Energy Technology Data Exchange (ETDEWEB)

Kurihara, Kenichi; Kawamata, Youichi [Japan Atomic Energy Research Inst., Naka, Ibaraki (Japan). Naka Fusion Research Establishment

1997-10-01

Long-time D-T burning operation in a tokamak requires that a magnetic sensor must work in an environment of 14-MeV intense neutron field, and that the measurement system must output precise magnetic field values. A method of time-integration of voltage produced in a simple pick-up coil seems to have preferable features of good time response, easy maintenance, and resistance to neutron irradiation. However, an inevitably-produced signal drift makes it difficult to apply the method to the long-time integral operation. To solve this problem, we have developed a new digital integrator (a voltage-to-frequency converter and an up-down counter) with testing the trial boards in the JT-60 magnetic measurements. This reports all of the problems and their measures through the development steps in details, and shows how to apply this method to the ITER operation. (author)

Analytical solutions for prediction of the ignition time of wood particles based on a time and space integral method

International Nuclear Information System (INIS)

Haseli, Y.; Oijen, J.A. van; Goey, L.P.H. de

2012-01-01

Highlights: ► A simple model for prediction of the ignition time of a wood particle is presented. ► The formulation is given for both thermally thin and thermally thick particles. ► Transition from thermally thin to thick regime occurs at a critical particle size. ► The model is validated against a numerical model and various experimental data. - Abstract: The main idea of this paper is to establish a simple approach for prediction of the ignition time of a wood particle assuming that the thermo-physical properties remain constant and ignition takes place at a characteristic ignition temperature. Using a time and space integral method, explicit relationships are derived for computation of the ignition time of particles of three common shapes (slab, cylinder and sphere), which may be characterized as thermally thin or thermally thick. It is shown through a dimensionless analysis that the dimensionless ignition time can be described as a function of non-dimensional ignition temperature, reactor temperature or external incident heat flux, and parameter K which represents the ratio of conduction heat transfer to the external radiation heat transfer. The numerical results reveal that for the dimensionless ignition temperature between 1.25 and 2.25 and for values of K up to 8000 (corresponding to woody materials), the variation of the ignition time of a thermally thin particle with K and the dimensionless ignition temperature is linear, whereas the dependence of the ignition time of a thermally thick particle on the above two parameters obeys a quadratic function. Furthermore, it is shown that the transition from the regime of thermally thin to the regime of thermally thick occurs at K cr (corresponding to a critical size of particle) which is found to be independent of the particle shape. The model is validated by comparing the predicted and the measured ignition time of several wood particles obtained from different sources. Good agreement is achieved which

An Accurate Integral Method for Vibration Signal Based on Feature Information Extraction

Directory of Open Access Journals (Sweden)

Yong Zhu

2015-01-01

Full Text Available After summarizing the advantages and disadvantages of current integral methods, a novel vibration signal integral method based on feature information extraction was proposed. This method took full advantage of the self-adaptive filter characteristic and waveform correction feature of ensemble empirical mode decomposition in dealing with nonlinear and nonstationary signals. This research merged the superiorities of kurtosis, mean square error, energy, and singular value decomposition on signal feature extraction. The values of the four indexes aforementioned were combined into a feature vector. Then, the connotative characteristic components in vibration signal were accurately extracted by Euclidean distance search, and the desired integral signals were precisely reconstructed. With this method, the interference problem of invalid signal such as trend item and noise which plague traditional methods is commendably solved. The great cumulative error from the traditional time-domain integral is effectively overcome. Moreover, the large low-frequency error from the traditional frequency-domain integral is successfully avoided. Comparing with the traditional integral methods, this method is outstanding at removing noise and retaining useful feature information and shows higher accuracy and superiority.

Fully implicit solution of large-scale non-equilibrium radiation diffusion with high order time integration

International Nuclear Information System (INIS)

Brown, Peter N.; Shumaker, Dana E.; Woodward, Carol S.

2005-01-01

We present a solution method for fully implicit radiation diffusion problems discretized on meshes having millions of spatial zones. This solution method makes use of high order in time integration techniques, inexact Newton-Krylov nonlinear solvers, and multigrid preconditioners. We explore the advantages and disadvantages of high order time integration methods for the fully implicit formulation on both two- and three-dimensional problems with tabulated opacities and highly nonlinear fusion source terms

Improvements of the integral transport theory method

International Nuclear Information System (INIS)

Kavenoky, A.; Lam-Hime, M.; Stankovski, Z.

1979-01-01

The integral transport theory is widely used in practical reactor design calculations however it is computer time consuming for two dimensional calculations of large media. In the first part of this report a new treatment is presented; it is based on the Galerkin method: inside each region the total flux is expanded over a three component basis. Numerical comparison shows that this method can considerably reduce the computing time. The second part of the this report is devoted to homogeneization theory: a straightforward calculation of the fundamental mode for an heterogeneous cell is presented. At first general presentation of the problem is given, then it is simplified to plane geometry and numerical results are presented

Adaptive time-stepping Monte Carlo integration of Coulomb collisions

Science.gov (United States)

Särkimäki, K.; Hirvijoki, E.; Terävä, J.

2018-01-01

We report an accessible and robust tool for evaluating the effects of Coulomb collisions on a test particle in a plasma that obeys Maxwell-Jüttner statistics. The implementation is based on the Beliaev-Budker collision integral which allows both the test particle and the background plasma to be relativistic. The integration method supports adaptive time stepping, which is shown to greatly improve the computational efficiency. The Monte Carlo method is implemented for both the three-dimensional particle momentum space and the five-dimensional guiding center phase space. Detailed description is provided for both the physics and implementation of the operator. The focus is in adaptive integration of stochastic differential equations, which is an overlooked aspect among existing Monte Carlo implementations of Coulomb collision operators. We verify that our operator converges to known analytical results and demonstrate that careless implementation of the adaptive time step can lead to severely erroneous results. The operator is provided as a self-contained Fortran 95 module and can be included into existing orbit-following tools that trace either the full Larmor motion or the guiding center dynamics. The adaptive time-stepping algorithm is expected to be useful in situations where the collision frequencies vary greatly over the course of a simulation. Examples include the slowing-down of fusion products or other fast ions, and the Dreicer generation of runaway electrons as well as the generation of fast ions or electrons with ion or electron cyclotron resonance heating.

New Approaches to Aluminum Integral Foam Production with Casting Methods

Directory of Open Access Journals (Sweden)

Ahmet Güner

2015-08-01

Full Text Available Integral foam has been used in the production of polymer materials for a long time. Metal integral foam casting systems are obtained by transferring and adapting polymer injection technology. Metal integral foam produced by casting has a solid skin at the surface and a foam core. Producing near-net shape reduces production expenses. Insurance companies nowadays want the automotive industry to use metallic foam parts because of their higher impact energy absorption properties. In this paper, manufacturing processes of aluminum integral foam with casting methods will be discussed.

Path integral density matrix dynamics: a method for calculating time-dependent properties in thermal adiabatic and non-adiabatic systems.

Science.gov (United States)

Habershon, Scott

2013-09-14

We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.

Path integral density matrix dynamics: A method for calculating time-dependent properties in thermal adiabatic and non-adiabatic systems

International Nuclear Information System (INIS)

Habershon, Scott

2013-01-01

We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency

A numerical integration-based yield estimation method for integrated circuits

International Nuclear Information System (INIS)

Liang Tao; Jia Xinzhang

2011-01-01

A novel integration-based yield estimation method is developed for yield optimization of integrated circuits. This method tries to integrate the joint probability density function on the acceptability region directly. To achieve this goal, the simulated performance data of unknown distribution should be converted to follow a multivariate normal distribution by using Box-Cox transformation (BCT). In order to reduce the estimation variances of the model parameters of the density function, orthogonal array-based modified Latin hypercube sampling (OA-MLHS) is presented to generate samples in the disturbance space during simulations. The principle of variance reduction of model parameters estimation through OA-MLHS together with BCT is also discussed. Two yield estimation examples, a fourth-order OTA-C filter and a three-dimensional (3D) quadratic function are used for comparison of our method with Monte Carlo based methods including Latin hypercube sampling and importance sampling under several combinations of sample sizes and yield values. Extensive simulations show that our method is superior to other methods with respect to accuracy and efficiency under all of the given cases. Therefore, our method is more suitable for parametric yield optimization. (semiconductor integrated circuits)

A numerical integration-based yield estimation method for integrated circuits

Energy Technology Data Exchange (ETDEWEB)

Liang Tao; Jia Xinzhang, E-mail: tliang@yahoo.cn [Key Laboratory of Ministry of Education for Wide Bandgap Semiconductor Materials and Devices, School of Microelectronics, Xidian University, Xi' an 710071 (China)

2011-04-15

A novel integration-based yield estimation method is developed for yield optimization of integrated circuits. This method tries to integrate the joint probability density function on the acceptability region directly. To achieve this goal, the simulated performance data of unknown distribution should be converted to follow a multivariate normal distribution by using Box-Cox transformation (BCT). In order to reduce the estimation variances of the model parameters of the density function, orthogonal array-based modified Latin hypercube sampling (OA-MLHS) is presented to generate samples in the disturbance space during simulations. The principle of variance reduction of model parameters estimation through OA-MLHS together with BCT is also discussed. Two yield estimation examples, a fourth-order OTA-C filter and a three-dimensional (3D) quadratic function are used for comparison of our method with Monte Carlo based methods including Latin hypercube sampling and importance sampling under several combinations of sample sizes and yield values. Extensive simulations show that our method is superior to other methods with respect to accuracy and efficiency under all of the given cases. Therefore, our method is more suitable for parametric yield optimization. (semiconductor integrated circuits)

SOCIAL INTEGRATION: TESTING ANTECEDENTS OF TIME SPENT ONLINE

Directory of Open Access Journals (Sweden)

Lily Suriani Mohd Arif

2013-07-01

Full Text Available The literature on the relationship of social integration and time spent onlineprovides conflicting evidence of the relationship of social integration with timespent online. The study identifies and highlightsthe controversy and attempts toclarify the relationship of social integration withtime spent online bydecomposing the construct social integration into its affective and behavioraldimensions . Thestudy tests antecedents and effects of time spent online in arandom sample of senior level undergraduate students at a public university inMalaysia. The findings indicated that while self-report measures of behavioralsocial integration did not predict time spent online, and, the affective socialintegration had an inverse relationship with time spent online.

A high-order boundary integral method for surface diffusions on elastically stressed axisymmetric rods.

Science.gov (United States)

Li, Xiaofan; Nie, Qing

2009-07-01

Many applications in materials involve surface diffusion of elastically stressed solids. Study of singularity formation and long-time behavior of such solid surfaces requires accurate simulations in both space and time. Here we present a high-order boundary integral method for an elastically stressed solid with axi-symmetry due to surface diffusions. In this method, the boundary integrals for isotropic elasticity in axi-symmetric geometry are approximated through modified alternating quadratures along with an extrapolation technique, leading to an arbitrarily high-order quadrature; in addition, a high-order (temporal) integration factor method, based on explicit representation of the mean curvature, is used to reduce the stability constraint on time-step. To apply this method to a periodic (in axial direction) and axi-symmetric elastically stressed cylinder, we also present a fast and accurate summation method for the periodic Green's functions of isotropic elasticity. Using the high-order boundary integral method, we demonstrate that in absence of elasticity the cylinder surface pinches in finite time at the axis of the symmetry and the universal cone angle of the pinching is found to be consistent with the previous studies based on a self-similar assumption. In the presence of elastic stress, we show that a finite time, geometrical singularity occurs well before the cylindrical solid collapses onto the axis of symmetry, and the angle of the corner singularity on the cylinder surface is also estimated.

Computing thermal Wigner densities with the phase integration method

International Nuclear Information System (INIS)

Beutier, J.; Borgis, D.; Vuilleumier, R.; Bonella, S.

2014-01-01

We discuss how the Phase Integration Method (PIM), recently developed to compute symmetrized time correlation functions [M. Monteferrante, S. Bonella, and G. Ciccotti, Mol. Phys. 109, 3015 (2011)], can be adapted to sampling/generating the thermal Wigner density, a key ingredient, for example, in many approximate schemes for simulating quantum time dependent properties. PIM combines a path integral representation of the density with a cumulant expansion to represent the Wigner function in a form calculable via existing Monte Carlo algorithms for sampling noisy probability densities. The method is able to capture highly non-classical effects such as correlation among the momenta and coordinates parts of the density, or correlations among the momenta themselves. By using alternatives to cumulants, it can also indicate the presence of negative parts of the Wigner density. Both properties are demonstrated by comparing PIM results to those of reference quantum calculations on a set of model problems

Computing thermal Wigner densities with the phase integration method.

Science.gov (United States)

Beutier, J; Borgis, D; Vuilleumier, R; Bonella, S

2014-08-28

We discuss how the Phase Integration Method (PIM), recently developed to compute symmetrized time correlation functions [M. Monteferrante, S. Bonella, and G. Ciccotti, Mol. Phys. 109, 3015 (2011)], can be adapted to sampling/generating the thermal Wigner density, a key ingredient, for example, in many approximate schemes for simulating quantum time dependent properties. PIM combines a path integral representation of the density with a cumulant expansion to represent the Wigner function in a form calculable via existing Monte Carlo algorithms for sampling noisy probability densities. The method is able to capture highly non-classical effects such as correlation among the momenta and coordinates parts of the density, or correlations among the momenta themselves. By using alternatives to cumulants, it can also indicate the presence of negative parts of the Wigner density. Both properties are demonstrated by comparing PIM results to those of reference quantum calculations on a set of model problems.

Unconditionally Energy Stable Implicit Time Integration: Application to Multibody System Analysis and Design

DEFF Research Database (Denmark)

Chen, Shanshin; Tortorelli, Daniel A.; Hansen, John Michael

1999-01-01

of ordinary diffferential equations is employed to avoid the instabilities associated with the direct integrations of differential-algebraic equations. To extend the unconditional stability of the implicit Newmark method to nonlinear dynamic systems, a discrete energy balance is enforced. This constraint......Advances in computer hardware and improved algorithms for multibody dynamics over the past decade have generated widespread interest in real-time simulations of multibody mechanics systems. At the heart of the widely used algorithms for multibody dynamics are a choice of coordinates which define...... the kinmatics of the system, and a choice of time integrations algorithms. The current approach uses a non-dissipative implict Newmark method to integrate the equations of motion defined in terms of the independent joint coordinates of the system. The reduction of the equations of motion to a minimal set...

A high-order boundary integral method for surface diffusions on elastically stressed axisymmetric rods

OpenAIRE

Li, Xiaofan; Nie, Qing

2009-01-01

Many applications in materials involve surface diffusion of elastically stressed solids. Study of singularity formation and long-time behavior of such solid surfaces requires accurate simulations in both space and time. Here we present a high-order boundary integral method for an elastically stressed solid with axi-symmetry due to surface diffusions. In this method, the boundary integrals for isotropic elasticity in axi-symmetric geometry are approximated through modified alternating quadratu...

Automatic numerical integration methods for Feynman integrals through 3-loop

International Nuclear Information System (INIS)

De Doncker, E; Olagbemi, O; Yuasa, F; Ishikawa, T; Kato, K

2015-01-01

We give numerical integration results for Feynman loop diagrams through 3-loop such as those covered by Laporta [1]. The methods are based on automatic adaptive integration, using iterated integration and extrapolation with programs from the QUADPACK package, or multivariate techniques from the ParInt package. The Dqags algorithm from QuadPack accommodates boundary singularities of fairly general types. PARINT is a package for multivariate integration layered over MPI (Message Passing Interface), which runs on clusters and incorporates advanced parallel/distributed techniques such as load balancing among processes that may be distributed over a network of nodes. Results are included for 3-loop self-energy diagrams without IR (infra-red) or UV (ultra-violet) singularities. A procedure based on iterated integration and extrapolation yields a novel method of numerical regularization for integrals with UV terms, and is applied to a set of 2-loop self-energy diagrams with UV singularities. (paper)

Issues in measure-preserving three dimensional flow integrators: Self-adjointness, reversibility, and non-uniform time stepping

International Nuclear Information System (INIS)

Finn, John M.

2015-01-01

Properties of integration schemes for solenoidal fields in three dimensions are studied, with a focus on integrating magnetic field lines in a plasma using adaptive time stepping. It is shown that implicit midpoint (IM) and a scheme we call three-dimensional leapfrog (LF) can do a good job (in the sense of preserving KAM tori) of integrating fields that are reversible, or (for LF) have a “special divergence-free” (SDF) property. We review the notion of a self-adjoint scheme, showing that such schemes are at least second order accurate and can always be formed by composing an arbitrary scheme with its adjoint. We also review the concept of reversibility, showing that a reversible but not exactly volume-preserving scheme can lead to a fractal invariant measure in a chaotic region, although this property may not often be observable. We also show numerical results indicating that the IM and LF schemes can fail to preserve KAM tori when the reversibility property (and the SDF property for LF) of the field is broken. We discuss extensions to measure preserving flows, the integration of magnetic field lines in a plasma and the integration of rays for several plasma waves. The main new result of this paper relates to non-uniform time stepping for volume-preserving flows. We investigate two potential schemes, both based on the general method of Feng and Shang [Numer. Math. 71, 451 (1995)], in which the flow is integrated in split time steps, each Hamiltonian in two dimensions. The first scheme is an extension of the method of extended phase space, a well-proven method of symplectic integration with non-uniform time steps. This method is found not to work, and an explanation is given. The second method investigated is a method based on transformation to canonical variables for the two split-step Hamiltonian systems. This method, which is related to the method of non-canonical generating functions of Richardson and Finn [Plasma Phys. Controlled Fusion 54, 014004 (2012

Assessing Backwards Integration as a Method of KBO Family Finding

Science.gov (United States)

Benfell, Nathan; Ragozzine, Darin

2018-04-01

The age of young asteroid collisional families can sometimes be determined by using backwards n-body integrations of the solar system. This method is not used for discovering young asteroid families and is limited by the unpredictable influence of the Yarkovsky effect on individual specific asteroids over time. Since these limitations are not as important for objects in the Kuiper belt, Marcus et al. 2011 suggested that backwards integration could be used to discover and characterize collisional families in the outer solar system. But various challenges present themselves when running precise and accurate 4+ Gyr integrations of Kuiper Belt objects. We have created simulated families of Kuiper Belt Objects with identical starting locations and velocity distributions, based on the Haumea Family. We then ran several long-term test integrations to observe the effect of various simulation parameters on integration results. These integrations were then used to investigate which parameters are of enough significance to require inclusion in the integration. Thereby we determined how to construct long-term integrations that both yield significant results and require manageable processing power. Additionally, we have tested the use of backwards integration as a method of discovery of potential young families in the Kuiper Belt.

Time-independent integral equation for Maxwell's system. Application of radar cross section computation

International Nuclear Information System (INIS)

Pujols, Agnes

1991-01-01

We prove that the scattering operator for the wave equation in the exterior of an non-homogeneous obstacle exists. Its distribution kernel is represented by a time-dependent boundary integral equation. A space-time integral variational formulation is developed for determining the current induced by the scattering of an electromagnetic wave by an homogeneous object. The discrete approximation of the variational problem using a finite element method in both space and time leads to stable convergent schemes, giving a numerical code for perfectly conducting cylinders. (author) [fr

Systematization of simplified J-integral evaluation method for flaw evaluation at high temperature

International Nuclear Information System (INIS)

Miura, Naoki; Takahashi, Yukio; Nakayama, Yasunari; Shimakawa, Takashi

2000-01-01

J-integral is an effective inelastic fracture parameter for the flaw evaluation of cracked components at high temperature. The evaluation of J-integral for an arbitrary crack configuration and an arbitrary loading condition can be generally accomplished by detailed numerical analysis such as finite element analysis, however, it is time-consuming and requires a high degree of expertise for its implementation. Therefore, it is important to develop simplified J-integral estimation techniques from the viewpoint of industrial requirements. In this study, a simplified J-integral evaluation method is proposed to estimate two types of J-integral parameters. One is the fatigue J-integral range to describe fatigue crack propagation behavior, and the other is the creep J-integral to describe creep crack propagation behavior. This paper presents the systematization of the simplified J-integral evaluation method incorporated with the reference stress method and the concept of elastic follow-up, and proposes a comprehensive evaluation procedure. The verification of the proposed method is presented in Part II of this paper. (author)

A comparison of response spectrum and direct integration analysis methods as applied to a nuclear component support structure

International Nuclear Information System (INIS)

Bryan, B.J.; Flanders, H.E. Jr.

1992-01-01

Seismic qualification of Class I nuclear components is accomplished using a variety of analytical methods. This paper compares the results of time history dynamic analyses of a heat exchanger support structure using response spectrum and time history direct integration analysis methods. Dynamic analysis is performed on the detailed component models using the two methods. A nonlinear elastic model is used for both the response spectrum and direct integration methods. A nonlinear model which includes friction and nonlinear springs, is analyzed using time history input by direct integration. The loads from the three cases are compared

Efficient Simulation of Compressible, Viscous Fluids using Multi-rate Time Integration

Science.gov (United States)

Mikida, Cory; Kloeckner, Andreas; Bodony, Daniel

2017-11-01

In the numerical simulation of problems of compressible, viscous fluids with single-rate time integrators, the global timestep used is limited to that of the finest mesh point or fastest physical process. This talk discusses the application of multi-rate Adams-Bashforth (MRAB) integrators to an overset mesh framework to solve compressible viscous fluid problems of varying scale with improved efficiency, with emphasis on the strategy of timescale separation and the application of the resulting numerical method to two sample problems: subsonic viscous flow over a cylinder and a viscous jet in crossflow. The results presented indicate the numerical efficacy of MRAB integrators, outline a number of outstanding code challenges, demonstrate the expected reduction in time enabled by MRAB, and emphasize the need for proper load balancing through spatial decomposition in order for parallel runs to achieve the predicted time-saving benefit. This material is based in part upon work supported by the Department of Energy, National Nuclear Security Administration, under Award Number DE-NA0002374.

Integral Equation Methods for Electromagnetic and Elastic Waves

CERN Document Server

Chew, Weng; Hu, Bin

2008-01-01

Integral Equation Methods for Electromagnetic and Elastic Waves is an outgrowth of several years of work. There have been no recent books on integral equation methods. There are books written on integral equations, but either they have been around for a while, or they were written by mathematicians. Much of the knowledge in integral equation methods still resides in journal papers. With this book, important relevant knowledge for integral equations are consolidated in one place and researchers need only read the pertinent chapters in this book to gain important knowledge needed for integral eq

Parallel Jacobi EVD Methods on Integrated Circuits

Directory of Open Access Journals (Sweden)

Chi-Chia Sun

2014-01-01

Full Text Available Design strategies for parallel iterative algorithms are presented. In order to further study different tradeoff strategies in design criteria for integrated circuits, A 10 × 10 Jacobi Brent-Luk-EVD array with the simplified μ-CORDIC processor is used as an example. The experimental results show that using the μ-CORDIC processor is beneficial for the design criteria as it yields a smaller area, faster overall computation time, and less energy consumption than the regular CORDIC processor. It is worth to notice that the proposed parallel EVD method can be applied to real-time and low-power array signal processing algorithms performing beamforming or DOA estimation.

Time-Dependent Heat Conduction Problems Solved by an Integral-Equation Approach

International Nuclear Information System (INIS)

Oberaigner, E.R.; Leindl, M.; Antretter, T.

2010-01-01

Full text: A classical task of mathematical physics is the formulation and solution of a time dependent thermoelastic problem. In this work we develop an algorithm for solving the time-dependent heat conduction equation c p ρ∂ t T-kT, ii =0 in an analytical, exact fashion for a two-component domain. By the Green's function approach the formal solution of the problem is obtained. As an intermediate result an integral-equation for the temperature history at the domain interface is formulated which can be solved analytically. This method is applied to a classical engineering problem, i.e. to a special case of a Stefan-Problem. The Green's function approach in conjunction with the integral-equation method is very useful in cases were strong discontinuities or jumps occur. The initial conditions and the system parameters of the investigated problem give rise to two jumps in the temperature field. Purely numerical solutions are obtained by using the FEM (finite element method) and the FDM (finite difference method) and compared with the analytical approach. At the domain boundary the analytical solution and the FEM-solution are in good agreement, but the FDM results show a signicant smearing effect. (author)

Subdomain Precise Integration Method for Periodic Structures

Directory of Open Access Journals (Sweden)

F. Wu

2014-01-01

Full Text Available A subdomain precise integration method is developed for the dynamical responses of periodic structures comprising many identical structural cells. The proposed method is based on the precise integration method, the subdomain scheme, and the repeatability of the periodic structures. In the proposed method, each structural cell is seen as a super element that is solved using the precise integration method, considering the repeatability of the structural cells. The computational efforts and the memory size of the proposed method are reduced, while high computational accuracy is achieved. Therefore, the proposed method is particularly suitable to solve the dynamical responses of periodic structures. Two numerical examples are presented to demonstrate the accuracy and efficiency of the proposed method through comparison with the Newmark and Runge-Kutta methods.

A Vision/Inertia Integrated Positioning Method Using Position and Orientation Matching

Directory of Open Access Journals (Sweden)

Xiaoyue Zhang

2017-01-01

Full Text Available A vision/inertia integrated positioning method using position and orientation matching which can be adopted on intelligent vehicle such as automated guided vehicle (AGV and mobile robot is proposed in this work. The method is introduced firstly. Landmarks are placed into the navigation field and camera and inertial measurement unit (IMU are installed on the vehicle. Vision processor calculates the azimuth and position information from the pictures which include artificial landmarks with the known direction and position. Inertial navigation system (INS calculates the azimuth and position of vehicle in real time and the calculated pixel position of landmark can be computed from the INS output position. Then the needed mathematical models are established and integrated navigation is implemented by Kalman filter with the observation of azimuth and the calculated pixel position of landmark. Navigation errors and IMU errors are estimated and compensated in real time so that high precision navigation results can be got. Finally, simulation and test are performed, respectively. Both simulation and test results prove that this vision/inertia integrated positioning method using position and orientation matching has feasibility and it can achieve centimeter-level autonomic continuous navigation.

Evolution equation of Lie-type for finite deformations, time-discrete integration, and incremental methods

Czech Academy of Sciences Publication Activity Database

Fiala, Zdeněk

2015-01-01

Roč. 226, č. 1 (2015), s. 17-35 ISSN 0001-5970 R&D Projects: GA ČR(CZ) GA103/09/2101 Institutional support: RVO:68378297 Keywords : solid mechanics * finite deformations * evolution equation of Lie-type * time-discrete integration Subject RIV: BA - General Mathematics OBOR OECD: Statistics and probability Impact factor: 1.694, year: 2015 http://link.springer.com/article/10.1007%2Fs00707-014-1162-9#page-1

Reliable Viscosity Calculation from Equilibrium Molecular Dynamics Simulations: A Time Decomposition Method.

Science.gov (United States)

Zhang, Yong; Otani, Akihito; Maginn, Edward J

2015-08-11

Equilibrium molecular dynamics is often used in conjunction with a Green-Kubo integral of the pressure tensor autocorrelation function to compute the shear viscosity of fluids. This approach is computationally expensive and is subject to a large amount of variability because the plateau region of the Green-Kubo integral is difficult to identify unambiguously. Here, we propose a time decomposition approach for computing the shear viscosity using the Green-Kubo formalism. Instead of one long trajectory, multiple independent trajectories are run and the Green-Kubo relation is applied to each trajectory. The averaged running integral as a function of time is fit to a double-exponential function with a weighting function derived from the standard deviation of the running integrals. Such a weighting function minimizes the uncertainty of the estimated shear viscosity and provides an objective means of estimating the viscosity. While the formal Green-Kubo integral requires an integration to infinite time, we suggest an integration cutoff time tcut, which can be determined by the relative values of the running integral and the corresponding standard deviation. This approach for computing the shear viscosity can be easily automated and used in computational screening studies where human judgment and intervention in the data analysis are impractical. The method has been applied to the calculation of the shear viscosity of a relatively low-viscosity liquid, ethanol, and relatively high-viscosity ionic liquid, 1-n-butyl-3-methylimidazolium bis(trifluoromethane-sulfonyl)imide ([BMIM][Tf2N]), over a range of temperatures. These test cases show that the method is robust and yields reproducible and reliable shear viscosity values.

Method and apparatus for in-system redundant array repair on integrated circuits

Science.gov (United States)

Bright, Arthur A.; Crumley, Paul G.; Dombrowa, Marc B.; Douskey, Steven M.; Haring, Rudolf A.; Oakland, Steven F.; Ouellette, Michael R.; Strissel, Scott A.

2007-12-18

Disclosed is a method of repairing an integrated circuit of the type comprising of a multitude of memory arrays and a fuse box holding control data for controlling redundancy logic of the arrays. The method comprises the steps of providing the integrated circuit with a control data selector for passing the control data from the fuse box to the memory arrays; providing a source of alternate control data, external of the integrated circuit; and connecting the source of alternate control data to the control data selector. The method comprises the further step of, at a given time, passing the alternate control data from the source thereof, through the control data selector and to the memory arrays to control the redundancy logic of the memory arrays.

Method and apparatus for in-system redundant array repair on integrated circuits

Science.gov (United States)

Bright, Arthur A [Croton-on-Hudson, NY; Crumley, Paul G [Yorktown Heights, NY; Dombrowa, Marc B [Bronx, NY; Douskey, Steven M [Rochester, MN; Haring, Rudolf A [Cortlandt Manor, NY; Oakland, Steven F [Colchester, VT; Ouellette, Michael R [Westford, VT; Strissel, Scott A [Byron, MN

2008-07-29

Disclosed is a method of repairing an integrated circuit of the type comprising of a multitude of memory arrays and a fuse box holding control data for controlling redundancy logic of the arrays. The method comprises the steps of providing the integrated circuit with a control data selector for passing the control data from the fuse box to the memory arrays; providing a source of alternate control data, external of the integrated circuit; and connecting the source of alternate control data to the control data selector. The method comprises the further step of, at a given time, passing the alternate control data from the source thereof, through the control data selector and to the memory arrays to control the redundancy logic of the memory arrays.

Integration of the time-dependent heat equation in the fuel rod performance program IAMBUS

International Nuclear Information System (INIS)

West, G.

1982-01-01

An iterative numerical method for integration of the time-dependent heat equation is described. No presuppositions are made for the dependency of the thermal conductivity and heat capacity on space, time and temperature. (orig.) [de

Parameters affecting temporal resolution of Time Resolved Integrative Optical Neutron Detector (TRION)

International Nuclear Information System (INIS)

Mor, I; Vartsky, D; Bar, D; Feldman, G; Goldberg, M B; Brandis, M; Dangendorf, V; Tittelmeier, K; Bromberger, B; Weierganz, M

2013-01-01

The Time-Resolved Integrative Optical Neutron (TRION) detector was developed for Fast Neutron Resonance Radiography (FNRR), a fast-neutron transmission imaging method that exploits characteristic energy-variations of the total scattering cross-section in the E n = 1–10 MeV range to detect specific elements within a radiographed object. As opposed to classical event-counting time of flight (ECTOF), it integrates the detector signal during a well-defined neutron Time of Flight window corresponding to a pre-selected energy bin, e.g., the energy-interval spanning a cross-section resonance of an element such as C, O and N. The integrative characteristic of the detector permits loss-free operation at very intense, pulsed neutron fluxes, at a cost however, of recorded temporal resolution degradation This work presents a theoretical and experimental evaluation of detector related parameters which affect temporal resolution of the TRION system

An integrative time-varying frequency detection and channel sounding method for dynamic plasma sheath

Science.gov (United States)

Shi, Lei; Yao, Bo; Zhao, Lei; Liu, Xiaotong; Yang, Min; Liu, Yanming

2018-01-01

The plasma sheath-surrounded hypersonic vehicle is a dynamic and time-varying medium and it is almost impossible to calculate time-varying physical parameters directly. The in-fight detection of the time-varying degree is important to understand the dynamic nature of the physical parameters and their effect on re-entry communication. In this paper, a constant envelope zero autocorrelation (CAZAC) sequence based on time-varying frequency detection and channel sounding method is proposed to detect the plasma sheath electronic density time-varying property and wireless channel characteristic. The proposed method utilizes the CAZAC sequence, which has excellent autocorrelation and spread gain characteristics, to realize dynamic time-varying detection/channel sounding under low signal-to-noise ratio in the plasma sheath environment. Theoretical simulation under a typical time-varying radio channel shows that the proposed method is capable of detecting time-variation frequency up to 200 kHz and can trace the channel amplitude and phase in the time domain well under -10 dB. Experimental results conducted in the RF modulation discharge plasma device verified the time variation detection ability in practical dynamic plasma sheath. Meanwhile, nonlinear phenomenon of dynamic plasma sheath on communication signal is observed thorough channel sounding result.

Numerical methods for engine-airframe integration

International Nuclear Information System (INIS)

Murthy, S.N.B.; Paynter, G.C.

1986-01-01

Various papers on numerical methods for engine-airframe integration are presented. The individual topics considered include: scientific computing environment for the 1980s, overview of prediction of complex turbulent flows, numerical solutions of the compressible Navier-Stokes equations, elements of computational engine/airframe integrations, computational requirements for efficient engine installation, application of CAE and CFD techniques to complete tactical missile design, CFD applications to engine/airframe integration, and application of a second-generation low-order panel methods to powerplant installation studies. Also addressed are: three-dimensional flow analysis of turboprop inlet and nacelle configurations, application of computational methods to the design of large turbofan engine nacelles, comparison of full potential and Euler solution algorithms for aeropropulsive flow field computations, subsonic/transonic, supersonic nozzle flows and nozzle integration, subsonic/transonic prediction capabilities for nozzle/afterbody configurations, three-dimensional viscous design methodology of supersonic inlet systems for advanced technology aircraft, and a user's technology assessment

Perfectly matched layer for the time domain finite element method

International Nuclear Information System (INIS)

Rylander, Thomas; Jin Jianming

2004-01-01

A new perfectly matched layer (PML) formulation for the time domain finite element method is described and tested for Maxwell's equations. In particular, we focus on the time integration scheme which is based on Galerkin's method with a temporally piecewise linear expansion of the electric field. The time stepping scheme is constructed by forming a linear combination of exact and trapezoidal integration applied to the temporal weak form, which reduces to the well-known Newmark scheme in the case without PML. Extensive numerical tests on scattering from infinitely long metal cylinders in two dimensions show good accuracy and no signs of instabilities. For a circular cylinder, the proposed scheme indicates the expected second order convergence toward the analytic solution and gives less than 2% root-mean-square error in the bistatic radar cross section (RCS) for resolutions with more than 10 points per wavelength. An ogival cylinder, which has sharp corners supporting field singularities, shows similar accuracy in the monostatic RCS

Reparametrization in the path integral over finite dimensional manifold with a time-dependent metric

International Nuclear Information System (INIS)

Storchak, S.N.

1988-01-01

The path reparametrization procedure in the path integral is considered using the methods of stochastic processes for diffusion on finite dimensional manifold with a time-dependent metric. the reparametrization Jacobian has been obtained. The formulas of reparametrization for a symbolic presentation of the path integral have been derived

A numerical method for resonance integral calculations

International Nuclear Information System (INIS)

Tanbay, Tayfun; Ozgener, Bilge

2013-01-01

A numerical method has been proposed for resonance integral calculations and a cubic fit based on least squares approximation to compute the optimum Bell factor is given. The numerical method is based on the discretization of the neutron slowing down equation. The scattering integral is approximated by taking into account the location of the upper limit in energy domain. The accuracy of the method has been tested by performing computations of resonance integrals for uranium dioxide isolated rods and comparing the results with empirical values. (orig.)

Review of Statistical Learning Methods in Integrated Omics Studies (An Integrated Information Science).

Science.gov (United States)

Zeng, Irene Sui Lan; Lumley, Thomas

2018-01-01

Integrated omics is becoming a new channel for investigating the complex molecular system in modern biological science and sets a foundation for systematic learning for precision medicine. The statistical/machine learning methods that have emerged in the past decade for integrated omics are not only innovative but also multidisciplinary with integrated knowledge in biology, medicine, statistics, machine learning, and artificial intelligence. Here, we review the nontrivial classes of learning methods from the statistical aspects and streamline these learning methods within the statistical learning framework. The intriguing findings from the review are that the methods used are generalizable to other disciplines with complex systematic structure, and the integrated omics is part of an integrated information science which has collated and integrated different types of information for inferences and decision making. We review the statistical learning methods of exploratory and supervised learning from 42 publications. We also discuss the strengths and limitations of the extended principal component analysis, cluster analysis, network analysis, and regression methods. Statistical techniques such as penalization for sparsity induction when there are fewer observations than the number of features and using Bayesian approach when there are prior knowledge to be integrated are also included in the commentary. For the completeness of the review, a table of currently available software and packages from 23 publications for omics are summarized in the appendix.

The general 2-D moments via integral transform method for acoustic radiation and scattering

Science.gov (United States)

Smith, Jerry R.; Mirotznik, Mark S.

2004-05-01

The moments via integral transform method (MITM) is a technique to analytically reduce the 2-D method of moments (MoM) impedance double integrals into single integrals. By using a special integral representation of the Green's function, the impedance integral can be analytically simplified to a single integral in terms of transformed shape and weight functions. The reduced expression requires fewer computations and reduces the fill times of the MoM impedance matrix. Furthermore, the resulting integral is analytic for nearly arbitrary shape and weight function sets. The MITM technique is developed for mixed boundary conditions and predictions with basic shape and weight function sets are presented. Comparisons of accuracy and speed between MITM and brute force are presented. [Work sponsored by ONR and NSWCCD ILIR Board.

Molecular radiotherapy: The NUKFIT software for calculating the time-integrated activity coefficient

Energy Technology Data Exchange (ETDEWEB)

Kletting, P.; Schimmel, S.; Luster, M. [Klinik für Nuklearmedizin, Universität Ulm, Ulm 89081 (Germany); Kestler, H. A. [Research Group Bioinformatics and Systems Biology, Institut für Neuroinformatik, Universität Ulm, Ulm 89081 (Germany); Hänscheid, H.; Fernández, M.; Lassmann, M. [Klinik für Nuklearmedizin, Universität Würzburg, Würzburg 97080 (Germany); Bröer, J. H.; Nosske, D. [Bundesamt für Strahlenschutz, Fachbereich Strahlenschutz und Gesundheit, Oberschleißheim 85764 (Germany); Glatting, G. [Medical Radiation Physics/Radiation Protection, Medical Faculty Mannheim, Heidelberg University, Mannheim 68167 (Germany)

2013-10-15

Purpose: Calculation of the time-integrated activity coefficient (residence time) is a crucial step in dosimetry for molecular radiotherapy. However, available software is deficient in that it is either not tailored for the use in molecular radiotherapy and/or does not include all required estimation methods. The aim of this work was therefore the development and programming of an algorithm which allows for an objective and reproducible determination of the time-integrated activity coefficient and its standard error.Methods: The algorithm includes the selection of a set of fitting functions from predefined sums of exponentials and the choice of an error model for the used data. To estimate the values of the adjustable parameters an objective function, depending on the data, the parameters of the error model, the fitting function and (if required and available) Bayesian information, is minimized. To increase reproducibility and user-friendliness the starting values are automatically determined using a combination of curve stripping and random search. Visual inspection, the coefficient of determination, the standard error of the fitted parameters, and the correlation matrix are provided to evaluate the quality of the fit. The functions which are most supported by the data are determined using the corrected Akaike information criterion. The time-integrated activity coefficient is estimated by analytically integrating the fitted functions. Its standard error is determined assuming Gaussian error propagation. The software was implemented using MATLAB.Results: To validate the proper implementation of the objective function and the fit functions, the results of NUKFIT and SAAM numerical, a commercially available software tool, were compared. The automatic search for starting values was successfully tested for reproducibility. The quality criteria applied in conjunction with the Akaike information criterion allowed the selection of suitable functions. Function fit

Mining method selection by integrated AHP and PROMETHEE method.

Science.gov (United States)

Bogdanovic, Dejan; Nikolic, Djordje; Ilic, Ivana

2012-03-01

Selecting the best mining method among many alternatives is a multicriteria decision making problem. The aim of this paper is to demonstrate the implementation of an integrated approach that employs AHP and PROMETHEE together for selecting the most suitable mining method for the "Coka Marin" underground mine in Serbia. The related problem includes five possible mining methods and eleven criteria to evaluate them. Criteria are accurately chosen in order to cover the most important parameters that impact on the mining method selection, such as geological and geotechnical properties, economic parameters and geographical factors. The AHP is used to analyze the structure of the mining method selection problem and to determine weights of the criteria, and PROMETHEE method is used to obtain the final ranking and to make a sensitivity analysis by changing the weights. The results have shown that the proposed integrated method can be successfully used in solving mining engineering problems.

Integrated circuit and method of arbitration in a network on an integrated circuit.

NARCIS (Netherlands)

2011-01-01

The invention relates to an integrated circuit and to a method of arbitration in a network on an integrated circuit. According to the invention, a method of arbitration in a network on an integrated circuit is provided, the network comprising a router unit, the router unit comprising a first input

An integration weighting method to evaluate extremum coordinates

International Nuclear Information System (INIS)

Ilyushchenko, V.I.

1990-01-01

The numerical version of the Laplace asymptotics has been used to evaluate the coordinates of extrema of multivariate continuous and discontinuous test functions. The performed computer experiments demonstrate the high efficiency of the integration method proposed. The saturating dependence of extremum coordinates on such parameters as a number of integration subregions and that of K going /theoretically/ to infinity has been studied in detail for the limitand being a ratio of two Laplace integrals with exponentiated K. The given method is an integral equivalent of that of weighted means. As opposed to the standard optimization methods of the zero, first and second order the proposed method can be successfully applied to optimize discontinuous objective functions, too. There are possibilities of applying the integration method in the cases, when the conventional techniques fail due to poor analytical properties of the objective functions near extremal points. The proposed method is efficient in searching for both local and global extrema of multimodal objective functions. 12 refs.; 4 tabs

Enhanced IMC based PID controller design for non-minimum phase (NMP) integrating processes with time delays.

Science.gov (United States)

Ghousiya Begum, K; Seshagiri Rao, A; Radhakrishnan, T K

2017-05-01

Internal model control (IMC) with optimal H 2 minimization framework is proposed in this paper for design of proportional-integral-derivative (PID) controllers. The controller design is addressed for integrating and double integrating time delay processes with right half plane (RHP) zeros. Blaschke product is used to derive the optimal controller. There is a single adjustable closed loop tuning parameter for controller design. Systematic guidelines are provided for selection of this tuning parameter based on maximum sensitivity. Simulation studies have been carried out on various integrating time delay processes to show the advantages of the proposed method. The proposed controller provides enhanced closed loop performances when compared to recently reported methods in the literature. Quantitative comparative analysis has been carried out using the performance indices, Integral Absolute Error (IAE) and Total Variation (TV). Copyright © 2017 ISA. Published by Elsevier Ltd. All rights reserved.

Hybrid perturbation methods based on statistical time series models

Science.gov (United States)

San-Juan, Juan Félix; San-Martín, Montserrat; Pérez, Iván; López, Rosario

2016-04-01

In this work we present a new methodology for orbit propagation, the hybrid perturbation theory, based on the combination of an integration method and a prediction technique. The former, which can be a numerical, analytical or semianalytical theory, generates an initial approximation that contains some inaccuracies derived from the fact that, in order to simplify the expressions and subsequent computations, not all the involved forces are taken into account and only low-order terms are considered, not to mention the fact that mathematical models of perturbations not always reproduce physical phenomena with absolute precision. The prediction technique, which can be based on either statistical time series models or computational intelligence methods, is aimed at modelling and reproducing missing dynamics in the previously integrated approximation. This combination results in the precision improvement of conventional numerical, analytical and semianalytical theories for determining the position and velocity of any artificial satellite or space debris object. In order to validate this methodology, we present a family of three hybrid orbit propagators formed by the combination of three different orders of approximation of an analytical theory and a statistical time series model, and analyse their capability to process the effect produced by the flattening of the Earth. The three considered analytical components are the integration of the Kepler problem, a first-order and a second-order analytical theories, whereas the prediction technique is the same in the three cases, namely an additive Holt-Winters method.

Integral equation approach to time-dependent kinematic dynamos in finite domains

International Nuclear Information System (INIS)

Xu Mingtian; Stefani, Frank; Gerbeth, Gunter

2004-01-01

The homogeneous dynamo effect is at the root of cosmic magnetic field generation. With only a very few exceptions, the numerical treatment of homogeneous dynamos is carried out in the framework of the differential equation approach. The present paper tries to facilitate the use of integral equations in dynamo research. Apart from the pedagogical value to illustrate dynamo action within the well-known picture of the Biot-Savart law, the integral equation approach has a number of practical advantages. The first advantage is its proven numerical robustness and stability. The second and perhaps most important advantage is its applicability to dynamos in arbitrary geometries. The third advantage is its intimate connection to inverse problems relevant not only for dynamos but also for technical applications of magnetohydrodynamics. The paper provides the first general formulation and application of the integral equation approach to time-dependent kinematic dynamos, with stationary dynamo sources, in finite domains. The time dependence is restricted to the magnetic field, whereas the velocity or corresponding mean-field sources of dynamo action are supposed to be stationary. For the spherically symmetric α 2 dynamo model it is shown how the general formulation is reduced to a coupled system of two radial integral equations for the defining scalars of the poloidal and toroidal field components. The integral equation formulation for spherical dynamos with general stationary velocity fields is also derived. Two numerical examples - the α 2 dynamo model with radially varying α and the Bullard-Gellman model - illustrate the equivalence of the approach with the usual differential equation method. The main advantage of the method is exemplified by the treatment of an α 2 dynamo in rectangular domains

Integrated Detection and Prediction of Influenza Activity for Real-Time Surveillance: Algorithm Design.

Science.gov (United States)

Spreco, Armin; Eriksson, Olle; Dahlström, Örjan; Cowling, Benjamin John; Timpka, Toomas

2017-06-15

Influenza is a viral respiratory disease capable of causing epidemics that represent a threat to communities worldwide. The rapidly growing availability of electronic "big data" from diagnostic and prediagnostic sources in health care and public health settings permits advance of a new generation of methods for local detection and prediction of winter influenza seasons and influenza pandemics. The aim of this study was to present a method for integrated detection and prediction of influenza virus activity in local settings using electronically available surveillance data and to evaluate its performance by retrospective application on authentic data from a Swedish county. An integrated detection and prediction method was formally defined based on a design rationale for influenza detection and prediction methods adapted for local surveillance. The novel method was retrospectively applied on data from the winter influenza season 2008-09 in a Swedish county (population 445,000). Outcome data represented individuals who met a clinical case definition for influenza (based on International Classification of Diseases version 10 [ICD-10] codes) from an electronic health data repository. Information from calls to a telenursing service in the county was used as syndromic data source. The novel integrated detection and prediction method is based on nonmechanistic statistical models and is designed for integration in local health information systems. The method is divided into separate modules for detection and prediction of local influenza virus activity. The function of the detection module is to alert for an upcoming period of increased load of influenza cases on local health care (using influenza-diagnosis data), whereas the function of the prediction module is to predict the timing of the activity peak (using syndromic data) and its intensity (using influenza-diagnosis data). For detection modeling, exponential regression was used based on the assumption that the beginning

Symmetric and arbitrarily high-order Birkhoff-Hermite time integrators and their long-time behaviour for solving nonlinear Klein-Gordon equations

Science.gov (United States)

Liu, Changying; Iserles, Arieh; Wu, Xinyuan

2018-03-01

The Klein-Gordon equation with nonlinear potential occurs in a wide range of application areas in science and engineering. Its computation represents a major challenge. The main theme of this paper is the construction of symmetric and arbitrarily high-order time integrators for the nonlinear Klein-Gordon equation by integrating Birkhoff-Hermite interpolation polynomials. To this end, under the assumption of periodic boundary conditions, we begin with the formulation of the nonlinear Klein-Gordon equation as an abstract second-order ordinary differential equation (ODE) and its operator-variation-of-constants formula. We then derive a symmetric and arbitrarily high-order Birkhoff-Hermite time integration formula for the nonlinear abstract ODE. Accordingly, the stability, convergence and long-time behaviour are rigorously analysed once the spatial differential operator is approximated by an appropriate positive semi-definite matrix, subject to suitable temporal and spatial smoothness. A remarkable characteristic of this new approach is that the requirement of temporal smoothness is reduced compared with the traditional numerical methods for PDEs in the literature. Numerical results demonstrate the advantage and efficiency of our time integrators in comparison with the existing numerical approaches.

Selective Integration in the Material-Point Method

DEFF Research Database (Denmark)

Andersen, Lars; Andersen, Søren; Damkilde, Lars

2009-01-01

The paper deals with stress integration in the material-point method. In order to avoid parasitic shear in bending, a formulation is proposed, based on selective integration in the background grid that is used to solve the governing equations. The suggested integration scheme is compared...... to a traditional material-point-method computation in which the stresses are evaluated at the material points. The deformation of a cantilever beam is analysed, assuming elastic or elastoplastic material behaviour....

Integration of Real-Time Data Into Building Automation Systems

Energy Technology Data Exchange (ETDEWEB)

Mark J. Stunder; Perry Sebastian; Brenda A. Chube; Michael D. Koontz

2003-04-16

The project goal was to investigate the possibility of using predictive real-time information from the Internet as an input to building management system algorithms. The objectives were to identify the types of information most valuable to commercial and residential building owners, managers, and system designers. To comprehensively investigate and document currently available electronic real-time information suitable for use in building management systems. Verify the reliability of the information and recommend accreditation methods for data and providers. Assess methodologies to automatically retrieve and utilize the information. Characterize equipment required to implement automated integration. Demonstrate the feasibility and benefits of using the information in building management systems. Identify evolutionary control strategies.

Integration of MDSplus in real-time systems

International Nuclear Information System (INIS)

Luchetta, A.; Manduchi, G.; Taliercio, C.

2006-01-01

RFX-mod makes extensive usage of real-time systems for feedback control and uses MDSplus to interface them to the main Data Acquisition system. For this purpose, the core of MDSplus has been ported to VxWorks, the operating system used for real-time control in RFX. Using this approach, it is possible to integrate real-time systems, but MDSplus is used only for non-real-time tasks, i.e. those tasks which are executed before and after the pulse and whose performance does not affect the system time constraints. More extensive use of MDSplus in real-time systems is foreseen, and a real-time layer for MDSplus is under development, which will provide access to memory-mapped pulse files, shared by the tasks running on the same CPU. Real-time communication will also be integrated in the MDSplus core to provide support for distributed memory-mapped pulse files

Multi-channel time-division integrator in HL-2A

International Nuclear Information System (INIS)

Yan Ji

2008-01-01

HL-2A is China's first Tokamak device with divertor configuration (magnetic confinement controlled nuclear fusion device). To find out the details of on-going fusion reaction at different times is of important significance in achieving controlled nuclear fusion. We developed a new type multi-channel time-division integrator for HL-2A. It has functions of automatic cutting off negative pulse of the input signals, optional integrating time division spacing 0.2-1 ms, TTL starting trigger signal, automatic regularly work 20 s, and integrating 10 channel at the same time. (authors)

The time-dependent density matrix renormalisation group method

Science.gov (United States)

Ma, Haibo; Luo, Zhen; Yao, Yao

2018-04-01

Substantial progress of the time-dependent density matrix renormalisation group (t-DMRG) method in the recent 15 years is reviewed in this paper. By integrating the time evolution with the sweep procedures in density matrix renormalisation group (DMRG), t-DMRG provides an efficient tool for real-time simulations of the quantum dynamics for one-dimensional (1D) or quasi-1D strongly correlated systems with a large number of degrees of freedom. In the illustrative applications, the t-DMRG approach is applied to investigate the nonadiabatic processes in realistic chemical systems, including exciton dissociation and triplet fission in polymers and molecular aggregates as well as internal conversion in pyrazine molecule.

Integration of RNA-Seq and RPPA data for survival time prediction in cancer patients.

Science.gov (United States)

Isik, Zerrin; Ercan, Muserref Ece

2017-10-01

Integration of several types of patient data in a computational framework can accelerate the identification of more reliable biomarkers, especially for prognostic purposes. This study aims to identify biomarkers that can successfully predict the potential survival time of a cancer patient by integrating the transcriptomic (RNA-Seq), proteomic (RPPA), and protein-protein interaction (PPI) data. The proposed method -RPBioNet- employs a random walk-based algorithm that works on a PPI network to identify a limited number of protein biomarkers. Later, the method uses gene expression measurements of the selected biomarkers to train a classifier for the survival time prediction of patients. RPBioNet was applied to classify kidney renal clear cell carcinoma (KIRC), glioblastoma multiforme (GBM), and lung squamous cell carcinoma (LUSC) patients based on their survival time classes (long- or short-term). The RPBioNet method correctly identified the survival time classes of patients with between 66% and 78% average accuracy for three data sets. RPBioNet operates with only 20 to 50 biomarkers and can achieve on average 6% higher accuracy compared to the closest alternative method, which uses only RNA-Seq data in the biomarker selection. Further analysis of the most predictive biomarkers highlighted genes that are common for both cancer types, as they may be driver proteins responsible for cancer progression. The novelty of this study is the integration of a PPI network with mRNA and protein expression data to identify more accurate prognostic biomarkers that can be used for clinical purposes in the future. Copyright © 2017 Elsevier Ltd. All rights reserved.

High-speed railway real-time localization auxiliary method based on deep neural network

Science.gov (United States)

Chen, Dongjie; Zhang, Wensheng; Yang, Yang

2017-11-01

High-speed railway intelligent monitoring and management system is composed of schedule integration, geographic information, location services, and data mining technology for integration of time and space data. Assistant localization is a significant submodule of the intelligent monitoring system. In practical application, the general access is to capture the image sequences of the components by using a high-definition camera, digital image processing technique and target detection, tracking and even behavior analysis method. In this paper, we present an end-to-end character recognition method based on a deep CNN network called YOLO-toc for high-speed railway pillar plate number. Different from other deep CNNs, YOLO-toc is an end-to-end multi-target detection framework, furthermore, it exhibits a state-of-art performance on real-time detection with a nearly 50fps achieved on GPU (GTX960). Finally, we realize a real-time but high-accuracy pillar plate number recognition system and integrate natural scene OCR into a dedicated classification YOLO-toc model.

Simulation electromagnetic scattering on bodies through integral equation and neural networks methods

Science.gov (United States)

Lvovich, I. Ya; Preobrazhenskiy, A. P.; Choporov, O. N.

2018-05-01

The paper deals with the issue of electromagnetic scattering on a perfectly conducting diffractive body of a complex shape. Performance calculation of the body scattering is carried out through the integral equation method. Fredholm equation of the second time was used for calculating electric current density. While solving the integral equation through the moments method, the authors have properly described the core singularity. The authors determined piecewise constant functions as basic functions. The chosen equation was solved through the moments method. Within the Kirchhoff integral approach it is possible to define the scattered electromagnetic field, in some way related to obtained electrical currents. The observation angles sector belongs to the area of the front hemisphere of the diffractive body. To improve characteristics of the diffractive body, the authors used a neural network. All the neurons contained a logsigmoid activation function and weighted sums as discriminant functions. The paper presents the matrix of weighting factors of the connectionist model, as well as the results of the optimized dimensions of the diffractive body. The paper also presents some basic steps in calculation technique of the diffractive bodies, based on the combination of integral equation and neural networks methods.

Hippocampal “Time Cells”: Time versus Path Integration

Science.gov (United States)

Kraus, Benjamin J.; Robinson, Robert J.; White, John A.; Eichenbaum, Howard; Hasselmo, Michael E.

2014-01-01

SUMMARY Recent studies have reported the existence of hippocampal “time cells,” neurons that fire at particular moments during periods when behavior and location are relatively constant. However, an alternative explanation of apparent time coding is that hippocampal neurons “path integrate” to encode the distance an animal has traveled. Here, we examined hippocampal neuronal firing patterns as rats ran in place on a treadmill, thus “clamping” behavior and location, while we varied the treadmill speed to distinguish time elapsed from distance traveled. Hippocampal neurons were strongly influenced by time and distance, and less so by minor variations in location. Furthermore, the activity of different neurons reflected integration over time and distance to varying extents, with most neurons strongly influenced by both factors and some significantly influenced by only time or distance. Thus, hippocampal neuronal networks captured both the organization of time and distance in a situation where these dimensions dominated an ongoing experience. PMID:23707613

Numerical functional integration method for studying the properties of the physical vacuum

International Nuclear Information System (INIS)

Lobanov, Yu.Yu.

1998-01-01

The new approach to investigate the physical vacuum in quantum theories including its nonperturbative topological structure is discussed. This approach is based on the representation of the matrix element of the evolution operator in Euclidean metrics in a form of the functional integral with a certain measure in the corresponding space and on the use of approximation formulas which we constructed for this kind of integral. No preliminary discretization of space and time is required, as well as no simplifying assumptions like semiclassical approximation, collective excitations, introduction of ''short-time'' propagators, etc. are necessary in this approach. The method allows to use the more preferable deterministic algorithms instead of the traditional stochastic technique. It has been proven that our approach has important advantages over the other known methods, including the higher efficiency of computations. Examples of application of the method to the numerical study of some potential nuclear models and to the computation of the topological susceptibility and the θ-vacua energy are presented. (author)

GENERIC Integrators: Structure Preserving Time Integration for Thermodynamic Systems

Science.gov (United States)

Öttinger, Hans Christian

2018-04-01

Thermodynamically admissible evolution equations for non-equilibrium systems are known to possess a distinct mathematical structure. Within the GENERIC (general equation for the non-equilibrium reversible-irreversible coupling) framework of non-equilibrium thermodynamics, which is based on continuous time evolution, we investigate the possibility of preserving all the structural elements in time-discretized equations. Our approach, which follows Moser's [1] construction of symplectic integrators for Hamiltonian systems, is illustrated for the damped harmonic oscillator. Alternative approaches are sketched.

Numerical time integration for air pollution models

NARCIS (Netherlands)

J.G. Verwer (Jan); W. Hundsdorfer (Willem); J.G. Blom (Joke)

1998-01-01

textabstractDue to the large number of chemical species and the three space dimensions, off-the-shelf stiff ODE integrators are not feasible for the numerical time integration of stiff systems of advection-diffusion-reaction equations [ fracpar{c{t + nabla cdot left( vu{u c right) = nabla cdot left(

A boundary integral method for two-dimensional (non)-Newtonian drops in slow viscous flow

NARCIS (Netherlands)

Toose, E.M.; Geurts, B.J.; Kuerten, J.G.M.

1995-01-01

A boundary integral method for the simulation of the time-dependent deformation of Newtonian or non-Newtonian drops suspended in a Newtonian fluid is developed. The boundary integral formulation for Stokes flow is used and the non-Newtonian stress is treated as a source term which yields an extra

Time-domain hybrid method for simulating large amplitude motions of ships advancing in waves

Directory of Open Access Journals (Sweden)

Shukui Liu

2011-03-01

Full Text Available Typical results obtained by a newly developed, nonlinear time domain hybrid method for simulating large amplitude motions of ships advancing with constant forward speed in waves are presented. The method is hybrid in the way of combining a time-domain transient Green function method and a Rankine source method. The present approach employs a simple double integration algorithm with respect to time to simulate the free-surface boundary condition. During the simulation, the diffraction and radiation forces are computed by pressure integration over the mean wetted surface, whereas the incident wave and hydrostatic restoring forces/moments are calculated on the instantaneously wetted surface of the hull. Typical numerical results of application of the method to the seakeeping performance of a standard containership, namely the ITTC S175, are herein presented. Comparisons have been made between the results from the present method, the frequency domain 3D panel method (NEWDRIFT of NTUA-SDL and available experimental data and good agreement has been observed for all studied cases between the results of the present method and comparable other data.

Achieving Integration in Mixed Methods Designs—Principles and Practices

Science.gov (United States)

Fetters, Michael D; Curry, Leslie A; Creswell, John W

2013-01-01

Mixed methods research offers powerful tools for investigating complex processes and systems in health and health care. This article describes integration principles and practices at three levels in mixed methods research and provides illustrative examples. Integration at the study design level occurs through three basic mixed method designs—exploratory sequential, explanatory sequential, and convergent—and through four advanced frameworks—multistage, intervention, case study, and participatory. Integration at the methods level occurs through four approaches. In connecting, one database links to the other through sampling. With building, one database informs the data collection approach of the other. When merging, the two databases are brought together for analysis. With embedding, data collection and analysis link at multiple points. Integration at the interpretation and reporting level occurs through narrative, data transformation, and joint display. The fit of integration describes the extent the qualitative and quantitative findings cohere. Understanding these principles and practices of integration can help health services researchers leverage the strengths of mixed methods. PMID:24279835

Atomization in graphite-furnace atomic absorption spectrometry. Peak-height method vs. integration method of measuring absorbance: carbon rod atomizer 63

International Nuclear Information System (INIS)

Sturgeon, R.E.; Chakrabarti, C.L.; Maines, I.S.; Bertels, P.C.

1975-01-01

Oscilloscopic traces of transient atomic absorption signals generated during continuous heating of a Carbon Rod Atomizer model 63 show features which are characteristic of the element being atomized. This research was undertaken to determine the significance and usefulness of the two analytically significant parameters, absorbance maximum and integrated absorbance. For measuring integrated absorbance, an electronic integrating control unit consisting of a timing circuit, a lock-in amplifier, and a digital voltmeter, which functions as a direct absorbance x second readout, has been designed, developed, and successfully tested. Oscilloscopic and recorder traces of the absorbance maximum and digital display of the integrated absorbance are simultaneously obtained. For the elements studied, Cd, Zn, Cu, Al, Sn, Mo, and V, the detection limits and the precision obtained are practically identical for both methods of measurements. The sensitivities by the integration method are about the same as, or less than, those obtained by the peak-height method, whereas the calibration curves by the former are generally linear over wider ranges of concentrations. (U.S.)

Integration of Small-Diameter Wood Harvesting in Early Thinnings using the Two pile Cutting Method

Energy Technology Data Exchange (ETDEWEB)

Kaerhae, Kalle (Metsaeteho Oy, P.O. Box 101, FI-00171 Helsinki (Finland))

2008-10-15

Metsaeteho Oy studied the integrated harvesting of industrial roundwood (pulpwood) and energy wood based on a two-pile cutting method, i.e. pulpwood and energy wood fractions are stacked into two separate piles when cutting a first-thinning stand. The productivity and cost levels of the integrated, two-pile cutting method were determined, and the harvesting costs of the two-pile method were compared with those of conventional separate wood harvesting methods. In the time study, when the size of removal was 50 dm3, the productivity in conventional whole-tree cutting was 6% higher than in integrated cutting. With a stem size of 100 dm3, the productivity in whole-tree cutting was 7% higher than in integrated cutting. The results indicated, however, that integrated harvesting based on the two-pile method enables harvesting costs to be decreased to below the current cost level of separate pulpwood harvesting in first thinning stands. The greatest cost-saving potential lies in small-sized first thinnings. The results showed that, when integrated wood harvesting based on the two-pile method is applied, the removals of both energy wood and pulpwood should be more than 15-20 m3/ha at the harvesting sites in order to achieve economically viable integrated procurement

Velocity time integral for right upper pulmonary vein in VLBW infants with patent ductus arteriosus

Directory of Open Access Journals (Sweden)

Gianluca Lista

Full Text Available OBJECTIVE: Early diagnosis of significant patent ductus arteriosus reduces the risk of clinical worsening in very low birth weight infants. Echocardiographic patent ductus arteriosus shunt flow pattern can be used to predict significant patent ductus arteriosus. Pulmonary venous flow, expressed as vein velocity time integral, is correlated to ductus arteriosus closure. The aim of this study is to investigate the relationship between significant reductions in vein velocity time integral and non-significant patent ductus arteriosus in the first week of life. METHODS: A multicenter, prospective, observational study was conducted to evaluate very low birth weight infants (<1500 g on respiratory support. Echocardiography was used to evaluate vein velocity time integral on days 1 and 4 of life. The relationship between vein velocity time integral and other parameters was studied. RESULTS: In total, 98 very low birth weight infants on respiratory support were studied. On day 1 of life, vein velocity time integral was similar in patients with open or closed ductus. The mean vein velocity time integral significantly reduced in the first four days of life. On the fourth day of life, there was less of a reduction in patients with patent ductus compared to those with closed patent ductus arteriosus and the difference was significant. CONCLUSIONS: A significant reduction in vein velocity time integral in the first days of life is associated with ductus closure. This parameter correlates well with other echocardiographic parameters and may aid in the diagnosis and management of patent ductus arteriosus.

Numerical method of singular problems on singular integrals

International Nuclear Information System (INIS)

Zhao Huaiguo; Mou Zongze

1992-02-01

As first part on the numerical research of singular problems, a numerical method is proposed for singular integrals. It is shown that the procedure is quite powerful for solving physics calculation with singularity such as the plasma dispersion function. Useful quadrature formulas for some class of the singular integrals are derived. In general, integrals with more complex singularities can be dealt by this method easily

Time step rescaling recovers continuous-time dynamical properties for discrete-time Langevin integration of nonequilibrium systems.

Science.gov (United States)

Sivak, David A; Chodera, John D; Crooks, Gavin E

2014-06-19

When simulating molecular systems using deterministic equations of motion (e.g., Newtonian dynamics), such equations are generally numerically integrated according to a well-developed set of algorithms that share commonly agreed-upon desirable properties. However, for stochastic equations of motion (e.g., Langevin dynamics), there is still broad disagreement over which integration algorithms are most appropriate. While multiple desiderata have been proposed throughout the literature, consensus on which criteria are important is absent, and no published integration scheme satisfies all desiderata simultaneously. Additional nontrivial complications stem from simulating systems driven out of equilibrium using existing stochastic integration schemes in conjunction with recently developed nonequilibrium fluctuation theorems. Here, we examine a family of discrete time integration schemes for Langevin dynamics, assessing how each member satisfies a variety of desiderata that have been enumerated in prior efforts to construct suitable Langevin integrators. We show that the incorporation of a novel time step rescaling in the deterministic updates of position and velocity can correct a number of dynamical defects in these integrators. Finally, we identify a particular splitting (related to the velocity Verlet discretization) that has essentially universally appropriate properties for the simulation of Langevin dynamics for molecular systems in equilibrium, nonequilibrium, and path sampling contexts.

Integral equation methods for vesicle electrohydrodynamics in three dimensions

Science.gov (United States)

Veerapaneni, Shravan

2016-12-01

In this paper, we develop a new boundary integral equation formulation that describes the coupled electro- and hydro-dynamics of a vesicle suspended in a viscous fluid and subjected to external flow and electric fields. The dynamics of the vesicle are characterized by a competition between the elastic, electric and viscous forces on its membrane. The classical Taylor-Melcher leaky-dielectric model is employed for the electric response of the vesicle and the Helfrich energy model combined with local inextensibility is employed for its elastic response. The coupled governing equations for the vesicle position and its transmembrane electric potential are solved using a numerical method that is spectrally accurate in space and first-order in time. The method uses a semi-implicit time-stepping scheme to overcome the numerical stiffness associated with the governing equations.

Calculations of the electromechanical transfer processes using implicit methods of numerical integration

Energy Technology Data Exchange (ETDEWEB)

Pogosyan, T A

1983-01-01

The article is dedicated to the solution of systems of differential equations which describe the transfer processes in an electric power system (EES) by implicit methods of numerical integration. The distinguishing feature of the implicit methods (Euler's reverse method and the trapeze method) is their absolute stability and, consequently, the relatively small accumulation of errors in each step of integration. Therefore, they are found to be very convenient for solving problems of electric power engineering, when the transfer processes are described by a rigid system of differential equations. The rigidity is associated with the range of values of the time constants considered. The advantage of the implicit methods over explicit are shown in a specific example (calculation of the dynamic stability of the simplest electric power system), along with the field of use of the implicit methods and the expedience of their use in power engineering problems.

[How timely are the methods taught in psychotherapy training and practice?].

Science.gov (United States)

Beutel, Manfred E; Michal, Matthias; Wiltink, Jörg; Subic-Wrana, Claudia

2015-01-01

Even though many psychotherapists consider themselves to be eclectic or integrative, training and reimbursement in the modern healthcare system are clearly oriented toward the model of distinct psychotherapy approaches. Prompted by the proposition to favor general, disorder-oriented psychotherapy, we investigate how timely distinctive methods are that are taught in training and practice. We reviewed the pertinent literature regarding general and specific factors, the effectiveness of integrative and eclectic treatments, orientation toward specific disorders, manualization and psychotherapeutic training. There is a lack of systematic studies on the efficacy of combining therapy methods from different approaches. The first empirical findings reveal that a superiority of combined versus single treatmentmethods has yet to be demonstrated. The development of transnosological manuals shows the limits of disorder-specific treatment.General factors such as therapeutic alliance or education about the model of disease and treatment rationale require specific definitions. Taking reference to a specific treatment approach provides important consistency of theory, training therapy and supervision, though this does not preclude an openness toward other therapy concepts. Current manualized examples show that methods and techniques can indeed be integrated from other approaches. Integrating different methods can also be seen as a developmental task for practitioners and researchers which may be mastered increasingly better with more experience.

Conservative fourth-order time integration of non-linear dynamic systems

DEFF Research Database (Denmark)

Krenk, Steen

2015-01-01

An energy conserving time integration algorithm with fourth-order accuracy is developed for dynamic systems with nonlinear stiffness. The discrete formulation is derived by integrating the differential state-space equations of motion over the integration time increment, and then evaluating...... the resulting time integrals of the inertia and stiffness terms via integration by parts. This process introduces the time derivatives of the state space variables, and these are then substituted from the original state-space differential equations. The resulting discrete form of the state-space equations...... is a direct fourth-order accurate representation of the original differential equations. This fourth-order form is energy conserving for systems with force potential in the form of a quartic polynomial in the displacement components. Energy conservation for a force potential of general form is obtained...

Investigation of the Adaptability of Transient Stability Assessment Methods to Real-Time Operation

OpenAIRE

Weckesser, Johannes Tilman Gabriel; Jóhannsson, Hjörtur; Sommer, Stefan; Østergaard, Jacob

2012-01-01

In this paper, an investigation of the adaptability of available transient stability assessment methods to real-time operation and their real-time performance is carried out. Two approaches based on Lyapunov’s method and the equal area criterion are analyzed. The results allow to determine the runtime of each method with respect to the number of inputs. Furthermore, it allows to identify, which method is preferable in case of changes in the power system such as the integration of distributed ...

Variational Symplectic Integrator for Long-Time Simulations of the Guiding-Center Motion of Charged Particles in General Magnetic Fields

International Nuclear Information System (INIS)

Qin Hong; Guan Xiaoyin

2008-01-01

A variational symplectic integrator for the guiding-center motion of charged particles in general magnetic fields is developed for long-time simulation studies of magnetized plasmas. Instead of discretizing the differential equations of the guiding-center motion, the action of the guiding-center motion is discretized and minimized to obtain the iteration rules for advancing the dynamics. The variational symplectic integrator conserves exactly a discrete Lagrangian symplectic structure, and has better numerical properties over long integration time, compared with standard integrators, such as the standard and variable time-step fourth order Runge-Kutta methods

Variational Symplectic Integrator for Long-Time Simulations of the Guiding-Center Motion of Charged Particles in General Magnetic Fields

International Nuclear Information System (INIS)

Qin, H.; Guan, X.

2008-01-01

A variational symplectic integrator for the guiding-center motion of charged particles in general magnetic fields is developed for long-time simulation studies of magnetized plasmas. Instead of discretizing the differential equations of the guiding-center motion, the action of the guiding-center motion is discretized and minimized to obtain the iteration rules for advancing the dynamics. The variational symplectic integrator conserves exactly a discrete Lagrangian symplectic structure, and has better numerical properties over long integration time, compared with standard integrators, such as the standard and variable time-step fourth order Runge-Kutta methods.

System Integration for Real-Time Mobile Manipulation

Directory of Open Access Journals (Sweden)

Reza Oftadeh

2014-03-01

Full Text Available Mobile manipulators are one of the most complicated types of mechatronics systems. The performance of these robots in performing complex manipulation tasks is highly correlated with the synchronization and integration of their low-level components. This paper discusses in detail the mechatronics design of a four wheel steered mobile manipulator. It presents the manipulator's mechanical structure and electrical interfaces, designs low-level software architecture based on embedded PC-based controls, and proposes a systematic solution based on code generation products of MATLAB and Simulink. The remote development environment described here is used to develop real-time controller software and modules for the mobile manipulator under a POSIX-compliant, real-time Linux operating system. Our approach enables developers to reliably design controller modules that meet the hard real-time constraints of the entire low-level system architecture. Moreover, it provides a systematic framework for the development and integration of hardware devices with various communication mediums and protocols, which facilitates the development and integration process of the software controller.

An Integrated Computational Materials Engineering Method for Woven Carbon Fiber Composites Preforming Process

Energy Technology Data Exchange (ETDEWEB)

Zhang, Weizhao; Ren, Huaqing; Wang, Zequn; Liu, Wing K.; Chen, Wei; Zeng, Danielle; Su, Xuming; Cao, Jian

2016-10-19

An integrated computational materials engineering method is proposed in this paper for analyzing the design and preforming process of woven carbon fiber composites. The goal is to reduce the cost and time needed for the mass production of structural composites. It integrates the simulation methods from the micro-scale to the macro-scale to capture the behavior of the composite material in the preforming process. In this way, the time consuming and high cost physical experiments and prototypes in the development of the manufacturing process can be circumvented. This method contains three parts: the micro-scale representative volume element (RVE) simulation to characterize the material; the metamodeling algorithm to generate the constitutive equations; and the macro-scale preforming simulation to predict the behavior of the composite material during forming. The results show the potential of this approach as a guidance to the design of composite materials and its manufacturing process.

Time for creative integration in medical sociology.

Science.gov (United States)

Levine, S

1995-01-01

The burgeoning of medical sociology has sometimes been accompanied by unfortunate parochialism and the presence of opposing intellectual camps that ignore and even impugn each other's work. We have lost opportunities to achieve creative discourse and integration of different perspectives, methods, and findings. At this stage we should consider how we can foster creative integration within our field.

Achieving integration in mixed methods designs-principles and practices.

Science.gov (United States)

Fetters, Michael D; Curry, Leslie A; Creswell, John W

2013-12-01

Mixed methods research offers powerful tools for investigating complex processes and systems in health and health care. This article describes integration principles and practices at three levels in mixed methods research and provides illustrative examples. Integration at the study design level occurs through three basic mixed method designs-exploratory sequential, explanatory sequential, and convergent-and through four advanced frameworks-multistage, intervention, case study, and participatory. Integration at the methods level occurs through four approaches. In connecting, one database links to the other through sampling. With building, one database informs the data collection approach of the other. When merging, the two databases are brought together for analysis. With embedding, data collection and analysis link at multiple points. Integration at the interpretation and reporting level occurs through narrative, data transformation, and joint display. The fit of integration describes the extent the qualitative and quantitative findings cohere. Understanding these principles and practices of integration can help health services researchers leverage the strengths of mixed methods. © Health Research and Educational Trust.

Local linearization methods for the numerical integration of ordinary differential equations: An overview

International Nuclear Information System (INIS)

Jimenez, J.C.

2009-06-01

Local Linearization (LL) methods conform a class of one-step explicit integrators for ODEs derived from the following primary and common strategy: the vector field of the differential equation is locally (piecewise) approximated through a first-order Taylor expansion at each time step, thus obtaining successive linear equations that are explicitly integrated. Hereafter, the LL approach may include some additional strategies to improve that basic affine approximation. Theoretical and practical results have shown that the LL integrators have a number of convenient properties. These include arbitrary order of convergence, A-stability, linearization preserving, regularity under quite general conditions, preservation of the dynamics of the exact solution around hyperbolic equilibrium points and periodic orbits, integration of stiff and high-dimensional equations, low computational cost, and others. In this paper, a review of the LL methods and their properties is presented. (author)

Space-time transformations in radial path integrals

International Nuclear Information System (INIS)

Steiner, F.

1984-09-01

Nonlinear space-time transformations in the radial path integral are discussed. A transformation formula is derived, which relates the original path integral to the Green's function of a new quantum system with an effective potential containing an observable quantum correction proportional(h/2π) 2 . As an example the formula is applied to spherical Brownian motion. (orig.)

Boundary-integral equation formulation for time-dependent inelastic deformation in metals

Energy Technology Data Exchange (ETDEWEB)

Kumar, V; Mukherjee, S

1977-01-01

The mathematical structure of various constitutive relations proposed in recent years for representing time-dependent inelastic deformation behavior of metals at elevated temperatues has certain features which permit a simple formulation of the three-dimensional inelasticity problem in terms of real time rates. A direct formulation of the boundary-integral equation method in terms of rates is discussed for the analysis of time-dependent inelastic deformation of arbitrarily shaped three-dimensional metallic bodies subjected to arbitrary mechanical and thermal loading histories and obeying constitutive relations of the kind mentioned above. The formulation is based on the assumption of infinitesimal deformations. Several illustrative examples involving creep of thick-walled spheres, long thick-walled cylinders, and rotating discs are discussed. The implementation of the method appears to be far easier than analogous BIE formulations that have been suggested for elastoplastic problems.

A simple method to adapt time sampling of the analog signal

International Nuclear Information System (INIS)

Kalinin, Yu.G.; Martyanov, I.S.; Sadykov, Kh.; Zastrozhnova, N.N.

2004-01-01

In this paper we briefly describe the time sampling method, which is adapted to the speed of the signal change. Principally, this method is based on a simple idea--the combination of discrete integration with differentiation of the analog signal. This method can be used in nuclear electronics research into the characteristics of detectors and the shape of the pulse signal, pulse and transitive characteristics of inertial systems of processing of signals, etc

Evaluation of integration methods for hybrid simulation of complex structural systems through collapse

Science.gov (United States)

Del Carpio R., Maikol; Hashemi, M. Javad; Mosqueda, Gilberto

2017-10-01

This study examines the performance of integration methods for hybrid simulation of large and complex structural systems in the context of structural collapse due to seismic excitations. The target application is not necessarily for real-time testing, but rather for models that involve large-scale physical sub-structures and highly nonlinear numerical models. Four case studies are presented and discussed. In the first case study, the accuracy of integration schemes including two widely used methods, namely, modified version of the implicit Newmark with fixed-number of iteration (iterative) and the operator-splitting (non-iterative) is examined through pure numerical simulations. The second case study presents the results of 10 hybrid simulations repeated with the two aforementioned integration methods considering various time steps and fixed-number of iterations for the iterative integration method. The physical sub-structure in these tests consists of a single-degree-of-freedom (SDOF) cantilever column with replaceable steel coupons that provides repeatable highlynonlinear behavior including fracture-type strength and stiffness degradations. In case study three, the implicit Newmark with fixed-number of iterations is applied for hybrid simulations of a 1:2 scale steel moment frame that includes a relatively complex nonlinear numerical substructure. Lastly, a more complex numerical substructure is considered by constructing a nonlinear computational model of a moment frame coupled to a hybrid model of a 1:2 scale steel gravity frame. The last two case studies are conducted on the same porotype structure and the selection of time steps and fixed number of iterations are closely examined in pre-test simulations. The generated unbalance forces is used as an index to track the equilibrium error and predict the accuracy and stability of the simulations.

Integral Methods in Science and Engineering

CERN Document Server

Constanda, Christian

2011-01-01

An enormous array of problems encountered by scientists and engineers are based on the design of mathematical models using many different types of ordinary differential, partial differential, integral, and integro-differential equations. Accordingly, the solutions of these equations are of great interest to practitioners and to science in general. Presenting a wealth of cutting-edge research by a diverse group of experts in the field, Integral Methods in Science and Engineering: Computational and Analytic Aspects gives a vivid picture of both the development of theoretical integral techniques

Multisensory integration: the case of a time window of gesture-speech integration.

Science.gov (United States)

Obermeier, Christian; Gunter, Thomas C

2015-02-01

This experiment investigates the integration of gesture and speech from a multisensory perspective. In a disambiguation paradigm, participants were presented with short videos of an actress uttering sentences like "She was impressed by the BALL, because the GAME/DANCE...." The ambiguous noun (BALL) was accompanied by an iconic gesture fragment containing information to disambiguate the noun toward its dominant or subordinate meaning. We used four different temporal alignments between noun and gesture fragment: the identification point (IP) of the noun was either prior to (+120 msec), synchronous with (0 msec), or lagging behind the end of the gesture fragment (-200 and -600 msec). ERPs triggered to the IP of the noun showed significant differences for the integration of dominant and subordinate gesture fragments in the -200, 0, and +120 msec conditions. The outcome of this integration was revealed at the target words. These data suggest a time window for direct semantic gesture-speech integration ranging from at least -200 up to +120 msec. Although the -600 msec condition did not show any signs of direct integration at the homonym, significant disambiguation was found at the target word. An explorative analysis suggested that gesture information was directly integrated at the verb, indicating that there are multiple positions in a sentence where direct gesture-speech integration takes place. Ultimately, this would implicate that in natural communication, where a gesture lasts for some time, several aspects of that gesture will have their specific and possibly distinct impact on different positions in an utterance.

Marching on-in-time solution of the time domain magnetic field integral equation using a predictor-corrector scheme

KAUST Repository

Ulku, Huseyin Arda; Bagci, Hakan; Michielssen, Eric

2013-01-01

An explicit marching on-in-time (MOT) scheme for solving the time-domain magnetic field integral equation (TD-MFIE) is presented. The proposed MOT-TD-MFIE solver uses Rao-Wilton-Glisson basis functions for spatial discretization and a PE(CE)m-type linear multistep method for time marching. Unlike previous explicit MOT-TD-MFIE solvers, the time step size can be chosen as large as that of the implicit MOT-TD-MFIE solvers without adversely affecting accuracy or stability. An algebraic stability analysis demonstrates the stability of the proposed explicit solver; its accuracy and efficiency are established via numerical examples. © 1963-2012 IEEE.

Marching on-in-time solution of the time domain magnetic field integral equation using a predictor-corrector scheme

KAUST Repository

Ulku, Huseyin Arda

2013-08-01

An explicit marching on-in-time (MOT) scheme for solving the time-domain magnetic field integral equation (TD-MFIE) is presented. The proposed MOT-TD-MFIE solver uses Rao-Wilton-Glisson basis functions for spatial discretization and a PE(CE)m-type linear multistep method for time marching. Unlike previous explicit MOT-TD-MFIE solvers, the time step size can be chosen as large as that of the implicit MOT-TD-MFIE solvers without adversely affecting accuracy or stability. An algebraic stability analysis demonstrates the stability of the proposed explicit solver; its accuracy and efficiency are established via numerical examples. © 1963-2012 IEEE.

Recent Advances in the Method of Forces: Integrated Force Method of Structural Analysis

Science.gov (United States)

Patnaik, Surya N.; Coroneos, Rula M.; Hopkins, Dale A.

1998-01-01

Stress that can be induced in an elastic continuum can be determined directly through the simultaneous application of the equilibrium equations and the compatibility conditions. In the literature, this direct stress formulation is referred to as the integrated force method. This method, which uses forces as the primary unknowns, complements the popular equilibrium-based stiffness method, which considers displacements as the unknowns. The integrated force method produces accurate stress, displacement, and frequency results even for modest finite element models. This version of the force method should be developed as an alternative to the stiffness method because the latter method, which has been researched for the past several decades, may have entered its developmental plateau. Stress plays a primary role in the development of aerospace and other products, and its analysis is difficult. Therefore, it is advisable to use both methods to calculate stress and eliminate errors through comparison. This paper examines the role of the integrated force method in analysis, animation and design.

ASPECTS OF INTEGRATION MANAGEMENT METHODS

Directory of Open Access Journals (Sweden)

Artemy Varshapetian

2015-10-01

Full Text Available For manufacturing companies to succeed in today's unstable economic environment, it is necessary to restructure the main components of its activities: designing innovative product, production using modern reconfigurable manufacturing systems, a business model that takes into account the global strategy and management methods using modern management models and tools. The first three components are discussed in numerous publications, for example, (Koren, 2010 and is therefore not considered in the article. A large number of publications devoted to the methods and tools of production management, for example (Halevi, 2007. On the basis of what was said in the article discusses the possibility of the integration of only three methods have received in recent years, the most widely used, namely: Six Sigma method - SS (George et al., 2005 and supplements its-Design for six sigm? - DFSS (Taguchi, 2003; Lean production transformed with the development to the "Lean management" and further to the "Lean thinking" - Lean (Hirano et al., 2006; Theory of Constraints, developed E.Goldratt - TOC (Dettmer, 2001. The article investigates some aspects of this integration: applications in diverse fields, positive features, changes in management structure, etc.

Retarded potentials and time domain boundary integral equations a road map

CERN Document Server

Sayas, Francisco-Javier

2016-01-01

This book offers a thorough and self-contained exposition of the mathematics of time-domain boundary integral equations associated to the wave equation, including applications to scattering of acoustic and elastic waves. The book offers two different approaches for the analysis of these integral equations, including a systematic treatment of their numerical discretization using Galerkin (Boundary Element) methods in the space variables and Convolution Quadrature in the time variable. The first approach follows classical work started in the late eighties, based on Laplace transforms estimates. This approach has been refined and made more accessible by tailoring the necessary mathematical tools, avoiding an excess of generality. A second approach contains a novel point of view that the author and some of his collaborators have been developing in recent years, using the semigroup theory of evolution equations to obtain improved results. The extension to electromagnetic waves is explained in one of the appendices...

A theoretical study for the real-time assessment of external gamma exposure using equivalent-volume numerical integration

International Nuclear Information System (INIS)

Han, Moon Hee

1995-02-01

An approximate method for estimating gamma external dose due to an arbitrary distribution of radioactive material has been developed. For the assessment of external gamma dose, the space over which radioactive material is distributed has been assumed to be composed of hexagonal cells. The evaluation of three-dimensional integration over the space is an extremely time-consuming task. Hence, a different approach has been used for the study, i.e., a equivalent-volume spherical approach in which a regular hexahedron is modeled as a equivalent-volume sphere to simplify the integration. For the justification of the current approach, two case studies have been performed: a comparison with a point source approximation and a comparison of external dose rate with the Monte Carlo integration. These comparisons show that the current approach gives reasonable results in a physical sense. Computing times of the developed and Monte Carlo integration method on VAX system have been compared as a function of the number of hexagonal cells. This comparison shows that CPU times for both methods are comparable in the region of small number of cells, but in the region of large number, Monte Carlo integration needs much more computing times. The proposed method is shown to have an accuracy equivalent to Monte Carlo method with an advantage of much shorter calculation time. Then, the method developed here evaluates early off-site consequences of a nuclear accident. An accident consequence assessment model has been integrated using Gaussian puff model which is used to obtain the distribution of radioactive material in the air and on the ground. For this work, the real meteorological data measured at Kori site for 10 years (1976 - 1985) have been statistically analyzed for obtaining site-specific conditions. The short-term external gamma exposures have been assessed for several site-specific meteorological conditions. The results show that the extent and the pattern of short-term external

Symplectic integrators for large scale molecular dynamics simulations: A comparison of several explicit methods

International Nuclear Information System (INIS)

Gray, S.K.; Noid, D.W.; Sumpter, B.G.

1994-01-01

We test the suitability of a variety of explicit symplectic integrators for molecular dynamics calculations on Hamiltonian systems. These integrators are extremely simple algorithms with low memory requirements, and appear to be well suited for large scale simulations. We first apply all the methods to a simple test case using the ideas of Berendsen and van Gunsteren. We then use the integrators to generate long time trajectories of a 1000 unit polyethylene chain. Calculations are also performed with two popular but nonsymplectic integrators. The most efficient integrators of the set investigated are deduced. We also discuss certain variations on the basic symplectic integration technique

Integrative real-time geographic visualization of energy resources

International Nuclear Information System (INIS)

Sorokine, A.; Shankar, M.; Stovall, J.; Bhaduri, B.; King, T.; Fernandez, S.; Datar, N.; Omitaomu, O.

2009-01-01

'Full text:' Several models forecast that climatic changes will increase the frequency of disastrous events like droughts, hurricanes, and snow storms. Responding to these events and also to power outages caused by system errors such as the 2003 North American blackout require an interconnect-wide real-time monitoring system for various energy resources. Such a system should be capable of providing situational awareness to its users in the government and energy utilities by dynamically visualizing the status of the elements of the energy grid infrastructure and supply chain in geographic contexts. We demonstrate an approach that relies on Google Earth and similar standard-based platforms as client-side geographic viewers with a data-dependent server component. The users of the system can view status information in spatial and temporal contexts. These data can be integrated with a wide range of geographic sources including all standard Google Earth layers and a large number of energy and environmental data feeds. In addition, we show a real-time spatio-temporal data sharing capability across the users of the system, novel methods for visualizing dynamic network data, and a fine-grain access to very large multi-resolution geographic datasets for faster delivery of the data. The system can be extended to integrate contingency analysis results and other grid models to assess recovery and repair scenarios in the case of major disruption. (author)

3-D electromagnetic modeling for very early time sounding of shallow targets using integral equations

International Nuclear Information System (INIS)

Xiong, Z.; Tripp, A.C.

1994-01-01

This paper presents an integral equation algorithm for 3D EM modeling at high frequencies for applications in engineering an environmental studies. The integral equation method remains the same for low and high frequencies, but the dominant roles of the displacements currents complicate both numerical treatments and interpretations. With singularity extraction technique they successively extended the application of the Hankel filtering technique to the computation of Hankel integrals occurring in high frequency EM modeling. Time domain results are calculated from frequency domain results via Fourier transforms. While frequency domain data are not obvious for interpretations, time domain data show wave-like pictures that resemble seismograms. Both 1D and 3D numerical results show clearly the layer interfaces

An arbitrary-order staggered time integrator for the linear acoustic wave equation

Science.gov (United States)

Lee, Jaejoon; Park, Hyunseo; Park, Yoonseo; Shin, Changsoo

2018-02-01

We suggest a staggered time integrator whose order of accuracy can arbitrarily be extended to solve the linear acoustic wave equation. A strategy to select the appropriate order of accuracy is also proposed based on the error analysis that quantitatively predicts the truncation error of the numerical solution. This strategy not only reduces the computational cost several times, but also allows us to flexibly set the modelling parameters such as the time step length, grid interval and P-wave speed. It is demonstrated that the proposed method can almost eliminate temporal dispersive errors during long term simulations regardless of the heterogeneity of the media and time step lengths. The method can also be successfully applied to the source problem with an absorbing boundary condition, which is frequently encountered in the practical usage for the imaging algorithms or the inverse problems.

A model of interval timing by neural integration.

Science.gov (United States)

Simen, Patrick; Balci, Fuat; de Souza, Laura; Cohen, Jonathan D; Holmes, Philip

2011-06-22

We show that simple assumptions about neural processing lead to a model of interval timing as a temporal integration process, in which a noisy firing-rate representation of time rises linearly on average toward a response threshold over the course of an interval. Our assumptions include: that neural spike trains are approximately independent Poisson processes, that correlations among them can be largely cancelled by balancing excitation and inhibition, that neural populations can act as integrators, and that the objective of timed behavior is maximal accuracy and minimal variance. The model accounts for a variety of physiological and behavioral findings in rodents, monkeys, and humans, including ramping firing rates between the onset of reward-predicting cues and the receipt of delayed rewards, and universally scale-invariant response time distributions in interval timing tasks. It furthermore makes specific, well-supported predictions about the skewness of these distributions, a feature of timing data that is usually ignored. The model also incorporates a rapid (potentially one-shot) duration-learning procedure. Human behavioral data support the learning rule's predictions regarding learning speed in sequences of timed responses. These results suggest that simple, integration-based models should play as prominent a role in interval timing theory as they do in theories of perceptual decision making, and that a common neural mechanism may underlie both types of behavior.

INTEGRAL EDUCATION, TIME AND SPACE: PROBLEMATIZING CONCEPTS

Directory of Open Access Journals (Sweden)

Ana Elisa Spaolonzi Queiroz Assis

2018-03-01

Full Text Available Integral Education, despite being the subject of public policy agenda for some decades, still carries disparities related to its concept. In this sense, this article aims to problematize not only the concepts of integral education but also the categories time and space contained in the magazines Em Aberto. They were organized and published by the National Institute of Educational Studies Anísio Teixeira (INEP, numbers 80 (2009 and 88 (2012, respectively entitled "Educação Integral e tempo integral" and " Políticas de educação integral em jornada ampliada". The methodology is based on Bardin’s content analysis, respecting the steps of pre-analysis (research corpus formed by the texts in the journals; material exploration (reading the texts encoding data choosing the registration units for categorization; and processing and interpretation of results, based on Saviani’s Historical-Critical Pedagogy. The work reveals convergent and divergent conceptual multiplicity, provoking a discussion about a critical conception of integral education. Keywords: Integral Education. Historical-Critical Pedagogy. Content Analysis.

Application of the heat-balance and refined integral methods to the Korteweg-de Vries equation

Directory of Open Access Journals (Sweden)

Myers Timothy G.

2009-01-01

Full Text Available In this paper we consider approximate travelling wave solutions to the Korteweg-de Vries equation. The heat-balance integral method is first applied to the problem, using two different quartic approximating functions, and then the refined integral method is investigated. We examine two types of solution, chosen by matching the wave speed to that of the exact solution and by imposing the same area. The first set of solutions is generally better with an error that is fixed in time. The second set of solutions has an error that grows with time. This is shown to be due to slight discrepancies in the wave speed.

Method for Determining the Time Parameter

Directory of Open Access Journals (Sweden)

K. P. Baslyk

2014-01-01

values of multipliers and parameter of penalty. The solving problem of the optimization without constraints is one step of the optimization problem with constraints.The proposed method was realized as a computational Pascal-language program. There is the multiple call of Runge-Kutta method procedure for integration of motion equations. To reduce operation-use time, normalization of motion equations is used, and the 4-th order RungeKutta method time step accuracy control is also applied.The program test results were compared with the solution, which was obtained using the existing software for the one stage missile design.The general results obtained in this paper are following: numerical determination method of the attack angle time parameter and maximum allowed value of the attack angle amplitude as functions of the projected and ballistic parameters; software implementation of this method. Deviation of angle trajectory value at the end of active trajectory leg as a function of the error of the time parameter was obtained as well. This paper can be used in the courses of learning such as “Introduction to rocketry” and “Launch vehicle design”.

Energy drift in reversible time integration

International Nuclear Information System (INIS)

McLachlan, R I; Perlmutter, M

2004-01-01

Energy drift is commonly observed in reversible integrations of systems of molecular dynamics. We show that this drift can be modelled as a diffusion and that the typical energy error after time T is O(√T). (letter to the editor)

Structural reliability calculation method based on the dual neural network and direct integration method.

Science.gov (United States)

Li, Haibin; He, Yun; Nie, Xiaobo

2018-01-01

Structural reliability analysis under uncertainty is paid wide attention by engineers and scholars due to reflecting the structural characteristics and the bearing actual situation. The direct integration method, started from the definition of reliability theory, is easy to be understood, but there are still mathematics difficulties in the calculation of multiple integrals. Therefore, a dual neural network method is proposed for calculating multiple integrals in this paper. Dual neural network consists of two neural networks. The neural network A is used to learn the integrand function, and the neural network B is used to simulate the original function. According to the derivative relationships between the network output and the network input, the neural network B is derived from the neural network A. On this basis, the performance function of normalization is employed in the proposed method to overcome the difficulty of multiple integrations and to improve the accuracy for reliability calculations. The comparisons between the proposed method and Monte Carlo simulation method, Hasofer-Lind method, the mean value first-order second moment method have demonstrated that the proposed method is an efficient and accurate reliability method for structural reliability problems.

Time-Lapse Measurement of Wellbore Integrity

Science.gov (United States)

Duguid, A.

2017-12-01

Well integrity is becoming more important as wells are used longer or repurposed. For CO2, shale gas, and other projects it has become apparent that wells represent the most likely unintended migration pathway for fluids out of the reservoir. Comprehensive logging programs have been employed to determine the condition of legacy wells in North America. These studies provide examples of assessment technologies. Logging programs have included pulsed neutron logging, ultrasonic well mapping, and cement bond logging. While these studies provide examples of what can be measured, they have only conducted a single round of logging and cannot show if the well has changed over time. Recent experience with time-lapse logging of three monitoring wells at a US Department of Energy sponsored CO2 project has shown the full value of similar tools. Time-lapse logging has shown that well integrity changes over time can be identified. It has also shown that the inclusion of and location of monitoring technologies in the well and the choice of construction materials must be carefully considered. Two of the wells were approximately eight years old at the time of study; they were constructed with steel and fiberglass casing sections and had lines on the outside of the casing running to the surface. The third well was 68 years old when it was studied and was originally constructed as a production well. Repeat logs were collected six or eight years after initial logging. Time-lapse logging showed the evolution of the wells. The results identified locations where cement degraded over time and locations that showed little change. The ultrasonic well maps show clearly that the lines used to connect the monitoring technology to the surface are visible and have a local effect on cement isolation. Testing and sampling was conducted along with logging. It provided insight into changes identified in the time-lapse log results. Point permeability testing was used to provide an in-situ point

A new integrated dual time-point amyloid PET/MRI data analysis method

International Nuclear Information System (INIS)

Cecchin, Diego; Zucchetta, Pietro; Turco, Paolo; Bui, Franco; Barthel, Henryk; Tiepolt, Solveig; Sabri, Osama; Poggiali, Davide; Cagnin, Annachiara; Gallo, Paolo; Frigo, Anna Chiara

2017-01-01

In the initial evaluation of patients with suspected dementia and Alzheimer's disease, there is no consensus on how to perform semiquantification of amyloid in such a way that it: (1) facilitates visual qualitative interpretation, (2) takes the kinetic behaviour of the tracer into consideration particularly with regard to at least partially correcting for blood flow dependence, (3) analyses the amyloid load based on accurate parcellation of cortical and subcortical areas, (4) includes partial volume effect correction (PVEC), (5) includes MRI-derived topographical indexes, (6) enables application to PET/MRI images and PET/CT images with separately acquired MR images, and (7) allows automation. A method with all of these characteristics was retrospectively tested in 86 subjects who underwent amyloid ( 18 F-florbetaben) PET/MRI in a clinical setting (using images acquired 90-110 min after injection, 53 were classified visually as amyloid-negative and 33 as amyloid-positive). Early images after tracer administration were acquired between 0 and 10 min after injection, and later images were acquired between 90 and 110 min after injection. PVEC of the PET data was carried out using the geometric transfer matrix method. Parametric images and some regional output parameters, including two innovative ''dual time-point'' indexes, were obtained. Subjects classified visually as amyloid-positive showed a sparse tracer uptake in the primary sensory, motor and visual areas in accordance with the isocortical stage of the topographic distribution of the amyloid plaque (Braak stages V/VI). In patients classified visually as amyloid-negative, the method revealed detectable levels of tracer uptake in the basal portions of the frontal and temporal lobes, areas that are known to be sites of early deposition of amyloid plaques that probably represented early accumulation (Braak stage A) that is typical of normal ageing. There was a strong correlation between age

A new integrated dual time-point amyloid PET/MRI data analysis method

Energy Technology Data Exchange (ETDEWEB)

Cecchin, Diego; Zucchetta, Pietro; Turco, Paolo; Bui, Franco [University Hospital of Padua, Nuclear Medicine Unit, Department of Medicine - DIMED, Padua (Italy); Barthel, Henryk; Tiepolt, Solveig; Sabri, Osama [Leipzig University, Department of Nuclear Medicine, Leipzig (Germany); Poggiali, Davide; Cagnin, Annachiara; Gallo, Paolo [University Hospital of Padua, Neurology, Department of Neurosciences (DNS), Padua (Italy); Frigo, Anna Chiara [University Hospital of Padua, Biostatistics, Epidemiology and Public Health Unit, Department of Cardiac, Thoracic and Vascular Sciences, Padua (Italy)

2017-11-15

In the initial evaluation of patients with suspected dementia and Alzheimer's disease, there is no consensus on how to perform semiquantification of amyloid in such a way that it: (1) facilitates visual qualitative interpretation, (2) takes the kinetic behaviour of the tracer into consideration particularly with regard to at least partially correcting for blood flow dependence, (3) analyses the amyloid load based on accurate parcellation of cortical and subcortical areas, (4) includes partial volume effect correction (PVEC), (5) includes MRI-derived topographical indexes, (6) enables application to PET/MRI images and PET/CT images with separately acquired MR images, and (7) allows automation. A method with all of these characteristics was retrospectively tested in 86 subjects who underwent amyloid ({sup 18}F-florbetaben) PET/MRI in a clinical setting (using images acquired 90-110 min after injection, 53 were classified visually as amyloid-negative and 33 as amyloid-positive). Early images after tracer administration were acquired between 0 and 10 min after injection, and later images were acquired between 90 and 110 min after injection. PVEC of the PET data was carried out using the geometric transfer matrix method. Parametric images and some regional output parameters, including two innovative ''dual time-point'' indexes, were obtained. Subjects classified visually as amyloid-positive showed a sparse tracer uptake in the primary sensory, motor and visual areas in accordance with the isocortical stage of the topographic distribution of the amyloid plaque (Braak stages V/VI). In patients classified visually as amyloid-negative, the method revealed detectable levels of tracer uptake in the basal portions of the frontal and temporal lobes, areas that are known to be sites of early deposition of amyloid plaques that probably represented early accumulation (Braak stage A) that is typical of normal ageing. There was a strong correlation between

Integral reinforcement learning for continuous-time input-affine nonlinear systems with simultaneous invariant explorations.

Science.gov (United States)

Lee, Jae Young; Park, Jin Bae; Choi, Yoon Ho

2015-05-01

This paper focuses on a class of reinforcement learning (RL) algorithms, named integral RL (I-RL), that solve continuous-time (CT) nonlinear optimal control problems with input-affine system dynamics. First, we extend the concepts of exploration, integral temporal difference, and invariant admissibility to the target CT nonlinear system that is governed by a control policy plus a probing signal called an exploration. Then, we show input-to-state stability (ISS) and invariant admissibility of the closed-loop systems with the policies generated by integral policy iteration (I-PI) or invariantly admissible PI (IA-PI) method. Based on these, three online I-RL algorithms named explorized I-PI and integral Q -learning I, II are proposed, all of which generate the same convergent sequences as I-PI and IA-PI under the required excitation condition on the exploration. All the proposed methods are partially or completely model free, and can simultaneously explore the state space in a stable manner during the online learning processes. ISS, invariant admissibility, and convergence properties of the proposed methods are also investigated, and related with these, we show the design principles of the exploration for safe learning. Neural-network-based implementation methods for the proposed schemes are also presented in this paper. Finally, several numerical simulations are carried out to verify the effectiveness of the proposed methods.

Simulating variable-density flows with time-consistent integration of Navier-Stokes equations

Science.gov (United States)

Lu, Xiaoyi; Pantano, Carlos

2017-11-01

In this talk, we present several features of a high-order semi-implicit variable-density low-Mach Navier-Stokes solver. A new formulation to solve pressure Poisson-like equation of variable-density flows is highlighted. With this formulation of the numerical method, we are able to solve all variables with a uniform order of accuracy in time (consistent with the time integrator being used). The solver is primarily designed to perform direct numerical simulations for turbulent premixed flames. Therefore, we also address other important elements, such as energy-stable boundary conditions, synthetic turbulence generation, and flame anchoring method. Numerical examples include classical non-reacting constant/variable-density flows, as well as turbulent premixed flames.

Numerical integration methods and layout improvements in the context of dynamic RNA visualization.

Science.gov (United States)

Shabash, Boris; Wiese, Kay C

2017-05-30

RNA visualization software tools have traditionally presented a static visualization of RNA molecules with limited ability for users to interact with the resulting image once it is complete. Only a few tools allowed for dynamic structures. One such tool is jViz.RNA. Currently, jViz.RNA employs a unique method for the creation of the RNA molecule layout by mapping the RNA nucleotides into vertexes in a graph, which we call the detailed graph, and then utilizes a Newtonian mechanics inspired system of forces to calculate a layout for the RNA molecule. The work presented here focuses on improvements to jViz.RNA that allow the drawing of RNA secondary structures according to common drawing conventions, as well as dramatic run-time performance improvements. This is done first by presenting an alternative method for mapping the RNA molecule into a graph, which we call the compressed graph, and then employing advanced numerical integration methods for the compressed graph representation. Comparing the compressed graph and detailed graph implementations, we find that the compressed graph produces results more consistent with RNA drawing conventions. However, we also find that employing the compressed graph method requires a more sophisticated initial layout to produce visualizations that would require minimal user interference. Comparing the two numerical integration methods demonstrates the higher stability of the Backward Euler method, and its resulting ability to handle much larger time steps, a high priority feature for any software which entails user interaction. The work in this manuscript presents the preferred use of compressed graphs to detailed ones, as well as the advantages of employing the Backward Euler method over the Forward Euler method. These improvements produce more stable as well as visually aesthetic representations of the RNA secondary structures. The results presented demonstrate that both the compressed graph representation, as well as the Backward

An integrated approach for facilities planning by ELECTRE method

Science.gov (United States)

Elbishari, E. M. Y.; Hazza, M. H. F. Al; Adesta, E. Y. T.; Rahman, Nur Salihah Binti Abdul

2018-01-01

Facility planning is concerned with the design, layout, and accommodation of people, machines and activities of a system. Most of the researchers try to investigate the production area layout and the related facilities. However, few of them try to investigate the relationship between the production space and its relationship with service departments. The aim of this research to is to integrate different approaches in order to evaluate, analyse and select the best facilities planning method that able to explain the relationship between the production area and other supporting departments and its effect on human efforts. To achieve the objective of this research two different approaches have been integrated: Apple’s layout procedure as one of the effective tools in planning factories, ELECTRE method as one of the Multi Criteria Decision Making methods (MCDM) to minimize the risk of getting poor facilities planning. Dalia industries have been selected as a case study to implement our integration the factory have been divided two main different area: the whole facility (layout A), and the manufacturing area (layout B). This article will be concerned with the manufacturing area layout (Layout B). After analysing the data gathered, the manufacturing area was divided into 10 activities. There are five factors that the alternative were compared upon which are: Inter department satisfactory level, total distance travelled for workers, total distance travelled for the product, total time travelled for the workers, and total time travelled for the product. Three different layout alternatives have been developed in addition to the original layouts. Apple’s layout procedure was used to study and evaluate the different alternatives layouts, the study and evaluation of the layouts was done by calculating scores for each of the factors. After obtaining the scores from evaluating the layouts, ELECTRE method was used to compare the proposed alternatives with each other and with

Nonperturbative time-convolutionless quantum master equation from the path integral approach

International Nuclear Information System (INIS)

Nan Guangjun; Shi Qiang; Shuai Zhigang

2009-01-01

The time-convolutionless quantum master equation is widely used to simulate reduced dynamics of a quantum system coupled to a bath. However, except for several special cases, applications of this equation are based on perturbative calculation of the dissipative tensor, and are limited to the weak system-bath coupling regime. In this paper, we derive an exact time-convolutionless quantum master equation from the path integral approach, which provides a new way to calculate the dissipative tensor nonperturbatively. Application of the new method is demonstrated in the case of an asymmetrical two-level system linearly coupled to a harmonic bath.

Optimization of the integration time of pulse shape analysis for dual-layer GSO detector with different amount of Ce

International Nuclear Information System (INIS)

Yamamoto, Seiichi

2008-01-01

For a multi-layer depth-of-interaction (DOI) detector using different decay times, pulse shape analysis based on two different integration times is often used to distinguish scintillators in DOI direction. This method measures a partial integration and a full integration, and calculates the ratio of these two to obtain the pulse shape distribution. The full integration time is usually set to integrate full width of the scintillation pulse. However, the optimum partial integration time is not obvious for obtaining the best separation of the pulse shape distribution. To make it clear, a theoretical analysis and experiments were conducted for pulse shape analysis by changing the partial integration time using a scintillation detector of GSOs with different amount of Ce. A scintillation detector with 1-in. round photomultiplier tube (PMT) optically coupled GSO of 1.5 mol% (decay time: 35 ns) and that of 0.5 mol% (decay time: 60 ns) was used for the experiments. The signal from PMT was digitally integrated with partial (50-150 ns) and full (160 ns) integration times and ratio of these two was calculated to obtain the pulse shape distribution. In the theoretical analysis, partial integration time of 50 ns showed largest distance between two peaks of the pulse shape distribution. In the experiments, it showed maximum at 70-80 ns of partial integration time. The peak to valley ratio showed the maximum at 120-130 ns. Because the separation of two peaks is determined by the peak to valley ratio, we conclude the optimum partial integration time for these combinations of GSOs is around 120-130 ns, relatively longer than the expected value

Time-dependent integral equations of neutron transport for calculating the kinetics of nuclear reactors by the Monte Carlo method

Energy Technology Data Exchange (ETDEWEB)

Davidenko, V. D., E-mail: Davidenko-VD@nrcki.ru; Zinchenko, A. S., E-mail: zin-sn@mail.ru; Harchenko, I. K. [National Research Centre Kurchatov Institute (Russian Federation)

2016-12-15

Integral equations for the shape functions in the adiabatic, quasi-static, and improved quasi-static approximations are presented. The approach to solving these equations by the Monte Carlo method is described.

Life cycle integrated thermoeconomic assessment method for energy conversion systems

International Nuclear Information System (INIS)

Kanbur, Baris Burak; Xiang, Liming; Dubey, Swapnil; Choo, Fook Hoong; Duan, Fei

2017-01-01

Highlights: • A new LCA integrated thermoeconomic approach is presented. • The new unit fuel cost is found 4.8 times higher than the classic method. • The new defined parameter increased the sustainability index by 67.1%. • The case studies are performed for countries with different CO 2 prices. - Abstract: Life cycle assessment (LCA) based thermoeconomic modelling has been applied for the evaluation of energy conversion systems since it provided more comprehensive and applicable assessment criteria. This study proposes an improved thermoeconomic method, named as life cycle integrated thermoeconomic assessment (LCiTA), which combines the LCA based enviroeconomic parameters in the production steps of the system components and fuel with the conventional thermoeconomic method for the energy conversion systems. A micro-cogeneration system is investigated and analyzed with the LCiTA method, the comparative studies show that the unit cost of fuel by using the LCiTA method is 3.8 times higher than the conventional thermoeconomic model. It is also realized that the enviroeconomic parameters during the operation of the system components do not have significant impacts on the system streams since the exergetic parameters are dominant in the thermoeconomic calculations. Moreover, the improved sustainability index is found roundly 67.2% higher than the previously defined sustainability index, suggesting that the enviroeconomic and thermoeconomic parameters decrease the impact of the exergy destruction in the sustainability index definition. To find the feasible operation conditions for the micro-cogeneration system, different assessment strategies are presented. Furthermore, a case study for Singapore is conducted to see the impact of the forecasted carbon dioxide prices on the thermoeconomic performance of the micro-cogeneration system.

Mutual Information Based Dynamic Integration of Multiple Feature Streams for Robust Real-Time LVCSR

Science.gov (United States)

Sato, Shoei; Kobayashi, Akio; Onoe, Kazuo; Homma, Shinichi; Imai, Toru; Takagi, Tohru; Kobayashi, Tetsunori

We present a novel method of integrating the likelihoods of multiple feature streams, representing different acoustic aspects, for robust speech recognition. The integration algorithm dynamically calculates a frame-wise stream weight so that a higher weight is given to a stream that is robust to a variety of noisy environments or speaking styles. Such a robust stream is expected to show discriminative ability. A conventional method proposed for the recognition of spoken digits calculates the weights front the entropy of the whole set of HMM states. This paper extends the dynamic weighting to a real-time large-vocabulary continuous speech recognition (LVCSR) system. The proposed weight is calculated in real-time from mutual information between an input stream and active HMM states in a searchs pace without an additional likelihood calculation. Furthermore, the mutual information takes the width of the search space into account by calculating the marginal entropy from the number of active states. In this paper, we integrate three features that are extracted through auditory filters by taking into account the human auditory system's ability to extract amplitude and frequency modulations. Due to this, features representing energy, amplitude drift, and resonant frequency drifts, are integrated. These features are expected to provide complementary clues for speech recognition. Speech recognition experiments on field reports and spontaneous commentary from Japanese broadcast news showed that the proposed method reduced error words by 9.2% in field reports and 4.7% in spontaneous commentaries relative to the best result obtained from a single stream.

The Green's matrix and the boundary integral equations for analysis of time-harmonic dynamics of elastic helical springs.

Science.gov (United States)

Sorokin, Sergey V

2011-03-01

Helical springs serve as vibration isolators in virtually any suspension system. Various exact and approximate methods may be employed to determine the eigenfrequencies of vibrations of these structural elements and their dynamic transfer functions. The method of boundary integral equations is a meaningful alternative to obtain exact solutions of problems of the time-harmonic dynamics of elastic springs in the framework of Bernoulli-Euler beam theory. In this paper, the derivations of the Green's matrix, of the Somigliana's identities, and of the boundary integral equations are presented. The vibrational power transmission in an infinitely long spring is analyzed by means of the Green's matrix. The eigenfrequencies and the dynamic transfer functions are found by solving the boundary integral equations. In the course of analysis, the essential features and advantages of the method of boundary integral equations are highlighted. The reported analytical results may be used to study the time-harmonic motion in any wave guide governed by a system of linear differential equations in a single spatial coordinate along its axis. © 2011 Acoustical Society of America

Approximations to the Probability of Failure in Random Vibration by Integral Equation Methods

DEFF Research Database (Denmark)

Nielsen, Søren R.K.; Sørensen, John Dalsgaard

Close approximations to the first passage probability of failure in random vibration can be obtained by integral equation methods. A simple relation exists between the first passage probability density function and the distribution function for the time interval spent below a barrier before...... passage probability density. The results of the theory agree well with simulation results for narrow banded processes dominated by a single frequency, as well as for bimodal processes with 2 dominating frequencies in the structural response....... outcrossing. An integral equation for the probability density function of the time interval is formulated, and adequate approximations for the kernel are suggested. The kernel approximation results in approximate solutions for the probability density function of the time interval, and hence for the first...

Mixed methods in psychotherapy research: A review of method(ology) integration in psychotherapy science.

Science.gov (United States)

Bartholomew, Theodore T; Lockard, Allison J

2018-06-13

Mixed methods can foster depth and breadth in psychological research. However, its use remains in development in psychotherapy research. Our purpose was to review the use of mixed methods in psychotherapy research. Thirty-one studies were identified via the PRISMA systematic review method. Using Creswell & Plano Clark's typologies to identify design characteristics, we assessed each study for rigor and how each used mixed methods. Key features of mixed methods designs and these common patterns were identified: (a) integration of clients' perceptions via mixing; (b) understanding group psychotherapy; (c) integrating methods with cases and small samples; (d) analyzing clinical data as qualitative data; and (e) exploring cultural identities in psychotherapy through mixed methods. The review is discussed with respect to the value of integrating multiple data in single studies to enhance psychotherapy research. © 2018 Wiley Periodicals, Inc.

Methods for enhancing numerical integration

International Nuclear Information System (INIS)

Doncker, Elise de

2003-01-01

We give a survey of common strategies for numerical integration (adaptive, Monte-Carlo, Quasi-Monte Carlo), and attempt to delineate their realm of applicability. The inherent accuracy and error bounds for basic integration methods are given via such measures as the degree of precision of cubature rules, the index of a family of lattice rules, and the discrepancy of uniformly distributed point sets. Strategies incorporating these basic methods often use paradigms to reduce the error by, e.g., increasing the number of points in the domain or decreasing the mesh size, locally or uniformly. For these processes the order of convergence of the strategy is determined by the asymptotic behavior of the error, and may be too slow in practice for the type of problem at hand. For certain problem classes we may be able to improve the effectiveness of the method or strategy by such techniques as transformations, absorbing a difficult part of the integrand into a weight function, suitable partitioning of the domain, transformations and extrapolation or convergence acceleration. Situations warranting the use of these techniques (possibly in an 'automated' way) are described and illustrated by sample applications

Verifying Real-Time Systems using Explicit-time Description Methods

Directory of Open Access Journals (Sweden)

Hao Wang

2009-12-01

Full Text Available Timed model checking has been extensively researched in recent years. Many new formalisms with time extensions and tools based on them have been presented. On the other hand, Explicit-Time Description Methods aim to verify real-time systems with general untimed model checkers. Lamport presented an explicit-time description method using a clock-ticking process (Tick to simulate the passage of time together with a group of global variables for time requirements. This paper proposes a new explicit-time description method with no reliance on global variables. Instead, it uses rendezvous synchronization steps between the Tick process and each system process to simulate time. This new method achieves better modularity and facilitates usage of more complex timing constraints. The two explicit-time description methods are implemented in DIVINE, a well-known distributed-memory model checker. Preliminary experiment results show that our new method, with better modularity, is comparable to Lamport's method with respect to time and memory efficiency.

Time resolved spectroscopy of GRB 030501 using INTEGRAL

DEFF Research Database (Denmark)

Beckmann, V.; Borkowski, J.; Courvoisier, T.J.L.

2003-01-01

The gamma-ray instruments on-board INTEGRAL offer an unique opportunity to perform time resolved analysis on GRBs. The imager IBIS allows accurate positioning of GRBs and broad band spectral analysis, while SPI provides high resolution spectroscopy. GRB 030501 was discovered by the INTEGRAL Burst...... the Ulysses and RHESSI experiments....

Momentum integral network method for thermal-hydraulic transient analysis

International Nuclear Information System (INIS)

Van Tuyle, G.J.

1983-01-01

A new momentum integral network method has been developed, and tested in the MINET computer code. The method was developed in order to facilitate the transient analysis of complex fluid flow and heat transfer networks, such as those found in the balance of plant of power generating facilities. The method employed in the MINET code is a major extension of a momentum integral method reported by Meyer. Meyer integrated the momentum equation over several linked nodes, called a segment, and used a segment average pressure, evaluated from the pressures at both ends. Nodal mass and energy conservation determined nodal flows and enthalpies, accounting for fluid compression and thermal expansion

Integrals of Frullani type and the method of brackets

Directory of Open Access Journals (Sweden)

Bravo Sergio

2017-01-01

Full Text Available The method of brackets is a collection of heuristic rules, some of which have being made rigorous, that provide a flexible, direct method for the evaluation of definite integrals. The present work uses this method to establish classical formulas due to Frullani which provide values of a specific family of integrals. Some generalizations are established.

Analysis of factors influencing the integrated bolus peak timing in contrast-enhanced brain computed tomographic angiography

International Nuclear Information System (INIS)

Son, Soon Yong; Choi, Kwan Woo; Jeong, Hoi Woun; Jang, Seo Goo; Jung, Jae Young; Yun, Jung Soo; Kim, Ki Won; Lee, Young Ah; Son, Jin Hyun; Min, Jung Whan

2016-01-01

The objective of this study was to analyze the factors influencing integrated bolus peak timing in contrast- enhanced computed tomographic angiography (CTA) and to determine a method of calculating personal peak time. The optimal time was calculated by performing multiple linear regression analysis, after finding the influence factors through correlation analysis between integrated peak time of contrast medium and personal measured value by monitoring CTA scans. The radiation exposure dose in CTA was 716.53 mGy·cm and the radiation exposure dose in monitoring scan was 15.52 mGy (2 - 34 mGy). The results were statistically significant (p < .01). Regression analysis revealed, a -0.160 times decrease with a one-step increase in heart rate in male, and -0.004, -0.174, and 0.006 times decrease with one-step in DBP, heart rate, and blood sugar, respectively, in female. In a consistency test of peak time by calculating measured peak time and peak time by using the regression equation, the consistency was determined to be very high for male and female. This study could prevent unnecessary dose exposure by encouraging in clinic calculation of personal integrated peak time of contrast medium prior to examination

Computational electrodynamics the finite-difference time-domain method

CERN Document Server

Taflove, Allen

2005-01-01

This extensively revised and expanded third edition of the Artech House bestseller, Computational Electrodynamics: The Finite-Difference Time-Domain Method, offers engineers the most up-to-date and definitive resource on this critical method for solving Maxwell's equations. The method helps practitioners design antennas, wireless communications devices, high-speed digital and microwave circuits, and integrated optical devices with unsurpassed efficiency. There has been considerable advancement in FDTD computational technology over the past few years, and the third edition brings professionals the very latest details with entirely new chapters on important techniques, major updates on key topics, and new discussions on emerging areas such as nanophotonics. What's more, to supplement the third edition, the authors have created a Web site with solutions to problems, downloadable graphics and videos, and updates, making this new edition the ideal textbook on the subject as well.

Approximation of the exponential integral (well function) using sampling methods

Science.gov (United States)

Baalousha, Husam Musa

2015-04-01

Exponential integral (also known as well function) is often used in hydrogeology to solve Theis and Hantush equations. Many methods have been developed to approximate the exponential integral. Most of these methods are based on numerical approximations and are valid for a certain range of the argument value. This paper presents a new approach to approximate the exponential integral. The new approach is based on sampling methods. Three different sampling methods; Latin Hypercube Sampling (LHS), Orthogonal Array (OA), and Orthogonal Array-based Latin Hypercube (OA-LH) have been used to approximate the function. Different argument values, covering a wide range, have been used. The results of sampling methods were compared with results obtained by Mathematica software, which was used as a benchmark. All three sampling methods converge to the result obtained by Mathematica, at different rates. It was found that the orthogonal array (OA) method has the fastest convergence rate compared with LHS and OA-LH. The root mean square error RMSE of OA was in the order of 1E-08. This method can be used with any argument value, and can be used to solve other integrals in hydrogeology such as the leaky aquifer integral.

A high-order relaxation method with projective integration for solving nonlinear systems of hyperbolic conservation laws

Science.gov (United States)

Lafitte, Pauline; Melis, Ward; Samaey, Giovanni

2017-07-01

We present a general, high-order, fully explicit relaxation scheme which can be applied to any system of nonlinear hyperbolic conservation laws in multiple dimensions. The scheme consists of two steps. In a first (relaxation) step, the nonlinear hyperbolic conservation law is approximated by a kinetic equation with stiff BGK source term. Then, this kinetic equation is integrated in time using a projective integration method. After taking a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution, the time derivative is estimated and used in an (outer) Runge-Kutta method of arbitrary order. We show that, with an appropriate choice of inner step size, the time step restriction on the outer time step is similar to the CFL condition for the hyperbolic conservation law. Moreover, the number of inner time steps is also independent of the stiffness of the BGK source term. We discuss stability and consistency, and illustrate with numerical results (linear advection, Burgers' equation and the shallow water and Euler equations) in one and two spatial dimensions.

Development of an integrated method for long-term water quality prediction using seasonal climate forecast

Directory of Open Access Journals (Sweden)

J. Cho

2016-10-01

Full Text Available The APEC Climate Center (APCC produces climate prediction information utilizing a multi-climate model ensemble (MME technique. In this study, four different downscaling methods, in accordance with the degree of utilizing the seasonal climate prediction information, were developed in order to improve predictability and to refine the spatial scale. These methods include: (1 the Simple Bias Correction (SBC method, which directly uses APCC's dynamic prediction data with a 3 to 6 month lead time; (2 the Moving Window Regression (MWR method, which indirectly utilizes dynamic prediction data; (3 the Climate Index Regression (CIR method, which predominantly uses observation-based climate indices; and (4 the Integrated Time Regression (ITR method, which uses predictors selected from both CIR and MWR. Then, a sampling-based temporal downscaling was conducted using the Mahalanobis distance method in order to create daily weather inputs to the Soil and Water Assessment Tool (SWAT model. Long-term predictability of water quality within the Wecheon watershed of the Nakdong River Basin was evaluated. According to the Korean Ministry of Environment's Provisions of Water Quality Prediction and Response Measures, modeling-based predictability was evaluated by using 3-month lead prediction data issued in February, May, August, and November as model input of SWAT. Finally, an integrated approach, which takes into account various climate information and downscaling methods for water quality prediction, was presented. This integrated approach can be used to prevent potential problems caused by extreme climate in advance.

Molecular radiotherapy: the NUKFIT software for calculating the time-integrated activity coefficient.

Science.gov (United States)

Kletting, P; Schimmel, S; Kestler, H A; Hänscheid, H; Luster, M; Fernández, M; Bröer, J H; Nosske, D; Lassmann, M; Glatting, G

2013-10-01

Calculation of the time-integrated activity coefficient (residence time) is a crucial step in dosimetry for molecular radiotherapy. However, available software is deficient in that it is either not tailored for the use in molecular radiotherapy and/or does not include all required estimation methods. The aim of this work was therefore the development and programming of an algorithm which allows for an objective and reproducible determination of the time-integrated activity coefficient and its standard error. The algorithm includes the selection of a set of fitting functions from predefined sums of exponentials and the choice of an error model for the used data. To estimate the values of the adjustable parameters an objective function, depending on the data, the parameters of the error model, the fitting function and (if required and available) Bayesian information, is minimized. To increase reproducibility and user-friendliness the starting values are automatically determined using a combination of curve stripping and random search. Visual inspection, the coefficient of determination, the standard error of the fitted parameters, and the correlation matrix are provided to evaluate the quality of the fit. The functions which are most supported by the data are determined using the corrected Akaike information criterion. The time-integrated activity coefficient is estimated by analytically integrating the fitted functions. Its standard error is determined assuming Gaussian error propagation. The software was implemented using MATLAB. To validate the proper implementation of the objective function and the fit functions, the results of NUKFIT and SAAM numerical, a commercially available software tool, were compared. The automatic search for starting values was successfully tested for reproducibility. The quality criteria applied in conjunction with the Akaike information criterion allowed the selection of suitable functions. Function fit parameters and their standard

Nodal integral method for the neutron diffusion equation in cylindrical geometry

International Nuclear Information System (INIS)

Azmy, Y.Y.

1987-01-01

The nodal methodology is based on retaining a higher a higher degree of analyticity in the process of deriving the discrete-variable equations compared to conventional numerical methods. As a result, extensive numerical testing of nodal methods developed for a wide variety of partial differential equations and comparison of the results to conventional methods have established the superior accuracy of nodal methods on coarse meshes. Moreover, these tests have shown that nodal methods are more computationally efficient than finite difference and finite-element methods in the sense that they require shorter CPU times to achieve comparable accuracy in the solutions. However, nodal formalisms and the final discrete-variable equations they produce are, in general, more complicated than their conventional counterparts. This, together with anticipated difficulties in applying the transverse-averaging procedure in curvilinear coordinates, has limited the applications of nodal methods, so far, to Cartesian geometry, and with additional approximations to hexagonal geometry. In this paper the authors report recent progress in deriving and numerically implementing a nodal integral method (NIM) for solving the neutron diffusion equation in cylindrical r-z geometry. Also, presented are comparisons of numerical solutions to two test problems with those obtained by the Exterminator-2 code, which indicate the superior accuracy of the nodal integral method solutions on much coarser meshes

A new approach to calculate charge carrier transport mobility in organic molecular crystals from imaginary time path integral simulations

International Nuclear Information System (INIS)

Song, Linze; Shi, Qiang

2015-01-01

We present a new non-perturbative method to calculate the charge carrier mobility using the imaginary time path integral approach, which is based on the Kubo formula for the conductivity, and a saddle point approximation to perform the analytic continuation. The new method is first tested using a benchmark calculation from the numerical exact hierarchical equations of motion method. Imaginary time path integral Monte Carlo simulations are then performed to explore the temperature dependence of charge carrier delocalization and mobility in organic molecular crystals (OMCs) within the Holstein and Holstein-Peierls models. The effects of nonlocal electron-phonon interaction on mobility in different charge transport regimes are also investigated

Explicit solution of Calderon preconditioned time domain integral equations

KAUST Repository

Ulku, Huseyin Arda

2013-07-01

An explicit marching on-in-time (MOT) scheme for solving Calderon-preconditioned time domain integral equations is proposed. The scheme uses Rao-Wilton-Glisson and Buffa-Christiansen functions to discretize the domain and range of the integral operators and a PE(CE)m type linear multistep to march on in time. Unlike its implicit counterpart, the proposed explicit solver requires the solution of an MOT system with a Gram matrix that is sparse and well-conditioned independent of the time step size. Numerical results demonstrate that the explicit solver maintains its accuracy and stability even when the time step size is chosen as large as that typically used by an implicit solver. © 2013 IEEE.

Explicit solution of the time domain volume integral equation using a stable predictor-corrector scheme

KAUST Repository

Al Jarro, Ahmed

2012-11-01

An explicit marching-on-in-time (MOT) scheme for solving the time domain volume integral equation is presented. The proposed method achieves its stability by employing, at each time step, a corrector scheme, which updates/corrects fields computed by the explicit predictor scheme. The proposedmethod is computationally more efficient when compared to the existing filtering techniques used for the stabilization of explicit MOT schemes. Numerical results presented in this paper demonstrate that the proposed method maintains its stability even when applied to the analysis of electromagnetic wave interactions with electrically large structures meshed using approximately half a million discretization elements.

Explicit solution of the time domain volume integral equation using a stable predictor-corrector scheme

KAUST Repository

Al Jarro, Ahmed; Salem, Mohamed; Bagci, Hakan; Benson, Trevor; Sewell, Phillip D.; Vuković, Ana

2012-01-01

An explicit marching-on-in-time (MOT) scheme for solving the time domain volume integral equation is presented. The proposed method achieves its stability by employing, at each time step, a corrector scheme, which updates/corrects fields computed by the explicit predictor scheme. The proposedmethod is computationally more efficient when compared to the existing filtering techniques used for the stabilization of explicit MOT schemes. Numerical results presented in this paper demonstrate that the proposed method maintains its stability even when applied to the analysis of electromagnetic wave interactions with electrically large structures meshed using approximately half a million discretization elements.

Real-time functional integral approach to the quantum disordered spin systems

International Nuclear Information System (INIS)

Kopec, T.K.

1989-01-01

In this paper the effect of randomness and frustration in the quantum Ising spin glass in a transverse field is studied by using the thermofield dynamics (TFD), the real time, finite temperature quantum field theory. It is shown that the method can be conveniently used for the averaging of the free energy of the system by completely avoiding the use of the n-replica trick. The effective dynamic Lagrangian for the disorder averaged causal, correlations and response Green functions is derived by functional integral approach. Furthermore, the properties of this Lagrangian are analyzed by the saddle point method which leads to the self-consistent equation for the spin glass order parameter

System and method for integrating hazard-based decision making tools and processes

Science.gov (United States)

Hodgin, C Reed [Westminster, CO

2012-03-20

A system and method for inputting, analyzing, and disseminating information necessary for identified decision-makers to respond to emergency situations. This system and method provides consistency and integration among multiple groups, and may be used for both initial consequence-based decisions and follow-on consequence-based decisions. The system and method in a preferred embodiment also provides tools for accessing and manipulating information that are appropriate for each decision-maker, in order to achieve more reasoned and timely consequence-based decisions. The invention includes processes for designing and implementing a system or method for responding to emergency situations.

A CMOS integrated timing discriminator circuit for fast scintillation counters

International Nuclear Information System (INIS)

Jochmann, M.W.

1998-01-01

Based on a zero-crossing discriminator using a CR differentiation network for pulse shaping, a new CMOS integrated timing discriminator circuit is proposed for fast (t r ≥ 2 ns) scintillation counters at the cooler synchrotron COSY-Juelich. By eliminating the input signal's amplitude information by means of an analog continuous-time divider, a normalized pulse shape at the zero-crossing point is gained over a wide dynamic input amplitude range. In combination with an arming comparator and a monostable multivibrator this yields in a highly precise timing discriminator circuit, that is expected to be useful in different time measurement applications. First measurement results of a CMOS integrated logarithmic amplifier, which is part of the analog continuous-time divider, agree well with the corresponding simulations. Moreover, SPICE simulations of the integrated discriminator circuit promise a time walk well below 200 ps (FWHM) over a 40 dB input amplitude dynamic range

An integrated technique for developing real-time systems

NARCIS (Netherlands)

Hooman, J.J.M.; Vain, J.

1995-01-01

The integration of conceptual modeling techniques, formal specification, and compositional verification is considered for real time systems within the knowledge engineering context. We define constructive transformations from a conceptual meta model to a real time specification language and give

Direct integration multiple collision integral transport analysis method for high energy fusion neutronics

International Nuclear Information System (INIS)

Koch, K.R.

1985-01-01

A new analysis method specially suited for the inherent difficulties of fusion neutronics was developed to provide detailed studies of the fusion neutron transport physics. These studies should provide a better understanding of the limitations and accuracies of typical fusion neutronics calculations. The new analysis method is based on the direct integration of the integral form of the neutron transport equation and employs a continuous energy formulation with the exact treatment of the energy angle kinematics of the scattering process. In addition, the overall solution is analyzed in terms of uncollided, once-collided, and multi-collided solution components based on a multiple collision treatment. Furthermore, the numerical evaluations of integrals use quadrature schemes that are based on the actual dependencies exhibited in the integrands. The new DITRAN computer code was developed on the Cyber 205 vector supercomputer to implement this direct integration multiple-collision fusion neutronics analysis. Three representative fusion reactor models were devised and the solutions to these problems were studied to provide suitable choices for the numerical quadrature orders as well as the discretized solution grid and to understand the limitations of the new analysis method. As further verification and as a first step in assessing the accuracy of existing fusion-neutronics calculations, solutions obtained using the new analysis method were compared to typical multigroup discrete ordinates calculations

Achieving Integration in Mixed Methods Designs—Principles and Practices

OpenAIRE

Fetters, Michael D; Curry, Leslie A; Creswell, John W

2013-01-01

Mixed methods research offers powerful tools for investigating complex processes and systems in health and health care. This article describes integration principles and practices at three levels in mixed methods research and provides illustrative examples. Integration at the study design level occurs through three basic mixed method designs—exploratory sequential, explanatory sequential, and convergent—and through four advanced frameworks—multistage, intervention, case study, and participato...

Characterization of HBV integration patterns and timing in liver cancer and HBV-infected livers.

Science.gov (United States)

Furuta, Mayuko; Tanaka, Hiroko; Shiraishi, Yuichi; Unida, Takuro; Imamura, Michio; Fujimoto, Akihiro; Fujita, Masahi; Sasaki-Oku, Aya; Maejima, Kazuhiro; Nakano, Kaoru; Kawakami, Yoshiiku; Arihiro, Koji; Aikata, Hiroshi; Ueno, Masaki; Hayami, Shinya; Ariizumi, Shun-Ichi; Yamamoto, Masakazu; Gotoh, Kunihito; Ohdan, Hideki; Yamaue, Hiroki; Miyano, Satoru; Chayama, Kazuaki; Nakagawa, Hidewaki

2018-05-18

Integration of Hepatitis B virus (HBV) into the human genome can cause genetic instability, leading to selective advantages for HBV-induced liver cancer. Despite the large number of studies for HBV integration into liver cancer, little is known about the mechanism of initial HBV integration events owing to the limitations of materials and detection methods. We conducted an HBV sequence capture, followed by ultra-deep sequencing, to screen for HBV integrations in 111 liver samples from human-hepatocyte chimeric mice with HBV infection and human clinical samples containing 42 paired samples from non-tumorous and tumorous liver tissues. The HBV infection model using chimeric mice verified the efficiency of our HBV-capture analysis and demonstrated that HBV integration could occur 23 to 49 days after HBV infection via microhomology-mediated end joining and predominantly in mitochondrial DNA. Overall HBV integration sites in clinical samples were significantly enriched in regions annotated as exhibiting open chromatin, a high level of gene expression, and early replication timing in liver cells. These data indicate that HBV integration in liver tissue was biased according to chromatin accessibility, with additional selection pressures in the gene promoters of tumor samples. Moreover, an integrative analysis using paired non-tumorous and tumorous samples and HBV-related transcriptional change revealed the involvement of TERT and MLL4 in clonal selection. We also found frequent and non-tumorous liver-specific HBV integrations in FN1 and HBV-FN1 fusion transcript. Extensive survey of HBV integrations facilitates and improves the understanding of the timing and biology of HBV integration during infection and HBV-related hepatocarcinogenesis.

A method for untriggered time-dependent searches for multiple flares from neutrino point sources

International Nuclear Information System (INIS)

Gora, D.; Bernardini, E.; Cruz Silva, A.H.

2011-04-01

A method for a time-dependent search for flaring astrophysical sources which can be potentially detected by large neutrino experiments is presented. The method uses a time-clustering algorithm combined with an unbinned likelihood procedure. By including in the likelihood function a signal term which describes the contribution of many small clusters of signal-like events, this method provides an effective way for looking for weak neutrino flares over different time-scales. The method is sensitive to an overall excess of events distributed over several flares which are not individually detectable. For standard cases (one flare) the discovery potential of the method is worse than a standard time-dependent point source analysis with unknown duration of the flare by a factor depending on the signal-to-background level. However, for flares sufficiently shorter than the total observation period, the method is more sensitive than a time-integrated analysis. (orig.)

A method for untriggered time-dependent searches for multiple flares from neutrino point sources

Energy Technology Data Exchange (ETDEWEB)

Gora, D. [Deutsches Elektronen-Synchrotron (DESY), Zeuthen (Germany); Institute of Nuclear Physics PAN, Cracow (Poland); Bernardini, E.; Cruz Silva, A.H. [Institute of Nuclear Physics PAN, Cracow (Poland)

2011-04-15

A method for a time-dependent search for flaring astrophysical sources which can be potentially detected by large neutrino experiments is presented. The method uses a time-clustering algorithm combined with an unbinned likelihood procedure. By including in the likelihood function a signal term which describes the contribution of many small clusters of signal-like events, this method provides an effective way for looking for weak neutrino flares over different time-scales. The method is sensitive to an overall excess of events distributed over several flares which are not individually detectable. For standard cases (one flare) the discovery potential of the method is worse than a standard time-dependent point source analysis with unknown duration of the flare by a factor depending on the signal-to-background level. However, for flares sufficiently shorter than the total observation period, the method is more sensitive than a time-integrated analysis. (orig.)

An Efficient Explicit-time Description Method for Timed Model Checking

Directory of Open Access Journals (Sweden)

Hao Wang

2009-12-01

Full Text Available Timed model checking, the method to formally verify real-time systems, is attracting increasing attention from both the model checking community and the real-time community. Explicit-time description methods verify real-time systems using general model constructs found in standard un-timed model checkers. Lamport proposed an explicit-time description method using a clock-ticking process (Tick to simulate the passage of time together with a group of global variables to model time requirements. Two methods, the Sync-based Explicit-time Description Method using rendezvous synchronization steps and the Semaphore-based Explicit-time Description Method using only one global variable were proposed; they both achieve better modularity than Lamport's method in modeling the real-time systems. In contrast to timed automata based model checkers like UPPAAL, explicit-time description methods can access and store the current time instant for future calculations necessary for many real-time systems, especially those with pre-emptive scheduling. However, the Tick process in the above three methods increments the time by one unit in each tick; the state spaces therefore grow relatively fast as the time parameters increase, a problem when the system's time period is relatively long. In this paper, we propose a more efficient method which enables the Tick process to leap multiple time units in one tick. Preliminary experimental results in a high performance computing environment show that this new method significantly reduces the state space and improves both the time and memory efficiency.

Numerical Simulation of Antennas with Improved Integral Equation Method

International Nuclear Information System (INIS)

Ma Ji; Fang Guang-You; Lu Wei

2015-01-01

Simulating antennas around a conducting object is a challenge task in computational electromagnetism, which is concerned with the behaviour of electromagnetic fields. To analyze this model efficiently, an improved integral equation-fast Fourier transform (IE-FFT) algorithm is presented in this paper. The proposed scheme employs two Cartesian grids with different size and location to enclose the antenna and the other object, respectively. On the one hand, IE-FFT technique is used to store matrix in a sparse form and accelerate the matrix-vector multiplication for each sub-domain independently. On the other hand, the mutual interaction between sub-domains is taken as the additional exciting voltage in each matrix equation. By updating integral equations several times, the whole electromagnetic system can achieve a stable status. Finally, the validity of the presented method is verified through the analysis of typical antennas in the presence of a conducting object. (paper)

Differential amplicons (ΔAmp)-a new molecular method to assess RNA integrity.

Science.gov (United States)

Björkman, J; Švec, D; Lott, E; Kubista, M; Sjöback, R

2016-01-01

Integrity of the mRNA in clinical samples has major impact on the quality of measured expression levels. This is independent of the measurement technique being next generation sequencing (NGS), Quantitative real-time PCR (qPCR) or microarray profiling. If mRNA is highly degraded or damaged, measured data will be very unreliable and the whole study is likely a waste of time and money. It is therefore common strategy to test the quality of RNA in samples before conducting large and costly studies. Most methods today to assess the quality of RNA are ignorant to the nature of the RNA and, therefore, reflect the integrity of ribosomal RNA, which is the dominant species, rather than of mRNAs, microRNAs and long non-coding RNAs, which usually are the species of interest. Here, we present a novel molecular approach to assess the quality of the targeted RNA species by measuring the differential amplification (ΔAmp) of an Endogenous RNase Resistant (ERR) marker relative to a reference gene, optionally combined with the measurement of two amplicons of different lengths. The combination reveals any mRNA degradation caused by ribonucleases as well as physical, chemical or UV damage. ΔAmp has superior sensitivity to common microfluidic electrophoretic methods, senses the integrity of the actual targeted RNA species, and allows for a smoother and more cost efficient workflow.

Differential amplicons (ΔAmp—a new molecular method to assess RNA integrity

Directory of Open Access Journals (Sweden)

J. Björkman

2016-01-01

Full Text Available Integrity of the mRNA in clinical samples has major impact on the quality of measured expression levels. This is independent of the measurement technique being next generation sequencing (NGS, Quantitative real-time PCR (qPCR or microarray profiling. If mRNA is highly degraded or damaged, measured data will be very unreliable and the whole study is likely a waste of time and money. It is therefore common strategy to test the quality of RNA in samples before conducting large and costly studies. Most methods today to assess the quality of RNA are ignorant to the nature of the RNA and, therefore, reflect the integrity of ribosomal RNA, which is the dominant species, rather than of mRNAs, microRNAs and long non-coding RNAs, which usually are the species of interest. Here, we present a novel molecular approach to assess the quality of the targeted RNA species by measuring the differential amplification (ΔAmp of an Endogenous RNase Resistant (ERR marker relative to a reference gene, optionally combined with the measurement of two amplicons of different lengths. The combination reveals any mRNA degradation caused by ribonucleases as well as physical, chemical or UV damage. ΔAmp has superior sensitivity to common microfluidic electrophoretic methods, senses the integrity of the actual targeted RNA species, and allows for a smoother and more cost efficient workflow.

Explicit solution of the time domain magnetic field integral equation using a predictor-corrector scheme

KAUST Repository

Ulku, Huseyin Arda; Bagci, Hakan; Michielssen, Eric

2012-01-01

An explicit yet stable marching-on-in-time (MOT) scheme for solving the time domain magnetic field integral equation (TD-MFIE) is presented. The stability of the explicit scheme is achieved via (i) accurate evaluation of the MOT matrix elements using closed form expressions and (ii) a PE(CE) m type linear multistep method for time marching. Numerical results demonstrate the accuracy and stability of the proposed explicit MOT-TD-MFIE solver. © 2012 IEEE.

Explicit solution of the time domain magnetic field integral equation using a predictor-corrector scheme

KAUST Repository

Ulku, Huseyin Arda

2012-09-01

An explicit yet stable marching-on-in-time (MOT) scheme for solving the time domain magnetic field integral equation (TD-MFIE) is presented. The stability of the explicit scheme is achieved via (i) accurate evaluation of the MOT matrix elements using closed form expressions and (ii) a PE(CE) m type linear multistep method for time marching. Numerical results demonstrate the accuracy and stability of the proposed explicit MOT-TD-MFIE solver. © 2012 IEEE.

Effects of integration time on in-water radiometric profiles.

Science.gov (United States)

D'Alimonte, Davide; Zibordi, Giuseppe; Kajiyama, Tamito

2018-03-05

This work investigates the effects of integration time on in-water downward irradiance E d , upward irradiance E u and upwelling radiance L u profile data acquired with free-fall hyperspectral systems. Analyzed quantities are the subsurface value and the diffuse attenuation coefficient derived by applying linear and non-linear regression schemes. Case studies include oligotrophic waters (Case-1), as well as waters dominated by Colored Dissolved Organic Matter (CDOM) and Non-Algal Particles (NAP). Assuming a 24-bit digitization, measurements resulting from the accumulation of photons over integration times varying between 8 and 2048ms are evaluated at depths corresponding to: 1) the beginning of each integration interval (Fst); 2) the end of each integration interval (Lst); 3) the averages of Fst and Lst values (Avg); and finally 4) the values weighted accounting for the diffuse attenuation coefficient of water (Wgt). Statistical figures show that the effects of integration time can bias results well above 5% as a function of the depth definition. Results indicate the validity of the Wgt depth definition and the fair applicability of the Avg one. Instead, both the Fst and Lst depths should not be adopted since they may introduce pronounced biases in E u and L u regression products for highly absorbing waters. Finally, the study reconfirms the relevance of combining multiple radiometric casts into a single profile to increase precision of regression products.

Analysis of electromagnetic wave interactions on nonlinear scatterers using time domain volume integral equations

KAUST Repository

Ulku, Huseyin Arda

2014-07-06

Effects of material nonlinearities on electromagnetic field interactions become dominant as field amplitudes increase. A typical example is observed in plasmonics, where highly localized fields “activate” Kerr nonlinearities. Naturally, time domain solvers are the method of choice when it comes simulating these nonlinear effects. Oftentimes, finite difference time domain (FDTD) method is used for this purpose. This is simply due to the fact that explicitness of the FDTD renders the implementation easier and the material nonlinearity can be easily accounted for using an auxiliary differential equation (J.H. Green and A. Taflove, Opt. Express, 14(18), 8305-8310, 2006). On the other hand, explicit marching on-in-time (MOT)-based time domain integral equation (TDIE) solvers have never been used for the same purpose even though they offer several advantages over FDTD (E. Michielssen, et al., ECCOMAS CFD, The Netherlands, Sep. 5-8, 2006). This is because explicit MOT solvers have never been stabilized until not so long ago. Recently an explicit but stable MOT scheme has been proposed for solving the time domain surface magnetic field integral equation (H.A. Ulku, et al., IEEE Trans. Antennas Propag., 61(8), 4120-4131, 2013) and later it has been extended for the time domain volume electric field integral equation (TDVEFIE) (S. B. Sayed, et al., Pr. Electromagn. Res. S., 378, Stockholm, 2013). This explicit MOT scheme uses predictor-corrector updates together with successive over relaxation during time marching to stabilize the solution even when time step is as large as in the implicit counterpart. In this work, an explicit MOT-TDVEFIE solver is proposed for analyzing electromagnetic wave interactions on scatterers exhibiting Kerr nonlinearity. Nonlinearity is accounted for using the constitutive relation between the electric field intensity and flux density. Then, this relation and the TDVEFIE are discretized together by expanding the intensity and flux - sing half

[Development method of healthcare information system integration based on business collaboration model].

Science.gov (United States)

Li, Shasha; Nie, Hongchao; Lu, Xudong; Duan, Huilong

2015-02-01

Integration of heterogeneous systems is the key to hospital information construction due to complexity of the healthcare environment. Currently, during the process of healthcare information system integration, people participating in integration project usually communicate by free-format document, which impairs the efficiency and adaptability of integration. A method utilizing business process model and notation (BPMN) to model integration requirement and automatically transforming it to executable integration configuration was proposed in this paper. Based on the method, a tool was developed to model integration requirement and transform it to integration configuration. In addition, an integration case in radiology scenario was used to verify the method.

Explicit integration of extremely stiff reaction networks: partial equilibrium methods

International Nuclear Information System (INIS)

Guidry, M W; Hix, W R; Billings, J J

2013-01-01

In two preceding papers (Guidry et al 2013 Comput. Sci. Disc. 6 015001 and Guidry and Harris 2013 Comput. Sci. Disc. 6 015002), we have shown that when reaction networks are well removed from equilibrium, explicit asymptotic and quasi-steady-state approximations can give algebraically stabilized integration schemes that rival standard implicit methods in accuracy and speed for extremely stiff systems. However, we also showed that these explicit methods remain accurate but are no longer competitive in speed as the network approaches equilibrium. In this paper, we analyze this failure and show that it is associated with the presence of fast equilibration timescales that neither asymptotic nor quasi-steady-state approximations are able to remove efficiently from the numerical integration. Based on this understanding, we develop a partial equilibrium method to deal effectively with the approach to equilibrium and show that explicit asymptotic methods, combined with the new partial equilibrium methods, give an integration scheme that can plausibly deal with the stiffest networks, even in the approach to equilibrium, with accuracy and speed competitive with that of implicit methods. Thus we demonstrate that such explicit methods may offer alternatives to implicit integration of even extremely stiff systems and that these methods may permit integration of much larger networks than have been possible before in a number of fields. (paper)

Optimal Real-time Dispatch for Integrated Energy Systems

DEFF Research Database (Denmark)

Anvari-Moghaddam, Amjad; Guerrero, Josep M.; Rahimi-Kian, Ashkan

2016-01-01

With the emerging of small-scale integrated energy systems (IESs), there are significant potentials to increase the functionality of a typical demand-side management (DSM) strategy and typical implementation of building-level distributed energy resources (DERs). By integrating DSM and DERs...... into a cohesive, networked package that fully utilizes smart energy-efficient end-use devices, advanced building control/automation systems, and integrated communications architectures, it is possible to efficiently manage energy and comfort at the end-use location. In this paper, an ontology-driven multi......-agent control system with intelligent optimizers is proposed for optimal real-time dispatch of an integrated building and microgrid system considering coordinated demand response (DR) and DERs management. The optimal dispatch problem is formulated as a mixed integer nonlinear programing problem (MINLP...

A simple reliability block diagram method for safety integrity verification

International Nuclear Information System (INIS)

Guo Haitao; Yang Xianhui

2007-01-01

IEC 61508 requires safety integrity verification for safety related systems to be a necessary procedure in safety life cycle. PFD avg must be calculated to verify the safety integrity level (SIL). Since IEC 61508-6 does not give detailed explanations of the definitions and PFD avg calculations for its examples, it is difficult for common reliability or safety engineers to understand when they use the standard as guidance in practice. A method using reliability block diagram is investigated in this study in order to provide a clear and feasible way of PFD avg calculation and help those who take IEC 61508-6 as their guidance. The method finds mean down times (MDTs) of both channel and voted group first and then PFD avg . The calculated results of various voted groups are compared with those in IEC61508 part 6 and Ref. [Zhang T, Long W, Sato Y. Availability of systems with self-diagnostic components-applying Markov model to IEC 61508-6. Reliab Eng System Saf 2003;80(2):133-41]. An interesting outcome can be realized from the comparison. Furthermore, although differences in MDT of voted groups exist between IEC 61508-6 and this paper, PFD avg of voted groups are comparatively close. With detailed description, the method of RBD presented can be applied to the quantitative SIL verification, showing a similarity of the method in IEC 61508-6

Ion-beam apparatus and method for analyzing and controlling integrated circuits

Science.gov (United States)

Campbell, Ann N.; Soden, Jerry M.

1998-01-01

An ion-beam apparatus and method for analyzing and controlling integrated circuits. The ion-beam apparatus comprises a stage for holding one or more integrated circuits (ICs); a source means for producing a focused ion beam; and a beam-directing means for directing the focused ion beam to irradiate a predetermined portion of the IC for sufficient time to provide an ion-beam-generated electrical input signal to a predetermined element of the IC. The apparatus and method have applications to failure analysis and developmental analysis of ICs and permit an alteration, control, or programming of logic states or device parameters within the IC either separate from or in combination with applied electrical stimulus to the IC for analysis thereof. Preferred embodiments of the present invention including a secondary particle detector and an electron floodgun further permit imaging of the IC by secondary ions or electrons, and allow at least a partial removal or erasure of the ion-beam-generated electrical input signal.

Improved Genetic Algorithm-Based Unit Commitment Considering Uncertainty Integration Method

Directory of Open Access Journals (Sweden)

Kyu-Hyung Jo

2018-05-01

Full Text Available In light of the dissemination of renewable energy connected to the power grid, it has become necessary to consider the uncertainty in the generation of renewable energy as a unit commitment (UC problem. A methodology for solving the UC problem is presented by considering various uncertainties, which are assumed to have a normal distribution, by using a Monte Carlo simulation. Based on the constructed scenarios for load, wind, solar, and generator outages, a combination of scenarios is found that meets the reserve requirement to secure the power balance of the power grid. In those scenarios, the uncertainty integration method (UIM identifies the best combination by minimizing the additional reserve requirements caused by the uncertainty of power sources. An integration process for uncertainties is formulated for stochastic unit commitment (SUC problems and optimized by the improved genetic algorithm (IGA. The IGA is composed of five procedures and finds the optimal combination of unit status at the scheduled time, based on the determined source data. According to the number of unit systems, the IGA demonstrates better performance than the other optimization methods by applying reserve repairing and an approximation process. To account for the result of the proposed method, various UC strategies are tested with a modified 24-h UC test system and compared.

A highly accurate boundary integral equation method for surfactant-laden drops in 3D

Science.gov (United States)

Sorgentone, Chiara; Tornberg, Anna-Karin

2018-05-01

The presence of surfactants alters the dynamics of viscous drops immersed in an ambient viscous fluid. This is specifically true at small scales, such as in applications of droplet based microfluidics, where the interface dynamics become of increased importance. At such small scales, viscous forces dominate and inertial effects are often negligible. Considering Stokes flow, a numerical method based on a boundary integral formulation is presented for simulating 3D drops covered by an insoluble surfactant. The method is able to simulate drops with different viscosities and close interactions, automatically controlling the time step size and maintaining high accuracy also when substantial drop deformation appears. To achieve this, the drop surfaces as well as the surfactant concentration on each surface are represented by spherical harmonics expansions. A novel reparameterization method is introduced to ensure a high-quality representation of the drops also under deformation, specialized quadrature methods for singular and nearly singular integrals that appear in the formulation are evoked and the adaptive time stepping scheme for the coupled drop and surfactant evolution is designed with a preconditioned implicit treatment of the surfactant diffusion.

A new integrated evaluation method of heavy metals pollution control during melting and sintering of MSWI fly ash.

Science.gov (United States)

Li, Rundong; Li, Yanlong; Yang, Tianhua; Wang, Lei; Wang, Weiyun

2015-05-30

Evaluations of technologies for heavy metal control mainly examine the residual and leaching rates of a single heavy metal, such that developed evaluation method have no coordination or uniqueness and are therefore unsuitable for hazard control effect evaluation. An overall pollution toxicity index (OPTI) was established in this paper, based on the developed index, an integrated evaluation method of heavy metal pollution control was established. Application of this method in the melting and sintering of fly ash revealed the following results: The integrated control efficiency of the melting process was higher in all instances than that of the sintering process. The lowest integrated control efficiency of melting was 56.2%, and the highest integrated control efficiency of sintering was 46.6%. Using the same technology, higher integrated control efficiency conditions were all achieved with lower temperatures and shorter times. This study demonstrated the unification and consistency of this method. Copyright © 2015 Elsevier B.V. All rights reserved.

Monolitic integrated circuit for the strobed charge-to-time converter

International Nuclear Information System (INIS)

Bel'skij, V.I.; Bushnin, Yu.B.; Zimin, S.A.; Punzhin, Yu.N.; Sen'ko, V.A.; Soldatov, M.M.; Tokarchuk, V.P.

1985-01-01

The developed and comercially produced semiconducting circuit - gating charge-to-time converter KR1101PD1 is described. The considered integrated circuit is a short pulse charge-to-time converter with integration of input current. The circuit is designed for construction of time-to-pulse analog-to-digital converters utilized in multichannel detection systems when studying complex topology processes. Input resistance of the circuit is 0.1 Ω permissible input current is 50 mA, maximum measured charge is 300-1000 pC

Improved parallel solution techniques for the integral transport matrix method

Energy Technology Data Exchange (ETDEWEB)

Zerr, R. Joseph, E-mail: rjz116@psu.edu [Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, University Park, PA (United States); Azmy, Yousry Y., E-mail: yyazmy@ncsu.edu [Department of Nuclear Engineering, North Carolina State University, Burlington Engineering Laboratories, Raleigh, NC (United States)

2011-07-01

Alternative solution strategies to the parallel block Jacobi (PBJ) method for the solution of the global problem with the integral transport matrix method operators have been designed and tested. The most straightforward improvement to the Jacobi iterative method is the Gauss-Seidel alternative. The parallel red-black Gauss-Seidel (PGS) algorithm can improve on the number of iterations and reduce work per iteration by applying an alternating red-black color-set to the subdomains and assigning multiple sub-domains per processor. A parallel GMRES(m) method was implemented as an alternative to stationary iterations. Computational results show that the PGS method can improve on the PBJ method execution time by up to 10´ when eight sub-domains per processor are used. However, compared to traditional source iterations with diffusion synthetic acceleration, it is still approximately an order of magnitude slower. The best-performing cases are optically thick because sub-domains decouple, yielding faster convergence. Further tests revealed that 64 sub-domains per processor was the best performing level of sub-domain division. An acceleration technique that improves the convergence rate would greatly improve the ITMM. The GMRES(m) method with a diagonal block pre conditioner consumes approximately the same time as the PBJ solver but could be improved by an as yet undeveloped, more efficient pre conditioner. (author)

Improved parallel solution techniques for the integral transport matrix method

International Nuclear Information System (INIS)

Zerr, R. Joseph; Azmy, Yousry Y.

2011-01-01

Alternative solution strategies to the parallel block Jacobi (PBJ) method for the solution of the global problem with the integral transport matrix method operators have been designed and tested. The most straightforward improvement to the Jacobi iterative method is the Gauss-Seidel alternative. The parallel red-black Gauss-Seidel (PGS) algorithm can improve on the number of iterations and reduce work per iteration by applying an alternating red-black color-set to the subdomains and assigning multiple sub-domains per processor. A parallel GMRES(m) method was implemented as an alternative to stationary iterations. Computational results show that the PGS method can improve on the PBJ method execution time by up to 10´ when eight sub-domains per processor are used. However, compared to traditional source iterations with diffusion synthetic acceleration, it is still approximately an order of magnitude slower. The best-performing cases are optically thick because sub-domains decouple, yielding faster convergence. Further tests revealed that 64 sub-domains per processor was the best performing level of sub-domain division. An acceleration technique that improves the convergence rate would greatly improve the ITMM. The GMRES(m) method with a diagonal block pre conditioner consumes approximately the same time as the PBJ solver but could be improved by an as yet undeveloped, more efficient pre conditioner. (author)

A Fully Integrated Discrete-Time Superheterodyne Receiver

NARCIS (Netherlands)

Tohidian, M.; Madadi, I.; Staszewski, R.B.

2017-01-01

The zero/low intermediate frequency (IF) receiver (RX) architecture has enabled full CMOS integration. As the technology scales and wireless standards become ever more challenging, the issues related to time-varying dc offsets, the second-order nonlinearity, and flicker noise become more critical.

A dynamic integrated fault diagnosis method for power transformers.

Science.gov (United States)

Gao, Wensheng; Bai, Cuifen; Liu, Tong

2015-01-01

In order to diagnose transformer fault efficiently and accurately, a dynamic integrated fault diagnosis method based on Bayesian network is proposed in this paper. First, an integrated fault diagnosis model is established based on the causal relationship among abnormal working conditions, failure modes, and failure symptoms of transformers, aimed at obtaining the most possible failure mode. And then considering the evidence input into the diagnosis model is gradually acquired and the fault diagnosis process in reality is multistep, a dynamic fault diagnosis mechanism is proposed based on the integrated fault diagnosis model. Different from the existing one-step diagnosis mechanism, it includes a multistep evidence-selection process, which gives the most effective diagnostic test to be performed in next step. Therefore, it can reduce unnecessary diagnostic tests and improve the accuracy and efficiency of diagnosis. Finally, the dynamic integrated fault diagnosis method is applied to actual cases, and the validity of this method is verified.

A Dynamic Integrated Fault Diagnosis Method for Power Transformers

Science.gov (United States)

Gao, Wensheng; Liu, Tong

2015-01-01

In order to diagnose transformer fault efficiently and accurately, a dynamic integrated fault diagnosis method based on Bayesian network is proposed in this paper. First, an integrated fault diagnosis model is established based on the causal relationship among abnormal working conditions, failure modes, and failure symptoms of transformers, aimed at obtaining the most possible failure mode. And then considering the evidence input into the diagnosis model is gradually acquired and the fault diagnosis process in reality is multistep, a dynamic fault diagnosis mechanism is proposed based on the integrated fault diagnosis model. Different from the existing one-step diagnosis mechanism, it includes a multistep evidence-selection process, which gives the most effective diagnostic test to be performed in next step. Therefore, it can reduce unnecessary diagnostic tests and improve the accuracy and efficiency of diagnosis. Finally, the dynamic integrated fault diagnosis method is applied to actual cases, and the validity of this method is verified. PMID:25685841

Optimal order and time-step criterion for Aarseth-type N-body integrators

International Nuclear Information System (INIS)

Makino, Junichiro

1991-01-01

How the selection of the time-step criterion and the order of the integrator change the efficiency of Aarseth-type N-body integrators is discussed. An alternative to Aarseth's scheme based on the direct calculation of the time derivative of the force using the Hermite interpolation is compared to Aarseth's scheme, which uses the Newton interpolation to construct the predictor and corrector. How the number of particles in the system changes the behavior of integrators is examined. The Hermite scheme allows a time step twice as large as that for the standard Aarseth scheme for the same accuracy. The calculation cost of the Hermite scheme per time step is roughly twice as much as that of the standard Aarseth scheme. The optimal order of the integrators depends on both the particle number and the accuracy required. The time-step criterion of the standard Aarseth scheme is found to be inapplicable to higher-order integrators, and a more uniformly reliable criterion is proposed. 18 refs

Time-optimal control of nuclear reactor power with adaptive proportional- integral-feedforward gains

International Nuclear Information System (INIS)

Park, Moon Ghu; Cho, Nam Zin

1993-01-01

A time-optimal control method which consists of coarse and fine control stages is described here. During the coarse control stage, the maximum control effort (time-optimal) is used to direct the system toward the switching boundary which is set near the desired power level. At this boundary, the controller is switched to the fine control stage in which an adaptive proportional-integral-feedforward (PIF) controller is used to compensate for any unmodeled reactivity feedback effects. This fine control is also introduced to obtain a constructive method for determining the (adaptive) feedback gains against the sampling effect. The feedforward control term is included to suppress the over-or undershoot. The estimation and feedback of the temperature-induced reactivity is also discussed

14th International Conference on Integral Methods in Science and Engineering

CERN Document Server

Riva, Matteo; Lamberti, Pier; Musolino, Paolo

2017-01-01

This contributed volume contains a collection of articles on the most recent advances in integral methods.  The first of two volumes, this work focuses on the construction of theoretical integral methods. Written by internationally recognized researchers, the chapters in this book are based on talks given at the Fourteenth International Conference on Integral Methods in Science and Engineering, held July 25-29, 2016, in Padova, Italy. A broad range of topics is addressed, such as: • Integral equations • Homogenization • Duality methods • Optimal design • Conformal techniques This collection will be of interest to researchers in applied mathematics, physics, and mechanical and electrical engineering, as well as graduate students in these disciplines, and to other professionals who use integration as an essential tool in their work.

Cultural adaptation and translation of measures: an integrated method.

Science.gov (United States)

Sidani, Souraya; Guruge, Sepali; Miranda, Joyal; Ford-Gilboe, Marilyn; Varcoe, Colleen

2010-04-01

Differences in the conceptualization and operationalization of health-related concepts may exist across cultures. Such differences underscore the importance of examining conceptual equivalence when adapting and translating instruments. In this article, we describe an integrated method for exploring conceptual equivalence within the process of adapting and translating measures. The integrated method involves five phases including selection of instruments for cultural adaptation and translation; assessment of conceptual equivalence, leading to the generation of a set of items deemed to be culturally and linguistically appropriate to assess the concept of interest in the target community; forward translation; back translation (optional); and pre-testing of the set of items. Strengths and limitations of the proposed integrated method are discussed. (c) 2010 Wiley Periodicals, Inc.

Interconnect rise time in superconducting integrating circuits

International Nuclear Information System (INIS)

Preis, D.; Shlager, K.

1988-01-01

The influence of resistive losses on the voltage rise time of an integrated-circuit interconnection is reported. A distribution-circuit model is used to present the interconnect. Numerous parametric curves are presented based on numerical evaluation of the exact analytical expression for the model's transient response. For the superconducting case in which the series resistance of the interconnect approaches zero, the step-response rise time is longer but signal strength increases significantly

ARE METHODS USED TO INTEGRATE STANDARDIZED MANAGEMENT SYSTEMS A CONDITIONING FACTOR OF THE LEVEL OF INTEGRATION? AN EMPIRICAL STUDY

Directory of Open Access Journals (Sweden)

Merce Bernardo

2011-09-01

Full Text Available Organizations are increasingly implementing multiple Management System Standards (M SSs and considering managing the related Management Systems (MSs as a single system.The aim of this paper is to analyze if methods us ed to integrate standardized MSs condition the level of integration of those MSs. A descriptive methodology has been applied to 343 Spanish organizations registered to, at least, ISO 9001 and ISO 14001. Seven groups of these organizations using different combinations of methods have been analyzed Results show that these organizations have a high level of integration of their MSs. The most common method used, was the process map. Organizations using a combination of different methods achieve higher levels of integration than those using a single method. However, no evidence has been found to confirm the relationship between the method used and the integration level achieved.

INTEGRATED FUSION METHOD FOR MULTIPLE TEMPORAL-SPATIAL-SPECTRAL IMAGES

Directory of Open Access Journals (Sweden)

H. Shen

2012-08-01

Full Text Available Data fusion techniques have been widely researched and applied in remote sensing field. In this paper, an integrated fusion method for remotely sensed images is presented. Differently from the existed methods, the proposed method has the performance to integrate the complementary information in multiple temporal-spatial-spectral images. In order to represent and process the images in one unified framework, two general image observation models are firstly presented, and then the maximum a posteriori (MAP framework is used to set up the fusion model. The gradient descent method is employed to solve the fused image. The efficacy of the proposed method is validated using simulated images.

Aspects for Run-time Component Integration

DEFF Research Database (Denmark)

Truyen, Eddy; Jørgensen, Bo Nørregaard; Joosen, Wouter

2000-01-01

Component framework technology has become the cornerstone of building a family of systems and applications. A component framework defines a generic architecture into which specialized components can be plugged. As such, the component framework leverages the glue that connects the different inserted...... to dynamically integrate into the architecture of middleware systems new services that support non-functional aspects such as security, transactions, real-time....

Deep Time Data Infrastructure: Integrating Our Current Geologic and Biologic Databases

Science.gov (United States)

Kolankowski, S. M.; Fox, P. A.; Ma, X.; Prabhu, A.

2016-12-01

As our knowledge of Earth's geologic and mineralogical history grows, we require more efficient methods of sharing immense amounts of data. Databases across numerous disciplines have been utilized to offer extensive information on very specific Epochs of Earth's history up to its current state, i.e. Fossil record, rock composition, proteins, etc. These databases could be a powerful force in identifying previously unseen correlations such as relationships between minerals and proteins. Creating a unifying site that provides a portal to these databases will aid in our ability as a collaborative scientific community to utilize our findings more effectively. The Deep-Time Data Infrastructure (DTDI) is currently being defined as part of a larger effort to accomplish this goal. DTDI will not be a new database, but an integration of existing resources. Current geologic and related databases were identified, documentation of their schema was established and will be presented as a stage by stage progression. Through conceptual modeling focused around variables from their combined records, we will determine the best way to integrate these databases using common factors. The Deep-Time Data Infrastructure will allow geoscientists to bridge gaps in data and further our understanding of our Earth's history.

Application of heat-balance integral method to conjugate thermal explosion

Directory of Open Access Journals (Sweden)

Novozhilov Vasily

2009-01-01

Full Text Available Conjugate thermal explosion is an extension of the classical theory, proposed and studied recently by the author. The paper reports application of heat-balance integral method for developing phase portraits for systems undergoing conjugate thermal explosion. The heat-balance integral method is used as an averaging method reducing partical differential equation problem to the set of first-order ordinary differential equations. The latter reduced problem allows natural interpretation in appropriately chosen phase space. It is shown that, with the help of heat-balance integral technique, conjugate thermal explosion problem can be described with a good accuracy by the set of non-linear first-order differential equations involving complex error function. Phase trajectories are presented for typical regimes emerging in conjugate thermal explosion. Use of heat-balance integral as a spatial averaging method allows efficient description of system evolution to be developed.

Uniform stable conformal convolutional perfectly matched layer for enlarged cell technique conformal finite-difference time-domain method

International Nuclear Information System (INIS)

Wang Yue; Wang Jian-Guo; Chen Zai-Gao

2015-01-01

Based on conformal construction of physical model in a three-dimensional Cartesian grid, an integral-based conformal convolutional perfectly matched layer (CPML) is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain (ECT-CFDTD) method is used to simulate the wave propagation inside a perfect electric conductor (PEC) waveguide. The algorithm has the same numerical stability as the ECT-CFDTD method. For the long-time propagation problems of an evanescent wave in a waveguide, several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML. Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide. (paper)

An integrated lean-methods approach to hospital facilities redesign.

Science.gov (United States)

Nicholas, John

2012-01-01

Lean production methods for eliminating waste and improving processes in manufacturing are now being applied in healthcare. As the author shows, the methods are appropriate for redesigning hospital facilities. When used in an integrated manner and employing teams of mostly clinicians, the methods produce facility designs that are custom-fit to patient needs and caregiver work processes, and reduce operational costs. The author reviews lean methods and an approach for integrating them in the redesign of hospital facilities. A case example of the redesign of an emergency department shows the feasibility and benefits of the approach.

Method of manufacturing Josephson junction integrated circuits

International Nuclear Information System (INIS)

Jillie, D.W. Jr.; Smith, L.N.

1985-01-01

Josephson junction integrated circuits of the current injection type and magnetically controlled type utilize a superconductive layer that forms both Josephson junction electrode for the Josephson junction devices on the integrated circuit as well as a ground plane for the integrated circuit. Large area Josephson junctions are utilized for effecting contact to lower superconductive layers and islands are formed in superconductive layers to provide isolation between the groudplane function and the Josephson junction electrode function as well as to effect crossovers. A superconductor-barrier-superconductor trilayer patterned by local anodization is also utilized with additional layers formed thereover. Methods of manufacturing the embodiments of the invention are disclosed

Criteria for quantitative and qualitative data integration: mixed-methods research methodology.

Science.gov (United States)

Lee, Seonah; Smith, Carrol A M

2012-05-01

Many studies have emphasized the need and importance of a mixed-methods approach for evaluation of clinical information systems. However, those studies had no criteria to guide integration of multiple data sets. Integrating different data sets serves to actualize the paradigm that a mixed-methods approach argues; thus, we require criteria that provide the right direction to integrate quantitative and qualitative data. The first author used a set of criteria organized from a literature search for integration of multiple data sets from mixed-methods research. The purpose of this article was to reorganize the identified criteria. Through critical appraisal of the reasons for designing mixed-methods research, three criteria resulted: validation, complementarity, and discrepancy. In applying the criteria to empirical data of a previous mixed methods study, integration of quantitative and qualitative data was achieved in a systematic manner. It helped us obtain a better organized understanding of the results. The criteria of this article offer the potential to produce insightful analyses of mixed-methods evaluations of health information systems.

High-order Path Integral Monte Carlo methods for solving strongly correlated fermion problems

Science.gov (United States)

Chin, Siu A.

2015-03-01

In solving for the ground state of a strongly correlated many-fermion system, the conventional second-order Path Integral Monte Carlo method is plagued with the sign problem. This is due to the large number of anti-symmetric free fermion propagators that are needed to extract the square of the ground state wave function at large imaginary time. In this work, I show that optimized fourth-order Path Integral Monte Carlo methods, which uses no more than 5 free-fermion propagators, in conjunction with the use of the Hamiltonian energy estimator, can yield accurate ground state energies for quantum dots with up to 20 polarized electrons. The correlations are directly built-in and no explicit wave functions are needed. This work is supported by the Qatar National Research Fund NPRP GRANT #5-674-1-114.

Integrative methods for analyzing big data in precision medicine.

Science.gov (United States)

Gligorijević, Vladimir; Malod-Dognin, Noël; Pržulj, Nataša

2016-03-01

We provide an overview of recent developments in big data analyses in the context of precision medicine and health informatics. With the advance in technologies capturing molecular and medical data, we entered the area of "Big Data" in biology and medicine. These data offer many opportunities to advance precision medicine. We outline key challenges in precision medicine and present recent advances in data integration-based methods to uncover personalized information from big data produced by various omics studies. We survey recent integrative methods for disease subtyping, biomarkers discovery, and drug repurposing, and list the tools that are available to domain scientists. Given the ever-growing nature of these big data, we highlight key issues that big data integration methods will face. © 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

New readout integrated circuit using continuous time fixed pattern noise correction

Science.gov (United States)

Dupont, Bertrand; Chammings, G.; Rapellin, G.; Mandier, C.; Tchagaspanian, M.; Dupont, Benoit; Peizerat, A.; Yon, J. J.

2008-04-01

LETI has been involved in IRFPA development since 1978; the design department (LETI/DCIS) has focused its work on new ROIC architecture since many years. The trend is to integrate advanced functions into the CMOS design to achieve cost efficient sensors production. Thermal imaging market is today more and more demanding of systems with instant ON capability and low power consumption. The purpose of this paper is to present the latest developments of fixed pattern noise continuous time correction. Several architectures are proposed, some are based on hardwired digital processing and some are purely analog. Both are using scene based algorithms. Moreover a new method is proposed for simultaneous correction of pixel offsets and sensitivities. In this scope, a new architecture of readout integrated circuit has been implemented; this architecture is developed with 0.18μm CMOS technology. The specification and the application of the ROIC are discussed in details.

A Multi-Objective Optimization Method to integrate Heat Pumps in Industrial Processes

OpenAIRE

Becker, Helen; Spinato, Giulia; Maréchal, François

2011-01-01

Aim of process integration methods is to increase the efficiency of industrial processes by using pinch analysis combined with process design methods. In this context, appropriate integrated utilities offer promising opportunities to reduce energy consumption, operating costs and pollutants emissions. Energy integration methods are able to integrate any type of predefined utility, but so far there is no systematic approach to generate potential utilities models based on their technology limit...

Hierarchical Matrices Method and Its Application in Electromagnetic Integral Equations

Directory of Open Access Journals (Sweden)

Han Guo

2012-01-01

Full Text Available Hierarchical (H- matrices method is a general mathematical framework providing a highly compact representation and efficient numerical arithmetic. When applied in integral-equation- (IE- based computational electromagnetics, H-matrices can be regarded as a fast algorithm; therefore, both the CPU time and memory requirement are reduced significantly. Its kernel independent feature also makes it suitable for any kind of integral equation. To solve H-matrices system, Krylov iteration methods can be employed with appropriate preconditioners, and direct solvers based on the hierarchical structure of H-matrices are also available along with high efficiency and accuracy, which is a unique advantage compared to other fast algorithms. In this paper, a novel sparse approximate inverse (SAI preconditioner in multilevel fashion is proposed to accelerate the convergence rate of Krylov iterations for solving H-matrices system in electromagnetic applications, and a group of parallel fast direct solvers are developed for dealing with multiple right-hand-side cases. Finally, numerical experiments are given to demonstrate the advantages of the proposed multilevel preconditioner compared to conventional “single level” preconditioners and the practicability of the fast direct solvers for arbitrary complex structures.

A generic, time-resolved, integrated digital image correlation, identification approach

NARCIS (Netherlands)

Hoefnagels, J.P.M.; Neggers, J.; Blaysat, Benoît; Hild, François; Geers, M.G.D.; Jin, H.; Sciammarella, C.; Yoshida, S.; Lamberti, L.

2015-01-01

A generic one-step Integrated Digital Image Correlation (I-DIC) inverse parameter identification approach is introduced that enables direct identification of constitutive model parameters by intimately integrating a Finite Elements Method (FEM) with Digital Image Correlation (DIC), directly

Simulation of three-dimensional, time-dependent, incompressible flows by a finite element method

International Nuclear Information System (INIS)

Chan, S.T.; Gresho, P.M.; Lee, R.L.; Upson, C.D.

1981-01-01

A finite element model has been developed for simulating the dynamics of problems encountered in atmospheric pollution and safety assessment studies. The model is based on solving the set of three-dimensional, time-dependent, conservation equations governing incompressible flows. Spatial discretization is performed via a modified Galerkin finite element method, and time integration is carried out via the forward Euler method (pressure is computed implicitly, however). Several cost-effective techniques (including subcycling, mass lumping, and reduced Gauss-Legendre quadrature) which have been implemented are discussed. Numerical results are presented to demonstrate the applicability of the model

Tau method approximation of the Hubbell rectangular source integral

International Nuclear Information System (INIS)

Kalla, S.L.; Khajah, H.G.

2000-01-01

The Tau method is applied to obtain expansions, in terms of Chebyshev polynomials, which approximate the Hubbell rectangular source integral:I(a,b)=∫ b 0 (1/(√(1+x 2 )) arctan(a/(√(1+x 2 )))) This integral corresponds to the response of an omni-directional radiation detector situated over a corner of a plane isotropic rectangular source. A discussion of the error in the Tau method approximation follows

Real-time GPS seismology using a single receiver: method comparison, error analysis and precision validation

Science.gov (United States)

Li, Xingxing

2014-05-01

Earthquake monitoring and early warning system for hazard assessment and mitigation has traditional been based on seismic instruments. However, for large seismic events, it is difficult for traditional seismic instruments to produce accurate and reliable displacements because of the saturation of broadband seismometers and problematic integration of strong-motion data. Compared with the traditional seismic instruments, GPS can measure arbitrarily large dynamic displacements without saturation, making them particularly valuable in case of large earthquakes and tsunamis. GPS relative positioning approach is usually adopted to estimate seismic displacements since centimeter-level accuracy can be achieved in real-time by processing double-differenced carrier-phase observables. However, relative positioning method requires a local reference station, which might itself be displaced during a large seismic event, resulting in misleading GPS analysis results. Meanwhile, the relative/network approach is time-consuming, particularly difficult for the simultaneous and real-time analysis of GPS data from hundreds or thousands of ground stations. In recent years, several single-receiver approaches for real-time GPS seismology, which can overcome the reference station problem of the relative positioning approach, have been successfully developed and applied to GPS seismology. One available method is real-time precise point positioning (PPP) relied on precise satellite orbit and clock products. However, real-time PPP needs a long (re)convergence period, of about thirty minutes, to resolve integer phase ambiguities and achieve centimeter-level accuracy. In comparison with PPP, Colosimo et al. (2011) proposed a variometric approach to determine the change of position between two adjacent epochs, and then displacements are obtained by a single integration of the delta positions. This approach does not suffer from convergence process, but the single integration from delta positions to

Generalized Runge-Kutta method for two- and three-dimensional space-time diffusion equations with a variable time step

International Nuclear Information System (INIS)

Aboanber, A.E.; Hamada, Y.M.

2008-01-01

An extensive knowledge of the spatial power distribution is required for the design and analysis of different types of current-generation reactors, and that requires the development of more sophisticated theoretical methods. Therefore, the need to develop new methods for multidimensional transient reactor analysis still exists. The objective of this paper is to develop a computationally efficient numerical method for solving the multigroup, multidimensional, static and transient neutron diffusion kinetics equations. A generalized Runge-Kutta method has been developed for the numerical integration of the stiff space-time diffusion equations. The method is fourth-order accurate, using an embedded third-order solution to arrive at an estimate of the truncation error for automatic time step control. In addition, the A(α)-stability properties of the method are investigated. The analyses of two- and three-dimensional benchmark problems as well as static and transient problems, demonstrate that very accurate solutions can be obtained with assembly-sized spatial meshes. Preliminary numerical evaluations using two- and three-dimensional finite difference codes showed that the presented generalized Runge-Kutta method is highly accurate and efficient when compared with other optimized iterative numerical and conventional finite difference methods

On the relationship between supplier integration and time-to-market

NARCIS (Netherlands)

Perols, J.; Zimmermann, C.; Kortmann, S.

2013-01-01

Recent operations management and innovation management research emphasizes the importance of supplier integration. However, the empirical results as to the relationship between supplier integration and time-to-market are ambivalent. To understand this important relationship, we incorporate two major

Integrals of Motion for Discrete-Time Optimal Control Problems

OpenAIRE

Torres, Delfim F. M.

2003-01-01

We obtain a discrete time analog of E. Noether's theorem in Optimal Control, asserting that integrals of motion associated to the discrete time Pontryagin Maximum Principle can be computed from the quasi-invariance properties of the discrete time Lagrangian and discrete time control system. As corollaries, results for first-order and higher-order discrete problems of the calculus of variations are obtained.

Variational method for integrating radial gradient field

Science.gov (United States)

Legarda-Saenz, Ricardo; Brito-Loeza, Carlos; Rivera, Mariano; Espinosa-Romero, Arturo

2014-12-01

We propose a variational method for integrating information obtained from circular fringe pattern. The proposed method is a suitable choice for objects with radial symmetry. First, we analyze the information contained in the fringe pattern captured by the experimental setup and then move to formulate the problem of recovering the wavefront using techniques from calculus of variations. The performance of the method is demonstrated by numerical experiments with both synthetic and real data.

First integral method for an oscillator system

Directory of Open Access Journals (Sweden)

Xiaoqian Gong

2013-04-01

Full Text Available In this article, we consider the nonlinear Duffing-van der Pol-type oscillator system by means of the first integral method. This system has physical relevance as a model in certain flow-induced structural vibration problems, which includes the van der Pol oscillator and the damped Duffing oscillator etc as particular cases. Firstly, we apply the Division Theorem for two variables in the complex domain, which is based on the ring theory of commutative algebra, to explore a quasi-polynomial first integral to an equivalent autonomous system. Then, through solving an algebraic system we derive the first integral of the Duffing-van der Pol-type oscillator system under certain parametric condition.

Time integration in the code Zgoubi and external usage of PTC's structures

International Nuclear Information System (INIS)

Forest, Etienne; Meot, F.

2006-06-01

The purpose of this note is to describe Zgoubi's integrator and to describe some pitfalls for time based integration when used in accelerators. We show why the convergence rate of an integrator can be affected by an improper treatment at the boundary when time is used as the integration variable. We also point out how the code PTC can be used as a container by other tracking engine. This work is not completed as far as incorporation of Zgoubi is concerned. (authors)

New method for calculation of integral characteristics of thermal plumes

DEFF Research Database (Denmark)

Zukowska, Daria; Popiolek, Zbigniew; Melikov, Arsen Krikor

2008-01-01

A method for calculation of integral characteristics of thermal plumes is proposed. The method allows for determination of the integral parameters of plumes based on speed measurements performed with omnidirectional low velocity thermoanemometers. The method includes a procedure for calculation...... of the directional velocity (upward component of the mean velocity). The method is applied for determination of the characteristics of an asymmetric thermal plume generated by a sitting person. The method was validated in full-scale experiments in a climatic chamber with a thermal manikin as a simulator of a sitting...

Parallel, explicit, and PWTD-enhanced time domain volume integral equation solver

KAUST Repository

Liu, Yang

2013-07-01

Time domain volume integral equations (TDVIEs) are useful for analyzing transient scattering from inhomogeneous dielectric objects in applications as varied as photonics, optoelectronics, and bioelectromagnetics. TDVIEs typically are solved by implicit marching-on-in-time (MOT) schemes [N. T. Gres et al., Radio Sci., 36, 379-386, 2001], requiring the solution of a system of equations at each and every time step. To reduce the computational cost associated with such schemes, [A. Al-Jarro et al., IEEE Trans. Antennas Propagat., 60, 5203-5215, 2012] introduced an explicit MOT-TDVIE method that uses a predictor-corrector technique to stably update field values throughout the scatterer. By leveraging memory-efficient nodal spatial discretization and scalable parallelization schemes [A. Al-Jarro et al., in 28th Int. Rev. Progress Appl. Computat. Electromagn., 2012], this solver has been successfully applied to the analysis of scattering phenomena involving 0.5 million spatial unknowns. © 2013 IEEE.

On the internal resonant modes in marching-on-in-time solution of the time domain electric field integral equation

KAUST Repository

Shi, Yifei; Bagci, Hakan; Lu, Mingyu

2013-01-01

Internal resonant modes are always observed in the marching-on-in-time (MOT) solution of the time domain electric field integral equation (EFIE), although 'relaxed initial conditions,' which are enforced at the beginning of time marching, should in theory prevent these spurious modes from appearing. It has been conjectured that, numerical errors built up during time marching establish the necessary initial conditions and induce the internal resonant modes. However, this conjecture has never been proved by systematic numerical experiments. Our numerical results in this communication demonstrate that, the internal resonant modes' amplitudes are indeed dictated by the numerical errors. Additionally, it is shown that in a few cases, the internal resonant modes can be made 'invisible' by significantly suppressing the numerical errors. These tests prove the conjecture that the internal resonant modes are induced by numerical errors when the time domain EFIE is solved by the MOT method. © 2013 IEEE.

On the internal resonant modes in marching-on-in-time solution of the time domain electric field integral equation

KAUST Repository

Shi, Yifei

2013-08-01

Internal resonant modes are always observed in the marching-on-in-time (MOT) solution of the time domain electric field integral equation (EFIE), although \\'relaxed initial conditions,\\' which are enforced at the beginning of time marching, should in theory prevent these spurious modes from appearing. It has been conjectured that, numerical errors built up during time marching establish the necessary initial conditions and induce the internal resonant modes. However, this conjecture has never been proved by systematic numerical experiments. Our numerical results in this communication demonstrate that, the internal resonant modes\\' amplitudes are indeed dictated by the numerical errors. Additionally, it is shown that in a few cases, the internal resonant modes can be made \\'invisible\\' by significantly suppressing the numerical errors. These tests prove the conjecture that the internal resonant modes are induced by numerical errors when the time domain EFIE is solved by the MOT method. © 2013 IEEE.

A Novel Multiple-Time Scale Integrator for the Hybrid Monte Carlo Algorithm

International Nuclear Information System (INIS)

Kamleh, Waseem

2011-01-01

Hybrid Monte Carlo simulations that implement the fermion action using multiple terms are commonly used. By the nature of their formulation they involve multiple integration time scales in the evolution of the system through simulation time. These different scales are usually dealt with by the Sexton-Weingarten nested leapfrog integrator. In this scheme the choice of time scales is somewhat restricted as each time step must be an exact multiple of the next smallest scale in the sequence. A novel generalisation of the nested leapfrog integrator is introduced which allows for far greater flexibility in the choice of time scales, as each scale now must only be an exact multiple of the smallest step size.

Integral-equation based methods for parameter estimation in output pulses of radiation detectors: Application in nuclear medicine and spectroscopy

Science.gov (United States)

Mohammadian-Behbahani, Mohammad-Reza; Saramad, Shahyar

2018-04-01

Model based analysis methods are relatively new approaches for processing the output data of radiation detectors in nuclear medicine imaging and spectroscopy. A class of such methods requires fast algorithms for fitting pulse models to experimental data. In order to apply integral-equation based methods for processing the preamplifier output pulses, this article proposes a fast and simple method for estimating the parameters of the well-known bi-exponential pulse model by solving an integral equation. The proposed method needs samples from only three points of the recorded pulse as well as its first and second order integrals. After optimizing the sampling points, the estimation results were calculated and compared with two traditional integration-based methods. Different noise levels (signal-to-noise ratios from 10 to 3000) were simulated for testing the functionality of the proposed method, then it was applied to a set of experimental pulses. Finally, the effect of quantization noise was assessed by studying different sampling rates. Promising results by the proposed method endorse it for future real-time applications.

Numerical counting ratemeter with variable time constant and integrated circuits

International Nuclear Information System (INIS)

Kaiser, J.; Fuan, J.

1967-01-01

We present here the prototype of a numerical counting ratemeter which is a special version of variable time-constant frequency meter (1). The originality of this work lies in the fact that the change in the time constant is carried out automatically. Since the criterion for this change is the accuracy in the annunciated result, the integration time is varied as a function of the frequency. For the prototype described in this report, the time constant varies from 1 sec to 1 millisec. for frequencies in the range 10 Hz to 10 MHz. This prototype is built entirely of MECL-type integrated circuits from Motorola and is thus contained in two relatively small boxes. (authors) [fr

Rapid expansion method (REM) for time‐stepping in reverse time migration (RTM)

KAUST Repository

Pestana, Reynam C.

2009-01-01

We show that the wave equation solution using a conventional finite‐difference scheme, derived commonly by the Taylor series approach, can be derived directly from the rapid expansion method (REM). After some mathematical manipulation we consider an analytical approximation for the Bessel function where we assume that the time step is sufficiently small. From this derivation we find that if we consider only the first two Chebyshev polynomials terms in the rapid expansion method we can obtain the second order time finite‐difference scheme that is frequently used in more conventional finite‐difference implementations. We then show that if we use more terms from the REM we can obtain a more accurate time integration of the wave field. Consequently, we have demonstrated that the REM is more accurate than the usual finite‐difference schemes and it provides a wave equation solution which allows us to march in large time steps without numerical dispersion and is numerically stable. We illustrate the method with post and pre stack migration results.

Explicit Time Integrators for Nonlinear Dynamics Derived from the Midpoint Rule

Directory of Open Access Journals (Sweden)

P. Krysl

2004-01-01

Full Text Available We address the design of time integrators for mechanical systems that are explicit in the forcing evaluations. Our starting point is the midpoint rule, either in the classical form for the vector space setting, or in the Lie form for the rotation group. By introducing discrete, concentrated impulses we can approximate the forcing impressed upon the system over the time step, and thus arrive at first-order integrators. These can then be composed to yield a second order integrator with very desirable properties: symplecticity and momentum conservation. 

Massively Parallel and Scalable Implicit Time Integration Algorithms for Structural Dynamics

Science.gov (United States)

Farhat, Charbel

1997-01-01

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

A Methodology for Conducting Integrative Mixed Methods Research and Data Analyses

Science.gov (United States)

Castro, Felipe González; Kellison, Joshua G.; Boyd, Stephen J.; Kopak, Albert

2011-01-01

Mixed methods research has gained visibility within the last few years, although limitations persist regarding the scientific caliber of certain mixed methods research designs and methods. The need exists for rigorous mixed methods designs that integrate various data analytic procedures for a seamless transfer of evidence across qualitative and quantitative modalities. Such designs can offer the strength of confirmatory results drawn from quantitative multivariate analyses, along with “deep structure” explanatory descriptions as drawn from qualitative analyses. This article presents evidence generated from over a decade of pilot research in developing an integrative mixed methods methodology. It presents a conceptual framework and methodological and data analytic procedures for conducting mixed methods research studies, and it also presents illustrative examples from the authors' ongoing integrative mixed methods research studies. PMID:22167325

Methods of assessing total doses integrated across pathways

International Nuclear Information System (INIS)

Grzechnik, M.; Camplin, W.; Clyne, F.; Allott, R.; Webbe-Wood, D.

2006-01-01

Calculated doses for comparison with limits resulting from discharges into the environment should be summed across all relevant pathways and food groups to ensure adequate protection. Current methodology for assessments used in the radioactivity in Food and the Environment (R.I.F.E.) reports separate doses from pathways related to liquid discharges of radioactivity to the environment from those due to gaseous releases. Surveys of local inhabitant food consumption and occupancy rates are conducted in the vicinity of nuclear sites. Information has been recorded in an integrated way, such that the data for each individual is recorded for all pathways of interest. These can include consumption of foods, such as fish, crustaceans, molluscs, fruit and vegetables, milk and meats. Occupancy times over beach sediments and time spent in close proximity to the site is also recorded for inclusion of external and inhalation radiation dose pathways. The integrated habits survey data may be combined with monitored environmental radionuclide concentrations to calculate total dose. The criteria for successful adoption of a method for this calculation were: Reproducibility can others easily use the approach and reassess doses? Rigour and realism how good is the match with reality?Transparency a measure of the ease with which others can understand how the calculations are performed and what they mean. Homogeneity is the group receiving the dose relatively homogeneous with respect to age, diet and those aspects that affect the dose received? Five methods of total dose calculation were compared and ranked according to their suitability. Each method was labelled (A to E) and given a short, relevant name for identification. The methods are described below; A) Individual doses to individuals are calculated and critical group selection is dependent on dose received. B) Individual Plus As in A, but consumption and occupancy rates for high dose is used to derive rates for application in

A travel time forecasting model based on change-point detection method

Science.gov (United States)

LI, Shupeng; GUANG, Xiaoping; QIAN, Yongsheng; ZENG, Junwei

2017-06-01

Travel time parameters obtained from road traffic sensors data play an important role in traffic management practice. A travel time forecasting model is proposed for urban road traffic sensors data based on the method of change-point detection in this paper. The first-order differential operation is used for preprocessing over the actual loop data; a change-point detection algorithm is designed to classify the sequence of large number of travel time data items into several patterns; then a travel time forecasting model is established based on autoregressive integrated moving average (ARIMA) model. By computer simulation, different control parameters are chosen for adaptive change point search for travel time series, which is divided into several sections of similar state.Then linear weight function is used to fit travel time sequence and to forecast travel time. The results show that the model has high accuracy in travel time forecasting.

A systematic and efficient method to compute multi-loop master integrals

Science.gov (United States)

Liu, Xiao; Ma, Yan-Qing; Wang, Chen-Yu

2018-04-01

We propose a novel method to compute multi-loop master integrals by constructing and numerically solving a system of ordinary differential equations, with almost trivial boundary conditions. Thus it can be systematically applied to problems with arbitrary kinematic configurations. Numerical tests show that our method can not only achieve results with high precision, but also be much faster than the only existing systematic method sector decomposition. As a by product, we find a new strategy to compute scalar one-loop integrals without reducing them to master integrals.

Computation of rectangular source integral by rational parameter polynomial method

International Nuclear Information System (INIS)

Prabha, Hem

2001-01-01

Hubbell et al. (J. Res. Nat Bureau Standards 64C, (1960) 121) have obtained a series expansion for the calculation of the radiation field generated by a plane isotropic rectangular source (plaque), in which leading term is the integral H(a,b). In this paper another integral I(a,b), which is related with the integral H(a,b) has been solved by the rational parameter polynomial method. From I(a,b), we compute H(a,b). Using this method the integral I(a,b) is expressed in the form of a polynomial of a rational parameter. Generally, a function f (x) is expressed in terms of x. In this method this is expressed in terms of x/(1+x). In this way, the accuracy of the expression is good over a wide range of x as compared to the earlier approach. The results for I(a,b) and H(a,b) are given for a sixth degree polynomial and are found to be in good agreement with the results obtained by numerically integrating the integral. Accuracy could be increased either by increasing the degree of the polynomial or by dividing the range of integration. The results of H(a,b) and I(a,b) are given for values of b and a up to 2.0 and 20.0, respectively

Integrated navigation method of a marine strapdown inertial navigation system using a star sensor

International Nuclear Information System (INIS)

Wang, Qiuying; Diao, Ming; Gao, Wei; Zhu, Minghong; Xiao, Shu

2015-01-01

This paper presents an integrated navigation method of the strapdown inertial navigation system (SINS) using a star sensor. According to the principle of SINS, its own navigation information contains an error that increases with time. Hence, the inertial attitude matrix from the star sensor is introduced as the reference information to correct the SINS increases error. For the integrated navigation method, the vehicle’s attitude can be obtained in two ways: one is calculated from SINS; the other, which we have called star sensor attitude, is obtained as the product between the SINS position and the inertial attitude matrix from the star sensor. Therefore, the SINS position error is introduced in the star sensor attitude error. Based on the characteristics of star sensor attitude error and the mathematical derivation, the SINS navigation errors can be obtained by the coupling calculation between the SINS attitude and the star sensor attitude. Unlike several current techniques, the navigation process of this method is non-radiating and invulnerable to jamming. The effectiveness of this approach was demonstrated by simulation and experimental study. The results show that this integrated navigation method can estimate the attitude error and the position error of SINS. Therefore, the SINS navigation accuracy is improved. (paper)

A study of compositional verification based IMA integration method

Science.gov (United States)

Huang, Hui; Zhang, Guoquan; Xu, Wanmeng

2018-03-01

The rapid development of avionics systems is driving the application of integrated modular avionics (IMA) systems. But meanwhile it is improving avionics system integration, complexity of system test. Then we need simplify the method of IMA system test. The IMA system supports a module platform that runs multiple applications, and shares processing resources. Compared with federated avionics system, IMA system is difficult to isolate failure. Therefore, IMA system verification will face the critical problem is how to test shared resources of multiple application. For a simple avionics system, traditional test methods are easily realizing to test a whole system. But for a complex system, it is hard completed to totally test a huge and integrated avionics system. Then this paper provides using compositional-verification theory in IMA system test, so that reducing processes of test and improving efficiency, consequently economizing costs of IMA system integration.

New design for photonic temporal integration with combined high processing speed and long operation time window.

Science.gov (United States)

Asghari, Mohammad H; Park, Yongwoo; Azaña, José

2011-01-17

We propose and experimentally prove a novel design for implementing photonic temporal integrators simultaneously offering a high processing bandwidth and a long operation time window, namely a large time-bandwidth product. The proposed scheme is based on concatenating in series a time-limited ultrafast photonic temporal integrator, e.g. implemented using a fiber Bragg grating (FBG), with a discrete-time (bandwidth limited) optical integrator, e.g. implemented using an optical resonant cavity. This design combines the advantages of these two previously demonstrated photonic integrator solutions, providing a processing speed as high as that of the time-limited ultrafast integrator and an operation time window fixed by the discrete-time integrator. Proof-of-concept experiments are reported using a uniform fiber Bragg grating (as the original time-limited integrator) connected in series with a bulk-optics coherent interferometers' system (as a passive 4-points discrete-time photonic temporal integrator). Using this setup, we demonstrate accurate temporal integration of complex-field optical signals with time-features as fast as ~6 ps, only limited by the processing bandwidth of the FBG integrator, over time durations as long as ~200 ps, which represents a 4-fold improvement over the operation time window (~50 ps) of the original FBG integrator.

Higher-Order Integral Equation Methods in Computational Electromagnetics

DEFF Research Database (Denmark)

Jørgensen, Erik; Meincke, Peter

Higher-order integral equation methods have been investigated. The study has focused on improving the accuracy and efficiency of the Method of Moments (MoM) applied to electromagnetic problems. A new set of hierarchical Legendre basis functions of arbitrary order is developed. The new basis...

Methods for determining time of death.

Science.gov (United States)

Madea, Burkhard

2016-12-01

Medicolegal death time estimation must estimate the time since death reliably. Reliability can only be provided empirically by statistical analysis of errors in field studies. Determining the time since death requires the calculation of measurable data along a time-dependent curve back to the starting point. Various methods are used to estimate the time since death. The current gold standard for death time estimation is a previously established nomogram method based on the two-exponential model of body cooling. Great experimental and practical achievements have been realized using this nomogram method. To reduce the margin of error of the nomogram method, a compound method was developed based on electrical and mechanical excitability of skeletal muscle, pharmacological excitability of the iris, rigor mortis, and postmortem lividity. Further increasing the accuracy of death time estimation involves the development of conditional probability distributions for death time estimation based on the compound method. Although many studies have evaluated chemical methods of death time estimation, such methods play a marginal role in daily forensic practice. However, increased precision of death time estimation has recently been achieved by considering various influencing factors (i.e., preexisting diseases, duration of terminal episode, and ambient temperature). Putrefactive changes may be used for death time estimation in water-immersed bodies. Furthermore, recently developed technologies, such as H magnetic resonance spectroscopy, can be used to quantitatively study decompositional changes. This review addresses the gold standard method of death time estimation in forensic practice and promising technological and scientific developments in the field.

The flow measurement methods for the primary system of integral reactors

International Nuclear Information System (INIS)

Lee, J.; Seo, J. K.; Lee, D. J.

2001-01-01

It is the common features of the integral reactors that the main components of the primary system are installed within the reactor vessel, and so there are no any flow pipes connecting the reactor coolant pumps or steam generators. Due to no any flow pipes, it is impossible to measure the differential pressure at the primary system of the integral reactors, and it also makes impossible measure the primary coolant flow rate. The objective of the study is to draw up the flow measurement methods for the primary system of integral reactors. As a result of the review, we have made a selection of the flow measurement method by pump speed, bt HBM, and by pump motor power as the flow measurement methods for the primary system of integral reactors. Peculiarly, we did not found out a precedent which the direct pump motor power-flow rate curve is used as the flow measurement method in the existing commercial nuclear power reactors. Therefore, to use this method for integral reactors, it is needed to bear the follow-up measures in mind. The follow-up measures is included in this report

Orientation, Evaluation, and Integration of Part-Time Nursing Faculty.

Science.gov (United States)

Carlson, Joanne S

2015-07-10

This study helps to quantify and describe orientation, evaluation, and integration practices pertaining to part-time clinical nursing faculty teaching in prelicensure nursing education programs. A researcher designed Web-based survey was used to collect information from a convenience sample of part-time clinical nursing faculty teaching in prelicensure nursing programs. Survey questions focused on the amount and type of orientation, evaluation, and integration practices. Descriptive statistics were used to analyze results. Respondents reported on average four hours of orientation, with close to half reporting no more than two hours. Evaluative feedback was received much more often from students than from full-time faculty. Most respondents reported receiving some degree of mentoring and that it was easy to get help from full-time faculty. Respondents reported being most informed about student evaluation procedures, grading, and the steps to take when students are not meeting course objectives, and less informed about changes to ongoing curriculum and policy.

A systematic and efficient method to compute multi-loop master integrals

Directory of Open Access Journals (Sweden)

Xiao Liu

2018-04-01

Full Text Available We propose a novel method to compute multi-loop master integrals by constructing and numerically solving a system of ordinary differential equations, with almost trivial boundary conditions. Thus it can be systematically applied to problems with arbitrary kinematic configurations. Numerical tests show that our method can not only achieve results with high precision, but also be much faster than the only existing systematic method sector decomposition. As a by product, we find a new strategy to compute scalar one-loop integrals without reducing them to master integrals.

Method of mechanical quadratures for solving singular integral equations of various types

Science.gov (United States)

Sahakyan, A. V.; Amirjanyan, H. A.

2018-04-01

The method of mechanical quadratures is proposed as a common approach intended for solving the integral equations defined on finite intervals and containing Cauchy-type singular integrals. This method can be used to solve singular integral equations of the first and second kind, equations with generalized kernel, weakly singular equations, and integro-differential equations. The quadrature rules for several different integrals represented through the same coefficients are presented. This allows one to reduce the integral equations containing integrals of different types to a system of linear algebraic equations.

A space-time mixed galerkin marching-on-in-time scheme for the time-domain combined field integral equation

KAUST Repository

Beghein, Yves

2013-03-01

The time domain combined field integral equation (TD-CFIE), which is constructed from a weighted sum of the time domain electric and magnetic field integral equations (TD-EFIE and TD-MFIE) for analyzing transient scattering from closed perfect electrically conducting bodies, is free from spurious resonances. The standard marching-on-in-time technique for discretizing the TD-CFIE uses Galerkin and collocation schemes in space and time, respectively. Unfortunately, the standard scheme is theoretically not well understood: stability and convergence have been proven for only one class of space-time Galerkin discretizations. Moreover, existing discretization schemes are nonconforming, i.e., the TD-MFIE contribution is tested with divergence conforming functions instead of curl conforming functions. We therefore introduce a novel space-time mixed Galerkin discretization for the TD-CFIE. A family of temporal basis and testing functions with arbitrary order is introduced. It is explained how the corresponding interactions can be computed efficiently by existing collocation-in-time codes. The spatial mixed discretization is made fully conforming and consistent by leveraging both Rao-Wilton-Glisson and Buffa-Christiansen basis functions and by applying the appropriate bi-orthogonalization procedures. The combination of both techniques is essential when high accuracy over a broad frequency band is required. © 2012 IEEE.

Feasibility of the integration of CRONOS, a 3-D neutronics code, into real-time simulators

International Nuclear Information System (INIS)

Ragusa, J.C.

2001-01-01

In its effort to contribute to nuclear power plant safety, CEA proposes the integration of an engineering grade 3-D neutronics code into a real-time plant analyser. This paper describes the capabilities of the neutronics code CRONOS to achieve a fast running performance. First, we will present current core models in simulators and explain their drawbacks. Secondly, the mean features of CRONOS's spatial-kinetics methods will be reviewed. We will then present an optimum core representation with respect to mesh size, choice of finite elements (FE) basis and execution time, for accurate results as well as the multi 1-D thermal-hydraulics (T/H) model developed to take into account 3-D effects in updating the cross-sections. A Main Steam Line Break (MSLB) End-of-Life (EOL) Hot-Zero-Power (HZP) accident will be used as an example, before we conclude with the perspectives of integrating CRONOS's 3-D core model into real-time simulators. (author)

Feasibility of the integration of CRONOS, a 3-D neutronics code, into real-time simulators

Energy Technology Data Exchange (ETDEWEB)

Ragusa, J.C. [CEA Saclay, Dept. de Mecanique et de Technologie, 91 - Gif-sur-Yvette (France)

2001-07-01

In its effort to contribute to nuclear power plant safety, CEA proposes the integration of an engineering grade 3-D neutronics code into a real-time plant analyser. This paper describes the capabilities of the neutronics code CRONOS to achieve a fast running performance. First, we will present current core models in simulators and explain their drawbacks. Secondly, the mean features of CRONOS's spatial-kinetics methods will be reviewed. We will then present an optimum core representation with respect to mesh size, choice of finite elements (FE) basis and execution time, for accurate results as well as the multi 1-D thermal-hydraulics (T/H) model developed to take into account 3-D effects in updating the cross-sections. A Main Steam Line Break (MSLB) End-of-Life (EOL) Hot-Zero-Power (HZP) accident will be used as an example, before we conclude with the perspectives of integrating CRONOS's 3-D core model into real-time simulators. (author)

Nonlinear structural analysis using integrated force method

Indian Academy of Sciences (India)

A new formulation termed the Integrated Force Method (IFM) was proposed by Patnaik ... nated ``Structure (nY m)'' where (nY m) are the force and displacement degrees of ..... Patnaik S N, Yadagiri S 1976 Frequency analysis of structures.

Deterministic factor analysis: methods of integro-differentiation of non-integral order

Directory of Open Access Journals (Sweden)

Valentina V. Tarasova

2016-12-01

Full Text Available Objective to summarize the methods of deterministic factor economic analysis namely the differential calculus and the integral method. nbsp Methods mathematical methods for integrodifferentiation of nonintegral order the theory of derivatives and integrals of fractional nonintegral order. Results the basic concepts are formulated and the new methods are developed that take into account the memory and nonlocality effects in the quantitative description of the influence of individual factors on the change in the effective economic indicator. Two methods are proposed for integrodifferentiation of nonintegral order for the deterministic factor analysis of economic processes with memory and nonlocality. It is shown that the method of integrodifferentiation of nonintegral order can give more accurate results compared with standard methods method of differentiation using the first order derivatives and the integral method using the integration of the first order for a wide class of functions describing effective economic indicators. Scientific novelty the new methods of deterministic factor analysis are proposed the method of differential calculus of nonintegral order and the integral method of nonintegral order. Practical significance the basic concepts and formulas of the article can be used in scientific and analytical activity for factor analysis of economic processes. The proposed method for integrodifferentiation of nonintegral order extends the capabilities of the determined factorial economic analysis. The new quantitative method of deterministic factor analysis may become the beginning of quantitative studies of economic agents behavior with memory hereditarity and spatial nonlocality. The proposed methods of deterministic factor analysis can be used in the study of economic processes which follow the exponential law in which the indicators endogenous variables are power functions of the factors exogenous variables including the processes

Undergraduate medical student's perception about an integrated method of teaching at a medical school in Oman

Directory of Open Access Journals (Sweden)

Harshal Sabane

2015-01-01

Full Text Available Objective In recent years, there has been a gradual but definitive shift in medical schools all over the globe to promote a more integrated way of teaching. Integration of medical disciplines promotes a holistic understanding of the medical curriculum in the students. This helps them better understand and appreciate the importance and role of each medical subject. Method The study was conducted among the 5th year Pre-clinical students. Questionnaire consisted of 4 questions on the level of integration, 5 questions on various aspects of the assessment and some questions which tested the level of awareness of the integrated method. Result Out of a total of 72 students present on the day of data collection, 65 participated in the study giving a response rate of 90.27 %. After primary data cleansing 4 questionnaires had to be omitted. Most of the students opined as “good” or “very good” for the questions on integration and its attributes. Only 27 (44 % were aware of integrated curriculum being taught in other medical schools in the gulf. Similar findings were observed regarding assessment related questions. Reduction in the number of block exams is unpopular among the students and only 6% have agreed for 3, 4, or 5 non-summative block assessments. Opinion regarding the help of integrated teaching in IFOM based OMSB entrance examination was mixed with a greater variance in the responses. 43% students have indicated that they would like to spend more time with PDCI. Conclusion The students of our institution seem to have a favourable opinion regarding the integrated system of teaching. The satisfaction with the conduct of examinations and its related variables is found to be high. A reduction in the number of block exams however is unpopular among the target group and they would appreciate a greater time allocation for subjects of PDCI and Pharmacology.

Analysis of noise in energy-dispersive spectrometers using time-domain methods

CERN Document Server

Goulding, F S

2002-01-01

This paper presents an integrated time domain approach to the optimization of the signal-to-noise ratio in all spectrometer systems that contain a detector that converts incoming quanta of radiation into electrical pulse signals that are amplified and shaped by an electronic pulse shaper. It allows analysis of normal passive pulse shapers as well as time-variant systems where switching of shaping elements occurs in synchronism with the signal. It also deals comfortably with microcalorimeters (sometimes referred to as bolometers), where noise-determining elements, such as the temperature-sensing element's resistance and temperature, change with time in the presence of a signal. As part of the purely time-domain approach, a new method of calculating the Johnson noise in resistors using only the statistics of electron motion is presented. The result is a time-domain analog of the Nyquist formula.

H∞ Control of Coronary Artery Input Time-Delay System via the Free-Matrix-Based Integral Inequality

Directory of Open Access Journals (Sweden)

Sha-sha Li

2018-01-01

Full Text Available The issue of H∞ control for the coronary artery input time-delay system with external disturbance is of concern. To further reduce conservation, we utilize the free-matrix-based integral inequality, Wirtinger-based integral inequality, and reciprocal convex combination approach to construct Lyapunov-Krasovskii function (LKF. Then a sufficient condition for controller design which can guarantee robust synchronization the coronary artery system is represented in terms of linear matrix inequality (LMI. Finally, a numerical example is exploited to show the effectiveness of the proposed methods.

Integrated response and transit time distributions of watersheds by combining hydrograph separation and long-term transit time modeling

Directory of Open Access Journals (Sweden)

M. C. Roa-García

2010-08-01

Full Text Available We present a new modeling approach analyzing and predicting the Transit Time Distribution (TTD and the Response Time Distribution (RTD from hourly to annual time scales as two distinct hydrological processes. The model integrates Isotope Hydrograph Separation (IHS and the Instantaneous Unit Hydrograph (IUH approach as a tool to provide a more realistic description of transit and response time of water in catchments. Individual event simulations and parameterizations were combined with long-term baseflow simulation and parameterizations; this provides a comprehensive picture of the catchment response for a long time span for the hydraulic and isotopic processes. The proposed method was tested in three Andean headwater catchments to compare the effects of land use on hydrological response and solute transport. Results show that the characteristics of events and antecedent conditions have a significant influence on TTD and RTD, but in general the RTD of the grassland dominated catchment is concentrated in the shorter time spans and has a higher cumulative TTD, while the forest dominated catchment has a relatively higher response distribution and lower cumulative TTD. The catchment where wetlands concentrate shows a flashier response, but wetlands also appear to prolong transit time.

Integrals of random fields treated by the model correction factor method

DEFF Research Database (Denmark)

Franchin, P.; Ditlevsen, Ove Dalager; Kiureghian, Armen Der

2002-01-01

The model correction factor method (MCFM) is used in conjunction with the first-order reliability method (FORM) to solve structural reliability problems involving integrals of non-Gaussian random fields. The approach replaces the limit-state function with an idealized one, in which the integrals ...

An Integrated Method for Airfoil Optimization

Science.gov (United States)

Okrent, Joshua B.

Design exploration and optimization is a large part of the initial engineering and design process. To evaluate the aerodynamic performance of a design, viscous Navier-Stokes solvers can be used. However this method can prove to be overwhelmingly time consuming when performing an initial design sweep. Therefore, another evaluation method is needed to provide accurate results at a faster pace. To accomplish this goal, a coupled viscous-inviscid method is used. This thesis proposes an integrated method for analyzing, evaluating, and optimizing an airfoil using a coupled viscous-inviscid solver along with a genetic algorithm to find the optimal candidate. The method proposed is different from prior optimization efforts in that it greatly broadens the design space, while allowing the optimization to search for the best candidate that will meet multiple objectives over a characteristic mission profile rather than over a single condition and single optimization parameter. The increased design space is due to the use of multiple parametric airfoil families, namely the NACA 4 series, CST family, and the PARSEC family. Almost all possible airfoil shapes can be created with these three families allowing for all possible configurations to be included. This inclusion of multiple airfoil families addresses a possible criticism of prior optimization attempts since by only focusing on one airfoil family, they were inherently limiting the number of possible airfoil configurations. By using multiple parametric airfoils, it can be assumed that all reasonable airfoil configurations are included in the analysis and optimization and that a global and not local maximum is found. Additionally, the method used is amenable to customization to suit any specific needs as well as including the effects of other physical phenomena or design criteria and/or constraints. This thesis found that an airfoil configuration that met multiple objectives could be found for a given set of nominal

Measuring fragmentation in dissociative identity disorder: the integration measure and relationship to switching and time in therapy

Directory of Open Access Journals (Sweden)

Margaret Rose Barlow

2014-01-01

Full Text Available Background: Some people with dissociative identity disorder (DID have very little communication or awareness among the parts of their identity, while others experience a great deal of cooperation among alternate identities. Previous research on this topic has been sparse. Currently, there is no empirical measure of integration versus fragmentation in a person with DID. In this study, we report the development of such a measure. Objective: The goal of this study was to pilot the integration measure (IM and to address its psychometric properties and relationships to other measures. The IM is the first standardized measure of integration in DID. Method: Eleven women with DID participated in an experiment that included a variety of tasks. They filled out questionnaires about trauma and dissociation as well as the IM. They also provided verbal results about switching among alternate identities during the study sessions. Results: Participants switched among identities an average of 5.8 times during the first session, and switching was highly correlated with trauma. Integration was related to switching, though this relationship may be non-linear. Integration was not related to time in psychotherapy. Conclusions: The IM provides a useful beginning to quantify and study integration and fragmentation in DID. Directions for future research are also discussed, including expanding the IM from this pilot. The IM may be useful in treatment settings to assess progress or change over time.

Permutation statistical methods an integrated approach

CERN Document Server

Berry, Kenneth J; Johnston, Janis E

2016-01-01

This research monograph provides a synthesis of a number of statistical tests and measures, which, at first consideration, appear disjoint and unrelated. Numerous comparisons of permutation and classical statistical methods are presented, and the two methods are compared via probability values and, where appropriate, measures of effect size. Permutation statistical methods, compared to classical statistical methods, do not rely on theoretical distributions, avoid the usual assumptions of normality and homogeneity of variance, and depend only on the data at hand. This text takes a unique approach to explaining statistics by integrating a large variety of statistical methods, and establishing the rigor of a topic that to many may seem to be a nascent field in statistics. This topic is new in that it took modern computing power to make permutation methods available to people working in the mainstream of research. This research monograph addresses a statistically-informed audience, and can also easily serve as a ...

Time delayed Ensemble Nudging Method

Science.gov (United States)

An, Zhe; Abarbanel, Henry

Optimal nudging method based on time delayed embedding theory has shows potentials on analyzing and data assimilation in previous literatures. To extend the application and promote the practical implementation, new nudging assimilation method based on the time delayed embedding space is presented and the connection with other standard assimilation methods are studied. Results shows the incorporating information from the time series of data can reduce the sufficient observation needed to preserve the quality of numerical prediction, making it a potential alternative in the field of data assimilation of large geophysical models.

Design of time-pulse coded optoelectronic neuronal elements for nonlinear transformation and integration

Science.gov (United States)

Krasilenko, Vladimir G.; Nikolsky, Alexander I.; Lazarev, Alexander A.; Lazareva, Maria V.

2008-03-01

In the paper the actuality of neurophysiologically motivated neuron arrays with flexibly programmable functions and operations with possibility to select required accuracy and type of nonlinear transformation and learning are shown. We consider neurons design and simulation results of multichannel spatio-time algebraic accumulation - integration of optical signals. Advantages for nonlinear transformation and summation - integration are shown. The offered circuits are simple and can have intellectual properties such as learning and adaptation. The integrator-neuron is based on CMOS current mirrors and comparators. The performance: consumable power - 100...500 μW, signal period- 0.1...1ms, input optical signals power - 0.2...20 μW time delays - less 1μs, the number of optical signals - 2...10, integration time - 10...100 of signal periods, accuracy or integration error - about 1%. Various modifications of the neuron-integrators with improved performance and for different applications are considered in the paper.

Integration of equations of parabolic type by the method of nets

CERN Document Server

Saul'Yev, V K; Stark, M; Ulam, S

1964-01-01

International Series of Monographs in Pure and Applied Mathematics, Volume 54: Integration of Equations of Parabolic Type by the Method of Nets deals with solving parabolic partial differential equations using the method of nets. The first part of this volume focuses on the construction of net equations, with emphasis on the stability and accuracy of the approximating net equations. The method of nets or method of finite differences (used to define the corresponding numerical method in ordinary differential equations) is one of many different approximate methods of integration of partial diff

Time series analysis of the developed financial markets' integration using visibility graphs

Science.gov (United States)

Zhuang, Enyu; Small, Michael; Feng, Gang

2014-09-01

A time series representing the developed financial markets' segmentation from 1973 to 2012 is studied. The time series reveals an obvious market integration trend. To further uncover the features of this time series, we divide it into seven windows and generate seven visibility graphs. The measuring capabilities of the visibility graphs provide means to quantitatively analyze the original time series. It is found that the important historical incidents that influenced market integration coincide with variations in the measured graphical node degree. Through the measure of neighborhood span, the frequencies of the historical incidents are disclosed. Moreover, it is also found that large "cycles" and significant noise in the time series are linked to large and small communities in the generated visibility graphs. For large cycles, how historical incidents significantly affected market integration is distinguished by density and compactness of the corresponding communities.

The 3D Lagrangian Integral Method. Henrik Koblitz Rasmussen

DEFF Research Database (Denmark)

Rasmussen, Henrik Koblitz

2003-01-01

. This are processes such as thermo-forming, gas-assisted injection moulding and all kind of simultaneous multi-component polymer processing operations. Though, in all polymer processing operations free surfaces (or interfaces) are present and the dynamic of these surfaces are of interest. In the "3D Lagrangian...... Integral Method" to simulate viscoelastic flow, the governing equations are solved for the particle positions (Lagrangian kinematics). Therefore, the transient motion of surfaces can be followed in a particularly simple fashion even in 3D viscoelastic flow. The "3D Lagrangian Integral Method" is described...

Time Varying Market Integration and Expected Rteurns in Emerging Markets

NARCIS (Netherlands)

de Jong, F.C.J.M.; de Roon, F.A.

2001-01-01

We use a simple model in which the expected returns in emerging markets depend on their systematic risk as measured by their beta relative to the world portfolio as well as on the level of integration in that market.The level of integration is a time-varying variable that depends on the market value

A hybrid data fusion method for GNSS/INS integration navigation system

Science.gov (United States)

Yang, Ling; Li, Bofeng; Shen, Yunzhong; Li, Haojun

2017-04-01

Although DGNSS is widely used and PPP-GNSS is nowadays a viable precise positioning technology option, the major disadvantage of GNSS still remains: signal blockage due to obstructions in urban and built up environments, and extreme power attenuation of the signals when operated indoors. The combination of GNSS with other sensors, such as a self-contained inertial navigation system (INS), provides an ideal position and attitude determination solution which can not only mitigate the weakness of GNSS, but also bound the INS error that otherwise would grow with time when the INS operates alone. However, the navigation accuracy provided by GNSS/INS strongly depends on the quality and geometry of the GNSS observations, the quality of the INS technology used, and the integration model applied. There are two main types of coupled schemes for integration systems: loosely coupled integration and tightly coupled integration. In loosely coupled integration, position measurements are taken from both systems and combined optimally, usually in a Kalman filter. Tightly coupled integration directly combines the raw pseudorange or carrier phase measurements of GNSS with inertial measurements in an extended Kalman filter. The latter technique improves the ability to resolve ambiguities, i.e. allows a quicker recovery from outage events such as a loss of signal under vegetation. In recent years, tightly coupled differential carrier phase GNSS/INS integration has become popular, because it has the advantage of providing accurate position information even when GPS measurements are rank-deficient in stand-alone processing and is theoretically optimal in a filtering sense, especially in urban navigation applications. However, the heavier computational burden and sensor communication usually complicate the tightly coupled integration and reduce the system efficiency, compared with the loosely coupled integration. In this paper, it has been proved that the loosely coupled and tightly

The comparative analysis of the current-meter method and the pressure-time method used for discharge measurements in the Kaplan turbine penstocks

Science.gov (United States)

Adamkowski, A.; Krzemianowski, Z.

2012-11-01

The paper presents experiences gathered during many years of utilizing the current-meter and pressure-time methods for flow rate measurements in many hydropower plants. The integration techniques used in these both methods are different from the recommendations contained in the relevant international standards, mainly from the graphical and arithmetical ones. The results of the comparative analysis of both methods applied at the same time during the hydraulic performance tests of two Kaplan turbines in one of the Polish hydropower plant are presented in the final part of the paper. In the case of the pressure-time method application, the concrete penstocks of the tested turbines required installing a special measuring instrumentation inside the penstock. The comparison has shown a satisfactory agreement between the results of discharge measurements executed using the both considered methods. Maximum differences between the discharge values have not exceeded 1.0 % and the average differences have not been greater than 0.5 %.

The comparative analysis of the current-meter method and the pressure-time method used for discharge measurements in the Kaplan turbine penstocks

International Nuclear Information System (INIS)

Adamkowski, A; Krzemianowski, Z

2012-01-01

The paper presents experiences gathered during many years of utilizing the current-meter and pressure-time methods for flow rate measurements in many hydropower plants. The integration techniques used in these both methods are different from the recommendations contained in the relevant international standards, mainly from the graphical and arithmetical ones. The results of the comparative analysis of both methods applied at the same time during the hydraulic performance tests of two Kaplan turbines in one of the Polish hydropower plant are presented in the final part of the paper. In the case of the pressure-time method application, the concrete penstocks of the tested turbines required installing a special measuring instrumentation inside the penstock. The comparison has shown a satisfactory agreement between the results of discharge measurements executed using the both considered methods. Maximum differences between the discharge values have not exceeded 1.0 % and the average differences have not been greater than 0.5 %.

Equilibrium and response properties of the integrate-and-fire neuron in discrete time

Directory of Open Access Journals (Sweden)

Moritz Helias

2010-01-01

Full Text Available The integrate-and-fire neuron with exponential postsynaptic potentials is a frequently employed model to study neural networks. Simulations in discrete time still have highest performance at moderate numerical errors, which makes them first choice for long-term simulations of plastic networks. Here we extend the population density approach to investigate how the equilibrium and response properties of the leaky integrate-and-fire neuron are affected by time discretization. We present a novel analytical treatment of the boundary condition at threshold, taking both discretization of time and finite synaptic weights into account. We uncover an increased membrane potential density just below threshold as the decisive property that explains the deviations found between simulations and the classical diffusion approximation. Temporal discretization and finite synaptic weights both contribute to this effect. Our treatment improves the standard formula to calculate the neuron’s equilibrium firing rate. Direct solution of the Markov process describing the evolution of the membrane potential density confirms our analysis and yields a method to calculate the firing rate exactly. Knowing the shape of the membrane potential distribution near threshold enables us to devise the transient response properties of the neuron model to synaptic input. We find a pronounced non-linear fast response component that has not been described by the prevailing continuous time theory for Gaussian white noise input.

Research on Multi - Person Parallel Modeling Method Based on Integrated Model Persistent Storage

Science.gov (United States)

Qu, MingCheng; Wu, XiangHu; Tao, YongChao; Liu, Ying

2018-03-01

This paper mainly studies the multi-person parallel modeling method based on the integrated model persistence storage. The integrated model refers to a set of MDDT modeling graphics system, which can carry out multi-angle, multi-level and multi-stage description of aerospace general embedded software. Persistent storage refers to converting the data model in memory into a storage model and converting the storage model into a data model in memory, where the data model refers to the object model and the storage model is a binary stream. And multi-person parallel modeling refers to the need for multi-person collaboration, the role of separation, and even real-time remote synchronization modeling.

Distributed finite-time containment control for double-integrator multiagent systems.

Science.gov (United States)

Wang, Xiangyu; Li, Shihua; Shi, Peng

2014-09-01

In this paper, the distributed finite-time containment control problem for double-integrator multiagent systems with multiple leaders and external disturbances is discussed. In the presence of multiple dynamic leaders, by utilizing the homogeneous control technique, a distributed finite-time observer is developed for the followers to estimate the weighted average of the leaders' velocities at first. Then, based on the estimates and the generalized adding a power integrator approach, distributed finite-time containment control algorithms are designed to guarantee that the states of the followers converge to the dynamic convex hull spanned by those of the leaders in finite time. Moreover, as a special case of multiple dynamic leaders with zero velocities, the proposed containment control algorithms also work for the case of multiple stationary leaders without using the distributed observer. Simulations demonstrate the effectiveness of the proposed control algorithms.

Integrated optical delay lines for time-division multiplexers

NARCIS (Netherlands)

Stopinski, S.T.; Malinowski, M.; Piramidowicz, R.; Kleijn, E.; Smit, M.K.; Leijtens, X.J.M.

2013-01-01

In this paper, we present a study of integrated optical delay lines (DLs) for application in optical time-division multiplexers. The investigated DLs are formed by spirally folded waveguides. The components were designed in a generic approach and fabricated in multi-project wafer runs on an

Integration and timing of basic and clinical sciences education.

Science.gov (United States)

Bandiera, Glen; Boucher, Andree; Neville, Alan; Kuper, Ayelet; Hodges, Brian

2013-05-01

Medical education has traditionally been compartmentalized into basic and clinical sciences, with the latter being viewed as the skillful application of the former. Over time, the relevance of basic sciences has become defined by their role in supporting clinical problem solving rather than being, of themselves, a defining knowledge base of physicians. As part of the national Future of Medical Education in Canada (FMEC MD) project, a comprehensive empirical environmental scan identified the timing and integration of basic sciences as a key pressing issue for medical education. Using the literature review, key informant interviews, stakeholder meetings, and subsequent consultation forums from the FMEC project, this paper details the empirical basis for focusing on the role of basic science, the evidentiary foundations for current practices, and the implications for medical education. Despite a dearth of definitive relevant studies, opinions about how best to integrate the sciences remain strong. Resource allocation, political power, educational philosophy, and the shift from a knowledge-based to a problem-solving profession all influence the debate. There was little disagreement that both sciences are important, that many traditional models emphasized deep understanding of limited basic science disciplines at the expense of other relevant content such as social sciences, or that teaching the sciences contemporaneously rather than sequentially has theoretical and practical merit. Innovations in integrated curriculum design have occurred internationally. Less clear are the appropriate balance of the sciences, the best integration model, and solutions to the political and practical challenges of integrated curricula. New curricula tend to emphasize integration, development of more diverse physician competencies, and preparation of physicians to adapt to evolving technology and patients' expectations. Refocusing the basic/clinical dichotomy to a foundational

Generalized Time-Limited Balanced Reduction Method

DEFF Research Database (Denmark)

Shaker, Hamid Reza; Shaker, Fatemeh

2013-01-01

In this paper, a new method for model reduction of bilinear systems is presented. The proposed technique is from the family of gramian-based model reduction methods. The method uses time-interval generalized gramians in the reduction procedure rather than the ordinary generalized gramians...... and in such a way it improves the accuracy of the approximation within the time-interval which the method is applied. The time-interval generalized gramians are the solutions to the generalized time-interval Lyapunov equations. The conditions for these equations to be solvable are derived and an algorithm...

Transient Thermal Analysis of 3-D Integrated Circuits Packages by the DGTD Method

KAUST Repository

Li, Ping

2017-03-11

Since accurate thermal analysis plays a critical role in the thermal design and management of the 3-D system-level integration, in this paper, a discontinuous Galerkin time-domain (DGTD) algorithm is proposed to achieve this purpose. Such as the parabolic partial differential equation (PDE), the transient thermal equation cannot be directly solved by the DGTD method. To address this issue, the heat flux, as an auxiliary variable, is introduced to reduce the Laplace operator to a divergence operator. The resulting PDE is hyperbolic, which can be further written into a conservative form. By properly choosing the definition of the numerical flux used for the information exchange between neighboring elements, the hyperbolic thermal PDE can be solved by the DGTD together with the auxiliary differential equation. The proposed algorithm is a kind of element-level domain decomposition method, which is suitable to deal with multiscale geometries in 3-D integrated systems. To verify the accuracy and robustness of the developed DGTD algorithm, several representative examples are benchmarked.

A Digitally Programmable Differential Integrator with Enlarged Time Constant

Directory of Open Access Journals (Sweden)

S. K. Debroy

1994-12-01

Full Text Available A new Operational Amplifier (OA-RC integrator network is described. The novelties of the design are used of single grounded capacitor, ideal integration function realization with dual-input capability and design flexibility for extremely large time constant involving an enlargement factor (K using product of resistor ratios. The aspect of the digital control of K through a programmable resistor array (PRA controlled by a microprocessor has also been implemented. The effect of the OA-poles has been analyzed which indicates degradation of the integrator-Q at higher frequencies. An appropriate Q-compensation design scheme exhibiting 1 : |A|2 order of Q-improvement has been proposed with supporting experimental observations.

Mathematical methods linear algebra normed spaces distributions integration

CERN Document Server

Korevaar, Jacob

1968-01-01

Mathematical Methods, Volume I: Linear Algebra, Normed Spaces, Distributions, Integration focuses on advanced mathematical tools used in applications and the basic concepts of algebra, normed spaces, integration, and distributions.The publication first offers information on algebraic theory of vector spaces and introduction to functional analysis. Discussions focus on linear transformations and functionals, rectangular matrices, systems of linear equations, eigenvalue problems, use of eigenvectors and generalized eigenvectors in the representation of linear operators, metric and normed vector

New methods for the numerical integration of ordinary differential equations and their application to the equations of motion of spacecraft

Science.gov (United States)

Banyukevich, A.; Ziolkovski, K.

1975-01-01

A number of hybrid methods for solving Cauchy problems are described on the basis of an evaluation of advantages of single and multiple-point numerical integration methods. The selection criterion is the principle of minimizing computer time. The methods discussed include the Nordsieck method, the Bulirsch-Stoer extrapolation method, and the method of recursive Taylor-Steffensen power series.

Application of Stochastic Sensitivity Analysis to Integrated Force Method

Directory of Open Access Journals (Sweden)

X. F. Wei

2012-01-01

Full Text Available As a new formulation in structural analysis, Integrated Force Method has been successfully applied to many structures for civil, mechanical, and aerospace engineering due to the accurate estimate of forces in computation. Right now, it is being further extended to the probabilistic domain. For the assessment of uncertainty effect in system optimization and identification, the probabilistic sensitivity analysis of IFM was further investigated in this study. A set of stochastic sensitivity analysis formulation of Integrated Force Method was developed using the perturbation method. Numerical examples are presented to illustrate its application. Its efficiency and accuracy were also substantiated with direct Monte Carlo simulations and the reliability-based sensitivity method. The numerical algorithm was shown to be readily adaptable to the existing program since the models of stochastic finite element and stochastic design sensitivity are almost identical.

The integrated circuit IC EMP transient state disturbance effect experiment method investigates

International Nuclear Information System (INIS)

Li Xiaowei

2004-01-01

Transient state disturbance characteristic study on the integrated circuit, IC, need from its coupling path outset. Through cable (aerial) coupling, EMP converts to an pulse current voltage and results in the impact to the integrated circuit I/O orifice passing the cable. Aiming at the armament system construction feature, EMP effect to the integrated circuit, IC inside the system is analyzed. The integrated circuit, IC EMP effect experiment current injection method is investigated and a few experiments method is given. (authors)

Nonlinear Fredholm Integral Equation of the Second Kind with Quadrature Methods

Directory of Open Access Journals (Sweden)

M. Jafari Emamzadeh

2010-06-01

Full Text Available In this paper, a numerical method for solving the nonlinear Fredholm integral equation is presented. We intend to approximate the solution of this equation by quadrature methods and by doing so, we solve the nonlinear Fredholm integral equation more accurately. Several examples are given at the end of this paper

Integrating external biological knowledge in the construction of regulatory networks from time-series expression data

Directory of Open Access Journals (Sweden)

Lo Kenneth

2012-08-01

Full Text Available Abstract Background Inference about regulatory networks from high-throughput genomics data is of great interest in systems biology. We present a Bayesian approach to infer gene regulatory networks from time series expression data by integrating various types of biological knowledge. Results We formulate network construction as a series of variable selection problems and use linear regression to model the data. Our method summarizes additional data sources with an informative prior probability distribution over candidate regression models. We extend the Bayesian model averaging (BMA variable selection method to select regulators in the regression framework. We summarize the external biological knowledge by an informative prior probability distribution over the candidate regression models. Conclusions We demonstrate our method on simulated data and a set of time-series microarray experiments measuring the effect of a drug perturbation on gene expression levels, and show that it outperforms leading regression-based methods in the literature.

An Efficient Method to Search Real-Time Bulk Data for an Information Processing System

International Nuclear Information System (INIS)

Kim, Seong Jin; Kim, Jong Myung; Suh, Yong Suk; Keum, Jong Yong; Park, Heui Youn

2005-01-01

The Man Machine Interface System (MMIS) of System-integrated Modular Advanced ReacTor (SMART) is designed with fully digitalized features. The Information Processing System (IPS) of the MMIS acquires and processes plant data from other systems. In addition, the IPS provides plant operation information to operators in the control room. The IPS is required to process bulky data in a real-time. So, it is necessary to consider a special processing method with regards to flexibility and performance because more than a few thousands of Plant Information converges on the IPS. Among other things, the processing time for searching for data from the bulk data consumes much more than other the processing times. Thus, this paper explores an efficient method for the search and examines its feasibility

Integrated vendor-buyer inventory models with inflation and time value of money in controllable lead time

Directory of Open Access Journals (Sweden)

Prashant Jindal

2016-01-01

Full Text Available In the global critical economic scenario, inflation plays a vital role in deciding optimal pricing of goods in any business entity. This article presents two single-vendor single-buyer integrated supply chain inventory models with inflation and time value of money. Shortage is allowed during the lead-time and it is partially backlogged. Lead time is controllable and can be reduced using crashing cost. In the first model, we consider the demand of lead time follows a normal distribution, and in the second model, it is considered distribution-free. For both cases, our objective is to minimize the integrated system cost by simultaneously optimizing the order quantity, safety factor, lead time and number of lots. The discounted cash flow and classical optimization technique are used to derive the optimal solution for both cases. Numerical examples including the sensitivity analysis of system parameters is provided to validate the results of the supply chain models.

Integral ceramic superstructure evaluation using time domain optical coherence tomography

Science.gov (United States)

Sinescu, Cosmin; Bradu, Adrian; Topala, Florin I.; Negrutiu, Meda Lavinia; Duma, Virgil-Florin; Podoleanu, Adrian G.

2014-02-01

Optical Coherence Tomography (OCT) is a non-invasive low coherence interferometry technique that includes several technologies (and the corresponding devices and components), such as illumination and detection, interferometry, scanning, adaptive optics, microscopy and endoscopy. From its large area of applications, we consider in this paper a critical aspect in dentistry - to be investigated with a Time Domain (TD) OCT system. The clinical situation of an edentulous mandible is considered; it can be solved by inserting 2 to 6 implants. On these implants a mesostructure will be manufactured and on it a superstructure is needed. This superstructure can be integral ceramic; in this case materials defects could be trapped inside the ceramic layers and those defects could lead to fractures of the entire superstructure. In this paper we demonstrate that a TD-OCT imaging system has the potential to properly evaluate the presence of the defects inside the ceramic layers and those defects can be fixed before inserting the prosthesis inside the oral cavity. Three integral ceramic superstructures were developed by using a CAD/CAM technology. After the milling, the ceramic layers were applied on the core. All the three samples were evaluated by a TD-OCT system working at 1300 nm. For two of the superstructures evaluated, no defects were found in the most stressed areas. The third superstructure presented four ceramic defects in the mentioned areas. Because of those defects the superstructure may fracture. The integral ceramic prosthesis was send back to the dental laboratory to fix the problems related to the material defects found. Thus, TD-OCT proved to be a valuable method for diagnosing the ceramic defects inside the integral ceramic superstructures in order to prevent fractures at this level.

Predictability of monthly temperature and precipitation using automatic time series forecasting methods

Science.gov (United States)

Papacharalampous, Georgia; Tyralis, Hristos; Koutsoyiannis, Demetris

2018-02-01

We investigate the predictability of monthly temperature and precipitation by applying automatic univariate time series forecasting methods to a sample of 985 40-year-long monthly temperature and 1552 40-year-long monthly precipitation time series. The methods include a naïve one based on the monthly values of the last year, as well as the random walk (with drift), AutoRegressive Fractionally Integrated Moving Average (ARFIMA), exponential smoothing state-space model with Box-Cox transformation, ARMA errors, Trend and Seasonal components (BATS), simple exponential smoothing, Theta and Prophet methods. Prophet is a recently introduced model inspired by the nature of time series forecasted at Facebook and has not been applied to hydrometeorological time series before, while the use of random walk, BATS, simple exponential smoothing and Theta is rare in hydrology. The methods are tested in performing multi-step ahead forecasts for the last 48 months of the data. We further investigate how different choices of handling the seasonality and non-normality affect the performance of the models. The results indicate that: (a) all the examined methods apart from the naïve and random walk ones are accurate enough to be used in long-term applications; (b) monthly temperature and precipitation can be forecasted to a level of accuracy which can barely be improved using other methods; (c) the externally applied classical seasonal decomposition results mostly in better forecasts compared to the automatic seasonal decomposition used by the BATS and Prophet methods; and (d) Prophet is competitive, especially when it is combined with externally applied classical seasonal decomposition.

Quadratic algebras in the noncommutative integration method of wave equation

International Nuclear Information System (INIS)

Varaksin, O.L.

1995-01-01

The paper deals with the investigation of applications of the method of noncommutative integration of linear differential equations by partial derivatives. Nontrivial example was taken for integration of three-dimensions wave equation with the use of non-Abelian quadratic algebras

FLRW cosmology in Weyl-integrable space-time

Energy Technology Data Exchange (ETDEWEB)

Gannouji, Radouane [Department of Physics, Faculty of Science, Tokyo University of Science, 1–3, Kagurazaka, Shinjuku-ku, Tokyo 162-8601 (Japan); Nandan, Hemwati [Department of Physics, Gurukula Kangri Vishwavidayalaya, Haridwar 249404 (India); Dadhich, Naresh, E-mail: gannouji@rs.kagu.tus.ac.jp, E-mail: hntheory@yahoo.co.in, E-mail: nkd@iucaa.ernet.in [IUCAA, Post Bag 4, Ganeshkhind, Pune 411 007 (India)

2011-11-01

We investigate the Weyl space-time extension of general relativity (GR) for studying the FLRW cosmology through focusing and defocusing of the geodesic congruences. We have derived the equations of evolution for expansion, shear and rotation in the Weyl space-time. In particular, we consider the Starobinsky modification, f(R) = R+βR{sup 2}−2Λ, of gravity in the Einstein-Palatini formalism, which turns out to reduce to the Weyl integrable space-time (WIST) with the Weyl vector being a gradient. The modified Raychaudhuri equation takes the form of the Hill-type equation which is then analysed to study the formation of the caustics. In this model, it is possible to have a Big Bang singularity free cyclic Universe but unfortunately the periodicity turns out to be extremely short.

Integral methods of solving boundary-value problems of nonstationary heat conduction and their comparative analysis

Science.gov (United States)

Kot, V. A.

2017-11-01

The modern state of approximate integral methods used in applications, where the processes of heat conduction and heat and mass transfer are of first importance, is considered. Integral methods have found a wide utility in different fields of knowledge: problems of heat conduction with different heat-exchange conditions, simulation of thermal protection, Stefantype problems, microwave heating of a substance, problems on a boundary layer, simulation of a fluid flow in a channel, thermal explosion, laser and plasma treatment of materials, simulation of the formation and melting of ice, inverse heat problems, temperature and thermal definition of nanoparticles and nanoliquids, and others. Moreover, polynomial solutions are of interest because the determination of a temperature (concentration) field is an intermediate stage in the mathematical description of any other process. The following main methods were investigated on the basis of the error norms: the Tsoi and Postol’nik methods, the method of integral relations, the Gudman integral method of heat balance, the improved Volkov integral method, the matched integral method, the modified Hristov method, the Mayer integral method, the Kudinov method of additional boundary conditions, the Fedorov boundary method, the method of weighted temperature function, the integral method of boundary characteristics. It was established that the two last-mentioned methods are characterized by high convergence and frequently give solutions whose accuracy is not worse that the accuracy of numerical solutions.

National Ignition Facility sub-system design requirements integrated timing system SSDR 1.5.3

International Nuclear Information System (INIS)

Wiedwald, J.; Van Aersau, P.; Bliss, E.

1996-01-01

This System Design Requirement document establishes the performance, design, development, and test requirements for the Integrated Timing System, WBS 1.5.3 which is part of the NIF Integrated Computer Control System (ICCS). The Integrated Timing System provides all temporally-critical hardware triggers to components and equipment in other NIF systems

Cochrane Qualitative and Implementation Methods Group guidance series-paper 5: methods for integrating qualitative and implementation evidence within intervention effectiveness reviews.

Science.gov (United States)

Harden, Angela; Thomas, James; Cargo, Margaret; Harris, Janet; Pantoja, Tomas; Flemming, Kate; Booth, Andrew; Garside, Ruth; Hannes, Karin; Noyes, Jane

2018-05-01

The Cochrane Qualitative and Implementation Methods Group develops and publishes guidance on the synthesis of qualitative and mixed-method evidence from process evaluations. Despite a proliferation of methods for the synthesis of qualitative research, less attention has focused on how to integrate these syntheses within intervention effectiveness reviews. In this article, we report updated guidance from the group on approaches, methods, and tools, which can be used to integrate the findings from quantitative studies evaluating intervention effectiveness with those from qualitative studies and process evaluations. We draw on conceptual analyses of mixed methods systematic review designs and the range of methods and tools that have been used in published reviews that have successfully integrated different types of evidence. We outline five key methods and tools as devices for integration which vary in terms of the levels at which integration takes place; the specialist skills and expertise required within the review team; and their appropriateness in the context of limited evidence. In situations where the requirement is the integration of qualitative and process evidence within intervention effectiveness reviews, we recommend the use of a sequential approach. Here, evidence from each tradition is synthesized separately using methods consistent with each tradition before integration takes place using a common framework. Reviews which integrate qualitative and process evaluation evidence alongside quantitative evidence on intervention effectiveness in a systematic way are rare. This guidance aims to support review teams to achieve integration and we encourage further development through reflection and formal testing. Copyright © 2017 Elsevier Inc. All rights reserved.

Integrating Expressive Methods in a Relational-Psychotherapy

Directory of Open Access Journals (Sweden)

Richard G. Erskine

2011-06-01

Full Text Available Therapeutic Involvement is an integral part of all effective psychotherapy.This article is written to illustrate the concept of Therapeutic Involvement in working within a therapeutic relationship – within the transference -- and with active expressive and experiential methods to resolve traumatic experiences, relational disturbances and life shaping decisions.

Integral method for the calculation of Hawking radiation in dispersive media. II. Asymmetric asymptotics.

Science.gov (United States)

Robertson, Scott