A multi scale approximation solution for the time dependent Boltzmann-transport equation

International Nuclear Information System (INIS)

Merk, B.

2004-03-01

The basis of all transient simulations for nuclear reactor cores is the reliable calculation of the power production. The local power distribution is generally calculated by solving the space, time, energy and angle dependent neutron transport equation known as Boltzmann equation. The computation of exact solutions of the Boltzmann equation is very time consuming. For practical numerical simulations approximated solutions are usually unavoidable. The objective of this work is development of an effective multi scale approximation solution for the Boltzmann equation. Most of the existing methods are based on separation of space and time. The new suggested method is performed without space-time separation. This effective approximation solution is developed on the basis of an expansion for the time derivative of different approximations to the Boltzmann equation. The method of multiple scale expansion is used for the expansion of the time derivative, because the problem of the stiff time behaviour can't be expressed by standard expansion methods. This multiple scale expansion is used in this work to develop approximation solutions for different approximations of the Boltzmann equation, starting from the expansion of the point kinetics equations. The resulting analytic functions are used for testing the applicability and accuracy of the multiple scale expansion method for an approximation solution with 2 delayed neutron groups. The results are tested versus the exact analytical results for the point kinetics equations. Very good agreement between both solutions is obtained. The validity of the solution with 2 delayed neutron groups to approximate the behaviour of the system with 6 delayed neutron groups is demonstrated in an additional analysis. A strategy for a solution with 4 delayed neutron groups is described. A multiple scale expansion is performed for the space-time dependent diffusion equation for one homogenized cell with 2 delayed neutron groups. The result is

Exact solutions to the supply chain equations for arbitrary, time-dependent demands

DEFF Research Database (Denmark)

Warburton, Roger D.H.; Hodgson, J.P.E.; Nielsen, Erland Hejn

2014-01-01

, so users can determine the inventory behavior to any desired precision. To illustrate, we solve the equations for a non-linear, quadratic time-dependence in the demand. For practical use, only a few terms in the series are required, a proposition illustrated by the For All Practical Purposes (FAPP......We study the impact on inventory of an unexpected, non-linear, time-dependent demand and present the exact solutions over time to the supply chain equations without requiring any approximations. We begin by imposing a boundary condition of stability at infinity, from which we derive expressions...... for the estimated demand and the target work in progress when the demand is time-dependent. The resulting inventory equation is solved in terms of the Lambert modes with all of the demand non-linearities confined to the pre-shape function. The series solution is exact, and all terms are reasonably easy to calculate...

Finite moments approach to the time-dependent neutron transport equation

International Nuclear Information System (INIS)

Kim, Sang Hyun

1994-02-01

Currently, nodal techniques are widely used in solving the multidimensional diffusion equation because of savings in computing time and storage. Thanks to the development of computer technology, one can now solve the transport equation instead of the diffusion equation to obtain more accurate solution. The finite moments method, one of the nodal methods, attempts to represent the fluxes in the cell and on cell surfaces more rigorously by retaining additional spatial moments. Generally, there are two finite moments schemes to solve the time-dependent transport equation. In one, the time variable is treated implicitly with finite moments method in space variable (implicit finite moments method), the other method uses finite moments method in both space and time (space-time finite moments method). In this study, these two schemes are applied to two types of time-dependent neutron transport problems. One is a fixed source problem, the other a heterogeneous fast reactor problem with delayed neutrons. From the results, it is observed that the two finite moments methods give almost the same solutions in both benchmark problems. However, the space-time finite moments method requires a little longer computing time than that of the implicit finite moments method. In order to reduce the longer computing time in the space-time finite moments method, a new iteration strategy is exploited, where a few time-stepwise calculation, in which original time steps are grouped into several coarse time divisions, is performed sequentially instead of performing iterations over the entire time steps. This strategy results in significant reduction of the computing time and we observe that 2-or 3-stepwise calculation is preferable. In addition, we propose a new finite moments method which is called mixed finite moments method in this thesis. Asymptotic analysis for the finite moments method shows that accuracy of the solution in a heterogeneous problem mainly depends on the accuracy of the

The large discretization step method for time-dependent partial differential equations

Science.gov (United States)

Haras, Zigo; Taasan, Shlomo

1995-01-01

A new method for the acceleration of linear and nonlinear time dependent calculations is presented. It is based on the Large Discretization Step (LDS) approximation, defined in this work, which employs an extended system of low accuracy schemes to approximate a high accuracy discrete approximation to a time dependent differential operator. Error bounds on such approximations are derived. These approximations are efficiently implemented in the LDS methods for linear and nonlinear hyperbolic equations, presented here. In these algorithms the high and low accuracy schemes are interpreted as the same discretization of a time dependent operator on fine and coarse grids, respectively. Thus, a system of correction terms and corresponding equations are derived and solved on the coarse grid to yield the fine grid accuracy. These terms are initialized by visiting the fine grid once in many coarse grid time steps. The resulting methods are very general, simple to implement and may be used to accelerate many existing time marching schemes.

Time-dependent field equations for paraxial relativistic electron beams: Beam Research Program

International Nuclear Information System (INIS)

Sharp, W.M.; Yu, S.S.; Lee, E.P.

1987-01-01

A simplified set of field equations for a paraxial relativistic electron beam is presented. These equations for the beam electrostatic potential phi and pinch potential Phi identical to A/sub z/ - phi retain previously neglected time-dependent terms and for axisymmetric beams reduce exactly to Maxwell's equations

Fortran programs for the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

Science.gov (United States)

Muruganandam, P.; Adhikari, S. K.

2009-10-01

Here we develop simple numerical algorithms for both stationary and non-stationary solutions of the time-dependent Gross-Pitaevskii (GP) equation describing the properties of Bose-Einstein condensates at ultra low temperatures. In particular, we consider algorithms involving real- and imaginary-time propagation based on a split-step Crank-Nicolson method. In a one-space-variable form of the GP equation we consider the one-dimensional, two-dimensional circularly-symmetric, and the three-dimensional spherically-symmetric harmonic-oscillator traps. In the two-space-variable form we consider the GP equation in two-dimensional anisotropic and three-dimensional axially-symmetric traps. The fully-anisotropic three-dimensional GP equation is also considered. Numerical results for the chemical potential and root-mean-square size of stationary states are reported using imaginary-time propagation programs for all the cases and compared with previously obtained results. Also presented are numerical results of non-stationary oscillation for different trap symmetries using real-time propagation programs. A set of convenient working codes developed in Fortran 77 are also provided for all these cases (twelve programs in all). In the case of two or three space variables, Fortran 90/95 versions provide some simplification over the Fortran 77 programs, and these programs are also included (six programs in all). Program summaryProgram title: (i) imagetime1d, (ii) imagetime2d, (iii) imagetime3d, (iv) imagetimecir, (v) imagetimesph, (vi) imagetimeaxial, (vii) realtime1d, (viii) realtime2d, (ix) realtime3d, (x) realtimecir, (xi) realtimesph, (xii) realtimeaxial Catalogue identifier: AEDU_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDU_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data

Spectral methods for time dependent partial differential equations

Science.gov (United States)

Gottlieb, D.; Turkel, E.

1983-01-01

The theory of spectral methods for time dependent partial differential equations is reviewed. When the domain is periodic Fourier methods are presented while for nonperiodic problems both Chebyshev and Legendre methods are discussed. The theory is presented for both hyperbolic and parabolic systems using both Galerkin and collocation procedures. While most of the review considers problems with constant coefficients the extension to nonlinear problems is also discussed. Some results for problems with shocks are presented.

Inverse scattering transform for the time dependent Schrödinger equation with applications to the KPI equation

Science.gov (United States)

Zhou, Xin

1990-03-01

For the direct-inverse scattering transform of the time dependent Schrödinger equation, rigorous results are obtained based on an opertor-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first order n x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution.

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.

Exact solutions of time-dependent Dirac equations and the quantum-classical correspondence

International Nuclear Information System (INIS)

Zhang Zhiguo

2006-01-01

Exact solutions to the Dirac equations with a time-dependent mass and a static magnetic field or a time-dependent linear potential are given. Matrix elements of the coordinate, momentum and velocity operator are calculated. In the large quantum number limit, these matrix elements give the classical solution

Energy decay of a variable-coefficient wave equation with nonlinear time-dependent localized damping

Directory of Open Access Journals (Sweden)

Jieqiong Wu

2015-09-01

Full Text Available We study the energy decay for the Cauchy problem of the wave equation with nonlinear time-dependent and space-dependent damping. The damping is localized in a bounded domain and near infinity, and the principal part of the wave equation has a variable-coefficient. We apply the multiplier method for variable-coefficient equations, and obtain an energy decay that depends on the property of the coefficient of the damping term.

Form-preserving Transformations for the Time-dependent Schroedinger Equation in (n + 1) Dimensions

International Nuclear Information System (INIS)

Schulze-Halberg, Axel

2006-01-01

We define a form-preserving transformation (also called point canonical transformation) for the time-dependent Schroedinger equation (TDSE) in (n+1) dimensions. The form-preserving transformation is shown to be invertible and to preserve L 2 -normalizability. We give a class of time-dependent TDSEs that can be mapped onto stationary Schroedinger equations by our form-preserving transformation. As an example, we generate a solvable, time-dependent potential of Coulombic ring-shaped type together with the corresponding exact solution of the TDSE in (3+1) dimensions. We further consider TDSEs with position-dependent (effective) masses and show that there is no form-preserving transformation between them and the conventional TDSEs, if the spatial dimension of the system is higher than one

Simulation of Time-Dependent P3 Equations Using a Semi-Analog Medium

International Nuclear Information System (INIS)

Hadad, K.; Pirouzmand, A.; Suh, Kune Y.

2010-01-01

A wide variety of numerical methods have been introduced to solve the neutron transport equation for reactor calculations. With the state-of-the-art computer technology, successful implementation of higher-order approximation of transport methods (P N , S N , MOC, etc.) may now be feasible. Although these methods have been adaptable to code parallelization techniques, the computational expense remains a significant obstacle and thwarts their implementation in a whole-core, time-dependent methodology. A novel method to remove this problem is based on the method of cellular neural networks (CNN) coupling with the PN method. Parallel data processing in CNN reduces the processing time and makes it possible to solve the time dependent models of neutron transport equation in real time

The master symmetry and time dependent symmetries of the differential–difference KP equation

International Nuclear Information System (INIS)

Khanizadeh, Farbod

2014-01-01

We first obtain the master symmetry of the differential–difference KP equation. Then we show how this master symmetry, through sl(2,C)-representation of the equation, can construct generators of time dependent symmetries. (paper)

Application of Trotter approximation for solving time dependent neutron transport equation

International Nuclear Information System (INIS)

Stancic, V.

1987-01-01

A method is proposed to solve multigroup time dependent neutron transport equation with arbitrary scattering anisotropy. The recurrence relation thus obtained is simple, numerically stable and especially suitable for treatment of complicated geometries. (author)

An elementary solution of the Maxwell equations for a time-dependent source

International Nuclear Information System (INIS)

Rivera, R; Villarroel, D

2002-01-01

We present an elementary solution of the Maxwell equations for a time-dependent source consisting of an infinite solenoid with a current density that increases linearly with time. The geometrical symmetries and the time dependence of the current density make possible a mathematical treatment that does not involve the usual technical difficulties, thus making this presentation suitable for students that are taking a first course in electromagnetism. We also show that the electric field generated by the solenoid can be used to construct an exact solution of the relativistic equation of motion of the electron that takes into account the effect of the radiation. In particular, we derive, in an almost trivial way, the formula for the radiation rate of an electron in circular motion

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

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

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

On time transformations for differential equations with state-dependent delay

Czech Academy of Sciences Publication Activity Database

Rezunenko, Oleksandr

2014-01-01

Roč. 12, č. 2 (2014), s. 298-307 ISSN 1895-1074 R&D Projects: GA ČR GAP103/12/2431 Institutional support: RVO:67985556 Keywords : differential equations * state-dependent delay * time transformations Subject RIV: BD - Theory of Information Impact factor: 0.578, year: 2014 http://library.utia.cas.cz/separaty/2014/AS/rezunenko-0429130.pdf

The time-dependent Hartree-Fock equations with Coulomb two-body interaction

International Nuclear Information System (INIS)

Chadam, J.M.; Glassey, R.T.

1975-06-01

The existence and uniqueness of global solutions to the Cauchy problem is proved in the space of ''smooth'' density matrices for the time-dependent Hartree-Fock equations describing the motion of finite Fermi systems interacting via a Coulomb two-body potential [fr

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)

Nonlinear Maxwell's and Schrodinger equations for describing the volumetric interaction of femtosecond laser pulses with transparent solid dielectrics: effect of the boundary conditions

Czech Academy of Sciences Publication Activity Database

Zhukov, V.P.; Bulgakova, Nadezhda M.; Fedoruk, M.P.

2017-01-01

Roč. 84, č. 7 (2017), s. 439-446 ISSN 1070-9762 R&D Projects: GA MŠk LO1602; GA ČR GA16-12960S Institutional support: RVO:68378271 Keywords : glass * femtosecond laser pulses * Maxwell's and Schrdinger equations Subject RIV: BH - Optics, Masers, Lasers OBOR OECD: Optics (including laser optics and quantum optics) Impact factor: 0.299, year: 2016

AIREK-MOD, Time Dependent Reactor Kinetics with Feedback Differential Equation

International Nuclear Information System (INIS)

Tamagnini, C.

1984-01-01

1 - Nature of physical problem solved: Solves the reactor kinetic equations with respect to time. A standard form for the reactivity behaviour has been introduced in which the reactivity is given by the sum of a polynomial, sine, cosine and exponential expansion. Tabular form is also included. The presence of feedback differential equations in which the dependence on variables different from the considered one is considered enables many heat-exchange problems to be dealt with. 2 - Method of solution: The method employed for the solution of the differential equations is the one developed by E.R. Cohen (Geneva Conference, 1958). 3 - Restrictions on the complexity of the problem: The maximum number of differential equations that can be solved simultaneously is 50. Within this limitation there may be n delayed neutron groups (n less than or equal to 25), on m other linear feedback equations (n+m less than or equal to 49). CDC 1604 version was offered by EIR (Institut Federal de Recherches en matiere de reacteurs, Switzerland)

1-Soliton solution of the generalized Zakharov-Kuznetsov equation with nonlinear dispersion and time-dependent coefficients

International Nuclear Information System (INIS)

Biswas, Anjan

2009-01-01

In this Letter, the 1-soliton solution of the Zakharov-Kuznetsov equation with power law nonlinearity and nonlinear dispersion along with time-dependent coefficients is obtained. There are two models for this kind of an equation that are studied. The constraint relation between these time-dependent coefficients is established for the solitons to exist. Subsequently, this equation is again analysed with generalized evolution. The solitary wave ansatz is used to carry out this investigation.

Density-Dependent Conformable Space-time Fractional Diffusion-Reaction Equation and Its Exact Solutions

Science.gov (United States)

Hosseini, Kamyar; Mayeli, Peyman; Bekir, Ahmet; Guner, Ozkan

2018-01-01

In this article, a special type of fractional differential equations (FDEs) named the density-dependent conformable fractional diffusion-reaction (DDCFDR) equation is studied. Aforementioned equation has a significant role in the modelling of some phenomena arising in the applied science. The well-organized methods, including the \\exp (-φ (\\varepsilon )) -expansion and modified Kudryashov methods are exerted to generate the exact solutions of this equation such that some of the solutions are new and have been reported for the first time. Results illustrate that both methods have a great performance in handling the DDCFDR equation.

Time dependent solutions of the Fokker-Planck equation for fast fusion ions

International Nuclear Information System (INIS)

Gnavi, G.; Gratton, F.T.; Heyn, M.

1990-01-01

Approximate time dependent solutions for the Fokker-Planck equation for fast fusion ions from an isotropic, monoenergetic source are presented, for the problem of D - T - He 3 reactions. The equations include the effect of diffusion, which is particularly noticeable in the distribution of particles of lower energy and in the formation of a tail of particles with energy higher than that of the source. (Author)

Estimation of time- and state-dependent delays and other parameters in functional differential equations

Science.gov (United States)

Murphy, K. A.

1990-01-01

A parameter estimation algorithm is developed which can be used to estimate unknown time- or state-dependent delays and other parameters (e.g., initial condition) appearing within a nonlinear nonautonomous functional differential equation. The original infinite dimensional differential equation is approximated using linear splines, which are allowed to move with the variable delay. The variable delays are approximated using linear splines as well. The approximation scheme produces a system of ordinary differential equations with nice computational properties. The unknown parameters are estimated within the approximating systems by minimizing a least-squares fit-to-data criterion. Convergence theorems are proved for time-dependent delays and state-dependent delays within two classes, which say essentially that fitting the data by using approximations will, in the limit, provide a fit to the data using the original system. Numerical test examples are presented which illustrate the method for all types of delay.

Time-dependent Hartree approximation and time-dependent harmonic oscillator model

International Nuclear Information System (INIS)

Blaizot, J.P.

1982-01-01

We present an analytically soluble model for studying nuclear collective motion within the framework of the time-dependent Hartree (TDH) approximation. The model reduces the TDH equations to the Schroedinger equation of a time-dependent harmonic oscillator. Using canonical transformations and coherent states we derive a few properties of the time-dependent harmonic oscillator which are relevant for applications. We analyse the role of the normal modes in the time evolution of a system governed by TDH equations. We show how these modes couple together due to the anharmonic terms generated by the non-linearity of the theory. (orig.)

Trial function method and exact solutions to the generalized nonlinear Schrödinger equation with time-dependent coefficient

International Nuclear Information System (INIS)

Cao Rui; Zhang Jian

2013-01-01

In this paper, the trial function method is extended to study the generalized nonlinear Schrödinger equation with time-dependent coefficients. On the basis of a generalized traveling wave transformation and a trial function, we investigate the exact envelope traveling wave solutions of the generalized nonlinear Schrödinger equation with time-dependent coefficients. Taking advantage of solutions to trial function, we successfully obtain exact solutions for the generalized nonlinear Schrödinger equation with time-dependent coefficients under constraint conditions. (general)

Time-dependent simplified PN approximation to the equations of radiative transfer

International Nuclear Information System (INIS)

Frank, Martin; Klar, Axel; Larsen, Edward W.; Yasuda, Shugo

2007-01-01

The steady-state simplified P N approximation to the radiative transport equation has been successfully applied to many problems involving radiation. This paper presents the derivation of time-dependent simplified P N (SP N ) equations (up to N = 3) via two different approaches. First, we use an asymptotic analysis, similar to the asymptotic derivation of the steady-state SP N equations. Second, we use an approach similar to the original derivation of the steady-state SP N equations and we show that both approaches lead to similar results. Special focus is put on the well-posedness of the equations and the question whether it can be guaranteed that the solution satisfies the correct physical bounds. Several numerical test cases are shown, including an analytical benchmark due to Su and Olson [B. Su, G.L. Olson, An analytical benchmark for non-equilibrium radiative transfer in an isotropically scattering medium, Ann. Nucl. Energy 24 (1997) 1035-1055.

Solution of the time-dependent, three-dimensional resistive magnetohydrodynamic equations

International Nuclear Information System (INIS)

Finan, C.H. III; Killeen, J.; California Univ., Davis

1981-01-01

Resistive magnetohydrodynamics (MHD) is described by a set of eight coupled, nonlinear, three-dimensional, time-dependent, partial differential equations. A computer code, IMP (Implicit MHD Program), has been developed to solve these equations numerically by the method of finite differences on an Eulerian mesh. In this model, the equations are expressed in orthogonal curvilinear coordinates, making the code applicable to a variety of coordinate systems. The Douglas-Gunn algorithm for Alternating-Direction Implicit (ADI) temporal advancement is used to avoid the limitations in timestep size imposed by explicit methods. The equations are expressed as conservation laws, the momentum and energy equations are nonconservative. This is to: (1) provide enhanced numerical stability by eliminating errors introduced by the nonvanishing of nabla x B on the finite difference mesh; and, (2) allow the simulation of low β plasmas. The resulting finite difference equations are a coupled system of nonlinear algebraic equations which are solved by the Newton-Raphson iteration technique. We apply our model to a number of problems of importance in magnetic fusion research. Ideal and resistive internal kink instabilities are simulated in a Cartesian geometry. Growth rates and nonlinear saturation amplitudes are found to be in agreement with previous analytic and numerical predictions. We also simulate these instabilities in a torus, which demonstrates the versatility of the orthogonal curvilinear coordinate representation. (orig.)

Similarity solutions of reaction–diffusion equation with space- and time-dependent diffusion and reaction terms

Energy Technology Data Exchange (ETDEWEB)

Ho, C.-L. [Department of Physics, Tamkang University, Tamsui 25137, Taiwan (China); Lee, C.-C., E-mail: chieh.no27@gmail.com [Center of General Education, Aletheia University, Tamsui 25103, Taiwan (China)

2016-01-15

We consider solvability of the generalized reaction–diffusion equation with both space- and time-dependent diffusion and reaction terms by means of the similarity method. By introducing the similarity variable, the reaction–diffusion equation is reduced to an ordinary differential equation. Matching the resulting ordinary differential equation with known exactly solvable equations, one can obtain corresponding exactly solvable reaction–diffusion systems. Several representative examples of exactly solvable reaction–diffusion equations are presented.

Initial-boundary value problems for multi-term time-fractional diffusion equations with x-dependent coefficients

OpenAIRE

Li, Zhiyuan; Huang, Xinchi; Yamamoto, Masahiro

2018-01-01

In this paper, we discuss an initial-boundary value problem (IBVP) for the multi-term time-fractional diffusion equation with x-dependent coefficients. By means of the Mittag-Leffler functions and the eigenfunction expansion, we reduce the IBVP to an equivalent integral equation to show the unique existence and the analyticity of the solution for the equation. Especially, in the case where all the coefficients of the time-fractional derivatives are non-negative, by the Laplace and inversion L...

On combined optical solitons of the one-dimensional Schrödinger’s equation with time dependent coefficients

Directory of Open Access Journals (Sweden)

Kilic Bulent

2016-01-01

Full Text Available This paper integrates dispersive optical solitons in special optical metamaterials with a time dependent coefficient. We obtained some optical solitons of the aforementioned equation. It is shown that the examined dependent coefficients are affected by the velocity of the wave. The first integral method (FIM and ansatz method are applied to reach the optical soliton solutions of the one-dimensional nonlinear Schrödinger’s equation (NLSE with time dependent coefficients.

Time Dependent Hartree Fock Equation: Gateway to Nonequilibrium Plasmas

International Nuclear Information System (INIS)

Dufty, James W.

2007-01-01

This is the Final Technical Report for DE-FG02-2ER54677 award 'Time Dependent Hartree Fock Equation - Gateway to Nonequilibrium Plasmas'. Research has focused on the nonequilibrium dynamics of electrons in the presence of ions, both via basic quantum theory and via semi-classical molecular dynamics (MD) simulation. In addition, fundamental notions of dissipative dynamics have been explored for models of grains and dust, and for scalar fields (temperature) in turbulent edge plasmas. The specific topics addressed were Quantum Kinetic Theory for Metallic Clusters, Semi-classical MD Simulation of Plasmas , and Effects of Dissipative Dynamics.

A Derivation of Source-based Kinetics Equation with Time Dependent Fission Kernel for Reactor Transient Analyses

International Nuclear Information System (INIS)

Kim, Song Hyun; Woo, Myeong Hyun; Shin, Chang Ho; Pyeon, Cheol Ho

2015-01-01

In this study, a new balance equation to overcome the problems generated by the previous methods is proposed using source-based balance equation. And then, a simple problem is analyzed with the proposed method. In this study, a source-based balance equation with the time dependent fission kernel was derived to simplify the kinetics equation. To analyze the partial variations of reactor characteristics, two representative methods were introduced in previous studies; (1) quasi-statics method and (2) multipoint technique. The main idea of quasistatics method is to use a low-order approximation for large integration times. To realize the quasi-statics method, first, time dependent flux is separated into the shape and amplitude functions, and shape function is calculated. It is noted that the method has a good accuracy; however, it can be expensive as a calculation cost aspect because the shape function should be fully recalculated to obtain accurate results. To improve the calculation efficiency, multipoint method was proposed. The multipoint method is based on the classic kinetics equation with using Green's function to analyze the flight probability from region r' to r. Those previous methods have been used to analyze the reactor kinetics analysis; however, the previous methods can have some limitations. First, three group variables (r g , E g , t g ) should be considered to solve the time dependent balance equation. This leads a big limitation to apply large system problem with good accuracy. Second, the energy group neutrons should be used to analyze reactor kinetics problems. In time dependent problem, neutron energy distribution can be changed at different time. It can affect the change of the group cross section; therefore, it can lead the accuracy problem. Third, the neutrons in a space-time region continually affect the other space-time regions; however, it is not properly considered in the previous method. Using birth history of the neutron sources

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.

Time-dependent embedding

OpenAIRE

Inglesfield, J. E.

2007-01-01

A method of solving the time-dependent Schr\\"odinger equation is presented, in which a finite region of space is treated explicitly, with the boundary conditions for matching the wave-functions on to the rest of the system replaced by an embedding term added on to the Hamiltonian. This time-dependent embedding term is derived from the Fourier transform of the energy-dependent embedding potential, which embeds the time-independent Schr\\"odinger equation. Results are presented for a one-dimensi...

Numerical method for solving the three-dimensional time-dependent neutron diffusion equation

International Nuclear Information System (INIS)

Khaled, S.M.; Szatmary, Z.

2005-01-01

A numerical time-implicit method has been developed for solving the coupled three-dimensional time-dependent multi-group neutron diffusion and delayed neutron precursor equations. The numerical stability of the implicit computation scheme and the convergence of the iterative associated processes have been evaluated. The computational scheme requires the solution of large linear systems at each time step. For this purpose, the point over-relaxation Gauss-Seidel method was chosen. A new scheme was introduced instead of the usual source iteration scheme. (author)

Stability of the Filter Equation for a Time-Dependent Signal on Rd

International Nuclear Information System (INIS)

Stannat, Wilhelm

2005-01-01

Stability of the pathwise filter equation for a time-dependent signal process induced by a d-dimensional stochastic differential equation and a linear observation is studied, using a variational approach. A lower bound for the rate of stability is identified in terms of the mass-gap of a parabolic ground state transform associated with the generator of the signal process and the square of the observation. The lower bound can be easily calculated a priori and provides hints on how precisely to measure the signal in order to reach a certain rate of stability. Ergodicity of the signal process is not needed

A Derivation of Source-based Kinetics Equation with Time Dependent Fission Kernel for Reactor Transient Analyses

Energy Technology Data Exchange (ETDEWEB)

Kim, Song Hyun; Woo, Myeong Hyun; Shin, Chang Ho [Hanyang University, Seoul (Korea, Republic of); Pyeon, Cheol Ho [Kyoto University, Osaka (Japan)

2015-10-15

In this study, a new balance equation to overcome the problems generated by the previous methods is proposed using source-based balance equation. And then, a simple problem is analyzed with the proposed method. In this study, a source-based balance equation with the time dependent fission kernel was derived to simplify the kinetics equation. To analyze the partial variations of reactor characteristics, two representative methods were introduced in previous studies; (1) quasi-statics method and (2) multipoint technique. The main idea of quasistatics method is to use a low-order approximation for large integration times. To realize the quasi-statics method, first, time dependent flux is separated into the shape and amplitude functions, and shape function is calculated. It is noted that the method has a good accuracy; however, it can be expensive as a calculation cost aspect because the shape function should be fully recalculated to obtain accurate results. To improve the calculation efficiency, multipoint method was proposed. The multipoint method is based on the classic kinetics equation with using Green's function to analyze the flight probability from region r' to r. Those previous methods have been used to analyze the reactor kinetics analysis; however, the previous methods can have some limitations. First, three group variables (r{sub g}, E{sub g}, t{sub g}) should be considered to solve the time dependent balance equation. This leads a big limitation to apply large system problem with good accuracy. Second, the energy group neutrons should be used to analyze reactor kinetics problems. In time dependent problem, neutron energy distribution can be changed at different time. It can affect the change of the group cross section; therefore, it can lead the accuracy problem. Third, the neutrons in a space-time region continually affect the other space-time regions; however, it is not properly considered in the previous method. Using birth history of the

Solution of the time-dependent inertial-frame equation of radiative transfer in moving media to O(v/c)

International Nuclear Information System (INIS)

Mihalas, D.; Klein, R.I.

1982-01-01

A stable and efficient mixed-frame method has been formulated for the solution of the time-dependent equation of radiative transfer with full retention of all velocity dependent terms to O(ν/c). The method retains the simplicity of the differential operator found in the inertial frame while transforming the absorption and emission coefficients to the comoving frame keeping them isotropic. The method is ideally suited to continuum calculations. To correctly treat the time dependence of the radiation field over fluid-flow time increments, the velocity-dependent terms on the right-hand side of both the transfer and moment equations must be retained for consistency

Pullback-Forward Dynamics for Damped Schrödinger Equations with Time-Dependent Forcing

Directory of Open Access Journals (Sweden)

Lianbing She

2018-01-01

Full Text Available This paper deals with pullback dynamics for the weakly damped Schrödinger equation with time-dependent forcing. An increasing, bounded, and pullback absorbing set is obtained if the forcing and its time-derivative are backward uniformly integrable. Also, we obtain the forward absorption, which is only used to deduce the backward compact-decay decomposition according to high and low frequencies. Based on a new existence theorem of a backward compact pullback attractor, we show that the nonautonomous Schrödinger equation has a pullback attractor which is compact in the past. The method of energy, high-low frequency decomposition, Sobolev embedding, and interpolation are quite involved in calculating a priori pullback or forward bound.

Asymptotic behaviors of solutions for viscoelastic wave equation with space-time dependent damping term

KAUST Repository

Said-Houari, Belkacem

2012-03-01

In this paper, we consider a viscoelastic wave equation with an absorbing term and space-time dependent damping term. Based on the weighted energy method, and by assuming that the kernel decaying exponentially, we obtain the L2 decay rates of the solutions. More precisely, we show that the decay rates are the same as those obtained in Lin et al. (2010) [15] for the semilinear wave equation with absorption term. © 2011 Elsevier Inc.

Asymptotic behaviors of solutions for viscoelastic wave equation with space-time dependent damping term

KAUST Repository

Said-Houari, Belkacem

2012-01-01

In this paper, we consider a viscoelastic wave equation with an absorbing term and space-time dependent damping term. Based on the weighted energy method, and by assuming that the kernel decaying exponentially, we obtain the L2 decay rates of the solutions. More precisely, we show that the decay rates are the same as those obtained in Lin et al. (2010) [15] for the semilinear wave equation with absorption term. © 2011 Elsevier Inc.

Time-dependent integral transport equation kernels, leakage rates and collision rates for plane and spherical geometry

International Nuclear Information System (INIS)

Henderson, D.L.

1987-01-01

Time-dependent integral transport equation flux and current kernels for plane and spherical geometry are derived for homogeneous media. Using the multiple collision formalism, isotropic sources that are delta distributions in time are considered for four different problems. The plane geometry flux kernel is applied to a uniformly distributed source within an infinite medium and to a surface source in a semi-infinite medium. The spherical flux kernel is applied to a point source in an infinite medium and to a point source at the origin of a finite sphere. The time-dependent first-flight leakage rates corresponding to the existing steady state first-flight escape probabilities are computed by the Laplace transform technique assuming a delta distribution source in time. The case of a constant source emitting neutrons over a time interval, Δt, for a spatially uniform source is obtained for a slab and a sphere. Time-dependent first-flight leakage rates are also determined for the general two region spherical medium problem for isotropic sources with a delta distribution in time uniformly distributed throughout both the inner and outer regions. The time-dependent collision rates due to the uncollided neutrons are computed for a slab and a sphere using the time-dependent first-flight leakage rates and the time-dependent continuity equation. The case of a constant source emitting neutrons over a time interval, Δt, is also considered

Application of the multigrid amplitude function method for time-dependent transport equation using MOC

International Nuclear Information System (INIS)

Tsujita, K.; Endo, T.; Yamamoto, A.

2013-01-01

An efficient numerical method for time-dependent transport equation, the mutigrid amplitude function (MAF) method, is proposed. The method of characteristics (MOC) is being widely used for reactor analysis thanks to the advances of numerical algorithms and computer hardware. However, efficient kinetic calculation method for MOC is still desirable since it requires significant computation time. Various efficient numerical methods for solving the space-dependent kinetic equation, e.g., the improved quasi-static (IQS) and the frequency transform methods, have been developed so far mainly for diffusion calculation. These calculation methods are known as effective numerical methods and they offer a way for faster computation. However, they have not been applied to the kinetic calculation method using MOC as the authors' knowledge. Thus, the MAF method is applied to the kinetic calculation using MOC aiming to reduce computation time. The MAF method is a unified numerical framework of conventional kinetic calculation methods, e.g., the IQS, the frequency transform, and the theta methods. Although the MAF method is originally developed for the space-dependent kinetic calculation based on the diffusion theory, it is extended to transport theory in the present study. The accuracy and computational time are evaluated though the TWIGL benchmark problem. The calculation results show the effectiveness of the MAF method. (authors)

Space-time Dependency of the Time and its Effect on the Relativistic Classical Equation of the String Theory

Science.gov (United States)

Gholibeigian, Hassan; Amirshahkarami, Abdolazim; Gholibeigian, Kazem

2017-01-01

In special relativity theory, time dilates in velocity of near light speed. Also based on ``Substantial motion'' theory of Sadra, relative time (time flux); R = f (mv , σ , τ) , for each atom is momentum of its involved fundamental particles, which is different from the other atoms. In this way, for modification of the relativistic classical equation of string theory and getting more precise results, we should use effect of dilation and contraction of time in equation. So we propose to add two derivatives of the time's flux to the equation as follows: n.tp∂/R ∂ τ +∂2Xμ/(σ , τ) ∂τ2 = n .tp (∂/R ∂ σ ) +c2∂2Xμ/(σ , τ) ∂σ2 In which, Xμ is space-time coordinates of the string, σ & τ are coordinates on the string world sheet, respectively space and time along the string, string's mass m , velocity of string's motion v , factor n depends on geometry of each hidden extra dimension which relates to its own flux time, and tp is Planck's time. AmirKabir University of Technology, Tehran, Iran.

Solution to the monoenergetic time-dependent neutron transport equation with a time-varying source

International Nuclear Information System (INIS)

Ganapol, B.D.

1986-01-01

Even though fundamental time-dependent neutron transport problems have existed since the inception of neutron transport theory, it has only been recently that a reliable numerical solution to one of the basic problems has been obtained. Experience in generating numerical solutions to time-dependent transport equations has indicated that the multiple collision formulation is the most versatile numerical technique for model problems. The formulation coupled with a moment reconstruction of each collided flux component has led to benchmark-quality (four- to five-digit accuracy) numerical evaluation of the neutron flux in plane infinite geometry for any degree of scattering anisotropy and for both pulsed isotropic and beam sources. As will be shown in this presentation, this solution can serve as a Green's function, thus extending the previous results to more complicated source situations. Here we will be concerned with a time-varying source at the center of an infinite medium. If accurate, such solutions have both pedagogical and practical uses as benchmarks against which other more approximate solutions designed for a wider class of problems can be compared

Self-consistent predictor/corrector algorithms for stable and efficient integration of the time-dependent Kohn-Sham equation

Science.gov (United States)

Zhu, Ying; Herbert, John M.

2018-01-01

The "real time" formulation of time-dependent density functional theory (TDDFT) involves integration of the time-dependent Kohn-Sham (TDKS) equation in order to describe the time evolution of the electron density following a perturbation. This approach, which is complementary to the more traditional linear-response formulation of TDDFT, is more efficient for computation of broad-band spectra (including core-excited states) and for systems where the density of states is large. Integration of the TDKS equation is complicated by the time-dependent nature of the effective Hamiltonian, and we introduce several predictor/corrector algorithms to propagate the density matrix, one of which can be viewed as a self-consistent extension of the widely used modified-midpoint algorithm. The predictor/corrector algorithms facilitate larger time steps and are shown to be more efficient despite requiring more than one Fock build per time step, and furthermore can be used to detect a divergent simulation on-the-fly, which can then be halted or else the time step modified.

Propagators for the Time-Dependent Kohn-Sham Equations: Multistep, Runge-Kutta, Exponential Runge-Kutta, and Commutator Free Magnus Methods.

Science.gov (United States)

Gómez Pueyo, Adrián; Marques, Miguel A L; Rubio, Angel; Castro, Alberto

2018-05-09

We examine various integration schemes for the time-dependent Kohn-Sham equations. Contrary to the time-dependent Schrödinger's equation, this set of equations is nonlinear, due to the dependence of the Hamiltonian on the electronic density. We discuss some of their exact properties, and in particular their symplectic structure. Four different families of propagators are considered, specifically the linear multistep, Runge-Kutta, exponential Runge-Kutta, and the commutator-free Magnus schemes. These have been chosen because they have been largely ignored in the past for time-dependent electronic structure calculations. The performance is analyzed in terms of cost-versus-accuracy. The clear winner, in terms of robustness, simplicity, and efficiency is a simplified version of a fourth-order commutator-free Magnus integrator. However, in some specific cases, other propagators, such as some implicit versions of the multistep methods, may be useful.

Angular distribution of scission neutrons studied with time-dependent Schrödinger equation

Science.gov (United States)

Wada, Takahiro; Asano, Tomomasa; Carjan, Nicolae

2018-03-01

We investigate the angular distribution of scission neutrons taking account of the effects of fission fragments. The time evolution of the wave function of the scission neutron is obtained by integrating the time-dependent Schrodinger equation numerically. The effects of the fission fragments are taken into account by means of the optical potentials. The angular distribution is strongly modified by the presence of the fragments. In the case of asymmetric fission, it is found that the heavy fragment has stronger effects. Dependence on the initial distribution and on the properties of fission fragments is discussed. We also discuss on the treatment of the boundary to avoid artificial reflections

Wave function for time-dependent harmonically confined electrons in a time-dependent electric field.

Science.gov (United States)

Li, Yu-Qi; Pan, Xiao-Yin; Sahni, Viraht

2013-09-21

The many-body wave function of a system of interacting particles confined by a time-dependent harmonic potential and perturbed by a time-dependent spatially homogeneous electric field is derived via the Feynman path-integral method. The wave function is comprised of a phase factor times the solution to the unperturbed time-dependent Schrödinger equation with the latter being translated by a time-dependent value that satisfies the classical driven equation of motion. The wave function reduces to that of the Harmonic Potential Theorem wave function for the case of the time-independent harmonic confining potential.

Time-dependent Schroedinger equations with effective mass in (2 + 1) dimensions: intertwining relations and Darboux operators

International Nuclear Information System (INIS)

Cobian, Hector; Schulze-Halberg, Axel

2011-01-01

We construct Darboux transformations for time-dependent Schroedinger equations with position-dependent mass in (2 + 1) dimensions. Several examples illustrate our results, which complement and generalize former findings for the constant mass case in two spatial variables (Schulze-Halberg 2010 J. Math. Phys. 51 033521).

Linear Scaling Solution of the Time-Dependent Self-Consistent-Field Equations

Directory of Open Access Journals (Sweden)

Matt Challacombe

2014-03-01

Full Text Available A new approach to solving the Time-Dependent Self-Consistent-Field equations is developed based on the double quotient formulation of Tsiper 2001 (J. Phys. B. Dual channel, quasi-independent non-linear optimization of these quotients is found to yield convergence rates approaching those of the best case (single channel Tamm-Dancoff approximation. This formulation is variational with respect to matrix truncation, admitting linear scaling solution of the matrix-eigenvalue problem, which is demonstrated for bulk excitons in the polyphenylene vinylene oligomer and the (4,3 carbon nanotube segment.

Determination of parameters for a stress-strain constitutive equation considering time-dependent behavior of Toki granite

International Nuclear Information System (INIS)

Hirano, Toru; Seno, Yasuhiro; Nakama, Shigeo; Okubo, Seisuke

2008-01-01

Toki granite was tested to obtain parameters for the constitutive equation. The testing method was uniaxial compressive loading at the moderate a constant strain rate that is decreased after yielding to obtain the complete stress-strain curve. In addition, two kinds of the strain rate were alternately switched to obtain the parameter n from one specimen. The n represents the strength time-dependence in the constitutive equation. The second parameter m can be obtained by fitting the experimental stress-strain curve to the calculated curve. The m accounts for the behavior after yielding. According to the results, Toki granite has n=52 and m=60, showing relatively weak time-dependence of creep failure. (author)

Quantum trajectories for time-dependent adiabatic master equations

Science.gov (United States)

Yip, Ka Wa; Albash, Tameem; Lidar, Daniel A.

2018-02-01

We describe a quantum trajectories technique for the unraveling of the quantum adiabatic master equation in Lindblad form. By evolving a complex state vector of dimension N instead of a complex density matrix of dimension N2, simulations of larger system sizes become feasible. The cost of running many trajectories, which is required to recover the master equation evolution, can be minimized by running the trajectories in parallel, making this method suitable for high performance computing clusters. In general, the trajectories method can provide up to a factor N advantage over directly solving the master equation. In special cases where only the expectation values of certain observables are desired, an advantage of up to a factor N2 is possible. We test the method by demonstrating agreement with direct solution of the quantum adiabatic master equation for 8-qubit quantum annealing examples. We also apply the quantum trajectories method to a 16-qubit example originally introduced to demonstrate the role of tunneling in quantum annealing, which is significantly more time consuming to solve directly using the master equation. The quantum trajectories method provides insight into individual quantum jump trajectories and their statistics, thus shedding light on open system quantum adiabatic evolution beyond the master equation.

Solution of time dependent atmospheric diffusion equation with a proposed diffusion coefficient

International Nuclear Information System (INIS)

Mayhoub, A.B.; Essa, KH.S.M.; Aly, SH.

2004-01-01

One-dimensional model for the dispersion of passive atmospheric contaminant (not included chemical reactions) in the atmospheric boundary layer is considered. On the basis of the gradient transfer theory (K-theory), the time dependent diffusion equation represents the dispersion of the pollutants is solved analytically. The solution depends on diffusion coefficient K', which is expressed in terms of the friction velocity 'u the vertical coordinate -L and the depth of the mixing layer 'h'. The solution is obtained to either the vertical coordinate 'z' is less or greater than the mixing height 'h'. The obtained solution may be applied to study the atmospheric dispersion of pollutants

A Solution of Time Dependent Schrodinger Equation by Quantum Walk

International Nuclear Information System (INIS)

Sekino, Hideo; Kawahata, Masayuki; Hamada, Shinji

2012-01-01

Time Dependent Schroedinger Equation (TDSE) with an initial Gaussian distribution, is solved by a discrete time/space Quantum Walk (QW) representing consecutive operations corresponding to a dot product of Pauli matrix and momentum operators. We call it as Schroedinger Walk (SW). Though an Hadamard Walk (HW) provides same dynamics of the probability distribution for delta-function-like initial distributions as that of the SW with a delta-function-like initial distribution, the former with a Gaussian initial distribution leads to a solution for advection of the probability distribution; the initial distribution splits into two distinctive distributions moving in opposite directions. Both mechanisms are analysed by investigating the evolution of the both amplitude components. Decoherence of the oscillating amplitudes in central region is found to be responsible for the splitting of the probability distribution in the HW.

Delay differential equations and the dose-time dependence of early radiotherapy reactions

International Nuclear Information System (INIS)

Fenwick, John D.

2006-01-01

The dose-time dependence of early radiotherapy reactions impacts on the design of accelerated fractionation schedules--oral mucositis, for example, can be dose limiting for short treatments designed to avoid tumor repopulation. In this paper a framework for modeling early reaction dose-time dependence is developed. Variation of stem cell number with time after the start of a radiation schedule is modeled using a first-order delay differential equation (DDE), motivated by experimental observations linking the speed of compensatory proliferation in early reacting tissues to the degree of tissue damage. The modeling suggests that two types of early reaction radiation response are possible, stem cell numbers either monotonically approaching equilibrium plateau levels or overshooting before returning to equilibrium. Several formulas have been derived from the delay differential equation, predicting changes in isoeffective total radiation dose with schedule duration for different types of fractionation scheme. The formulas have been fitted to a wide range of published animal early reaction data, the fits all implying a degree of overshoot. Results are presented illustrating the scope of the delay differential model: most of the data are fitted well, although the model struggles with a few datasets measured for schedules with distinctive dose-time patterns. Ways of extending the current model to cope with these particular dose-time patterns are briefly discussed. The DDE approach is conceptually more complex than earlier descriptive dose-time models but potentially more powerful. It can be used to study issues not addressed by simpler models, such as the likely effects of increasing or decreasing the dose-per-day over time, or of splitting radiation courses into intense segments separated by gaps. It may also prove useful for modeling the effects of chemoirradiation

Delay differential equations and the dose-time dependence of early radiotherapy reactions.

Science.gov (United States)

Fenwick, John D

2006-09-01

The dose-time dependence of early radiotherapy reactions impacts on the design of accelerated fractionation schedules--oral mucositis, for example, can be dose limiting for short treatments designed to avoid tumor repopulation. In this paper a framework for modeling early reaction dose-time dependence is developed. Variation of stem cell number with time after the start of a radiation schedule is modeled using a first-order delay differential equation (DDE), motivated by experimental observations linking the speed of compensatory proliferation in early reacting tissues to the degree of tissue damage. The modeling suggests that two types of early reaction radiation response are possible, stem cell numbers either monotonically approaching equilibrium plateau levels or overshooting before returning to equilibrium. Several formulas have been derived from the delay differential equation, predicting changes in isoeffective total radiation dose with schedule duration for different types of fractionation scheme. The formulas have been fitted to a wide range of published animal early reaction data, the fits all implying a degree of overshoot. Results are presented illustrating the scope of the delay differential model: most of the data are fitted well, although the model struggles with a few datasets measured for schedules with distinctive dose-time patterns. Ways of extending the current model to cope with these particular dose-time patterns are briefly discussed. The DDE approach is conceptually more complex than earlier descriptive dose-time models but potentially more powerful. It can be used to study issues not addressed by simpler models, such as the likely effects of increasing or decreasing the dose-per-day over time, or of splitting radiation courses into intense segments separated by gaps. It may also prove useful for modeling the effects of chemoirradiation.

Introducing Relativity into Quantum Chemistry

Science.gov (United States)

Li, Wai-Kee; Blinder, S. M.

2011-01-01

It is not often realized by chemists that the special theory of relativity is behind several aspects of quantum chemistry. The Schrdinger equation itself is based on relations between space-time and energy-momentum four vectors. Electron spin is, of course, the most obvious manifestation of relativity. The chemistry of some heavy elements is…

Alternating-direction implicit numerical solution of the time-dependent, three-dimensional, single fluid, resistive magnetohydrodynamic equations

Energy Technology Data Exchange (ETDEWEB)

Finan, C.H. III

1980-12-01

Resistive magnetohydrodynamics (MHD) is described by a set of eight coupled, nonlinear, three-dimensional, time-dependent, partial differential equations. A computer code, IMP (Implicit MHD Program), has been developed to solve these equations numerically by the method of finite differences on an Eulerian mesh. In this model, the equations are expressed in orthogonal curvilinear coordinates, making the code applicable to a variety of coordinate systems. The Douglas-Gunn algorithm for Alternating-Direction Implicit (ADI) temporal advancement is used to avoid the limitations in timestep size imposed by explicit methods. The equations are solved simultaneously to avoid syncronization errors.

Darboux transformations for the time-dependent nonhomogeneous Burgers equation in (1+1) dimensions

International Nuclear Information System (INIS)

Schulze-Halberg, Axel; Manuel Carballo Jimenez, Juan

2009-01-01

We extend the formalism of nth order Darboux transformations to the time-dependent nonhomogeneous Burgers equation (NBE) in (1+1) dimensions. Similar to the Schroedinger case, our Darboux transformation retains the form of the NBE, while changing the nonhomogeneous term. The transformed solution of the NBE and the corresponding transformed nonhomogeneity are given in closed form. Furthermore, properties of the transformation are discussed and an application is given.

Time dependent resonating Hartree-Bogoliubov theory

International Nuclear Information System (INIS)

Nishiyama, Seiya; Fukutome, Hideo.

1989-01-01

Very recently, we have developed a theory of excitations in superconducting Fermion systems with large quantum fluctuations that can be described by resonance of time dependent non-orthogonal Hartree-Bogoliubov (HB) wave functions with different correlation structures. We have derived a new kind of variation equation called the time dependent Resonating HB equation, in order to determine both the time dependent Resonating HB wave functions and coefficients of a superposition of the HB wave functions. Further we have got a new approximation for excitations from time dependent small fluctuations of the Resonating HB ground state, i.e., the Resonating HB RPA. The Res HB RPA equation is represented in a given single particle basis. It, however, has drawbacks that the constraints for the Res HB RPA amplitudes are not taken into account and the equation contains equations which are not independent. We shall derive another form of the Res HB RPA equation eliminating these drawbacks. The Res HB RPA gives a unified description of the vibrons and resonons and their interactions. (author)

Uniform in N global well-posedness of the time-dependent Hartree-Fock-Bogoliubov equations in R^{1+1}

Science.gov (United States)

Chong, Jacky Jia Wei

2018-04-01

We prove the global well-posedness of the time-dependent Hartree-Fock-Bogoliubov (TDHFB) equations in R^{1+1} with two-body interaction potential of the form N^{-1}v_N(x) = N^{β -1} v(N^β x) where v≥0 is a sufficiently regular radial function, i.e., v \\in L^1(R)\\cap C^∞ (R) . In particular, using methods of dispersive PDEs similar to the ones used in Grillakis and Machedon (Commun Partial Differ Equ 42:24-67, 2017), we are able to show for any scaling parameter β >0 the TDHFB equations are globally well-posed in some Strichartz-type spaces independent of N, cf. (Bach et al. in The time-dependent Hartree-Fock-Bogoliubov equations for Bosons, 2016. arXiv:1602.05171).

Global existence of solutions to the Cauchy problem for time-dependent Hartree equations

International Nuclear Information System (INIS)

Chadam, J.M.; Glassey, R.T.

1975-01-01

The existence of global solutions to the Cauchy problem for time-dependent Hartree equations for N electrons is established. The solution is shown to have a uniformly bounded H 1 (R 3 ) norm and to satisfy an estimate of the form two parallel PSI (t) two parallel/sub H 2 ; less than or equal to c exp(kt). It is shown that ''negative energy'' solutions do not converge uniformly to zero as t → infinity. (U.S.)

Dispersion Decay and Scattering Theory

CERN Document Server

Komech, Alexander

2012-01-01

A simplified, yet rigorous treatment of scattering theory methods and their applications Dispersion Decay and Scattering Theory provides thorough, easy-to-understand guidance on the application of scattering theory methods to modern problems in mathematics, quantum physics, and mathematical physics. Introducing spectral methods with applications to dispersion time-decay and scattering theory, this book presents, for the first time, the Agmon-Jensen-Kato spectral theory for the Schr?dinger equation, extending the theory to the Klein-Gordon equation. The dispersion decay plays a crucial role i

Time-dependent Gross-Pitaevskii equation for composite bosons as the strong-coupling limit of the fermionic broken-symmetry random-phase approximation

International Nuclear Information System (INIS)

Strinati, G.C.; Pieri, P.

2004-01-01

The linear response to a space- and time-dependent external disturbance of a system of dilute condensed composite bosons at zero temperature, as obtained from the linearized version of the time-dependent Gross-Pitaevskii equation, is shown to result also from the strong-coupling limit of the time-dependent BCS (or broken-symmetry random-phase) approximation for the constituent fermions subject to the same external disturbance. In this way, it is possible to connect excited-state properties of the bosonic and fermionic systems by placing the Gross-Pitaevskii equation in perspective with the corresponding fermionic approximations

Rhodium SPND's Error Reduction using Extended Kalman Filter combined with Time Dependent Neutron Diffusion Equation

International Nuclear Information System (INIS)

Lee, Jeong Hun; Park, Tong Kyu; Jeon, Seong Su

2014-01-01

The Rhodium SPND is accurate in steady-state conditions but responds slowly to changes in neutron flux. The slow response time of Rhodium SPND precludes its direct use for control and protection purposes specially when nuclear power plant is used for load following. To shorten the response time of Rhodium SPND, there were some acceleration methods but they could not reflect neutron flux distribution in reactor core. On the other hands, some methods for core power distribution monitoring could not consider the slow response time of Rhodium SPND and noise effect. In this paper, time dependent neutron diffusion equation is directly used to estimate reactor power distribution and extended Kalman filter method is used to correct neutron flux with Rhodium SPND's and to shorten the response time of them. Extended Kalman filter is effective tool to reduce measurement error of Rhodium SPND's and even simple FDM to solve time dependent neutron diffusion equation can be an effective measure. This method reduces random errors of detectors and can follow reactor power level without cross-section change. It means monitoring system may not calculate cross-section at every time steps and computing time will be shorten. To minimize delay of Rhodium SPND's conversion function h should be evaluated in next study. Neutron and Rh-103 reaction has several decay chains and half-lives over 40 seconds causing delay of detection. Time dependent neutron diffusion equation will be combined with decay chains. Power level and distribution change corresponding movement of control rod will be tested with more complicated reference code as well as xenon effect. With these efforts, final result is expected to be used as a powerful monitoring tool of nuclear reactor core

Rhodium SPND's Error Reduction using Extended Kalman Filter combined with Time Dependent Neutron Diffusion Equation

Energy Technology Data Exchange (ETDEWEB)

Lee, Jeong Hun; Park, Tong Kyu; Jeon, Seong Su [FNC Technology Co., Ltd., Yongin (Korea, Republic of)

2014-05-15

The Rhodium SPND is accurate in steady-state conditions but responds slowly to changes in neutron flux. The slow response time of Rhodium SPND precludes its direct use for control and protection purposes specially when nuclear power plant is used for load following. To shorten the response time of Rhodium SPND, there were some acceleration methods but they could not reflect neutron flux distribution in reactor core. On the other hands, some methods for core power distribution monitoring could not consider the slow response time of Rhodium SPND and noise effect. In this paper, time dependent neutron diffusion equation is directly used to estimate reactor power distribution and extended Kalman filter method is used to correct neutron flux with Rhodium SPND's and to shorten the response time of them. Extended Kalman filter is effective tool to reduce measurement error of Rhodium SPND's and even simple FDM to solve time dependent neutron diffusion equation can be an effective measure. This method reduces random errors of detectors and can follow reactor power level without cross-section change. It means monitoring system may not calculate cross-section at every time steps and computing time will be shorten. To minimize delay of Rhodium SPND's conversion function h should be evaluated in next study. Neutron and Rh-103 reaction has several decay chains and half-lives over 40 seconds causing delay of detection. Time dependent neutron diffusion equation will be combined with decay chains. Power level and distribution change corresponding movement of control rod will be tested with more complicated reference code as well as xenon effect. With these efforts, final result is expected to be used as a powerful monitoring tool of nuclear reactor core.

An implicit fast Fourier transform method for integration of the time dependent Schrodinger or diffusion equation

International Nuclear Information System (INIS)

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

1997-06-01

The authors have found that the conventional exponentiated split operator procedure is subject to difficulties in energy conservation when solving the time-dependent Schrodinger equation for Coulombic systems. By rearranging the kinetic and potential energy terms in the temporal propagator of the finite difference equations, one can find a propagation algorithm for three dimensions that looks much like the Crank-Nicholson and alternating direction implicit methods for one- and two-space-dimensional partial differential equations. They report comparisons of this novel implicit split operator procedure with the conventional exponentiated split operator procedure on hydrogen atom solutions. The results look promising for a purely numerical approach to certain electron quantum mechanical problems

State-dependent differential Riccati equation to track control of time-varying systems with state and control nonlinearities.

Science.gov (United States)

Korayem, M H; Nekoo, S R

2015-07-01

This work studies an optimal control problem using the state-dependent Riccati equation (SDRE) in differential form to track for time-varying systems with state and control nonlinearities. The trajectory tracking structure provides two nonlinear differential equations: the state-dependent differential Riccati equation (SDDRE) and the feed-forward differential equation. The independence of the governing equations and stability of the controller are proven along the trajectory using the Lyapunov approach. Backward integration (BI) is capable of solving the equations as a numerical solution; however, the forward solution methods require the closed-form solution to fulfill the task. A closed-form solution is introduced for SDDRE, but the feed-forward differential equation has not yet been obtained. Different ways of solving the problem are expressed and analyzed. These include BI, closed-form solution with corrective assumption, approximate solution, and forward integration. Application of the tracking problem is investigated to control robotic manipulators possessing rigid or flexible joints. The intention is to release a general program for automatic implementation of an SDDRE controller for any manipulator that obeys the Denavit-Hartenberg (D-H) principle when only D-H parameters are received as input data. Copyright © 2015 ISA. Published by Elsevier Ltd. All rights reserved.

A stable computational scheme for stiff time-dependent constitutive equations

International Nuclear Information System (INIS)

Shih, C.F.; Delorenzi, H.G.; Miller, A.K.

1977-01-01

Viscoplasticity and creep type constitutive equations are increasingly being employed in finite element codes for evaluating the deformation of high temperature structural members. These constitutive equations frequently exhibit stiff regimes which makes an analytical assessment of the structure very costly. A computational scheme for handling deformation in stiff regimes is proposed in this paper. By the finite element discretization, the governing partial differential equations in the spatial (x) and time (t) variables are reduced to a system of nonlinear ordinary differential equations in the independent variable t. The constitutive equations are expanded in a Taylor's series about selected values of t. The resulting system of differential equations are then integrated by an implicit scheme which employs a predictor technique to initiate the Newton-Raphson procedure. To examine the stability and accuracy of the computational scheme, a series of calculations were carried out for uniaxial specimens and thick wall tubes subjected to mechanical and thermal loading. (Auth.)

Spectral fitting method for the solution of time-dependent Schroedinger equations: Applications to atoms in intense laser fields

International Nuclear Information System (INIS)

Qiao Haoxue; Cai Qingyu; Rao Jianguo; Li Baiwen

2002-01-01

A spectral fitting method for solving the time-dependent Schroedinger equation has been developed and applied to the atom in intense laser fields. This method allows us to obtain a highly accurate time-dependent wave function with a contribution from the high-order term of Δt. Moreover, the time-dependent wave function is determined on a small number of discrete mesh points, thus making calculations simple and accurate. This method is illustrated by computing wave functions and harmonic generation spectra of a model atom in laser fields

From Feshbach-resonance managed Bose-Einstein condensates to anisotropic universes: Applications of the Ermakov-Pinney equation with time-dependent nonlinearity

International Nuclear Information System (INIS)

Herring, G.; Kevrekidis, P.G.; Williams, F.; Christodoulakis, T.; Frantzeskakis, D.J.

2008-01-01

In this work we revisit the topic of two-dimensional Bose-Einstein condensates under the influence of time-dependent magnetic confinement and time-dependent scattering length. A moment approach reduces the examination of moments of the wavefunction (in particular, of its width) to an Ermakov-Pinney (EP) ordinary differential equation (ODE). We use the well-known structure of the solutions of this nonlinear ODE to 'engineer' trapping and interatomic interaction conditions that lead to condensates dispersing, breathing or even collapsing. The advantage of the approach is that it is fully tractable analytically, in excellent agreement with our numerical observations. As an aside, we also discuss how similar time-dependent EP equations may arise in the description of anisotropic scalar field cosmologies

From Feshbach-resonance managed Bose-Einstein condensates to anisotropic universes: Applications of the Ermakov-Pinney equation with time-dependent nonlinearity

International Nuclear Information System (INIS)

Herring, G.; Kevrekidis, P.G.; Williams, F.; Christodoulakis, T.; Frantzeskakis, D.J.

2007-01-01

In this work we revisit the topic of two-dimensional Bose-Einstein condensates under the influence of time-dependent magnetic confinement and time-dependent scattering length. A moment approach reduces the examination of moments of the wavefunction (in particular, of its width) to an Ermakov-Pinney (EP) ordinary differential equation (ODE). We use the well-known structure of the solutions of this nonlinear ODE to 'engineer' trapping and interatomic interaction conditions that lead to condensates dispersing, breathing or even collapsing. The advantage of the approach is that it is fully tractable analytically, in excellent agreement with our numerical observations. As an aside, we also discuss how similar time-dependent EP equations may arise in the description of anisotropic scalar field cosmologies

LUCKY-TD code for solving the time-dependent transport equation with the use of parallel computations

Energy Technology Data Exchange (ETDEWEB)

Moryakov, A. V., E-mail: sailor@orc.ru [National Research Centre Kurchatov Institute (Russian Federation)

2016-12-15

An algorithm for solving the time-dependent transport equation in the P{sub m}S{sub n} group approximation with the use of parallel computations is presented. The algorithm is implemented in the LUCKY-TD code for supercomputers employing the MPI standard for the data exchange between parallel processes.

Derivation and solution of a time-dependent, nonlinear, Schrodinger-like equation for the superconductivity order parameter

International Nuclear Information System (INIS)

Esrick, M.A.

1981-01-01

A time-dependent, nonlinear, Schrodinger-like equation for the superconductivity order parameter is derived from the Gor'kov equations. Three types of traveling wave solutions of the equation are discussed. The phases and amplitudes of these solutions propagate at different speeds. The first type of solution has an amplitude that propagates as a soliton and it is suggested that this solution might correspond to the recently observed propagating collective modes of the order parameter. The amplitude of the second type of solution propagates as a periodic disturbance in space and time. It is suggested that this type of solution might explain the recently observed multiple values of the superconductor energy gap as well as the spatially inhomogenous superconducting state. The third type of solution, which is of a more general character, might provide some insight into non-periodic, inhomogeneous states occuring in superconductors. It is also proposed that quasiparticle injection and microwave irradiation might generate soliton-like disturbances in superconductors

Exact Time-Dependent Wave Functions of a Confined Time-Dependent Harmonic Oscillator with Two Moving Boundaries

International Nuclear Information System (INIS)

Lo, C.F.

2009-01-01

By applying the standard analytical techniques of solving partial differential equations, we have obtained the exact solution in terms of the Fourier sine series to the time-dependent Schroedinger equation describing a quantum one-dimensional harmonic oscillator of time-dependent frequency confined in an infinite square well with the two walls moving along some parametric trajectories. Based upon the orthonormal basis of quasi-stationary wave functions, the exact propagator of the system has also been analytically derived. Special cases like (i) a confined free particle, (ii) a confined time-independent harmonic oscillator, and (iii) an aging oscillator are examined, and the corresponding time-dependent wave functions are explicitly determined. Besides, the approach has been extended to solve the case of a confined generalized time-dependent harmonic oscillator for some parametric moving boundaries as well. (general)

Time-dependent massless Dirac fermions in graphene

Energy Technology Data Exchange (ETDEWEB)

Khantoul, Boubakeur, E-mail: bobphys@gmail.com [Department of Mathematics, City University London, Northampton Square, London EC1V 0HB (United Kingdom); Department of Physics, University of Jijel, BP 98, Ouled Aissa, 18000 Jijel (Algeria); Fring, Andreas, E-mail: a.fring@city.ac.uk [Department of Mathematics, City University London, Northampton Square, London EC1V 0HB (United Kingdom)

2015-10-30

Using the Lewis–Riesenfeld method of invariants we construct explicit analytical solutions for the massless Dirac equation in 2+1 dimensions describing quasi-particles in graphene. The Hamiltonian of the system considered contains some explicit time-dependence in addition to one resulting from being minimally coupled to a time-dependent vector potential. The eigenvalue equations for the two spinor components of the Lewis–Riesenfeld invariant are found to decouple into a pair of supersymmetric invariants in a similar fashion as the known decoupling for the time-independent Dirac Hamiltonians. - Highlights: • An explicit analytical solution for a massless 2+1 dimensional time-dependent Dirac equation is found. • All steps of the Lewis–Riesenfeld method have been carried out.

Time-dependent magnetization of a type-II superconductor numerically calculated by using the flux-creep equation

International Nuclear Information System (INIS)

Lee, J. H.; Park, I. S.; Ahmad, D.; Kim, D.; Kim, Y. C.; Ko, R. K.; Jeong, D. Y.

2012-01-01

The macroscopic magnetic behaviors of a type-II superconductor, such as the field- or the temperature-dependent magnetization, have been described by using critical state models. However, because the models are time-independent, the magnetic relaxation in a type-II superconductor cannot be described by them, and the time dependence of the magnetization can affect the field or the temperature-dependent magnetization curve described by the models. In order to avoid the time independence of critical state models, we try the numerical calculation used by Qin et al., who mainly calculated the temperature dependence of the ac susceptibility χ(T). Their calculation showed that the frequency-dependent χ(T) could be obtained by using the flux-creep equation. We calculated the field-dependent magnetization and magnetic relaxation by using a numerical method. The calculated field-dependent magnetization M(H) curves shows the shapes of a typical type-II superconductor. The calculated magnetic relaxation do not show a logarithmic decay of the magnetization, but the addition of a surface barrier to the relaxation calculation caused a clear logarithmic decay of the magnetization, producing a crossover at a mid-time. This means that the logarithmic magnetic relaxation is caused by not only flux creep but also a combination of flux creep and a surface barrier.

Numerical solution of the time dependent neutron transport equation by the method of the characteristics

International Nuclear Information System (INIS)

Talamo, Alberto

2013-01-01

This study presents three numerical algorithms to solve the time dependent neutron transport equation by the method of the characteristics. The algorithms have been developed taking into account delayed neutrons and they have been implemented into the novel MCART code, which solves the neutron transport equation for two-dimensional geometry and an arbitrary number of energy groups. The MCART code uses regular mesh for the representation of the spatial domain, it models up-scattering, and takes advantage of OPENMP and OPENGL algorithms for parallel computing and plotting, respectively. The code has been benchmarked with the multiplication factor results of a Boiling Water Reactor, with the analytical results for a prompt jump transient in an infinite medium, and with PARTISN and TDTORT results for cross section and source transients. The numerical simulations have shown that only two numerical algorithms are stable for small time steps

Numerical solution of the time dependent neutron transport equation by the method of the characteristics

Energy Technology Data Exchange (ETDEWEB)

Talamo, Alberto, E-mail: alby@anl.gov [Nuclear Engineering Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439 (United States)

2013-05-01

This study presents three numerical algorithms to solve the time dependent neutron transport equation by the method of the characteristics. The algorithms have been developed taking into account delayed neutrons and they have been implemented into the novel MCART code, which solves the neutron transport equation for two-dimensional geometry and an arbitrary number of energy groups. The MCART code uses regular mesh for the representation of the spatial domain, it models up-scattering, and takes advantage of OPENMP and OPENGL algorithms for parallel computing and plotting, respectively. The code has been benchmarked with the multiplication factor results of a Boiling Water Reactor, with the analytical results for a prompt jump transient in an infinite medium, and with PARTISN and TDTORT results for cross section and source transients. The numerical simulations have shown that only two numerical algorithms are stable for small time steps.

Simulating transient dynamics of the time-dependent time fractional Fokker-Planck systems

Science.gov (United States)

Kang, Yan-Mei

2016-09-01

For a physically realistic type of time-dependent time fractional Fokker-Planck (FP) equation, derived as the continuous limit of the continuous time random walk with time-modulated Boltzmann jumping weight, a semi-analytic iteration scheme based on the truncated (generalized) Fourier series is presented to simulate the resultant transient dynamics when the external time modulation is a piece-wise constant signal. At first, the iteration scheme is demonstrated with a simple time-dependent time fractional FP equation on finite interval with two absorbing boundaries, and then it is generalized to the more general time-dependent Smoluchowski-type time fractional Fokker-Planck equation. The numerical examples verify the efficiency and accuracy of the iteration method, and some novel dynamical phenomena including polarized motion orientations and periodic response death are discussed.

Solitary wave dynamics in time-dependent potentials

International Nuclear Information System (INIS)

Abou Salem, Walid K.

2008-01-01

The long time dynamics of solitary wave solutions of the nonlinear Schroedinger equation in time-dependent external potentials is rigorously studied. To set the stage, the well-posedness of the Cauchy problem for a generalized nonautonomous nonlinear Schroedinger equation with time-dependent nonlinearities and potential is established. Afterward, the dynamics of NLS solitary waves in time-dependent potentials is studied. It is shown that in the space-adiabatic regime where the external potential varies slowly in space compared to the size of the soliton, the dynamics of the center of the soliton is described by Hamilton's equations, plus terms due to radiation damping. Finally, two physical applications are discussed: the first is adiabatic transportation of solitons and the second is the Mathieu instability of trapped solitons due to time-periodic perturbations

A Generalized Estimating Equations Approach to Model Heterogeneity and Time Dependence in Capture-Recapture Studies

Directory of Open Access Journals (Sweden)

Akanda Md. Abdus Salam

2017-03-01

Full Text Available Individual heterogeneity in capture probabilities and time dependence are fundamentally important for estimating the closed animal population parameters in capture-recapture studies. A generalized estimating equations (GEE approach accounts for linear correlation among capture-recapture occasions, and individual heterogeneity in capture probabilities in a closed population capture-recapture individual heterogeneity and time variation model. The estimated capture probabilities are used to estimate animal population parameters. Two real data sets are used for illustrative purposes. A simulation study is carried out to assess the performance of the GEE estimator. A Quasi-Likelihood Information Criterion (QIC is applied for the selection of the best fitting model. This approach performs well when the estimated population parameters depend on the individual heterogeneity and the nature of linear correlation among capture-recapture occasions.

Simulating transient dynamics of the time-dependent time fractional Fokker–Planck systems

Energy Technology Data Exchange (ETDEWEB)

Kang, Yan-Mei, E-mail: ymkang@mail.xjtu.edu.cn

2016-09-16

For a physically realistic type of time-dependent time fractional Fokker–Planck (FP) equation, derived as the continuous limit of the continuous time random walk with time-modulated Boltzmann jumping weight, a semi-analytic iteration scheme based on the truncated (generalized) Fourier series is presented to simulate the resultant transient dynamics when the external time modulation is a piece-wise constant signal. At first, the iteration scheme is demonstrated with a simple time-dependent time fractional FP equation on finite interval with two absorbing boundaries, and then it is generalized to the more general time-dependent Smoluchowski-type time fractional Fokker–Planck equation. The numerical examples verify the efficiency and accuracy of the iteration method, and some novel dynamical phenomena including polarized motion orientations and periodic response death are discussed. - Highlights: • An iteration method is proposed for the transient dynamics of time-dependent time fractional Fokker–Planck equations. • The method is based on Fourier Series solution and the multi-step transition probability formula. • With the time-modulated subdiffusion on finite interval as example, the polarized motion orientation is disclosed. • With the time-modulated subdiffusion within a confined potential as example, the death of dynamic response is observed.

Simulating transient dynamics of the time-dependent time fractional Fokker–Planck systems

International Nuclear Information System (INIS)

Kang, Yan-Mei

2016-01-01

For a physically realistic type of time-dependent time fractional Fokker–Planck (FP) equation, derived as the continuous limit of the continuous time random walk with time-modulated Boltzmann jumping weight, a semi-analytic iteration scheme based on the truncated (generalized) Fourier series is presented to simulate the resultant transient dynamics when the external time modulation is a piece-wise constant signal. At first, the iteration scheme is demonstrated with a simple time-dependent time fractional FP equation on finite interval with two absorbing boundaries, and then it is generalized to the more general time-dependent Smoluchowski-type time fractional Fokker–Planck equation. The numerical examples verify the efficiency and accuracy of the iteration method, and some novel dynamical phenomena including polarized motion orientations and periodic response death are discussed. - Highlights: • An iteration method is proposed for the transient dynamics of time-dependent time fractional Fokker–Planck equations. • The method is based on Fourier Series solution and the multi-step transition probability formula. • With the time-modulated subdiffusion on finite interval as example, the polarized motion orientation is disclosed. • With the time-modulated subdiffusion within a confined potential as example, the death of dynamic response is observed.

Use of a theoretical equation of state to interpret time-dependent free volume in polymer glasses

International Nuclear Information System (INIS)

Curro, J.G.; Lagasse, R.R.; Simha, R.

1981-01-01

Many physical properties of polymer glasses change spontaneously during isothermal aging by a process commonly modeled as collapse of free volume. The model has not been verified rigorously because free volume cannot be unambiguously measured. In the present investigation we tentatively identify the free-volume fraction with the fraction of empty sites in the equation of state of Simha and Somcynsky. With this theory, volume recovery measurements can be analyzed to yield directly the time-dependent, free-volume fraction. Using this approach, recent volume measurements on poly(methyl methacrylate) are analyzed. The resulting free-volume fractions are then used in the Doolittle equation to predict the shift in stress relaxation curves at 23 0 C. These predicted shift factors agree with the experimental measurements of Cizmecioglu et al. In addition, it is shown that previous assumptions concerning temperature dependence of free volume are inconsistent with the theory

Solution of the atmospheric diffusion equation with a realistic diffusion coefficient and time dependent mixing height

International Nuclear Information System (INIS)

Mayhoub, A.B.; Etman, S.M.

1997-01-01

One dimensional model for the dispersion of a passive atmospheric contaminant (neglecting chemical reactions) in the atmospheric boundary layer is introduced. The differential equation representing the dispersion of pollutants is solved on the basis of gradient-transfer theory (K- theory). The present approach deals with a more appropriate and realistic profile for the diffusion coefficient K, which is expressed in terms of the friction velocity U, the vertical coordinate z and the depth of the mixing layer h, which is taken time dependent. After some mathematical simplification, the equation analytic obtained solution can be easily applied to case study concerning atmospheric dispersion of pollutants

Optimized effective potential in real time: Problems and prospects in time-dependent density-functional theory

International Nuclear Information System (INIS)

Mundt, Michael; Kuemmel, Stephan

2006-01-01

The integral equation for the time-dependent optimized effective potential (TDOEP) in time-dependent density-functional theory is transformed into a set of partial-differential equations. These equations only involve occupied Kohn-Sham orbitals and orbital shifts resulting from the difference between the exchange-correlation potential and the orbital-dependent potential. Due to the success of an analog scheme in the static case, a scheme that propagates orbitals and orbital shifts in real time is a natural candidate for an exact solution of the TDOEP equation. We investigate the numerical stability of such a scheme. An approximation beyond the Krieger-Li-Iafrate approximation for the time-dependent exchange-correlation potential is analyzed

The semi-classical limit of the time dependent Hartree-Fock equation. II. The Wick symbol of the solution

OpenAIRE

Amour, Laurent; Khodja, Mohamed; Nourrigat, Jean

2011-01-01

We study the Wick symbol of a solution of the time dependent Hartree Fock equation, under weaker hypotheses than those needed for the Weyl symbol in the first paper with thesame title. With similar, we prove some kind of Ehrenfest theorem for observables that are not pseudo-differential operators.

Exact norm-conserving stochastic time-dependent Hartree-Fock

International Nuclear Information System (INIS)

Tessieri, Luca; Wilkie, Joshua; Cetinbas, Murat

2005-01-01

We derive an exact single-body decomposition of the time-dependent Schroedinger equation for N pairwise interacting fermions. Each fermion obeys a stochastic time-dependent norm-preserving wave equation. As a first test of the method, we calculate the low energy spectrum of helium. An extension of the method to bosons is outlined

Soliton solutions of some nonlinear evolution equations with time ...

Indian Academy of Sciences (India)

Abstract. In this paper, we obtain exact soliton solutions of the modified KdV equation, inho- mogeneous nonlinear Schrödinger equation and G(m, n) equation with variable coefficients using solitary wave ansatz. The constraint conditions among the time-dependent coefficients turn out as necessary conditions for the ...

Separation of massive field equation of arbitrary spin in Robertson-Walker space-time

International Nuclear Information System (INIS)

Zecca, A.

2006-01-01

The massive spin-(3/2) field equation is explicitly integrated in the Robertson-Walker space-time by the Newman Penrose formalism. The solution is obtained by extending a separation procedure previously used to solve the spin-1 equation. The separated time dependence results in two coupled equations depending on the cosmological background evolution. The separated angular equations are explicitly integrated and the eigenvalues determined. The separated radial equations are integrated in the flat space-time case. The separation method of solution is then generalized, by induction, to prove the main result, that is the separability of the massive field equations of arbitrary spin in the Robertson-Walker space-time

Exact solution of a quantum forced time-dependent harmonic oscillator

Science.gov (United States)

Yeon, Kyu Hwang; George, Thomas F.; Um, Chung IN

1992-01-01

The Schrodinger equation is used to exactly evaluate the propagator, wave function, energy expectation values, uncertainty values, and coherent state for a harmonic oscillator with a time dependent frequency and an external driving time dependent force. These quantities represent the solution of the classical equation of motion for the time dependent harmonic oscillator.

Finite difference solution of the time dependent neutron group diffusion equations

International Nuclear Information System (INIS)

Hendricks, J.S.; Henry, A.F.

1975-08-01

In this thesis two unrelated topics of reactor physics are examined: the prompt jump approximation and alternating direction checkerboard methods. In the prompt jump approximation it is assumed that the prompt and delayed neutrons in a nuclear reactor may be described mathematically as being instantaneously in equilibrium with each other. This approximation is applied to the spatially dependent neutron diffusion theory reactor kinetics model. Alternating direction checkerboard methods are a family of finite difference alternating direction methods which may be used to solve the multigroup, multidimension, time-dependent neutron diffusion equations. The reactor mesh grid is not swept line by line or point by point as in implicit or explicit alternating direction methods; instead, the reactor mesh grid may be thought of as a checkerboard in which all the ''red squares'' and '' black squares'' are treated successively. Two members of this family of methods, the ADC and NSADC methods, are at least as good as other alternating direction methods. It has been found that the accuracy of implicit and explicit alternating direction methods can be greatly improved by the application of an exponential transformation. This transformation is incompatible with checkerboard methods. Therefore, a new formulation of the exponential transformation has been developed which is compatible with checkerboard methods and at least as good as the former transformation for other alternating direction methods

State-dependent neutral delay equations from population dynamics.

Science.gov (United States)

Barbarossa, M V; Hadeler, K P; Kuttler, C

2014-10-01

A novel class of state-dependent delay equations is derived from the balance laws of age-structured population dynamics, assuming that birth rates and death rates, as functions of age, are piece-wise constant and that the length of the juvenile phase depends on the total adult population size. The resulting class of equations includes also neutral delay equations. All these equations are very different from the standard delay equations with state-dependent delay since the balance laws require non-linear correction factors. These equations can be written as systems for two variables consisting of an ordinary differential equation (ODE) and a generalized shift, a form suitable for numerical calculations. It is shown that the neutral equation (and the corresponding ODE--shift system) is a limiting case of a system of two standard delay equations.

On an nth-order infinitesimal generator and time-dependent operator differential equation with a strongly almost periodic solution

Directory of Open Access Journals (Sweden)

Aribindi Satyanarayan Rao

2002-01-01

Full Text Available In a Banach space, if u is a Stepanov almost periodic solution of a certain nth-order infinitesimal generator and time-dependent operator differential equation with a Stepanov almost periodic forcing function, then u,u′,…,u (n−2 are all strongly almost periodic and u (n−1 is weakly almost periodic.

On the time evolution operator for time-dependent quadratic Hamiltonians

International Nuclear Information System (INIS)

Fernandez, F.M.

1989-01-01

The Schroedinger equation with a time-dependent quadratic Hamiltonian is investigated. The time-evolution operator is written as a product of exponential operators determined by the Heisenberg equations of motion. This product operator is shown to be global in the occupation number representation when the Hamiltonian is Hermitian. The success of some physical applications of the product-form representation is explained

A quasi-Lagrangian finite element method for the Navier-Stokes equations in a time-dependent domain

Science.gov (United States)

Lozovskiy, Alexander; Olshanskii, Maxim A.; Vassilevski, Yuri V.

2018-05-01

The paper develops a finite element method for the Navier-Stokes equations of incompressible viscous fluid in a time-dependent domain. The method builds on a quasi-Lagrangian formulation of the problem. The paper provides stability and convergence analysis of the fully discrete (finite-difference in time and finite-element in space) method. The analysis does not assume any CFL time-step restriction, it rather needs mild conditions of the form $\\Delta t\\le C$, where $C$ depends only on problem data, and $h^{2m_u+2}\\le c\\,\\Delta t$, $m_u$ is polynomial degree of velocity finite element space. Both conditions result from a numerical treatment of practically important non-homogeneous boundary conditions. The theoretically predicted convergence rate is confirmed by a set of numerical experiments. Further we apply the method to simulate a flow in a simplified model of the left ventricle of a human heart, where the ventricle wall dynamics is reconstructed from a sequence of contrast enhanced Computed Tomography images.

Generation of matter-wave solitons of the Gross-Pitaevskii equation with a time-dependent complicated potential

Energy Technology Data Exchange (ETDEWEB)

Mohamadou, Alidou [Condensed Matter Laboratory, Department of Physics, Faculty of Science, University of Douala, P.O. Box 24157, Douala (Cameroon); Abdus Salam International Centre for Theoretical Physics, P.O. Box 538, Strada Costiera 11, I-34014 Trieste (Italy); Wamba, Etienne; Kofane, Timoleon C. [Laboratory of Mechanics, Department of Physics, Faculty of Science, University of Yaounde I, P.O. Box 812, Yaounde (Cameroon); Doka, Serge Y. [Higher Teacher Training College, University of Maroua, P.O. Box 55, Maroua (Cameroon); Ekogo, Thierry B. [Departement de Physique, Universite des Sciences et Techniques de Masuku, B.P. 943, Franceville (Gabonese Republic)

2011-08-15

We examine the generation of bright matter-wave solitons in the Gross-Pitaevskii equation describing Bose-Einstein condensates with a time-dependent complex potential, which is composed of a repulsive parabolic background potential and a gravitational field. By performing a modified lens-type transformation, an explicit expression for the growth rate of a purely growing modulational instability is presented and analyzed. We point out the effects of the gravitational field, as well as of the parameter related to the feeding or loss of atoms in the condensate, on the instability growth rate. It is evident from numerical simulations that the feeding with atoms and the magnetic trap have opposite effects on the dynamics of the system. It is shown that the feeding or loss parameter can be well used to control the instability domain. Our study shows that the gravitational field changes the condensate trail of the soliton trains during the propagation. We also perform a numerical analysis to solve the Gross-Pitaevskii equation with a time-dependent complicated potential. The numerical results on the effect of both the gravitational field and the parameter of feeding or loss of atoms in the condensate agree well with predictions of the linear stability analysis. Another result of the present work is the modification of the background wave function in the Thomas-Fermi approximation during the numerical simulations.

On the initial conditions of time-dependent mean-field equations of evolution. Pt. 2

International Nuclear Information System (INIS)

Troudet, T.; Paris-11 Univ., 91 - Orsay

1986-01-01

We analyze the problem so far untouched of determining the initial mean-field wavefunction in the context of zero-temperature mean-field descriptions of time-dependent expectation values and quantum fluctuations of nuclear observables. The nucleus, at zero temperature, is taken to be in a low-lying excited many-body eigenstate and is approximated by the corresponding RPA wavefunction as a continuous superposition of coherent states (i.e. Slater determinants). A generating function Gsub(A)(lambda) for time-dependent expectation values and quantum fluctuations is constructed within the formalism of functional integration. By applying the saddle-point method to the functional action of Gsub(A)(lambda) and then taking its lambda-derivatives, we recover the well-known TDHF theory and propose a simple determination of the initial Slater determinant for an appropriate mean-field description of time-dependent expectation values. The analog mean-field description of quadratic-quantum fluctuations proceeds similarly and in addition includes the contribution of the uncorrelated TDHF-RPA phonons coupled to collective excitations of the initial (static) mean-field configuration. When the collective TDHF-RPA excitations are solely taken into account, we obtain an improved version of the Balian-Veneroni dispersion formula by showing how to determine the initial mean-field wavefunction. By first taking the lambda-derivatives of Gsub(A)(lambda) before applying the saddle-point method, the initial mean-field wavefunction is found to be non-linearly coupled to the mean-field dynamics themselves. In return, and in contrast to the first quantization scheme, these both depend non-trivially upon the observable A being measured so that approximations must be proposed to simplify the resulting mean-field equations. (orig.)

Time evolution of some quantum-mechanical systems. Wavefunction cloning in evolving rotating systems. Finite range boundary conditions for time dependent Schroedinger Equation

International Nuclear Information System (INIS)

Arvieu, R.; Carbonell, J.; Gignoux, C.; Mangin-Brinet, M.; Rozmej, P.

1997-01-01

The time evolution of coherent rotational wave packets associated to a diatomic molecule or to a deformed nucleus has been studied. Assuming a rigid body dynamics the J(J+1) law leads to a mechanism of cloning: the way function is divided into wave packets identical to the initial one at specific time. Applications are studied for a nuclear wave packed formed by Coulomb excitation. Exact boundary conditions at finite distance for the solution of the time-dependent Schroedinger equation are derived. A numerical scheme based on Crank-Nicholson method is proposed to illustrate its applicability in several examples. (authors)

Cellular neural networks (CNN) simulation for the TN approximation of the time dependent neutron transport equation in slab geometry

International Nuclear Information System (INIS)

Hadad, Kamal; Pirouzmand, Ahmad; Ayoobian, Navid

2008-01-01

This paper describes the application of a multilayer cellular neural network (CNN) to model and solve the time dependent one-speed neutron transport equation in slab geometry. We use a neutron angular flux in terms of the Chebyshev polynomials (T N ) of the first kind and then we attempt to implement the equations in an equivalent electrical circuit. We apply this equivalent circuit to analyze the T N moments equation in a uniform finite slab using Marshak type vacuum boundary condition. The validity of the CNN results is evaluated with numerical solution of the steady state T N moments equations by MATLAB. Steady state, as well as transient simulations, shows a very good comparison between the two methods. We used our CNN model to simulate space-time response of total flux and its moments for various c (where c is the mean number of secondary neutrons per collision). The complete algorithm could be implemented using very large-scale integrated circuit (VLSI) circuitry. The efficiency of the calculation method makes it useful for neutron transport calculations

Solution of the linearised Vlasov equation for collisionless plasmas evolving in external fields of arbitrary spatial and time dependence: Pt. 2

International Nuclear Information System (INIS)

Skarka, V.; Coveney, P.V.

1990-01-01

We solve perturbatively the linearised Vlasov equation describing inhomogeneous collisionless plasmas evolving in time-dependent external fields. The method employs an explicitly time-dependent formalism and is facilitated by the used of diagrammatic techniques. It leads to a straightforward algorithm for computing the contribution to the solution, order by order in the external field. In the previous paper we provided the solution to first order; higher orders are described in the present paper. (author)

Dynamics of second order in time evolution equations with state-dependent delay

Czech Academy of Sciences Publication Activity Database

Chueshov, I.; Rezunenko, Oleksandr

123-124, č. 1 (2015), s. 126-149 ISSN 0362-546X R&D Projects: GA ČR GAP103/12/2431 Institutional support: RVO:67985556 Keywords : Second order evolution equations * State dependent delay * Nonlinear plate * Finite-dimensional attractor Subject RIV: BD - Theory of Information Impact factor: 1.125, year: 2015 http://library.utia.cas.cz/separaty/2015/AS/rezunenko-0444708.pdf

Application of Trotter approximation for solving time dependent neutron transport equation; Primena Trotterove aproksimacije za resavanje vremenski zavisne transportne jednacine neutrona

Energy Technology Data Exchange (ETDEWEB)

Stancic, V [Institut za nuklearne nauke Boris Kidric, Vinca, Beograd (Yugoslavia)

1987-07-01

A method is proposed to solve multigroup time dependent neutron transport equation with arbitrary scattering anisotropy. The recurrence relation thus obtained is simple, numerically stable and especially suitable for treatment of complicated geometries. (author)

The constrained discrete-time state-dependent Riccati equation technique for uncertain nonlinear systems

Science.gov (United States)

Chang, Insu

The objective of the thesis is to introduce a relatively general nonlinear controller/estimator synthesis framework using a special type of the state-dependent Riccati equation technique. The continuous time state-dependent Riccati equation (SDRE) technique is extended to discrete-time under input and state constraints, yielding constrained (C) discrete-time (D) SDRE, referred to as CD-SDRE. For the latter, stability analysis and calculation of a region of attraction are carried out. The derivation of the D-SDRE under state-dependent weights is provided. Stability of the D-SDRE feedback system is established using Lyapunov stability approach. Receding horizon strategy is used to take into account the constraints on D-SDRE controller. Stability condition of the CD-SDRE controller is analyzed by using a switched system. The use of CD-SDRE scheme in the presence of constraints is then systematically demonstrated by applying this scheme to problems of spacecraft formation orbit reconfiguration under limited performance on thrusters. Simulation results demonstrate the efficacy and reliability of the proposed CD-SDRE. The CD-SDRE technique is further investigated in a case where there are uncertainties in nonlinear systems to be controlled. First, the system stability under each of the controllers in the robust CD-SDRE technique is separately established. The stability of the closed-loop system under the robust CD-SDRE controller is then proven based on the stability of each control system comprising switching configuration. A high fidelity dynamical model of spacecraft attitude motion in 3-dimensional space is derived with a partially filled fuel tank, assumed to have the first fuel slosh mode. The proposed robust CD-SDRE controller is then applied to the spacecraft attitude control system to stabilize its motion in the presence of uncertainties characterized by the first fuel slosh mode. The performance of the robust CD-SDRE technique is discussed. Subsequently

Boundary layer phenomena for differential-delay equations with state-dependent time lags: III

Science.gov (United States)

Mallet-Paret, John; Nussbaum, Roger D.

We consider a class of singularly perturbed delay-differential equations of the form ɛ ẋ(t)=f(x(t),x(t-r)), where r= r( x( t)) is a state-dependent delay. We study the asymptotic shape, as ɛ→0, of slowly oscillating periodic solutions. In particular, we show that the limiting shape of such solutions can be explicitly described by the solution of a pair of so-called max-plus equations. We are able thereby to characterize both the regular parts of the solution graph and the internal transition layers arising from the singular perturbation structure.

Velocity form of the Kohn-Sham frequency-dependent polarizability equations

International Nuclear Information System (INIS)

Bartolotti, L.J.

1987-01-01

A single equation is derived for the determination of the first-order correction to the frequency-dependent density, due to the perturbation of a time-varying electric field. This new expression for the first-order correction to the frequency-dependent Kohn-Sham amplitudes depends explicitly upon the velocity form of the dipole-moment operator and the square of the Kohn-Sham Hamiltonian

A time dependent mixing model to close PDF equations for transport in heterogeneous aquifers

Science.gov (United States)

Schüler, L.; Suciu, N.; Knabner, P.; Attinger, S.

2016-10-01

Probability density function (PDF) methods are a promising alternative to predicting the transport of solutes in groundwater under uncertainty. They make it possible to derive the evolution equations of the mean concentration and the concentration variance, used in moment methods. The mixing model, describing the transport of the PDF in concentration space, is essential for both methods. Finding a satisfactory mixing model is still an open question and due to the rather elaborate PDF methods, a difficult undertaking. Both the PDF equation and the concentration variance equation depend on the same mixing model. This connection is used to find and test an improved mixing model for the much easier to handle concentration variance. Subsequently, this mixing model is transferred to the PDF equation and tested. The newly proposed mixing model yields significantly improved results for both variance modelling and PDF modelling.

Analytic solutions of the multigroup space-time reactor kinetics equations

International Nuclear Information System (INIS)

Lee, C.E.; Rottler, S.

1986-01-01

The development of analytical and numerical solutions to the reactor kinetics equations is reviewed. Analytic solutions of the multigroup space-time reactor kinetics equations are developed for bare and reflected slabs and spherical reactors for zero flux, zero current and extrapolated endpoint boundary conditions. The material properties of the reactors are assumed constant in space and time, but spatially-dependent source terms and initial conditions are investigated. The system of partial differential equations is reduced to a set of linear ordinary differential equations by the Laplace transform method. These equations are solved by matrix Green's functions yielding a general matrix solution for the neutron flux and precursor concentration in the Laplace transform space. The detailed pole structure of the Laplace transform matrix solutions is investigated. The temporally- and spatially-dependent solutions are determined from the inverse Laplace transform using the Cauchy residue theorem, the theorem of Frobenius, a knowledge of the detailed pole structure and matrix operators. (author)

Time-Dependent Mean-Field Games in the Subquadratic Case

KAUST Repository

Gomes, Diogo A.; Pimentel, Edgard A.; Sá nchez-Morgado, Hé ector

2014-01-01

In this paper we consider time-dependent mean-field games with subquadratic Hamiltonians and power-like local dependence on the measure. We establish existence of classical solutions under a certain set of conditions depending on both the growth of the Hamiltonian and the dimension. This is done by combining regularity estimates for the Hamilton-Jacobi equation based on the Gagliardo-Nirenberg interpolation inequality with polynomial estimates for the Fokker-Planck equation. This technique improves substantially the previous results on the regularity of time-dependent mean-field games.

Time-Dependent Mean-Field Games in the Subquadratic Case

KAUST Repository

Gomes, Diogo A.

2014-10-14

In this paper we consider time-dependent mean-field games with subquadratic Hamiltonians and power-like local dependence on the measure. We establish existence of classical solutions under a certain set of conditions depending on both the growth of the Hamiltonian and the dimension. This is done by combining regularity estimates for the Hamilton-Jacobi equation based on the Gagliardo-Nirenberg interpolation inequality with polynomial estimates for the Fokker-Planck equation. This technique improves substantially the previous results on the regularity of time-dependent mean-field games.

Numerical Simulations of Light Bullets, Using The Full Vector, Time Dependent, Nonlinear Maxwell Equations

Science.gov (United States)

Goorjian, Peter M.; Silberberg, Yaron; Kwak, Dochan (Technical Monitor)

1995-01-01

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that we currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Karr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Time dependent drift Hamiltonian

International Nuclear Information System (INIS)

Boozer, A.H.

1982-04-01

The motion of individual charged particles in a given magnetic and an electric fields is discussed. An idea of a guiding center distribution function f is introduced. The guiding center distribution function is connected to the asymptotic Hamiltonian through the drift kinetic equation. The general non-stochastic magnetic field can be written in a contravariant and a covariant forms. The drift Hamiltonian is proposed, and the canonical gyroradius is presented. The proposed drift Hamiltonian agrees with Alfven's drift velocity to lowest non-vanishing order in the gyroradius. The relation between the exact, time dependent equations of motion and the guiding center equation is clarified by a Lagrangian analysis. The deduced Lagrangian represents the drift motion. (Kato, T.)

On solving wave equations on fixed bounded intervals involving Robin boundary conditions with time-dependent coefficients

Science.gov (United States)

van Horssen, Wim T.; Wang, Yandong; Cao, Guohua

2018-06-01

In this paper, it is shown how characteristic coordinates, or equivalently how the well-known formula of d'Alembert, can be used to solve initial-boundary value problems for wave equations on fixed, bounded intervals involving Robin type of boundary conditions with time-dependent coefficients. A Robin boundary condition is a condition that specifies a linear combination of the dependent variable and its first order space-derivative on a boundary of the interval. Analytical methods, such as the method of separation of variables (SOV) or the Laplace transform method, are not applicable to those types of problems. The obtained analytical results by applying the proposed method, are in complete agreement with those obtained by using the numerical, finite difference method. For problems with time-independent coefficients in the Robin boundary condition(s), the results of the proposed method also completely agree with those as for instance obtained by the method of separation of variables, or by the finite difference method.

Functional approach to a time-dependent self-consistent field theory

International Nuclear Information System (INIS)

Reinhardt, H.

1979-01-01

The time-dependent Hartree-Fock approximation is formulated within the path integral approach. It is shown that by a suitable choice of the collective field the classical equation of motion of the collective field coincides with the time-dependent Hartree (TDH) equation. The consideration is restricted to the TDH equation, since the exchange terms do not appear in the functional approach on the same footing as the direct terms

Exact wavefunctions for a time-dependent Coulomb potential

International Nuclear Information System (INIS)

Menouar, S; Maamache, M; Saadi, Y; Choi, J R

2008-01-01

The one-dimensional Schroedinger equation associated with a time-dependent Coulomb potential is studied. The invariant operator method (Lewis and Riesenfeld) and unitary transformation approach are employed to derive quantum solutions of the system. We obtain an ordinary second-order differential equation whose analytical exact solution has been unknown. It is confirmed that the form of this equation is similar to the radial Schroedinger equation for the hydrogen atom in a (arbitrary) strong magnetic field. The qualitative properties for the eigenstates spectrum are described separately for the different values of the parameter ω 0 appearing in the x 2 term, x being the position, i.e., ω 0 > 0, ω 0 0 = 0. For the ω 0 = 0 case, the eigenvalue equation of invariant operator reduces to a solvable form and, consequently, we have provided exact eigenstates of the time-dependent Hamiltonian system

Relativistic extension of a charge-conservative finite element solver for time-dependent Maxwell-Vlasov equations

Science.gov (United States)

Na, D.-Y.; Moon, H.; Omelchenko, Y. A.; Teixeira, F. L.

2018-01-01

Accurate modeling of relativistic particle motion is essential for physical predictions in many problems involving vacuum electronic devices, particle accelerators, and relativistic plasmas. A local, explicit, and charge-conserving finite-element time-domain (FETD) particle-in-cell (PIC) algorithm for time-dependent (non-relativistic) Maxwell-Vlasov equations on irregular (unstructured) meshes was recently developed by Moon et al. [Comput. Phys. Commun. 194, 43 (2015); IEEE Trans. Plasma Sci. 44, 1353 (2016)]. Here, we extend this FETD-PIC algorithm to the relativistic regime by implementing and comparing three relativistic particle-pushers: (relativistic) Boris, Vay, and Higuera-Cary. We illustrate the application of the proposed relativistic FETD-PIC algorithm for the analysis of particle cyclotron motion at relativistic speeds, harmonic particle oscillation in the Lorentz-boosted frame, and relativistic Bernstein modes in magnetized charge-neutral (pair) plasmas.

On the time-dependent Aharonov–Bohm effect

Directory of Open Access Journals (Sweden)

Jian Jing

2017-11-01

Full Text Available The Aharonov–Bohm effect in the background of a time-dependent vector potential is re-examined for both non-relativistic and relativistic cases. Based on the solutions to the Schrodinger and Dirac equations which contain the time-dependent magnetic vector potential, we find that contrary to the conclusions in a recent paper (Singleton and Vagenas 2013 [4], the interference pattern will be altered with respect to time because of the time-dependent vector potential.

Attainable conditions and exact invariant for the time-dependent harmonic oscillator

Energy Technology Data Exchange (ETDEWEB)

Guasti, Manuel Fernandez [Lab. de Optica Cuantica, Dep. de Fisica, Universidad A. Metropolitana, Unidad Iztapalapa, Mexico DF, Ap. Post. 55-534 (Mexico)

2006-09-22

The time-dependent oscillator equation is solved numerically for various trajectories in amplitude and phase variables. The solutions exhibit a finite time-dependent parameter whenever the squared amplitude times the derivative of the phase is invariant. If the invariant relationship does not hold, the time-dependent parameter has divergent singularities. These observations lead to the proposition that the harmonic oscillator equation with finite time-dependent parameter must have amplitude and phase solutions fulfilling the invariant relationship. Since the time-dependent parameter or the potential must be finite for any real oscillator implementation, the invariant must hold for any such physically realizable system.

Attainable conditions and exact invariant for the time-dependent harmonic oscillator

International Nuclear Information System (INIS)

Guasti, Manuel Fernandez

2006-01-01

The time-dependent oscillator equation is solved numerically for various trajectories in amplitude and phase variables. The solutions exhibit a finite time-dependent parameter whenever the squared amplitude times the derivative of the phase is invariant. If the invariant relationship does not hold, the time-dependent parameter has divergent singularities. These observations lead to the proposition that the harmonic oscillator equation with finite time-dependent parameter must have amplitude and phase solutions fulfilling the invariant relationship. Since the time-dependent parameter or the potential must be finite for any real oscillator implementation, the invariant must hold for any such physically realizable system

Construction of an exact solution of time-dependent Ginzburg ...

Indian Academy of Sciences (India)

time-dependent Ginzburg–Landau (TDGL) equations we have calculated the ... The prototype of such equations is the parabolic reaction diffusion equation [7,8] ..... It may be possible to compare the above results with suitable experiments, ...

Time-dependent problems and difference methods

CERN Document Server

Gustafsson, Bertil; Oliger, Joseph

2013-01-01

Praise for the First Edition "". . . fills a considerable gap in the numerical analysis literature by providing a self-contained treatment . . . this is an important work written in a clear style . . . warmly recommended to any graduate student or researcher in the field of the numerical solution of partial differential equations."" -SIAM Review Time-Dependent Problems and Difference Methods, Second Edition continues to provide guidance for the analysis of difference methods for computing approximate solutions to partial differential equations for time-de

Time-dependent deterministic transport on parallel architectures using PARTISN

International Nuclear Information System (INIS)

Alcouffe, R.E.; Baker, R.S.

1998-01-01

In addition to the ability to solve the static transport equation, the authors have also incorporated time dependence into the parallel S N code PARTISN. Using a semi-implicit scheme, PARTISN is capable of performing time-dependent calculations for both fissioning and pure source driven problems. They have applied this to various types of problems such as shielding and prompt fission experiments. This paper describes the form of the time-dependent equations implemented, their solution strategies in PARTISN including iteration acceleration, and the strategies used for time-step control. Results are presented for a iron-water shielding calculation and a criticality excursion in a uranium solution configuration

Introducing time-dependent molecular fields: a new derivation of the wave equations

Science.gov (United States)

Baer, Michael

2018-02-01

This article is part of a series of articles trying to establish the concept molecular field. The theory that induced us to introduce this novel concept is based on the Born-Huang expansion as applied to the Schroedinger equation that describes the interaction of a molecular system with an external electric field. Assuming the molecular system is made up of two coupled adiabatic states the theory leads from a single spatial curl equation, two space-time curl equations and one single space-time divergent equation to a pair of decoupled wave equations usually encountered within the theory of fields. In the present study, just like in the previous study [see Baer et al., Mol. Phys. 114, 227 (2016)] the wave equations are derived for an electric field having two features: (a) its intensity is high enough; (b) its duration is short enough. Although not all the findings are new the derivation, in the present case, is new, straightforward, fluent and much friendlier as compared to the previous one and therefore should be presented again. For this situation the study reveals that the just described interaction creates two fields that coexist within a molecule: one is a novel vectorial field formed via the interaction of the electric field with the Born-Huang non-adiabatic coupling terms (NACTs) and the other is an ordinary, scalar, electric field essentially identical to the original electric field. Section 4 devoted to the visualization of the outcomes via two intersecting Jahn-Teller cones which contain NACTs that become singular at the intersection point of these cones. Finally, the fact that eventually we are facing a kind of a cosmic situation may bring us to speculate that singular NACTs are a result of cosmic phenomena. Thus, if indeed this singularity is somehow connected to reality then, like other singularities in physics, it is formed at (or immediately after) the Big Bang and consequently, guarantees the formation of molecules.

Time-dependent quantum transport through an interacting quantum dot beyond sequential tunneling: second-order quantum rate equations

International Nuclear Information System (INIS)

Dong, B; Ding, G H; Lei, X L

2015-01-01

A general theoretical formulation for the effect of a strong on-site Coulomb interaction on the time-dependent electron transport through a quantum dot under the influence of arbitrary time-varying bias voltages and/or external fields is presented, based on slave bosons and the Keldysh nonequilibrium Green's function (GF) techniques. To avoid the difficulties of computing double-time GFs, we generalize the propagation scheme recently developed by Croy and Saalmann to combine the auxiliary-mode expansion with the celebrated Lacroix's decoupling approximation in dealing with the second-order correlated GFs and then establish a closed set of coupled equations of motion, called second-order quantum rate equations (SOQREs), for an exact description of transient dynamics of electron correlated tunneling. We verify that the stationary solution of our SOQREs is able to correctly describe the Kondo effect on a qualitative level. Moreover, a comparison with other methods, such as the second-order von Neumann approach and Hubbard-I approximation, is performed. As illustrations, we investigate the transient current behaviors in response to a step voltage pulse and a harmonic driving voltage, and linear admittance as well, in the cotunneling regime. (paper)

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.

Motivation for Using Generalized Geometry in the Time Dependent Transport Code TDKENO

Energy Technology Data Exchange (ETDEWEB)

Dustin Popp; Zander Mausolff; Sedat Goluoglu

2016-04-01

We are proposing to use the code, TDKENO, to model TREAT. TDKENO solves the time dependent, three dimensional Boltzmann transport equation with explicit representation of delayed neutrons. Instead of directly integrating this equation, the neutron flux is factored into two components – a rapidly varying amplitude equation and a slowly varying shape equation and each is solved separately on different time scales. The shape equation is solved using the 3D Monte Carlo transport code KENO, from Oak Ridge National Laboratory’s SCALE code package. Using the Monte Carlo method to solve the shape equation is still computationally intensive, but the operation is only performed when needed. The amplitude equation is solved deterministically and frequently, so the solution gives an accurate time-dependent solution without having to repeatedly We have modified TDKENO to incorporate KENO-VI so that we may accurately represent the geometries within TREAT. This paper explains the motivation behind using generalized geometry, and provides the results of our modifications. TDKENO uses the Improved Quasi-Static method to accomplish this. In this method, the neutron flux is factored into two components. One component is a purely time-dependent and rapidly varying amplitude function, which is solved deterministically and very frequently (small time steps). The other is a slowly varying flux shape function that weakly depends on time and is only solved when needed (significantly larger time steps).

Time-dependent 2-stream particle transport

International Nuclear Information System (INIS)

Corngold, Noel

2015-01-01

Highlights: • We consider time-dependent transport in the 2-stream or “rod” model via an attractive matrix formalism. • After reviewing some classical problems in homogeneous media we discuss transport in materials with whose density may vary. • There we achieve a significant contraction of the underlying Telegrapher’s equation. • We conclude with a discussion of stochastics, treated by the “first-order smoothing approximation.” - Abstract: We consider time-dependent transport in the 2-stream or “rod” model via an attractive matrix formalism. After reviewing some classical problems in homogeneous media we discuss transport in materials whose density may vary. There we achieve a significant contraction of the underlying Telegrapher’s equation. We conclude with a discussion of stochastics, treated by the “first-order smoothing approximation.”

Time dependent mean-field games

KAUST Repository

Gomes, Diogo A.

2014-01-06

We consider time dependent mean-field games (MFG) with a local power-like dependence on the measure and Hamiltonians satisfying both sub and superquadratic growth conditions. We establish existence of smooth solutions under a certain set of conditions depending both on the growth of the Hamiltonian as well as on the dimension. In the subquadratic case this is done by combining a Gagliardo-Nirenberg type of argument with a new class of polynomial estimates for solutions of the Fokker-Planck equation in terms of LrLp- norms of DpH. These techniques do not apply to the superquadratic case. In this setting we recur to a delicate argument that combines the non-linear adjoint method with polynomial estimates for solutions of the Fokker-Planck equation in terms of L1L1-norms of DpH. Concerning the subquadratic case, we substantially improve and extend the results previously obtained. Furthermore, to the best of our knowledge, the superquadratic case has not been addressed in the literature yet. In fact, it is likely that our estimates may also add to the current understanding of Hamilton-Jacobi equations with superquadratic Hamiltonians.

Time-dependent, Bianchi II, rotating universe

International Nuclear Information System (INIS)

Reboucas, M.J.

1981-01-01

An exact cosmological solution of Einstein's equations which has time-dependent rotation is presented. The t-constant sections are of Bianchi type II. The source of this geometry is a fluid which has not been thermalized. (Author) [pt

Variational derivation of a time-dependent Hartree-Fock Hamiltonian

International Nuclear Information System (INIS)

Lichtner, P.C.; Griffin, J.J.; Schultheis, H.; Schultheis, R.; Volkov, A.B.

1979-01-01

The variational derivation of the time-dependent Hartree-Fock equation is reviewed. When norm-violating variations are included, a unique time-dependent Hartree-Fock Hamiltonian, which differs from that customarily used in time-dependent Hartree-Fock analyses, is implied. This variationally ''true'' Hartree-Fock Hamiltonian has the same expectation value as the exact Hamiltonian, equal to the average energy of the system. Since this quantity remains constant under time-dependent Hartree-Fock time evolution, we suggest the label ''constant '' for this form of time-dependent Hartree-Fock theory

A NUMERICAL SCHEME FOR SPECIAL RELATIVISTIC RADIATION MAGNETOHYDRODYNAMICS BASED ON SOLVING THE TIME-DEPENDENT RADIATIVE TRANSFER EQUATION

Energy Technology Data Exchange (ETDEWEB)

Ohsuga, Ken; Takahashi, Hiroyuki R. [National Astronomical Observatory of Japan, Osawa, Mitaka, Tokyo 181-8588 (Japan)

2016-02-20

We develop a numerical scheme for solving the equations of fully special relativistic, radiation magnetohydrodynamics (MHDs), in which the frequency-integrated, time-dependent radiation transfer equation is solved to calculate the specific intensity. The radiation energy density, the radiation flux, and the radiation stress tensor are obtained by the angular quadrature of the intensity. In the present method, conservation of total mass, momentum, and energy of the radiation magnetofluids is guaranteed. We treat not only the isotropic scattering but also the Thomson scattering. The numerical method of MHDs is the same as that of our previous work. The advection terms are explicitly solved, and the source terms, which describe the gas–radiation interaction, are implicitly integrated. Our code is suitable for massive parallel computing. We present that our code shows reasonable results in some numerical tests for propagating radiation and radiation hydrodynamics. Particularly, the correct solution is given even in the optically very thin or moderately thin regimes, and the special relativistic effects are nicely reproduced.

An “Airy gun”: Self-accelerating solutions of the time-dependent Schrödinger equation in vacuum

International Nuclear Information System (INIS)

Mahalov, Alex; Suslov, Sergei K.

2012-01-01

We consider generalizations of the Berry and Balazs nonspreading and accelerating solution of the time-dependent Schrödinger equation in empty space, which has been experimentally demonstrated in paraxial optics. In particular, we show that the original nonspreading wave packet is unstable. An explicit variation of the initial Airy-state evolves into the self-accelerating and self-compressing solution presented here. Quasi-diffraction-free finite energy Airy beams that are more realistic for experimental study are obtained by analytic continuation and their Wigner function is evaluated. Nonlinear generalizations related to second Painlevé transcendents are briefly discussed.

Blow-Up Time for Nonlinear Heat Equations with Transcendental Nonlinearity

Directory of Open Access Journals (Sweden)

Hee Chul Pak

2012-01-01

Full Text Available A blow-up time for nonlinear heat equations with transcendental nonlinearity is investigated. An upper bound of the first blow-up time is presented. It is pointed out that the upper bound of the first blow-up time depends on the support of the initial datum.

Time dependent variational method in quantum mechanics

International Nuclear Information System (INIS)

Torres del Castillo, G.F.

1987-01-01

Using the fact that the solutions to the time-dependent Schodinger equation can be obtained from a variational principle, by restricting the evolution of the state vector to some surface in the corresponding Hilbert space, approximations to the exact solutions can be obtained, which are determined by equations similar to Hamilton's equations. It is shown that, in order for the approximate evolution to be well defined on a given surface, the imaginary part of the inner product restricted to the surface must be non-singular. (author)

Space-Time Fractional Diffusion-Advection Equation with Caputo Derivative

Directory of Open Access Journals (Sweden)

José Francisco Gómez Aguilar

2014-01-01

Full Text Available An alternative construction for the space-time fractional diffusion-advection equation for the sedimentation phenomena is presented. The order of the derivative is considered as 0<β, γ≤1 for the space and time domain, respectively. The fractional derivative of Caputo type is considered. In the spatial case we obtain the fractional solution for the underdamped, undamped, and overdamped case. In the temporal case we show that the concentration has amplitude which exhibits an algebraic decay at asymptotically large times and also shows numerical simulations where both derivatives are taken in simultaneous form. In order that the equation preserves the physical units of the system two auxiliary parameters σx and σt are introduced characterizing the existence of fractional space and time components, respectively. A physical relation between these parameters is reported and the solutions in space-time are given in terms of the Mittag-Leffler function depending on the parameters β and γ. The generalization of the fractional diffusion-advection equation in space-time exhibits anomalous behavior.

Isostable reduction with applications to time-dependent partial differential equations.

Science.gov (United States)

Wilson, Dan; Moehlis, Jeff

2016-07-01

Isostables and isostable reduction, analogous to isochrons and phase reduction for oscillatory systems, are useful in the study of nonlinear equations which asymptotically approach a stationary solution. In this work, we present a general method for isostable reduction of partial differential equations, with the potential power to reduce the dimensionality of a nonlinear system from infinity to 1. We illustrate the utility of this reduction by applying it to two different models with biological relevance. In the first example, isostable reduction of the Fokker-Planck equation provides the necessary framework to design a simple control strategy to desynchronize a population of pathologically synchronized oscillatory neurons, as might be relevant to Parkinson's disease. Another example analyzes a nonlinear reaction-diffusion equation with relevance to action potential propagation in a cardiac system.

Subdiffusive master equation with space-dependent anomalous exponent and structural instability

Science.gov (United States)

Fedotov, Sergei; Falconer, Steven

2012-03-01

We derive the fractional master equation with space-dependent anomalous exponent. We analyze the asymptotic behavior of the corresponding lattice model both analytically and by Monte Carlo simulation. We show that the subdiffusive fractional equations with constant anomalous exponent μ in a bounded domain [0,L] are not structurally stable with respect to the nonhomogeneous variations of parameter μ. In particular, the Gibbs-Boltzmann distribution is no longer the stationary solution of the fractional Fokker-Planck equation whatever the space variation of the exponent might be. We analyze the random distribution of μ in space and find that in the long-time limit, the probability distribution is highly intermediate in space and the behavior is completely dominated by very unlikely events. We show that subdiffusive fractional equations with the nonuniform random distribution of anomalous exponent is an illustration of a “Black Swan,” the low probability event of the small value of the anomalous exponent that completely dominates the long-time behavior of subdiffusive systems.

The Harmonic Potential Theorem for a Quantum System with Time-Dependent Effective Mass

International Nuclear Information System (INIS)

Lai Meng-Yun; Xiao Duan-Liang; Pan Xiao-Yin

2015-01-01

We investigate the many-body wave function of a quantum system with time-dependent effective mass, confined by a harmonic potential with time-dependent frequency, and perturbed by a time-dependent spatially homogeneous electric field. It is found that the wave function is comprised of a phase factor times the solution to the unperturbed time-dependent Schrödinger equation with the latter being translated by a time-dependent value that satisfies the classical driven equation of motion. The wave function reduces to that of the harmonic potential theorem wave function when both the effective mass and frequency are static. An example of application is also given. (paper)

A time-dependent neutron transport method of characteristics formulation with time derivative propagation

Energy Technology Data Exchange (ETDEWEB)

Hoffman, Adam J., E-mail: adamhoff@umich.edu; Lee, John C., E-mail: jcl@umich.edu

2016-02-15

A new time-dependent Method of Characteristics (MOC) formulation for nuclear reactor kinetics was developed utilizing angular flux time-derivative propagation. This method avoids the requirement of storing the angular flux at previous points in time to represent a discretized time derivative; instead, an equation for the angular flux time derivative along 1D spatial characteristics is derived and solved concurrently with the 1D transport characteristic equation. This approach allows the angular flux time derivative to be recast principally in terms of the neutron source time derivatives, which are approximated to high-order accuracy using the backward differentiation formula (BDF). This approach, called Source Derivative Propagation (SDP), drastically reduces the memory requirements of time-dependent MOC relative to methods that require storing the angular flux. An SDP method was developed for 2D and 3D applications and implemented in the computer code DeCART in 2D. DeCART was used to model two reactor transient benchmarks: a modified TWIGL problem and a C5G7 transient. The SDP method accurately and efficiently replicated the solution of the conventional time-dependent MOC method using two orders of magnitude less memory.

A time-dependent neutron transport method of characteristics formulation with time derivative propagation

International Nuclear Information System (INIS)

Hoffman, Adam J.; Lee, John C.

2016-01-01

A new time-dependent Method of Characteristics (MOC) formulation for nuclear reactor kinetics was developed utilizing angular flux time-derivative propagation. This method avoids the requirement of storing the angular flux at previous points in time to represent a discretized time derivative; instead, an equation for the angular flux time derivative along 1D spatial characteristics is derived and solved concurrently with the 1D transport characteristic equation. This approach allows the angular flux time derivative to be recast principally in terms of the neutron source time derivatives, which are approximated to high-order accuracy using the backward differentiation formula (BDF). This approach, called Source Derivative Propagation (SDP), drastically reduces the memory requirements of time-dependent MOC relative to methods that require storing the angular flux. An SDP method was developed for 2D and 3D applications and implemented in the computer code DeCART in 2D. DeCART was used to model two reactor transient benchmarks: a modified TWIGL problem and a C5G7 transient. The SDP method accurately and efficiently replicated the solution of the conventional time-dependent MOC method using two orders of magnitude less memory.

Exactly solvable position dependent mass schroedinger equation

International Nuclear Information System (INIS)

Koc, R.; Tuetuencueler, H.; Koercuek, E.

2002-01-01

Exact solution of the Schrodinger equation with a variable mass is presented. We have derived general expressions for the eigenstates and eigenvalues of the position dependent mass systems. We provide supersymmetric and Lie algebraic methods to discuss the position dependent mass systems

''Free-space'' boundary conditions for the time-dependent wave equation

International Nuclear Information System (INIS)

Lindman, E.L.

1975-01-01

Boundary conditions for the discrete wave equation which act like an infinite region of free space in contact with the computational region can be constructed using projection operators. Propagating and evanescent waves coming from within the computational region generate no reflected waves as they cross the boundary. At the same time arbitrary waves may be launched into the computational region. Well known projection operators for one-dimensional waves may be used for this purpose in one dimension. Extensions of these operators to higher dimensions along with numerically efficient approximations to them are described for higher-dimensional problems. The separation of waves into ingoing and outgoing waves inherent in these boundary conditions greatly facilitates diagnostics

Equation for disentangling time-ordered exponentials with arbitrary quadratic generators

International Nuclear Information System (INIS)

Budanov, V.G.

1987-01-01

In many quantum-mechanical constructions, it is necessary to disentangle an operator-valued time-ordered exponential with time-dependent generators quadratic in the creation and annihilation operators. By disentangling, one understands the finding of the matrix elements of the time-ordered exponential or, in a more general formulation. The solution of the problem can also be reduced to calculation of a matrix time-ordered exponential that solves the corresponding classical problem. However, in either case the evolution equations in their usual form do not enable one to take into account explicitly the symmetry of the system. In this paper the methods of Weyl analysis are used to find an ordinary differential equation on a matrix Lie algebra that is invariant with respect to the adjoint action of the dynamical symmetry group of a quadratic Hamiltonian and replaces the operator evolution equation for the Green's function

First passage time problems in time-dependent fields

International Nuclear Information System (INIS)

Fletcher, J.E.; Havlin, S.; Weiss, G.H.

1988-01-01

This paper discusses the simplest first passage time problems for random walks and diffusion processes on a line segment. When a diffusing particle moves in a time-varying field, use of the adjoint equation does not lead to any simplification in the calculation of moments of the first passage time as is the case for diffusion in a time-invariant field. We show that for a discrete random walk in the presence of a sinusoidally varying field there is a resonant frequency omega* for which the mean residence time on the line segment in a minimum. It is shown that for a random walk on a line segment of length L the mean residence time goes like L 2 for large L when omega omega*, but when omega = omega* the dependence is proportional to L. The results of our simulation are numerical, but can be regarded as exact. Qualitatively similar results are shown to hold for diffusion processes by a perturbation expansion in powers of a dimensionless velocity. These results are extended to higher values of this parameter by a numerical solution of the forward equation

Time-dependent tumour repopulation factors in linear-quadratic equations

International Nuclear Information System (INIS)

Dale, R.G.

1989-01-01

Tumour proliferation effects can be tentatively quantified in the linear-quadratic (LQ) method by the incorporation of a time-dependent factor, the magnitude of which is related both to the value of α in the tumour α/β ratio, and to the tumour doubling time. The method, the principle of which has been suggested by a numbre of other workers for use in fractionated therapy, is here applied to both fractionated and protracted radiotherapy treatments, and examples of its uses are given. By assuming that repopulation of late-responding tissues is significant during normal treatment strategies in terms of the behaviour of the Extrapolated Response Dose (ERD). Although the numerical credibility of the analysis used here depends on the reliability of the LQ model, and on the assumption that the rate of repopulation is constant throughout treatment, the predictions are consistent with other lines of reasoning which point to the advantages of accelerated hyperfractionation. In particular, it is demonstrated that accelerated fractionation represents a relatively 'foregiving' treatment which enables tumours of a variety of sensitivities and clonogenic growth rates to be treated moderately successfully, even though the critical cellular parameters may not be known in individual cases. The analysis also suggests that tumours which combine low intrinsic sensitivity with a very short doubling time might be bettter controlled by low dose-rate continuous therapy than by almost any form of accelerated hyperfractionation. (author). 24 refs.; 5 figs

An explicit marching on-in-time solver for the time domain volume magnetic field integral equation

KAUST Repository

Sayed, Sadeed Bin

2014-07-01

Transient scattering from inhomogeneous dielectric objects can be modeled using time domain volume integral equations (TDVIEs). TDVIEs are oftentimes solved using marching on-in-time (MOT) techniques. Classical MOT-TDVIE solvers expand the field induced on the scatterer using local spatio-temporal basis functions. Inserting this expansion into the TDVIE and testing the resulting equation in space and time yields a system of equations that is solved by time marching. Depending on the type of the basis and testing functions and the time step, the time marching scheme can be implicit (N. T. Gres, et al., Radio Sci., 36(3), 379-386, 2001) or explicit (A. Al-Jarro, et al., IEEE Trans. Antennas Propag., 60(11), 5203-5214, 2012). Implicit MOT schemes are known to be more stable and accurate. However, under low-frequency excitation, i.e., when the time step size is large, they call for inversion of a full matrix system at very time step.

An explicit marching on-in-time solver for the time domain volume magnetic field integral equation

KAUST Repository

Sayed, Sadeed Bin; Ulku, Huseyin Arda; Bagci, Hakan

2014-01-01

Transient scattering from inhomogeneous dielectric objects can be modeled using time domain volume integral equations (TDVIEs). TDVIEs are oftentimes solved using marching on-in-time (MOT) techniques. Classical MOT-TDVIE solvers expand the field induced on the scatterer using local spatio-temporal basis functions. Inserting this expansion into the TDVIE and testing the resulting equation in space and time yields a system of equations that is solved by time marching. Depending on the type of the basis and testing functions and the time step, the time marching scheme can be implicit (N. T. Gres, et al., Radio Sci., 36(3), 379-386, 2001) or explicit (A. Al-Jarro, et al., IEEE Trans. Antennas Propag., 60(11), 5203-5214, 2012). Implicit MOT schemes are known to be more stable and accurate. However, under low-frequency excitation, i.e., when the time step size is large, they call for inversion of a full matrix system at very time step.

Time-dependent fatigue--phenomenology and life prediction

International Nuclear Information System (INIS)

Coffin, L.F.

1979-01-01

The time-dependent fatigue behavior of materials used or considered for use in present and advanced systems for power generation is outlined. A picture is first presented to show how basic mechanisms and phenomenological information relate to the performance of the component under consideration through the so-called local strain approach. By this means life prediction criteria and design rules can be formulated utilizing laboratory test information which is directly translated to predicting the performance of a component. The body of phenomenological information relative to time-dependent fatigue is reviewed. Included are effects of strain range, strain rate and frequency, environment and wave shape, all of which are shown to be important in developing both an understanding and design base for time dependent fatigue. Using this information, some of the current methods being considered for the life prediction of components are reviewed. These include the current ASME code case, frequency-modified fatigue equations, strain range partitioning, the damage function method, frequency separation and damage rate equations. From this review, it is hoped that a better perspective on future directions for basic material science at high temperature can be achieved

Efficient exact-exchange time-dependent density-functional theory methods and their relation to time-dependent Hartree-Fock.

Science.gov (United States)

Hesselmann, Andreas; Görling, Andreas

2011-01-21

A recently introduced time-dependent exact-exchange (TDEXX) method, i.e., a response method based on time-dependent density-functional theory that treats the frequency-dependent exchange kernel exactly, is reformulated. In the reformulated version of the TDEXX method electronic excitation energies can be calculated by solving a linear generalized eigenvalue problem while in the original version of the TDEXX method a laborious frequency iteration is required in the calculation of each excitation energy. The lowest eigenvalues of the new TDEXX eigenvalue equation corresponding to the lowest excitation energies can be efficiently obtained by, e.g., a version of the Davidson algorithm appropriate for generalized eigenvalue problems. Alternatively, with the help of a series expansion of the new TDEXX eigenvalue equation, standard eigensolvers for large regular eigenvalue problems, e.g., the standard Davidson algorithm, can be used to efficiently calculate the lowest excitation energies. With the help of the series expansion as well, the relation between the TDEXX method and time-dependent Hartree-Fock is analyzed. Several ways to take into account correlation in addition to the exact treatment of exchange in the TDEXX method are discussed, e.g., a scaling of the Kohn-Sham eigenvalues, the inclusion of (semi)local approximate correlation potentials, or hybrids of the exact-exchange kernel with kernels within the adiabatic local density approximation. The lowest lying excitations of the molecules ethylene, acetaldehyde, and pyridine are considered as examples.

Time-dependent potential-functional embedding theory

International Nuclear Information System (INIS)

Huang, Chen; Libisch, Florian; Peng, Qing; Carter, Emily A.

2014-01-01

We introduce a time-dependent potential-functional embedding theory (TD-PFET), in which atoms are grouped into subsystems. In TD-PFET, subsystems can be propagated by different suitable time-dependent quantum mechanical methods and their interactions can be treated in a seamless, first-principles manner. TD-PFET is formulated based on the time-dependent quantum mechanics variational principle. The action of the total quantum system is written as a functional of the time-dependent embedding potential, i.e., a potential-functional formulation. By exploiting the Runge-Gross theorem, we prove the uniqueness of the time-dependent embedding potential under the constraint that all subsystems share a common embedding potential. We derive the integral equation that such an embedding potential needs to satisfy. As proof-of-principle, we demonstrate TD-PFET for a Na 4 cluster, in which each Na atom is treated as one subsystem and propagated by time-dependent Kohn-Sham density functional theory (TDDFT) using the adiabatic local density approximation (ALDA). Our results agree well with a direct TDDFT calculation on the whole Na 4 cluster using ALDA. We envision that TD-PFET will ultimately be useful for studying ultrafast quantum dynamics in condensed matter, where key regions are solved by highly accurate time-dependent quantum mechanics methods, and unimportant regions are solved by faster, less accurate methods

Davydov–Chaban Hamiltonian in presence of time-dependent potential

Energy Technology Data Exchange (ETDEWEB)

Sobhani, Hadi; Hassanabadi, Hassan, E-mail: h.hassanabadi@shahroodut.ac.ir

2016-09-10

In this article, we have investigated collective effects of atomic nuclei in presence of a time-dependent potential in Davydov–Chaban Hamiltonian. Since such potential has an explicit time-dependency, in order to obtain the wave function of considered system, we should face with time-dependent Schrödinger equation. Obtaining the wave function could be possible using Lewis–Riesenfeld dynamical invariant method. Appropriate dynamical invariant has been constructed after determining the wave functions and values, the wave function will obtain.

Time-dependent mean-field games in the superquadratic case

KAUST Repository

Gomes, Diogo A.

2016-04-06

We investigate time-dependent mean-field games with superquadratic Hamiltonians and a power dependence on the measure. Such problems pose substantial mathematical challenges as key techniques used in the subquadratic case, which was studied in a previous publication of the authors, do not extend to the superquadratic setting. The main objective of the present paper is to address these difficulties. Because of the superquadratic structure of the Hamiltonian, Lipschitz estimates for the solutions of the Hamilton−Jacobi equation are obtained here through a novel set of techniques. These explore the parabolic nature of the problem through the nonlinear adjoint method. Well-posedness is proven by combining Lipschitz regularity for the Hamilton−Jacobi equation with polynomial estimates for solutions of the Fokker−Planck equation. Existence of classical solutions is then established under conditions depending only on the growth of the Hamiltonian and the dimension. Our results also add to current understanding of superquadratic Hamilton−Jacobi equations.

Time-dependent mean-field games in the superquadratic case

KAUST Repository

Gomes, Diogo A.; Pimentel, Edgard; Sá nchez-Morgado, Hé ctor

2016-01-01

We investigate time-dependent mean-field games with superquadratic Hamiltonians and a power dependence on the measure. Such problems pose substantial mathematical challenges as key techniques used in the subquadratic case, which was studied in a previous publication of the authors, do not extend to the superquadratic setting. The main objective of the present paper is to address these difficulties. Because of the superquadratic structure of the Hamiltonian, Lipschitz estimates for the solutions of the Hamilton−Jacobi equation are obtained here through a novel set of techniques. These explore the parabolic nature of the problem through the nonlinear adjoint method. Well-posedness is proven by combining Lipschitz regularity for the Hamilton−Jacobi equation with polynomial estimates for solutions of the Fokker−Planck equation. Existence of classical solutions is then established under conditions depending only on the growth of the Hamiltonian and the dimension. Our results also add to current understanding of superquadratic Hamilton−Jacobi equations.

Preliminary Formulation of Finite Element Solution for the 1-D, 1-G Time Dependent Neutron Diffusion Equation without Consideration about Delay Neutron

Energy Technology Data Exchange (ETDEWEB)

Ryu, Eun Hyun; Song, Yong Mann; Park, Joo Hwan [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of)

2013-05-15

If time-dependent equation is solved with the FEM, the limitation of the input geometry will disappear. It has often been pointed out that the numerical methods implemented in the RFSP code are not state-of-the-art. Although an acceleration method such as the Coarse Mesh Finite Difference (CMFD) for Finite Difference Method (FDM) does not exist for the FEM, one should keep in mind that the number of time steps for the transient simulation is not large. The rigorous formulation in this study will richen the theoretical basis of the FEM and lead to an extension of the dynamics code to deal with a more complicated problem. In this study, the formulation for the 1-D, 1-G Time Dependent Neutron Diffusion Equation (TDNDE) without consideration of the delay neutron will first be done. A problem including one multiplying medium will be solved. Also several conclusions from a comparison between the numerical and analytic solutions, a comparison between solutions with various element orders, and a comparison between solutions with different time differencing will be made to be certain about the formulation and FEM solution. By investigating various cases with different values of albedo, theta, and the order of elements, it can be concluded that the finite element solution is agree well with the analytic solution. The higher the element order used, the higher the accuracy improvements are obtained.

Numerical integration of the Teukolsky equation in the time domain

International Nuclear Information System (INIS)

Pazos-Avalos, Enrique; Lousto, Carlos O.

2005-01-01

We present a fourth-order convergent (2+1)-dimensional, numerical formalism to solve the Teukolsky equation in the time domain. Our approach is first to rewrite the Teukolsky equation as a system of first-order differential equations. In this way we get a system that has the form of an advection equation. This is then used in combination with a series expansion of the solution in powers of time. To obtain a fourth-order scheme we kept terms up to fourth derivative in time and use the advectionlike system of differential equations to substitute the temporal derivatives by spatial derivatives. This scheme is applied to evolve gravitational perturbations in the Schwarzschild and Kerr backgrounds. Our numerical method proved to be stable and fourth-order convergent in r* and θ directions. The correct power-law tail, ∼1/t 2l+3 , for general initial data, and ∼1/t 2l+4 , for time-symmetric data, was found in our runs. We noted that it is crucial to resolve accurately the angular dependence of the mode at late times in order to obtain these values of the exponents in the power-law decay. In other cases, when the decay was too fast and round-off error was reached before a tail was developed, then the quasinormal modes frequencies provided a test to determine the validity of our code

Analysis of the dual phase lag bio-heat transfer equation with constant and time-dependent heat flux conditions on skin surface

Directory of Open Access Journals (Sweden)

Ziaei Poor Hamed

2016-01-01

Full Text Available This article focuses on temperature response of skin tissue due to time-dependent surface heat fluxes. Analytical solution is constructed for DPL bio-heat transfer equation with constant, periodic and pulse train heat flux conditions on skin surface. Separation of variables and Duhamel’s theorem for a skin tissue as a finite domain are employed. The transient temperature responses for constant and time-dependent boundary conditions are obtained and discussed. The results show that there is major discrepancy between the predicted temperature of parabolic (Pennes bio-heat transfer, hyperbolic (thermal wave and DPL bio-heat transfer models when high heat flux accidents on the skin surface with a short duration or propagation speed of thermal wave is finite. The results illustrate that the DPL model reduces to the hyperbolic model when τT approaches zero and the classic Fourier model when both thermal relaxations approach zero. However for τq = τT the DPL model anticipates different temperature distribution with that predicted by the Pennes model. Such discrepancy is due to the blood perfusion term in energy equation. It is in contrast to results from the literature for pure conduction material, where the DPL model approaches the Fourier heat conduction model when τq = τT . The burn injury is also investigated.

Simulation of time-dependent free-surface Navier-Stokes flows

International Nuclear Information System (INIS)

Muldowney, G.P.

1989-01-01

Two numerical methods for simulation of time-dependent free-surface Navier-Stokes flows are developed. Both techniques are based on semi-implicit time advancement of the momentum equations, integral formulation of the spatial problem at each timestep, and spectral-element discretization to solve the resulting integral equation. Central to each algorithm is a boundary-specific solution step which permits the spatial treatment in two dimensions to be performed in O(N 3 ) operations per timestep despite the presence of deforming geometry. The first approach is a domain-integral formulation involving integrals over the entire flow domain of kernel functions which arise in time-differencing the Navier-Stokes equations. The second is a particular-solution formulation which replaces domain integration with an iterative scheme to generate particular velocity and pressure fields on individual elements, followed by a patching step to produce a particular solution continuous over the full domain. Two of the most difficult aspects of viscous free-surface flow simulations, namely time-dependent geometry and nontrivial boundary conditions, are well accommodated by these integral equation techniques. In addition the methods offer spectral accuracy in space and admit arbitrarily high-order discretization in time. For large-scale computations and/or long-term time advancement the domain-integral algorithm must be executed on a supercomputer to deliver results in reasonable processing time. A detailed simulation of gas liquid flow with full resolution of the free phase boundary requires approximately five CPU hours at 80 megaflops

Effect of the time spent by the photon in the absorbed state on the time-dependent transfer of radiation

International Nuclear Information System (INIS)

Rao, D.M.; Rangarajan, K.E.; Peraiah, A.

1990-01-01

The time-dependent transfer equation is derived for a two-level atomic model which takes both bound-bound and bound-free transitions into account. A numerical scheme is proposed for solving the monochromatic time-dependent transfer equation when the time spent by the photon in the absorbed state is significant. The method can be easily extended to solve the problem of time-dependent line formation of the bound-free continuum. It is used here to study three types of boundary conditions of the incident radiation incident on a scattering atmosphere. The quantitative results show that the relaxation of the radiation field depends on the optical depth of the medium and on the ray's angle of emergence. 21 refs

The recovery of a time-dependent point source in a linear transport equation: application to surface water pollution

International Nuclear Information System (INIS)

Hamdi, Adel

2009-01-01

The aim of this paper is to localize the position of a point source and recover the history of its time-dependent intensity function that is both unknown and constitutes the right-hand side of a 1D linear transport equation. Assuming that the source intensity function vanishes before reaching the final control time, we prove that recording the state with respect to the time at two observation points framing the source region leads to the identification of the source position and the recovery of its intensity function in a unique manner. Note that at least one of the two observation points should be strategic. We establish an identification method that determines quasi-explicitly the source position and transforms the task of recovering its intensity function into solving directly a well-conditioned linear system. Some numerical experiments done on a variant of the water pollution BOD model are presented

almaBTE : A solver of the space-time dependent Boltzmann transport equation for phonons in structured materials

Science.gov (United States)

Carrete, Jesús; Vermeersch, Bjorn; Katre, Ankita; van Roekeghem, Ambroise; Wang, Tao; Madsen, Georg K. H.; Mingo, Natalio

2017-11-01

almaBTE is a software package that solves the space- and time-dependent Boltzmann transport equation for phonons, using only ab-initio calculated quantities as inputs. The program can predictively tackle phonon transport in bulk crystals and alloys, thin films, superlattices, and multiscale structures with size features in the nm- μm range. Among many other quantities, the program can output thermal conductances and effective thermal conductivities, space-resolved average temperature profiles, and heat-current distributions resolved in frequency and space. Its first-principles character makes almaBTE especially well suited to investigate novel materials and structures. This article gives an overview of the program structure and presents illustrative examples for some of its uses. PROGRAM SUMMARY Program Title:almaBTE Program Files doi:http://dx.doi.org/10.17632/8tfzwgtp73.1 Licensing provisions: Apache License, version 2.0 Programming language: C++ External routines/libraries: BOOST, MPI, Eigen, HDF5, spglib Nature of problem: Calculation of temperature profiles, thermal flux distributions and effective thermal conductivities in structured systems where heat is carried by phonons Solution method: Solution of linearized phonon Boltzmann transport equation, Variance-reduced Monte Carlo

TIMEX: a time-dependent explicit discrete ordinates program for the solution of multigroup transport equations with delayed neutrons

International Nuclear Information System (INIS)

Hill, T.R.; Reed, W.H.

1976-01-01

TIMEX solves the time-dependent, one-dimensional multigroup transport equation with delayed neutrons in plane, cylindrical, spherical, and two-angle plane geometries. Both regular and adjoint, inhomogeneous and homogeneous problems subject to vacuum, reflective, periodic, white, albedo or inhomogeneous boundary flux conditions are solved. General anisotropic scattering is allowed and anisotropic inhomogeneous sources are permitted. The discrete ordinates approximation for the angular variable is used with the diamond (central) difference approximation for the angular extrapolation in curved geometries. A linear discontinuous finite element representation for the angular flux in each spatial mesh cell is used. The time variable is differenced by an explicit technique that is unconditionally stable so that arbitrarily large time steps can be taken. Because no iteration is performed the method is exceptionally fast in terms of computing time per time step. Two acceleration methods, exponential extrapolation and rebalance, are utilized to improve the accuracy of the time differencing scheme. Variable dimensioning is used so that any combination of problem parameters leading to a container array less than MAXCOR can be accommodated. The running time for TIMEX is highly problem-dependent, but varies almost linearly with the total number of unknowns and time steps. Provision is made for creation of standard interface output files for angular fluxes and angle-integrated fluxes. Five interface units (use of interface units is optional), five output units, and two system input/output units are required. A large bulk memory is desirable, but may be replaced by disk, drum, or tape storage. 13 tables, 9 figures

A maximum principle for time dependent transport in systems with voids

International Nuclear Information System (INIS)

Schofield, S.L.; Ackroyd, R.T.

1996-01-01

A maximum principle is developed for the first-order time dependent Boltzmann equation. The maximum principle is a generalization of Schofield's κ(θ) principle for the first-order steady state Boltzmann equation, and provides a treatment of time dependent transport in systems with void regions. The formulation comprises a direct least-squares minimization allied with a suitable choice of bilinear functional, and gives rise to a maximum principle whose functional is free of terms that have previously led to difficulties in treating void regions. (Author)

Exact time-dependent exchange-correlation potentials for strong-field electron dynamics

International Nuclear Information System (INIS)

Lein, Manfred; Kuemmel, Stephan

2005-01-01

By solving the time-dependent Schroedinger equation and inverting the time-dependent Kohn-Sham scheme we obtain the exact time-dependent exchange-correlation potential of density-functional theory for the strong-field dynamics of a correlated system. We demonstrate that essential features of the exact exchange-correlation potential can be related to derivative discontinuities in stationary density-functional theory. Incorporating the discontinuity in a time-dependent density-functional calculation greatly improves the description of the ionization process

The time dependent Hartree-Fock-theory for collective nuclear motions

International Nuclear Information System (INIS)

Goeke, K.

1976-11-01

The time-dependent Hartree-Fock theory (TDHF) approximately solves the Schroedinger equation by a variational method in the space of the time-dependent Slater determinants. As the TDHF wave function, similar to the exact solution has the property of being determined completely for all times by the nucleon-nucleon interaction and by assuming initial conditions. TDHF is expected to describe collective motion of nuclei with large amplitudes, too. The subject of this paper is to formulate the TDHF theory and its adiabatic limiting case (ATDHF) suited for setting up a collective Schroedinger equation, to investigate the relations with other theories, and to show the applicability for solving practical problems. (orig./WL) [de

FAST: a three-dimensional time-dependent FEL simulation code

International Nuclear Information System (INIS)

Saldin, E.L.; Schneidmiller, E.A.; Yurkov, M.V.

1999-01-01

In this report we briefly describe the three-dimensional, time-dependent FEL simulation code FAST. The equations of motion of the particles and Maxwell's equations are solved simultaneously taking into account the slippage effect. Radiation fields are calculated using an integral solution of Maxwell's equations. A special technique has been developed for fast calculations of the radiation field, drastically reducing the required CPU time. As a result, the developed code allows one to use a personal computer for time-dependent simulations. The code allows one to simulate the radiation from the electron bunch of any transverse and longitudinal bunch shape; to simulate simultaneously an external seed with superimposed noise in the electron beam; to take into account energy spread in the electron beam and the space charge fields; and to simulate a high-gain, high-efficiency FEL amplifier with a tapered undulator. It is important to note that there are no significant memory limitations in the developed code and an electron bunch of any length can be simulated

QUIPS: Time-dependent properties of quasi-invariant self-gravitating polytropes

International Nuclear Information System (INIS)

Munier, A.; Feix, M.R.

1983-01-01

Quasi-invariance, a method based on group tranformations, is used to obtain time-dependent solutions for the expansion and/or contraction of a self-gravitating sphere of perfect gas with polytopic index n. Quasi-invariance transforms the equations of hydrodynamics into ''dual equations'' exhibiting extra terms such as a friction, a mass source or sink term, and a centripetal/centrifugal force. The search for stationary solutions in this ''dual space'' leads to a new class of time-dependent solutions, the QUIP (for Quasi-invariant polytrope), which generalizes Emden's static model and introduces a characteristic frequency a related to Jean's frequency. The second order differential equation describing the solution is integrated numerically. A critical point is seen always to exist for nnot =3. Solutions corresponding in the ''dual space'' to a time-dependent generalization of Eddington's standard model (n = 3) are discussed. These solutions conserve both the total mass and the energy. A transition between closed and open structures is seen to take place at a particular frequency a/sub c/. For n = 3, no critical point arises in the ''dual space'' due to the self-similar motion of the fluid. A new time-dependent mass-radius relation and a generalized Betti-Ritter relation are obtained. Conclusions about the existence of a minimum Q-factor are presented

Exact solutions of Feinberg–Horodecki equation for time-dependent ...

Indian Academy of Sciences (India)

proposed to find the exact solutions of the Feinberg–Horodecki equation for the .... is a polynomial of degree n and satisfies the Rodrigues relation [24] ... The discriminant of expression (15) under the square root has to be zero, so that the.

Numerical studies of time-independent and time-dependent scattering by several elliptical cylinders

Science.gov (United States)

Nigsch, Martin

2007-07-01

A numerical solution to the problem of time-dependent scattering by an array of elliptical cylinders with parallel axes is presented. The solution is an exact one, based on the separation-of-variables technique in the elliptical coordinate system, the addition theorem for Mathieu functions, and numerical integration. Time-independent solutions are described by a system of linear equations of infinite order which are truncated for numerical computations. Time-dependent solutions are obtained by numerical integration involving a large number of these solutions. First results of a software package generating these solutions are presented: wave propagation around three impenetrable elliptical scatterers. As far as we know, this method described has never been used for time-dependent multiple scattering.

Sensitivity and uncertainty analysis for functionals of the time-dependent nuclide density field

International Nuclear Information System (INIS)

Williams, M.L.; Weisbin, C.R.

1978-04-01

An approach to extend the present ORNL sensitivity program to include functionals of the time-dependent nuclide density field is developed. An adjoint equation for the nuclide field was derived previously by using generalized perturbation theory; the present derivation makes use of a variational principle and results in the same equation. The physical significance of this equation is discussed and compared to that of the time-dependent neutron adjoint equation. Computational requirements for determining sensitivity profiles and uncertainties for functionals of the time-dependent nuclide density vector are developed within the framework of the existing FORSS system; in this way the current capability is significantly extended. The development, testing, and use of an adjoint version of the ORIGEN isotope generation and depletion code are documented. Finally, a sample calculation is given which estimates the uncertainty in the plutonium inventory at shutdown of a PWR due to assumed uncertainties in uranium and plutonium cross sections. 8 figures, 4 tables

Two-Dimensional Space-Time Dependent Multi-group Diffusion Equation with SLOR Method

International Nuclear Information System (INIS)

Yulianti, Y.; Su'ud, Z.; Waris, A.; Khotimah, S. N.

2010-01-01

The research of two-dimensional space-time diffusion equations with SLOR (Successive-Line Over Relaxation) has been done. SLOR method is chosen because this method is one of iterative methods that does not required to defined whole element matrix. The research is divided in two cases, homogeneous case and heterogeneous case. Homogeneous case has been inserted by step reactivity. Heterogeneous case has been inserted by step reactivity and ramp reactivity. In general, the results of simulations are agreement, even in some points there are differences.

Integrable time-dependent Hamiltonians, solvable Landau-Zener models and Gaudin magnets

Science.gov (United States)

Yuzbashyan, Emil A.

2018-05-01

We solve the non-stationary Schrödinger equation for several time-dependent Hamiltonians, such as the BCS Hamiltonian with an interaction strength inversely proportional to time, periodically driven BCS and linearly driven inhomogeneous Dicke models as well as various multi-level Landau-Zener tunneling models. The latter are Demkov-Osherov, bow-tie, and generalized bow-tie models. We show that these Landau-Zener problems and their certain interacting many-body generalizations map to Gaudin magnets in a magnetic field. Moreover, we demonstrate that the time-dependent Schrödinger equation for the above models has a similar structure and is integrable with a similar technique as Knizhnik-Zamolodchikov equations. We also discuss applications of our results to the problem of molecular production in an atomic Fermi gas swept through a Feshbach resonance and to the evaluation of the Landau-Zener transition probabilities.

Surface harmonics method for two-dimensional time-dependent neutron transport problems of square-lattice nuclear reactors

Energy Technology Data Exchange (ETDEWEB)

Boyarinov, V. F.; Kondrushin, A. E.; Fomichenko, P. A. [National Research Centre Kurchatov Institute, Kurchatov Sq. 1, Moscow (Russian Federation)

2013-07-01

Time-dependent equations of the Surface Harmonics Method (SHM) have been derived from the time-dependent neutron transport equation with explicit representation of delayed neutrons for solving the two-dimensional time-dependent problems. These equations have been realized in the SUHAM-TD code. The TWIGL benchmark problem has been used for verification of the SUHAM-TD code. The results of the study showed that computational costs required to achieve necessary accuracy of the solution can be an order of magnitude less than with the use of the conventional finite difference method (FDM). (authors)

Finite-difference solution of the space-angle-lethargy-dependent slowing-down transport equation

Energy Technology Data Exchange (ETDEWEB)

Matausek, M V [Boris Kidric Vinca Institute of Nuclear Sciences, Vinca, Belgrade (Yugoslavia)

1972-07-01

A procedure has been developed for solving the slowing-down transport equation for a cylindrically symmetric reactor system. The anisotropy of the resonance neutron flux is treated by the spherical harmonics formalism, which reduces the space-angle-Iethargy-dependent transport equation to a matrix integro-differential equation in space and lethargy. Replacing further the lethargy transfer integral by a finite-difference form, a set of matrix ordinary differential equations is obtained, with lethargy-and space dependent coefficients. If the lethargy pivotal points are chosen dense enough so that the difference correction term can be ignored, this set assumes a lower block triangular form and can be solved directly by forward block substitution. As in each step of the finite-difference procedure a boundary value problem has to be solved for a non-homogeneous system of ordinary differential equations with space-dependent coefficients, application of any standard numerical procedure, for example, the finite-difference method or the method of adjoint equations, is too cumbersome and would make the whole procedure practically inapplicable. A simple and efficient approximation is proposed here, allowing analytical solution for the space dependence of the spherical-harmonics flux moments, and hence the derivation of the recurrence relations between the flux moments at successive lethargy pivotal points. According to the procedure indicated above a computer code has been developed for the CDC -3600 computer, which uses the KEDAK nuclear data file. The space and lethargy distribution of the resonance neutrons can be computed in such a detailed fashion as the neutron cross-sections are known for the reactor materials considered. The computing time is relatively short so that the code can be efficiently used, either autonomously, or as part of some complex modular scheme. Typical results will be presented and discussed in order to prove and illustrate the applicability of the

A comparison of the Method of Lines to finite difference techniques in solving time-dependent partial differential equations. [with applications to Burger equation and stream function-vorticity problem

Science.gov (United States)

Kurtz, L. A.; Smith, R. E.; Parks, C. L.; Boney, L. R.

1978-01-01

Steady state solutions to two time dependent partial differential systems have been obtained by the Method of Lines (MOL) and compared to those obtained by efficient standard finite difference methods: (1) Burger's equation over a finite space domain by a forward time central space explicit method, and (2) the stream function - vorticity form of viscous incompressible fluid flow in a square cavity by an alternating direction implicit (ADI) method. The standard techniques were far more computationally efficient when applicable. In the second example, converged solutions at very high Reynolds numbers were obtained by MOL, whereas solution by ADI was either unattainable or impractical. With regard to 'set up' time, solution by MOL is an attractive alternative to techniques with complicated algorithms, as much of the programming difficulty is eliminated.

The roles of time and displacement in velocity-dependent volumetric strain of fault zones

Science.gov (United States)

Beeler, N.M.; Tullis, T.E.

1997-01-01

The relationship between measured friction??A and volumetric strain during frictional sliding was determined using a rate and state variable dependent friction constitutive equation, a common work balance relating friction and volume change, and two types of experimental faults: initially bare surfaces of Westerly granite and rock surfaces separated by a 1 mm layer of derivative of fault normal displacement with respect shear displacement, d??n ld??s. An implication of this relationship is that the rate dependence of d??n ld??s contributes to the rate dependence of ??A. Experiments show changes in sliding velocity lead to changes in both fault strength and volume. Analysis of data with the rate and state equations combined with the work balance relationship preclude the conventional interpretation of the direct effect in the rate and state variable constitutive equations. Consideration of a model bare surface fault consisting of an undeformable indentor sliding on a deformable surface reveals a serious flaw in the work balance relationship if volume change is time-dependent. For the model, at zero slip rate indentation creep under the normal load leads to time-dependent strengthening of the fault surface but, according to the work balance relationship, no work is done because compaction or dilatancy can only be induced by shearing. Additional tests on initially bare surfaces and gouges show that fault normal strain in experiments is time-dependent, consistent with the model. This time-dependent fault normal strain, which is not accounted for in the work balance relationship, explains the inconsistency between the constitutive equations and the work balance. For initially bare surface faults, all rate dependence of volume change is due to time dependence. Similar results are found for gouge. We conclude that ??A reflects the frictional resistance that results in shear heating, and no correction needs to be made for the volume changes. The result that time-dependent

Coupled electron and atomic kinetics through the solution of the Boltzmann equation for generating time-dependent X-ray spectra

International Nuclear Information System (INIS)

Sherrill, M.E.; Abdallah, J. Jr.; Csanak, G.; Kilcrease, D.P.; Dodd, E.S.; Fukuda, Y.; Akahane, Y.; Aoyama, M.; Inoue, N.; Ueda, H.; Yamakawa, K.; Faenov, A.Ya.; Magunov, A.I.; Pikuz, T.A.; Skobelev, I.Yu.

2006-01-01

In this work, we present a model that solves self-consistently the electron and atomic kinetics to characterize highly non-equilibrium plasmas, in particular for those systems where both the electron distribution function is far from Maxwellian and the evolution of the ion level populations are dominated by time-dependent atomic kinetics. In this model, level populations are obtained from a detailed collisional-radiative model where collision rates are computed from a time varying electron distribution function obtained from the solution of the zero-dimensional Boltzmann equation. The Boltzmann collision term includes the effects of electron-electron collisions, electron collisional ionization, excitation and de-excitation. An application for He α spectra from a short pulse laser irradiated argon cluster target will be shown to illustrate the results of our model

Coupled electron and atomic kinetics through the solution of the Boltzmann equation for generating time-dependent X-ray spectra

Energy Technology Data Exchange (ETDEWEB)

Sherrill, M.E. [Los Alamos National Laboratory, T-4, Los Alamos, NM 87545 (United States)]. E-mail: manolo@t4.lanl.gov; Abdallah, J. Jr. [Los Alamos National Laboratory, T-4, Los Alamos, NM 87545 (United States); Csanak, G. [Los Alamos National Laboratory, T-4, Los Alamos, NM 87545 (United States); Kilcrease, D.P. [Los Alamos National Laboratory, T-4, Los Alamos, NM 87545 (United States); Dodd, E.S. [Los Alamos National Laboratory, X-1, Los Alamos, NM 87545 (United States); Fukuda, Y. [Advanced Photon Research Center, JAERI, Kyoto 619-0215 (Japan); Akahane, Y. [Advanced Photon Research Center, JAERI, Kyoto 619-0215 (Japan); Aoyama, M. [Advanced Photon Research Center, JAERI, Kyoto 619-0215 (Japan); Inoue, N. [Advanced Photon Research Center, JAERI, Kyoto 619-0215 (Japan); Ueda, H. [Advanced Photon Research Center, JAERI, Kyoto 619-0215 (Japan); Yamakawa, K. [Advanced Photon Research Center, JAERI, Kyoto 619-0215 (Japan); Faenov, A.Ya. [Multicharged Ions Spectra Data Center of VNIIFTRI, Mendeleevo, Moscow Region 141570 (Russian Federation); Magunov, A.I. [Multicharged Ions Spectra Data Center of VNIIFTRI, Mendeleevo, Moscow Region 141570 (Russian Federation); Pikuz, T.A. [Multicharged Ions Spectra Data Center of VNIIFTRI, Mendeleevo, Moscow Region 141570 (Russian Federation); Skobelev, I.Yu. [Multicharged Ions Spectra Data Center of VNIIFTRI, Mendeleevo, Moscow Region 141570 (Russian Federation)

2006-05-15

In this work, we present a model that solves self-consistently the electron and atomic kinetics to characterize highly non-equilibrium plasmas, in particular for those systems where both the electron distribution function is far from Maxwellian and the evolution of the ion level populations are dominated by time-dependent atomic kinetics. In this model, level populations are obtained from a detailed collisional-radiative model where collision rates are computed from a time varying electron distribution function obtained from the solution of the zero-dimensional Boltzmann equation. The Boltzmann collision term includes the effects of electron-electron collisions, electron collisional ionization, excitation and de-excitation. An application for He{sub {alpha}} spectra from a short pulse laser irradiated argon cluster target will be shown to illustrate the results of our model.

Stability analysis of implicit time discretizations for the Compton-scattering Fokker-Planck equation

International Nuclear Information System (INIS)

Densmore, Jeffery D.; Warsa, James S.; Lowrie, Robert B.; Morel, Jim E.

2009-01-01

The Fokker-Planck equation is a widely used approximation for modeling the Compton scattering of photons in high energy density applications. In this paper, we perform a stability analysis of three implicit time discretizations for the Compton-Scattering Fokker-Planck equation. Specifically, we examine (i) a Semi-Implicit (SI) scheme that employs backward-Euler differencing but evaluates temperature-dependent coefficients at their beginning-of-time-step values, (ii) a Fully Implicit (FI) discretization that instead evaluates temperature-dependent coefficients at their end-of-time-step values, and (iii) a Linearized Implicit (LI) scheme, which is developed by linearizing the temperature dependence of the FI discretization within each time step. Our stability analysis shows that the FI and LI schemes are unconditionally stable and cannot generate oscillatory solutions regardless of time-step size, whereas the SI discretization can suffer from instabilities and nonphysical oscillations for sufficiently large time steps. With the results of this analysis, we present time-step limits for the SI scheme that prevent undesirable behavior. We test the validity of our stability analysis and time-step limits with a set of numerical examples.

Stability analysis of implicit time discretizations for the Compton-scattering Fokker-Planck equation

Energy Technology Data Exchange (ETDEWEB)

Densmore, Jeffery D [Los Alamos National Laboratory; Warsa, James S [Los Alamos National Laboratory; Lowrie, Robert B [Los Alamos National Laboratory; Morel, Jim E [TEXAS A& M UNIV

2008-01-01

The Fokker-Planck equation is a widely used approximation for modeling the Compton scattering of photons in high energy density applications. In this paper, we perform a stability analysis of three implicit time discretizations for the Compton-Scattering Fokker-Planck equation. Specifically, we examine (i) a Semi-Implicit (SI) scheme that employs backward-Euler differencing but evaluates temperature-dependent coefficients at their beginning-of-time-step values, (ii) a Fully Implicit (FI) discretization that instead evaluates temperature-dependent coefficients at their end-of-time-step values, and (iii) a Linearized Implicit (LI) scheme, which is developed by linearizing the temperature dependence of the FI discretization within each time step. Our stability analysis shows that the FI and LI schemes are unconditionally stable and cannot generate oscillatory solutions regardless of time-step size, whereas the SI discretization can suffer from instabilities and nonphysical oscillations for sufficiently large time steps. With the results of this analysis, we present time-step limits for the SI scheme that prevent undesirable behavior. We test the validity of our stability analysis and time-step limits with a set of numerical examples.

Stability analysis of implicit time discretizations for the Compton-scattering Fokker-Planck equation

Science.gov (United States)

Densmore, Jeffery D.; Warsa, James S.; Lowrie, Robert B.; Morel, Jim E.

2009-09-01

The Fokker-Planck equation is a widely used approximation for modeling the Compton scattering of photons in high energy density applications. In this paper, we perform a stability analysis of three implicit time discretizations for the Compton-Scattering Fokker-Planck equation. Specifically, we examine (i) a Semi-Implicit (SI) scheme that employs backward-Euler differencing but evaluates temperature-dependent coefficients at their beginning-of-time-step values, (ii) a Fully Implicit (FI) discretization that instead evaluates temperature-dependent coefficients at their end-of-time-step values, and (iii) a Linearized Implicit (LI) scheme, which is developed by linearizing the temperature dependence of the FI discretization within each time step. Our stability analysis shows that the FI and LI schemes are unconditionally stable and cannot generate oscillatory solutions regardless of time-step size, whereas the SI discretization can suffer from instabilities and nonphysical oscillations for sufficiently large time steps. With the results of this analysis, we present time-step limits for the SI scheme that prevent undesirable behavior. We test the validity of our stability analysis and time-step limits with a set of numerical examples.

Stability of time-dependent particle-like solutions of some wave equations

International Nuclear Information System (INIS)

Voronov, N.A.

1978-01-01

The proof of the nonstability of the one-dimensional periodical localized solutions of the equation with a spontaneously broken symmetry is given. The stability of the one-dimensional oscillating solutions of the sine-Gordon equation was also considered with regard to such perturbations. As it was expected these solutions proved to be stable

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

OpenMP GNU and Intel Fortran programs for solving the time-dependent Gross-Pitaevskii equation

Science.gov (United States)

Young-S., Luis E.; Muruganandam, Paulsamy; Adhikari, Sadhan K.; Lončar, Vladimir; Vudragović, Dušan; Balaž, Antun

2017-11-01

We present Open Multi-Processing (OpenMP) version of Fortran 90 programs for solving the Gross-Pitaevskii (GP) equation for a Bose-Einstein condensate in one, two, and three spatial dimensions, optimized for use with GNU and Intel compilers. We use the split-step Crank-Nicolson algorithm for imaginary- and real-time propagation, which enables efficient calculation of stationary and non-stationary solutions, respectively. The present OpenMP programs are designed for computers with multi-core processors and optimized for compiling with both commercially-licensed Intel Fortran and popular free open-source GNU Fortran compiler. The programs are easy to use and are elaborated with helpful comments for the users. All input parameters are listed at the beginning of each program. Different output files provide physical quantities such as energy, chemical potential, root-mean-square sizes, densities, etc. We also present speedup test results for new versions of the programs. Program files doi:http://dx.doi.org/10.17632/y8zk3jgn84.2 Licensing provisions: Apache License 2.0 Programming language: OpenMP GNU and Intel Fortran 90. Computer: Any multi-core personal computer or workstation with the appropriate OpenMP-capable Fortran compiler installed. Number of processors used: All available CPU cores on the executing computer. Journal reference of previous version: Comput. Phys. Commun. 180 (2009) 1888; ibid.204 (2016) 209. Does the new version supersede the previous version?: Not completely. It does supersede previous Fortran programs from both references above, but not OpenMP C programs from Comput. Phys. Commun. 204 (2016) 209. Nature of problem: The present Open Multi-Processing (OpenMP) Fortran programs, optimized for use with commercially-licensed Intel Fortran and free open-source GNU Fortran compilers, solve the time-dependent nonlinear partial differential (GP) equation for a trapped Bose-Einstein condensate in one (1d), two (2d), and three (3d) spatial dimensions for

A perfectly matched layer for the time-dependent wave equation in heterogeneous and layered media

KAUST Repository

Duru, Kenneth

2014-01-01

A mathematical analysis of the perfectly matched layer (PML) for the time-dependent wave equation in heterogeneous and layered media is presented. We prove the stability of the PML for discontinuous media with piecewise constant coefficients, and derive energy estimates for discontinuous media with piecewise smooth coefficients. We consider a computational setup consisting of smaller structured subdomains that are discretized using high order accurate finite difference operators for approximating spatial derivatives. The subdomains are then patched together into a global domain by a weak enforcement of interface conditions using penalties. In order to ensure the stability of the discrete PML, it is necessary to transform the interface conditions to include the auxiliary variables. In the discrete setting, the transformed interface conditions are crucial in deriving discrete energy estimates analogous to the continuous energy estimates, thus proving stability and convergence of the numerical method. Finally, we present numerical experiments demonstrating the stability of the PML in a layered medium and high order accuracy of the proposed interface conditions. © 2013 Elsevier Inc.

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.

Solving Algebraic Riccati Equation Real Time for Integrated Vehicle Dynamics Control

NARCIS (Netherlands)

Kunnappillil Madhusudhanan, A; Corno, M.; Bonsen, B.; Holweg, E.

2012-01-01

In this paper we present a comparison study of different computational methods to implement State Dependent Riccati Equation (SDRE) based control in real time for a vehicle dynamics control application. Vehicles are mechatronic systems with nonlinear dynamics. One of the promising nonlinear control

Time-dependent weak values and their intrinsic phases of evolution

International Nuclear Information System (INIS)

Parks, A D

2008-01-01

The equation of motion for a time-dependent weak value of a quantum-mechanical observable is known to contain a complex valued energy factor (the weak energy of evolution) that is defined by the dynamics of the pre-selected and post-selected states which specify the observable's weak value. In this paper, the mechanism responsible for the creation of this energy is identified and it is shown that the cumulative effect over time of this energy is manifested as dynamical phases and pure geometric phases (the intrinsic phases of evolution) which govern the evolution of the weak value during its measurement process. These phases are simply related to a Pancharatnam phase and Fubini-Study metric distance defined by the Hilbert space evolution of the associated pre-selected and post-selected states. A characterization of time-dependent weak value evolution as Pancharatnam phase angle rotations and Fubini-Study distance scalings of a vector in the Argand plane is discussed as an application of this relationship. The theory of weak values is also reviewed and simple 'gedanken experiments' are used to illustrate both the time-independent and the time-dependent versions of the theory. It is noted that the direct experimental observation of the weak energy of evolution would strongly support the time-symmetric paradigm of quantum mechanics and it is suggested that weak value equations of motion represent a new category of nonlocal equations of motion

Computational Diffusion Magnetic Resonance Imaging Based on Time-Dependent Bloch NMR Flow Equation and Bessel Functions.

Science.gov (United States)

Awojoyogbe, Bamidele O; Dada, Michael O; Onwu, Samuel O; Ige, Taofeeq A; Akinwande, Ninuola I

2016-04-01

Magnetic resonance imaging (MRI) uses a powerful magnetic field along with radio waves and a computer to produce highly detailed "slice-by-slice" pictures of virtually all internal structures of matter. The results enable physicians to examine parts of the body in minute detail and identify diseases in ways that are not possible with other techniques. For example, MRI is one of the few imaging tools that can see through bones, making it an excellent tool for examining the brain and other soft tissues. Pulsed-field gradient experiments provide a straightforward means of obtaining information on the translational motion of nuclear spins. However, the interpretation of the data is complicated by the effects of restricting geometries as in the case of most cancerous tissues and the mathematical concept required to account for this becomes very difficult. Most diffusion magnetic resonance techniques are based on the Stejskal-Tanner formulation usually derived from the Bloch-Torrey partial differential equation by including additional terms to accommodate the diffusion effect. Despite the early success of this technique, it has been shown that it has important limitations, the most of which occurs when there is orientation heterogeneity of the fibers in the voxel of interest (VOI). Overcoming this difficulty requires the specification of diffusion coefficients as function of spatial coordinate(s) and such a phenomenon is an indication of non-uniform compartmental conditions which can be analyzed accurately by solving the time-dependent Bloch NMR flow equation analytically. In this study, a mathematical formulation of magnetic resonance flow sequence in restricted geometry is developed based on a general second order partial differential equation derived directly from the fundamental Bloch NMR flow equations. The NMR signal is obtained completely in terms of NMR experimental parameters. The process is described based on Bessel functions and properties that can make it

Extended rate equations

International Nuclear Information System (INIS)

Shore, B.W.

1981-01-01

The equations of motion are discussed which describe time dependent population flows in an N-level system, reviewing the relationship between incoherent (rate) equations, coherent (Schrodinger) equations, and more general partially coherent (Bloch) equations. Approximations are discussed which replace the elaborate Bloch equations by simpler rate equations whose coefficients incorporate long-time consequences of coherence

Semiclassical approximation to time-dependent Hartree--Fock theory

International Nuclear Information System (INIS)

Dworzecka, M.; Poggioli, R.

1976-01-01

Working within a time-dependent Hartree-Fock framework, one develops a semiclassical approximation appropriate for large systems. It is demonstrated that the standard semiclassical approach, the Thomas-Fermi approximation, is inconsistent with Hartree-Fock theory when the basic two-body interaction is short-ranged (as in nuclear systems, for example). However, by introducing a simple extension of the Thomas-Fermi approximation, one overcomes this problem. One also discusses the infinite nuclear matter problem and point out that time-dependent Hartree-Fock theory yields collective modes of the zero sound variety instead of ordinary hydrodynamic (first) sound. One thus emphasizes that one should be extremely circumspect when attempting to cast the equations of motion of time-dependent Hartree-Fock theory into a hydrodynamic-like form

Anomalous transport in fluid field with random waiting time depending on the preceding jump length

International Nuclear Information System (INIS)

Zhang Hong; Li Guo-Hua

2016-01-01

Anomalous (or non-Fickian) transport behaviors of particles have been widely observed in complex porous media. To capture the energy-dependent characteristics of non-Fickian transport of a particle in flow fields, in the present paper a generalized continuous time random walk model whose waiting time probability distribution depends on the preceding jump length is introduced, and the corresponding master equation in Fourier–Laplace space for the distribution of particles is derived. As examples, two generalized advection-dispersion equations for Gaussian distribution and lévy flight with the probability density function of waiting time being quadratic dependent on the preceding jump length are obtained by applying the derived master equation. (paper)

FEM for time-fractional diffusion equations, novel optimal error analyses

OpenAIRE

Mustapha, Kassem

2016-01-01

A semidiscrete Galerkin finite element method applied to time-fractional diffusion equations with time-space dependent diffusivity on bounded convex spatial domains will be studied. The main focus is on achieving optimal error results with respect to both the convergence order of the approximate solution and the regularity of the initial data. By using novel energy arguments, for each fixed time $t$, optimal error bounds in the spatial $L^2$- and $H^1$-norms are derived for both cases: smooth...

The development of the time dependence of the nuclear EMP electric field

International Nuclear Information System (INIS)

Eng, C.

2009-01-01

The nuclear electromagnetic pulse (EMP) electric field calculated with the legacy code CHAP is compared with the field given by an integral solution of Maxwell's equations, also known as the Jefimenko equation, to aid our current understanding on the factors that affect the time dependence of the EMP. For a fair comparison the CHAP current density is used as a source in the Jefimenko equation. At first, the comparison is simplified by neglecting the conduction current and replacing the standard atmosphere with a constant density air slab. The simplicity of the resultant current density aids in determining the factors that affect the rise, peak and tail of the EMP electric field versus time. The three dimensional nature of the radiating source, i.e. sources off the line-of-sight, and the time dependence of the derivative of the current density with respect to time are found to play significant roles in shaping the EMP electric field time dependence. These results are found to hold even when the conduction current and the standard atmosphere are properly accounted for. Comparison of the CHAP electric field with the Jefimenko electric field offers a direct validation of the high-frequency/outgoing wave approximation.

Global existence and large time asymptotic behavior of strong solutions to the Cauchy problem of 2D density-dependent Navier–Stokes equations with vacuum

Science.gov (United States)

Lü, Boqiang; Shi, Xiaoding; Zhong, Xin

2018-06-01

We are concerned with the Cauchy problem of the two-dimensional (2D) nonhomogeneous incompressible Navier–Stokes equations with vacuum as far-field density. It is proved that if the initial density decays not too slow at infinity, the 2D Cauchy problem of the density-dependent Navier–Stokes equations on the whole space admits a unique global strong solution. Note that the initial data can be arbitrarily large and the initial density can contain vacuum states and even have compact support. Furthermore, we also obtain the large time decay rates of the spatial gradients of the velocity and the pressure, which are the same as those of the homogeneous case.

Continuous-time random walk as a guide to fractional Schroedinger equation

International Nuclear Information System (INIS)

Lenzi, E. K.; Ribeiro, H. V.; Mukai, H.; Mendes, R. S.

2010-01-01

We argue that the continuous-time random walk approach may be a useful guide to extend the Schroedinger equation in order to incorporate nonlocal effects, avoiding the inconsistencies raised by Jeng et al. [J. Math. Phys. 51, 062102 (2010)]. As an application, we work out a free particle in a half space, obtaining the time dependent solution by considering an arbitrary initial condition.

Simulation of Near-Edge X-ray Absorption Fine Structure with Time-Dependent Equation-of-Motion Coupled-Cluster Theory.

Science.gov (United States)

Nascimento, Daniel R; DePrince, A Eugene

2017-07-06

An explicitly time-dependent (TD) approach to equation-of-motion (EOM) coupled-cluster theory with single and double excitations (CCSD) is implemented for simulating near-edge X-ray absorption fine structure in molecular systems. The TD-EOM-CCSD absorption line shape function is given by the Fourier transform of the CCSD dipole autocorrelation function. We represent this transform by its Padé approximant, which provides converged spectra in much shorter simulation times than are required by the Fourier form. The result is a powerful framework for the blackbox simulation of broadband absorption spectra. K-edge X-ray absorption spectra for carbon, nitrogen, and oxygen in several small molecules are obtained from the real part of the absorption line shape function and are compared with experiment. The computed and experimentally obtained spectra are in good agreement; the mean unsigned error in the predicted peak positions is only 1.2 eV. We also explore the spectral signatures of protonation in these molecules.

From interatomic interaction potentials via Einstein field equation techniques to time dependent contact mechanics

International Nuclear Information System (INIS)

Schwarzer, N

2014-01-01

In order to understand the principle differences between rheological or simple stress tests like the uniaxial tensile test to contact mechanical tests and the differences between quasistatic contact experiments and oscillatory ones, this study resorts to effective first principles. This study will show how relatively simple models simulating bond interactions in solids using effective potentials like Lennard-Jones and Morse can be used to investigate the effect of time dependent stress-induced softening or stiffening of these solids. The usefulness of the current study is in the possibility of deriving relatively simple dependences of the bulk-modulus B on time, shear and pressure P with time t. In cases where it is possible to describe, or at least partially describe a material by Lennard-Jones potential approaches, the above- mentioned dependences are even completely free of microscopic material parameters. Instead of bond energies and length, only specific integral parameters like Young’s modulus and Poisson’s ratio are required. However, in the case of time dependent (viscose) material behavior the parameters are not constants anymore. They themselves depend on time and the actual stress field, especially the shear field. A body completely consisting of so called standard linear solid interacting particles will then phenomenologically show a completely different and usually much more complicated mechanical behavior. The influence of the time dependent pressure-shear-induced Young’s modulus change is discussed with respect to mechanical contact experiments and their analysis in the case of viscose materials. (papers)

Dissipative time-dependent quantum transport theory.

Science.gov (United States)

Zhang, Yu; Yam, Chi Yung; Chen, GuanHua

2013-04-28

A dissipative time-dependent quantum transport theory is developed to treat the transient current through molecular or nanoscopic devices in presence of electron-phonon interaction. The dissipation via phonon is taken into account by introducing a self-energy for the electron-phonon coupling in addition to the self-energy caused by the electrodes. Based on this, a numerical method is proposed. For practical implementation, the lowest order expansion is employed for the weak electron-phonon coupling case and the wide-band limit approximation is adopted for device and electrodes coupling. The corresponding hierarchical equation of motion is derived, which leads to an efficient and accurate time-dependent treatment of inelastic effect on transport for the weak electron-phonon interaction. The resulting method is applied to a one-level model system and a gold wire described by tight-binding model to demonstrate its validity and the importance of electron-phonon interaction for the quantum transport. As it is based on the effective single-electron model, the method can be readily extended to time-dependent density functional theory.

Existence of solitary waves in dipolar quantum gases

KAUST Repository

Antonelli, Paolo; Sparber, Christof

2011-01-01

We study a nonlinear Schrdinger equation arising in the mean field description of dipolar quantum gases. Under the assumption of sufficiently strong dipolar interactions, the existence of standing waves, and hence solitons, is proved together with some of their properties. This gives a rigorous argument for the possible existence of solitary waves in BoseEinstein condensates, which originate solely due to the dipolar interaction between the particles. © 2010 Elsevier B.V. All rights reserved.

Existence of solitary waves in dipolar quantum gases

KAUST Repository

Antonelli, Paolo

2011-02-01

We study a nonlinear Schrdinger equation arising in the mean field description of dipolar quantum gases. Under the assumption of sufficiently strong dipolar interactions, the existence of standing waves, and hence solitons, is proved together with some of their properties. This gives a rigorous argument for the possible existence of solitary waves in BoseEinstein condensates, which originate solely due to the dipolar interaction between the particles. © 2010 Elsevier B.V. All rights reserved.

Monoenergetic time-dependent neutron transport in an infinite medium with time-varying cross sections

International Nuclear Information System (INIS)

Ganapol, B.D.

1987-01-01

For almost 20 yr, the main thrust of the author's research has been the generation of as many benchmark solutions to the time-dependent monoenergetic neutron transport equation as possible. The major motivation behind this effort has been to provide code developers with highly accurate numerical solutions to serve as standards in the assessment of numerical transport algorithms. In addition, these solutions provide excellent educational tools since the important physical features of neutron transport are still present even though the problems solved are idealized. A secondary motivation, though of equal importance, is the intellectual stimulation and understanding provided by the combination of the analytical, numerical, and computational techniques required to obtain these solutions. Therefore, to further the benchmark development, the added complication of time-dependent cross sections in the one-group transport equation is considered here

Nonlinear time-dependent simulation of helix traveling wave tubes

International Nuclear Information System (INIS)

Peng Wei-Feng; Yang Zhong-Hai; Hu Yu-Lu; Li Jian-Qing; Lu Qi-Ru; Li Bin

2011-01-01

A one-dimensional nonlinear time-dependent theory for helix traveling wave tubes is studied. A generalized electromagnetic field is applied to the expression of the radio frequency field. To simulate the variations of the high frequency structure, such as the pitch taper and the effect of harmonics, the spatial average over a wavelength is substituted by a time average over a wave period in the equation of the radio frequency field. Under this assumption, the space charge field of the electron beam can be treated by a space charge wave model along with the space charge coefficient. The effects of the radio frequency and the space charge fields on the electrons are presented by the equations of the electron energy and the electron phase. The time-dependent simulation is compared with the frequency-domain simulation for a helix TWT, which validates the availability of this theory. (interdisciplinary physics and related areas of science and technology)

Applications of alternating direction methods to the solution of the time-dependent heat conduction equation with source and in transients stage

International Nuclear Information System (INIS)

Gebrin, A.N.

1981-10-01

Various types and also some variants of alternating direction methods, (A.D.M.), were applied to the solution of the time-dependent heat conduction equation, with source, in region with axial symetry. The results shown that some of the variants perform consistently better than the Classical Cranck-Nicolson method. Having in mind a combination of accuracy, ability to support larg time steps and computational efficiency, the 'alternating direction explicit', (A.D.E.) method appears as the best choice, being the 'alternating direction checkerboard', (A.D.C), method the second best. Additional operations like the exponential transformation or the truncation pos-correction don't seem to be worth, excect for some special cases. (Author) [pt

Time evolution of the wave equation using rapid expansion method

KAUST Repository

Pestana, Reynam C.; Stoffa, Paul L.

2010-01-01

Forward modeling of seismic data and reverse time migration are based on the time evolution of wavefields. For the case of spatially varying velocity, we have worked on two approaches to evaluate the time evolution of seismic wavefields. An exact solution for the constant-velocity acoustic wave equation can be used to simulate the pressure response at any time. For a spatially varying velocity, a one-step method can be developed where no intermediate time responses are required. Using this approach, we have solved for the pressure response at intermediate times and have developed a recursive solution. The solution has a very high degree of accuracy and can be reduced to various finite-difference time-derivative methods, depending on the approximations used. Although the two approaches are closely related, each has advantages, depending on the problem being solved. © 2010 Society of Exploration Geophysicists.

Time evolution of the wave equation using rapid expansion method

KAUST Repository

Pestana, Reynam C.

2010-07-01

Forward modeling of seismic data and reverse time migration are based on the time evolution of wavefields. For the case of spatially varying velocity, we have worked on two approaches to evaluate the time evolution of seismic wavefields. An exact solution for the constant-velocity acoustic wave equation can be used to simulate the pressure response at any time. For a spatially varying velocity, a one-step method can be developed where no intermediate time responses are required. Using this approach, we have solved for the pressure response at intermediate times and have developed a recursive solution. The solution has a very high degree of accuracy and can be reduced to various finite-difference time-derivative methods, depending on the approximations used. Although the two approaches are closely related, each has advantages, depending on the problem being solved. © 2010 Society of Exploration Geophysicists.

Ordinary differential equation for local accumulation time.

Science.gov (United States)

Berezhkovskii, Alexander M

2011-08-21

Cell differentiation in a developing tissue is controlled by the concentration fields of signaling molecules called morphogens. Formation of these concentration fields can be described by the reaction-diffusion mechanism in which locally produced molecules diffuse through the patterned tissue and are degraded. The formation kinetics at a given point of the patterned tissue can be characterized by the local accumulation time, defined in terms of the local relaxation function. Here, we show that this time satisfies an ordinary differential equation. Using this equation one can straightforwardly determine the local accumulation time, i.e., without preliminary calculation of the relaxation function by solving the partial differential equation, as was done in previous studies. We derive this ordinary differential equation together with the accompanying boundary conditions and demonstrate that the earlier obtained results for the local accumulation time can be recovered by solving this equation. © 2011 American Institute of Physics

Time-dependent non-equilibrium dielectric response in QM/continuum approaches

Energy Technology Data Exchange (ETDEWEB)

Ding, Feizhi; Lingerfelt, David B.; Li, Xiaosong, E-mail: benedetta.mennucci@unipi.it, E-mail: li@chem.washington.edu [Department of Chemistry, University of Washington, Seattle, Washington 98195 (United States); Mennucci, Benedetta, E-mail: benedetta.mennucci@unipi.it, E-mail: li@chem.washington.edu [Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, 56126 Pisa (Italy)

2015-01-21

The Polarizable Continuum Models (PCMs) are some of the most inexpensive yet successful methods for including the effects of solvation in quantum-mechanical calculations of molecular systems. However, when applied to the electronic excitation process, these methods are restricted to dichotomously assuming either that the solvent has completely equilibrated with the excited solute charge density (infinite-time limit), or that it retains the configuration that was in equilibrium with the solute prior to excitation (zero-time limit). This renders the traditional PCMs inappropriate for resolving time-dependent solvent effects on non-equilibrium solute electron dynamics like those implicated in the instants following photoexcitation of a solvated molecular species. To extend the existing methods to this non-equilibrium regime, we herein derive and apply a new formalism for a general time-dependent continuum embedding method designed to be propagated alongside the solute’s electronic degrees of freedom in the time domain. Given the frequency-dependent dielectric constant of the solvent, an equation of motion for the dielectric polarization is derived within the PCM framework and numerically integrated simultaneously with the time-dependent Hartree fock/density functional theory equations. Results for small molecular systems show the anticipated dipole quenching and electronic state dephasing/relaxation resulting from out-of-phase charge fluctuations in the dielectric and embedded quantum system.

Excitonic effects in solids : time-dependent density functional theory versus the Bethe-Salpeter equation

International Nuclear Information System (INIS)

Sagmeister, S.

2009-01-01

The aim of this work is to compare two state-of-the-art methods for the investigation of excitonic effects in solids, namely Time-Dependent Density Functional Theory (TDDFT) and Many-Body Perturbation Theory (MBPT), for selected simple gap systems as well as semiconducting polymers. Within TDDFT, the linear response framework is used and the Dyson equation for the density-density response function is solved, whereas within MBPT, the Bethe-Salpeter equation (BSE) for the electron-hole correlation function is solved. The dielectric function is obtained as a last step. Both techniques take into account the excitonic effects caused by the interaction of electron-hole pairs. In the former these effects are included in the exchange-correlation (xc) kernel, whereas in the latter they are located in the interaction kernel of the BSE. Kohn-Sham single-particle wave functions obtained from Density Functional Theory within the linearized augmented planewave (LAPW) method are used to calculate all relevant quantities of the formalism. For the simple systems GaAs, Si and LiF are chosen. The role of several approximations to the xc kernel is studied and it is found that for GaAs and Si simple semi-empirical models provide a dielectric function in accordance with the BSE. For the case of LiF, being a system with a weak screening and a strongly bound exciton, only an xc kernel derived from MBPT yields reasonable results but still a slight discrepancy to the BSE is observed. Finally, the semiconducting polymers poly-acetylene and poly(phenylene-vinylene) (PPV) are studied. For both materials the concept of semi-empirical approximations to the xc kernel turns out to be ambiguous due to their low-dimensional character. In the case of poly-acetylene, the xc kernel derived from MBPT yields a dielectric function which is in close but not exact agreement with the one obtained from the BSE. (author) [de

An application of transverse momentum dependent evolution equations in QCD

International Nuclear Information System (INIS)

Ceccopieri, Federico A.; Trentadue, Luca

2008-01-01

The properties and behaviour of the solutions of the recently obtained k t -dependent QCD evolution equations are investigated. When used to reproduce transverse momentum spectra of hadrons in Semi-Inclusive DIS, an encouraging agreement with data is found. The present analysis also supports at the phenomenological level the factorization properties of the Semi-Inclusive DIS cross-sections in terms of k t -dependent distributions. Further improvements and possible developments of the proposed evolution equations are envisaged

Induced voltage due to time-dependent magnetisation textures

International Nuclear Information System (INIS)

Kudtarkar, Santosh Kumar; Dhadwal, Renu

2010-01-01

We determine the induced voltage generated by spatial and temporal magnetisation textures (inhomogeneities) in metallic ferromagnets due to the spin diffusion of non-equilibrium electrons. Using time dependent semi-classical theory as formulated in Zhang and Li and the drift-diffusion model of transport it is shown that the voltage generated depends critically on the difference in the diffusion constants of up and down spins. Including spin relaxation results in a crucial contribution to the induced voltage. We also show that the presence of magnetisation textures results in the modification of the conductivity of the system. As an illustration, we calculate the voltage generated due to a time dependent field driven helimagnet by solving the Landau-Lifshitz equation with Gilbert damping and explicitly calculate the dependence on the relaxation and damping parameters.

Existence results for impulsive evolution differential equations with state-dependent delay

OpenAIRE

Eduardo Hernandez M.; Rathinasamy Sakthivel; Sueli Tanaka Aki

2008-01-01

We study the existence of mild solution for impulsive evolution abstract differential equations with state-dependent delay. A concrete application to partial delayed differential equations is considered.

Watching excitons move: the time-dependent transition density matrix

Science.gov (United States)

Ullrich, Carsten

2012-02-01

Time-dependent density-functional theory allows one to calculate excitation energies and the associated transition densities in principle exactly. The transition density matrix (TDM) provides additional information on electron-hole localization and coherence of specific excitations of the many-body system. We have extended the TDM concept into the real-time domain in order to visualize the excited-state dynamics in conjugated molecules. The time-dependent TDM is defined as an implicit density functional, and can be approximately obtained from the time-dependent Kohn-Sham orbitals. The quality of this approximation is assessed in simple model systems. A computational scheme for real molecular systems is presented: the time-dependent Kohn-Sham equations are solved with the OCTOPUS code and the time-dependent Kohn-Sham TDM is calculated using a spatial partitioning scheme. The method is applied to show in real time how locally created electron-hole pairs spread out over neighboring conjugated molecular chains. The coupling mechanism, electron-hole coherence, and the possibility of charge separation are discussed.

Long-time behavior for suspension bridge equations with time delay

Science.gov (United States)

Park, Sun-Hye

2018-04-01

In this paper, we consider suspension bridge equations with time delay of the form u_{tt}(x,t) + Δ ^2 u (x,t) + k u^+ (x,t) + a_0 u_t (x,t) + a_1 u_t (x, t- τ ) + f(u(x,t)) = g(x). Many researchers have studied well-posedness, decay rates of energy, and existence of attractors for suspension bridge equations without delay effects. But, as far as we know, there is no work about suspension equations with time delay. In addition, there are not many studies on attractors for other delayed systems. Thus we first provide well-posedness for suspension equations with time delay. And then show the existence of global attractors and the finite dimensionality of the attractors by establishing energy functionals which are related to the norm of the phase space to our problem.

SYNTH-C, Steady-State and Time-Dependent 3-D Neutron Diffusion with Thermohydraulic Feedback

Energy Technology Data Exchange (ETDEWEB)

Brega, E [ENEL-CRTN, Bastioni di Porta Volta 10, Milan (Italy); Salina, E [A.R.S. Spa, Viale Maino 35, Milan (Italy)

1980-04-01

1 - Description of problem or function: SYNTH-C-STEADY and SYNTH-C- TRANS solve respectively the steady-state and time-dependent few- group neutron diffusion equations in three dimensions x,y,z in the presence of fuel temperature and thermal-hydraulic feedback. The neutron diffusion and delayed precursor equations are approximated by a space-time (z,t) synthesis method with axially discontinuous trial functions. Three thermal-hydraulic and fuel heat transfer models are available viz. COBRA-3C/MIT model, lumped parameter (WIGL) model and adiabatic fuel heat-up model. 2 - Method of solution: The steady-state and time-dependent synthesis equations are solved respectively by the Wielandt's power method and by the theta-difference method (in time), both coupled with a block factorization technique and double precision arithmetic. The thermal-hydraulic model equations are solved by fully implicit finite differences (WIGL) or explicit-implicit difference techniques with iterations (COBRA-EC/MIT). 3 - Restrictions on the complexity of the problem: Except for the few- group limitation, the programs have no other fixed limitation so the ability to run a problem depends only on the available computer storage.

Simulating time-dependent energy transfer between crossed laser beams in an expanding plasma

International Nuclear Information System (INIS)

Hittinger, J.A.F.; Dorr, M.R.; Berger, R.L.; Williams, E.A.

2005-01-01

A coupled mode system is derived to investigate a three-wave parametric instability leading to energy transfer between co-propagating laser beams crossing in a plasma flow. The model includes beams of finite width refracting in a prescribed transverse plasma flow with spatial and temporal gradients in velocity and density. The resulting paraxial light equations are discretized spatially with a Crank-Nicholson-type scheme, and these algebraic constraints are nonlinearly coupled with ordinary differential equations in time that describe the ion acoustic response. The entire nonlinear differential-algebraic system is solved using an adaptive, backward-differencing method coupled with Newton's method. A numerical study is conducted in two dimensions that compares the intensity gain of the fully time-dependent coupled mode system with the gain computed under the further assumption of a strongly damped ion acoustic response. The results demonstrate a time-dependent gain suppression when the beam diameter is commensurate with the velocity gradient scale length. The gain suppression is shown to depend on time-dependent beam refraction and is interpreted as a time-dependent frequency shift

Multiconfiguration hartree-fock theory for pseudorelativistic systems: The time-dependent case

KAUST Repository

Hajaiej, Hichem

2014-03-01

In [Setting and analysis of the multi-configuration time-dependent Hartree-Fock equations, Arch. Ration. Mech. Anal. 198 (2010) 273-330] the third author has studied in collaboration with Bardos, Catto and Mauser the nonrelativistic multiconfiguration time-dependent Hartree-Fock system of equations arising in the modeling of molecular dynamics. In this paper, we extend the previous work to the case of pseudorelativistic atoms. We show the existence and the uniqueness of global-in-time solution to the underlying system under technical assumptions on the energy of the initial data and the charge of the nucleus. Moreover, we prove that the result can be extended to the case of neutron stars when the number of electrons is less than a critical number N cr. © 2014 World Scientific Publishing Company.

Time-dependent theoretical treatments of the dynamics of electrons and nuclei in molecular systems

International Nuclear Information System (INIS)

Deumens, E.; Diz, A.; Longo, R.; Oehrn, Y.

1994-01-01

An overview is presented of methods for time-dependent treatments of molecules as systems of electrons and nuclei. The theoretical details of these methods are reviewed and contrasted in the light of a recently developed time-dependent method called electron-nuclear dynamics. Electron-nuclear dynamics (END) is a formulation of the complete dynamics of electrons and nuclei of a molecular system that eliminates the necessity of constructing potential-energy surfaces. Because of its general formulation, it encompasses many aspects found in other formulations and can serve as a didactic device for clarifying many of the principles and approximations relevant in time-dependent treatments of molecular systems. The END equations are derived from the time-dependent variational principle applied to a chosen family of efficiently parametrized approximate state vectors. A detailed analysis of the END equations is given for the case of a single-determinantal state for the electrons and a classical treatment of the nuclei. The approach leads to a simple formulation of the fully nonlinear time-dependent Hartree-Fock theory including nuclear dynamics. The nonlinear END equations with the ab initio Coulomb Hamiltonian have been implemented at this level of theory in a computer program, ENDyne, and have been shown feasible for the study of small molecular systems. Implementation of the Austin Model 1 semiempirical Hamiltonian is discussed as a route to large molecular systems. The linearized END equations at this level of theory are shown to lead to the random-phase approximation for the coupled system of electrons and nuclei. The qualitative features of the general nonlinear solution are analyzed using the results of the linearized equations as a first approximation. Some specific applications of END are presented, and the comparison with experiment and other theoretical approaches is discussed

Communication: Time-dependent optimized coupled-cluster method for multielectron dynamics

Science.gov (United States)

Sato, Takeshi; Pathak, Himadri; Orimo, Yuki; Ishikawa, Kenichi L.

2018-02-01

Time-dependent coupled-cluster method with time-varying orbital functions, called time-dependent optimized coupled-cluster (TD-OCC) method, is formulated for multielectron dynamics in an intense laser field. We have successfully derived the equations of motion for CC amplitudes and orthonormal orbital functions based on the real action functional, and implemented the method including double excitations (TD-OCCD) and double and triple excitations (TD-OCCDT) within the optimized active orbitals. The present method is size extensive and gauge invariant, a polynomial cost-scaling alternative to the time-dependent multiconfiguration self-consistent-field method. The first application of the TD-OCC method of intense-laser driven correlated electron dynamics in Ar atom is reported.

Similarity transformed coupled cluster response (ST-CCR) theory--a time-dependent similarity transformed equation-of-motion coupled cluster (STEOM-CC) approach.

Science.gov (United States)

Landau, Arie

2013-07-07

This paper presents a new method for calculating spectroscopic properties in the framework of response theory utilizing a sequence of similarity transformations (STs). The STs are preformed using the coupled cluster (CC) and Fock-space coupled cluster operators. The linear and quadratic response functions of the new similarity transformed CC response (ST-CCR) method are derived. The poles of the linear response yield excitation-energy (EE) expressions identical to the ones in the similarity transformed equation-of-motion coupled cluster (STEOM-CC) approach. ST-CCR and STEOM-CC complement each other, in analogy to the complementarity of CC response (CCR) and equation-of-motion coupled cluster (EOM-CC). ST-CCR/STEOM-CC and CCR/EOM-CC yield size-extensive and size-intensive EEs, respectively. Other electronic-properties, e.g., transition dipole strengths, are also size-extensive within ST-CCR, in contrast to STEOM-CC. Moreover, analysis suggests that in comparison with CCR, the ST-CCR expressions may be confined to a smaller subspace, however, the precise scope of the truncation can only be determined numerically. In addition, reformulation of the time-independent STEOM-CC using the same parameterization as in ST-CCR, as well as an efficient truncation scheme, is presented. The shown convergence of the time-dependent and time-independent expressions displays the completeness of the presented formalism.

Transforming parts of a differential equations system to difference equations as a method for run-time savings in NONMEM.

Science.gov (United States)

Petersson, K J F; Friberg, L E; Karlsson, M O

2010-10-01

Computer models of biological systems grow more complex as computing power increase. Often these models are defined as differential equations and no analytical solutions exist. Numerical integration is used to approximate the solution; this can be computationally intensive, time consuming and be a large proportion of the total computer runtime. The performance of different integration methods depend on the mathematical properties of the differential equations system at hand. In this paper we investigate the possibility of runtime gains by calculating parts of or the whole differential equations system at given time intervals, outside of the differential equations solver. This approach was tested on nine models defined as differential equations with the goal to reduce runtime while maintaining model fit, based on the objective function value. The software used was NONMEM. In four models the computational runtime was successfully reduced (by 59-96%). The differences in parameter estimates, compared to using only the differential equations solver were less than 12% for all fixed effects parameters. For the variance parameters, estimates were within 10% for the majority of the parameters. Population and individual predictions were similar and the differences in OFV were between 1 and -14 units. When computational runtime seriously affects the usefulness of a model we suggest evaluating this approach for repetitive elements of model building and evaluation such as covariate inclusions or bootstraps.

Wavelet and adaptive methods for time dependent problems and applications in aerosol dynamics

Science.gov (United States)

Guo, Qiang

Time dependent partial differential equations (PDEs) are widely used as mathematical models of environmental problems. Aerosols are now clearly identified as an important factor in many environmental aspects of climate and radiative forcing processes, as well as in the health effects of air quality. The mathematical models for the aerosol dynamics with respect to size distribution are nonlinear partial differential and integral equations, which describe processes of condensation, coagulation and deposition. Simulating the general aerosol dynamic equations on time, particle size and space exhibits serious difficulties because the size dimension ranges from a few nanometer to several micrometer while the spatial dimension is usually described with kilometers. Therefore, it is an important and challenging task to develop efficient techniques for solving time dependent dynamic equations. In this thesis, we develop and analyze efficient wavelet and adaptive methods for the time dependent dynamic equations on particle size and further apply them to the spatial aerosol dynamic systems. Wavelet Galerkin method is proposed to solve the aerosol dynamic equations on time and particle size due to the fact that aerosol distribution changes strongly along size direction and the wavelet technique can solve it very efficiently. Daubechies' wavelets are considered in the study due to the fact that they possess useful properties like orthogonality, compact support, exact representation of polynomials to a certain degree. Another problem encountered in the solution of the aerosol dynamic equations results from the hyperbolic form due to the condensation growth term. We propose a new characteristic-based fully adaptive multiresolution numerical scheme for solving the aerosol dynamic equation, which combines the attractive advantages of adaptive multiresolution technique and the characteristics method. On the aspect of theoretical analysis, the global existence and uniqueness of

Some aspects of the description of relativistic particles in external fields. [Time dependent and time independent potentials

Energy Technology Data Exchange (ETDEWEB)

Labonte, G

1973-01-01

We study the time description of the motion of relativistic particles in both the dependent and time independent potentials. The differential equations of motion considered are the standard linear spin zero and one half equations. They are always meaningful in the sense that, at all times, unique well defined operator valued distributions in the three space variables are determined. We discuss the problem of determining which set of creation and annihilation operators is relevant in a given problem. We examine the implementation of certain simple requirements which seem to be necessary in order for the mathematical formalism to be able to describe a physical system. We show that whenever the equation of motion is homogeneous, the study of all physical requirements reduces to studying Bogoliubov transformations between creation and annihilation operators. We study such transformations where we obtain some new important results concerning their general properties. We examine in detail a quantized field in presence of an external source, electrons and positrons acted upon by a plane electromagnetic wave, Dirac fields acted upon by potentials of the form A(x) delta (t) and A(x) THETA (t-t/sub 0/). We study Dirac fields in presence of potentials which have time dependences which can be represented by sequences of step functions. We then discuss the limiting case where the time dependence is continuous. We prove that the requirements that there exists a unitary evolution operator or that physical particles can be described are exactly equivalent. (auth)

Dynamics of Bose-Einstein condensates in a time-dependent trap

International Nuclear Information System (INIS)

Kumar, V. Ramesh; Radha, R.; Panigrahi, Prasanta K.

2008-01-01

In this paper, we generate the Lax pair for the one-dimensional Gross-Pitaevskii equation with time-dependent scattering length in the presence of a confining or expulsive harmonic time-dependent trap. We then exploit the Lax pair profitably to construct multisoliton solutions using gauge transformation from a trivial input solution. In particular, we have investigated the effect of both expulsive and confining traps on soliton interaction. Even though we find that the amplitude of the bright soliton relies upon the time-dependent scattering length and the external time-dependent trap with the velocity being dictated by the external trap alone, the observation of interdependence of the scattering length on the trap shows that the bright solitons not only can be compressed into a desirable width and amplitude but also can be remote controlled and driven anywhere in the plane by suitably maneuvering the external time-dependent trap alone

Approximate solution of space and time fractional higher order phase field equation

Science.gov (United States)

Shamseldeen, S.

2018-03-01

This paper is concerned with a class of space and time fractional partial differential equation (STFDE) with Riesz derivative in space and Caputo in time. The proposed STFDE is considered as a generalization of a sixth-order partial phase field equation. We describe the application of the optimal homotopy analysis method (OHAM) to obtain an approximate solution for the suggested fractional initial value problem. An averaged-squared residual error function is defined and used to determine the optimal convergence control parameter. Two numerical examples are studied, considering periodic and non-periodic initial conditions, to justify the efficiency and the accuracy of the adopted iterative approach. The dependence of the solution on the order of the fractional derivative in space and time and model parameters is investigated.

Construction of time-dependent dynamical invariants: A new approach

International Nuclear Information System (INIS)

Bertin, M. C.; Pimentel, B. M.; Ramirez, J. A.

2012-01-01

We propose a new way to obtain polynomial dynamical invariants of the classical and quantum time-dependent harmonic oscillator from the equations of motion. We also establish relations between linear and quadratic invariants, and discuss how the quadratic invariant can be related to the Ermakov invariant.

Time-dependent crashworthiness of polyurethane foam

Science.gov (United States)

Basit, Munshi Mahbubul; Cheon, Seong Sik

2018-05-01

Time-dependent stress-strain relationship as well as crashworthiness of polyurethane foam was investigated under constant impact energy with different velocities, considering inertia and strain-rate effects simultaneously during the impact testing. Even though the impact energies were same, the percentage in increase in densification strain due to higher impact velocities was found, which yielded the wider plateau region, i.e. growth in crashworthiness. This phenomenon is analyzed by the microstructure of polyurethane foam obtained from scanning electron microscopy. The equations, coupled with the Sherwood-Frost model and the impulse-momentum theory, were employed to build the constitutive equation of the polyurethane foam and calculate energy absorption capacity of the foam. The nominal stress-strain curves obtained from the constitutive equation were compared with results from impact tests and were found to be in good agreement. This study is dedicated to guiding designer use polyurethane foam in crashworthiness structures such as an automotive bumper system by providing crashworthiness data, determining the crush mode, and addressing a mathematical model of the crashworthiness.

Analog computing for a new nuclear reactor dynamic model based on a time-dependent second order form of the neutron transport equation

International Nuclear Information System (INIS)

Pirouzmand, Ahmad; Hadad, Kamal; Suh, Kune Y.

2011-01-01

This paper considers the concept of analog computing based on a cellular neural network (CNN) paradigm to simulate nuclear reactor dynamics using a time-dependent second order form of the neutron transport equation. Instead of solving nuclear reactor dynamic equations numerically, which is time-consuming and suffers from such weaknesses as vulnerability to transient phenomena, accumulation of round-off errors and floating-point overflows, use is made of a new method based on a cellular neural network. The state-of-the-art shows the CNN as being an alternative solution to the conventional numerical computation method. Indeed CNN is an analog computing paradigm that performs ultra-fast calculations and provides accurate results. In this study use is made of the CNN model to simulate the space-time response of scalar flux distribution in steady state and transient conditions. The CNN model also is used to simulate step perturbation in the core. The accuracy and capability of the CNN model are examined in 2D Cartesian geometry for two fixed source problems, a mini-BWR assembly, and a TWIGL Seed/Blanket problem. We also use the CNN model concurrently for a typical small PWR assembly to simulate the effect of temperature feedback, poisons, and control rods on the scalar flux distribution

Analytical Solution of Heat Conduction for Hollow Cylinders with Time-Dependent Boundary Condition and Time-Dependent Heat Transfer Coefficient

Directory of Open Access Journals (Sweden)

Te-Wen Tu

2015-01-01

Full Text Available An analytical solution for the heat transfer in hollow cylinders with time-dependent boundary condition and time-dependent heat transfer coefficient at different surfaces is developed for the first time. The methodology is an extension of the shifting function method. By dividing the Biot function into a constant plus a function and introducing two specially chosen shifting functions, the system is transformed into a partial differential equation with homogenous boundary conditions only. The transformed system is thus solved by series expansion theorem. Limiting cases of the solution are studied and numerical results are compared with those in the literature. The convergence rate of the present solution is fast and the analytical solution is simple and accurate. Also, the influence of physical parameters on the temperature distribution of a hollow cylinder along the radial direction is investigated.

A time-dependent dusty gas dynamic model of axisymmetric cometary jets

International Nuclear Information System (INIS)

Korosmezey, A.; Gombosi, T.I.

1990-01-01

The present time-dependent, axisymmetric dusty gas dynamical model of inner cometary atmospheres solves the coupled and time-dependent equations of continuity, momentum, and energy for a gas-dust mixture between the surface of the nucleus and 100 km, using an axisymmetric 40 x 40 grid structure. A novel numerical method employing a second-order accurate Godunov-type scheme with dimensional splitting is used to solve the time-dependent pde system. It is established that a subsolar dust spike not predicted by previous calculations is generated by narrow axisymmetric jets, together with a jet cone whose opening angle depends on the jet length. 28 refs

Time-dependent angularly averaged inverse transport

International Nuclear Information System (INIS)

Bal, Guillaume; Jollivet, Alexandre

2009-01-01

This paper concerns the reconstruction of the absorption and scattering parameters in a time-dependent linear transport equation from knowledge of angularly averaged measurements performed at the boundary of a domain of interest. Such measurement settings find applications in medical and geophysical imaging. We show that the absorption coefficient and the spatial component of the scattering coefficient are uniquely determined by such measurements. We obtain stability results on the reconstruction of the absorption and scattering parameters with respect to the measured albedo operator. The stability results are obtained by a precise decomposition of the measurements into components with different singular behavior in the time domain

Time-dependent quantum fluid density functional theory of hydrogen ...

Indian Academy of Sciences (India)

A time-dependent generalized non-linear Schrödinger equation (GNLSE) of motion was earlier derived in our laboratory by combining density functional theory and quantum fluid dynamics in threedimensional space. In continuation of the work reported previously, the GNLSE is applied to provide additional knowledge on ...

Time-dependent dilatancy for brittle rocks

Directory of Open Access Journals (Sweden)

Jie Li

2017-12-01

Full Text Available This paper presents a theoretical study on time-dependent dilatancy behaviors for brittle rocks. The theory employs a well-accepted postulation that macroscopically observed dilatancy originates from the expansion of microcracks. The mechanism and dynamic process that microcracks initiate from local stress concentration and grow due to localized tensile stress are analyzed. Then, by generalizing the results from the analysis of single cracks, a parameter and associated equations for its evolution are developed to describe the behaviors of the microcracks. In this circumstance, the relationship between microcracking and dilatancy can be established, and the theoretical equations for characterizing the process of rock dilatancy behaviors are derived. Triaxial compression and creep tests are conducted to validate the developed theory. With properly chosen model parameters, the theory yields a satisfactory accuracy in comparison with the experimental results.

Goedel, Penrose, anti-Mach: extra supersymmetries of time-dependent plane waves

International Nuclear Information System (INIS)

Blau, Matthias; O'Loughlin, Martin; Meessen, Patrick

2003-01-01

We prove that M-theory plane waves with extra supersymmetries are necessarily homogeneous (but possibly time-dependent), and we show by explicit construction that such time-dependent plane waves can admit extra supersymmetries. To that end we study the Penrose limits of Goedel-like metrics, show that the Penrose limit of the M-theory Goedel metric (with 20 supercharges) is generically a time-dependent homogeneous plane wave of the anti-Mach type, and display the four extra Killings spinors in that case. We conclude with some general remarks on the Killing spinor equations for homogeneous plane waves. (author)

Goedel, Penrose, anti-Mach: extra supersymmetries of time-dependent plane waves

Energy Technology Data Exchange (ETDEWEB)

Blau, Matthias; O' Loughlin, Martin; Meessen, Patrick [SISSA/ISAS, Via Beirut 2-4, 34014 Trieste (Italy)]. E-mail: meessen@sissa.it

2003-09-01

We prove that M-theory plane waves with extra supersymmetries are necessarily homogeneous (but possibly time-dependent), and we show by explicit construction that such time-dependent plane waves can admit extra supersymmetries. To that end we study the Penrose limits of Goedel-like metrics, show that the Penrose limit of the M-theory Goedel metric (with 20 supercharges) is generically a time-dependent homogeneous plane wave of the anti-Mach type, and display the four extra Killings spinors in that case. We conclude with some general remarks on the Killing spinor equations for homogeneous plane waves. (author)

Sensitivity analysis of time-dependent laminar flows

International Nuclear Information System (INIS)

Hristova, H.; Etienne, S.; Pelletier, D.; Borggaard, J.

2004-01-01

This paper presents a general sensitivity equation method (SEM) for time dependent incompressible laminar flows. The SEM accounts for complex parameter dependence and is suitable for a wide range of problems. The formulation is verified on a problem with a closed form solution obtained by the method of manufactured solution. Systematic grid convergence studies confirm the theoretical rates of convergence in both space and time. The methodology is then applied to pulsatile flow around a square cylinder. Computations show that the flow starts with symmetrical vortex shedding followed by a transition to the traditional Von Karman street (alternate vortex shedding). Simulations show that the transition phase manifests itself earlier in the sensitivity fields than in the flow field itself. Sensitivities are then demonstrated for fast evaluation of nearby flows and uncertainty analysis. (author)

Time-dependent shock acceleration of energetic electrons including synchrotron losses

International Nuclear Information System (INIS)

Fritz, K.; Webb, G.M.

1990-01-01

The present investigation of the time-dependent particle acceleration problem in strong shocks, including synchrotron radiation losses, solves the transport equation analytically by means of Laplace transforms. The particle distribution thus obtained is then transformed numerically into real space for the cases of continuous and impulsive injections of particles at the shock. While in the continuous case the steady-state spectrum undergoes evolution, impulsive injection is noted to yield such unpredicted features as a pile-up of high-energy particles or a steep power-law with time-dependent spectral index. The time-dependent calculations reveal varying spectral shapes and more complex features for the higher energies which may be useful in the interpretation of outburst spectra. 33 refs

On the algebraic approach to the time-dependent quadratic Hamiltonian

Energy Technology Data Exchange (ETDEWEB)

Urdaneta, Ines; Palma, Alejandro [Instituto de Fisica, Benemerita Universidad Autonoma de Puebla, Puebla (Mexico); Sandoval, Lourdes, E-mail: urdaneta@sirio.ifuap.buap.m [Facultad de Ciencias de la Computacion, Benemerita Universidad Autonoma de Puebla, Puebla (Mexico)

2010-09-24

The unitary operator V(t) that diagonalizes the time-dependent quadratic Hamiltonian (TDQH) into a time-dependent harmonic oscillator (TDHO) is obtained using a Lie algebra. The method involves a factorization of the TDQH into a TDHO through a unitary Bogoliubov transformation in terms of creation and annihilation operators with time-dependent coefficients. It is shown that this operator can be easily achieved by means of the factorization, together with the commonly known Wei-Norman theorem. We discuss the conditions under which this unitary operator converges to the evolution operator U(t) of the Schroedinger equation for the TDQH, giving then a straightforward calculation of the evolution operator with respect to the procedures published in the literature.

Lie symmetry analysis, conservation laws and exact solutions of the seventh-order time fractional Sawada–Kotera–Ito equation

Directory of Open Access Journals (Sweden)

Emrullah Yaşar

Full Text Available In this paper Lie symmetry analysis of the seventh-order time fractional Sawada–Kotera–Ito (FSKI equation with Riemann–Liouville derivative is performed. Using the Lie point symmetries of FSKI equation, it is shown that it can be transformed into a nonlinear ordinary differential equation of fractional order with a new dependent variable. In the reduced equation the derivative is in Erdelyi–Kober sense. Furthermore, adapting the Ibragimov’s nonlocal conservation method to time fractional partial differential equations, we obtain conservation laws of the underlying equation. In addition, we construct some exact travelling wave solutions for the FSKI equation using the sub-equation method. Keywords: Fractional Sawada–Kotera–Ito equation, Lie symmetry, Riemann–Liouville fractional derivative, Conservation laws, Exact solutions

Anomalous transport in fluid field with random waiting time depending on the preceding jump length

Science.gov (United States)

Zhang, Hong; Li, Guo-Hua

2016-11-01

Anomalous (or non-Fickian) transport behaviors of particles have been widely observed in complex porous media. To capture the energy-dependent characteristics of non-Fickian transport of a particle in flow fields, in the present paper a generalized continuous time random walk model whose waiting time probability distribution depends on the preceding jump length is introduced, and the corresponding master equation in Fourier-Laplace space for the distribution of particles is derived. As examples, two generalized advection-dispersion equations for Gaussian distribution and lévy flight with the probability density function of waiting time being quadratic dependent on the preceding jump length are obtained by applying the derived master equation. Project supported by the Foundation for Young Key Teachers of Chengdu University of Technology, China (Grant No. KYGG201414) and the Opening Foundation of Geomathematics Key Laboratory of Sichuan Province, China (Grant No. scsxdz2013009).

Time-Dependent Mean-Field Games with Logarithmic Nonlinearities

KAUST Repository

Gomes, Diogo A.

2015-10-06

In this paper, we prove the existence of classical solutions for time-dependent mean-field games with a logarithmic nonlinearity and subquadratic Hamiltonians. Because the logarithm is unbounded from below, this nonlinearity poses substantial mathematical challenges that have not been addressed in the literature. Our result is proven by recurring to a delicate argument which combines Lipschitz regularity for the Hamilton-Jacobi equation with estimates for the nonlinearity in suitable Lebesgue spaces. Lipschitz estimates follow from an application of the nonlinear adjoint method. These are then combined with a priori bounds for solutions of the Fokker-Planck equation and a concavity argument for the nonlinearity.

Time-Dependent Mean-Field Games with Logarithmic Nonlinearities

KAUST Repository

Gomes, Diogo A.; Pimentel, Edgard

2015-01-01

In this paper, we prove the existence of classical solutions for time-dependent mean-field games with a logarithmic nonlinearity and subquadratic Hamiltonians. Because the logarithm is unbounded from below, this nonlinearity poses substantial mathematical challenges that have not been addressed in the literature. Our result is proven by recurring to a delicate argument which combines Lipschitz regularity for the Hamilton-Jacobi equation with estimates for the nonlinearity in suitable Lebesgue spaces. Lipschitz estimates follow from an application of the nonlinear adjoint method. These are then combined with a priori bounds for solutions of the Fokker-Planck equation and a concavity argument for the nonlinearity.

Singular Hopf bifurcation in a differential equation with large state-dependent delay.

Science.gov (United States)

Kozyreff, G; Erneux, T

2014-02-08

We study the onset of sustained oscillations in a classical state-dependent delay (SDD) differential equation inspired by control theory. Owing to the large delays considered, the Hopf bifurcation is singular and the oscillations rapidly acquire a sawtooth profile past the instability threshold. Using asymptotic techniques, we explicitly capture the gradual change from nearly sinusoidal to sawtooth oscillations. The dependence of the delay on the solution can be either linear or nonlinear, with at least quadratic dependence. In the former case, an asymptotic connection is made with the Rayleigh oscillator. In the latter, van der Pol's equation is derived for the small-amplitude oscillations. SDD differential equations are currently the subject of intense research in order to establish or amend general theorems valid for constant-delay differential equation, but explicit analytical construction of solutions are rare. This paper illustrates the use of singular perturbation techniques and the unusual way in which solvability conditions can arise for SDD problems with large delays.

Analytical representation for solution of the neutron point kinetics equation with time-dependent reactivity and free of the stiffness character

International Nuclear Information System (INIS)

Silva, Milena Wollmann da

2013-01-01

In this work, we report a genuine analytical representation for the solution of the neutron point kinetics equation free of the stiffness character, assuming that the reactivity is a continuous and sectionally continuous function of time. To this end, we initially cast the point kinetics equation in a first order linear differential equation. Next, we split the corresponding matrix as a sum of a diagonal matrix with a matrix, whose components contain the off-diagonal elements. Next, expanding the neutron density and the delayed neutron precursors concentrations in a truncated series, and replacing these expansions in the matrix equation, we come out with an equation, which allows to construct a recursive system, a first order matrix differential equation with source. The fundamental characteristic of this system relies on the fact that the corresponding matrix is diagonal, meanwhile the source term is written in terms of the matrix with the off-diagonal components. Further, the first equation of the recursive system has no source and satisfies the initial conditions. On the other hand, the remaining equations satisfy the null initial condition. Due to the diagonal feature of the matrix, we attain analytical solutions for these recursive equations. We also mention that we evaluate the results for any time value, without the analytical continuity because the purposed solution is free on the stiffness character. Finally, we present numerical simulations and comparisons against literature results, considering specific the applications for the following reactivity functions: constant, step, ramp, and sine. (author)

THE WIGNER–FOKKER–PLANCK EQUATION: STATIONARY STATES AND LARGE TIME BEHAVIOR

KAUST Repository

ARNOLD, ANTON

2012-11-01

We consider the linear WignerFokkerPlanck equation subject to confining potentials which are smooth perturbations of the harmonic oscillator potential. For a certain class of perturbations we prove that the equation admits a unique stationary solution in a weighted Sobolev space. A key ingredient of the proof is a new result on the existence of spectral gaps for FokkerPlanck type operators in certain weighted L 2-spaces. In addition we show that the steady state corresponds to a positive density matrix operator with unit trace and that the solutions of the time-dependent problem converge towards the steady state with an exponential rate. © 2012 World Scientific Publishing Company.

Numerical Solution of Time-Dependent Problems with a Fractional-Power Elliptic Operator

Science.gov (United States)

Vabishchevich, P. N.

2018-03-01

A time-dependent problem in a bounded domain for a fractional diffusion equation is considered. The first-order evolution equation involves a fractional-power second-order elliptic operator with Robin boundary conditions. A finite-element spatial approximation with an additive approximation of the operator of the problem is used. The time approximation is based on a vector scheme. The transition to a new time level is ensured by solving a sequence of standard elliptic boundary value problems. Numerical results obtained for a two-dimensional model problem are presented.

Implicit time-dependent finite different algorithm for quench simulation

International Nuclear Information System (INIS)

Koizumi, Norikiyo; Takahashi, Yoshikazu; Tsuji, Hiroshi

1994-12-01

A magnet in a fusion machine has many difficulties in its application because of requirement of a large operating current, high operating field and high breakdown voltage. A cable-in-conduit (CIC) conductor is the best candidate to overcome these difficulties. However, there remained uncertainty in a quench event in the cable-in-conduit conductor because of a difficulty to analyze a fluid dynamics equation. Several scientists, then, developed the numerical code for the quench simulation. However, most of them were based on an explicit time-dependent finite difference scheme. In this scheme, a discrete time increment is strictly restricted by CFL (Courant-Friedrichs-Lewy) condition. Therefore, long CPU time was consumed for the quench simulation. Authors, then, developed a new quench simulation code, POCHI1, which is based on an implicit time dependent scheme. In POCHI1, the fluid dynamics equation is linearlized according to a procedure applied by Beam and Warming and then, a tridiagonal system can be offered. Therefore, no iteration is necessary to solve the fluid dynamics equation. This leads great reduction of the CPU time. Also, POCHI1 can cope with non-linear boundary condition. In this study, comparison with experimental results was carried out. The normal zone propagation behavior was investigated in two samples of CIC conductors which had different hydraulic diameters. The measured and simulated normal zone propagation length showed relatively good agreement. However, the behavior of the normal voltage shows a little disagreement. These results indicate necessity to improve the treatment of the heat transfer coefficient in the turbulent flow region and the electric resistivity of the copper stabilizer in high temperature and high field region. (author)

A numerical study of time-dependent Schrödinger equation for ...

Indian Academy of Sciences (India)

Unknown

Theoretical Chemistry Group, Department of Chemistry, Panjab University,. Chandigarh 160 ... probability, potential energy curve and dipole moment. ... quantum Monte Carlo (DQMC)-type equation.23 The system is then evolved in imaginary.

Coherent states for the time dependent harmonic oscillator: the step function

International Nuclear Information System (INIS)

Moya-Cessa, Hector; Fernandez Guasti, Manuel

2003-01-01

We study the time evolution for the quantum harmonic oscillator subjected to a sudden change of frequency. It is based on an approximate analytic solution to the time dependent Ermakov equation for a step function. This approach allows for a continuous treatment that differs from former studies that involve the matching of two time independent solutions at the time when the step occurs

A numerical study of adaptive space and time discretisations for Gross–Pitaevskii equations

Science.gov (United States)

Thalhammer, Mechthild; Abhau, Jochen

2012-01-01

As a basic principle, benefits of adaptive discretisations are an improved balance between required accuracy and efficiency as well as an enhancement of the reliability of numerical computations. In this work, the capacity of locally adaptive space and time discretisations for the numerical solution of low-dimensional nonlinear Schrödinger equations is investigated. The considered model equation is related to the time-dependent Gross–Pitaevskii equation arising in the description of Bose–Einstein condensates in dilute gases. The performance of the Fourier-pseudo spectral method constrained to uniform meshes versus the locally adaptive finite element method and of higher-order exponential operator splitting methods with variable time stepsizes is studied. Numerical experiments confirm that a local time stepsize control based on a posteriori local error estimators or embedded splitting pairs, respectively, is effective in different situations with an enhancement either in efficiency or reliability. As expected, adaptive time-splitting schemes combined with fast Fourier transform techniques are favourable regarding accuracy and efficiency when applied to Gross–Pitaevskii equations with a defocusing nonlinearity and a mildly varying regular solution. However, the numerical solution of nonlinear Schrödinger equations in the semi-classical regime becomes a demanding task. Due to the highly oscillatory and nonlinear nature of the problem, the spatial mesh size and the time increments need to be of the size of the decisive parameter 0Fourier pseudo-spectral and the finite element method. Nevertheless, for smaller parameter values locally adaptive time discretisations facilitate to determine the time stepsizes sufficiently small in order that the numerical approximation captures correctly the behaviour of the analytical solution. Further illustrations for Gross–Pitaevskii equations with a focusing nonlinearity or a sharp Gaussian as initial condition, respectively

Application of the nodal method RTN-0 for the solution of the neutron diffusion equation dependent of time in hexagonal-Z geometry

International Nuclear Information System (INIS)

Esquivel E, J.; Alonso V, G.; Del Valle G, E.

2015-09-01

The solution of the neutron diffusion equation either for reactors in steady state or time dependent, is obtained through approximations generated by implementing of nodal methods such as RTN-0 (Raviart-Thomas-Nedelec of zero index), which is used in this study. Since the nodal methods are applied in quadrangular geometries, in this paper a technique in which the hexagonal geometry through the transfinite interpolation of Gordon-Hall becomes the appropriate geometry to make use of the nodal method RTN-0 is presented. As a result, a computer program was developed, whereby is possible to obtain among other results the neutron multiplication effective factor (k eff ), and the distribution of radial and/or axial power. To verify the operation of the code, was applied to three benchmark problems: in the first two reactors VVER and FBR, results k eff and power distribution are obtained, considering the steady state case of reactor; while the third problem a type VVER is analyzed, in its case dependent of time, which qualitative results are presented on the behavior of the reactor power. (Author)

Time-dependent quantum many-body systems. Linear response, electronic transport, and reduced density matrices

International Nuclear Information System (INIS)

Appel, H.

2007-05-01

In part I of this work we present a double-pole approximation (DPA) to the response equations of time-dependent density functional theory (TDDFT). The double-pole approximation provides an exact description of systems with two strongly coupled excitations which are isolated from the rest of the spectrum. In contrast to the traditional single-pole approximation of TDDFT the DPA also yields corrections to the Kohn-Sham oscillator strengths. We also demonstrate how to invert the double-pole solution which allows us to predict matrix elements of the exchange-correlation kernel f xc from experimental input. We attempt some first steps towards a time-dependent generalization of reduced density matrix functional theory (RDMFT). In part II we derive equations of motion for natural orbitals and occupation numbers. Using the equation of motion for the occupation numbers we show that an adiabatic extension of presently known ground-state functionals of static RDMFT always leads to occupation numbers which are constant in time. From the stationary conditions of the equations of motion for the N-body correlations (correlated parts of the N-body matrices) we derive a new class of ground-state functionals which can be used in static RDMFT. Applications are presented for a one-dimensional model system where the time-dependent many-body Schroedinger equation can be propagated numerically. We use optimal control theory to find optimized laser pulses for transitions in a model for atomic Helium. From the numerically exact correlated wavefunction we extract the exact time evolution of natural orbitals and occupation numbers for (i) laser-driven Helium and (ii) electron-ion scattering. Part III of this work considers time-dependent quantum transport within TDDFT. We present an algorithm for the calculation of extended eigenstates of single-particle Hamiltonians which is especially tailored to a finite-difference discretization of the Schroedinger equation. We consider the propagation

Time-dependent quantum many-body systems. Linear response, electronic transport, and reduced density matrices

Energy Technology Data Exchange (ETDEWEB)

Appel, H.

2007-05-15

In part I of this work we present a double-pole approximation (DPA) to the response equations of time-dependent density functional theory (TDDFT). The double-pole approximation provides an exact description of systems with two strongly coupled excitations which are isolated from the rest of the spectrum. In contrast to the traditional single-pole approximation of TDDFT the DPA also yields corrections to the Kohn-Sham oscillator strengths. We also demonstrate how to invert the double-pole solution which allows us to predict matrix elements of the exchange-correlation kernel f{sub xc} from experimental input. We attempt some first steps towards a time-dependent generalization of reduced density matrix functional theory (RDMFT). In part II we derive equations of motion for natural orbitals and occupation numbers. Using the equation of motion for the occupation numbers we show that an adiabatic extension of presently known ground-state functionals of static RDMFT always leads to occupation numbers which are constant in time. From the stationary conditions of the equations of motion for the N-body correlations (correlated parts of the N-body matrices) we derive a new class of ground-state functionals which can be used in static RDMFT. Applications are presented for a one-dimensional model system where the time-dependent many-body Schroedinger equation can be propagated numerically. We use optimal control theory to find optimized laser pulses for transitions in a model for atomic Helium. From the numerically exact correlated wavefunction we extract the exact time evolution of natural orbitals and occupation numbers for (i) laser-driven Helium and (ii) electron-ion scattering. Part III of this work considers time-dependent quantum transport within TDDFT. We present an algorithm for the calculation of extended eigenstates of single-particle Hamiltonians which is especially tailored to a finite-difference discretization of the Schroedinger equation. We consider the

A Sequential, Implicit, Wavelet-Based Solver for Multi-Scale Time-Dependent Partial Differential Equations

Directory of Open Access Journals (Sweden)

Donald A. McLaren

2013-04-01

Full Text Available This paper describes and tests a wavelet-based implicit numerical method for solving partial differential equations. Intended for problems with localized small-scale interactions, the method exploits the form of the wavelet decomposition to divide the implicit system created by the time-discretization into multiple smaller systems that can be solved sequentially. Included is a test on a basic non-linear problem, with both the results of the test, and the time required to calculate them, compared with control results based on a single system with fine resolution. The method is then tested on a non-trivial problem, its computational time and accuracy checked against control results. In both tests, it was found that the method requires less computational expense than the control. Furthermore, the method showed convergence towards the fine resolution control results.

Tracer kinetics: Modelling by partial differential equations of inhomogeneous compartments with age-dependent elimination rates. Pt. 1

International Nuclear Information System (INIS)

Winkler, E.

1991-01-01

Mathematical models in tracer kinetics are usually based on ordinary differential equations which correspond to a system of kinetically homogeneous compartments (standard compartments). A generalization is possible by the admission of inhomogeneities in the behaviour of the elements belonging to a compartment. The important special case of the age-dependence of elimination rates is treated in its deterministic version. It leads to partial different equations (i.e., systems with distributed coefficients) with the 'age' or the 'residence time' of an element of the compartment as a variable additional to 'time'. The basic equations for one generalized compartment and for systems of such compartments are given together with their general solutions. (orig.) [de

Long-time behaviour of discretizations of the simple pendulum equation

Energy Technology Data Exchange (ETDEWEB)

Cieslinski, Jan L [Uniwersytet w Bialymstoku, Wydzial Fizyki, ul. Lipowa 41, 15-424 Bialystok (Poland); Ratkiewicz, Boguslaw [Doctoral Studies, Wydzial Fizyki, Uniwersytet Adama Mickiewicza, Poznan (Poland)], E-mail: janek@alpha.uwb.edu.pl, E-mail: bograt@poczta.onet.pl

2009-03-13

We compare several discretizations of the simple pendulum equation in a series of numerical experiments. The stress is put on the long-time behaviour. The chosen numerical schemes are either symplectic maps or integrable (energy-preserving) maps, or both. Therefore, they preserve qualitative features of solutions (such as periodicity). We describe characteristic periodic time dependences of numerical estimates of the period and the amplitude, and explain them as systematic numerical by-effects produced by any method. Finally, we propose a new numerical scheme which is a modification of the discrete gradient method. This modified discrete gradient method preserves (almost exactly) the period of small oscillations for any time step.

Long-time behaviour of discretizations of the simple pendulum equation

International Nuclear Information System (INIS)

Cieslinski, Jan L; Ratkiewicz, Boguslaw

2009-01-01

We compare several discretizations of the simple pendulum equation in a series of numerical experiments. The stress is put on the long-time behaviour. The chosen numerical schemes are either symplectic maps or integrable (energy-preserving) maps, or both. Therefore, they preserve qualitative features of solutions (such as periodicity). We describe characteristic periodic time dependences of numerical estimates of the period and the amplitude, and explain them as systematic numerical by-effects produced by any method. Finally, we propose a new numerical scheme which is a modification of the discrete gradient method. This modified discrete gradient method preserves (almost exactly) the period of small oscillations for any time step

MINARET: Towards a time-dependent neutron transport parallel solver

International Nuclear Information System (INIS)

Baudron, A.M.; Lautard, J.J.; Maday, Y.; Mula, O.

2013-01-01

We present the newly developed time-dependent 3D multigroup discrete ordinates neutron transport solver that has recently been implemented in the MINARET code. The solver is the support for a study about computing acceleration techniques that involve parallel architectures. In this work, we will focus on the parallelization of two of the variables involved in our equation: the angular directions and the time. This last variable has been parallelized by a (time) domain decomposition method called the para-real in time algorithm. (authors)

Time-dependent shell-model theory of dissipative heavy-ion collisions

International Nuclear Information System (INIS)

Ayik, S.; Noerenberg, W.

1982-01-01

A transport theory is formulated within a time-dependent shell-model approach. Time averaging of the equations for macroscopic quantities lead to irreversibility and justifies weak-coupling limit and Markov approximation for the (energy-conserving) one- and two-body collision terms. Two coupled equations for the occupation probabilities of dynamical single-particle states and for the collective variable are derived and explicit formulas for transition rates, dynamical forces, mass parameters and friction coefficients are given. The applicability of the formulation in terms of characteristic quantities of nuclear systems is considered in detail and some peculiarities due to memory effects in the initial equilibration process of heavy-ion collisions are discussed. (orig.)

ON PARTIAL DIFFERENTIAL AND DIFFERENCE EQUATIONS WITH SYMMETRIES DEPENDING ON ARBITRARY FUNCTIONS

Directory of Open Access Journals (Sweden)

Giorgio Gubbiotti

2016-06-01

Full Text Available In this note we present some ideas on when Lie symmetries, both point and generalized, can depend on arbitrary functions. We show a few examples, both in partial differential and partial difference equations where this happens. Moreover we show that the infinitesimal generators of generalized symmetries depending on arbitrary functions, both for continuous and discrete equations, effectively play the role of master symmetries.

Time-dependent approach to collisional ionization using exterior complex scaling

International Nuclear Information System (INIS)

McCurdy, C. William; Horner, Daniel A.; Rescigno, Thomas N.

2002-01-01

We present a time-dependent formulation of the exterior complex scaling method that has previously been used to treat electron-impact ionization of the hydrogen atom accurately at low energies. The time-dependent approach solves a driven Schroedinger equation, and scales more favorably with the number of electrons than the original formulation. The method is demonstrated in calculations for breakup processes in two dimensions (2D) and three dimensions for systems involving short-range potentials and in 2D for electron-impact ionization in the Temkin-Poet model for electron-hydrogen atom collisions

Algebraic time-dependent variational approach to dynamical calculations

International Nuclear Information System (INIS)

Shi, S.; Rabitz, H.

1988-01-01

A set of time-dependent basis states is obtained with a group of unitary transformations generated by a Lie algebra. Applying the time-dependent variational principle to the trial function subspace constructed from the linear combination of the time-dependent basis states gives rise to a set of ''classical'' equations of motion for the group parameters and the expansion coefficients from which the time evolution of the system state can be determined. The formulation is developed for a general Lie algebra as well as for the commonly encountered algebra containing homogeneous polynominal products of the coordinate Q and momentum P operators (or equivalently the boson creation a/sup dagger/ and annihilation a operators) of order 0, 1, and 2. Explicit expressions for the transition amplitudes are derived by virtue of the cannonical transformation properties of the unitary transformation. The applicability of the present formalism in a variety of problems is implied by two illustrative examples: (a) a parametric amplifier; (b) the collinear collision of an atom with a Morse oscillator

Relativistic time-dependent Fermion-mass renormalization using statistical regularization

Science.gov (United States)

Kutnink, Timothy; McMurray, Christian; Santrach, Amelia; Hockett, Sarah; Barcus, Scott; Petridis, Athanasios

2017-09-01

The time-dependent electromagnetically self-coupled Dirac equation is solved numerically by means of the staggered-leap-frog algorithm with reflecting boundary conditions. The stability region of the method versus the interaction strength and the spatial-grid size over time-step ratio is established. The expectation values of several dynamic operators are then evaluated as functions of time. These include the fermion and electromagnetic energies and the fermion dynamic mass. There is a characteristic, non-exponential, oscillatory dependence leading to asymptotic constants of these expectation values. In the case of the fermion mass this amounts to renormalization. The dependence of the expectation values on the spatial-grid size is evaluated in detail. Furthermore, the contribution of positive and negative energy states to the asymptotic values and the gauge fields is analyzed. Statistical regularization, employing a canonical ensemble whose temperature is the inverse of the grid size, is used to remove the grid-size and momentum-dependence and produce a finite result in the continuum limit.

Exact results on diffusion in a piecewise linear potential with a time-dependent sink

Energy Technology Data Exchange (ETDEWEB)

Diwaker, E-mail: diwakerphysics@gmail.com [Central University of Himachal Pradesh, School of Physical and Astronomical Sciences (India); Chakraborty, Aniruddha [Indian Institute of Technology Mandi (India)

2016-02-15

The Smoluchowski equation with a time-dependent sink term is solved exactly. In this method, knowing the probability distribution P(0, s) at the origin, allows deriving the probability distribution P(x, s) at all positions. Exact solutions of the Smoluchowski equation are also provided in different cases where the sink term has linear, constant, inverse, and exponential variation in time.

Darboux Transformations for Energy-Dependent Potentials and the Klein–Gordon Equation

International Nuclear Information System (INIS)

Schulze-Halberg, Axel

2013-01-01

We construct explicit Darboux transformations for a generalized Schrödinger-type equation with energy-dependent potential, a special case of which is the stationary Klein–Gordon equation. Our results complement and generalize former findings (Lin et al., Phys Lett A 362:212–214, 2007).

Reparametrization invariance and the Schroedinger equation

International Nuclear Information System (INIS)

Tkach, V.I.; Pashnev, A.I.; Rosales, J.J.

1999-01-01

A time-dependent Schroedinger equation for systems invariant under the reparametrization of time is considered. We develop the two-stage procedure of construction such systems from a given initial ones, which are not invariant under the time reparametrization. One of the first-class constraints of the systems in such description becomes the time-dependent Schroedinger equation. The procedure is applicable in the supersymmetric theories as well. The n = 2 supersymmetric quantum mechanics is coupled to world-line supergravity, and the local supersymmetric action is constructed leading to the square root representation of the time-dependent Schroedinger equation

Solving the Sea-Level Equation in an Explicit Time Differencing Scheme

Science.gov (United States)

Klemann, V.; Hagedoorn, J. M.; Thomas, M.

2016-12-01

In preparation of coupling the solid-earth to an ice-sheet compartment in an earth-system model, the dependency of initial topography on the ice-sheet history and viscosity structure has to be analysed. In this study, we discuss this dependency and how it influences the reconstruction of former sea level during a glacial cycle. The modelling is based on the VILMA code in which the field equations are solved in the time domain applying an explicit time-differencing scheme. The sea-level equation is solved simultaneously in the same explicit scheme as the viscoleastic field equations (Hagedoorn et al., 2007). With the assumption of only small changes, we neglect the iterative solution at each time step as suggested by e.g. Kendall et al. (2005). Nevertheless, the prediction of the initial paleo topography in case of moving coastlines remains to be iterated by repeated integration of the whole load history. The sensitivity study sketched at the beginning is accordingly motivated by the question if the iteration of the paleo topography can be replaced by a predefined one. This study is part of the German paleoclimate modelling initiative PalMod. Lit:Hagedoorn JM, Wolf D, Martinec Z, 2007. An estimate of global mean sea-level rise inferred from tide-gauge measurements using glacial-isostatic models consistent with the relative sea-level record. Pure appl. Geophys. 164: 791-818, doi:10.1007/s00024-007-0186-7Kendall RA, Mitrovica JX, Milne GA, 2005. On post-glacial sea level - II. Numerical formulation and comparative reesults on spherically symmetric models. Geophys. J. Int., 161: 679-706, doi:10.1111/j.365-246.X.2005.02553.x

The Laplace transformation of adjoint transport equations

International Nuclear Information System (INIS)

Hoogenboom, J.E.

1985-01-01

A clarification is given of the difference between the equation adjoint to the Laplace-transformed time-dependent transport equation and the Laplace-transformed time-dependent adjoint transport equation. Proper procedures are derived to obtain the Laplace transform of the instantaneous detector response. (author)

Towards a frequency-dependent discrete maximum principle for the implicit Monte Carlo equations

International Nuclear Information System (INIS)

Wollaber, Allan B.; Larsen, Edward W.; Densmore, Jeffery D.

2011-01-01

It has long been known that temperature solutions of the Implicit Monte Carlo (IMC) equations can exceed the external boundary temperatures, a so-called violation of the 'maximum principle'. Previous attempts at prescribing a maximum value of the time-step size Δ t that is sufficient to eliminate these violations have recommended a Δ t that is typically too small to be used in practice and that appeared to be much too conservative when compared to numerical solutions of the IMC equations for practical problems. In this paper, we derive a new estimator for the maximum time-step size that includes the spatial-grid size Δ x . This explicitly demonstrates that the effect of coarsening Δ x is to reduce the limitation on Δ t , which helps explain the overly conservative nature of the earlier, grid-independent results. We demonstrate that our new time-step restriction is a much more accurate means of predicting violations of the maximum principle. We discuss how the implications of the new, grid-dependent time-step restriction can impact IMC solution algorithms. (author)

Elliptic and solitary wave solutions for Bogoyavlenskii equations system, couple Boiti-Leon-Pempinelli equations system and Time-fractional Cahn-Allen equation

Directory of Open Access Journals (Sweden)

Mostafa M.A. Khater

Full Text Available In this article and for the first time, we introduce and describe Khater method which is a new technique for solving nonlinear partial differential equations (PDEs.. We apply this method for each of the following models Bogoyavlenskii equation, couple Boiti-Leon-Pempinelli system and Time-fractional Cahn-Allen equation. Khater method is very powerful, Effective, felicitous and fabulous method to get exact and solitary wave solution of (PDEs.. Not only just like that but it considers too one of the general methods for solving that kind of equations since it involves some methods as we will see in our discuss of the results. We make a comparison between the results of this new method and another method. Keywords: Bogoyavlenskii equations system, Couple Boiti-Leon-Pempinelli equations system, Time-fractional Cahn-Allen equation, Khater method, Traveling wave solutions, Solitary wave solutions

On time-dependent diffusion coefficients arising from stochastic processes with memory

Science.gov (United States)

Carpio-Bernido, M. Victoria; Barredo, Wilson I.; Bernido, Christopher C.

2017-08-01

Time-dependent diffusion coefficients arise from anomalous diffusion encountered in many physical systems such as protein transport in cells. We compare these coefficients with those arising from analysis of stochastic processes with memory that go beyond fractional Brownian motion. Facilitated by the Hida white noise functional integral approach, diffusion propagators or probability density functions (pdf) are obtained and shown to be solutions of modified diffusion equations with time-dependent diffusion coefficients. This should be useful in the study of complex transport processes.

Time dependent non-extinction probability for prompt critical systems

International Nuclear Information System (INIS)

Gregson, M. W.; Prinja, A. K.

2009-01-01

The time dependent non-extinction probability equation is presented for slab geometry. Numerical solutions are provided for a nested inner/outer iteration routine where the fission terms (both linear and non-linear) are updated and then held fixed over the inner scattering iteration. Time dependent results are presented highlighting the importance of the injection position and angle. The iteration behavior is also described as the steady state probability of initiation is approached for both small and large time steps. Theoretical analysis of the nested iteration scheme is shown and highlights poor numerical convergence for marginally prompt critical systems. An acceleration scheme for the outer iterations is presented to improve convergence of such systems. Theoretical analysis of the acceleration scheme is also provided and the associated decrease in computational run time addressed. (authors)

Equal-Time and Equal-Space Poisson Brackets of the N -Component Coupled NLS Equation

International Nuclear Information System (INIS)

Zhou Ru-Guang; Li Pei-Yao; Gao Yuan

2017-01-01

Two Poisson brackets for the N-component coupled nonlinear Schrödinger (NLS) equation are derived by using the variantional principle. The first one is called the equal-time Poisson bracket which does not depend on time but only on the space variable. Actually it is just the usual one describing the time evolution of system in the traditional theory of integrable Hamiltonian systems. The second one is equal-space and new. It is shown that the spatial part of Lax pair with respect to the equal-time Poisson bracket and temporal part of Lax pair with respect to the equal-space Poisson bracket share the same r-matrix formulation. These properties are similar to that of the NLS equation. (paper)

Derivation of the RPA (Random Phase Approximation) Equation of ATDDFT (Adiabatic Time Dependent Density Functional Ground State Response Theory) from an Excited State Variational Approach Based on the Ground State Functional.

Science.gov (United States)

Ziegler, Tom; Krykunov, Mykhaylo; Autschbach, Jochen

2014-09-09

The random phase approximation (RPA) equation of adiabatic time dependent density functional ground state response theory (ATDDFT) has been used extensively in studies of excited states. It extracts information about excited states from frequency dependent ground state response properties and avoids, thus, in an elegant way, direct Kohn-Sham calculations on excited states in accordance with the status of DFT as a ground state theory. Thus, excitation energies can be found as resonance poles of frequency dependent ground state polarizability from the eigenvalues of the RPA equation. ATDDFT is approximate in that it makes use of a frequency independent energy kernel derived from the ground state functional. It is shown in this study that one can derive the RPA equation of ATDDFT from a purely variational approach in which stationary states above the ground state are located using our constricted variational DFT (CV-DFT) method and the ground state functional. Thus, locating stationary states above the ground state due to one-electron excitations with a ground state functional is completely equivalent to solving the RPA equation of TDDFT employing the same functional. The present study is an extension of a previous work in which we demonstrated the equivalence between ATDDFT and CV-DFT within the Tamm-Dancoff approximation.

Time-domain representation of frequency-dependent foundation impedance functions

Science.gov (United States)

Safak, E.

2006-01-01

Foundation impedance functions provide a simple means to account for soil-structure interaction (SSI) when studying seismic response of structures. Impedance functions represent the dynamic stiffness of the soil media surrounding the foundation. The fact that impedance functions are frequency dependent makes it difficult to incorporate SSI in standard time-history analysis software. This paper introduces a simple method to convert frequency-dependent impedance functions into time-domain filters. The method is based on the least-squares approximation of impedance functions by ratios of two complex polynomials. Such ratios are equivalent, in the time-domain, to discrete-time recursive filters, which are simple finite-difference equations giving the relationship between foundation forces and displacements. These filters can easily be incorporated into standard time-history analysis programs. Three examples are presented to show the applications of the method.

Two-dimensional time dependent Riemann solvers for neutron transport

International Nuclear Information System (INIS)

Brunner, Thomas A.; Holloway, James Paul

2005-01-01

A two-dimensional Riemann solver is developed for the spherical harmonics approximation to the time dependent neutron transport equation. The eigenstructure of the resulting equations is explored, giving insight into both the spherical harmonics approximation and the Riemann solver. The classic Roe-type Riemann solver used here was developed for one-dimensional problems, but can be used in multidimensional problems by treating each face of a two-dimensional computation cell in a locally one-dimensional way. Several test problems are used to explore the capabilities of both the Riemann solver and the spherical harmonics approximation. The numerical solution for a simple line source problem is compared to the analytic solution to both the P 1 equation and the full transport solution. A lattice problem is used to test the method on a more challenging problem

Time-independent limit of a creep-recovery constitutive equation

International Nuclear Information System (INIS)

Chang, S.J.

1984-01-01

The effect of strain recovery is taken into consideration in ORNL efforts to establish unified constitutive equations for time-dependent plastic deformation for metals at elevated temperatures. Representation by internal state variables and Rice's flow potential are under consideration. Here the growth law for the internal state variables is discussed and interpreted in terms of a generalized form of the kinematic hardening condition of Prager. The yield condition is obtained from the flow potential representation of the inelastic strain rate. A consistency condition is derived from the yield condition and leads to a flow rule which assumes a slightly general form as compared with that of the classical plasticity due to the effect of strain recovery and the time-dependent property of the yield condition. Based on this representation, the time-independent limit is discussed. From a vanishing effect of recovery and a rate-independent limit for the yield condition at low temperature, this flow rule reduces to the well-known form of time-independent plasticity with a kinematic hardening condition. The duration of time (the characteristic time) required for the inelastic strain to reach its saturated value is defined for the inelastic loading condition. It provides the measure of a minimum duration of time which is required for a valid approximation made by the time-independent plasticity model

Microscopic description of fission dynamics: Toward a 3D computation of the time dependent GCM equation

Directory of Open Access Journals (Sweden)

Regnier D.

2017-01-01

Full Text Available Accurate knowledge of fission fragment yields is an essential ingredient of numerous applications ranging from the formation of elements in the r-process to fuel cycle optimization in nuclear energy. The need for a predictive theory applicable where no data is available, together with the variety of potential applications, is an incentive to develop a fully microscopic approach to fission dynamics. One of the most promising theoretical frameworks is the time dependent generator coordinate method (TDGCM applied under the Gaussian overlap approximation (GOA. However, the computational cost of this method makes it difficult to perform calculations with more than two collective degree of freedom. Meanwhile, it is well-known from both semi-phenomenological and fully microscopic approaches that at least four or five dimensions may play a role in the dynamics of fission. To overcome this limitation, we develop the code FELIX aiming to solve the TDGCM+GOA equation for an arbitrary number of collective variables. In this talk, we report the recent progress toward this enriched description of fission dynamics. We will briefly present the numerical methods adopted as well as the status of the latest version of FELIX. Finally, we will discuss fragments yields obtained within this approach for the low energy fission of major actinides.

Microscopic description of fission dynamics: Toward a 3D computation of the time dependent GCM equation

Science.gov (United States)

Regnier, D.; Dubray, N.; Schunck, N.; Verrière, M.

2017-09-01

Accurate knowledge of fission fragment yields is an essential ingredient of numerous applications ranging from the formation of elements in the r-process to fuel cycle optimization in nuclear energy. The need for a predictive theory applicable where no data is available, together with the variety of potential applications, is an incentive to develop a fully microscopic approach to fission dynamics. One of the most promising theoretical frameworks is the time dependent generator coordinate method (TDGCM) applied under the Gaussian overlap approximation (GOA). However, the computational cost of this method makes it difficult to perform calculations with more than two collective degree of freedom. Meanwhile, it is well-known from both semi-phenomenological and fully microscopic approaches that at least four or five dimensions may play a role in the dynamics of fission. To overcome this limitation, we develop the code FELIX aiming to solve the TDGCM+GOA equation for an arbitrary number of collective variables. In this talk, we report the recent progress toward this enriched description of fission dynamics. We will briefly present the numerical methods adopted as well as the status of the latest version of FELIX. Finally, we will discuss fragments yields obtained within this approach for the low energy fission of major actinides.

The nonlinear heat equation with state–dependent parameters and its connection to the Burgers’ and the potential Burgers’ equation

DEFF Research Database (Denmark)

Backi, Christoph Josef; Bendtsen, Jan Dimon; Leth, John-Josef

2014-01-01

In this work the stability properties of a nonlinear partial differential equation (PDE) with state–dependent parameters is investigated. Among other things, the PDE describes freezing of foodstuff, and is closely related to the (Potential) Burgers’ Equation. We show that for certain forms of coe...

Time-dependent resonant tunnelling for parallel-coupled double quantum dots

International Nuclear Information System (INIS)

Dong Bing; Djuric, Ivana; Cui, H L; Lei, X L

2004-01-01

We derive the quantum rate equations for an Aharonov-Bohm interferometer with two vertically coupled quantum dots embedded in each of two arms by means of the nonequilibrium Green function in the sequential tunnelling regime. Based on these equations, we investigate time-dependent resonant tunnelling under a small amplitude irradiation and find that the resonant photon-assisted tunnelling peaks in photocurrent demonstrate a combination behaviour of Fano and Lorentzian resonances due to the interference effect between the two pathways in this parallel configuration, which is controllable by threading the magnetic flux inside this device

Quadratic time dependent Hamiltonians and separation of variables

International Nuclear Information System (INIS)

Anzaldo-Meneses, A.

2017-01-01

Time dependent quantum problems defined by quadratic Hamiltonians are solved using canonical transformations. The Green’s function is obtained and a comparison with the classical Hamilton–Jacobi method leads to important geometrical insights like exterior differential systems, Monge cones and time dependent Gaussian metrics. The Wei–Norman approach is applied using unitary transformations defined in terms of generators of the associated Lie groups, here the semi-direct product of the Heisenberg group and the symplectic group. A new explicit relation for the unitary transformations is given in terms of a finite product of elementary transformations. The sequential application of adequate sets of unitary transformations leads naturally to a new separation of variables method for time dependent Hamiltonians, which is shown to be related to the Inönü–Wigner contraction of Lie groups. The new method allows also a better understanding of interacting particles or coupled modes and opens an alternative way to analyze topological phases in driven systems. - Highlights: • Exact unitary transformation reducing time dependent quadratic quantum Hamiltonian to zero. • New separation of variables method and simultaneous uncoupling of modes. • Explicit examples of transformations for one to four dimensional problems. • New general evolution equation for quadratic form in the action, respectively Green’s function.

Time-dependent, many-body scattering theory and nuclear reaction applications

International Nuclear Information System (INIS)

Levin, F.S.

1977-01-01

The channel component state form of the channel coupling array theory of many-body scattering is briefly reviewed. These states obey a non-hermitian matrix equation whose exact solution yields the Schroedinger eigenstates, eigenvalues and scattering amplitudes. A time-dependent formulation of the theory is introduced in analogy to the time-dependent Schrodinger equation and several consequences of the development are noted. These include an interaction picture, a single (matrix) S operator, and the usual connection between the t = 0 time-dependent and the time-independent scattering states. Finally, the channel component states (psi/sub j/) are shown to have the useful property that only psi/sub j/ has (two-body) outgoing waves in channel j: psi/sub m/, m not equal to j, is asymptotically zero in two-body channel j. This formalism is then considered as a means for direct nuclear reaction analysis. Typical bound state approximations are introduced and it is shown that a DWBA amplitude occurs in only one channel. The non-time-reversal invariance of the approximate theory is noted. Results of calculations based on a realistic model for two sets of light-ion induced, one-particle transfer reactions are discussed and compared with the coupled reaction channel (CRC) results using the CRC procedure of Cotanch and Vincent. Angular distributions for the two calculational methods are found to be similar in shape and magnitude. Higher ordercorrections are small as are time-reversal non-invariant effects. Post- and prior-type CRC calculations are seen to differ; the latter are closer to the full CRC results

A Structural Equation Approach to Models with Spatial Dependence

NARCIS (Netherlands)

Oud, Johan H. L.; Folmer, Henk

We introduce the class of structural equation models (SEMs) and corresponding estimation procedures into a spatial dependence framework. SEM allows both latent and observed variables within one and the same (causal) model. Compared with models with observed variables only, this feature makes it

A structural equation approach to models with spatial dependence

NARCIS (Netherlands)

Oud, J.H.L.; Folmer, H.

2008-01-01

We introduce the class of structural equation models (SEMs) and corresponding estimation procedures into a spatial dependence framework. SEM allows both latent and observed variables within one and the same (causal) model. Compared with models with observed variables only, this feature makes it

A Structural Equation Approach to Models with Spatial Dependence

NARCIS (Netherlands)

Oud, J.H.L.; Folmer, H.

2008-01-01

We introduce the class of structural equation models (SEMs) and corresponding estimation procedures into a spatial dependence framework. SEM allows both latent and observed variables within one and the same (causal) model. Compared with models with observed variables only, this feature makes it

Theoretical information measurement in nonrelativistic time-dependent approach

Science.gov (United States)

Najafizade, S. A.; Hassanabadi, H.; Zarrinkamar, S.

2018-02-01

The information-theoretic measures of time-dependent Schrödinger equation are investigated via the Shannon information entropy, variance and local Fisher quantities. In our calculations, we consider the two first states n = 0,1 and obtain the position Sx (t) and momentum Sp (t) Shannon entropies as well as Fisher information Ix (t) in position and momentum Ip (t) spaces. Using the Fourier transformed wave function, we obtain the results in momentum space. Some interesting features of the information entropy densities ρs (x,t) and γs (p,t), as well as the probability densities ρ (x,t) and γ (p,t) for time-dependent states are demonstrated. We establish a general relation between variance and Fisher's information. The Bialynicki-Birula-Mycielski inequality is tested and verified for the states n = 0,1.

On shallow water waves in a medium with time-dependent

Directory of Open Access Journals (Sweden)

Hamdy I. Abdel-Gawad

2015-07-01

Full Text Available In this paper, we studied the progression of shallow water waves relevant to the variable coefficient Korteweg–de Vries (vcKdV equation. We investigated two kinds of cases: when the dispersion and nonlinearity coefficients are proportional, and when they are not linearly dependent. In the first case, it was shown that the progressive waves have some geometric structures as in the case of KdV equation with constant coefficients but the waves travel with time dependent speed. In the second case, the wave structure is maintained when the nonlinearity balances the dispersion. Otherwise, water waves collapse. The objectives of the study are to find a wide class of exact solutions by using the extended unified method and to present a new algorithm for treating the coupled nonlinear PDE’s.

Boundary layer phenomena for differential-delay equations with state-dependent time lags, I.

Science.gov (United States)

Mallet-Paret, John; Nussbaum, Roger D.

1992-11-01

In this paper we begin a study of the differential-delay equation \\varepsilon x'(t) = - x(t) + f(x(t - r)), r = r(x(t)) . We prove the existence of periodic solutions for 0equations. In a companion paper these results will be used to investigate the limiting profile and corresponding boundary layer phenomena for periodic solutions as ɛ approaches zero.

A numerical study of adaptive space and time discretisations for Gross-Pitaevskii equations.

Science.gov (United States)

Thalhammer, Mechthild; Abhau, Jochen

2012-08-15

As a basic principle, benefits of adaptive discretisations are an improved balance between required accuracy and efficiency as well as an enhancement of the reliability of numerical computations. In this work, the capacity of locally adaptive space and time discretisations for the numerical solution of low-dimensional nonlinear Schrödinger equations is investigated. The considered model equation is related to the time-dependent Gross-Pitaevskii equation arising in the description of Bose-Einstein condensates in dilute gases. The performance of the Fourier-pseudo spectral method constrained to uniform meshes versus the locally adaptive finite element method and of higher-order exponential operator splitting methods with variable time stepsizes is studied. Numerical experiments confirm that a local time stepsize control based on a posteriori local error estimators or embedded splitting pairs, respectively, is effective in different situations with an enhancement either in efficiency or reliability. As expected, adaptive time-splitting schemes combined with fast Fourier transform techniques are favourable regarding accuracy and efficiency when applied to Gross-Pitaevskii equations with a defocusing nonlinearity and a mildly varying regular solution. However, the numerical solution of nonlinear Schrödinger equations in the semi-classical regime becomes a demanding task. Due to the highly oscillatory and nonlinear nature of the problem, the spatial mesh size and the time increments need to be of the size of the decisive parameter [Formula: see text], especially when it is desired to capture correctly the quantitative behaviour of the wave function itself. The required high resolution in space constricts the feasibility of numerical computations for both, the Fourier pseudo-spectral and the finite element method. Nevertheless, for smaller parameter values locally adaptive time discretisations facilitate to determine the time stepsizes sufficiently small in order that

Towards a frequency-dependent discrete maximum principle for the implicit Monte Carlo equations

Energy Technology Data Exchange (ETDEWEB)

Wollaber, Allan B [Los Alamos National Laboratory; Larsen, Edward W [Los Alamos National Laboratory; Densmore, Jeffery D [Los Alamos National Laboratory

2010-12-15

It has long been known that temperature solutions of the Implicit Monte Carlo (IMC) equations can exceed the external boundary temperatures, a so-called violation of the 'maximum principle.' Previous attempts at prescribing a maximum value of the time-step size {Delta}{sub t} that is sufficient to eliminate these violations have recommended a {Delta}{sub t} that is typically too small to be used in practice and that appeared to be much too conservative when compared to numerical solutions of the IMC equations for practical problems. In this paper, we derive a new estimator for the maximum time-step size that includes the spatial-grid size {Delta}{sub x}. This explicitly demonstrates that the effect of coarsening {Delta}{sub x} is to reduce the limitation on {Delta}{sub t}, which helps explain the overly conservative nature of the earlier, grid-independent results. We demonstrate that our new time-step restriction is a much more accurate means of predicting violations of the maximum principle. We discuss how the implications of the new, grid-dependent timestep restriction can impact IMC solution algorithms.

Parametric Resonance in a Time-Dependent Harmonic Oscillator

Directory of Open Access Journals (Sweden)

P. N. Nesterov

2013-01-01

Full Text Available In this paper, we study the phenomenon of appearance of new resonances in a timedependent harmonic oscillator under an oscillatory decreasing force. The studied equation belongs to the class of adiabatic oscillators and arises in connection with the spectral problem for the one-dimensional Schr¨odinger equation with Wigner–von Neumann type potential. We use a specially developed method for asymptotic integration of linear systems of differential equations with oscillatory decreasing coefficients. This method uses the ideas of the averaging method to simplify the initial system. Then we apply Levinson’s fundamental theorem to get the asymptotics for its solutions. Finally, we analyze the features of a parametric resonance phenomenon. The resonant frequencies of perturbation are found and the pointwise type of the parametric resonance phenomenon is established. In conclusion, we construct an example of a time-dependent harmonic oscillator (adiabatic oscillator in which the parametric resonances, mentioned in the paper, may occur.

Stability theory for dynamic equations on time scales

CERN Document Server

Martynyuk, Anatoly A

2016-01-01

This monograph is a first in the world to present three approaches for stability analysis of solutions of dynamic equations. The first approach is based on the application of dynamic integral inequalities and the fundamental matrix of solutions of linear approximation of dynamic equations. The second is based on the generalization of the direct Lyapunovs method for equations on time scales, using scalar, vector and matrix-valued auxiliary functions. The third approach is the application of auxiliary functions (scalar, vector, or matrix-valued ones) in combination with differential dynamic inequalities. This is an alternative comparison method, developed for time continuous and time discrete systems. In recent decades, automatic control theory in the study of air- and spacecraft dynamics and in other areas of modern applied mathematics has encountered problems in the analysis of the behavior of solutions of time continuous-discrete linear and/or nonlinear equations of perturbed motion. In the book “Men of Ma...

Time-domain analysis of frequency dependent inertial wave forces on cylinders

DEFF Research Database (Denmark)

Krenk, Steen

2013-01-01

a simple time-domain procedure for the inertial force, in which the frequency dependence is represented via a simple explicit time filter on the wave particle acceleration or velocity. The frequency dependence of the inertia coefficient is known analytically as a function of the wave......-number, and the relevant range of waves shorter than about six times the diameter typically corresponds to deep water waves. This permits a universal non-dimensional frequency representation, that is converted to rational form to provide the relevant filter equation. Simple time-domain simulations demonstrate...... the reduction of the resonant part of the response for natural structural frequencies above the dominating wave frequency....

Empty space-times with separable Hamilton-Jacobi equation

International Nuclear Information System (INIS)

Collinson, C.D.; Fugere, J.

1977-01-01

All empty space-times admitting a one-parameter group of motions and in which the Hamilton-Jacobi equation is (partially) separable are obtained. Several different cases of such empty space-times exist and the Riemann tensor is found to be either type D or N. The results presented here complete the search for empty space-times with separable Hamilton-Jacobi equation. (author)

Applying Boundary Conditions Using a Time-Dependent Lagrangian for Modeling Laser-Plasma Interactions

Science.gov (United States)

Reyes, Jonathan; Shadwick, B. A.

2016-10-01

Modeling the evolution of a short, intense laser pulse propagating through an underdense plasma is of particular interest in the physics of laser-plasma interactions. Numerical models are typically created by first discretizing the equations of motion and then imposing boundary conditions. Using the variational principle of Chen and Sudan, we spatially discretize the Lagrangian density to obtain discrete equations of motion and a discrete energy conservation law which is exactly satisfied regardless of the spatial grid resolution. Modifying the derived equations of motion (e.g., enforcing boundary conditions) generally ruins energy conservation. However, time-dependent terms can be added to the Lagrangian which force the equations of motion to have the desired boundary conditions. Although some foresight is needed to choose these time-dependent terms, this approach provides a mechanism for energy to exit the closed system while allowing the conservation law to account for the loss. An appropriate time discretization scheme is selected based on stability analysis and resolution requirements. We present results using this variational approach in a co-moving coordinate system and compare such results to those using traditional second-order methods. This work was supported by the U. S. Department of Energy under Contract No. DE-SC0008382 and by the National Science Foundation under Contract No. PHY- 1104683.

Equations of motion of test particles for solving the spin-dependent Boltzmann–Vlasov equation

Energy Technology Data Exchange (ETDEWEB)

Xia, Yin [Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800 (China); University of Chinese Academy of Science, Beijing 100049 (China); Xu, Jun, E-mail: xujun@sinap.ac.cn [Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800 (China); Li, Bao-An [Department of Physics and Astronomy, Texas A& M University-Commerce, Commerce, TX 75429-3011 (United States); Department of Applied Physics, Xi' an Jiao Tong University, Xi' an 710049 (China); Shen, Wen-Qing [Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800 (China)

2016-08-10

A consistent derivation of the equations of motion (EOMs) of test particles for solving the spin-dependent Boltzmann–Vlasov equation is presented. The resulting EOMs in phase space are similar to the canonical equations in Hamiltonian dynamics, and the EOM of spin is the same as that in the Heisenburg picture of quantum mechanics. Considering further the quantum nature of spin and choosing the direction of total angular momentum in heavy-ion reactions as a reference of measuring nucleon spin, the EOMs of spin-up and spin-down nucleons are given separately. The key elements affecting the spin dynamics in heavy-ion collisions are identified. The resulting EOMs provide a solid foundation for using the test-particle approach in studying spin dynamics in heavy-ion collisions at intermediate energies. Future comparisons of model simulations with experimental data will help to constrain the poorly known in-medium nucleon spin–orbit coupling relevant for understanding properties of rare isotopes and their astrophysical impacts.

The time-dependent development of electric double-layers in saline solutions

International Nuclear Information System (INIS)

Morrow, R; McKenzie, D R; Bilek, M M M

2006-01-01

We have studied the time-dependent development of electric double-layers (ionic sheaths) in saline solutions by simultaneously solving the sodium and chlorine ion continuity equations coupled with Poisson's equation in one dimension. The study of the effects of time-varying electric fields in solution is relevant to the possible health effect of radio-frequency electric fields on cells in the human body and to assessing the potential of using external electric fields to orient proteins for attachment to surfaces for biosensing applications. Our calculations, for applied voltages of 10-175 mV between the electrode and the solution, predict time scales of ∼0.1-110 μs for the formation of double-layers in solutions of concentration between 0.001 and 1.0 M. We develop an empirical equation that can predict the double-layer formation time to within 10% over this wide parameter range. The method has been validated by comparing the solutions obtained, once the program has run to a steady state, with the standard non-linear Poisson-Boltzmann equations. Excellent agreement is found with the Gouy-Chapman solution of the non-linear Poisson-Boltzmann equation. Thus the method is not restricted in accuracy and applicability as is the case for the linear Poisson-Boltzmann equation. The method can also provide solutions for cases where there are orders of magnitude changes in the ion densities; this has not been the case for previous studies where small perturbation analysis has been employed. The method developed here can readily be extended to two and three dimensions using time-splitting methods

Semianalytic Solution of Space-Time Fractional Diffusion Equation

Directory of Open Access Journals (Sweden)

A. Elsaid

2016-01-01

Full Text Available We study the space-time fractional diffusion equation with spatial Riesz-Feller fractional derivative and Caputo fractional time derivative. The continuation of the solution of this fractional equation to the solution of the corresponding integer order equation is proved. The series solution of this problem is obtained via the optimal homotopy analysis method (OHAM. Numerical simulations are presented to validate the method and to show the effect of changing the fractional derivative parameters on the solution behavior.

A new iterative solver for the time-harmonic wave equation

NARCIS (Netherlands)

Riyanti, C.D.; Erlangga, Y.A.; Plessix, R.E.; Mulder, W.A.; Vuik, C.; Oosterlee, C.

2006-01-01

The time-harmonic wave equation, also known as the Helmholtz equation, is obtained if the constant-density acoustic wave equation is transformed from the time domain to the frequency domain. Its discretization results in a large, sparse, linear system of equations. In two dimensions, this system can

Reactive scattering theory for molecular transitions in time-dependent fields

International Nuclear Information System (INIS)

Peskin, U.; Miller, W.H.

1995-01-01

A new approach is introduced for computing probabilities of molecular transitions in time-dependent fields. The method is based on the stationary (t,t') representation of the Schroedinger equation and is shown to be equivalent to infinite order time-dependent perturbation theory. Bound-to-bound (i.e., photoexcitation) and bound-to-continuum (i.e., photoreaction) transitions are regarded as reactive collisions with the ''time coordinate'' as the reaction coordinate in an extended Hilbert space. A numerical method based on imposing absorbing boundary conditions for the time coordinate in a discrete variable representation framework is introduced. A single operation of the Green's operator provides all the state-specific transition probabilities as well as partial state-resolved (inclusive) reaction probabilities. Illustrative numerical applications are given for model systems

Approximate method for solving the velocity dependent transport equation in a slab lattice

International Nuclear Information System (INIS)

Ferrari, A.

1966-01-01

A method is described that is intended to provide an approximate solution of the transport equation in a medium simulating a water-moderated plate filled reactor core. This medium is constituted by a periodic array of water channels and absorbing plates. The velocity dependent transport equation in slab geometry is included. The computation is performed in a water channel: the absorbing plates are accounted for by the boundary conditions. The scattering of neutrons in water is assumed isotropic, which allows the use of a double Pn approximation to deal with the angular dependence. This method is able to represent the discontinuity of the angular distribution at the channel boundary. The set of equations thus obtained is dependent only on x and v and the coefficients are independent on x. This solution suggests to try solutions involving Legendre polynomials. This scheme leads to a set of equations v dependent only. To obtain an explicit solution, a thermalization model must now be chosen. Using the secondary model of Cadilhac a solution of this set is easy to get. The numerical computations were performed with a particular secondary model, the well-known model of Wigner and Wilkins. (author) [fr

Two-dimensional nonlinear dynamics of an axially moving viscoelastic beam with time-dependent axial speed

International Nuclear Information System (INIS)

Ghayesh, Mergen H.; Amabili, Marco; Farokhi, Hamed

2013-01-01

In the present study, the coupled nonlinear dynamics of an axially moving viscoelastic beam with time-dependent axial speed is investigated employing a numerical technique. The equations of motion for both the transverse and longitudinal motions are obtained using Newton’s second law of motion and the constitutive relations. A two-parameter rheological model of the Kelvin–Voigt energy dissipation mechanism is employed in the modelling of the viscoelastic beam material, in which the material time derivative is used in the viscoelastic constitutive relation. The Galerkin method is then applied to the coupled nonlinear equations, which are in the form of partial differential equations, resulting in a set of nonlinear ordinary differential equations (ODEs) with time-dependent coefficients due to the axial acceleration. A change of variables is then introduced to this set of ODEs to transform them into a set of first-order ordinary differential equations. A variable step-size modified Rosenbrock method is used to conduct direct time integration upon this new set of first-order nonlinear ODEs. The mean axial speed and the amplitude of the speed variations, which are taken as bifurcation parameters, are varied, resulting in the bifurcation diagrams of Poincaré maps of the system. The dynamical characteristics of the system are examined more precisely via plotting time histories, phase-plane portraits, Poincaré sections, and fast Fourier transforms (FFTs)

Integrated vehicle dynamics control using State Dependent Riccati Equations

NARCIS (Netherlands)

Bonsen, B.; Mansvelders, R.; Vermeer, E.

2010-01-01

In this paper we discuss a State Dependent Riccati Equations (SDRE) solution for Integrated Vehicle Dynamics Control (IVDC). The SDRE approach is a nonlinear variant of the well known Linear Quadratic Regulator (LQR) and implements a quadratic cost function optimization. A modified version of this

An efficient explicit marching on in time solver for magnetic field volume integral equation

KAUST Repository

Sayed, Sadeed Bin

2015-07-25

An efficient explicit marching on in time (MOT) scheme for solving the magnetic field volume integral equation is proposed. The MOT system is cast in the form of an ordinary differential equation and is integrated in time using a PE(CE)m multistep scheme. At each time step, a system with a Gram matrix is solved for the predicted/corrected field expansion coefficients. Depending on the type of spatial testing scheme Gram matrix is sparse or consists of blocks with only diagonal entries regardless of the time step size. Consequently, the resulting MOT scheme is more efficient than its implicit counterparts, which call for inversion of fuller matrix system at lower frequencies. Numerical results, which demonstrate the efficiency, accuracy, and stability of the proposed MOT scheme, are presented.

Parsimonious wave-equation travel-time inversion for refraction waves

KAUST Repository

Fu, Lei

2017-02-14

We present a parsimonious wave-equation travel-time inversion technique for refraction waves. A dense virtual refraction dataset can be generated from just two reciprocal shot gathers for the sources at the endpoints of the survey line, with N geophones evenly deployed along the line. These two reciprocal shots contain approximately 2N refraction travel times, which can be spawned into O(N2) refraction travel times by an interferometric transformation. Then, these virtual refraction travel times are used with a source wavelet to create N virtual refraction shot gathers, which are the input data for wave-equation travel-time inversion. Numerical results show that the parsimonious wave-equation travel-time tomogram has about the same accuracy as the tomogram computed by standard wave-equation travel-time inversion. The most significant benefit is that a reciprocal survey is far less time consuming than the standard refraction survey where a source is excited at each geophone location.

Continuous-time random walks on networks with vertex- and time-dependent forcing.

Science.gov (United States)

Angstmann, C N; Donnelly, I C; Henry, B I; Langlands, T A M

2013-08-01

We have investigated the transport of particles moving as random walks on the vertices of a network, subject to vertex- and time-dependent forcing. We have derived the generalized master equations for this transport using continuous time random walks, characterized by jump and waiting time densities, as the underlying stochastic process. The forcing is incorporated through a vertex- and time-dependent bias in the jump densities governing the random walking particles. As a particular case, we consider particle forcing proportional to the concentration of particles on adjacent vertices, analogous to self-chemotactic attraction in a spatial continuum. Our algebraic and numerical studies of this system reveal an interesting pair-aggregation pattern formation in which the steady state is composed of a high concentration of particles on a small number of isolated pairs of adjacent vertices. The steady states do not exhibit this pair aggregation if the transport is random on the vertices, i.e., without forcing. The manifestation of pair aggregation on a transport network may thus be a signature of self-chemotactic-like forcing.

Simultaneous inversion for the space-dependent diffusion coefficient and the fractional order in the time-fractional diffusion equation

International Nuclear Information System (INIS)

Li, Gongsheng; Zhang, Dali; Jia, Xianzheng; Yamamoto, Masahiro

2013-01-01

This paper deals with an inverse problem of simultaneously identifying the space-dependent diffusion coefficient and the fractional order in the 1D time-fractional diffusion equation with smooth initial functions by using boundary measurements. The uniqueness results for the inverse problem are proved on the basis of the inverse eigenvalue problem, and the Lipschitz continuity of the solution operator is established. A modified optimal perturbation algorithm with a regularization parameter chosen by a sigmoid-type function is put forward for the discretization of the minimization problem. Numerical inversions are performed for the diffusion coefficient taking on different functional forms and the additional data having random noise. Several factors which have important influences on the realization of the algorithm are discussed, including the approximate space of the diffusion coefficient, the regularization parameter and the initial iteration. The inversion solutions are good approximations to the exact solutions with stability and adaptivity demonstrating that the optimal perturbation algorithm with the sigmoid-type regularization parameter is efficient for the simultaneous inversion. (paper)

On the XFEL Schrödinger Equation: Highly Oscillatory Magnetic Potentials and Time Averaging

KAUST Repository

Antonelli, Paolo

2014-01-14

We analyse a nonlinear Schrödinger equation for the time-evolution of the wave function of an electron beam, interacting selfconsistently through a Hartree-Fock nonlinearity and through the repulsive Coulomb interaction of an atomic nucleus. The electrons are supposed to move under the action of a time dependent, rapidly periodically oscillating electromagnetic potential. This can be considered a simplified effective single particle model for an X-ray free electron laser. We prove the existence and uniqueness for the Cauchy problem and the convergence of wave-functions to corresponding solutions of a Schrödinger equation with a time-averaged Coulomb potential in the high frequency limit for the oscillations of the electromagnetic potential. © 2014 Springer-Verlag Berlin Heidelberg.

Novel Exponentially Fitted Two-Derivative Runge-Kutta Methods with Equation-Dependent Coefficients for First-Order Differential Equations

Directory of Open Access Journals (Sweden)

Yanping Yang

2016-01-01

Full Text Available The construction of exponentially fitted two-derivative Runge-Kutta (EFTDRK methods for the numerical solution of first-order differential equations is investigated. The revised EFTDRK methods proposed, with equation-dependent coefficients, take into consideration the errors produced in the internal stages to the update. The local truncation errors and stability of the new methods are analyzed. The numerical results are reported to show the accuracy of the new methods.

Universal time-dependence of the mean-square displacement in extremely rugged energy landscapes with equal minima

DEFF Research Database (Denmark)

Dyre, Jeppe; Jacobsen, Jacob M.

1995-01-01

This paper presents a calculation of the time dependence of the mean-square displacement for symmetric random energy barrier hopping models at low temperatures, where the frequency dependence of the normalized diffusion constant D-tilde becomes universal, i.e., independent of the energy barrier...... probability distribution [J. C. Dyre, Phys. Rev. B 49, 11 709 (1994)]. The universal time dependence of the mean-square displacement is calculated from the effective medium approximation (EMA) universality equation, D-tilde lnD-tilde=s-tilde, where s-tilde is the dimensionless imaginary frequency, as well...... as for the approximation to the EMA universality equation D-tilde~=s-tilde/ln(1+s-tilde). At long times the universal mean-square displacement is linear in time, corresponding to ordinary diffusion, whereas the mean-square displacement at short times t in dimensionless units varies as 2/ln(t-1)....

Solution of the one-dimensional time-dependent discrete ordinates problem in a slab by the spectral and LTSN methods

International Nuclear Information System (INIS)

Oliveira, J.V.P. de; Cardona, A.V.; Vilhena, M.T.M.B. de

2002-01-01

In this work, we present a new approach to solve the one-dimensional time-dependent discrete ordinates problem (S N problem) in a slab. The main idea is based upon the application of the spectral method to the set of S N time-dependent differential equations and solution of the resulting coupling equations by the LTS N method. We report numerical simulations

Photodissociation of NaH using time-dependent Fourier grid method

Indian Academy of Sciences (India)

We have solved the time dependent Schrödinger equation by using the Chebyshev polynomial scheme and Fourier grid Hamiltonian method to calculate the dissociation cross section of NaH molecule by 1-photon absorption from the 1+ state to the 1 state. We have found that the results differ significantly from an ...

Electronic structure and time-dependent description of rotational predissociation of LiH

DEFF Research Database (Denmark)

Jasik, P.; Sienkiewicz, J. E.; Domsta, J.

2017-01-01

parameters and the experimental ones. The dynamics of the rotational predissociation process of the 11Π state were studied by solving the time-dependent Schrödinger equation. The classical experiment of Velasco [Can. J. Phys., 1957, 35, 1204] on dissociation in the 11Π state is explained for the first time...

Quaternion wave equations in curved space-time

Science.gov (United States)

Edmonds, J. D., Jr.

1974-01-01

The quaternion formulation of relativistic quantum theory is extended to include curvilinear coordinates and curved space-time in order to provide a framework for a unified quantum/gravity theory. Six basic quaternion fields are identified in curved space-time, the four-vector basis quaternions are identified, and the necessary covariant derivatives are obtained. Invariant field equations are derived, and a general invertable coordinate transformation is developed. The results yield a way of writing quaternion wave equations in curvilinear coordinates and curved space-time as well as a natural framework for solving the problem of second quantization for gravity.

Continuous dependence estimates for viscosity solutions of fully nonlinear degenerate elliptic equations

Directory of Open Access Journals (Sweden)

Espen R. Jakobsen

2002-05-01

Full Text Available Using the maximum principle for semicontinuous functions [3,4], we prove a general ``continuous dependence on the nonlinearities'' estimate for bounded Holder continuous viscosity solutions of fully nonlinear degenerate elliptic equations. Furthermore, we provide existence, uniqueness, and Holder continuity results for bounded viscosity solutions of such equations. Our results are general enough to encompass Hamilton-Jacobi-Bellman-Isaacs's equations of zero-sum, two-player stochastic differential games. An immediate consequence of the results obtained herein is a rate of convergence for the vanishing viscosity method for fully nonlinear degenerate elliptic equations.

The fast algorithm solving the one-dimensional time-dependent Schroedinger equation for teaching purposes

International Nuclear Information System (INIS)

Skoczen, A.; Machowski, W.; Kaprzyk, S.

1990-07-01

Computer program aiming at application in quantum mechanics didactics has been proposed. This program can generate the moving pictures of one-dimensional quantum mechanics scattering phenomena. Constructions of this program provide two options. In the first option the wave packet is generated in infinite one-dimensional well which has walls on the borders of graphic window. In the second option the square potential barrier is located in this well and transmission and reflection of wave packet are shown. We have selected a Gaussian wave packet to represent the initial state of the particle. The wave equation is solved numerically by a method discussed in detail. Solutions for the succesive time moments are graphically presented on the monitor screen. In this way observer can watch whole time-development of physical system. Graphically presented results are physically realistic when program parameters satisfy conditions discussed in this paper. (author)

Tracer kinetics: Modelling by partial differential equations of inhomogeneous compartments with age-dependent elimination rates. Pt. 2

International Nuclear Information System (INIS)

Winkler, E.

1991-01-01

The general theory of inhomogeneous compartments with age-dependent elimination rates is illustrated by examples. Mathematically, it turns out that models consisting of partial differential equations include ordinary, delayed and integro-differential equations, a general fact which is treated here in the context of linear tracer kinetics. The examples include standard compartments as a degenerate case, systems of standard compartments (compartment blocks), models resulting in special residence time distributions, models with pipes, and systems with heterogeneous particles. (orig./BBR) [de

Time-dependent theoretical model of the polar wind: Preliminary results

International Nuclear Information System (INIS)

Gombosi, T.I.; Cravens, T.E.; Nagy, A.F.

1985-01-01

The coupled time dependent continuity, momentum and energy equations of a two ion (O + and H + ) quasineutral plasma were solved in order to extend our understanding of polar wind behavior. This numerical code allows studies of the time dependent behavior of polar wind-type flows into and out of the ionosphere. Initial studies indicate that the typical time constants for electron and ion temperature changes are of the order of minutes and tens of minutes, respectively. The response time of the minor high altitude ion O + is less than an hour, whereas that of the major ion, H + , is many hours. The initial test runs also demonstrate the fact that temporary supersonic flows of both O + and H + are possible, especially in the presence of significant ion heating

Nonlinear diffusion in the presence of a time-dependent external electric field

International Nuclear Information System (INIS)

Lima e Silva, T. de; Galvao, R.M.O.

1987-09-01

The influence of a time-dependent external electric field on the nonlinear diffusion process of weakly ionized plasmas is investigated. A new solution of the diffusion equation is obtained for the case when electron-ion collisions can be neglected. (author) [pt

Time-dependent Hartree-Fock dynamics and phase transition in Lipkin-Meshkov-Glick model

International Nuclear Information System (INIS)

Kan, K.; Lichtner, P.C.; Dworzecka, M.; Griffin, J.J.

1980-01-01

The time-dependent Hartree-Fock solutions of the two-level Lipkin-Meshkov-Glick model are studied by transforming the time-dependent Hartree-Fock equations into Hamilton's canonical form and analyzing the qualitative structure of the Hartree-Fock energy surface in the phase space. It is shown that as the interaction strength increases these time-dependent Hartree-Fock solutions undergo a qualitative change associated with the ground state phase transition previously studied in terms of coherent states. For two-body interactions stronger than the critical value, two types of time-dependent Hartree-Fock solutions (the ''librations'' and ''rotations'' in Hamilton's mechanics) exist simultaneously, while for weaker interactions only the rotations persist. It is also shown that the coherent states with the maximum total pseudospin value are determinants, so that time-dependent Hartree-Fock analysis is equivalent to the coherent state method

Analytical solutions of time-fractional models for homogeneous Gardner equation and non-homogeneous differential equations

Directory of Open Access Journals (Sweden)

Olaniyi Samuel Iyiola

2014-09-01

Full Text Available In this paper, we obtain analytical solutions of homogeneous time-fractional Gardner equation and non-homogeneous time-fractional models (including Buck-master equation using q-Homotopy Analysis Method (q-HAM. Our work displays the elegant nature of the application of q-HAM not only to solve homogeneous non-linear fractional differential equations but also to solve the non-homogeneous fractional differential equations. The presence of the auxiliary parameter h helps in an effective way to obtain better approximation comparable to exact solutions. The fraction-factor in this method gives it an edge over other existing analytical methods for non-linear differential equations. Comparisons are made upon the existence of exact solutions to these models. The analysis shows that our analytical solutions converge very rapidly to the exact solutions.

Quantum features of a charged particle in ionized plasma controlled by a time-dependent magnetic field

Directory of Open Access Journals (Sweden)

Jeong Ryeol eChoi

2014-08-01

Full Text Available Quantum characteristics of a charged particle traveling under the influence of an external time-dependent magnetic field in ionized plasma are investigated using the invariant operator method. The Hamiltonian that gives the radial part of the classical equation of motion for the charged particle is dependent on time. The corresponding invariant operator that satisfies Liouville-von Neumann equation is constructed using fundamental relations. The exact radial wave functions are derived by taking advantage of the eigenstates of the invariant operator. Quantum properties of the system is studied using these wave functions. Especially, the time behavior of the radial component of the quantized energy is addressed in detail.

Coupled kinetic equations for fermions and bosons in the relaxation-time approximation

Science.gov (United States)

Florkowski, Wojciech; Maksymiuk, Ewa; Ryblewski, Radoslaw

2018-02-01

Kinetic equations for fermions and bosons are solved numerically in the relaxation-time approximation for the case of one-dimensional boost-invariant geometry. Fermions are massive and carry baryon number, while bosons are massless. The conservation laws for the baryon number, energy, and momentum lead to two Landau matching conditions, which specify the coupling between the fermionic and bosonic sectors and determine the proper-time dependence of the effective temperature and baryon chemical potential of the system. The numerical results illustrate how a nonequilibrium mixture of fermions and bosons approaches hydrodynamic regime described by the Navier-Stokes equations with appropriate forms of the kinetic coefficients. The shear viscosity of a mixture is the sum of the shear viscosities of fermion and boson components, while the bulk viscosity is given by the formula known for a gas of fermions, however, with the thermodynamic variables characterising the mixture. Thus, we find that massless bosons contribute in a nontrivial way to the bulk viscosity of a mixture, provided fermions are massive. We further observe the hydrodynamization effect, which takes place earlier in the shear sector than in the bulk one. The numerical studies of the ratio of the longitudinal and transverse pressures show, to a good approximation, that it depends on the ratio of the relaxation and proper times only. This behavior is connected with the existence of an attractor solution for conformal systems.

Lagrangian finite element method for 3D time-dependent non-isothermal flow of K-BKZ fluids

DEFF Research Database (Denmark)

Román Marín, José Manuel; Rasmussen, Henrik K.

2009-01-01

equation is replaced with a temperature dependent pseudo time. The spatial coordinate system attached to the particles is discretized by 10-node quadratic tetrahedral elements using Cartesian coordinates. The temperature and the pressure are discretized by 10-node quadratic and linear interpolation...... utilizing an implicit variable step backwards differencing (BDF2) scheme, obtaining second order convergence of the temperature in time. A quadratic interpolation in time is applied to approximate the time integral in the K-BKZ equation. This type of scheme ensures third order accuracy with respect......, respectively, in the tetrahedral particle elements. The spatial discretization of the governing equations follows a mixed Galerkin finite element method. This type of scheme ensures third order accuracy with respect to the discretization of spatial dimension. The temperature equation is solved in time...

Time-domain representation of frequency dependent inertial forces on offshore structures

DEFF Research Database (Denmark)

Krenk, Steen

2013-01-01

dependence is then approximated by a rational function, corresponding to a set of ordinary differential equations in the time domain. The MacCamy-Fuchs solution leads to a representation of the inertial force coefficient as a complex function with argument mainly corresponding to a 'phase lead', in contrast...... history of the inertial force is determined by processing the stable part of the transformation by a forward time integration, followed by an integration in the negative time-direction to obtain the final inertial force time history. The differential equations of the local inertial force at a cross......The inertial wave force on a vertical cylinder decreases with decreasing wave length, when the wave length is less than about six times the diameter of the diameter of the cylinder. In structures with a largediameter component like mono-towers the resonance frequency of the structure is typically...

The Limit Behavior of a Stochastic Logistic Model with Individual Time-Dependent Rates

Directory of Open Access Journals (Sweden)

Yilun Shang

2013-01-01

Full Text Available We investigate a variant of the stochastic logistic model that allows individual variation and time-dependent infection and recovery rates. The model is described as a heterogeneous density dependent Markov chain. We show that the process can be approximated by a deterministic process defined by an integral equation as the population size grows.

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

Unifying treatment of nonequilibrium and unstable dynamics of cold bosonic atom system with time-dependent order parameter in Thermo Field Dynamics

International Nuclear Information System (INIS)

Nakamura, Y.; Yamanaka, Y.

2011-01-01

Research highlights: → Cold atoms with time-dependent condensate in nonequilibrium Thermo Field Dynamics. → Coupled equations which describe the temporal evolution of the system are derived. → They are not the naive assemblages of presumable equations, but the self-consistently ones. → Valid even for systems with Landau or dynamical instability, and describing decays. → Transport equation has new collision term that is important in Landau instability. - Abstract: The coupled equations which describe the temporal evolution of the Bose-Einstein condensed system are derived in the framework of nonequilibrium Thermo Field Dynamics. The key element is that they are not the naive assemblages of assumed equations, but are the self-consistent ones derived by appropriate renormalization conditions. While the order parameter is time-dependent, an explicit quasiparticle picture is constructed by a time-dependent expansion. Our formulation is valid even for the system with a unstable condensate, and describes the condensate decay caused by the Landau instability as well as by the dynamical one.

A discussion of the relativistic equal-time equation

International Nuclear Information System (INIS)

Chengrui, Q.; Danhua, Q.

1981-03-01

Ruan Tu-nan et al have proposed an equal-time equation for composite particles which is derived from Bethe-Salpeter (B-S) equation. Its advantage is that the kernel of this equation is a completely definite single rearrangement of the B-S irreducible kernel without any artificial assumptions. In this paper we shall give a further discussion of the properties of this equation. We discuss the behaviour of this equation as the mass of one of the two particles approaches the limit M 2 → infinite in the ladder approximation of single photon exchange. We show that up to order O(α 4 ) this equation is consistent with the Dirac equation. If the crossed two photon exchange diagrams are taken into account the difference between them is of order O(α 6 ). (author)

Stochastic TDHF and the Boltzman-Langevin equation

International Nuclear Information System (INIS)

Suraud, E.; Reinhard, P.G.

1991-01-01

Outgoing from a time-dependent theory of correlations, we present a stochastic differential equation for the propagation of ensembles of Slater determinants, called Stochastic Time-Dependent Hartree-Fock (Stochastic TDHF). These ensembles are allowed to develop large fluctuations in the Hartree-Fock mean fields. An alternative stochastic differential equation, the Boltzmann-Langevin equation, can be derived from Stochastic TDHF by averaging over subensembles with small fluctuations

The Dirac equation in the Lobachevsky space-time

International Nuclear Information System (INIS)

Paramonov, D.V.; Paramonova, N.N.; Shavokhina, N.S.

2000-01-01

The product of the Lobachevsky space and the time axis is termed the Lobachevsky space-time. The Lobachevsky space is considered as a hyperboloid's sheet in the four-dimensional pseudo-Euclidean space. The Dirac-Fock-Ivanenko equation is reduced to the Dirac equation in two special forms by passing from Lame basis in the Lobachevsky space to the Cartesian basis in the enveloping pseudo-Euclidean space

Time-dependent generalized Gibbs ensembles in open quantum systems

Science.gov (United States)

Lange, Florian; Lenarčič, Zala; Rosch, Achim

2018-04-01

Generalized Gibbs ensembles have been used as powerful tools to describe the steady state of integrable many-particle quantum systems after a sudden change of the Hamiltonian. Here, we demonstrate numerically that they can be used for a much broader class of problems. We consider integrable systems in the presence of weak perturbations which break both integrability and drive the system to a state far from equilibrium. Under these conditions, we show that the steady state and the time evolution on long timescales can be accurately described by a (truncated) generalized Gibbs ensemble with time-dependent Lagrange parameters, determined from simple rate equations. We compare the numerically exact time evolutions of density matrices for small systems with a theory based on block-diagonal density matrices (diagonal ensemble) and a time-dependent generalized Gibbs ensemble containing only a small number of approximately conserved quantities, using the one-dimensional Heisenberg model with perturbations described by Lindblad operators as an example.

Existence of a time-dependent heat flux-related ponderomotive effect

International Nuclear Information System (INIS)

Schamel, H.; Sack, C.

1980-01-01

The existence of a new ponderomotive effect associated with high-frequency waves is pointed out. It originates when time-dependency, mean velocities, or divergent heat fluxes are involved and it supplements the two effects known previously, namely, the ponderomotive force and fake heating. Two proofs are presented; the first is obtained by establishing the momentum equations generalized by including radiation effects and the second by solving the quasi-linear-type diffusion equation explicitly. For a time-dependent wave packet the solution exhibits a new contribution in terms of an integral over previous states. Owing to this term, the plasma has a memory which leads to a breaking of the time symmetry of the plasma response. The range, influenced by the localized wave packet, expands during the course of time due to streamers emanating from the wave active region. Perturbations, among which is the heat flux, are carried to remote positions and, consequently, the region accessible to wave heating is increased. The density dip appears to be less pronounced at the center, and its generation and decay are delayed. The analysis includes a self-consistent action of high-frequency waves as well as the case of traveling wave packets. In order to establish the existence of this new effect, the analytical results are compared with recent microwave experiments. The possibility of generating fast particles by this new ponderomotive effect is emphasized

Cosmologies with a time dependent vacuum

International Nuclear Information System (INIS)

Sola, Joan

2011-01-01

The idea that the cosmological term Λ should be a time dependent quantity in cosmology is a most natural one. It is difficult to conceive an expanding universe with a strictly constant vacuum energy density, ρ Λ = Λ/(8π G), namely one that has remained immutable since the origin of time. A smoothly evolving vacuum energy density ρ Λ = ρ Λ (ξ(t)) that inherits its time-dependence from cosmological functions ξ = ξ(t), such as the Hubble rate H(t) or the scale factor a(t), is not only a qualitatively more plausible and intuitive idea, but is also suggested by fundamental physics, in particular by quantum field theory (QFT) in curved space-time. To implement this notion, is not strictly necessary to resort to ad hoc scalar fields, as usually done in the literature (e.g. in quintessence formulations and the like). A 'running' Λ term can be expected on very similar grounds as one expects (and observes) the running of couplings and masses with a physical energy scale in QFT. Furthermore, the experimental evidence that the equation of state (EOS) of the dark energy (DE) could be evolving with time/redshift (including the possibility that it might currently behave phantom-like) suggests that a time-variable Λ = Λ(t) term (possibly accompanied by a variable Newton's gravitational coupling too, G = G(t)) could account in a natural way for all these features. Remarkably enough, a class of these models (the 'new cosmon') could even be the clue for solving the old cosmological constant problem, including the coincidence problem.

Fractional dynamic calculus and fractional dynamic equations on time scales

CERN Document Server

Georgiev, Svetlin G

2018-01-01

Pedagogically organized, this monograph introduces fractional calculus and fractional dynamic equations on time scales in relation to mathematical physics applications and problems. Beginning with the definitions of forward and backward jump operators, the book builds from Stefan Hilger’s basic theories on time scales and examines recent developments within the field of fractional calculus and fractional equations. Useful tools are provided for solving differential and integral equations as well as various problems involving special functions of mathematical physics and their extensions and generalizations in one and more variables. Much discussion is devoted to Riemann-Liouville fractional dynamic equations and Caputo fractional dynamic equations.  Intended for use in the field and designed for students without an extensive mathematical background, this book is suitable for graduate courses and researchers looking for an introduction to fractional dynamic calculus and equations on time scales. .

Time-dependent coupled harmonic oscillators: classical and quantum solutions

International Nuclear Information System (INIS)

Macedo, D.X.; Guedes, I.

2014-01-01

In this work we present the classical and quantum solutions for an arbitrary system of time-dependent coupled harmonic oscillators, where the masses (m), frequencies (ω) and coupling parameter (k) are functions of time. To obtain the classical solutions, we use a coordinate and momentum transformations along with a canonical transformation to write the original Hamiltonian as the sum of two Hamiltonians of uncoupled harmonic oscillators with modified time-dependent frequencies and unitary masses. To obtain the exact quantum solutions we use a unitary transformation and the Lewis and Riesenfeld (LR) invariant method. The exact wave functions are obtained by solving the respective Milne–Pinney (MP) equation for each system. We obtain the solutions for the system with m 1 = m 2 = m 0 e γt , ω 1 = ω 01 e -γt/2 , ω 2 = ω 02 e -γt/2 and k = k 0 . (author)

Time-dependent local-to-normal mode transition in triatomic molecules

Science.gov (United States)

Cruz, Hans; Bermúdez-Montaña, Marisol; Lemus, Renato

2018-01-01

Time-evolution of the vibrational states of two interacting harmonic oscillators in the local mode scheme is presented. A local-to-normal mode transition (LNT) is identified and studied from temporal perspective through time-dependent frequencies of the oscillators. The LNT is established as a polyad-breaking phenomenon from the local standpoint for the stretching degrees of freedom in a triatomic molecule. This study is carried out in the algebraic representation of bosonic operators. The dynamics of the states are determined via the solutions of the corresponding nonlinear Ermakov equation and a local time-dependent polyad is obtained as a tool to identify the LNT. Applications of this formalism to H2O, CO2, O3 and NO2 molecules in the adiabatic, sudden and linear regime are considered.

Computational micromagnetics: prediction of time dependent and thermal properties

International Nuclear Information System (INIS)

Schrefl, T.; Scholz, W.; Suess, Dieter; Fidler, J.

2001-01-01

Finite element modeling treats magnetization processes on a length scale of several nanometers and thus gives a quantitative correlation between the microstructure and the magnetic properties of ferromagnetic materials. This work presents a novel finite element/boundary element micro-magnetics solver that combines a wavelet-based matrix compression technique for magnetostatic field calculations with a BDF/GMRES method for the time integration of the Gilbert equation of motion. The simulations show that metastable energy minima and nonuniform magnetic states within the grains are important factors in the reversal dynamics at finite temperature. The numerical solution of the Gilbert equation shows how reversed domains nucleate and expand. The switching time of submicron magnetic elements depends on the shape of the elements. Elements with slanted ends decrease the overall reversal time, as a transverse demagnetizing field suppresses oscillations of the magnetization. Thermal activated processes can be included adding a random thermal field to the effective magnetic field. Thermally assisted reversal was studied for CoCrPtTa thin-film media

How a dependent's variable non-randomness affects taper equation ...

African Journals Online (AJOL)

In order to apply the least squares method in regression analysis, the values of the dependent variable Y should be random. In an example of regression analysis linear and nonlinear taper equations, which estimate the diameter of the tree dhi at any height of the tree hi, were compared. For each tree the diameter at the ...

International Conference in Memory of Professor Kôsaku Yosida held at RIMS

CERN Document Server

1993-01-01

In these proceedings of the international conference held in Kyoto in memoryof the late Professor K saku Yosida, twenty six invited speakers display in their many facets of functional analysis and its applications in the research tradition of Yosida's school. Many of the topics are related tolinear and non-linear partial differential equations, including the Schr|dinger equations, the Navier-Stokes equations and quasilinear hyperbolic equations. Several of the papers are survey articles, the others are original (unpublished) and refereed research articles. Also included is a full listing of the publications of K. Yosida. Recommendedto students and research workers looking for a bird's-eye view of current research activity in functional analysis and its applications. FROM THE CONTENTS: K. Ito: Semigroups in probability theory.- T. Kato: Abstract evolution equations, linear and quasilinear, revisited.- J.L. Lions: Remarkson systems with incompletely given initial data and incompletely given part of the boundary...

Josephson-like currents in graphene for arbitrary time-dependent potential barriers

OpenAIRE

Savel'ev, Sergey E.; Hausler, Wolfgang; Hanggi, Peter

2011-01-01

From the exact solution of the Dirac-Weyl equation we find unusual currents j_y running in y-direction parallel to a time-dependent scalar potential barrier W(x,t) placed upon a monolayer of graphene, even for vanishing momentum component p_y. In their sine-like dependence on the phase difference of wave functions, describing left and right moving Dirac fermions, these currents resemble Josephson currents in superconductors, including the occurance of Shapiro steps at certain frequencies of p...

Time-Dependent Close-Coupling Methods for Electron-Atom/Molecule Scattering

International Nuclear Information System (INIS)

Colgan, James

2014-01-01

The time-dependent close-coupling (TDCC) method centers on an accurate representation of the interaction between two outgoing electrons moving in the presence of a Coulomb field. It has been extensively applied to many problems of electrons, photons, and ions scattering from light atomic targets. Theoretical Description: The TDCC method centers on a solution of the time-dependent Schrödinger equation for two interacting electrons. The advantages of a time-dependent approach are two-fold; one treats the electron-electron interaction essentially in an exact manner (within numerical accuracy) and a time-dependent approach avoids the difficult boundary condition encountered when two free electrons move in a Coulomb field (the classic three-body Coulomb problem). The TDCC method has been applied to many fundamental atomic collision processes, including photon-, electron- and ion-impact ionization of light atoms. For application to electron-impact ionization of atomic systems, one decomposes the two-electron wavefunction in a partial wave expansion and represents the subsequent two-electron radial wavefunctions on a numerical lattice. The number of partial waves required to converge the ionization process depends on the energy of the incoming electron wavepacket and on the ionization threshold of the target atom or ion.

Experimental test of depth dependence of solutions for time-resolved diffusion equation

Energy Technology Data Exchange (ETDEWEB)

Laidevant, A.; Da Silva, A.; Moy, J.P.; Berger, M.; Dinten, J.M

2004-07-01

The determination of optical properties of a semi-infinite medium such as biological tissue has been widely investigated by many authors. Reflectance formulas can be derived from the diffusion equation for different boundary conditions at the medium-air interface. This quantity can be measured at the medium surface. For realistic objects, such as a mouse, tissue optical properties can realistically only be determined at the object surface. However, near the surface diffusion approximation is weak and boundary models have to be considered. In order to investigate the validity of the time resolved reflectance approach at the object boundary, we have estimated optical properties of a liquid semi-infinite medium by this method for different boundary conditions and different fiber's position beneath the surface. The time-correlated single photon counting (TCSPC) technique is used to measure the reflectance curve. Our liquid phantoms are made of water, Intra-lipid and Ink. Laser light is delivered by a pulsed laser diode. Measurements are then fitted to theoretical solutions expressed as a function of source and detector's depth and distance. By taking as reference the optical properties obtained from the infinite model for fibers deeply immersed, influence of the different boundary conditions and bias induced are established for different fibers' depth and a variety of solutions. This influence is analysed by comparing evolution of the reflectance models, as well as estimations of absorption and scattering coefficients. According to this study we propose a strategy for determining optical properties of a solid phantom where measurements can only be realized at the surface. (authors)

Time-dependent patterns in quasivertical cylindrical binary convection

Science.gov (United States)

Alonso, Arantxa; Mercader, Isabel; Batiste, Oriol

2018-02-01

This paper reports on numerical investigations of the effect of a slight inclination α on pattern formation in a shallow vertical cylindrical cell heated from below for binary mixtures with a positive value of the Soret coefficient. By using direct numerical simulation of the three-dimensional Boussinesq equations with Soret effect in cylindrical geometry, we show that a slight inclination of the cell in the range α ≈0.036 rad =2∘ strongly influences pattern selection. The large-scale shear flow (LSSF) induced by the small tilt of gravity overcomes the squarelike arrangements observed in noninclined cylinders in the Soret regime, stratifies the fluid along the direction of inclination, and produces an enhanced separation of the two components of the mixture. The competition between shear effects and horizontal and vertical buoyancy alters significantly the dynamics observed in noninclined convection. Additional unexpected time-dependent patterns coexist with the basic LSSF. We focus on an unsual periodic state recently discovered in an experiment, the so-called superhighway convection state (SHC), in which ascending and descending regions of fluid move in opposite directions. We provide numerical confirmation that Boussinesq Navier-Stokes equations with standard boundary conditions contain the essential ingredients that allow for the existence of such a state. Also, we obtain a persistent heteroclinic structure where regular oscillations between a SHC pattern and a state of nearly stationary longitudinal rolls take place. We characterize numerically these time-dependent patterns and investigate the dynamics around the threshold of convection.

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.

Angle-dependent strong-field molecular ionization rates with tuned range-separated time-dependent density functional theory

Energy Technology Data Exchange (ETDEWEB)

Sissay, Adonay [Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803 (United States); Abanador, Paul; Mauger, François; Gaarde, Mette; Schafer, Kenneth J. [Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803 (United States); Lopata, Kenneth, E-mail: klopata@lsu.edu [Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803 (United States); Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana 70803 (United States)

2016-09-07

Strong-field ionization and the resulting electronic dynamics are important for a range of processes such as high harmonic generation, photodamage, charge resonance enhanced ionization, and ionization-triggered charge migration. Modeling ionization dynamics in molecular systems from first-principles can be challenging due to the large spatial extent of the wavefunction which stresses the accuracy of basis sets, and the intense fields which require non-perturbative time-dependent electronic structure methods. In this paper, we develop a time-dependent density functional theory approach which uses a Gaussian-type orbital (GTO) basis set to capture strong-field ionization rates and dynamics in atoms and small molecules. This involves propagating the electronic density matrix in time with a time-dependent laser potential and a spatial non-Hermitian complex absorbing potential which is projected onto an atom-centered basis set to remove ionized charge from the simulation. For the density functional theory (DFT) functional we use a tuned range-separated functional LC-PBE*, which has the correct asymptotic 1/r form of the potential and a reduced delocalization error compared to traditional DFT functionals. Ionization rates are computed for hydrogen, molecular nitrogen, and iodoacetylene under various field frequencies, intensities, and polarizations (angle-dependent ionization), and the results are shown to quantitatively agree with time-dependent Schrödinger equation and strong-field approximation calculations. This tuned DFT with GTO method opens the door to predictive all-electron time-dependent density functional theory simulations of ionization and ionization-triggered dynamics in molecular systems using tuned range-separated hybrid functionals.

Angle-dependent strong-field molecular ionization rates with tuned range-separated time-dependent density functional theory

International Nuclear Information System (INIS)

Sissay, Adonay; Abanador, Paul; Mauger, François; Gaarde, Mette; Schafer, Kenneth J.; Lopata, Kenneth

2016-01-01

Strong-field ionization and the resulting electronic dynamics are important for a range of processes such as high harmonic generation, photodamage, charge resonance enhanced ionization, and ionization-triggered charge migration. Modeling ionization dynamics in molecular systems from first-principles can be challenging due to the large spatial extent of the wavefunction which stresses the accuracy of basis sets, and the intense fields which require non-perturbative time-dependent electronic structure methods. In this paper, we develop a time-dependent density functional theory approach which uses a Gaussian-type orbital (GTO) basis set to capture strong-field ionization rates and dynamics in atoms and small molecules. This involves propagating the electronic density matrix in time with a time-dependent laser potential and a spatial non-Hermitian complex absorbing potential which is projected onto an atom-centered basis set to remove ionized charge from the simulation. For the density functional theory (DFT) functional we use a tuned range-separated functional LC-PBE*, which has the correct asymptotic 1/r form of the potential and a reduced delocalization error compared to traditional DFT functionals. Ionization rates are computed for hydrogen, molecular nitrogen, and iodoacetylene under various field frequencies, intensities, and polarizations (angle-dependent ionization), and the results are shown to quantitatively agree with time-dependent Schrödinger equation and strong-field approximation calculations. This tuned DFT with GTO method opens the door to predictive all-electron time-dependent density functional theory simulations of ionization and ionization-triggered dynamics in molecular systems using tuned range-separated hybrid functionals.

Critique of the foundations of time-dependent density-functional theory

International Nuclear Information System (INIS)

Schirmer, J.; Dreuw, A.

2007-01-01

The general expectation that, in principle, the time-dependent density-functional theory (TDDFT) is an exact formulation of the time evolution of an interacting N-electron system is critically reexamined. It is demonstrated that the previous TDDFT foundation, resting on four theorems by Runge and Gross (RG) [Phys. Rev. Lett. 52, 997 (1984)], is invalid because undefined phase factors corrupt the RG action integral functionals. Our finding confirms much of a previous analysis by van Leeuwen [Int. J. Mod. Phys. B 15, 1969 (2001)]. To analyze the RG theorems and other aspects of TDDFT, an utmost simplification of the Kohn-Sham (KS) concept has been introduced, in which the ground-state density is obtained from a single KS equation for one spatial (spinless) orbital. The time-dependent (TD) form of this radical Kohn-Sham (rKS) scheme, which has the same validity status as the ordinary KS version, has proved to be a valuable tool for analysis. The rKS concept is used to clarify also the alternative nonvariational formulation of TD KS theory. We argue that it is just a formal theory, allowing one to reproduce but not predict the time development of the exact density of the interacting N-electron system. Besides the issue of the formal exactness of TDDFT, it is shown that both the static and time-dependent KS linear response equations neglect the particle-particle (p-p) and hole-hole (h-h) matrix elements of the perturbing operator. For a local (multiplicative) operator this does not lead to a loss of information due to a remarkable general property of local operators. Accordingly, no logical inconsistency arises with respect to DFT, because DFT requires any external potential to be local. For a general nonlocal operator the error resulting from the neglected matrix elements is of second order in the electronic repulsion

Similarity Solutions for Multiterm Time-Fractional Diffusion Equation

OpenAIRE

Elsaid, A.; Abdel Latif, M. S.; Maneea, M.

2016-01-01

Similarity method is employed to solve multiterm time-fractional diffusion equation. The orders of the fractional derivatives belong to the interval (0,1] and are defined in the Caputo sense. We illustrate how the problem is reduced from a multiterm two-variable fractional partial differential equation to a multiterm ordinary fractional differential equation. Power series solution is obtained for the resulting ordinary problem and the convergence of the series solution is discussed. Based on ...

Time-Dependent Stress Rupture Strength Degradation of Hi-Nicalon Fiber-Reinforced Silicon Carbide Composites at Intermediate Temperatures

Science.gov (United States)

Sullivan, Roy M.

2016-01-01

The stress rupture strength of silicon carbide fiber-reinforced silicon carbide composites with a boron nitride fiber coating decreases with time within the intermediate temperature range of 700 to 950 degree Celsius. Various theories have been proposed to explain the cause of the time-dependent stress rupture strength. The objective of this paper is to investigate the relative significance of the various theories for the time-dependent strength of silicon carbide fiber-reinforced silicon carbide composites. This is achieved through the development of a numerically based progressive failure analysis routine and through the application of the routine to simulate the composite stress rupture tests. The progressive failure routine is a time-marching routine with an iterative loop between a probability of fiber survival equation and a force equilibrium equation within each time step. Failure of the composite is assumed to initiate near a matrix crack and the progression of fiber failures occurs by global load sharing. The probability of survival equation is derived from consideration of the strength of ceramic fibers with randomly occurring and slow growing flaws as well as the mechanical interaction between the fibers and matrix near a matrix crack. The force equilibrium equation follows from the global load sharing presumption. The results of progressive failure analyses of the composite tests suggest that the relationship between time and stress-rupture strength is attributed almost entirely to the slow flaw growth within the fibers. Although other mechanisms may be present, they appear to have only a minor influence on the observed time-dependent behavior.

Concept of frequency separation in life prediction for time-dependent fatigue

International Nuclear Information System (INIS)

Coffin, L.F.

1976-01-01

Two methods are described to improve the predictive capability of the frequency-modified fatigue equations for time-dependent fatigue. These built-on the earlier approach and are applicable to severely unbalanced hysteresis loops. Comparisons made with various wave-shape investigations show favorable results. A new form of hysteresis loop is introduced utilizing frequency separation concepts. 4 tables, 11 fig

Time dependence of entropy flux and entropy production for a dynamical system driven by noises with coloured cross-correlation

Institute of Scientific and Technical Information of China (English)

Xie Wen-Xian; Xu Wei; Cai Li

2007-01-01

This paper shows the Fokker-Planck equation of a dynamical system driven by coloured cross-correlated white noises in the absence and presence of a small external force. Based on the Fokker-Planck equation and the definition of Shannon's information entropy, the time dependence of entropy flux and entropy production can be calculated. The present results can be used to explain the extremal behaviour of time dependence of entropy flux and entropy production in view of the dissipative parameter γ of the system, coloured cross-correlation time τ and coloured cross-correlation strength λ.

Time-dependent density functional theory of open quantum systems in the linear-response regime.

Science.gov (United States)

Tempel, David G; Watson, Mark A; Olivares-Amaya, Roberto; Aspuru-Guzik, Alán

2011-02-21

Time-dependent density functional theory (TDDFT) has recently been extended to describe many-body open quantum systems evolving under nonunitary dynamics according to a quantum master equation. In the master equation approach, electronic excitation spectra are broadened and shifted due to relaxation and dephasing of the electronic degrees of freedom by the surrounding environment. In this paper, we develop a formulation of TDDFT linear-response theory (LR-TDDFT) for many-body electronic systems evolving under a master equation, yielding broadened excitation spectra. This is done by mapping an interacting open quantum system onto a noninteracting open Kohn-Sham system yielding the correct nonequilibrium density evolution. A pseudoeigenvalue equation analogous to the Casida equations of the usual LR-TDDFT is derived for the Redfield master equation, yielding complex energies and Lamb shifts. As a simple demonstration, we calculate the spectrum of a C(2 +) atom including natural linewidths, by treating the electromagnetic field vacuum as a photon bath. The performance of an adiabatic exchange-correlation kernel is analyzed and a first-order frequency-dependent correction to the bare Kohn-Sham linewidth based on the Görling-Levy perturbation theory is calculated.

Dealing with Dependent Uncertainty in Modelling: A Comparative Study Case through the Airy Equation

Directory of Open Access Journals (Sweden)

J.-C. Cortés

2013-01-01

Full Text Available The consideration of uncertainty in differential equations leads to the emergent area of random differential equations. Under this approach, inputs become random variables and/or stochastic processes. Often one assumes that inputs are independent, a hypothesis that simplifies the mathematical treatment although it could not be met in applications. In this paper, we analyse, through the Airy equation, the influence of statistical dependence of inputs on the output, computing its expectation and standard deviation by Fröbenius and Polynomial Chaos methods. The results are compared with Monte Carlo sampling. The analysis is conducted by the Airy equation since, as in the deterministic scenario its solutions are highly oscillatory, it is expected that differences will be better highlighted. To illustrate our study, and motivated by the ubiquity of Gaussian random variables in numerous practical problems, we assume that inputs follow a multivariate Gaussian distribution throughout the paper. The application of Fröbenius method to solve Airy equation is based on an extension of the method to the case where inputs are dependent. The numerical results show that the existence of statistical dependence among the inputs and its magnitude entails changes on the variability of the output.

Parallel processing for a 1-D time-dependent solution to impurity rate equations for fusion plasma simulations

International Nuclear Information System (INIS)

Veerasingam, R.

1990-01-01

In fusion plasmas impurities such as carbon, oxygen or nickel can contaminate the plasma and cause degradation of the performance of a fusion device through radiation. However, impurities can also be used as diagnostics to obtain information about a plasma through spectroscopic experiments which can then be used in plasma modeling and simulations. In the past, serial algorithms have been described for either the time dependent or steady state problem. In this paper, we describe a parallel procedure adopted to solve the time-dependent problem. It can be shown that for the steady state problem a parallel procedure would not be a useful application of parallelization because a few seconds of the Central Processing Unit time on a CRAY-XMP or IBM 3090/600S would suffice to obtain the solution, while this is not the case for the time-dependent problem. In order to study the effects of low Z and high Z impurities on the final state of a plasma, time-dependent solutions are necessary. For purposes of diagnostics and comparisons with experiments, a fast turn around time of the simulations would be advantageous. We have implemented a parallel algorithm on and IBM 3090/600S and tested its performance for a typical set of fusion plasma parameters. 4 refs., 1 tab

Light pressure of time-dependent fields in plasmas

International Nuclear Information System (INIS)

Zeidler, A.; Schnabl, H.; Mulser, P.

1985-01-01

An expression of the light pressure Pi is derived for the case of a nearly monochromatic electromagnetic wave with arbitrarily time-dependent amplitude. Thereby Pi is defined as the time-averaged force density exerted on a plasma by the wave. The resulting equations are valid for both transverse and longitudinal waves. The light pressure turns out to consist of two components: the well-known gradient-type term and a new nonstationary solenoidal term. This is true for warm as well as cold plasmas. The importance of the new term for the generation of static magnetic fields is shown, and a model in which shear forces may result is given. Formulas for the nonstationary light pressure developed previously are discussed

Parallel solution of the time-dependent Ginzburg-Landau equations and other experiences using BlockComm-Chameleon and PCN on the IBM SP, Intel iPSC/860, and clusters of workstations

International Nuclear Information System (INIS)

Coskun, E.

1995-09-01

Time-dependent Ginzburg-Landau (TDGL) equations are considered for modeling a thin-film finite size superconductor placed under magnetic field. The problem then leads to the use of so-called natural boundary conditions. Computational domain is partitioned into subdomains and bond variables are used in obtaining the corresponding discrete system of equations. An efficient time-differencing method based on the Forward Euler method is developed. Finally, a variable strength magnetic field resulting in a vortex motion in Type II High T c superconducting films is introduced. The authors tackled the problem using two different state-of-the-art parallel computing tools: BlockComm/Chameleon and PCN. They had access to two high-performance distributed memory supercomputers: the Intel iPSC/860 and IBM SP1. They also tested the codes using, as a parallel computing environment, a cluster of Sun Sparc workstations

Numerical resolution of Navier-Stokes equations coupled to the heat equation

International Nuclear Information System (INIS)

Zenouda, Jean-Claude

1970-08-01

The author proves a uniqueness theorem for the time dependent Navier-Stokes equations coupled with heat flow in the two-dimensional case. He studies stability and convergence of several finite - difference schemes to solve these equations. Numerical experiments are done in the case of a square domain. (author) [fr

Advanced Time Approach of FW-H Equations for Predicting Noise

DEFF Research Database (Denmark)

Haiqing, Si; Yan, Shi; Shen, Wen Zhong

2013-01-01

An advanced time approach of Ffowcs Williams-Hawkings (FW-H) acoustic analogy is developed, and the integral equations and integral solution of FW-H acoustic analogy are derived. Compared with the retarded time approach, the transcendental equation need not to be solved in the advanced time...

Regularity of the Rotation Number for the One-Dimensional Time-Continuous Schroedinger Equation

Energy Technology Data Exchange (ETDEWEB)

Amor, Sana Hadj, E-mail: sana_hadjamor@yahoo.fr [Ecole Nationale d' Ingenieurs de Monastir (Tunisia)

2012-12-15

Starting from results already obtained for quasi-periodic co-cycles in SL(2, R), we show that the rotation number of the one-dimensional time-continuous Schroedinger equation with Diophantine frequencies and a small analytic potential has the behavior of a 1/2-Hoelder function. We give also a sub-exponential estimate of the length of the gaps which depends on its label given by the gap-labeling theorem.

General exact solution for homogeneous time-dependent self-gravitating perfect fluids

International Nuclear Information System (INIS)

Gaete, P.; Hojman, R.

1988-01-01

A procedure to obtain the general exact solution of Einstein equations for a self-gravitating spherically-symmetric static perfect fluid obeying an arbitrary equation of state, is applied to time-dependent Kantowsky-Sachs line elements (with spherical, planar and hyperbolic symmetry). As in the static case, the solution is generated by an arbitrary function of the independent variable and its first derivative. To illustrate the results, the whole family of (plane-symmetric) solutions with a ''gamma-law'' equation of state is explicity obtained in terms of simple known functions. It is also shown that, while in the static plane-symmtric line elements, every metric is in one to one correspondence with a ''partner-metric'' (both originated from the same generatrix function), in this case every generatrix function univocally determines one metric. (author) [pt

Exact solutions of space-time fractional EW and modified EW equations

International Nuclear Information System (INIS)

Korkmaz, Alper

2017-01-01

The bright soliton solutions and singular solutions are constructed for the space-time fractional EW and the space-time fractional modified EW (MEW) equations. Both equations are reduced to ordinary differential equations by the use of fractional complex transform (FCT) and properties of modified Riemann–Liouville derivative. Then, various ansatz method are implemented to construct the solutions for both equations.

Operator of Time and Generalized Schrödinger Equation

Directory of Open Access Journals (Sweden)

Slobodan Prvanović

2018-01-01

Full Text Available The equation describing the change of the state of the quantum system with respect to energy is introduced within the framework of the self-adjoint operator of time in nonrelativistic quantum mechanics. In this proposal, the operator of time appears to be the generator of the change of the energy, while the operator of energy that is conjugate to the operator of time generates the time evolution. Two examples, one with discrete time and the other with continuous one, are given and the generalization of Schrödinger equation is proposed.

Wave-packet revival for the Schroedinger equation with position-dependent mass

International Nuclear Information System (INIS)

Schmidt, Alexandre G.M.

2006-01-01

We study the temporal evolution of solutions of 1D Schroedinger equation with position-dependent mass inside an infinite well. Revival of wave-packet is shown to exist and partial revivals are different from the usual ones

Flavored quantum Boltzmann equations

International Nuclear Information System (INIS)

Cirigliano, Vincenzo; Lee, Christopher; Ramsey-Musolf, Michael J.; Tulin, Sean

2010-01-01

We derive from first principles, using nonequilibrium field theory, the quantum Boltzmann equations that describe the dynamics of flavor oscillations, collisions, and a time-dependent mass matrix in the early universe. Working to leading nontrivial order in ratios of relevant time scales, we study in detail a toy model for weak-scale baryogenesis: two scalar species that mix through a slowly varying time-dependent and CP-violating mass matrix, and interact with a thermal bath. This model clearly illustrates how the CP asymmetry arises through coherent flavor oscillations in a nontrivial background. We solve the Boltzmann equations numerically for the density matrices, investigating the impact of collisions in various regimes.

Pomarning-eddington approximation for time-dependent radiation transfer in finite slab media

International Nuclear Information System (INIS)

El-Wakil, S.A.; Degheidy, A.R.; Sallah, M.

2005-01-01

The time-dependent monoenergetic radiation transfer equation with linear anisotropic scattering is proposed. Pomraning-Eddington approximation is used to calculate the radiation intensity in finite plane-parallel media. Numerical results are done for the isotropic media. Shielding calculations are shown for reflectivity and transmissivity at different times. The medium is assumed to have specular-reflecting boundaries. Two different weight functions are introduced to force the boundary conditions to be fulfilled

Solving the Schroedinger equation using the finite difference time domain method

International Nuclear Information System (INIS)

Sudiarta, I Wayan; Geldart, D J Wallace

2007-01-01

In this paper, we solve the Schroedinger equation using the finite difference time domain (FDTD) method to determine energies and eigenfunctions. In order to apply the FDTD method, the Schroedinger equation is first transformed into a diffusion equation by the imaginary time transformation. The resulting time-domain diffusion equation is then solved numerically by the FDTD method. The theory and an algorithm are provided for the procedure. Numerical results are given for illustrative examples in one, two and three dimensions. It is shown that the FDTD method accurately determines eigenfunctions and energies of these systems

Time-Dependent Damage Investigation of Rock Mass in an In Situ Experimental Tunnel

Science.gov (United States)

Jiang, Quan; Cui, Jie; Chen, Jing

2012-01-01

In underground tunnels or caverns, time-dependent deformation or failure of rock mass, such as extending cracks, gradual rock falls, etc., are a costly irritant and a major safety concern if the time-dependent damage of surrounding rock is serious. To understand the damage evolution of rock mass in underground engineering, an in situ experimental testing was carried out in a large belowground tunnel with a scale of 28.5 m in width, 21 m in height and 352 m in length. The time-dependent damage of rock mass was detected in succession by an ultrasonic wave test after excavation. The testing results showed that the time-dependent damage of rock mass could last a long time, i.e., nearly 30 days. Regression analysis of damage factors defined by wave velocity, resulted in the time-dependent evolutional damage equation of rock mass, which corresponded with logarithmic format. A damage viscoelastic-plastic model was developed to describe the exposed time-dependent deterioration of rock mass by field test, such as convergence of time-dependent damage, deterioration of elastic modules and logarithmic format of damage factor. Furthermore, the remedial measures for damaged surrounding rock were discussed based on the measured results and the conception of damage compensation, which provides new clues for underground engineering design.

Exact solutions to the time-fractional differential equations via local fractional derivatives

Science.gov (United States)

Guner, Ozkan; Bekir, Ahmet

2018-01-01

This article utilizes the local fractional derivative and the exp-function method to construct the exact solutions of nonlinear time-fractional differential equations (FDEs). For illustrating the validity of the method, it is applied to the time-fractional Camassa-Holm equation and the time-fractional-generalized fifth-order KdV equation. Moreover, the exact solutions are obtained for the equations which are formed by different parameter values related to the time-fractional-generalized fifth-order KdV equation. This method is an reliable and efficient mathematical tool for solving FDEs and it can be applied to other non-linear FDEs.

Time-dependent solutions for stochastic systems with delays: Perturbation theory and applications to financial physics

International Nuclear Information System (INIS)

Frank, T.D.

2006-01-01

First-order approximations of time-dependent solutions are determined for stochastic systems perturbed by time-delayed feedback forces. To this end, the theory of delay Fokker-Planck equations is applied in combination with Bayes' theorem. Applications to a time-delayed Ornstein-Uhlenbeck process and the geometric Brownian walk of financial physics are discussed

Geometric Implications of Maxwell's Equations

Science.gov (United States)

Smith, Felix T.

2015-03-01

Maxwell's synthesis of the varied results of the accumulated knowledge of electricity and magnetism, based largely on the searching insights of Faraday, still provide new issues to explore. A case in point is a well recognized anomaly in the Maxwell equations: The laws of electricity and magnetism require two 3-vector and two scalar equations, but only six dependent variables are available to be their solutions, the 3-vectors E and B. This leaves an apparent redundancy of two degrees of freedom (J. Rosen, AJP 48, 1071 (1980); Jiang, Wu, Povinelli, J. Comp. Phys. 125, 104 (1996)). The observed self-consistency of the eight equations suggests that they contain additional information. This can be sought as a previously unnoticed constraint connecting the space and time variables, r and t. This constraint can be identified. It distorts the otherwise Euclidean 3-space of r with the extremely slight, time dependent curvature k (t) =Rcurv-2 (t) of the 3-space of a hypersphere whose radius has the time dependence dRcurv / dt = +/- c nonrelativistically, or dRcurvLor / dt = +/- ic relativistically. The time dependence is exactly that of the Hubble expansion. Implications of this identification will be explored.

The adaptive CCCG({eta}) method for efficient solution of time dependent partial differential equations

Energy Technology Data Exchange (ETDEWEB)

Campos, F.F. [Universidade Federal de Minas Gerais, Belo Horizonte (Brazil); Birkett, N.R.C. [Oxford Univ. Computing Lab. (United Kingdom)

1996-12-31

The Controlled Cholesky factorisation has been shown to be a robust preconditioner for the Conjugate Gradient method. In this scheme the amount of fill-in is defined in terms of a parameter {eta}, the number of extra elements allowed per column. It is demonstrated how an optimum value of {eta} can be automatically determined when solving time dependent p.d.e.`s using an implicit time step method. A comparison between CCCG({eta}) and the standard ICCG solving parabolic problems on general grids shows CCCG({eta}) to be an efficient general purpose solver.

A variable-order time-dependent neutron transport method for nuclear reactor kinetics using analytically-integrated space-time characteristics

International Nuclear Information System (INIS)

Hoffman, A. J.; Lee, J. C.

2013-01-01

A new time-dependent neutron transport method based on the method of characteristics (MOC) has been developed. Whereas most spatial kinetics methods treat time dependence through temporal discretization, this new method treats time dependence by defining the characteristics to span space and time. In this implementation regions are defined in space-time where the thickness of the region in time fulfills an analogous role to the time step in discretized methods. The time dependence of the local source is approximated using a truncated Taylor series expansion with high order derivatives approximated using backward differences, permitting the solution of the resulting space-time characteristic equation. To avoid a drastic increase in computational expense and memory requirements due to solving many discrete characteristics in the space-time planes, the temporal variation of the boundary source is similarly approximated. This allows the characteristics in the space-time plane to be represented analytically rather than discretely, resulting in an algorithm comparable in implementation and expense to one that arises from conventional time integration techniques. Furthermore, by defining the boundary flux time derivative in terms of the preceding local source time derivative and boundary flux time derivative, the need to store angularly-dependent data is avoided without approximating the angular dependence of the angular flux time derivative. The accuracy of this method is assessed through implementation in the neutron transport code DeCART. The method is employed with variable-order local source representation to model a TWIGL transient. The results demonstrate that this method is accurate and more efficient than the discretized method. (authors)

Time-dependent reduced density matrix functional theory applied to laser-driven, correlated two-electron dynamics

Energy Technology Data Exchange (ETDEWEB)

Brics, Martins; Kapoor, Varun; Bauer, Dieter [Institut fuer Physik, Universitaet Rostock, 18051 Rostock (Germany)

2013-07-01

Time-dependent density functional theory (TDDFT) with known and practicable exchange-correlation potentials does not capture highly correlated electron dynamics such as single-photon double ionization, autoionization, or nonsequential ionization. Time-dependent reduced density matrix functional theory (TDRDMFT) may remedy these problems. The key ingredients in TDRDMFT are the natural orbitals (NOs), i.e., the eigenfunctions of the one-body reduced density matrix (1-RDM), and the occupation numbers (OCs), i.e., the respective eigenvalues. The two-body reduced density matrix (2-RDM) is then expanded in NOs, and equations of motion for the NOs can be derived. If the expansion coefficients of the 2-RDM were known exactly, the problem at hand would be solved. In practice, approximations have to be made. We study the prospects of TDRDMFT following a top-down approach. We solve the exact two-electron time-dependent Schroedinger equation for a model Helium atom in intense laser fields in order to study highly correlated phenomena such as the population of autoionizing states or single-photon double ionization. From the exact wave function we calculate the exact NOs, OCs, the exact expansion coefficients of the 2-RDM, and the exact potentials in the equations of motion. In that way we can identify how many NOs and which level of approximations are necessary to capture such phenomena.

Rotating Hele-Shaw cell with a time-dependent angular velocity

Science.gov (United States)

Anjos, Pedro H. A.; Alvarez, Victor M. M.; Dias, Eduardo O.; Miranda, José A.

2017-12-01

Despite the large number of existing studies of viscous flows in rotating Hele-Shaw cells, most investigations analyze rotational motion with a constant angular velocity, under vanishing Reynolds number conditions in which inertial effects can be neglected. In this work, we examine the linear and weakly nonlinear dynamics of the interface between two immiscible fluids in a rotating Hele-Shaw cell, considering the action of a time-dependent angular velocity, and taking into account the contribution of inertia. By using a generalized Darcy's law, we derive a second-order mode-coupling equation which describes the time evolution of the interfacial perturbation amplitudes. For arbitrary values of viscosity and density ratios, and for a range of values of a rotational Reynolds number, we investigate how the time-dependent angular velocity and inertia affect the important finger competition events that traditionally arise in rotating Hele-Shaw flows.

Time-dependent theory of double ionization of helium under XUV radiation

International Nuclear Information System (INIS)

Nikolopoulos, L A A; Lambropoulos, P

2007-01-01

We present non-perturbative time-dependent calculations of single and double ionization of helium, under XUV radiation of photon energy ranging from 40 to 45 eV, through the direct propagation of the time-dependent Schroedinger equation. The time-dependent wavefunction of the atom under the field is expanded in terms of correlated multichannel states normalized with incoming-wave boundary conditions. In addition to presenting a new non-perturbative approach to the three-body problem, in a fully correlated scheme, capable of providing in the same calculation photoelectron energy and angularly resolved spectra, as well as cross sections through the lowest non-vanishing order transition amplitude, we also present a detailed comparison of the values of certain key quantities that have been obtained through a variety of other methods. The degree of agreement we find, while lending credence to the approach and its versatility, also highlights the remaining open questions in this novel context of double ionization

A note on Burgers' equation with time delay: Instability via finite-time blow-up

International Nuclear Information System (INIS)

Jordan, P.M.

2008-01-01

Burgers' equation with time delay is considered. Using the Cole-Hopf transformation, the exact solution of this nonlinear partial differential equation (PDE) is determined in the context of a (seemingly) well-posed initial-boundary value problem (IBVP) involving homogeneous Dirichlet data. The solution obtained, however, is shown to exhibit a delay-induced instability, suffering blow-up in finite-time

Theory of coherent time-dependent transport in one-dimensional multiband semiconductor super-lattices

DEFF Research Database (Denmark)

Rotvig, J.; Smith, H.; Jauho, Antti-Pekka

1996-01-01

We present an analytical study of one-dimensional semiconductor superlattices in external electric fields, which may be time dependent. A number of general results for the (quasi)energies and eigenstates are derived. An equation of motion for the density matrix is obtained for a two-band model...

Travelling wave solutions for some time-delayed equations through factorizations

International Nuclear Information System (INIS)

Fahmy, E.S.

2008-01-01

In this work, we use factorization method to find explicit particular travelling wave solutions for the following important nonlinear second-order partial differential equations: The generalized time-delayed Burgers-Huxley, time-delayed convective Fishers, and the generalized time-delayed Burgers-Fisher. Using the particular solutions for these equations we find the general solutions, two-parameter solution, as special cases

Time-dependent simulation of organic light-emitting diodes

International Nuclear Information System (INIS)

Sharifi, M J

2009-01-01

Several methods to simulate the behavior of organic light-emitting diodes (OLEDs) have been proposed in the past. In this paper, we develop a previous method, based on the master equation, in order to allow the simulation of time-dependent behavior and transient states. The calculation algorithm of the program that we have written is described. The time-dependent behaviors of two simple monolayer devices and of a more complicated three-layer device were simulated by means of this program, and the results are discussed. The results show that the turn-off speed of an OLED might be very slow, especially in the case of a multilayer device. This behavior is related to the low mobility of the organic material in weak electric fields. An interesting feature of the time behavior is pointed out, whereby the recombination rate may become considerably larger after the falling edge of an applied voltage pulse. Moreover, the validity of the transient electro-luminescent method for measuring carrier mobility in organic material has been examined by means of simulation. The results show that there is some inconsistency especially in high electric fields

Space-Time Discrete KPZ Equation

Science.gov (United States)

Cannizzaro, G.; Matetski, K.

2018-03-01

We study a general family of space-time discretizations of the KPZ equation and show that they converge to its solution. The approach we follow makes use of basic elements of the theory of regularity structures (Hairer in Invent Math 198(2):269-504, 2014) as well as its discrete counterpart (Hairer and Matetski in Discretizations of rough stochastic PDEs, 2015. arXiv:1511.06937). Since the discretization is in both space and time and we allow non-standard discretization for the product, the methods mentioned above have to be suitably modified in order to accommodate the structure of the models under study.

Similarity Solutions for Multiterm Time-Fractional Diffusion Equation

Directory of Open Access Journals (Sweden)

A. Elsaid

2016-01-01

Full Text Available Similarity method is employed to solve multiterm time-fractional diffusion equation. The orders of the fractional derivatives belong to the interval (0,1] and are defined in the Caputo sense. We illustrate how the problem is reduced from a multiterm two-variable fractional partial differential equation to a multiterm ordinary fractional differential equation. Power series solution is obtained for the resulting ordinary problem and the convergence of the series solution is discussed. Based on the obtained results, we propose a definition for a multiterm error function with generalized coefficients.

Solution of Moving Boundary Space-Time Fractional Burger’s Equation

Directory of Open Access Journals (Sweden)

E. A.-B. Abdel-Salam

2014-01-01

Full Text Available The fractional Riccati expansion method is used to solve fractional differential equations with variable coefficients. To illustrate the effectiveness of the method, the moving boundary space-time fractional Burger’s equation is studied. The obtained solutions include generalized trigonometric and hyperbolic function solutions. Among these solutions, some are found for the first time. The linear and periodic moving boundaries for the kink solution of the Burger’s equation are presented graphically and discussed.

Remarks on time-dependent [current]-density functional theory for open quantum systems.

Science.gov (United States)

Yuen-Zhou, Joel; Aspuru-Guzik, Alán

2013-08-14

Time-dependent [current]-density functional theory for open quantum systems (OQS) has emerged as a formalism that can incorporate dissipative effects in the dynamics of many-body quantum systems. Here, we review and clarify some formal aspects of these theories that have been recently questioned in the literature. In particular, we provide theoretical support for the following conclusions: (1) contrary to what we and others had stated before, within the master equation framework, there is in fact a one-to-one mapping between vector potentials and current densities for fixed initial state, particle-particle interaction, and memory kernel; (2) regardless of the first conclusion, all of our recently suggested Kohn-Sham (KS) schemes to reproduce the current and particle densities of the original OQS, and in particular, the use of a KS closed driven system, remains formally valid; (3) the Lindblad master equation maintains the positivity of the density matrix regardless of the time-dependence of the Hamiltonian or the dissipation operators; (4) within the stochastic Schrödinger equation picture, a one-to-one mapping from stochastic vector potential to stochastic current density for individual trajectories has not been proven so far, except in the case where the vector potential is the same for every member of the ensemble, in which case, it reduces to the Lindblad master equation picture; (5) master equations may violate certain desired properties of the density matrix, such as positivity, but they remain as one of the most useful constructs to study OQS when the environment is not easily incorporated explicitly in the calculation. The conclusions support our previous work as formally rigorous, offer new insights into it, and provide a common ground to discuss related theories.

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.

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.

Large-Time Behavior of Solutions to Vlasov-Poisson-Fokker-Planck Equations: From Evanescent Collisions to Diffusive Limit

Science.gov (United States)

Herda, Maxime; Rodrigues, L. Miguel

2018-03-01

The present contribution investigates the dynamics generated by the two-dimensional Vlasov-Poisson-Fokker-Planck equation for charged particles in a steady inhomogeneous background of opposite charges. We provide global in time estimates that are uniform with respect to initial data taken in a bounded set of a weighted L^2 space, and where dependencies on the mean-free path τ and the Debye length δ are made explicit. In our analysis the mean free path covers the full range of possible values: from the regime of evanescent collisions τ → ∞ to the strongly collisional regime τ → 0. As a counterpart, the largeness of the Debye length, that enforces a weakly nonlinear regime, is used to close our nonlinear estimates. Accordingly we pay a special attention to relax as much as possible the τ -dependent constraint on δ ensuring exponential decay with explicit τ -dependent rates towards the stationary solution. In the strongly collisional limit τ → 0, we also examine all possible asymptotic regimes selected by a choice of observation time scale. Here also, our emphasis is on strong convergence, uniformity with respect to time and to initial data in bounded sets of a L^2 space. Our proofs rely on a detailed study of the nonlinear elliptic equation defining stationary solutions and a careful tracking and optimization of parameter dependencies of hypocoercive/hypoelliptic estimates.

Similarity solutions of time-dependent relativistic radiation-hydrodynamical plane-parallel flows

Science.gov (United States)

Fukue, Jun

2018-04-01

Similarity solutions are examined for the frequency-integrated relativistic radiation-hydrodynamical flows, which are described by the comoving quantities. The flows are vertical plane-parallel time-dependent ones with a gray opacity coefficient. For adequate boundary conditions, the flows are accelerated in a somewhat homologous manner, but terminate at some singular locus, which originates from the pathological behavior in relativistic radiation moment equations truncated in finite orders.

Finite element approximation for time-dependent diffusion with measure-valued source

Czech Academy of Sciences Publication Activity Database

Seidman, T.; Gobbert, M.; Trott, D.; Kružík, Martin

2012-01-01

Roč. 122, č. 4 (2012), s. 709-723 ISSN 0029-599X R&D Projects: GA AV ČR IAA100750802 Institutional support: RVO:67985556 Keywords : measure-valued source * diffusion equation Subject RIV: BA - General Mathematics Impact factor: 1.329, year: 2012 http://library.utia.cas.cz/separaty/2012/MTR/kruzik-finite element approximation for time - dependent diffusion with measure-valued source.pdf

Direct Yaw Control of Vehicle using State Dependent Riccati Equation with Integral Terms

Directory of Open Access Journals (Sweden)

SANDHU, F.

2016-05-01

Full Text Available Direct yaw control of four-wheel vehicles using optimal controllers such as the linear quadratic regulator (LQR and the sliding mode controller (SMC either considers only certain parameters constant in the nonlinear equations of vehicle model or totally neglect their effects to obtain simplified models, resulting in loss of states for the system. In this paper, a modified state-dependent Ricatti equation method obtained by the simplification of the vehicle model is proposed. This method overcomes the problem of the lost states by including state integrals. The results of the proposed system are compared with the sliding mode slip controller and state-dependent Ricatti equation method using high fidelity vehicle model in the vehicle simulation software package, Carsim. Results show 38% reduction in the lateral velocity, 34% reduction in roll and 16% reduction in excessive yaw by only increasing the fuel consumption by 6.07%.

Communication: On the calculation of time-dependent electron flux within the Born-Oppenheimer approximation: A flux-flux reflection principle

Science.gov (United States)

Albert, Julian; Hader, Kilian; Engel, Volker

2017-12-01

It is commonly assumed that the time-dependent electron flux calculated within the Born-Oppenheimer (BO) approximation vanishes. This is not necessarily true if the flux is directly determined from the continuity equation obeyed by the electron density. This finding is illustrated for a one-dimensional model of coupled electronic-nuclear dynamics. There, the BO flux is in perfect agreement with the one calculated from a solution of the time-dependent Schrödinger equation for the coupled motion. A reflection principle is derived where the nuclear BO flux is mapped onto the electronic flux.

Time-dependent anisotropic external sources in transient 3-D transport code TORT-TD

International Nuclear Information System (INIS)

Seubert, A.; Pautz, A.; Becker, M.; Dagan, R.

2009-01-01

This paper describes the implementation of a time-dependent distributed external source in TORT-TD by explicitly considering the external source in the ''fixed-source'' term of the implicitly time-discretised 3-D discrete ordinates transport equation. Anisotropy of the external source is represented by a spherical harmonics series expansion similar to the angular fluxes. The YALINA-Thermal subcritical assembly serves as a test case. The configuration with 280 fuel rods has been analysed with TORT-TD using cross sections in 18 energy groups and P1 scattering order generated by the KAPROS code system. Good agreement is achieved concerning the multiplication factor. The response of the system to an artificial time-dependent source consisting of two square-wave pulses demonstrates the time-dependent external source capability of TORT-TD. The result is physically plausible as judged from validation calculations. (orig.)

Two-dimensional time dependent calculations for the training reactor of Budapest University of Technology and Economics

International Nuclear Information System (INIS)

Mahmoud, K.S.; Szatmary, Z.

2005-01-01

An iterative method was developed for the numerical solution of the coupled two-dimensional time dependent multigroup diffusion equation and delayed precursor equations. Both forward (Explicit) and backward (Implicit) schemes were used. The second scheme was found to be numerically stable, while the first scheme requires that Δt -10 sec. for stability. An example is given for the second method. (authors)