Knapp, Jenny; Umhang, Gérald; Poulle, Marie-Lazarine; Millon, Laurence
2016-05-15
Studying the environmental occurrence of parasites of concern for humans and animals based on coprosamples is an expanding field of work in epidemiology and the ecology of health. Detecting and quantifying Toxocara spp. and Echinococcus multilocularis, two predominant zoonotic helminths circulating in European carnivores, in feces may help to better target measures for prevention. A rapid, sensitive, and one-step quantitative PCR (qPCR) allowing detection of E. multilocularis and Toxocara spp. was developed in the present study, combined with a host fecal test based on the identification of three carnivores (red fox, dog, and cat) involved in the life cycles of these parasites. A total of 68 coprosamples were collected from identified specimens from Vulpes vulpes, Canis lupus familiaris, Canis lupus, Felis silvestris catus, Meles meles, Martes foina, and Martes martes With DNA coprosamples, real-time PCR was performed in duplex with a qPCR inhibitor control specifically designed for this study. All the coprosample host identifications were confirmed by qPCR combined with sequencing, and parasites were detected and confirmed (E. multilocularis in red foxes and Toxocara cati in cats; 16% of samples presented inhibition). By combining parasite detection and quantification, the host fecal test, and a new qPCR inhibitor control, we created a technique with a high sensitivity that may considerably improve environmental studies of pathogens. PMID:26969697
Simulating system dynamics with arbitrary time step
NASA Astrophysics Data System (ADS)
Kantorovich, L.
2007-02-01
We suggest a dynamic simulation method that allows efficient and realistic modeling of kinetic processes, such as atomic diffusion, in which time has its actual meaning. Our method is similar in spirit to widely used kinetic Monte Carlo (KMC) techniques; however, in our approach, the time step can be chosen arbitrarily. This has an advantage in some cases, e.g., when the transition rates change with time sufficiently fast over the period of the KMC time step (e.g., due to time dependence of some external factors influencing kinetics such as moving scanning probe microscopy tip or external time-dependent field) or when the clock time is set by some external conditions, and it is convenient to use equal time steps instead of the random choice of the KMC algorithm in order to build up probability distribution functions. We show that an arbitrary choice of the time step can be afforded by building up the complete list of events including the “residence site” and multihop transitions. The idea of the method is illustrated in a simple “toy” model of a finite one-dimensional lattice of potential wells with unequal jump rates to either side, which can be studied analytically. We show that for any choice of the time step, our general kinetics method reproduces exactly the solution of the corresponding master equations for any choice of the time steps. The final kinetics also matches the standard KMC, and this allows better understanding of this algorithm, in which the time step is chosen in a certain way and the system always advances by a single hop.
Grief: Difficult Times, Simple Steps.
ERIC Educational Resources Information Center
Waszak, Emily Lane
This guide presents techniques to assist others in coping with the loss of a loved one. Using the language of 9 layperson, the book contains more than 100 tips for caregivers or loved ones. A simple step is presented on each page, followed by reasons and instructions for each step. Chapters include: "What to Say"; "Helpful Things to Do"; "Dealing…
Extrapolated implicit-explicit time stepping.
Constantinescu, E. M.; Sandu, A.; Mathematics and Computer Science; Virginia Polytechnic Inst. and State Univ.
2010-01-01
This paper constructs extrapolated implicit-explicit time stepping methods that allow one to efficiently solve problems with both stiff and nonstiff components. The proposed methods are based on Euler steps and can provide very high order discretizations of ODEs, index-1 DAEs, and PDEs in the method-of-lines framework. Implicit-explicit schemes based on extrapolation are simple to construct, easy to implement, and straightforward to parallelize. This work establishes the existence of perturbed asymptotic expansions of global errors, explains the convergence orders of these methods, and studies their linear stability properties. Numerical results with stiff ODE, DAE, and PDE test problems confirm the theoretical findings and illustrate the potential of these methods to solve multiphysics multiscale problems.
Obtaining Runge-Kutta Solutions Between Time Steps
NASA Technical Reports Server (NTRS)
Horn, M. K.
1984-01-01
New interpolation method used with existing Runge-Kutta algorithms. Algorithm evaluates solution at intermediate point within integration step. Only few additional computations required to produce intermediate solution data. Runge-Kutta method provides accurate solution with larger time steps than allowable in other methods.
40 CFR 264.113 - Closure; time allowed for closure.
Code of Federal Regulations, 2011 CFR
2011-07-01
... 40 Protection of Environment 26 2011-07-01 2011-07-01 false Closure; time allowed for closure. 264... (CONTINUED) STANDARDS FOR OWNERS AND OPERATORS OF HAZARDOUS WASTE TREATMENT, STORAGE, AND DISPOSAL FACILITIES Closure and Post-Closure § 264.113 Closure; time allowed for closure. (a) Within 90 days after...
40 CFR 265.113 - Closure; time allowed for closure.
Code of Federal Regulations, 2011 CFR
2011-07-01
... 40 Protection of Environment 26 2011-07-01 2011-07-01 false Closure; time allowed for closure. 265... (CONTINUED) INTERIM STATUS STANDARDS FOR OWNERS AND OPERATORS OF HAZARDOUS WASTE TREATMENT, STORAGE, AND DISPOSAL FACILITIES Closure and Post-Closure § 265.113 Closure; time allowed for closure. (a) Within...
Progress Report on Alloy 617 Time Dependent Allowables
Wright, Julie Knibloe
2015-06-01
Time dependent allowable stresses are required in the ASME Boiler and Pressure Vessel Code for design of components in the temperature range where time dependent deformation (i.e., creep) is expected to become significant. There are time dependent allowable stresses in Section IID of the Code for use in the non-nuclear construction codes, however, there are additional criteria that must be considered in developing time dependent allowables for nuclear components. These criteria are specified in Section III NH. St is defined as the lesser of three quantities: 100% of the average stress required to obtain a total (elastic, plastic, primary and secondary creep) strain of 1%; 67% of the minimum stress to cause rupture; and 80% of the minimum stress to cause the initiation of tertiary creep. The values are reported for a range of temperatures and for time increments up to 100,000 hours. These values are determined from uniaxial creep tests, which involve the elevated temperature application of a constant load which is relatively small, resulting in deformation over a long time period prior to rupture. The stress which is the minimum resulting from these criteria is the time dependent allowable stress St. In this report data from a large number of creep and creep-rupture tests on Alloy 617 are analyzed using the ASME Section III NH criteria. Data which are used in the analysis are from the ongoing DOE sponsored high temperature materials program, form Korea Atomic Energy Institute through the Generation IV VHTR Materials Program and historical data from previous HTR research and vendor data generated in developing the alloy. It is found that the tertiary creep criterion determines St at highest temperatures, while the stress to cause 1% total strain controls at low temperatures. The ASME Section III Working Group on Allowable Stress Criteria has recommended that the uncertainties associated with determining the onset of tertiary creep and the lack of significant
Simulating stochastic dynamics using large time steps.
Corradini, O; Faccioli, P; Orland, H
2009-12-01
We present an approach to investigate the long-time stochastic dynamics of multidimensional classical systems, in contact with a heat bath. When the potential energy landscape is rugged, the kinetics displays a decoupling of short- and long-time scales and both molecular dynamics or Monte Carlo (MC) simulations are generally inefficient. Using a field theoretic approach, we perform analytically the average over the short-time stochastic fluctuations. This way, we obtain an effective theory, which generates the same long-time dynamics of the original theory, but has a lower time-resolution power. Such an approach is used to develop an improved version of the MC algorithm, which is particularly suitable to investigate the dynamics of rare conformational transitions. In the specific case of molecular systems at room temperature, we show that elementary integration time steps used to simulate the effective theory can be chosen a factor approximately 100 larger than those used in the original theory. Our results are illustrated and tested on a simple system, characterized by a rugged energy landscape. PMID:20365123
Margul, Daniel T; Tuckerman, Mark E
2016-05-10
the polarizable AMOEBA water model. As was seen with fixed-charge models, we are able to obtain large time steps exceeding 100 fs, allowing calculations to be performed 10 to 20 times faster than standard thermostated molecular dynamics. PMID:27054809
Time to pause before the next step
Siemon, R.E.
1998-12-31
Many scientists, who have staunchly supported ITER for years, are coming to realize it is time to further rethink fusion energy`s development strategy. Specifically, as was suggested by Grant Logan and Dale Meade, and in keeping with the restructuring of 1996, a theme of better, cheaper, faster fusion would serve the program more effectively than ``demonstrating controlled ignition...and integrated testing of the high-heat-flux and nuclear components required to utilize fusion energy...`` which are the important ingredients of ITER`s objectives. The author has personally shifted his view for a mixture of technical and political reasons. On the technical side, he senses that through advanced tokamak research, spherical tokamak research, and advanced stellarator work, scientists are coming to a new understanding that might make a burning-plasma device significantly smaller and less expensive. Thus waiting for a few years, even ten years, seems prudent. Scientifically, there is fascinating physics to be learned through studies of burning plasma on a tokamak. And clearly if one wishes to study burning plasma physics in a sustained plasma, there is no other configuration with an adequate database on which to proceed. But what is the urgency of moving towards an ITER-like step focused on burning plasma? Some of the arguments put forward and the counter arguments are discussed here.
Adaptive time steps in trajectory surface hopping simulations.
Spörkel, Lasse; Thiel, Walter
2016-05-21
Trajectory surface hopping (TSH) simulations are often performed in combination with active-space multi-reference configuration interaction (MRCI) treatments. Technical problems may arise in such simulations if active and inactive orbitals strongly mix and switch in some particular regions. We propose to use adaptive time steps when such regions are encountered in TSH simulations. For this purpose, we present a computational protocol that is easy to implement and increases the computational effort only in the critical regions. We test this procedure through TSH simulations of a GFP chromophore model (OHBI) and a light-driven rotary molecular motor (F-NAIBP) on semiempirical MRCI potential energy surfaces, by comparing the results from simulations with adaptive time steps to analogous ones with constant time steps. For both test molecules, the number of successful trajectories without technical failures rises significantly, from 53% to 95% for OHBI and from 25% to 96% for F-NAIBP. The computed excited-state lifetime remains essentially the same for OHBI and increases somewhat for F-NAIBP, and there is almost no change in the computed quantum efficiency for internal rotation in F-NAIBP. We recommend the general use of adaptive time steps in TSH simulations with active-space CI methods because this will help to avoid technical problems, increase the overall efficiency and robustness of the simulations, and allow for a more complete sampling. PMID:27208937
Adaptive time steps in trajectory surface hopping simulations
NASA Astrophysics Data System (ADS)
Spörkel, Lasse; Thiel, Walter
2016-05-01
Trajectory surface hopping (TSH) simulations are often performed in combination with active-space multi-reference configuration interaction (MRCI) treatments. Technical problems may arise in such simulations if active and inactive orbitals strongly mix and switch in some particular regions. We propose to use adaptive time steps when such regions are encountered in TSH simulations. For this purpose, we present a computational protocol that is easy to implement and increases the computational effort only in the critical regions. We test this procedure through TSH simulations of a GFP chromophore model (OHBI) and a light-driven rotary molecular motor (F-NAIBP) on semiempirical MRCI potential energy surfaces, by comparing the results from simulations with adaptive time steps to analogous ones with constant time steps. For both test molecules, the number of successful trajectories without technical failures rises significantly, from 53% to 95% for OHBI and from 25% to 96% for F-NAIBP. The computed excited-state lifetime remains essentially the same for OHBI and increases somewhat for F-NAIBP, and there is almost no change in the computed quantum efficiency for internal rotation in F-NAIBP. We recommend the general use of adaptive time steps in TSH simulations with active-space CI methods because this will help to avoid technical problems, increase the overall efficiency and robustness of the simulations, and allow for a more complete sampling.
Multiple time step integrators in ab initio molecular dynamics
Luehr, Nathan; Martínez, Todd J.; Markland, Thomas E.
2014-02-28
Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we present two schemes that enable efficient time-scale separation in ab initio calculations: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy.
Dependence of aqua-planet simulations on time step
NASA Astrophysics Data System (ADS)
Williamson, David L.; Olson, Jerry G.
2003-04-01
Aqua-planet simulations with Eulerian and semi-Lagrangian dynamical cores coupled to the NCAR CCM3 parametrization suite produce very different zonal average precipitation patterns. The model with the Eulerian core forms a narrow single precipitation peak centred on the sea surface temperature (SST) maximum. The one with the semi-Lagrangian core forms a broad structure often with a double peak straddling the SST maximum with a precipitation minimum centred on the SST maximum. The different structure is shown to be caused primarily by the different time step adopted by each core and its effect on the parametrizations rather than by different truncation errors introduced by the dynamical cores themselves. With a longer discrete time step, the surface exchange parametrization deposits more moisture in the atmosphere in a single time step, resulting in convection being initiated farther from the equator, closer to the maximum source. Different diffusive smoothing associated with different spectral resolutions is a secondary effect influencing the strength of the double structure. When the semi-Lagrangian core is configured to match the Eulerian with the same time step, a three-time-level formulation and same spectral truncation it produces precipitation fields similar to those from the Eulerian. It is argued that the broad and double structure forms in this model with the longer time step because more water is put into the atmosphere over a longer discrete time step, the evaporation rate being the same. The additional water vapour in the region of equatorial moisture convergence results in more convective available potential energy farther from the equator which allows convection to initiate farther from the equator.The resulting heating drives upward vertical motion and low-level convergence away from the equator, resulting in much weaker upward motion at the equator. The feedback between the convective heating and dynamics reduces the instability at the equator and
Seven Steps to On-Time Delivery.
ERIC Educational Resources Information Center
Konchar, Mark; Sanvido, Victor
1999-01-01
Describes seven steps to consider when making project-delivery decisions that include defining the school district's goals and profile, selecting the project-delivery system and procurement method, selecting the project team and contract type, and developing and confirming the facility program. Concluding comments address the district review of…
Assessment Data at Your Fingertips: Advances Allow for Timely Reporting
ERIC Educational Resources Information Center
Bolch, Matt
2009-01-01
The ever-increasing standards of No Child Left Behind regulations and various state assessments have put more pressure on teachers and administrators to monitor the learning process. Fortunately, the advent of technology is allowing teachers to test more often to prepare students for high-stakes tests and for districts to understand results for…
Collocation and Galerkin Time-Stepping Methods
NASA Technical Reports Server (NTRS)
Huynh, H. T.
2011-01-01
We study the numerical solutions of ordinary differential equations by one-step methods where the solution at tn is known and that at t(sub n+1) is to be calculated. The approaches employed are collocation, continuous Galerkin (CG) and discontinuous Galerkin (DG). Relations among these three approaches are established. A quadrature formula using s evaluation points is employed for the Galerkin formulations. We show that with such a quadrature, the CG method is identical to the collocation method using quadrature points as collocation points. Furthermore, if the quadrature formula is the right Radau one (including t(sub n+1)), then the DG and CG methods also become identical, and they reduce to the Radau IIA collocation method. In addition, we present a generalization of DG that yields a method identical to CG and collocation with arbitrary collocation points. Thus, the collocation, CG, and generalized DG methods are equivalent, and the latter two methods can be formulated using the differential instead of integral equation. Finally, all schemes discussed can be cast as s-stage implicit Runge-Kutta methods.
Time scaling relations for step bunches from models with step-step attractions (B1-type models)
NASA Astrophysics Data System (ADS)
Krasteva, A.; Popova, H.; Akutsu, N.; Tonchev, V.
2016-03-01
The step bunching instability is studied in three models of step motion defined in terms of ordinary differential equations (ODE). The source of instability in these models is step-step attraction, it is opposed by step-step repulsion and the developing surface patterns reflect the balance between the two. The first model, TE2, is a generalization of the seminal model of Tersoff et al. (1995). The second one, LW2, is obtained from the model of Liu and Weeks (1998) using the repulsions term to construct the attractions one with retained possibility to change the parameters in the two independently. The third model, MM2, is a minimal one constructed ad hoc and in this article it plays a central role. New scheme for scaling the ODE in vicinal studies is applied towards deciphering the pre-factors in the time-scaling relations. In all these models the patterned surface is self-similar - only one length scale is necessary to describe its evolution (hence B1-type). The bunches form finite angles with the terraces. Integrating numerically the equations for step motion and changing systematically the parameters we obtain the overall dependence of time-scaling exponent β on the power of step-step attractions p as β = 1/(3+p) for MM2 and hypothesize based on restricted set of data that it is β = 1/(5+p) for LW2 and TE2.
Improving quality: one step at a time.
Blankson-seck, N; Butta, P
1999-01-01
The notion that health care workers have the power to improve the quality of their services is a key to AVSC's efforts worldwide. The COPE process, AVSC's low-cost intervention for improving quality at service sites, brings together supervisors and staff at all levels to identify barriers to quality services and helps them find solutions they can implement with their own resources. For example, a hospital in Tanzania had tried unsuccessfully to obtain the funds to repair or replace broken equipment. Using the COPE process, the hospital used available funds to send a technician for training in maintenance and repair. Now everything from blood pressure equipment to bedsprings is repaired promptly, and quality has improved. Another hospital in Tanzania coped with the problem of broken bedsprings (patients were putting mattresses on the floor) by using readily available wire mesh to make repairs. In Kenya, the lack of running water forced staff to collect water from a cistern, taking time from their other responsibilities. During a COPE meeting to resolve the problem the staff bemoaned the fact that they did not have the funds to replace the water system. Then the gardener told the group that all they needed to do was fix a broken pipe. The repair was made at minimal cost, and the water supply was restored. The COPE process reveals that health care staff not only can identify obstacles to quality, they often know the cause of the problem and can offer the best solutions. PMID:12295155
NASA Astrophysics Data System (ADS)
Ficchi, Andrea; Perrin, Charles; Andréassian, Vazken
2015-04-01
We investigate the operational utility of fine time step hydro-climatic information using a large catchment data set. The originality of this data set lies in the availability of precipitation data from the 6-minute rain gauges of Météo-France, and in the size of the catchment set (217 French catchments in total). The rainfall-runoff model used (GR4) has been adapted to hourly and sub-hourly time steps (up to 6-minute) from the daily time step version (Perrin et al., 2003). The model is applied at different time steps ranging from 6-minute to 1 day (6-, 12-, 30-minute, 1-, 3-, 6-, 12-hour and 1 day) and the evolution of model performance for each catchment is evaluated at the daily time step by aggregation of model outputs. Three classes of behavior are found according to the trend of model performance as the time step becomes finer: (i) catchments presenting an improvement of model performance; (ii) catchments with a model performance insensitive to the time step; (iii) catchments for which the performance even deteriorates as the time step becomes finer. The reasons behind these different trends are investigated from a hydrological point of view, by relating the model sensitivity to data at finer time step to catchment descriptors. References: Perrin, C., C. Michel and V. Andréassian (2003), "Improvement of a parsimonious model for streamflow simulation", Journal of Hydrology, 279(1-4): 275-289.
Multiple-time-stepping generalized hybrid Monte Carlo methods
NASA Astrophysics Data System (ADS)
Escribano, Bruno; Akhmatskaya, Elena; Reich, Sebastian; Azpiroz, Jon M.
2015-01-01
Performance of the generalized shadow hybrid Monte Carlo (GSHMC) method [1], which proved to be superior in sampling efficiency over its predecessors [2-4], molecular dynamics and hybrid Monte Carlo, can be further improved by combining it with multi-time-stepping (MTS) and mollification of slow forces. We demonstrate that the comparatively simple modifications of the method not only lead to better performance of GSHMC itself but also allow for beating the best performed methods, which use the similar force splitting schemes. In addition we show that the same ideas can be successfully applied to the conventional generalized hybrid Monte Carlo method (GHMC). The resulting methods, MTS-GHMC and MTS-GSHMC, provide accurate reproduction of thermodynamic and dynamical properties, exact temperature control during simulation and computational robustness and efficiency. MTS-GHMC uses a generalized momentum update to achieve weak stochastic stabilization to the molecular dynamics (MD) integrator. MTS-GSHMC adds the use of a shadow (modified) Hamiltonian to filter the MD trajectories in the HMC scheme. We introduce a new shadow Hamiltonian formulation adapted to force-splitting methods. The use of such Hamiltonians improves the acceptance rate of trajectories and has a strong impact on the sampling efficiency of the method. Both methods were implemented in the open-source MD package ProtoMol and were tested on a water and a protein systems. Results were compared to those obtained using a Langevin Molly (LM) method [5] on the same systems. The test results demonstrate the superiority of the new methods over LM in terms of stability, accuracy and sampling efficiency. This suggests that putting the MTS approach in the framework of hybrid Monte Carlo and using the natural stochasticity offered by the generalized hybrid Monte Carlo lead to improving stability of MTS and allow for achieving larger step sizes in the simulation of complex systems.
Multiple-time-stepping generalized hybrid Monte Carlo methods
Escribano, Bruno; Akhmatskaya, Elena; Reich, Sebastian; Azpiroz, Jon M.
2015-01-01
Performance of the generalized shadow hybrid Monte Carlo (GSHMC) method [1], which proved to be superior in sampling efficiency over its predecessors [2–4], molecular dynamics and hybrid Monte Carlo, can be further improved by combining it with multi-time-stepping (MTS) and mollification of slow forces. We demonstrate that the comparatively simple modifications of the method not only lead to better performance of GSHMC itself but also allow for beating the best performed methods, which use the similar force splitting schemes. In addition we show that the same ideas can be successfully applied to the conventional generalized hybrid Monte Carlo method (GHMC). The resulting methods, MTS-GHMC and MTS-GSHMC, provide accurate reproduction of thermodynamic and dynamical properties, exact temperature control during simulation and computational robustness and efficiency. MTS-GHMC uses a generalized momentum update to achieve weak stochastic stabilization to the molecular dynamics (MD) integrator. MTS-GSHMC adds the use of a shadow (modified) Hamiltonian to filter the MD trajectories in the HMC scheme. We introduce a new shadow Hamiltonian formulation adapted to force-splitting methods. The use of such Hamiltonians improves the acceptance rate of trajectories and has a strong impact on the sampling efficiency of the method. Both methods were implemented in the open-source MD package ProtoMol and were tested on a water and a protein systems. Results were compared to those obtained using a Langevin Molly (LM) method [5] on the same systems. The test results demonstrate the superiority of the new methods over LM in terms of stability, accuracy and sampling efficiency. This suggests that putting the MTS approach in the framework of hybrid Monte Carlo and using the natural stochasticity offered by the generalized hybrid Monte Carlo lead to improving stability of MTS and allow for achieving larger step sizes in the simulation of complex systems.
Time step and shadow Hamiltonian in molecular dynamics simulations
NASA Astrophysics Data System (ADS)
Kim, Sangrak
2015-08-01
We examine the time step and the shadow Hamiltonian of symplectic algorithms for a bound system of a simple harmonic oscillator as a specific example. The phase space trajectory moves on the hyperplane of a constant shadow Hamiltonian. We find a stationary condition for the time step τ n with which the motion repeats itself on the phase space with a period n. Interestingly, that the time steps satisfying the stationary condition turn out to be independent of the symplectic algorithms chosen. Furthermore, the phase volume enclosed by the phase trajectory is given by n τ n Ẽ n , where Ẽ n is the initial shadow energy of the corresponding symplectic algorithm.
Improving Leadership and Management Practices: One Step at a Time
ERIC Educational Resources Information Center
Bella, Jill
2008-01-01
Taking small steps toward change is a sensible way to improve the leadership and management practices in an early care and education program. A director must be able to make continuous improvements without alienating staff by asking them to make drastic changes that seem overwhelming and unachievable. Taking on change one step at a time is a way…
Accuracy-based time step criteria for solving parabolic equations
Mohtar, R.; Segerlind, L.
1995-12-31
Parabolic equations govern many transient engineering problems. Space integration using finite element or finite difference methods changes the parabolic partial differential equation into an ordinary differential equation. Time integration schemes are needed to solve the later equation. In order to accurately perform the later integration a proper time step must be provided. Time step estimates based on a stability criteria have been prescribed in the literature. The following paper presents time step estimates that satisfy stability as well as accuracy criteria. These estimates were correlated to the Froude and Courant Numbers. The later criteria were found to be overly conservative for some integration schemes. Suggestions as to which time integration scheme is the best to use are also presented.
IMPROVEMENTS TO THE TIME STEPPING ALGORITHM OF RELAP5-3D
Cumberland, R.; Mesina, G.
2009-01-01
The RELAP5-3D time step method is used to perform thermo-hydraulic and neutronic simulations of nuclear reactors and other devices. It discretizes time and space by numerically solving several differential equations. Previously, time step size was controlled by halving or doubling the size of a previous time step. This process caused the code to run slower than it potentially could. In this research project, the RELAP5-3D time step method was modifi ed to allow a new method of changing time steps to improve execution speed and to control error. The new RELAP5-3D time step method being studied involves making the time step proportional to the material courant limit (MCL), while insuring that the time step does not increase by more than a factor of two between advancements. As before, if a step fails or mass error is excessive, the time step is cut in half. To examine performance of the new method, a measure of run time and a measure of error were plotted against a changing MCL proportionality constant (m) in seven test cases. The removal of the upper time step limit produced a small increase in error, but a large decrease in execution time. The best value of m was found to be 0.9. The new algorithm is capable of producing a signifi cant increase in execution speed, with a relatively small increase in mass error. The improvements made are now under consideration for inclusion as a special option in the RELAP5-3D production code.
40 CFR 97.41 - Timing requirements for NOX allowance allocations.
Code of Federal Regulations, 2010 CFR
2010-07-01
... 40 Protection of Environment 20 2010-07-01 2010-07-01 false Timing requirements for NOX allowance... PROGRAMS (CONTINUED) FEDERAL NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS NOX Allowance Allocations § 97.41 Timing requirements for NOX allowance allocations. (a) The NOX allowance...
40 CFR 97.41 - Timing requirements for NOX allowance allocations.
Code of Federal Regulations, 2011 CFR
2011-07-01
... 40 Protection of Environment 21 2011-07-01 2011-07-01 false Timing requirements for NOX allowance... PROGRAMS (CONTINUED) FEDERAL NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS NOX Allowance Allocations § 97.41 Timing requirements for NOX allowance allocations. (a) The NOX allowance...
Sensitivity of a thermodynamic sea ice model with leads to time step size
NASA Technical Reports Server (NTRS)
Ledley, T. S.
1985-01-01
The characteristics of sea ice models, developed to study the physics of the growth and melt of ice at the ocean surface and the variations in ice extent, depend on the size of the time step. Thus, to study longer-term variations within a reasonable computer budget, a model with a scheme allowing longer time steps has been constructed. However, the results produced by the model can definitely depend on the length of the time step. The sensitivity of a model to time-step size can be reduced by appropriate approaches. The present investigation is concerned with experiments which use a formulation of a lead parameterization that can be considered as a first step toward the development of a lead parameterization suitable for a use in long-term climate studies.
40 CFR 96.41 - Timing requirements for NOX allowance allocations.
Code of Federal Regulations, 2011 CFR
2011-07-01
... 40 Protection of Environment 21 2011-07-01 2011-07-01 false Timing requirements for NOX allowance... PROGRAMS (CONTINUED) NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS FOR STATE IMPLEMENTATION PLANS NOX Allowance Allocations § 96.41 Timing requirements for NOX allowance allocations. (a)...
40 CFR 96.41 - Timing requirements for NOX allowance allocations.
Code of Federal Regulations, 2010 CFR
2010-07-01
... 40 Protection of Environment 20 2010-07-01 2010-07-01 false Timing requirements for NOX allowance... PROGRAMS (CONTINUED) NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS FOR STATE IMPLEMENTATION PLANS NOX Allowance Allocations § 96.41 Timing requirements for NOX allowance allocations. (a)...
Short-term Time Step Convergence in a Climate Model
Wan, Hui; Rasch, Philip J.; Taylor, Mark; Jablonowski, Christiane
2015-02-11
A testing procedure is designed to assess the convergence property of a global climate model with respect to time step size, based on evaluation of the root-mean-square temperature difference at the end of very short (1 h) simulations with time step sizes ranging from 1 s to 1800 s. A set of validation tests conducted without sub-grid scale parameterizations confirmed that the method was able to correctly assess the convergence rate of the dynamical core under various configurations. The testing procedure was then applied to the full model, and revealed a slow convergence of order 0.4 in contrast to the expected first-order convergence. Sensitivity experiments showed without ambiguity that the time stepping errors in the model were dominated by those from the stratiform cloud parameterizations, in particular the cloud microphysics. This provides a clear guidance for future work on the design of more accurate numerical methods for time stepping and process coupling in the model.
Short-term Time Step Convergence in a Climate Model
Wan, Hui; Rasch, Philip J.; Taylor, Mark; Jablonowski, Christiane
2015-02-11
A testing procedure is designed to assess the convergence property of a global climate model with respect to time step size, based on evaluation of the root-mean-square temperature difference at the end of very short (1 h) simulations with time step sizes ranging from 1 s to 1800 s. A set of validation tests conducted without sub-grid scale parameterizations confirmed that the method was able to correctly assess the convergence rate of the dynamical core under various configurations. The testing procedure was then applied to the full model, and revealed a slow convergence of order 0.4 in contrast to themore » expected first-order convergence. Sensitivity experiments showed without ambiguity that the time stepping errors in the model were dominated by those from the stratiform cloud parameterizations, in particular the cloud microphysics. This provides a clear guidance for future work on the design of more accurate numerical methods for time stepping and process coupling in the model.« less
Accuracy of Pedometer Steps and Time for Youth with Disabilities
ERIC Educational Resources Information Center
Beets, Michael W.; Combs, Cindy; Pitetti, Kenneth H.; Morgan, Melinda; Bryan, Rebecca R.; Foley, John T.
2007-01-01
The purpose of the study was to examine the accuracy of pedometer steps and activity time (Walk4Life, WL) for youth with developmental disabilities. Eighteen youth (11 girls, 7 boys) 4-14 years completed six 80-meter self-paced walking trials while wearing a pedometer at five waist locations (front right, front left, back right, back left, middle…
Modified Chebyshev pseudospectral method with O(N exp -1) time step restriction
NASA Technical Reports Server (NTRS)
Kosloff, Dan; Tal-Ezer, Hillel
1989-01-01
The extreme eigenvalues of the Chebyshev pseudospectral differentiation operator are O(N exp 2) where N is the number of grid points. As a result of this, the allowable time step in an explicit time marching algorithm is O(N exp -2) which, in many cases, is much below the time step dictated by the physics of the partial differential equation. A new set of interpolating points is introduced such that the eigenvalues of the differentiation operator are O(N) and the allowable time step is O(N exp -1). The properties of the new algorithm are similar to those of the Fourier method. The new algorithm also provides a highly accurate solution for non-periodic boundary value problems.
40 CFR 96.341 - Timing requirements for CAIR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2012 CFR
2012-07-01
... Season allowance allocations. 96.341 Section 96.341 Protection of Environment ENVIRONMENTAL PROTECTION... PROGRAMS FOR STATE IMPLEMENTATION PLANS CAIR NOX Ozone Season Allowance Allocations § 96.341 Timing requirements for CAIR NOX Ozone Season allowance allocations. (a) By October 31, 2006, the permitting...
40 CFR 97.341 - Timing requirements for CAIR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2013 CFR
2013-07-01
... Season allowance allocations. 97.341 Section 97.341 Protection of Environment ENVIRONMENTAL PROTECTION... TRADING PROGRAMS CAIR NOX Ozone Season Allowance Allocations § 97.341 Timing requirements for CAIR NOX Ozone Season allowance allocations. (a) The Administrator will determine by order the CAIR NOX...
40 CFR 97.341 - Timing requirements for CAIR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2012 CFR
2012-07-01
... Season allowance allocations. 97.341 Section 97.341 Protection of Environment ENVIRONMENTAL PROTECTION... TRADING PROGRAMS CAIR NOX Ozone Season Allowance Allocations § 97.341 Timing requirements for CAIR NOX Ozone Season allowance allocations. (a) The Administrator will determine by order the CAIR NOX...
40 CFR 96.341 - Timing requirements for CAIR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2013 CFR
2013-07-01
... Season allowance allocations. 96.341 Section 96.341 Protection of Environment ENVIRONMENTAL PROTECTION... PROGRAMS FOR STATE IMPLEMENTATION PLANS CAIR NOX Ozone Season Allowance Allocations § 96.341 Timing requirements for CAIR NOX Ozone Season allowance allocations. (a) By October 31, 2006, the permitting...
40 CFR 97.341 - Timing requirements for CAIR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2014 CFR
2014-07-01
... Season allowance allocations. 97.341 Section 97.341 Protection of Environment ENVIRONMENTAL PROTECTION... TRADING PROGRAMS CAIR NOX Ozone Season Allowance Allocations § 97.341 Timing requirements for CAIR NOX Ozone Season allowance allocations. (a) The Administrator will determine by order the CAIR NOX...
40 CFR 96.341 - Timing requirements for CAIR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2014 CFR
2014-07-01
... Season allowance allocations. 96.341 Section 96.341 Protection of Environment ENVIRONMENTAL PROTECTION... PROGRAMS FOR STATE IMPLEMENTATION PLANS CAIR NOX Ozone Season Allowance Allocations § 96.341 Timing requirements for CAIR NOX Ozone Season allowance allocations. (a) By October 31, 2006, the permitting...
40 CFR 97.411 - Timing requirements for TR NOX Annual allowance allocations.
Code of Federal Regulations, 2013 CFR
2013-07-01
... 40 Protection of Environment 22 2013-07-01 2013-07-01 false Timing requirements for TR NOX Annual... TR NOX Annual Trading Program § 97.411 Timing requirements for TR NOX Annual allowance allocations. (a) Existing units. (1) TR NOX Annual allowances are allocated, for the control periods in 2012...
40 CFR 97.711 - Timing requirements for TR SO2 Group 2 allowance allocations.
Code of Federal Regulations, 2012 CFR
2012-07-01
... 40 Protection of Environment 22 2012-07-01 2012-07-01 false Timing requirements for TR SO2 Group 2... TR SO2 Group 2 Trading Program § 97.711 Timing requirements for TR SO2 Group 2 allowance allocations. (a) Existing units. (1) TR SO2 Group 2 allowances are allocated, for the control periods in 2012...
40 CFR 97.711 - Timing requirements for TR SO2 Group 2 allowance allocations.
Code of Federal Regulations, 2014 CFR
2014-07-01
... 40 Protection of Environment 21 2014-07-01 2014-07-01 false Timing requirements for TR SO2 Group 2... TR SO2 Group 2 Trading Program § 97.711 Timing requirements for TR SO2 Group 2 allowance allocations. (a) Existing units. (1) TR SO2 Group 2 allowances are allocated, for the control periods in 2012...
40 CFR 97.611 - Timing requirements for TR SO2 Group 1 allowance allocations.
Code of Federal Regulations, 2013 CFR
2013-07-01
... 40 Protection of Environment 22 2013-07-01 2013-07-01 false Timing requirements for TR SO2 Group 1... TR SO2 Group 1 Trading Program § 97.611 Timing requirements for TR SO2 Group 1 allowance allocations. (a) Existing units. (1) TR SO2 Group 1 allowances are allocated, for the control periods in 2012...
40 CFR 97.511 - Timing requirements for TR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2013 CFR
2013-07-01
... 40 Protection of Environment 22 2013-07-01 2013-07-01 false Timing requirements for TR NOX Ozone... TRADING PROGRAMS TR NOX Ozone Season Trading Program § 97.511 Timing requirements for TR NOX Ozone Season allowance allocations. (a) Existing units. (1) TR NOX Ozone Season allowances are allocated, for the...
40 CFR 97.511 - Timing requirements for TR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2012 CFR
2012-07-01
... 40 Protection of Environment 22 2012-07-01 2012-07-01 false Timing requirements for TR NOX Ozone... TRADING PROGRAMS TR NOX Ozone Season Trading Program § 97.511 Timing requirements for TR NOX Ozone Season allowance allocations. (a) Existing units. (1) TR NOX Ozone Season allowances are allocated, for the...
40 CFR 97.611 - Timing requirements for TR SO2 Group 1 allowance allocations.
Code of Federal Regulations, 2014 CFR
2014-07-01
... 40 Protection of Environment 21 2014-07-01 2014-07-01 false Timing requirements for TR SO2 Group 1... TR SO2 Group 1 Trading Program § 97.611 Timing requirements for TR SO2 Group 1 allowance allocations. (a) Existing units. (1) TR SO2 Group 1 allowances are allocated, for the control periods in 2012...
40 CFR 97.411 - Timing requirements for TR NOX Annual allowance allocations.
Code of Federal Regulations, 2014 CFR
2014-07-01
... 40 Protection of Environment 21 2014-07-01 2014-07-01 false Timing requirements for TR NOX Annual... TR NOX Annual Trading Program § 97.411 Timing requirements for TR NOX Annual allowance allocations. (a) Existing units. (1) TR NOX Annual allowances are allocated, for the control periods in 2012...
40 CFR 97.411 - Timing requirements for TR NOX Annual allowance allocations.
Code of Federal Regulations, 2012 CFR
2012-07-01
... 40 Protection of Environment 22 2012-07-01 2012-07-01 false Timing requirements for TR NOX Annual... TR NOX Annual Trading Program § 97.411 Timing requirements for TR NOX Annual allowance allocations. (a) Existing units. (1) TR NOX Annual allowances are allocated, for the control periods in 2012...
40 CFR 97.511 - Timing requirements for TR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2014 CFR
2014-07-01
... 40 Protection of Environment 21 2014-07-01 2014-07-01 false Timing requirements for TR NOX Ozone... TRADING PROGRAMS TR NOX Ozone Season Trading Program § 97.511 Timing requirements for TR NOX Ozone Season allowance allocations. (a) Existing units. (1) TR NOX Ozone Season allowances are allocated, for the...
40 CFR 97.711 - Timing requirements for TR SO2 Group 2 allowance allocations.
Code of Federal Regulations, 2013 CFR
2013-07-01
... 40 Protection of Environment 22 2013-07-01 2013-07-01 false Timing requirements for TR SO2 Group 2... TR SO2 Group 2 Trading Program § 97.711 Timing requirements for TR SO2 Group 2 allowance allocations. (a) Existing units. (1) TR SO2 Group 2 allowances are allocated, for the control periods in 2012...
40 CFR 97.611 - Timing requirements for TR SO2 Group 1 allowance allocations.
Code of Federal Regulations, 2012 CFR
2012-07-01
... 40 Protection of Environment 22 2012-07-01 2012-07-01 false Timing requirements for TR SO2 Group 1... TR SO2 Group 1 Trading Program § 97.611 Timing requirements for TR SO2 Group 1 allowance allocations. (a) Existing units. (1) TR SO2 Group 1 allowances are allocated, for the control periods in 2012...
40 CFR 97.141 - Timing requirements for CAIR NOX allowance allocations.
Code of Federal Regulations, 2010 CFR
2010-07-01
... 40 Protection of Environment 20 2010-07-01 2010-07-01 false Timing requirements for CAIR NOX... (CONTINUED) AIR PROGRAMS (CONTINUED) FEDERAL NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS CAIR NOX Allowance Allocations § 97.141 Timing requirements for CAIR NOX allowance allocations. (a)...
40 CFR 97.141 - Timing requirements for CAIR NOX allowance allocations.
Code of Federal Regulations, 2011 CFR
2011-07-01
... 40 Protection of Environment 21 2011-07-01 2011-07-01 false Timing requirements for CAIR NOX... (CONTINUED) AIR PROGRAMS (CONTINUED) FEDERAL NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS CAIR NOX Allowance Allocations § 97.141 Timing requirements for CAIR NOX allowance allocations. (a)...
A method for improving time-stepping numerics
NASA Astrophysics Data System (ADS)
Williams, P. D.
2012-04-01
In contemporary numerical simulations of the atmosphere, evidence suggests that time-stepping errors may be a significant component of total model error, on both weather and climate time-scales. This presentation will review the available evidence, and will then suggest a simple but effective method for substantially improving the time-stepping numerics at no extra computational expense. The most common time-stepping method is the leapfrog scheme combined with the Robert-Asselin (RA) filter. This method is used in the following atmospheric models (and many more): ECHAM, MAECHAM, MM5, CAM, MESO-NH, HIRLAM, KMCM, LIMA, SPEEDY, IGCM, PUMA, COSMO, FSU-GSM, FSU-NRSM, NCEP-GFS, NCEP-RSM, NSEAM, NOGAPS, RAMS, and CCSR/NIES-AGCM. Although the RA filter controls the time-splitting instability in these models, it also introduces non-physical damping and reduces the accuracy. This presentation proposes a simple modification to the RA filter. The modification has become known as the RAW filter (Williams 2011). When used in conjunction with the leapfrog scheme, the RAW filter eliminates the non-physical damping and increases the amplitude accuracy by two orders, yielding third-order accuracy. (The phase accuracy remains second-order.) The RAW filter can easily be incorporated into existing models, typically via the insertion of just a single line of code. Better simulations are obtained at no extra computational expense. Results will be shown from recent implementations of the RAW filter in various atmospheric models, including SPEEDY and COSMO. For example, in SPEEDY, the skill of weather forecasts is found to be significantly improved. In particular, in tropical surface pressure predictions, five-day forecasts made using the RAW filter have approximately the same skill as four-day forecasts made using the RA filter (Amezcua, Kalnay & Williams 2011). These improvements are encouraging for the use of the RAW filter in other models.
Accurate Monotonicity - Preserving Schemes With Runge-Kutta Time Stepping
NASA Technical Reports Server (NTRS)
Suresh, A.; Huynh, H. T.
1997-01-01
A new class of high-order monotonicity-preserving schemes for the numerical solution of conservation laws is presented. The interface value in these schemes is obtained by limiting a higher-order polynominal reconstruction. The limiting is designed to preserve accuracy near extrema and to work well with Runge-Kutta time stepping. Computational efficiency is enhanced by a simple test that determines whether the limiting procedure is needed. For linear advection in one dimension, these schemes are shown as well as the Euler equations also confirm their high accuracy, good shock resolution, and computational efficiency.
30 CFR 203.66 - What happens if MMS does not act in the time allowed?
Code of Federal Regulations, 2011 CFR
2011-07-01
... Expansion Projects § 203.66 What happens if MMS does not act in the time allowed? If we do not act within.... (b) An expansion project You get a royalty suspension for the first year of production Abide...
Enabling fast, stable and accurate peridynamic computations using multi-time-step integration
Lindsay, P.; Parks, M. L.; Prakash, A.
2016-04-13
Peridynamics is a nonlocal extension of classical continuum mechanics that is well-suited for solving problems with discontinuities such as cracks. This paper extends the peridynamic formulation to decompose a problem domain into a number of smaller overlapping subdomains and to enable the use of different time steps in different subdomains. This approach allows regions of interest to be isolated and solved at a small time step for increased accuracy while the rest of the problem domain can be solved at a larger time step for greater computational efficiency. Lastly, performance of the proposed method in terms of stability, accuracy, andmore » computational cost is examined and several numerical examples are presented to corroborate the findings.« less
40 CFR 96.141 - Timing requirements for CAIR NOX allowance allocations.
Code of Federal Regulations, 2010 CFR
2010-07-01
... 40 Protection of Environment 20 2010-07-01 2010-07-01 false Timing requirements for CAIR NOX... (CONTINUED) AIR PROGRAMS (CONTINUED) NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS FOR STATE IMPLEMENTATION PLANS CAIR NOX Allowance Allocations § 96.141 Timing requirements for CAIR...
40 CFR 96.141 - Timing requirements for CAIR NOX allowance allocations.
Code of Federal Regulations, 2011 CFR
2011-07-01
... 40 Protection of Environment 21 2011-07-01 2011-07-01 false Timing requirements for CAIR NOX... (CONTINUED) AIR PROGRAMS (CONTINUED) NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS FOR STATE IMPLEMENTATION PLANS CAIR NOX Allowance Allocations § 96.141 Timing requirements for CAIR...
40 CFR 97.341 - Timing requirements for CAIR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2010 CFR
2010-07-01
... 40 Protection of Environment 20 2010-07-01 2010-07-01 false Timing requirements for CAIR NOX Ozone... AGENCY (CONTINUED) AIR PROGRAMS (CONTINUED) FEDERAL NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS CAIR NOX Ozone Season Allowance Allocations § 97.341 Timing requirements for CAIR...
40 CFR 97.341 - Timing requirements for CAIR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2011 CFR
2011-07-01
... 40 Protection of Environment 21 2011-07-01 2011-07-01 false Timing requirements for CAIR NOX Ozone... AGENCY (CONTINUED) AIR PROGRAMS (CONTINUED) FEDERAL NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS CAIR NOX Ozone Season Allowance Allocations § 97.341 Timing requirements for CAIR...
Viral DNA Packaging: One Step at a Time
NASA Astrophysics Data System (ADS)
Bustamante, Carlos; Moffitt, Jeffrey R.
During its life-cycle the bacteriophage φ29 actively packages its dsDNA genome into a proteinacious capsid, compressing its genome to near crystalline densities against large electrostatic, elastic, and entropic forces. This remarkable process is accomplished by a nano-scale, molecular DNA pump - a complex assembly of three protein and nucleic acid rings which utilizes the free energy released in ATP hydrolysis to perform the mechanical work necessary to overcome these large energetic barriers. We have developed a single molecule optical tweezers assay which has allowed us to probe the detailed mechanism of this packaging motor. By following the rate of packaging of a single bacteriophage as the capsid is filled with genome and as a function of optically applied load, we find that the compression of the genome results in the build-up of an internal force, on the order of ˜ 55 pN, due to the compressed genome. The ability to work against such large forces makes the packaging motor one of the strongest known molecular motors. By titrating the concentration of ATP, ADP, and inorganic phosphate at different opposing load, we are able to determine features of the mechanochemistry of this motor - the coupling between the mechanical and chemical cycles. We find that force is generated not upon binding of ATP, but rather upon release of hydrolysis products. Finally, by improving the resolution of the optical tweezers assay, we are able to observe the discrete increments of DNA encapsidated each cycle of the packaging motor. We find that DNA is packaged in 10-bp increments preceded by the binding of multiple ATPs. The application of large external forces slows the packaging rate of the motor, revealing that the 10-bp steps are actually composed of four 2.5-bp steps which occur in rapid succession. These data show that the individual subunits of the pentameric ring-ATPase at the core of the packaging motor are highly coordinated, with the binding of ATP and the
Accelerating spectral-element simulations of seismic wave propagation using local time stepping
NASA Astrophysics Data System (ADS)
Peter, D. B.; Rietmann, M.; Galvez, P.; Nissen-Meyer, T.; Grote, M.; Schenk, O.
2013-12-01
Seismic tomography using full-waveform inversion requires accurate simulations of seismic wave propagation in complex 3D media. However, finite element meshing in complex media often leads to areas of local refinement, generating small elements that accurately capture e.g. strong topography and/or low-velocity sediment basins. For explicit time schemes, this dramatically reduces the global time-step for wave-propagation problems due to numerical stability conditions, ultimately making seismic inversions prohibitively expensive. To alleviate this problem, local time stepping (LTS) algorithms allow an explicit time-stepping scheme to adapt the time-step to the element size, allowing near-optimal time-steps everywhere in the mesh. Numerical simulations are thus liberated of global time-step constraints potentially speeding up simulation runtimes significantly. We present here a new, efficient multi-level LTS-Newmark scheme for general use with spectral-element methods (SEM) with applications in seismic wave propagation. We fit the implementation of our scheme onto the package SPECFEM3D_Cartesian, which is a widely used community code, simulating seismic and acoustic wave propagation in earth-science applications. Our new LTS scheme extends the 2nd-order accurate Newmark time-stepping scheme, and leads to an efficient implementation, producing real-world speedup of multi-resolution seismic applications. Furthermore, we generalize the method to utilize many refinement levels with a design specifically for continuous finite elements. We demonstrate performance speedup using a state-of-the-art dynamic earthquake rupture model for the Tohoku-Oki event, which is currently limited by small elements along the rupture fault. Utilizing our new algorithmic LTS implementation together with advances in exploiting graphic processing units (GPUs), numerical seismic wave propagation simulations in complex media will dramatically reduce computation times, empowering high
40 CFR 96.141 - Timing requirements for CAIR NOX allowance allocations.
Code of Federal Regulations, 2013 CFR
2013-07-01
... 40 Protection of Environment 22 2013-07-01 2013-07-01 false Timing requirements for CAIR NOX... allowance allocations. (a) By October 31, 2006, the permitting authority will submit to the Administrator... § 96.142(a) and (b), for the control periods in 2009, 2010, 2011, 2012, 2013, and 2014. (b) By...
40 CFR 96.141 - Timing requirements for CAIR NOX allowance allocations.
Code of Federal Regulations, 2012 CFR
2012-07-01
... 40 Protection of Environment 22 2012-07-01 2012-07-01 false Timing requirements for CAIR NOX... allowance allocations. (a) By October 31, 2006, the permitting authority will submit to the Administrator... § 96.142(a) and (b), for the control periods in 2009, 2010, 2011, 2012, 2013, and 2014. (b) By...
40 CFR 60.4141 - Timing requirements for Hg allowance allocations.
Code of Federal Regulations, 2010 CFR
2010-07-01
... 40 Protection of Environment 6 2010-07-01 2010-07-01 false Timing requirements for Hg allowance... accordance with § 60.4142(a) and (b), for the control periods in 2010, 2011, 2012, 2013, and 2014. (b)(1) By October 31, 2008 and October 31 of each year thereafter, the permitting authority will submit to...
40 CFR 60.4141 - Timing requirements for Hg allowance allocations.
Code of Federal Regulations, 2011 CFR
2011-07-01
... 40 Protection of Environment 6 2011-07-01 2011-07-01 false Timing requirements for Hg allowance... accordance with § 60.4142(a) and (b), for the control periods in 2010, 2011, 2012, 2013, and 2014. (b)(1) By October 31, 2008 and October 31 of each year thereafter, the permitting authority will submit to...
40 CFR 96.341 - Timing requirements for CAIR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2011 CFR
2011-07-01
... 40 Protection of Environment 21 2011-07-01 2011-07-01 false Timing requirements for CAIR NOX Ozone... AGENCY (CONTINUED) AIR PROGRAMS (CONTINUED) NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS FOR STATE IMPLEMENTATION PLANS CAIR NOX Ozone Season Allowance Allocations § 96.341...
40 CFR 96.341 - Timing requirements for CAIR NOX Ozone Season allowance allocations.
Code of Federal Regulations, 2010 CFR
2010-07-01
... 40 Protection of Environment 20 2010-07-01 2010-07-01 false Timing requirements for CAIR NOX Ozone... AGENCY (CONTINUED) AIR PROGRAMS (CONTINUED) NOX BUDGET TRADING PROGRAM AND CAIR NOX AND SO2 TRADING PROGRAMS FOR STATE IMPLEMENTATION PLANS CAIR NOX Ozone Season Allowance Allocations § 96.341...
Rules Allowing Extended Time on Graduation: Advocates Debate Effects of Change in Regulations
ERIC Educational Resources Information Center
Gewertz, Catherine
2009-01-01
Federal regulations have opened a door that allows schools to get credit under the No Child Left Behind Act for students who take longer than four years to earn a high school diploma. That option worries some education advocates, who fear it could relieve valuable pressure on high schools to graduate students on time. Under the law's…
NASA Astrophysics Data System (ADS)
Karimi, S.; Nakshatrala, K. B.
2014-12-01
Advection-Diffusion-Reaction (ADR) equations play a crucial role in simulating numerous geo- physical phenomena. It is well-known that the solution to these equations exhibit disparate spatial and temporal scales. These mathematical scales occur due to relative dominance of either advec- tion, diffusion, or reaction processes. Hence, in a careful simulation, one has to choose appropriate time-integrators, time-steps, and numerical formulations for spatial discretization. Multi-time-step coupling methods allow specific choice of integration methods (either temporal or spatial) in dif- ferent regions of the spatial domain. In recent years, most of the attempts to design monolithic multi-time-step frameworks favored second-order transient systems in structural dynamics. In this presentation, we will introduce monolithic multi-time-step computational frameworks for ADR equations. These methods are based on the theory of differential/algebraic equations. We shall also provide an overview of results from stability analysis, study of drift from compatibility con- straints, and analysis of influence of perturbations. Several benchmark problems will be utilized to demonstrate the theoretical findings and features of the proposed frameworks. Finally, application of the proposed methods to fast bimolecular reactive systems will be shown.
A new time-stepping method for regional climate models
NASA Astrophysics Data System (ADS)
Williams, P. D.
2010-12-01
The dynamical cores of many regional climate models use the Robert-Asselin filter to suppress the spurious computational mode of the leapfrog scheme. Unfortunately, whilst successfully eliminating the unwanted mode, the Robert-Asselin filter also weakly suppresses the physical solution and degrades the numerical accuracy. These two concomitant problems occur because the filter does not conserve the mean state, averaged over the three time slices on which it operates. This presentation proposes a simple modification to the Robert-Asselin filter, which does conserve the three-time-level mean state. When used in conjunction with the leapfrog scheme, the modification vastly reduces the artificial damping of the physical solution. Correspondingly, the modification increases the numerical accuracy for amplitude errors by two orders, yielding third-order accuracy. The modified filter may easily be incorporated into existing regional climate models, via the addition of only a few lines of code that are computationally very inexpensive. Results will be shown from recent implementations of the modified filter in various models. The modification will be shown to reduce model biases and to significantly improve the predictive skill. Magnitude of the complex amplification factor as a function of the non-dimensional time step, for leapfrog integrations. This quantity would be identical to 1 for a perfect numerical scheme. Clearly, the filter proposed here (case α=0.53) has much smaller numerical errors than the original Robert-Asselin filter (case α=1). Moreover, the proposed filter is trivial to implement and is no more computationally expensive. Taken from Williams (2009; Monthly Weather Review).
Watching Proteins Direct Crystal Growth One Step at a Time
2009-01-01
Researchers at Berkeley Labs Molecular Foundry use an atomic force microscope to record this movie of a peptide being adsorbed to a crystal surface while two successive crystal steps interact, then progress beyond the peptide. The peptide temporarily slows the step before transferring up to the next atomic layer. The lattice pattern on the surface corresponds to the molecular structure of the underlying crystal.
The USMLE Step 2 CS: Time for a change.
Alvin, Matthew D
2016-08-01
The United States Medical Licensing Examination (USMLE(®)) Steps are a series of mandatory licensing assessments for all allopathic (MD degree) medical students in their transition from student to intern to resident physician. Steps 1, 2 Clinical Knowledge (CK), and 3 are daylong multiple-choice exams that quantify a medical student's basic science and clinical knowledge as well as their application of that knowledge using a three-digit score. In doing so, these Steps provide a standardized assessment that residency programs use to differentiate applicants and evaluate their competitiveness. Step 2 Clinical Skills (CS), the only other Step exam and the second component of Step 2, was created in 2004 to test clinical reasoning and patient-centered skills. As a Pass/Fail exam without a numerical scoring component, Step 2 CS provides minimal differentiation among applicants for residency programs. In this personal view article, it is argued that the current Step 2 CS exam should be eliminated for US medical students and propose an alternative consistent with the mission and purpose of the exam that imposes less of a burden on medical students. PMID:27007882
Constrained Density Functional Theory by Imaginary Time-Step Method
NASA Astrophysics Data System (ADS)
Kidd, Daniel
Constrained Density Functional Theory (CDFT) has been a popular choice within the last decade for sidestepping the self interaction problem within long-range charge transfer calculations. Typically an inner constraint loop is added within the self-consistent field iterations of DFT in order to enforce this charge transfer state by means of a Lagrange multiplier method. In this work, an alternate implementation of CDFT is introduced, that of the imaginary time-step method, which lends itself more readily to real space calculations in the ability to solve numerically for 3D local external potentials which enforce arbitrary given densities. This method has been shown to reproduce the proper 1 / R dependence of charge transfer systems in real space calculations as well as the ability to generate useful constraint potentials. As an example application, this method is shown to be capable of describing defects within periodic systems using finite calculations by constraining the 3D density to that of the periodically calculated perfect system at the boundaries.
DNA walks one step at a time in electrophoresis
NASA Astrophysics Data System (ADS)
Guan, Juan; Wang, Bo; Granick, Steve
2011-03-01
Testing the classical view that in DNA gel electrophoresis, long polymer chains navigate through their gel environment via reptation, we reach a different conclusion: this driven motion proceeds by stick-slip. Our single-molecule experiments visualize fluorescent-labeled lambda-DNA, whose intramolecular conformations are resolved with 30 ms resolution using home-written software. Combining hundreds to thousands of trajectories under amplitudes of electric field ranging from zero to large, we quantify the full statistical distribution of motion with unprecedented statistics. Pauses are seen between steps of driven motion, probably reflecting that the chain is trapped inside the gel matrix. The pausing time is exponentially distributed and decreases with increasing electric field strength, suggesting that the jerky behavior is an activated process, facilitated by electric field. We propose a stretch-assisted mechanism: that the energy barrier to move through the gel environment is first overcome by a leading segment, the ensuing intramolecular stress from stretching causing lagging segments to recoil and follow along.
Calvo-Lorenzo, M S; Hulbert, L E; Fowler, A L; Louie, A; Gershwin, L J; Pinkerton, K E; Ballou, M A; Klasing, K C; Mitloehner, F M
2016-06-01
Dairy calves in the western United States are commonly raised individually in wooden hutches with a space allowance of 1.23m(2)/calf. Recent legislative initiatives in California and across the United States were passed regarding concern over space allowance for farm animals. The objective of this study was to determine if rearing male Holstein calves in wooden hutches modified to increase space allowance would influence measures of performance, lying time per day, health, and respiratory immunocompetence. At 4d of age, 60 calves were randomly assigned to 1 of 3housing treatments: (1) conventional housing (CONV; 1.23m(2)/calf), (2) 1.5 × CONV (MOD; 1.85m(2)/calf), or (3) 3 × CONV (MAX; 3.71m(2)/calf). Intakes of milk and solid feed were recorded daily and body weight was measured at 0, 3, 6, 10, and 12 wk of age. For the first 3 wk of the trial, calves were scored daily for fecal consistency, hydration, and hide cleanliness. In addition, calves were scored for respiratory health (i.e., nasal and eye discharge, ear position) until 7 wk of age. The total lying duration per day was recorded using data loggers at 3, 6, and 10 wk of age. Eight clinically healthy calves from each treatment were sensitized with subcutaneous ovalbumin (OVA) and then challenged with aerosolized OVA to assess calf respiratory immunity at 11 wk of age. Bronchoalveolar lavage fluid (BALF) was collected 4d after the OVA challenge and analyzed for leukocyte differentials and OVA-specific IgG, IgG1, IgA, and IgE. Calf average daily gain and body weight were positively associated with space allowance at approximately 3 wk before weaning and throughout postweaning, respectively. A greater space allowance decreased lying time after 46d. Space allowance did not influence fecal consistency, but there was a tendency for MAX calves to take 1d longer to recover from loose feces than MOD calves. The MAX calves had the fewest (%) observations with feces on their body compared with CONV or MOD. At 3 wk of
Two-stepping through time: mammals and viruses.
Meyerson, Nicholas R; Sawyer, Sara L
2011-06-01
Recent studies have identified ancient virus genomes preserved as fossils within diverse animal genomes. These fossils have led to the revelation that a broad range of mammalian virus families are older and more ubiquitous than previously appreciated. Long-term interactions between viruses and their hosts often develop into genetic arms races where both parties continually jockey for evolutionary dominance. It is difficult to imagine how mammalian hosts have kept pace in the evolutionary race against rapidly evolving viruses over large expanses of time, given their much slower evolutionary rates. However, recent data has begun to reveal the evolutionary strategy of slowly-evolving hosts. We review these data and suggest a modified arms race model where the evolutionary possibilities of viruses are relatively constrained. Such a model could allow more accurate forecasting of virus evolution. PMID:21531564
NASA Astrophysics Data System (ADS)
Ramadan, Omar
2014-12-01
Systematic split-step finite difference time domain (SS-FDTD) formulations, based on the general Lie-Trotter-Suzuki product formula, are presented for solving the time-dependent Maxwell equations in double-dispersive electromagnetic materials. The proposed formulations provide a unified tool for constructing a family of unconditionally stable algorithms such as the first order split-step FDTD (SS1-FDTD), the second order split-step FDTD (SS2-FDTD), and the second order alternating direction implicit FDTD (ADI-FDTD) schemes. The theoretical stability of the formulations is included and it has been demonstrated that the formulations are unconditionally stable by construction. Furthermore, the dispersion relation of the formulations is derived and it has been found that the proposed formulations are best suited for those applications where a high space resolution is needed. Two-dimensional (2-D) and 3-D numerical examples are included and it has been observed that the SS1-FDTD scheme is computationally more efficient than the ADI-FDTD counterpart, while maintaining approximately the same numerical accuracy. Moreover, the SS2-FDTD scheme allows using larger time step than the SS1-FDTD or ADI-FDTD and therefore necessitates less CPU time, while giving approximately the same numerical accuracy.
Ruprecht, Verena; Wieser, Stefan; Marguet, Didier; Schütz, Gerhard J.
2011-01-01
Resolving the dynamical interplay of proteins and lipids in the live-cell plasma membrane represents a central goal in current cell biology. Superresolution concepts have introduced a means of capturing spatial heterogeneity at a nanoscopic length scale. Similar concepts for detecting dynamical transitions (superresolution chronoscopy) are still lacking. Here, we show that recently introduced spot-variation fluorescence correlation spectroscopy allows for sensing transient confinement times of membrane constituents at dramatically improved resolution. Using standard diffraction-limited optics, spot-variation fluorescence correlation spectroscopy captures signatures of single retardation events far below the transit time of the tracer through the focal spot. We provide an analytical description of special cases of transient binding of a tracer to pointlike traps, or association of a tracer with nanodomains. The influence of trap mobility and the underlying binding kinetics are quantified. Experimental approaches are suggested that allow for gaining quantitative mechanistic insights into the interaction processes of membrane constituents. PMID:21641330
Error correction in short time steps during the application of quantum gates
NASA Astrophysics Data System (ADS)
de Castro, L. A.; Napolitano, R. d. J.
2016-04-01
We propose a modification of the standard quantum error-correction method to enable the correction of errors that occur due to the interaction with a noisy environment during quantum gates without modifying the codification used for memory qubits. Using a perturbation treatment of the noise that allows us to separate it from the ideal evolution of the quantum gate, we demonstrate that in certain cases it is necessary to divide the logical operation in short time steps intercalated by correction procedures. A prescription of how these gates can be constructed is provided, as well as a proof that, even for the cases when the division of the quantum gate in short time steps is not necessary, this method may be advantageous for reducing the total duration of the computation.
The Semi-implicit Time-stepping Algorithm in MH4D
NASA Astrophysics Data System (ADS)
Vadlamani, Srinath; Shumlak, Uri; Marklin, George; Meier, Eric; Lionello, Roberto
2006-10-01
The Plasma Science and Innovation Center (PSI Center) at the University of Washington is developing MHD codes to accurately model Emerging Concept (EC) devices. Examination of the semi-implicit time stepping algorithm implemented in the tetrahedral mesh MHD simulation code, MH4D, is presented. The time steps for standard explicit methods, which are constrained by the Courant-Friedrichs-Lewy (CFL) condition, are typically small for simulations of EC experiments due to the large Alfven speed. The CFL constraint is more severe with a tetrahedral mesh because of the irregular cell geometry. The semi-implicit algorithm [1] removes the fast waves constraint, thus allowing for larger time steps. We will present the implementation method of this algorithm, and numerical results for test problems in simple geometry. Also, we will present the effectiveness in simulations of complex geometry, similar to the ZaP [2] experiment at the University of Washington. References: [1]Douglas S. Harned and D. D. Schnack, Semi-implicit method for long time scale magnetohy drodynamic computations in three dimensions, JCP, Volume 65, Issue 1, July 1986, Pages 57-70. [2]U. Shumlak, B. A. Nelson, R. P. Golingo, S. L. Jackson, E. A. Crawford, and D. J. Den Hartog, Sheared flow stabilization experiments in the ZaP flow Zpinch, Phys. Plasmas 10, 1683 (2003).
Automatic multirate methods for ordinary differential equations. [Adaptive time steps
Gear, C.W.
1980-01-01
A study is made of the application of integration methods in which different step sizes are used for different members of a system of equations. Such methods can result in savings if the cost of derivative evaluation is high or if a system is sparse; however, the estimation and control of errors is very difficult and can lead to high overheads. Three approaches are discussed, and it is shown that the least intuitive is the most promising. 2 figures.
Risk-based evaluation of Allowed Outage Times (AOTs) considering risk of shutdown
Mankamo, T.; Kim, I.S.; Samanta, P.K.
1992-12-31
When safety systems fail during power operation, Technical Specifications (TS) usually limit the repair within Allowed Outage Time (AOT). If the repair cannot be completed within the AOT, or no AOT is allowed, the plant is required to be shut down for the repair. However, if the capability to remove decay heat is degraded, shutting down the plant with the need to operate the affected decay-heat removal systems may impose a substantial risk compared to continued power operation over a usual repair time. Thus, defining a proper AOT in such situations can be considered as a risk-comparison between the repair in frill power state with a temporarily increased level of risk, and the altemative of shutting down the plant for the repair in zero power state with a specific associated risk. The methodology of the risk-comparison approach, with a due consideration of the shutdown risk, has been further developed and applied to the AOT considerations of residual heat removal and standby service water systems of a boiling water reactor (BWR) plant. Based on the completed work, several improvements to the TS requirements for the systems studied can be suggested.
Risk-based evaluation of Allowed Outage Times (AOTs) considering risk of shutdown
Mankamo, T. ); Kim, I.S.; Samanta, P.K. )
1992-01-01
When safety systems fail during power operation, Technical Specifications (TS) usually limit the repair within Allowed Outage Time (AOT). If the repair cannot be completed within the AOT, or no AOT is allowed, the plant is required to be shut down for the repair. However, if the capability to remove decay heat is degraded, shutting down the plant with the need to operate the affected decay-heat removal systems may impose a substantial risk compared to continued power operation over a usual repair time. Thus, defining a proper AOT in such situations can be considered as a risk-comparison between the repair in frill power state with a temporarily increased level of risk, and the altemative of shutting down the plant for the repair in zero power state with a specific associated risk. The methodology of the risk-comparison approach, with a due consideration of the shutdown risk, has been further developed and applied to the AOT considerations of residual heat removal and standby service water systems of a boiling water reactor (BWR) plant. Based on the completed work, several improvements to the TS requirements for the systems studied can be suggested.
McCrorie, P Rw; Duncan, E; Granat, M H; Stansfield, B W
2012-11-01
Evidence suggests that behaviours such as standing are beneficial for our health. Unfortunately, little is known of the prevalence of this state, its importance in relation to time spent stepping or variation across seasons. The aim of this study was to quantify, in young adolescents, the prevalence and seasonal changes in time spent upright and not stepping (UNSt(time)) as well as time spent upright and stepping (USt(time)), and their contribution to overall upright time (U(time)). Thirty-three adolescents (12.2 ± 0.3 y) wore the activPAL activity monitor during four school days on two occasions: November/December (winter) and May/June (summer). UNSt(time) contributed 60% of daily U(time) at winter (Mean = 196 min) and 53% at summer (Mean = 171 min); a significant seasonal effect, p < 0.001. USt(time) was significantly greater in summer compared to winter (153 min versus 131 min, p < 0.001). The effects in UNSt(time) could be explained through significant seasonal differences during the school hours (09:00-16:00), whereas the effects in USt(time) could be explained through significant seasonal differences in the evening period (16:00-22:00). Adolescents spent a greater amount of time upright and not stepping than they did stepping, in both winter and summer. The observed seasonal effects for both UNSt(time) and USt(time) provide important information for behaviour change intervention programs. PMID:23111187
One step at a time: endoplasmic reticulum-associated degradation
Vembar, Shruthi S.; Brodsky, Jeffrey L.
2009-01-01
Protein folding in the endoplasmic reticulum (ER) is monitored by ER quality control (ERQC) mechanisms. Proteins that pass ERQC criteria traffic to their final destinations through the secretory pathway, whereas non-native and unassembled subunits of multimeric proteins are degraded by the ER-associated degradation (ERAD) pathway. During ERAD, molecular chaperones and associated factors recognize and target substrates for retrotranslocation to the cytoplasm, where they are degraded by the ubiquitin–proteasome machinery. The discovery of diseases that are associated with ERAD substrates highlights the importance of this pathway. Here, we summarize our current understanding of each step during ERAD, with emphasis on the factors that catalyse distinct activities. PMID:19002207
Kim, Taewon; Rhee, Joohyun; Wright, David L
2016-01-01
Two experiments were conducted to examine the efficacy of random (RP) and blocked practice (BP) for enhancing later motor learning. Each experiment involved practicing three unique seven key serial reaction time (SRT) tasks in either a blocked or random format followed by practice of a novel SRT task either 2-min (Experiment 1) or 24-h (Experiment 2) later. While the expected benefit of RP for retention was present in both experiments, in Experiment 1 there was no advantage from prior RP for new learning. Experiment 2 explored the possibility that increasing the interval, from 2-min to 24-h, between BP or RP and practice of the novel motor task might allow consolidation of sequence knowledge acquired during BP or RP which in turn might facilitate new learning. As a result of the additional time between training bouts RP facilitated the rate at which the novel motor task was acquired. Interestingly, when this additional time was provided, both BP and RP supported (a) a performance saving for the first trial with the novel task, and (b) an offline improvement in performance across a 24-h interval not present when only the novel motor task was practiced. The latter benefits for new learning may have resulted from exposure to prior physical practice per se. or practice variability. These data are discussed with respect to (a) future learning benefits from prior experience training with greater CI, and (b) the importance of memory consolidation for motor learning. PMID:26686835
Evaluation of allowed outage times (AOTs) from a risk and reliability standpoint
Vesely, W.E. )
1989-08-01
This report describes the basic risks which are associated with allowed outage times (AOTs), defines strategies for selecting the risks to be quantified, and describes how the risks can be quantified. The report furthermore describes criteria considerations in determining the acceptability of calculated AOT risks, and discusses the merits of relative risk criteria versus absolute risk criteria. The detailed evaluations which are involved in calculating AOT risks, including uncertainty considerations are also discussed. The report also describes the proper ways that risks from multiple AOTs should be considered so that risks are properly accumulated from proposed multiple AOT changes, but are not double-counted. Generally, average AOT risks which include the frequency of occurrence of the AOT need to be accumulated but single downtime risks don't since they apply to individual AOTs. 8 refs., 22 tabs.
PIC Algorithm with Multiple Poisson Equation Solves During One Time Step
NASA Astrophysics Data System (ADS)
Ren, Junxue; Godar, Trenton; Menart, James; Mahalingam, Sudhakar; Choi, Yongjun; Loverich, John; Stoltz, Peter H.
2015-09-01
In order to reduce the overall computational time of a PIC (particle-in-cell) computer simulation, an attempt was made to utilize larger time step sizes by implementing multiple solutions of Poisson's equation within one time step. The hope was this would make the PIC simulation stable at larger time steps than an explicit technique can use, and using larger time steps would reduce the overall computational time, even though the computational time per time step would increase. A three-dimensional PIC code that tracks electrons and ions throughout a three-dimensional Cartesian computational domain is used to perform this study. The results of altering the number of times Poisson's equation is solved during a single time step are presented. Also, the size of the time that can be used and still maintain a stable solution is surveyed. The results indicate that using multiple Poisson solves during one time step provides some ability to use larger time steps in PIC simulations, but the increase in time step size is not significant and the overall simulation run time is not reduced
Space-time variability of floods across Germany: Gradual trends, step changes and fluctuations
NASA Astrophysics Data System (ADS)
Merz, Bruno; Vorogushyn, Sergiy; Viet Dung, Nguyen; Schröter, Kai
2015-04-01
The space-time variability of flood magnitude and frequency across Germany at the interannual and decadal time scale is analyzed and interpreted. The analyses are based on flood time series of 68 catchments for a joint period of 74 years. The catchments are distributed across Germany and show different flood regimes. Different statistical tests are applied to investigate different types of flood changes: gradual trends, step changes and fluctuations. In addition, changes in the mean behavior and in the variability are studied. A focus is placed on the spatial stability of changes, i.e. answering the question to which extent flood changes are coherent across Germany. The joint analysis of changes for a large number of catchments allows interpreting the causes of the observed changes. For instance, climate-related flood changes are expected to show a different behavior than changes caused by river training or land-use change.
Empirical versus time stepping with embedded error control for density-driven flow in porous media
NASA Astrophysics Data System (ADS)
Younes, Anis; Ackerer, Philippe
2010-08-01
Modeling density-driven flow in porous media may require very long computational time due to the nonlinear coupling between flow and transport equations. Time stepping schemes are often used to adapt the time step size in order to reduce the computational cost of the simulation. In this work, the empirical time stepping scheme which adapts the time step size according to the performance of the iterative nonlinear solver is compared to an adaptive time stepping scheme where the time step length is controlled by the temporal truncation error. Results of the simulations of the Elder problem show that (1) the empirical time stepping scheme can lead to inaccurate results even with a small convergence criterion, (2) accurate results are obtained when the time step size selection is based on the truncation error control, (3) a non iterative scheme with proper time step management can be faster and leads to more accurate solution than the standard iterative procedure with the empirical time stepping and (4) the temporal truncation error can have a significant effect on the results and can be considered as one of the reasons for the differences observed in the Elder numerical results.
Competencies for Part-Time Faculty--the First Step.
ERIC Educational Resources Information Center
Haddad, Margaret; Dickens, Mary Ellen
1978-01-01
Discusses hiring, evaluation, involvement, and competencies of the increasing number of part-time teachers in colleges throughout the country, and the unclear expectations placed on them. Includes a competencies questionnaire for part-time instructors developed at Caldwell Community College and Technical Institute.
Terahertz time domain spectroscopy allows contactless monitoring of grapevine water status
Santesteban, Luis G.; Palacios, Inés; Miranda, Carlos; Iriarte, Juan C.; Royo, José B.; Gonzalo, Ramón
2015-01-01
Agriculture is the sector with the greatest water consumption, since food production is frequently based on crop irrigation. Proper irrigation management requires reliable information on plant water status, but all the plant-based methods to determine it suffer from several inconveniences, mainly caused by the necessity of destructive sampling or of alteration of the plant organ due to contact installation. The aim of this work is to test if terahertz (THz) time domain reflectance measurements made on the grapevine trunk allows contactless monitoring of plant status. The experiments were performed on a potted 14-years-old plant, using a general purpose THz emitter receiver head. Trunk THz time-domain reflection signal proved to be very sensitive to changes in plant water availability, as its pattern follows the trend of soil water content and trunk growth variations. Therefore, it could be used to contactless monitor plant water status. Apart from that, THz reflection signal was observed to respond to light conditions which, according to a specifically designed girdling experiment, was caused by changes in the phloem. This latter results opens a promising field of research for contactless monitoring of phloem activity. PMID:26082791
Terahertz time domain spectroscopy allows contactless monitoring of grapevine water status.
Santesteban, Luis G; Palacios, Inés; Miranda, Carlos; Iriarte, Juan C; Royo, José B; Gonzalo, Ramón
2015-01-01
Agriculture is the sector with the greatest water consumption, since food production is frequently based on crop irrigation. Proper irrigation management requires reliable information on plant water status, but all the plant-based methods to determine it suffer from several inconveniences, mainly caused by the necessity of destructive sampling or of alteration of the plant organ due to contact installation. The aim of this work is to test if terahertz (THz) time domain reflectance measurements made on the grapevine trunk allows contactless monitoring of plant status. The experiments were performed on a potted 14-years-old plant, using a general purpose THz emitter receiver head. Trunk THz time-domain reflection signal proved to be very sensitive to changes in plant water availability, as its pattern follows the trend of soil water content and trunk growth variations. Therefore, it could be used to contactless monitor plant water status. Apart from that, THz reflection signal was observed to respond to light conditions which, according to a specifically designed girdling experiment, was caused by changes in the phloem. This latter results opens a promising field of research for contactless monitoring of phloem activity. PMID:26082791
NASA Astrophysics Data System (ADS)
Cavalcanti, José Rafael; Dumbser, Michael; Motta-Marques, David da; Fragoso Junior, Carlos Ruberto
2015-12-01
In this article we propose a new conservative high resolution TVD (total variation diminishing) finite volume scheme with time-accurate local time stepping (LTS) on unstructured grids for the solution of scalar transport problems, which are typical in the context of water quality simulations. To keep the presentation of the new method as simple as possible, the algorithm is only derived in two space dimensions and for purely convective transport problems, hence neglecting diffusion and reaction terms. The new numerical method for the solution of the scalar transport is directly coupled to the hydrodynamic model of Casulli and Walters (2000) that provides the dynamics of the free surface and the velocity vector field based on a semi-implicit discretization of the shallow water equations. Wetting and drying is handled rigorously by the nonlinear algorithm proposed by Casulli (2009). The new time-accurate LTS algorithm allows a different time step size for each element of the unstructured grid, based on an element-local Courant-Friedrichs-Lewy (CFL) stability condition. The proposed method does not need any synchronization between different time steps of different elements and is by construction locally and globally conservative. The LTS scheme is based on a piecewise linear polynomial reconstruction in space-time using the MUSCL-Hancock method, to obtain second order of accuracy in both space and time. The new algorithm is first validated on some classical test cases for pure advection problems, for which exact solutions are known. In all cases we obtain a very good level of accuracy, showing also numerical convergence results; we furthermore confirm mass conservation up to machine precision and observe an improved computational efficiency compared to a standard second order TVD scheme for scalar transport with global time stepping (GTS). Then, the new LTS method is applied to some more complex problems, where the new scalar transport scheme has also been coupled to
Schönauer, Ria; Kaiser, Anette; Holze, Cathleen; Babilon, Stefanie; Köbberling, Johannes; Riedl, Bernd; Beck-Sickinger, Annette G
2015-12-01
The human adrenomedullin (ADM) is a 52 amino acid peptide hormone belonging to the calcitonin family of peptides, which plays a major role in the development and regulation of cardiovascular and lymphatic systems. For potential use in clinical applications, we aimed to investigate the fate of the peptide ligand after binding and activation of the adrenomedullin receptor (AM1), a heterodimer consisting of the calcitonin receptor-like receptor (CLR), a G protein-coupled receptor, associated with the receptor activity-modifying protein 2 (RAMP2). Full length and N-terminally shortened ADM peptides were synthesized using Fmoc/tBu solid phase peptide synthesis and site-specifically labeled with the fluorophore carboxytetramethylrhodamine (Tam) either by amide bond formation or copper(I)-catalyzed azide alkyne cycloaddition. For the first time, Tam-labeled ligands allowed the observation of co-internalization of the whole ligand-receptor complex in living cells co-transfected with fluorescent fusion proteins of CLR and RAMP2. Application of a fluorescent probe to track lysosomal compartments revealed that ADM together with the CLR/RAMP2-complex is routed to the degradative pathway. Moreover, we found that the N-terminus of ADM is not a crucial component of the peptide sequence in terms of AM1 internalization behavior. PMID:26767744
Write a Research Paper One Step at a Time: Research Writing Guide.
ERIC Educational Resources Information Center
Evans, Helen, Ed.
Intended to supplement the textbook series "Houghton Mifflin English Grammar and Composition" and to offer students and classroom teachers in the secondary schools a review of research writing, this guide outlines a step-by-step process allowing for thorough student comprehension and comfort with the application of basic research and writing…
Time-step limits for a Monte Carlo Compton-scattering method
Densmore, Jeffery D; Warsa, James S; Lowrie, Robert B
2009-01-01
We perform a stability analysis of a Monte Carlo method for simulating the Compton scattering of photons by free electron in high energy density applications and develop time-step limits that avoid unstable and oscillatory solutions. Implementing this Monte Carlo technique in multi physics problems typically requires evaluating the material temperature at its beginning-of-time-step value, which can lead to this undesirable behavior. With a set of numerical examples, we demonstrate the efficacy of our time-step limits.
NASA Technical Reports Server (NTRS)
Chao, W. C.
1982-01-01
With appropriate modifications, a recently proposed explicit-multiple-time-step scheme (EMTSS) is incorporated into the UCLA model. In this scheme, the linearized terms in the governing equations that generate the gravity waves are split into different vertical modes. Each mode is integrated with an optimal time step, and at periodic intervals these modes are recombined. The other terms are integrated with a time step dictated by the CFL condition for low-frequency waves. This large time step requires a special modification of the advective terms in the polar region to maintain stability. Test runs for 72 h show that EMTSS is a stable, efficient and accurate scheme.
Importance of variable time-step algorithms in spatial kinetics calculations
Aviles, B.N.
1994-12-31
The use of spatial kinetics codes in conjunction with advanced thermal-hydraulics codes is becoming more widespread as better methods and faster computers appear. The integrated code packages are being used for routine nuclear power plant design and analysis, including simulations with instrumentation and control systems initiating system perturbations such as rod motion and scrams. As a result, it is important to include a robust variable time-step algorithm that can accurately and efficiently follow widely varying plant neutronic behavior. This paper describes the variable time-step algorithm in SPANDEX and compares the automatic time-step scheme with a more traditional fixed time-step scheme.
Large time-step stability of explicit one-dimensional advection schemes
NASA Technical Reports Server (NTRS)
Leonard, B. P.
1993-01-01
There is a wide-spread belief that most explicit one-dimensional advection schemes need to satisfy the so-called 'CFL condition' - that the Courant number, c = udelta(t)/delta(x), must be less than or equal to one, for stability in the von Neumann sense. This puts severe limitations on the time-step in high-speed, fine-grid calculations and is an impetus for the development of implicit schemes, which often require less restrictive time-step conditions for stability, but are more expensive per time-step. However, it turns out that, at least in one dimension, if explicit schemes are formulated in a consistent flux-based conservative finite-volume form, von Neumann stability analysis does not place any restriction on the allowable Courant number. Any explicit scheme that is stable for c is less than 1, with a complex amplitude ratio, G(c), can be easily extended to arbitrarily large c. The complex amplitude ratio is then given by exp(- (Iota)(Nu)(Theta)) G(delta(c)), where N is the integer part of c, and delta(c) = c - N (less than 1); this is clearly stable. The CFL condition is, in fact, not a stability condition at all, but, rather, a 'range restriction' on the 'pieces' in a piece-wise polynomial interpolation. When a global view is taken of the interpolation, the need for a CFL condition evaporates. A number of well-known explicit advection schemes are considered and thus extended to large delta(t). The analysis also includes a simple interpretation of (large delta(t)) total-variation-diminishing (TVD) constraints.
Convergence Acceleration for Multistage Time-Stepping Schemes
NASA Technical Reports Server (NTRS)
Swanson, R. C.; Turkel, Eli L.; Rossow, C-C; Vasta, V. N.
2006-01-01
The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 could be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. Numerical dissipation operators (based on the Roe scheme, a matrix formulation, and the CUSP scheme) as well as the number of RK stages are considered in evaluating the RK/implicit scheme. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. In two dimensions, turbulent flows over an airfoil at subsonic and transonic conditions are computed. The effects of mesh cell aspect ratio on convergence are investigated for Reynolds numbers between 5.7 x 10(exp 6) and 100.0 x 10(exp 6). Results are also obtained for a transonic wing flow. For both 2-D and 3-D problems, the computational time of a well-tuned standard RK scheme is reduced at least a factor of four.
A Dynamic Era-Based Time-Symmetric Block Time-Step Algorithm with Parallel Implementations
NASA Astrophysics Data System (ADS)
Kaplan, Murat; Saygin, Hasan
2012-06-01
The time-symmetric block time-step (TSBTS) algorithm is a newly developed efficient scheme for N-body integrations. It is constructed on an era-based iteration. In this work, we re-designed the TSBTS integration scheme with a dynamically changing era size. A number of numerical tests were performed to show the importance of choosing the size of the era, especially for long-time integrations. Our second aim was to show that the TSBTS scheme is as suitable as previously known schemes for developing parallel N-body codes. In this work, we relied on a parallel scheme using the copy algorithm for the time-symmetric scheme. We implemented a hybrid of data and task parallelization for force calculation to handle load balancing problems that can appear in practice. Using the Plummer model initial conditions for different numbers of particles, we obtained the expected efficiency and speedup for a small number of particles. Although parallelization of the direct N-body codes is negatively affected by the communication/calculation ratios, we obtained good load-balanced results. Moreover, we were able to conserve the advantages of the algorithm (e.g., energy conservation for long-term simulations).
NASA Astrophysics Data System (ADS)
Hirthe, E. M.; Graf, T.
2012-04-01
Fluid density variations occur due to changes in the solute concentration, temperature and pressure of groundwater. Examples are interaction between freshwater and seawater, radioactive waste disposal, groundwater contamination, and geothermal energy production. The physical coupling between flow and transport introduces non-linearity in the governing mathematical equations, such that solving variable-density flow problems typically requires very long computational time. Computational efficiency can be attained through the use of adaptive time-stepping schemes. The aim of this work is therefore to apply a non-iterative adaptive time-stepping scheme based on local truncation error in variable-density flow problems. That new scheme is implemented into the code of the HydroGeoSphere model (Therrien et al., 2011). The new time-stepping scheme is applied to the Elder (1967) and the Shikaze et al. (1998) problem of free convection in porous and fractured-porous media, respectively. Numerical simulations demonstrate that non-iterative time-stepping based on local truncation error control fully automates the time step size and efficiently limits the temporal discretization error to the user-defined tolerance. Results of the Elder problem show that the new time-stepping scheme presented here is significantly more efficient than uniform time-stepping when high accuracy is required. Results of the Shikaze problem reveal that the new scheme is considerably faster than conventional time-stepping where time step sizes are either constant or controlled by absolute head/concentration changes. Future research will focus on the application of the new time-stepping scheme to variable-density flow in complex real-world fractured-porous rock.
An adaptive time-stepping strategy for solving the phase field crystal model
Zhang, Zhengru; Ma, Yuan; Qiao, Zhonghua
2013-09-15
In this work, we will propose an adaptive time step method for simulating the dynamics of the phase field crystal (PFC) model. The numerical simulation of the PFC model needs long time to reach steady state, and then large time-stepping method is necessary. Unconditionally energy stable schemes are used to solve the PFC model. The time steps are adaptively determined based on the time derivative of the corresponding energy. It is found that the use of the proposed time step adaptivity cannot only resolve the steady state solution, but also the dynamical development of the solution efficiently and accurately. The numerical experiments demonstrate that the CPU time is significantly saved for long time simulations.
An Explicit Super-Time-Stepping Scheme for Non-Symmetric Parabolic Problems
NASA Astrophysics Data System (ADS)
Gurski, K. F.; O'Sullivan, S.
2010-09-01
Explicit numerical methods for the solution of a system of differential equations may suffer from a time step size that approaches zero in order to satisfy stability conditions. When the differential equations are dominated by a skew-symmetric component, the problem is that the real eigenvalues are dominated by imaginary eigenvalues. We compare results for stable time step limits for the super-time-stepping method of Alexiades, Amiez, and Gremaud (super-time-stepping methods belong to the Runge-Kutta-Chebyshev class) and a new method modeled on a predictor-corrector scheme with multiplicative operator splitting. This new explicit method increases stability of the original super-time-stepping whenever the skew-symmetric term is nonzero, which occurs in particular convection-diffusion problems and more generally when the iteration matrix represents a nonlinear operator. The new method is stable for skew symmetric dominated systems where the regular super-time-stepping scheme fails. This method is second order in time (may be increased by Richardson extrapolation) and the spatial order is determined by the user's choice of discretization scheme. We present a comparison between the two super-time-stepping methods to show the speed up available for any non-symmetric system using the nearly symmetric Black-Scholes equation as an example.
NASA Astrophysics Data System (ADS)
Tsekouras, Konstantinos; Presse, Steve
Counting of photobleaching steps is of importance in the investigation of many open problems in biophysics. Current methods of counting photo- bleaching steps cannot directly account for fluorophore photophysical behaviors such as fluorophore self-quenching, blinking and flickering. Our Bayesian approach to the counting problem allows for fluorophore blinking and reactivation as well as for multiple simultaneous photobleaching events and is neither computational resource- nor time- heavy. We detail the method's applicability and limitations and present examples of application in photobleach event counting.
Entosis allows timely elimination of the luminal epithelial barrier for embryo implantation.
Li, Yingju; Sun, Xiaofei; Dey, Sudhansu K
2015-04-21
During implantation, uterine luminal epithelial (LE) cells first interact with the blastocyst trophectoderm. Within 30 hr after the initiation of attachment, LE cells surrounding the blastocyst in the implantation chamber (crypt) disappear, allowing trophoblast cells to make direct physical contact with the underneath stroma for successful implantation. The mechanism for the extraction of LE cells was thought to be mediated by apoptosis. Here, we show that LE cells in direct contact with the blastocyst are endocytosed by trophoblast cells by adopting the nonapoptotic cell-in-cell invasion process (entosis) in the absence of caspase 3 activation. Our in vivo observations were reinforced by the results of co-culture experiments with primary uterine epithelial cells with trophoblast stem cells or blastocysts showing internalization of epithelial cells by trophoblasts. We have identified entosis as a mechanism to remove LE cells by trophoblast cells in implantation, conferring a role for entosis in an important physiological process. PMID:25865893
Consistency of internal fluxes in a hydrological model running at multiple time steps
NASA Astrophysics Data System (ADS)
Ficchi, Andrea; Perrin, Charles; Andréassian, Vazken
2016-04-01
Improving hydrological models remains a difficult task and many ways can be explored, among which one can find the improvement of spatial representation, the search for more robust parametrization, the better formulation of some processes or the modification of model structures by trial-and-error procedure. Several past works indicate that model parameters and structure can be dependent on the modelling time step, and there is thus some rationale in investigating how a model behaves across various modelling time steps, to find solutions for improvements. Here we analyse the impact of data time step on the consistency of the internal fluxes of a rainfall-runoff model run at various time steps, by using a large data set of 240 catchments. To this end, fine time step hydro-climatic information at sub-hourly resolution is used as input of a parsimonious rainfall-runoff model (GR) that is run at eight different model time steps (from 6 minutes to one day). The initial structure of the tested model (i.e. the baseline) corresponds to the daily model GR4J (Perrin et al., 2003), adapted to be run at variable sub-daily time steps. The modelled fluxes considered are interception, actual evapotranspiration and intercatchment groundwater flows. Observations of these fluxes are not available, but the comparison of modelled fluxes at multiple time steps gives additional information for model identification. The joint analysis of flow simulation performance and consistency of internal fluxes at different time steps provides guidance to the identification of the model components that should be improved. Our analysis indicates that the baseline model structure is to be modified at sub-daily time steps to warrant the consistency and realism of the modelled fluxes. For the baseline model improvement, particular attention is devoted to the interception model component, whose output flux showed the strongest sensitivity to modelling time step. The dependency of the optimal model
Automatic Time Stepping with Global Error Control for Groundwater Flow Models
Tang, Guoping
2008-09-01
An automatic time stepping with global error control is proposed for the time integration of the diffusion equation to simulate groundwater flow in confined aquifers. The scheme is based on an a posteriori error estimate for the discontinuous Galerkin (dG) finite element methods. A stability factor is involved in the error estimate and it is used to adapt the time step and control the global temporal error for the backward difference method. The stability factor can be estimated by solving a dual problem. The stability factor is not sensitive to the accuracy of the dual solution and the overhead computational cost can be minimized by solving the dual problem using large time steps. Numerical experiments are conducted to show the application and the performance of the automatic time stepping scheme. Implementation of the scheme can lead to improvement in accuracy and efficiency for groundwater flow models.
40 CFR 97.141 - Timing requirements for CAIR NOX allowance allocations.
Code of Federal Regulations, 2012 CFR
2012-07-01
...), for the control periods in 2009, 2010, 2011, 2012, 2013, and 2014. (b) By July 31, 2011 and July 31 of... 40 Protection of Environment 22 2012-07-01 2012-07-01 false Timing requirements for CAIR NOX... applicable deadline for determination under this paragraph. (c) By July 31, 2009 and July 31 of each...
40 CFR 97.141 - Timing requirements for CAIR NOX allowance allocations.
Code of Federal Regulations, 2013 CFR
2013-07-01
...), for the control periods in 2009, 2010, 2011, 2012, 2013, and 2014. (b) By July 31, 2011 and July 31 of... 40 Protection of Environment 22 2013-07-01 2013-07-01 false Timing requirements for CAIR NOX... applicable deadline for determination under this paragraph. (c) By July 31, 2009 and July 31 of each...
NASA Astrophysics Data System (ADS)
Pal, Robert; Beeby, Andrew
2014-09-01
An inverted microscope has been adapted to allow time-gated imaging and spectroscopy to be carried out on samples containing responsive lanthanide probes. The adaptation employs readily available components, including a pulsed light source, time-gated camera, spectrometer and photon counting detector, allowing imaging, emission spectroscopy and lifetime measurements. Each component is controlled by a suite of software written in LabVIEW and is powered via conventional USB ports.
Real-time quantum trajectories for classically allowed dynamics in strong laser fields
NASA Astrophysics Data System (ADS)
Plimak, L. I.; Ivanov, Misha Yu.
2015-10-01
Both the physical picture of the dynamics of atoms and molecules in intense infrared fields and its theoretical description use the concept of electron trajectories. Here, we address a key question which arises in this context: Are distinctly quantum features of these trajectories, such as the complex-valued coordinates, physically relevant in the classically allowed region of phase space, and what is their origin? First, we argue that solutions of classical equations of motion can account for quantum effects. To this end, we construct an exact solution to the classical Hamilton-Jacobi equation which accounts for dynamics of the wave packet, and show that this solution is physically correct in the limit ?. Second, we show that imaginary components of classical trajectories are directly linked to the finite size of the initial wave packet in momentum space. This way, if the electronic wave packet produced by optical tunnelling in strong infrared fields is localised both in coordinate and momentum, its motion after tunnelling ipso facto cannot be described with purely classical trajectories - in contrast to popular models in the literature.
Time Gating of Chloroplast Autofluorescence Allows Clearer Fluorescence Imaging In Planta
Kodama, Yutaka
2016-01-01
Chloroplast, an organelle facilitating photosynthesis, exhibits strong autofluorescence, which is an undesired background signal that restricts imaging experiments with exogenous fluorophore in plants. In this study, the autofluorescence was characterized in planta under confocal laser microscopy, and it was found that the time-gated imaging technique completely eliminates the autofluorescence. As a demonstration of the technique, a clearer signal of fluorescent protein-tagged phototropin, a blue-light photoreceptor localized at the chloroplast periphery, was visualized in planta. PMID:27027881
Time Gating of Chloroplast Autofluorescence Allows Clearer Fluorescence Imaging In Planta.
Kodama, Yutaka
2016-01-01
Chloroplast, an organelle facilitating photosynthesis, exhibits strong autofluorescence, which is an undesired background signal that restricts imaging experiments with exogenous fluorophore in plants. In this study, the autofluorescence was characterized in planta under confocal laser microscopy, and it was found that the time-gated imaging technique completely eliminates the autofluorescence. As a demonstration of the technique, a clearer signal of fluorescent protein-tagged phototropin, a blue-light photoreceptor localized at the chloroplast periphery, was visualized in planta. PMID:27027881
Cherenkov Video Imaging Allows for the First Visualization of Radiation Therapy in Real Time
Jarvis, Lesley A.; Zhang, Rongxiao; Gladstone, David J.; Jiang, Shudong; Hitchcock, Whitney; Friedman, Oscar D.; Glaser, Adam K.; Jermyn, Michael; Pogue, Brian W.
2014-07-01
Purpose: To determine whether Cherenkov light imaging can visualize radiation therapy in real time during breast radiation therapy. Methods and Materials: An intensified charge-coupled device (CCD) camera was synchronized to the 3.25-μs radiation pulses of the clinical linear accelerator with the intensifier set × 100. Cherenkov images were acquired continuously (2.8 frames/s) during fractionated whole breast irradiation with each frame an accumulation of 100 radiation pulses (approximately 5 monitor units). Results: The first patient images ever created are used to illustrate that Cherenkov emission can be visualized as a video during conditions typical for breast radiation therapy, even with complex treatment plans, mixed energies, and modulated treatment fields. Images were generated correlating to the superficial dose received by the patient and potentially the location of the resulting skin reactions. Major blood vessels are visible in the image, providing the potential to use these as biological landmarks for improved geometric accuracy. The potential for this system to detect radiation therapy misadministrations, which can result from hardware malfunction or patient positioning setup errors during individual fractions, is shown. Conclusions: Cherenkoscopy is a unique method for visualizing surface dose resulting in real-time quality control. We propose that this system could detect radiation therapy errors in everyday clinical practice at a time when these errors can be corrected to result in improved safety and quality of radiation therapy.
Stability analysis and time-step limits for a Monte Carlo Compton-scattering method
Densmore, Jeffery D. Warsa, James S. Lowrie, Robert B.
2010-05-20
A Monte Carlo method for simulating Compton scattering in high energy density applications has been presented that models the photon-electron collision kinematics exactly [E. Canfield, W.M. Howard, E.P. Liang, Inverse Comptonization by one-dimensional relativistic electrons, Astrophys. J. 323 (1987) 565]. However, implementing this technique typically requires an explicit evaluation of the material temperature, which can lead to unstable and oscillatory solutions. In this paper, we perform a stability analysis of this Monte Carlo method and develop two time-step limits that avoid undesirable behavior. The first time-step limit prevents instabilities, while the second, more restrictive time-step limit avoids both instabilities and nonphysical oscillations. With a set of numerical examples, we demonstrate the efficacy of these time-step limits.
Stance time and step width variability have unique contributing impairments in older persons.
Brach, Jennifer S; Studenski, Stephanie; Perera, Subashan; VanSwearingen, Jessie M; Newman, Anne B
2008-04-01
Gait variability may have multiple causes. We hypothesized that central nervous system (CNS) impairments would affect motor control and be manifested as increased stance time and step length variability, while sensory impairments would affect balance and be manifested as increased step width variability. Older adults (mean+/-standard deviation (S.D.) age=79.4+/-4.1, n=558) from the Pittsburgh site of the Cardiovascular Health Study participated. The S.D. across steps was the indicator of gait variability, determined for three gait measures, step length, stance time and step width, using a computerized walkway. Impairment measures included CNS function (modified mini-mental state examination, Trails A and B, Digit Symbol Substitution, finger tapping), sensory function (lower extremity (LE) vibration, vision), strength (grip strength, repeated chair stands), mood, and LE pain. Linear regression models were fit for the three gait variability characteristics using impairment measures as independent variables, adjusted for age, race, gender, and height. Analyses were repeated stratified by gait speed. All measures of CNS impairment were directly related to stance time variability (p<0.01), with increased CNS impairment associated with increased stance time variability. CNS impairments were not related to step length or width variability. Both sensory impairments were inversely related to step width (p<0.01) but not step length or stance time variability. CNS impairments affected stance time variability especially in slow walkers while sensory impairments affected step width variability in fast walkers. Specific patterns of gait variability may imply different underlying causes. Types of gait variability should be specified. Interventions may be targeted at specific types of gait variability. PMID:17632004
Egea, Emilie; Mérigot, Bastien; Mahé-Bézac, Chantal; Féral, Jean-Pierre; Chenuil, Anne
2011-01-01
Echinocardium cordatum had long been considered as cosmopolitan, but molecular data revealed it is a complex of cryptic species, with two non-hybridizing species (B1 & B2) in the Mediterranean Sea living in syntopy with Echinocardium mediterraneum. Histological analyses of the gonads from a 17-month sampling period revealed a statistically significant time lag between the Maturity Indices of E. cordatum and E. mediterraneum. The main environmental stimulus may be different for the two nominal species, possibly seawater temperature for E. cordatum and chlorophyll a concentration for E. mediterraneum. Within the E. cordatum complex, spawning timing and synchrony are different according to major geographic areas (Atlantic/Pacific/Mediterranean) and/or the corresponding genetic subdivision [A/P/(B1 & B2)]. In contrast, the effects of temperature on the reproductive cycle seem rather to mirror the genetic lineages than environmental similarities of the different localities. Between the sister species (B1 & B2) no differences could be detected, maybe due to small sample sizes. PMID:21262482
Hecht, Nils; Woitzik, Johannes; König, Susanne; Horn, Peter; Vajkoczy, Peter
2013-01-01
Currently, there is no adequate technique for intraoperative monitoring of cerebral blood flow (CBF). To evaluate laser speckle imaging (LSI) for assessment of relative CBF, LSI was performed in 30 patients who underwent direct surgical revascularization for treatment of arteriosclerotic cerebrovascular disease (ACVD), Moyamoya disease (MMD), or giant aneurysms, and in 8 control patients who underwent intracranial surgery for reasons other than hemodynamic compromise. The applicability and sensitivity of LSI was investigated through baseline perfusion and CO2 reactivity testing. The dynamics of LSI were assessed during bypass test occlusion and flow initiation procedures. Laser speckle imaging permitted robust (pseudo-) quantitative assessment of relative microcirculatory flow and standard bypass grafting resulted in significantly higher postoperative baseline perfusion values in ACVD and MMD. The applicability and sensitivity of LSI was shown by a significantly reduced CO2 reactivity in ACVD (9.6±9%) and MMD (8.5±8%) compared with control (31.2±5% P<0.0001). In high- and intermediate-flow bypass patients, LSI was characterized by a dynamic real-time response to acute perfusion changes and ultimately confirmed a sufficient flow substitution through the bypass graft. Thus, LSI can be used for sensitive and continuous, non-invasive real-time visualization and measurement of relative cortical CBF in excellent spatial–temporal resolution. PMID:23512134
An Adaptive Fourier Filter for Relaxing Time Stepping Constraints for Explicit Solvers
Gelb, Anne; Archibald, Richard K
2015-01-01
Filtering is necessary to stabilize piecewise smooth solutions. The resulting diffusion stabilizes the method, but may fail to resolve the solution near discontinuities. Moreover, high order filtering still requires cost prohibitive time stepping. This paper introduces an adaptive filter that controls spurious modes of the solution, but is not unnecessarily diffusive. Consequently we are able to stabilize the solution with larger time steps, but also take advantage of the accuracy of a high order filter.
Database Integration: An Intial Step Towards the Deep-Time Data Infrastructure
NASA Astrophysics Data System (ADS)
Kolankowski, S. M.; Fox, P. A.; Ma, X.
2015-12-01
As our knowledge of Earth's geologic history grows, we require more robust methods of sharing immense amounts of data. Various databases across numerous disciplines have been widely utilized to offer extensive information on very specific pieces of both Earth's history and its current state, ie. fossil record, rock composition, proteins, etc. In order to gain a deeper understanding of our planet's past we must combine the resources present in our online communities. These databases could be a powerful force in identifying previously unseen correlations if used in tandem rather than as separate entities. Creating a unifying site that provides links to these databases will aid in our ability as a collaborative scientific community to utilize our findings on a larger scale. The Deep-Time Data Infrastructure is currently underway as part of a larger effort to accomplish this goal. DTDI works not to build a new database, but to integrate existing resources. This research is the beginning step in the DTDI program. To create this infrastructure, all current geologic and related databases had to be identified and their schema recorded. Using variables from their combined records, we are able to determine the best way to integrate them using common factors. The Deep-Time Data Infrastructure will allow geoscientists to bridge gaps in data and further our understanding of our planet's history.
NASA Astrophysics Data System (ADS)
Clark, Martyn P.; Kavetski, Dmitri
2010-10-01
A major neglected weakness of many current hydrological models is the numerical method used to solve the governing model equations. This paper thoroughly evaluates several classes of time stepping schemes in terms of numerical reliability and computational efficiency in the context of conceptual hydrological modeling. Numerical experiments are carried out using 8 distinct time stepping algorithms and 6 different conceptual rainfall-runoff models, applied in a densely gauged experimental catchment, as well as in 12 basins with diverse physical and hydroclimatic characteristics. Results show that, over vast regions of the parameter space, the numerical errors of fixed-step explicit schemes commonly used in hydrology routinely dwarf the structural errors of the model conceptualization. This substantially degrades model predictions, but also, disturbingly, generates fortuitously adequate performance for parameter sets where numerical errors compensate for model structural errors. Simply running fixed-step explicit schemes with shorter time steps provides a poor balance between accuracy and efficiency: in some cases daily-step adaptive explicit schemes with moderate error tolerances achieved comparable or higher accuracy than 15 min fixed-step explicit approximations but were nearly 10 times more efficient. From the range of simple time stepping schemes investigated in this work, the fixed-step implicit Euler method and the adaptive explicit Heun method emerge as good practical choices for the majority of simulation scenarios. In combination with the companion paper, where impacts on model analysis, interpretation, and prediction are assessed, this two-part study vividly highlights the impact of numerical errors on critical performance aspects of conceptual hydrological models and provides practical guidelines for robust numerical implementation.
Mailey, Emily L.; Gothe, Neha P.; Wójcicki, Thomas R.; Szabo, Amanda N.; Olson, Erin A.; Mullen, Sean P.; Fanning, Jason T.; Motl, Robert W.; McAuley, Edward
2013-01-01
The criteria one uses to reduce accelerometer data can profoundly influence the interpretation of research outcomes. The purpose of this study was to examine the influence of three different interruption periods (i.e., 20, 30, and 60 minutes) on the amount of data retained for analyses and estimates of sedentary time among older adults. Older adults (N=311; Mage=71.1) wore an accelerometer for seven days and reported wear time on an accelerometer log. Accelerometer data were downloaded and scored using 20, 30, and 60-minute interruption periods. Estimates of wear time derived using each interruption period were compared to self-reported wear time, and descriptive statistics were used to compare estimates of sedentary time. Results showed a longer interruption period (i.e., 60 minutes) yields the largest sample size and the closest approximation of self-reported wear time. A short interruption period (i.e., 20 minutes) is likely to underestimate sedentary time among older adults. PMID:23752299
Modeling solute transport in distribution networks with variable demand and time step sizes.
Peyton, Chad E.; Bilisoly, Roger Lee; Buchberger, Steven G.; McKenna, Sean Andrew; Yarrington, Lane
2004-06-01
The effect of variable demands at short time scales on the transport of a solute through a water distribution network has not previously been studied. We simulate flow and transport in a small water distribution network using EPANET to explore the effect of variable demand on solute transport across a range of hydraulic time step scales from 1 minute to 2 hours. We show that variable demands at short time scales can have the following effects: smoothing of a pulse of tracer injected into a distribution network and increasing the variability of both the transport pathway and transport timing through the network. Variable demands are simulated for these different time step sizes using a previously developed Poisson rectangular pulse (PRP) demand generator that considers demand at a node to be a combination of exponentially distributed arrival times with log-normally distributed intensities and durations. Solute is introduced at a tank and at three different network nodes and concentrations are modeled through the system using the Lagrangian transport scheme within EPANET. The transport equations within EPANET assume perfect mixing of the solute within a parcel of water and therefore physical dispersion cannot occur. However, variation in demands along the solute transport path contribute to both removal and distortion of the injected pulse. The model performance measures examined are the distribution of the Reynolds number, the variation in the center of mass of the solute across time, and the transport path and timing of the solute through the network. Variation in all three performance measures is greatest at the shortest time step sizes. As the scale of the time step increases, the variability in these performance measures decreases. The largest time steps produce results that are inconsistent with the results produced by the smaller time steps.
Halsey, Lewis G; Watkins, David A R; Duggan, Brendan M
2012-01-01
Stairway climbing provides a ubiquitous and inconspicuous method of burning calories. While typically two strategies are employed for climbing stairs, climbing one stair step per stride or two steps per stride, research to date has not clarified if there are any differences in energy expenditure between them. Fourteen participants took part in two stair climbing trials whereby measures of heart rate were used to estimate energy expenditure during stairway ascent at speeds chosen by the participants. The relationship between rate of oxygen consumption ([Formula: see text]) and heart rate was calibrated for each participant using an inclined treadmill. The trials involved climbing up and down a 14.05 m high stairway, either ascending one step per stride or ascending two stair steps per stride. Single-step climbing used 8.5±0.1 kcal min(-1), whereas double step climbing used 9.2±0.1 kcal min(-1). These estimations are similar to equivalent measures in all previous studies, which have all directly measured [Formula: see text] The present study findings indicate that (1) treadmill-calibrated heart rate recordings can be used as a valid alternative to respirometry to ascertain rate of energy expenditure during stair climbing; (2) two step climbing invokes a higher rate of energy expenditure; however, one step climbing is energetically more expensive in total over the entirety of a stairway. Therefore to expend the maximum number of calories when climbing a set of stairs the single-step strategy is better. PMID:23251455
Halsey, Lewis G.; Watkins, David A. R.; Duggan, Brendan M.
2012-01-01
Stairway climbing provides a ubiquitous and inconspicuous method of burning calories. While typically two strategies are employed for climbing stairs, climbing one stair step per stride or two steps per stride, research to date has not clarified if there are any differences in energy expenditure between them. Fourteen participants took part in two stair climbing trials whereby measures of heart rate were used to estimate energy expenditure during stairway ascent at speeds chosen by the participants. The relationship between rate of oxygen consumption () and heart rate was calibrated for each participant using an inclined treadmill. The trials involved climbing up and down a 14.05 m high stairway, either ascending one step per stride or ascending two stair steps per stride. Single-step climbing used 8.5±0.1 kcal min−1, whereas double step climbing used 9.2±0.1 kcal min−1. These estimations are similar to equivalent measures in all previous studies, which have all directly measured The present study findings indicate that (1) treadmill-calibrated heart rate recordings can be used as a valid alternative to respirometry to ascertain rate of energy expenditure during stair climbing; (2) two step climbing invokes a higher rate of energy expenditure; however, one step climbing is energetically more expensive in total over the entirety of a stairway. Therefore to expend the maximum number of calories when climbing a set of stairs the single-step strategy is better. PMID:23251455
Lutz, Barry; Liang, Tinny; Fu, Elain; Ramachandran, Sujatha; Kauffman, Peter; Yager, Paul
2013-01-01
Lateral flow tests (LFTs) are an ingenious format for rapid and easy-to-use diagnostics, but they are fundamentally limited to assay chemistries that can be reduced to a single chemical step. In contrast, most laboratory diagnostic assays rely on multiple timed steps carried out by a human or a machine. Here, we use dissolvable sugar applied to paper to create programmable flow delays and present a paper network topology that uses these time delays to program automated multi-step fluidic protocols. Solutions of sucrose at different concentrations (10-70% of saturation) were added to paper strips and dried to create fluidic time delays spanning minutes to nearly an hour. A simple folding card format employing sugar delays was shown to automate a four-step fluidic process initiated by a single user activation step (folding the card); this device was used to perform a signal-amplified sandwich immunoassay for a diagnostic biomarker for malaria. The cards are capable of automating multi-step assay protocols normally used in laboratories, but in a rapid, low-cost, and easy-to-use format. PMID:23685876
2015-01-01
When simulating molecular systems using deterministic equations of motion (e.g., Newtonian dynamics), such equations are generally numerically integrated according to a well-developed set of algorithms that share commonly agreed-upon desirable properties. However, for stochastic equations of motion (e.g., Langevin dynamics), there is still broad disagreement over which integration algorithms are most appropriate. While multiple desiderata have been proposed throughout the literature, consensus on which criteria are important is absent, and no published integration scheme satisfies all desiderata simultaneously. Additional nontrivial complications stem from simulating systems driven out of equilibrium using existing stochastic integration schemes in conjunction with recently developed nonequilibrium fluctuation theorems. Here, we examine a family of discrete time integration schemes for Langevin dynamics, assessing how each member satisfies a variety of desiderata that have been enumerated in prior efforts to construct suitable Langevin integrators. We show that the incorporation of a novel time step rescaling in the deterministic updates of position and velocity can correct a number of dynamical defects in these integrators. Finally, we identify a particular splitting (related to the velocity Verlet discretization) that has essentially universally appropriate properties for the simulation of Langevin dynamics for molecular systems in equilibrium, nonequilibrium, and path sampling contexts. PMID:24555448
Multi time-step wavefront reconstruction for tomographic adaptive-optics systems.
Ono, Yoshito H; Akiyama, Masayuki; Oya, Shin; Lardiére, Olivier; Andersen, David R; Correia, Carlos; Jackson, Kate; Bradley, Colin
2016-04-01
In tomographic adaptive-optics (AO) systems, errors due to tomographic wavefront reconstruction limit the performance and angular size of the scientific field of view (FoV), where AO correction is effective. We propose a multi time-step tomographic wavefront reconstruction method to reduce the tomographic error by using measurements from both the current and previous time steps simultaneously. We further outline the method to feed the reconstructor with both wind speed and direction of each turbulence layer. An end-to-end numerical simulation, assuming a multi-object AO (MOAO) system on a 30 m aperture telescope, shows that the multi time-step reconstruction increases the Strehl ratio (SR) over a scientific FoV of 10 arc min in diameter by a factor of 1.5-1.8 when compared to the classical tomographic reconstructor, depending on the guide star asterism and with perfect knowledge of wind speeds and directions. We also evaluate the multi time-step reconstruction method and the wind estimation method on the RAVEN demonstrator under laboratory setting conditions. The wind speeds and directions at multiple atmospheric layers are measured successfully in the laboratory experiment by our wind estimation method with errors below 2 ms^{-1}. With these wind estimates, the multi time-step reconstructor increases the SR value by a factor of 1.2-1.5, which is consistent with a prediction from the end-to-end numerical simulation. PMID:27140785
Ejupi, Andreas; Brodie, Matthew; Gschwind, Yves J; Schoene, Daniel; Lord, Stephen; Delbaere, Kim
2014-01-01
Accidental falls remain an important problem in older people. Stepping is a common task to avoid a fall and requires good interplay between sensory functions, central processing and motor execution. Increased choice stepping reaction time has been associated with recurrent falls in older people. The aim of this study was to examine if a sensor-based Exergame Choice Stepping Reaction Time test can successfully discriminate older fallers from non-fallers. The stepping test was conducted in a cohort of 104 community-dwelling older people (mean age: 80.7 ± 7.0 years). Participants were asked to step laterally as quickly as possible after a light stimulus appeared on a TV screen. Spatial and temporal measurements of the lower and upper body were derived from a low-cost and portable 3D-depth sensor (i.e. Microsoft Kinect) and 3D-accelerometer. Fallers had a slower stepping reaction time (970 ± 228 ms vs. 858 ± 123 ms, P = 0.001) and a slower reaction of their upper body (719 ± 289 ms vs. 631 ± 166 ms, P = 0.052) compared to non-fallers. It took fallers significantly longer than non-fallers to recover their balance after initiating the step (2147 ± 800 ms vs. 1841 ± 591 ms, P = 0.029). This study demonstrated that a sensor-based, low-cost and easy to administer stepping test, with the potential to be used in clinical practice or regular unsupervised home assessments, was able to identify significant differences between performances by fallers and non-fallers. PMID:25571596
Time-step limits for a Monte Carlo Compton-scattering method
Densmore, Jeffery D; Warsa, James S; Lowrie, Robert B
2008-01-01
Compton scattering is an important aspect of radiative transfer in high energy density applications. In this process, the frequency and direction of a photon are altered by colliding with a free electron. The change in frequency of a scattered photon results in an energy exchange between the photon and target electron and energy coupling between radiation and matter. Canfield, Howard, and Liang have presented a Monte Carlo method for simulating Compton scattering that models the photon-electron collision kinematics exactly. However, implementing their technique in multiphysics problems that include the effects of radiation-matter energy coupling typically requires evaluating the material temperature at its beginning-of-time-step value. This explicit evaluation can lead to unstable and oscillatory solutions. In this paper, we perform a stability analysis of this Monte Carlo method and present time-step limits that avoid instabilities and nonphysical oscillations by considering a spatially independent, purely scattering radiative-transfer problem. Examining a simplified problem is justified because it isolates the effects of Compton scattering, and existing Monte Carlo techniques can robustly model other physics (such as absorption, emission, sources, and photon streaming). Our analysis begins by simplifying the equations that are solved via Monte Carlo within each time step using the Fokker-Planck approximation. Next, we linearize these approximate equations about an equilibrium solution such that the resulting linearized equations describe perturbations about this equilibrium. We then solve these linearized equations over a time step and determine the corresponding eigenvalues, quantities that can predict the behavior of solutions generated by a Monte Carlo simulation as a function of time-step size and other physical parameters. With these results, we develop our time-step limits. This approach is similar to our recent investigation of time discretizations for the
ERIC Educational Resources Information Center
Burkman, Ernest; And Others
1981-01-01
The simultaneous effects of two instructional variables of degree of teacher direction (student- versus teacher-directed) and time allowed for study and the student variable of general ability are described for student achievement (N=912) in three life-science instructional modules. (Author)
NASA Technical Reports Server (NTRS)
Grantham, William D.; Smith, Paul M.; Person, Lee H., Jr.; Meyer, Robert T.; Tingas, Stephen A.
1987-01-01
A piloted simulation study was conducted to determine the permissible time delay in the flight control system of a 10-percent statically unstable transport airplane during cruise flight conditions. The math model used for the simulation was a derivative Lockheed L-1011 wide-body jet transport. Data were collected and analyzed from a total of 137 cruising flights in both calm- and turbulent-air conditions. Results of this piloted simulation study verify previous findings that show present military specifications for allowable control-system time delay may be too stringent when applied to transport-size airplanes. Also, the degree of handling-qualities degradation due to time delay is shown to be strongly dependent on the source of the time delay in an advanced flight control system. Maximum allowable time delay for each source of time delay in the control system, in addition to a less stringent overall maximum level of time delay, should be considered for large aircraft. Preliminary results also suggest that adverse effects of control-system time delay may be at least partially offset by variations in control gearing. It is recommended that the data base include different airplane baselines, control systems, and piloting tasks with many pilots participating, so that a reasonable set of limits for control-system time delay can be established to replace the military specification limits currently being used.
Exponential time-differencing with embedded Runge–Kutta adaptive step control
Whalen, P.; Brio, M.; Moloney, J.V.
2015-01-01
We have presented the first embedded Runge–Kutta exponential time-differencing (RKETD) methods of fourth order with third order embedding and fifth order with third order embedding for non-Rosenbrock type nonlinear systems. A procedure for constructing RKETD methods that accounts for both order conditions and stability is outlined. In our stability analysis, the fast time scale is represented by a full linear operator in contrast to particular scalar cases considered before. An effective time-stepping strategy based on reducing both ETD function evaluations and rejected steps is described. Comparisons of performance with adaptive-stepping integrating factor (IF) are carried out on a set of canonical partial differential equations: the shock-fronts of Burgers equation, interacting KdV solitons, KS controlled chaos, and critical collapse of two-dimensional NLS.
Muscle contraction: the step-size distance and the impulse-time per ATP.
Worthington, C R; Elliott, G F
1996-02-01
We derive the step-size distance, and the impulse time per ATP split, from a consideration of Hill's energy rate equation coupled with the enthalpy available per ATP split. This definition of step-size distance is model-independent, and is calculated to have a maximum of 17 A at no load and to reduce to zero at isometric tension, since it will depend on the velocity of shortening. We revisit a derivation of Hill's force-velocity equation depending on impulsive forces working against frictional forces and show that this gives a physical meaning to Hill's constants a and b. This is particularly elegant for Hill's constant b, which is directly related to the impulse time; the value of this impulse time is 1/2 ms. The question that muscle contraction may involve overlapping interactions is considered. However, we find that the step-size distance is not dependent on the possibility of overlapping interactions. PMID:8852761
NASA Astrophysics Data System (ADS)
Wang, Yue
A new variable grid-size and time-step finite-difference (FD) method is developed and applied to three different geophysical problems: simulation of tube waves in boreholes, three-dimensional (3-D) ground-motion simulation in sedimentary basin models, and reverse-time migration of multicomponent data. Unlike the conventional FD method, which uses a fixed grid-size and time-step for the entire model region, spatially variable grid-sizes and time-steps are used to achieve the optimal computational efficiency. For tube wave simulations, a fine grid-spacing is used for simulation inside the borehole region, while a coarse grid is used in the exterior region. While the stability condition requires a very fine time step for the fine grid, a variable time-step method provides coarse time steps for simulation in the coarse grid. Variable grid-size and time-step changes are used to achieve both accuracy and efficiency in the simulations. Numerical tests are performed for the Bayou Choctaw salt-flank model with different borehole models. The results show the important borehole effects on the seismic wavefield for a realistic source bandwidth. The combination of variable grid-size and time-step methods reduces computational costs by several orders of magnitude for the borehole models. Viscoelastic 3-D simulations are performed for a three-layer Salt Lake basin model. The near-surface unconsolidated layer is modeled with a fine grid, and the deep part of the model is modeled by a coarse grid. Simulation results show that the 3-D basin features and the shallow layer significantly affect the amplitude and duration time of the ground motion. In the elastic case, the approximation by 2-D modeling is insufficient to simulate the 3-D ground motion response. A basin model without a shallow low-velocity layer underestimates the ground motion duration and cumulative kinetic energy by 50% or more. The simulation of a Bingham Mine blast suggests that a lower S-velocity should be used to
NASA Technical Reports Server (NTRS)
Jameson, A.; Schmidt, Wolfgang; Turkel, Eli
1981-01-01
A new combination of a finite volume discretization in conjunction with carefully designed dissipative terms of third order, and a Runge Kutta time stepping scheme, is shown to yield an effective method for solving the Euler equations in arbitrary geometric domains. The method has been used to determine the steady transonic flow past an airfoil using an O mesh. Convergence to a steady state is accelerated by the use of a variable time step determined by the local Courant member, and the introduction of a forcing term proportional to the difference between the local total enthalpy and its free stream value.
A step in time: Changes in standard-frequency and time-signal broadcasts, 1 January 1972
NASA Technical Reports Server (NTRS)
Chi, A. R.; Fosque, H. S.
1973-01-01
An improved coordinated universal time (UTC) system has been adopted by the International Radio Consultative Committee. It was implemented internationally by the standard-frequency and time-broadcast stations on 1 Jan. 1972. The new UTC system eliminates the frequency offset of 300 parts in 10 to the 10th power between the old UTC and atomic time, thus making the broadcast time interval (the UTC second) constant and defined by the resonant frequency of cesium atoms. The new time scale is kept in synchronism with the rotation of the Earth within plus or minus 0.7 s by step-time adjustments of exactly 1 s, when needed. A time code has been added to the disseminated time signals to permit universal time to be obtained from the broadcasts to the nearest 0.1 s for users requiring such precision. The texts of the International Radio Consultative Committee recommendation and report to implement the new UTC system are given. The coding formats used by various standard time broadcast services to transmit the difference between the universal time (UT1) and the UTC are also given. For users' convenience, worldwide primary VLF and HF transmissions stations, frequencies, and schedules of time emissions are also included. Actual time-step adjustments made by various stations on 1 Jan. 1972, are provided for future reference.
Time-step Considerations in Particle Simulation Algorithms for Coulomb Collisions in Plasmas
Cohen, B I; Dimits, A; Friedman, A; Caflisch, R
2009-10-29
The accuracy of first-order Euler and higher-order time-integration algorithms for grid-based Langevin equations collision models in a specific relaxation test problem is assessed. We show that statistical noise errors can overshadow time-step errors and argue that statistical noise errors can be conflated with time-step effects. Using a higher-order integration scheme may not achieve any benefit in accuracy for examples of practical interest. We also investigate the collisional relaxation of an initial electron-ion relative drift and the collisional relaxation to a resistive steady-state in which a quasi-steady current is driven by a constant applied electric field, as functions of the time step used to resolve the collision processes using binary and grid-based, test-particle Langevin equations models. We compare results from two grid-based Langevin equations collision algorithms to results from a binary collision algorithm for modeling electronion collisions. Some guidance is provided regarding how large a time step can be used compared to the inverse of the characteristic collision frequency for specific relaxation processes.
The large discretization step method for time-dependent partial differential equations
NASA Technical Reports Server (NTRS)
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.
Simulating diffusion processes in discontinuous media: A numerical scheme with constant time steps
Lejay, Antoine; Pichot, Geraldine
2012-08-30
In this article, we propose new Monte Carlo techniques for moving a diffusive particle in a discontinuous media. In this framework, we characterize the stochastic process that governs the positions of the particle. The key tool is the reduction of the process to a Skew Brownian motion (SBM). In a zone where the coefficients are locally constant on each side of the discontinuity, the new position of the particle after a constant time step is sampled from the exact distribution of the SBM process at the considered time. To do so, we propose two different but equivalent algorithms: a two-steps simulation with a stop at the discontinuity and a one-step direct simulation of the SBM dynamic. Some benchmark tests illustrate their effectiveness.
Suggestions for CAP-TSD mesh and time-step input parameters
NASA Technical Reports Server (NTRS)
Bland, Samuel R.
1991-01-01
Suggestions for some of the input parameters used in the CAP-TSD (Computational Aeroelasticity Program-Transonic Small Disturbance) computer code are presented. These parameters include those associated with the mesh design and time step. The guidelines are based principally on experience with a one-dimensional model problem used to study wave propagation in the vertical direction.
Adaptive time stepping algorithm for Lagrangian transport models: Theory and idealised test cases
NASA Astrophysics Data System (ADS)
Shah, Syed Hyder Ali Muttaqi; Heemink, Arnold Willem; Gräwe, Ulf; Deleersnijder, Eric
2013-08-01
Random walk simulations have an excellent potential in marine and oceanic modelling. This is essentially due to their relative simplicity and their ability to represent advective transport without being plagued by the deficiencies of the Eulerian methods. The physical and mathematical foundations of random walk modelling of turbulent diffusion have become solid over the years. Random walk models rest on the theory of stochastic differential equations. Unfortunately, the latter and the related numerical aspects have not attracted much attention in the oceanic modelling community. The main goal of this paper is to help bridge the gap by developing an efficient adaptive time stepping algorithm for random walk models. Its performance is examined on two idealised test cases of turbulent dispersion; (i) pycnocline crossing and (ii) non-flat isopycnal diffusion, which are inspired by shallow-sea dynamics and large-scale ocean transport processes, respectively. The numerical results of the adaptive time stepping algorithm are compared with the fixed-time increment Milstein scheme, showing that the adaptive time stepping algorithm for Lagrangian random walk models is more efficient than its fixed step-size counterpart without any loss in accuracy.
Dependence of Hurricane intensity and structures on vertical resolution and time-step size
NASA Astrophysics Data System (ADS)
Zhang, Da-Lin; Wang, Xiaoxue
2003-09-01
In view of the growing interests in the explicit modeling of clouds and precipitation, the effects of varying vertical resolution and time-step sizes on the 72-h explicit simulation of Hurricane Andrew (1992) are studied using the Pennsylvania State University/National Center for Atmospheric Research (PSU/NCAR) mesoscale model (i.e., MM5) with the finest grid size of 6 km. It is shown that changing vertical resolution and time-step size has significant effects on hurricane intensity and inner-core cloud/precipitation, but little impact on the hurricane track. In general, increasing vertical resolution tends to produce a deeper storm with lower central pressure and stronger three-dimensional winds, and more precipitation. Similar effects, but to a less extent, occur when the time-step size is reduced. It is found that increasing the low-level vertical resolution is more efficient in intensifying a hurricane, whereas changing the upper-level vertical resolution has little impact on the hurricane intensity. Moreover, the use of a thicker surface layer tends to produce higher maximum surface winds. It is concluded that the use of higher vertical resolution, a thin surface layer, and smaller time-step sizes, along with higher horizontal resolution, is desirable to model more realistically the intensity and inner-core structures and evolution of tropical storms as well as the other convectively driven weather systems.
Causal-Path Local Time-Stepping in the discontinuous Galerkin method for Maxwell's equations
NASA Astrophysics Data System (ADS)
Angulo, L. D.; Alvarez, J.; Teixeira, F. L.; Pantoja, M. F.; Garcia, S. G.
2014-01-01
We introduce a novel local time-stepping technique for marching-in-time algorithms. The technique is denoted as Causal-Path Local Time-Stepping (CPLTS) and it is applied for two time integration techniques: fourth-order low-storage explicit Runge-Kutta (LSERK4) and second-order Leap-Frog (LF2). The CPLTS method is applied to evolve Maxwell's curl equations using a Discontinuous Galerkin (DG) scheme for the spatial discretization. Numerical results for LF2 and LSERK4 are compared with analytical solutions and the Montseny's LF2 technique. The results show that the CPLTS technique improves the dispersive and dissipative properties of LF2-LTS scheme.
Development of a variable time-step transient NEW code: SPANDEX
Aviles, B.N. )
1993-01-01
This paper describes a three-dimensional, variable time-step transient multigroup diffusion theory code, SPANDEX (space-time nodal expansion method). SPANDEX is based on the static nodal expansion method (NEM) code, NODEX (Ref. 1), and employs a nonlinear algorithm and a fifth-order expansion of the transverse-integrated fluxes. The time integration scheme in SPANDEX is a fourth-order implicit generalized Runge-Kutta method (GRK) with on-line error control and variable time-step selection. This Runge-Kutta method has been applied previously to point kinetics and one-dimensional finite difference transient analysis. This paper describes the application of the Runge-Kutta method to three-dimensional reactor transient analysis in a multigroup NEM code.
Promoting rest using a quiet time innovation in an adult neuroscience step down unit.
Bergner, Tara
2014-01-01
Sleep and rest are fundamental for the restoration of energy needed to recuperate from illness, trauma and surgery. At present hospitals are too noisy to promote rest for patients. A literature search produced research that described how quiet time interventions addressing noise levels have met with positive patient and staff satisfaction, as well as creating a more peaceful and healing environment. In this paper, a description of the importance of quiet time and how a small butfeasible innovation was carried out in an adult neuroscience step down unit in a large tertiary health care facility in Canada is provided. Anecdotal evidence from patients, families, and staff suggests that quiet time may have positive effects for patients, their families, and the adult neuroscience step down unit staff Future research examining the effect of quiet time on patient, family and staff satisfaction and patient healing is necessary. PMID:25638912
NASA Astrophysics Data System (ADS)
Mulder, W. A.; Zhebel, E.; Minisini, S.
2014-02-01
We analyse the time-stepping stability for the 3-D acoustic wave equation, discretized on tetrahedral meshes. Two types of methods are considered: mass-lumped continuous finite elements and the symmetric interior-penalty discontinuous Galerkin method. Combining the spatial discretization with the leap-frog time-stepping scheme, which is second-order accurate and conditionally stable, leads to a fully explicit scheme. We provide estimates of its stability limit for simple cases, namely, the reference element with Neumann boundary conditions, its distorted version of arbitrary shape, the unit cube that can be partitioned into six tetrahedra with periodic boundary conditions and its distortions. The Courant-Friedrichs-Lewy stability limit contains an element diameter for which we considered different options. The one based on the sum of the eigenvalues of the spatial operator for the first-degree mass-lumped element gives the best results. It resembles the diameter of the inscribed sphere but is slightly easier to compute. The stability estimates show that the mass-lumped continuous and the discontinuous Galerkin finite elements of degree 2 have comparable stability conditions, whereas the mass-lumped elements of degree one and three allow for larger time steps.
NASA Astrophysics Data System (ADS)
Zhang, Peng; Zhang, Na; Deng, Yuefan; Bluestein, Danny
2015-03-01
We developed a multiple time-stepping (MTS) algorithm for multiscale modeling of the dynamics of platelets flowing in viscous blood plasma. This MTS algorithm improves considerably the computational efficiency without significant loss of accuracy. This study of the dynamic properties of flowing platelets employs a combination of the dissipative particle dynamics (DPD) and the coarse-grained molecular dynamics (CGMD) methods to describe the dynamic microstructures of deformable platelets in response to extracellular flow-induced stresses. The disparate spatial scales between the two methods are handled by a hybrid force field interface. However, the disparity in temporal scales between the DPD and CGMD that requires time stepping at microseconds and nanoseconds respectively, represents a computational challenge that may become prohibitive. Classical MTS algorithms manage to improve computing efficiency by multi-stepping within DPD or CGMD for up to one order of magnitude of scale differential. In order to handle 3-4 orders of magnitude disparity in the temporal scales between DPD and CGMD, we introduce a new MTS scheme hybridizing DPD and CGMD by utilizing four different time stepping sizes. We advance the fluid system at the largest time step, the fluid-platelet interface at a middle timestep size, and the nonbonded and bonded potentials of the platelet structural system at two smallest timestep sizes. Additionally, we introduce parameters to study the relationship of accuracy versus computational complexities. The numerical experiments demonstrated 3000x reduction in computing time over standard MTS methods for solving the multiscale model. This MTS algorithm establishes a computationally feasible approach for solving a particle-based system at multiple scales for performing efficient multiscale simulations.
Zhang, Peng; Zhang, Na; Deng, Yuefan; Bluestein, Danny
2015-01-01
We developed a multiple time-stepping (MTS) algorithm for multiscale modeling of the dynamics of platelets flowing in viscous blood plasma. This MTS algorithm improves considerably the computational efficiency without significant loss of accuracy. This study of the dynamic properties of flowing platelets employs a combination of the dissipative particle dynamics (DPD) and the coarse-grained molecular dynamics (CGMD) methods to describe the dynamic microstructures of deformable platelets in response to extracellular flow-induced stresses. The disparate spatial scales between the two methods are handled by a hybrid force field interface. However, the disparity in temporal scales between the DPD and CGMD that requires time stepping at microseconds and nanoseconds respectively, represents a computational challenge that may become prohibitive. Classical MTS algorithms manage to improve computing efficiency by multi-stepping within DPD or CGMD for up to one order of magnitude of scale differential. In order to handle 3–4 orders of magnitude disparity in the temporal scales between DPD and CGMD, we introduce a new MTS scheme hybridizing DPD and CGMD by utilizing four different time stepping sizes. We advance the fluid system at the largest time step, the fluid-platelet interface at a middle timestep size, and the nonbonded and bonded potentials of the platelet structural system at two smallest timestep sizes. Additionally, we introduce parameters to study the relationship of accuracy versus computational complexities. The numerical experiments demonstrated 3000x reduction in computing time over standard MTS methods for solving the multiscale model. This MTS algorithm establishes a computationally feasible approach for solving a particle-based system at multiple scales for performing efficient multiscale simulations. PMID:25641983
An implicit time-stepping scheme for rigid body dynamics with Coulomb friction
STEWART,DAVID; TRINKLE,JEFFREY C.
2000-02-15
In this paper a new time-stepping method for simulating systems of rigid bodies is given. Unlike methods which take an instantaneous point of view, the method is based on impulse-momentum equations, and so does not need to explicitly resolve impulsive forces. On the other hand, the method is distinct from previous impulsive methods in that it does not require explicit collision checking and it can handle simultaneous impacts. Numerical results are given for one planar and one three-dimensional example, which demonstrate the practicality of the method, and its convergence as the step size becomes small.
Inertial stochastic dynamics. I. Long-time-step methods for Langevin dynamics
NASA Astrophysics Data System (ADS)
Beard, Daniel A.; Schlick, Tamar
2000-05-01
Two algorithms are presented for integrating the Langevin dynamics equation with long numerical time steps while treating the mass terms as finite. The development of these methods is motivated by the need for accurate methods for simulating slow processes in polymer systems such as two-site intermolecular distances in supercoiled DNA, which evolve over the time scale of milliseconds. Our new approaches refine the common Brownian dynamics (BD) scheme, which approximates the Langevin equation in the highly damped diffusive limit. Our LTID ("long-time-step inertial dynamics") method is based on an eigenmode decomposition of the friction tensor. The less costly integrator IBD ("inertial Brownian dynamics") modifies the usual BD algorithm by the addition of a mass-dependent correction term. To validate the methods, we evaluate the accuracy of LTID and IBD and compare their behavior to that of BD for the simple example of a harmonic oscillator. We find that the LTID method produces the expected correlation structure for Langevin dynamics regardless of the level of damping. In fact, LTID is the only consistent method among the three, with error vanishing as the time step approaches zero. In contrast, BD is accurate only for highly overdamped systems. For cases of moderate overdamping, and for the appropriate choice of time step, IBD is significantly more accurate than BD. IBD is also less computationally expensive than LTID (though both are the same order of complexity as BD), and thus can be applied to simulate systems of size and time scale ranges previously accessible to only the usual BD approach. Such simulations are discussed in our companion paper, for long DNA molecules modeled as wormlike chains.
High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs
NASA Technical Reports Server (NTRS)
Nielsen, Tanner B.; Carpenter, Mark H.; Fisher, Travis C.; Frankel, Steven H.
2014-01-01
This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.
Daily Time Step Refinement of Optimized Flood Control Rule Curves for a Global Warming Scenario
NASA Astrophysics Data System (ADS)
Lee, S.; Fitzgerald, C.; Hamlet, A. F.; Burges, S. J.
2009-12-01
Pacific Northwest temperatures have warmed by 0.8 °C since 1920 and are predicted to further increase in the 21st century. Simulated streamflow timing shifts associated with climate change have been found in past research to degrade water resources system performance in the Columbia River Basin when using existing system operating policies. To adapt to these hydrologic changes, optimized flood control operating rule curves were developed in a previous study using a hybrid optimization-simulation approach which rebalanced flood control and reservoir refill at a monthly time step. For the climate change scenario, use of the optimized flood control curves restored reservoir refill capability without increasing flood risk. Here we extend the earlier studies using a detailed daily time step simulation model applied over a somewhat smaller portion of the domain (encompassing Libby, Duncan, and Corra Linn dams, and Kootenai Lake) to evaluate and refine the optimized flood control curves derived from monthly time step analysis. Moving from a monthly to daily analysis, we found that the timing of flood control evacuation needed adjustment to avoid unintended outcomes affecting Kootenai Lake. We refined the flood rule curves derived from monthly analysis by creating a more gradual evacuation schedule, but kept the timing and magnitude of maximum evacuation the same as in the monthly analysis. After these refinements, the performance at monthly time scales reported in our previous study proved robust at daily time scales. Due to a decrease in July storage deficits, additional benefits such as more revenue from hydropower generation and more July and August outflow for fish augmentation were observed when the optimized flood control curves were used for the climate change scenario.
NASA Astrophysics Data System (ADS)
Kuraz, Michal
2016-06-01
Modelling the transport processes in a vadose zone, e.g. modelling contaminant transport or the effect of the soil water regime on changes in soil structure and composition, plays an important role in predicting the reactions of soil biotopes to anthropogenic activity. Water flow is governed by the quasilinear Richards equation. The paper concerns the implementation of a multi-time-step approach for solving a nonlinear Richards equation. When modelling porous media flow with a Richards equation, due to a possible convection dominance and a convergence of a nonlinear solver, a stable finite element approximation requires accurate temporal and spatial integration. The method presented here enables adaptive domain decomposition algorithm together with a multi-time-step treatment of actively changing subdomains.
NASA Astrophysics Data System (ADS)
Aronoff, H. I.; Leslie, J. J.; Mittleman, A. N.; Holt, S.
1983-11-01
This manual describes a Shared Time Engineering Program (STEP) conducted by the New England Apparel Manufacturers Association (NEAMA) headquartered in Fall River Massachusetts, and funded by the Office of Trade Adjustment Assistance of the U.S. Department of Commerce. It is addressed to industry association executives, industrial engineers and others interested in examining an innovative model of industrial engineering assistance to small plants which might be adapted to their particular needs.
ERIC Educational Resources Information Center
Pitetti, Kenneth H.; Beets, Michael W.; Flaming, Judy
2009-01-01
Pedometer accuracy for steps and activity time during dynamic movement for youth with intellectual disabilities (ID) were examined. Twenty-four youth with ID (13 girls, 13.1 [plus or minus] 3.2 yrs; 11 boys, 14.7 [plus or minus] 2.7 yrs) were videotaped during adapted physical education class while wearing a Walk4Life 2505 pedometer in five…
Chen, Yunjie; Kale, Seyit; Weare, Jonathan; Dinner, Aaron R; Roux, Benoît
2016-04-12
A multiple time-step integrator based on a dual Hamiltonian and a hybrid method combining molecular dynamics (MD) and Monte Carlo (MC) is proposed to sample systems in the canonical ensemble. The Dual Hamiltonian Multiple Time-Step (DHMTS) algorithm is based on two similar Hamiltonians: a computationally expensive one that serves as a reference and a computationally inexpensive one to which the workload is shifted. The central assumption is that the difference between the two Hamiltonians is slowly varying. Earlier work has shown that such dual Hamiltonian multiple time-step schemes effectively precondition nonlinear differential equations for dynamics by reformulating them into a recursive root finding problem that can be solved by propagating a correction term through an internal loop, analogous to RESPA. Of special interest in the present context, a hybrid MD-MC version of the DHMTS algorithm is introduced to enforce detailed balance via a Metropolis acceptance criterion and ensure consistency with the Boltzmann distribution. The Metropolis criterion suppresses the discretization errors normally associated with the propagation according to the computationally inexpensive Hamiltonian, treating the discretization error as an external work. Illustrative tests are carried out to demonstrate the effectiveness of the method. PMID:26918826
Weare, Jonathan; Dinner, Aaron R.; Roux, Benoît
2016-01-01
A multiple time-step integrator based on a dual Hamiltonian and a hybrid method combining molecular dynamics (MD) and Monte Carlo (MC) is proposed to sample systems in the canonical ensemble. The Dual Hamiltonian Multiple Time-Step (DHMTS) algorithm is based on two similar Hamiltonians: a computationally expensive one that serves as a reference and a computationally inexpensive one to which the workload is shifted. The central assumption is that the difference between the two Hamiltonians is slowly varying. Earlier work has shown that such dual Hamiltonian multiple time-step schemes effectively precondition nonlinear differential equations for dynamics by reformulating them into a recursive root finding problem that can be solved by propagating a correction term through an internal loop, analogous to RESPA. Of special interest in the present context, a hybrid MD-MC version of the DHMTS algorithm is introduced to enforce detailed balance via a Metropolis acceptance criterion and ensure consistency with the Boltzmann distribution. The Metropolis criterion suppresses the discretization errors normally associated with the propagation according to the computationally inexpensive Hamiltonian, treating the discretization error as an external work. Illustrative tests are carried out to demonstrate the effectiveness of the method. PMID:26918826
Real-Time, Single-Step Bioassay Using Nanoplasmonic Resonator With Ultra-High Sensitivity
NASA Technical Reports Server (NTRS)
Zhang, Xiang (Inventor); Ellman, Jonathan A. (Inventor); Chen, Fanqing Frank (Inventor); Su, Kai-Hang (Inventor); Wei, Qi-Huo (Inventor); Sun, Cheng (Inventor)
2014-01-01
A nanoplasmonic resonator (NPR) comprising a metallic nanodisk with alternating shielding layer(s), having a tagged biomolecule conjugated or tethered to the surface of the nanoplasmonic resonator for highly sensitive measurement of enzymatic activity. NPRs enhance Raman signals in a highly reproducible manner, enabling fast detection of protease and enzyme activity, such as Prostate Specific Antigen (paPSA), in real-time, at picomolar sensitivity levels. Experiments on extracellular fluid (ECF) from paPSA-positive cells demonstrate specific detection in a complex bio-fluid background in real-time single-step detection in very small sample volumes.
Real-time, single-step bioassay using nanoplasmonic resonator with ultra-high sensitivity
Zhang, Xiang; Ellman, Jonathan A; Chen, Fanqing Frank; Su, Kai-Hang; Wei, Qi-Huo; Sun, Cheng
2014-04-01
A nanoplasmonic resonator (NPR) comprising a metallic nanodisk with alternating shielding layer(s), having a tagged biomolecule conjugated or tethered to the surface of the nanoplasmonic resonator for highly sensitive measurement of enzymatic activity. NPRs enhance Raman signals in a highly reproducible manner, enabling fast detection of protease and enzyme activity, such as Prostate Specific Antigen (paPSA), in real-time, at picomolar sensitivity levels. Experiments on extracellular fluid (ECF) from paPSA-positive cells demonstrate specific detection in a complex bio-fluid background in real-time single-step detection in very small sample volumes.
Moment tensor inversion of waveforms: a two-step time-frequency approach
NASA Astrophysics Data System (ADS)
Vavryčuk, Václav; Kühn, Daniela
2012-09-01
We present a moment tensor inversion of waveforms, which is more robust and yields more stable and more accurate results than standard approaches. The inversion is performed in two steps and combines inversions in time and frequency domains. First, the inversion for the source-time function is performed in the frequency domain using complex spectra. Second, the time-domain inversion for the moment tensor is performed using the source-time function calculated in the first step. In this way, we can consider a realistic, complex source-time function and still keep the final moment tensor inversion linear. Using numerical modelling, we compare the efficiency and accuracy of the proposed approach with standard waveform inversions. We study the sensitivity of the retrieved double-couple and non-double-couple components of the moment tensors to noise in the data, to inaccuracies of the location and of the velocity model, and to the type of the focal mechanism. Finally, the proposed moment tensor inversion is tested on real data observed in a complex 3-D inhomogeneous geological environment: a production blast and a rockburst in the Pyhäsalmi ore mine, Finland.
Toggweiler, Matthias; Adelmann, Andreas; Arbenz, Peter; Yang, Jianjun
2014-09-15
We show that adaptive time stepping in particle accelerator simulation is an enhancement for certain problems. The new algorithm has been implemented in the OPAL (Object Oriented Parallel Accelerator Library) framework. The idea is to adjust the frequency of costly self-field calculations, which are needed to model Coulomb interaction (space charge) effects. In analogy to a Kepler orbit simulation that requires a higher time step resolution at the close encounter, we propose to choose the time step based on the magnitude of the space charge forces. Inspired by geometric integration techniques, our algorithm chooses the time step proportional to a function of the current phase space state instead of calculating a local error estimate like a conventional adaptive procedure. Building on recent work, a more profound argument is given on how exactly the time step should be chosen. An intermediate algorithm, initially built to allow a clearer analysis by introducing separate time steps for external field and self-field integration, turned out to be useful by its own, for a large class of problems.
Garrido, Jesús A; Luque, Niceto R; Tolu, Silvia; D'Angelo, Egidio
2016-08-01
The majority of operations carried out by the brain require learning complex signal patterns for future recognition, retrieval and reuse. Although learning is thought to depend on multiple forms of long-term synaptic plasticity, the way this latter contributes to pattern recognition is still poorly understood. Here, we have used a simple model of afferent excitatory neurons and interneurons with lateral inhibition, reproducing a network topology found in many brain areas from the cerebellum to cortical columns. When endowed with spike-timing dependent plasticity (STDP) at the excitatory input synapses and at the inhibitory interneuron-interneuron synapses, the interneurons rapidly learned complex input patterns. Interestingly, induction of plasticity required that the network be entrained into theta-frequency band oscillations, setting the internal phase-reference required to drive STDP. Inhibitory plasticity effectively distributed multiple patterns among available interneurons, thus allowing the simultaneous detection of multiple overlapping patterns. The addition of plasticity in intrinsic excitability made the system more robust allowing self-adjustment and rescaling in response to a broad range of input patterns. The combination of plasticity in lateral inhibitory connections and homeostatic mechanisms in the inhibitory interneurons optimized mutual information (MI) transfer. The storage of multiple complex patterns in plastic interneuron networks could be critical for the generation of sparse representations of information in excitatory neuron populations falling under their control. PMID:27079422
NASA Astrophysics Data System (ADS)
Yu, Chunxue; Yin, Xin'an; Yang, Zhifeng; Cai, Yanpeng; Sun, Tao
2016-09-01
The time step used in the operation of eco-friendly reservoirs has decreased from monthly to daily, and even sub-daily. The shorter time step is considered a better choice for satisfying downstream environmental requirements because it more closely resembles the natural flow regime. However, little consideration has been given to the influence of different time steps on the ability to simultaneously meet human and environmental flow requirements. To analyze this influence, we used an optimization model to explore the relationships among the time step, environmental flow (e-flow) requirements, and human water needs for a wide range of time steps and e-flow scenarios. We used the degree of hydrologic alteration to evaluate the regime's ability to satisfy the e-flow requirements of riverine ecosystems, and used water supply reliability to evaluate the ability to satisfy human needs. We then applied the model to a case study of China's Tanghe Reservoir. We found four efficient time steps (2, 3, 4, and 5 days), with a remarkably high water supply reliability (around 80%) and a low alteration of the flow regime (<35%). Our analysis of the hydrologic alteration revealed the smallest alteration at time steps ranging from 1 to 7 days. However, longer time steps led to higher water supply reliability to meet human needs under several e-flow scenarios. Our results show that adjusting the time step is a simple way to improve reservoir operation performance to balance human and e-flow needs.
Sensitivity of The High-resolution Wam Model With Respect To Time Step
NASA Astrophysics Data System (ADS)
Kasemets, K.; Soomere, T.
The northern part of the Baltic Proper and its subbasins (Bothnian Sea, the Gulf of Finland, Moonsund) serve as a challenge for wave modellers. In difference from the southern and the eastern parts of the Baltic Sea, their coasts are highly irregular and contain many peculiarities with the characteristic horizontal scale of the order of a few kilometres. For example, the northern coast of the Gulf of Finland is extremely ragged and contains a huge number of small islands. Its southern coast is more or less regular but has up to 50m high cliff that is frequently covered by high forests. The area also contains numerous banks that have water depth a couple of meters and that may essentially modify wave properties near the banks owing to topographical effects. This feature suggests that a high-resolution wave model should be applied for the region in question, with a horizontal resolution of an order of 1 km or even less. According to the Courant-Friedrich-Lewy criterion, the integration time step for such models must be of the order of a few tens of seconds. A high-resolution WAM model turns out to be fairly sensitive with respect to the particular choice of the time step. In our experiments, a medium-resolution model for the whole Baltic Sea was used, with the horizontal resolution 3 miles (3' along latitudes and 6' along longitudes) and the angular resolution 12 directions. The model was run with steady wind blowing 20 m/s from different directions and with two time steps (1 and 3 minutes). For most of the wind directions, the rms. difference of significant wave heights calculated with differ- ent time steps did not exceed 10 cm and typically was of the order of a few per cents. The difference arose within a few tens of minutes and generally did not increase in further computations. However, in the case of the north wind, the difference increased nearly monotonously and reached 25-35 cm (10-15%) within three hours of integra- tion whereas mean of significant wave
BIOMAP A Daily Time Step, Mechanistic Model for the Study of Ecosystem Dynamics
NASA Astrophysics Data System (ADS)
Wells, J. R.; Neilson, R. P.; Drapek, R. J.; Pitts, B. S.
2010-12-01
BIOMAP simulates competition between two Plant Functional Types (PFT) at any given point in the conterminous U.S. using a time series of daily temperature (mean, minimum, maximum), precipitation, humidity, light and nutrients, with PFT-specific rooting within a multi-layer soil. The model employs a 2-layer canopy biophysics, Farquhar photosynthesis, the Beer-Lambert Law for light attenuation and a mechanistic soil hydrology. In essence, BIOMAP is a re-built version of the biogeochemistry model, BIOME-BGC, into the form of the MAPSS biogeography model. Specific enhancements are: 1) the 2-layer canopy biophysics of Dolman (1993); 2) the unique MAPSS-based hydrology, which incorporates canopy evaporation, snow dynamics, infiltration and saturated and unsaturated percolation with ‘fast’ flow and base flow and a ‘tunable aquifer’ capacity, a metaphor of D’Arcy’s Law; and, 3) a unique MAPSS-based stomatal conductance algorithm, which simultaneously incorporates vapor pressure and soil water potential constraints, based on physiological information and many other improvements. Over small domains the PFTs can be parameterized as individual species to investigate fundamental vs. potential niche theory; while, at more coarse scales the PFTs can be rendered as more general functional groups. Since all of the model processes are intrinsically leaf to plot scale (physiology to PFT competition), it essentially has no ‘intrinsic’ scale and can be implemented on a grid of any size, taking on the characteristics defined by the homogeneous climate of each grid cell. Currently, the model is implemented on the VEMAP 1/2 degree, daily grid over the conterminous U.S. Although both the thermal and water-limited ecotones are dynamic, following climate variability, the PFT distributions remain fixed. Thus, the model is currently being fitted with a ‘reproduction niche’ to allow full dynamic operation as a Dynamic General Vegetation Model (DGVM). While global simulations
NASA Astrophysics Data System (ADS)
Sheridan, J. A.; Bloom, D. M.; Solomon, P. M.
1995-03-01
We have built a system capable of measuring the step response of III-V electronic devices on the picosecond time scale, with no alteration in device design or epitaxy. To switch on the device under test (DUT), we have designed and fabricated a new type of photoconductor, the recessed-ohmic photoconductor, which swings 0.45 V with a 2-ps rise time and maintains constant output voltage for 100 ps. This switch is monolithically integrated with the DUT. To measure the output current of the DUT, we have built a Ti:sapphire-laser-based pump-probe direct electro-optic sampling system that has a minimum detectable voltage of 70 mu V / \\radical Hz \\end-radical and a measurement bandwidth of 750 GHz. The overall system, comprised of the recessed ohmic photoconductor and the electro-optic sampling system, can be used to measure the step response of III-V electronic devices on the picosecond time scale.
A class of large time step Godunov schemes for hyperbolic conservation laws and applications
NASA Astrophysics Data System (ADS)
Qian, ZhanSen; Lee, Chun-Hian
2011-08-01
A large time step (LTS) Godunov scheme firstly proposed by LeVeque is further developed in the present work and applied to Euler equations. Based on the analysis of the computational performances of LeVeque's linear approximation on wave interactions, a multi-wave approximation on rarefaction fan is proposed to avoid the occurrences of rarefaction shocks in computations. The developed LTS scheme is validated using 1-D test cases, manifesting high resolution for discontinuities and the capability of maintaining computational stability when large CFL numbers are imposed. The scheme is then extended to multidimensional problems using dimensional splitting technique; the treatment of boundary condition for this multidimensional LTS scheme is also proposed. As for demonstration problems, inviscid flows over NACA0012 airfoil and ONERA M6 wing with given swept angle are simulated using the developed LTS scheme. The numerical results reveal the high resolution nature of the scheme, where the shock can be captured within 1-2 grid points. The resolution of the scheme would improve gradually along with the increasing of CFL number under an upper bound where the solution becomes severely oscillating across the shock. Computational efficiency comparisons show that the developed scheme is capable of reducing the computational time effectively with increasing the time step (CFL number).
A multistage time-stepping scheme for the Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Swanson, R. C.; Turkel, E.
1985-01-01
A class of explicit multistage time-stepping schemes is used to construct an algorithm for solving the compressible Navier-Stokes equations. Flexibility in treating arbitrary geometries is obtained with a finite-volume formulation. Numerical efficiency is achieved by employing techniques for accelerating convergence to steady state. Computer processing is enhanced through vectorization of the algorithm. The scheme is evaluated by solving laminar and turbulent flows over a flat plate and an NACA 0012 airfoil. Numerical results are compared with theoretical solutions or other numerical solutions and/or experimental data.
Imaginary Time Step Method to Solve the Dirac Equation with Nonlocal Potential
Zhang Ying; Liang Haozhao; Meng Jie
2009-08-26
The imaginary time step (ITS) method is applied to solve the Dirac equation with nonlocal potentials in coordinate space. Taking the nucleus {sup 12}C as an example, even with nonlocal potentials, the direct ITS evolution for the Dirac equation still meets the disaster of the Dirac sea. However, following the recipe in our former investigation, the disaster can be avoided by the ITS evolution for the corresponding Schroedinger-like equation without localization, which gives the convergent results exactly the same with those obtained iteratively by the shooting method with localized effective potentials.
Omelyan, Igor E-mail: omelyan@icmp.lviv.ua; Kovalenko, Andriy
2013-12-28
steered by effective solvation forces allows huge outer time steps up to tens of picoseconds without affecting the equilibrium and conformational properties, and thus provides a 100- to 500-fold effective speedup in comparison to conventional MD with explicit solvent. With the statistical-mechanical 3D-RISM-KH account for effective solvation forces, the method provides efficient sampling of biomolecular processes with slow and/or rare solvation events such as conformational transitions of hydrated alanine dipeptide with the mean life times ranging from 30 ps up to 10 ns for “flip-flop” conformations, and is particularly beneficial for biomolecular systems with exchange and localization of solvent and ions, ligand binding, and molecular recognition.
A simple method for improving the time-stepping accuracy in atmosphere and ocean models
NASA Astrophysics Data System (ADS)
Williams, P. D.
2012-12-01
In contemporary numerical simulations of the atmosphere and ocean, evidence suggests that time-stepping errors may be a significant component of total model error, on both weather and climate time-scales. This presentation will review the available evidence, and will then suggest a simple but effective method for substantially improving the time-stepping numerics at no extra computational expense. A common time-stepping method in atmosphere and ocean models is the leapfrog scheme combined with the Robert-Asselin (RA) filter. This method is used in the following models (and many more): ECHAM, MAECHAM, MM5, CAM, MESO-NH, HIRLAM, KMCM, LIMA, SPEEDY, IGCM, PUMA, COSMO, FSU-GSM, FSU-NRSM, NCEP-GFS, NCEP-RSM, NSEAM, NOGAPS, RAMS, and CCSR/NIES-AGCM. Although the RA filter controls the time-splitting instability, it also introduces non-physical damping and reduces the accuracy. This presentation proposes a simple modification to the RA filter, which has become known as the RAW filter (Williams 2009, 2011). When used in conjunction with the leapfrog scheme, the RAW filter eliminates the non-physical damping and increases the amplitude accuracy by two orders, yielding third-order accuracy. (The phase accuracy remains second-order.) The RAW filter can easily be incorporated into existing models, typically via the insertion of just a single line of code. Better simulations are obtained at no extra computational expense. Results will be shown from recent implementations of the RAW filter in various models, including SPEEDY and COSMO. For example, in SPEEDY, the skill of weather forecasts is found to be significantly improved. In particular, in tropical surface pressure predictions, five-day forecasts made using the RAW filter have approximately the same skill as four-day forecasts made using the RA filter (Amezcua, Kalnay & Williams 2011). These improvements are encouraging for the use of the RAW filter in other atmosphere and ocean models. References PD Williams (2009) A
NASA Astrophysics Data System (ADS)
Ho, C. Y.; Leung, R. C. K.; Zhou, K.; Lam, G. C. Y.; Jiang, Z.
2011-09-01
One-step direct aeroacoustic simulation (DAS) has received attention from aerospace and mechanical high-pressure fluid-moving system manufacturers for quite some time. They aim to simulate the unsteady flow and acoustic field in the duct simultaneously in order to investigate the aeroacoustic generation mechanisms. Because of the large length and energy scale disparities between the acoustic far field and the aerodynamic near field, highly accurate and high-resolution simulation scheme is required. This involves the use of high order compact finite difference and time advancement schemes in simulation. However, in this situation, large buffer zones are always needed to suppress the spurious numerical waves emanating from computational boundaries. This further increases the computational resources to yield accurate results. On the other hand, for such problem as supersonic jet noise, the numerical scheme should be able to resolve both strong shock waves and weak acoustic waves simultaneously. Usually numerical aeroa-coustic scheme that is good for low Mach number flow is not able to give satisfactory simulation results for shock wave. Therefore, the aeroacoustic research community has been looking for a more efficient one-step DAS scheme that has the comparable accuracy to the finite-difference approach with smaller buffer regions, yet is able to give accurate solutions from subsonic to supersonic flows. The conservation element and solution element (CE/SE) scheme is one of the possible schemes satisfying the above requirements. This paper aims to report the development of a CE/SE scheme for one-step DAS and illustrate its robustness and effectiveness with two selected benchmark problems.
Detection of Zika virus by SYBR green one-step real-time RT-PCR.
Xu, Ming-Yue; Liu, Si-Qing; Deng, Cheng-Lin; Zhang, Qiu-Yan; Zhang, Bo
2016-10-01
The ongoing Zika virus (ZIKV) outbreak has rapidly spread to new areas of Americas, which were the first transmissions outside its traditional endemic areas in Africa and Asia. Due to the link with newborn defects and neurological disorder, numerous infected cases throughout the world and various mosquito vectors, the virus has been considered to be an international public health emergency. In the present study, we developed a SYBR Green based one-step real-time RT-PCR assay for rapid detection of ZIKV. Our results revealed that the real-time assay is highly specific and sensitive in detection of ZIKV in cell samples. Importantly, the replication of ZIKV at different time points in infected cells could be rapidly monitored by the real-time RT-PCR assay. Specifically, the real-time RT-PCR showed acceptable performance in measurement of infectious ZIKV RNA. This assay could detect ZIKV at a titer as low as 1PFU/mL. The real-time RT-PCR assay could be a useful tool for further virology surveillance and diagnosis of ZIKV. PMID:27444120
Multiple ``time step'' Monte Carlo simulations: Application to charged systems with Ewald summation
NASA Astrophysics Data System (ADS)
Bernacki, Katarzyna; Hetényi, Balázs; Berne, B. J.
2004-07-01
Recently, we have proposed an efficient scheme for Monte Carlo simulations, the multiple "time step" Monte Carlo (MTS-MC) [J. Chem. Phys. 117, 8203 (2002)] based on the separation of the potential interactions into two additive parts. In this paper, the structural and thermodynamic properties of the simple point charge water model combined with the Ewald sum are compared for the MTS-MC real-/reciprocal-space split of the Ewald summation and the common Metropolis Monte Carlo method. We report a number of observables as a function of CPU time calculated using MC and MTS-MC. The correlation functions indicate that speedups on the order of 4.5-7.5 can be obtained for systems of 108-500 waters for n=10 splitting parameter.
The multiple time step r-RESPA procedure and polarizable potentials based on induced dipole moments
NASA Astrophysics Data System (ADS)
Masella, Michel
In the present study, we present an accelerating scheme based on the reversible multiple time step r-RESPA method to be used in molecular dynamics simulations with polarizable potentials based on induced dipole moments. Even if the induced dipoles are estimated with an iterative self-consistent procedure, this scheme significantly reduces the CPU time needed to perform a molecular dynamics simulation, up to a factor 2, as compared to the Car-Parrinello method where additional dynamical variables are introduced for the treatment of the induced dipoles. The tests show that stable and reliable molecular dynamics trajectories can be generated with that scheme, and that the physical properties derived from the trajectories are equivalent to those computed with the classical all atom iterative approach and the Car-Parrinello one.
An Efficient Time-Stepping Scheme for Ab Initio Molecular Dynamics Simulations
NASA Astrophysics Data System (ADS)
Tsuchida, Eiji
2016-08-01
In ab initio molecular dynamics simulations of real-world problems, the simple Verlet method is still widely used for integrating the equations of motion, while more efficient algorithms are routinely used in classical molecular dynamics. We show that if the Verlet method is used in conjunction with pre- and postprocessing, the accuracy of the time integration is significantly improved with only a small computational overhead. We also propose several extensions of the algorithm required for use in ab initio molecular dynamics. The validity of the processed Verlet method is demonstrated in several examples including ab initio molecular dynamics simulations of liquid water. The structural properties obtained from the processed Verlet method are found to be sufficiently accurate even for large time steps close to the stability limit. This approach results in a 2× performance gain over the standard Verlet method for a given accuracy. We also show how to generate a canonical ensemble within this approach.
NASA Astrophysics Data System (ADS)
Gupta, Shubhangi; Wohlmuth, Barbara; Helmig, Rainer
2016-05-01
We present an extrapolation-based semi-implicit multi-rate time stepping (MRT) scheme and a compound-fast MRT scheme for a naturally partitioned, multi-time-scale hydro-geomechanical hydrate reservoir model. We evaluate the performance of the two MRT methods compared to an iteratively coupled solution scheme and discuss their advantages and disadvantages. The performance of the two MRT methods is evaluated in terms of speed-up and accuracy by comparison to an iteratively coupled solution scheme. We observe that the extrapolation-based semi-implicit method gives a higher speed-up but is strongly dependent on the relative time scales of the latent (slow) and active (fast) components. On the other hand, the compound-fast method is more robust and less sensitive to the relative time scales, but gives lower speed up as compared to the semi-implicit method, especially when the relative time scales of the active and latent components are comparable.
NASA Technical Reports Server (NTRS)
Chang, Chau-Lyan; Venkatachari, Balaji Shankar; Cheng, Gary
2013-01-01
With the wide availability of affordable multiple-core parallel supercomputers, next generation numerical simulations of flow physics are being focused on unsteady computations for problems involving multiple time scales and multiple physics. These simulations require higher solution accuracy than most algorithms and computational fluid dynamics codes currently available. This paper focuses on the developmental effort for high-fidelity multi-dimensional, unstructured-mesh flow solvers using the space-time conservation element, solution element (CESE) framework. Two approaches have been investigated in this research in order to provide high-accuracy, cross-cutting numerical simulations for a variety of flow regimes: 1) time-accurate local time stepping and 2) highorder CESE method. The first approach utilizes consistent numerical formulations in the space-time flux integration to preserve temporal conservation across the cells with different marching time steps. Such approach relieves the stringent time step constraint associated with the smallest time step in the computational domain while preserving temporal accuracy for all the cells. For flows involving multiple scales, both numerical accuracy and efficiency can be significantly enhanced. The second approach extends the current CESE solver to higher-order accuracy. Unlike other existing explicit high-order methods for unstructured meshes, the CESE framework maintains a CFL condition of one for arbitrarily high-order formulations while retaining the same compact stencil as its second-order counterpart. For large-scale unsteady computations, this feature substantially enhances numerical efficiency. Numerical formulations and validations using benchmark problems are discussed in this paper along with realistic examples.
The two-step shape and timing of the last deglaciation in Antarctica
Jouzel, J.; Petit, J.R. |; Duclos, Y.
1995-04-01
The two-step character of the last deglaciation is well recognized in Western Europe, in Greenland and in the North Atlantic. For example, in Greenland, a gradual temperature decrease started at the Boelling (B) around 14.5 ky BP, spanned through the Alleroed (A) and was followed by the cold Younger Dryas (YD) event which terminated abruptly around 11.5 ky BP. Recent results suggest that this BA/YD sequence may have extended throughout all the Northern Hemisphere but the evidence of a late transition cooling is still poor for the Southern Hemisphere. Here we present a detailed isotopic record analyzed in a new ice core drilled at Dome B in East Antarctica that fully demonstrates the existence of an Antarctic cold reversal (ACR). These results suggest that the two-step shape of the last deglaciation has a worldwide character but they also point to noticeable interhemispheric differences. Thus. the coldest part of the ACR. which shows a temperature drop about three times weaker than that recorded during the YD in Greenland, may have preceded the YD. Antarctica did not experienced abrupt changes and the two warming periods started there before they started in Greenland. The links between Southern and Northern Hemisphere climates throughout this period are discussed in the light of additional information derived from the Antarctic dust record. 87 refs., 5 figs.
[Photodissociation of Acetylene and Acetone using Step-Scan Time-Resolved FTIR Emission Spectroscopy
NASA Technical Reports Server (NTRS)
McLaren, Ian A.; Wrobel, Jacek D.
1997-01-01
The photodissociation of acetylene and acetone was investigated as a function of added quenching gas pressures using step-scan time-resolved FTIR emission spectroscopy. Its main components consist of Bruker IFS88, step-scan Fourier Transform Infrared (FTIR) spectrometer coupled to a flow cell equipped with Welsh collection optics. Vibrationally excited C2H radicals were produced from the photodissociation of acetylene in the unfocused experiments. The infrared (IR) emission from these excited C2H radicals was investigated as a function of added argon pressure. Argon quenching rate constants for all C2H emission bands are of the order of 10(exp -13)cc/molecule.sec. Quenching of these radicals by acetylene is efficient, with a rate constant in the range of 10(exp -11) cc/molecule.sec. The relative intensity of the different C2H emission bands did not change with the increasing argon or acetylene pressure. However, the overall IR emission intensity decreased, for example, by more than 50% when the argon partial pressure was raised from 0.2 to 2 Torr at fixed precursor pressure of 160mTorr. These observations provide evidence for the formation of a metastable C2H2 species, which are collisionally quenched by argon or acetylene. Problems encountered in the course of the experimental work are also described.
Structural damage evolution assessment using the regularised time step integration method
NASA Astrophysics Data System (ADS)
Chen, Hua-Peng; Maung, Than Soe
2014-09-01
This paper presents an approach to identify both the location and severity evolution of damage in engineering structures directly from measured dynamic response data. A relationship between the change in structural parameters such as stiffness caused by structural damage development and the measured dynamic response data such as accelerations is proposed, on the basis of the governing equations of motion for the original and damaged structural systems. Structural damage parameters associated with time are properly chosen to reflect both the location and severity development over time of damage in a structure. Basic equations are provided to solve the chosen time-dependent damage parameters, which are constructed by using the Newmark time step integration method without requiring a modal analysis procedure. The Tikhonov regularisation method incorporating the L-curve criterion for determining the regularisation parameter is then employed to reduce the influence of measurement errors in dynamic response data and then to produce stable solutions for structural damage parameters. Results for two numerical examples with various simulated damage scenarios show that the proposed method can accurately identify the locations of structural damage and correctly assess the evolution of damage severity from information on vibration measurements with uncertainties.
NASA Astrophysics Data System (ADS)
Roland, Teboh; Mavroidis, Panayiotis; Shi, Chengyu; Papanikolaou, Nikos
2010-05-01
System latency introduces geometric errors in the course of real-time target tracking radiotherapy. This effect can be minimized, for example by the use of predictive filters, but cannot be completely avoided. In this work, we present a convolution technique that can incorporate the effect as part of the treatment planning process. The method can be applied independently or in conjunction with the predictive filters to compensate for residual latency effects. The implementation was performed on TrackBeam (Initia Ltd, Israel), a prototype real-time target tracking system assembled and evaluated at our Cancer Institute. For the experimental system settings examined, a Gaussian distribution attributable to the TrackBeam latency was derived with σ = 3.7 mm. The TrackBeam latency, expressed as an average response time, was deduced to be 172 ms. Phantom investigations were further performed to verify the convolution technique. In addition, patient studies involving 4DCT volumes of previously treated lung cancer patients were performed to incorporate the latency effect in the dose prediction step. This also enabled us to effectively quantify the dosimetric and radiobiological impact of the TrackBeam and other higher latency effects on the clinical outcome of a real-time target tracking delivery.
NASA Astrophysics Data System (ADS)
Murthi, A.; Menon, S.; Sednev, I.
2011-12-01
An inherent difficulty in the ability of global climate models to accurately simulate precipitation lies in the use of a large time step, Δt (usually 30 minutes), to solve the governing equations. Since microphysical processes are characterized by small time scales compared to Δt, finite difference approximations used to advance microphysics equations suffer from numerical instability and large time truncation errors. With this in mind, the sensitivity of precipitation simulated by the atmospheric component of CESM, namely the Community Atmosphere Model (CAM 5.1), to the microphysics time step (τ) is investigated. Model integrations are carried out for a period of five years with a spin up time of about six months for a horizontal resolution of 2.5 × 1.9 degrees and 30 levels in the vertical, with Δt = 1800 s. The control simulation with τ = 900 s is compared with one using τ = 300 s for accumulated precipitation and radi- ation budgets at the surface and top of the atmosphere (TOA), while keeping Δt fixed. Our choice of τ = 300 s is motivated by previous work on warm rain processes wherein it was shown that a value of τ around 300 s was necessary, but not sufficient, to ensure positive definiteness and numerical stability of the explicit time integration scheme used to integrate the microphysical equations. However, since the entire suite of microphysical processes are represented in our case, we suspect that this might impose additional restrictions on τ. The τ = 300 s case produces differences in large-scale accumulated rainfall from the τ = 900 s case by as large as 200 mm, over certain regions of the globe. The spatial patterns of total accumulated precipitation using τ = 300 s are in closer agreement with satellite observed precipitation, when compared to the τ = 900 s case. Differences are also seen in the radiation budget with the τ = 300 (900) s cases producing surpluses that range between 1-3 W/m2 at both the TOA and surface in the global
Driancourt, M A; Cox, P; Rubion, S; Harnois-Milon, G; Kemp, B; Soede, N M
2013-09-01
The aim of this study was to demonstrate successful breeding of sows with a single fixed-time insemination following ovulation induction by buserelin, a GnRH analogue. In a first step, the optimal dose of buserelin (6, 10, or 16 μg) injected at 77 hours after weaning was determined in weaned sows (N = 15, 11, and 12, respectively) using its ability to induce an LH surge of similar magnitude as in control sows (N = 15) and induce ovulation. In 29/38 treated sows (76%), ovulation was induced and synchronized between 32 and 44 hours after injection, and the proportion of females ovulating during this time window was similar between groups at 73%, 73%, and 83% (6, 10, or 16 μg, respectively). Interestingly, whereas ovulation of 100% multiparous sows was induced and synchronized in the 32 to 44 hours posttreatment time window, successful induction was achieved in a lower proportion of primiparous sows (50%, 50%, and 67% following 6, 10, or 16 μg, respectively), the dose effect being nonsignificant. The magnitude of the LH surge was similar between control and treated sows, irrespective of the buserelin dose injected. Neither ovulation rate nor the number of good embryos on Day 5 postovulation differed between groups. Interestingly, the frequency of follicular cysts at slaughter was significantly affected by treatment (P < 0.05), being minimal and maximal in sows treated with 10 or 6 μg buserelin, respectively. In a second step, 419 sows from commercial herds in Spain, Germany, and France were randomly allocated to a control or treated group. The control sows were inseminated twice 12 ± 4 hours apart once estrus was detected. Treated sows received 10 μg buserelin at 86 ± 3 hours after weaning and were inseminated once 30 to 33 hours later. Farrowing rate of treated sows (87%, 166/192) was similar to that of control sows (84.5%, 169/200). Litter size was also similar between treated and control sows (13.6 ± 3.8 vs. 13.7 ± 3.2). In multiparous sows, neither
Detection and Correction of Step Discontinuities in Kepler Flux Time Series
NASA Technical Reports Server (NTRS)
Kolodziejczak, J. J.; Morris, R. L.
2011-01-01
PDC 8.0 includes an implementation of a new algorithm to detect and correct step discontinuities appearing in roughly one of every 20 stellar light curves during a given quarter. The majority of such discontinuities are believed to result from high-energy particles (either cosmic or solar in origin) striking the photometer and causing permanent local changes (typically -0.5%) in quantum efficiency, though a partial exponential recovery is often observed [1]. Since these features, dubbed sudden pixel sensitivity dropouts (SPSDs), are uncorrelated across targets they cannot be properly accounted for by the current detrending algorithm. PDC detrending is based on the assumption that features in flux time series are due either to intrinsic stellar phenomena or to systematic errors and that systematics will exhibit measurable correlations across targets. SPSD events violate these assumptions and their successful removal not only rectifies the flux values of affected targets, but demonstrably improves the overall performance of PDC detrending [1].
Electric and hybrid electric vehicle study utilizing a time-stepping simulation
NASA Technical Reports Server (NTRS)
Schreiber, Jeffrey G.; Shaltens, Richard K.; Beremand, Donald G.
1992-01-01
The applicability of NASA's advanced power technologies to electric and hybrid vehicles was assessed using a time-stepping computer simulation to model electric and hybrid vehicles operating over the Federal Urban Driving Schedule (FUDS). Both the energy and power demands of the FUDS were taken into account and vehicle economy, range, and performance were addressed simultaneously. Results indicate that a hybrid electric vehicle (HEV) configured with a flywheel buffer energy storage device and a free-piston Stirling convertor fulfills the emissions, fuel economy, range, and performance requirements that would make it acceptable to the consumer. It is noted that an assessment to determine which of the candidate technologies are suited for the HEV application has yet to be made. A proper assessment should take into account the fuel economy and range, along with the driveability and total emissions produced.
Comparison of Fixed and Variable Time Step Trajectory Integration Methods for Cislunar Trajectories
NASA Technical Reports Server (NTRS)
Weeks, ichael W.; Thrasher, Stephen W.
2007-01-01
Due to the nonlinear nature of the Earth-Moon-Sun three-body problem and non-spherical gravity, CEV cislunar targeting algorithms will require many propagations in their search for a desired trajectory. For on-board targeting especially, the algorithm must have a simple, fast, and accurate propagator to calculate a trajectory with reasonable computation time, and still be robust enough to remain stable in the various flight regimes that the CEV will experience. This paper compares Cowell s method with a fourth-order Runge- Kutta integrator (RK4), Encke s method with a fourth-order Runge-Kutta- Nystr m integrator (RKN4), and a method known as Multi-Conic. Additionally, the study includes the Bond-Gottlieb 14-element method (BG14) and extends the investigation of Encke-Nystrom methods to integrators of higher order and with variable step size.
Nikzad, Nasim; Sahari, Mohammad A; Vanak, Zahra Piravi; Safafar, Hamed; Boland-nazar, Seyed A
2013-08-01
Weight, oil, fatty acids, tocopherol, polyphenol, and sterol properties of 5 olive cultivars (Zard, Fishomi, Ascolana, Amigdalolia, and Conservalia) during crude, lye treatment, washing, fermentation, and pasteurization steps were studied. Results showed: oil percent was higher and lower in Ascolana (crude step) and in Fishomi (pasteurization step), respectively; during processing steps, in all cultivars, oleic, palmitic, linoleic, and stearic acids were higher; the highest changes in saturated and unsaturated fatty acids were in fermentation step; the highest and the lowest ratios of ω3 / ω6 were in Ascolana (washing step) and in Zard (pasteurization step), respectively; the highest and the lowest tocopherol were in Amigdalolia and Fishomi, respectively, and major damage occurred in lye step; the highest and the lowest polyphenols were in Ascolana (crude step) and in Zard and Ascolana (pasteurization step), respectively; the major damage among cultivars occurred during lye step, in which the polyphenol reduced to 1/10 of first content; sterol did not undergo changes during steps. Reviewing of olive patents shows that many compositions of fruits such as oil quality, fatty acids, quantity and its fraction can be changed by alteration in cultivar and process. PMID:23688142
NASA Astrophysics Data System (ADS)
Tassotto, Michael
2001-08-01
Liquid surfaces are very abundant in nature. Despite the importance of the liquid interface in general, experimental molecular-level data was almost completely lacking prior to the last decade and a half. The intent of this work is to provide a means by which experimental data on the composition of liquid surfaces and the average orientation of their constituent molecules can be obtained in order to supplement data from molecular dynamics and related computational techniques. To this end, a unique time-of-flight (TOF) spectrometer, which constitutes the backbone of a new method to study liquid surfaces, was constructed and commissioned. The performance of the spectrometer is demonstrated in a number of exemplary TOF spectra obtained from liquid glycerol. Moving from mere qualitative to quantitative surface analysis necessitates the ability to relate physical quantities such as detection efficiencies, accurate signal intensities, and interaction cross-sections for all elements to one another. As a first step, the absolute detection efficiency of a channel electron multiplier, used as particle detector in the spectrometer, was measured for the noble gas ions He+, Ar+, and Xe +. The data obtained led to an empirically derived, general expression of the detection efficiency that is applicable to particles of any atomic number. The results also show that the threshold velocity, below which a multiplier does not respond to impinging ions, cannot be regarded as independent of the ion's atomic number as previously reported in the literature. The second step involved a comprehensive investigation of ion-atom interactions and spectral features that are crucial for the processing of experimental signal intensities for quantitative analysis. For this purpose, the binary collision code MARLOWE was used in extensive trajectory calculations simulating TOF spectra. The simulation results confirm the high surface sensitivity of the technique and reveal the strong dependence of the
The role of the time step and overshooting in the modelling of PMS evolution: The case of EK Cephei
NASA Astrophysics Data System (ADS)
Marques, J. P.; Fernandes, J.; Monteiro, M. J. P. F. G.
2004-07-01
EK Cephei (HD 206821) is a unique candidate to test predictions based on stellar evolutionary models. It is a double-lined detached eclipsing binary system with accurate absolute dimensions available and a precise determination of the metallicity. Most importantly for our work, its low mass (1.12 Msun) component appears to be in the pre-main sequence (PMS) phase. We have produced detailed evolutionary models of the binary EK Cep using the CESAM stellar evolution code (Morel \\cite{morel97}). A χ2-minimisation was performed to derive the most reliable set of modelling parameters (age, αA, αB and Yi). We have found that an evolutionary age of about 26.8 Myr fits both components in the same isochrone. The positions of EK Cep A and B in the HR diagram are consistent (within the observational uncertainties) with our results. Our revised calibration shows clearly that EK Cep A is in the beginning of the main sequence, while EK Cep B is indeed a PMS star. Such a combination allows for a precise age determination of the binary, and provides a strict test of the modelling. In particular we have found that the definition of the time step in calculating the PMS evolution is crucial to reproduce the observations. A discussion of the optimal time step for calculating PMS evolution is presented. The fitting to the radii of both components is a more difficult task; although we managed to do it for EK Cep B, EK Cep A has a lower radius than our best models. We further studied the effect of the inclusion of a moderate convective overshooting; the calibration of the binary is not significantly altered, but the effect of the inclusion of overshooting can be dramatic in the approach to the main sequence of stars with masses high enough to burn hydrogen through the CNO cycle on the main sequence.
ERIC Educational Resources Information Center
Svetcov, Eric
2005-01-01
This article provides a list of the essential steps to keeping a school's or district's network safe and sound. It describes how to establish a security architecture and approach that will continually evolve as the threat environment changes over time. The article discusses the methodology for implementing this approach and then discusses the…
Halleroed, Tomas Rylander, Thomas
2008-04-20
A stable hybridization of the finite-element method (FEM) and the finite-difference time-domain (FDTD) scheme for Maxwell's equations with electric and magnetic losses is presented for two-dimensional problems. The hybrid method combines the flexibility of the FEM with the efficiency of the FDTD scheme and it is based directly on Ampere's and Faraday's law. The electric and magnetic losses can be treated implicitly by the FEM on an unstructured mesh, which allows for local mesh refinement in order to resolve rapid variations in the material parameters and/or the electromagnetic field. It is also feasible to handle larger homogeneous regions with losses by the explicit FDTD scheme connected to an implicitly time-stepped and lossy FEM region. The hybrid method shows second-order convergence for smooth scatterers. The bistatic radar cross section (RCS) for a circular metal cylinder with a lossy coating converges to the analytical solution and an accuracy of 2% is achieved for about 20 points per wavelength. The monostatic RCS for an airfoil that features sharp corners yields a lower order of convergence and it is found to agree well with what can be expected for singular fields at the sharp corners. A careful convergence study with resolutions from 20 to 140 points per wavelength provides accurate extrapolated results for this non-trivial test case, which makes it possible to use as a reference problem for scattering codes that model both electric and magnetic losses.
Outward Bound to the Galaxies--One Step at a Time
ERIC Educational Resources Information Center
Ward, R. Bruce; Miller-Friedmann, Jaimie; Sienkiewicz, Frank; Antonucci, Paul
2012-01-01
Less than a century ago, astronomers began to unlock the cosmic distances within and beyond the Milky Way. Understanding the size and scale of the universe is a continuing, step-by-step process that began with the remarkably accurate measurement of the distance to the Moon made by early Greeks. In part, the authors have ITEAMS (Innovative…
ERIC Educational Resources Information Center
Berkeley, Viv
2008-01-01
Disability equality and improving life chances for disabled people are high on the U.K. Government's agenda at present. As someone who was once a practitioner and a manager in a busy further education environment, the author knows how difficult it is to find time to think about, and discuss with colleagues, the implications of policy. This article…
Herzog, W; Zatsiorsky, V; Prilutsky, B I; Leonard, T R
1994-06-01
Force-sharing among muscles during locomotion has been studied experimentally using 'representative' or 'average' step cycles. Mathematical approaches aimed at predicting individual muscle forces during locomotion are based on the assumption that force-sharing among muscles occurs in a consistent and unique way. In this study, we quantify normal variations in muscular force-time histories for step cycles executed at a given nominal speed, so that we can appreciate what it means to analyze 'representative' or 'average' step cycles and can evaluate whether these normal variations in muscular force-time histories are random or may be associated with variations in the kinematics of consecutive step cycles. Forces in gastrocnemius, soleus and plantaris muscles were measured for step cycles performed at a constant nominal speed in freely moving cats. Gastrocnemius forces were always larger than peak plantaris or soleus forces. Also, peak gastrocnemius forces typically occurred first after paw contact, followed by peak soleus and then peak plantaris forces. Furthermore, it was found that variations in muscular force-time histories were substantial and were systematically related to step-cycle durations. The results of this study suggest that findings based on 'representative' or 'average' step cycles for a given nominal speed of locomotion should be viewed cautiously and that variations in force-sharing among muscles are systematically related to variations in locomotor kinematics. PMID:7931035
Singh, N; Maniscalco, S J; Fisher, H F
1993-01-01
We introduce a novel transient-state kinetic approach which can resolve proton and product time courses into a series of individual steps that comprise the reaction path. We have applied this approach to the oxidative deamination reaction catalyzed by bovine liver glutamate dehydrogenase, measuring both the product (NADPH) and proton time courses at various pH values. The global treatment (over all pH values) resolves the very early portion of this reaction quantitatively and provides a continuous time course for each of the six protonic species. We propose the following mechanism: L-glutamate binds to an open conformation of the enzyme-NADP complex, forming salt bridges between its alpha- and gamma-carboxyl groups and the protonated forms of enzyme lysine residues 114 and 90, respectively. In this position, the alpha-H atom of the substrate is too far from the nicotinamide ring for hydride transfer to occur. In the next step, three events occur in a concerted manner: lysine 126 loses a proton and acquires a single water molecule; the active site cleft closes; bulk water is expelled; the substrate and coenzyme are forced closer together and remain in a nonaqueous environment during the ensuing chemical events, returning to an open conformation only in time to allow the product release steps to occur. Thus, substrate binding accomplishes a number of important tasks which are themselves an integral part of the catalytic mechanism. Combining the novel transient state approach developed here with steady-state kinetic information can produce a detailed mechanistic resolution of otherwise hidden steps. PMID:8093240
NASA Astrophysics Data System (ADS)
Jordan, D. M.; Hill, J. D.; Uman, M. A.; Dwyer, J. R.; Rassoul, H. K.
2012-12-01
Time-of-arrival (TOA) techniques were used to determine the three-dimensional locations and emission times of x-ray and dE/dt sources measured at ground level in association with dart-stepped leader steps in natural and rocket-and-wire triggered lightning discharges recorded during summer 2011 at Camp Blanding, FL. The measurement network consisted of ten flat plate dE/dt antennas approximately co-located with eight plastic and two Lanthanum Bromide scintillation detectors arrayed around the launching facility over an area of about 0.25 square kilometers. For two triggered lightning dart-stepped leaders, x-ray sources were emitted from locations separated by average distances of 22.7 m and 29 m, respectively, from the locations of the associated dE/dt pulse peaks. The x-ray sources occurred beneath the dE/dt sources in 88% of the cases. X-rays were emitted from 20 ns to 2.16 μs following the dE/dt pulse peaks, with average temporal separations of 150 ns and 290 ns, respectively, for the two triggered lightning events. For one natural lightning dart-stepped leader, x-ray sources were emitted an average total distance of 39.2 m from the associated dE/dt pulse peak, and occurred beneath the location of the dE/dt source in 86% of the cases. The x-rays were emitted from 10 ns to 1.76 μs following the dE/dt pulse peak with an average temporal separation of 280 ns. In each of the three events, the altitude displacement between the dE/dt and x-ray sources dominated the total separation, accounting for 90%, 63%, and 72%, respectively, of the total separation. X-ray sources were distributed randomly in the lateral directions about the lightning channel in each event. For the triggered lightning events, x-rays were located from 2.5-83.5 μs prior to the return stroke at altitudes ranging from 24-336 m. For the natural lightning event, x-rays were located from 40.4-222.3 μs prior to the return stroke at altitudes ranging from 99-394 m. Cumulatively, 67% of the located x
Code of Federal Regulations, 2014 CFR
2014-07-01
... that we must consider for allowing an employee to complete all aspects of relocation? 302-2.110 Section... Are there time factors that we must consider for allowing an employee to complete all aspects of... travel is approved and inform employees that they must complete all aspects of relocation within a...
Code of Federal Regulations, 2011 CFR
2011-07-01
... that we must consider for allowing an employee to complete all aspects of relocation? 302-2.110 Section... Are there time factors that we must consider for allowing an employee to complete all aspects of... travel is approved and inform employees that they must complete all aspects of relocation within a...
Code of Federal Regulations, 2010 CFR
2010-07-01
... that we must consider for allowing an employee to complete all aspects of relocation? 302-2.110 Section... Are there time factors that we must consider for allowing an employee to complete all aspects of... travel is approved and inform employees that they must complete all aspects of relocation within a...
Code of Federal Regulations, 2013 CFR
2013-07-01
... that we must consider for allowing an employee to complete all aspects of relocation? 302-2.110 Section... Are there time factors that we must consider for allowing an employee to complete all aspects of... travel is approved and inform employees that they must complete all aspects of relocation within a...
Code of Federal Regulations, 2012 CFR
2012-07-01
... that we must consider for allowing an employee to complete all aspects of relocation? 302-2.110 Section... Are there time factors that we must consider for allowing an employee to complete all aspects of... travel is approved and inform employees that they must complete all aspects of relocation within a...
Haworth, Kevin J.; Fowlkes, J. Brian; Carson, Paul L.; Kripfgans, Oliver D.
2009-01-01
A theoretical shot noise model to describe the output of a time-reversal experiment in a multiple-scattering medium is developed. This (non-wave equation based) model describes the following process. An arbitrary waveform is transmitted through a high-order multiple-scattering environment and recorded. The recorded signal is arbitrarily windowed and then time-reversed. The processed signal is retransmitted into the environment and the resulting signal recorded. The temporal and spatial signal and noise of this process is predicted statistically. It is found that the time when the noise is largest depends on the arbitrary windowing and this noise peak can occur at times outside the main lobe. To determine further trends, a common set of parameters is applied to the general result. It is seen that as the duration of the input function increases, the signal-to-noise ratio (SNR) decreases (independent of signal bandwidth). It is also seen that longer persisting impulse responses result in increased main lobe amplitudes and SNR. Assumptions underpinning the generalized shot noise model are compared to an experimental realization of a multiple-scattering medium (a time-reversal chaotic cavity). Results from the model are compared to random number numerical simulation. PMID:19425655
Kutscher, Marika; Rosenberger, Manuel; Schmauss, Bernhard; Meinel, Lorenz; Lorenz, Udo; Ohlsen, Knut; Hellmann, Ralf; Germershaus, Oliver
2016-07-01
In this study, sensor surface functionalization allowing the repetitive use of a sensing device was evaluated for antibody-based detection of living bacteria using an optical planar Bragg grating sensor. To achieve regenerable immobilization of bacteria specific antibodies, the heterobifunctional cross-linker N-succinimidyl 3-(2-pyridyldithio) propionate (SPDP) was linked to an aminosilanized sensor surface and subsequently reduced to expose sulfhydryl groups enabling the covalent conjugation of SPDP-activated antibodies via disulfide bonds. The immobilization of a capture antibody specific for Staphylococcus aureus on the sensor surface as well as specific binding of S. aureus could be monitored, highlighting the applicability of optical sensors for the specific detection of large biological structures. Reusability of bacteria saturated sensors was successfully demonstrated by cleaving the antibody along with bound bacteria through reduction of disulfide bonds and subsequent re-functionalization with activated antibody, resulting in comparable sensitivity towards S. aureus. PMID:26486822
NASA Astrophysics Data System (ADS)
Taylor, N. E.; Williamson, D. M.; Jardine, A. P.
2014-05-01
Equations of state can be used to predict the relationship between pressure, volume and temperature. However, in shock physics, they are usually only constrained by experimental observations of pressure and volume. Direct observation of temperature in a shock is therefore valuable in constraining equations of state. Bloomquist and Sheffield (1980, 1981) and Rosenberg and Partom (1984) have attempted such observations in poly(methyl methacrylate) (PMMA). However, their results disagree strongly above 2GPa shock pressure. Here we present an improved fabrication technique, to examine this outstanding issue. We make use of the fact that the electrical resistivity of most metals is a known function of both pressure and temperature. If the change in resistance of a thin metal thermistor gauge is measured during a shock experiment of known pressure, the temperature can be calculated directly. The time response is limited by the time taken for the gauge to reach thermal equilibrium with the medium in which it is embedded. Gold gauges of thickness up to 200 nm have been produced by thermal evaporation, and fully embedded in PMMA. These reach thermal equilibrium with the host material in under 1 us, allowing temperature measurement within the duration of a plate impact experiment.
NASA Astrophysics Data System (ADS)
Hirthe, Eugenia M.; Graf, Thomas
2012-12-01
The automatic non-iterative second-order time-stepping scheme based on the temporal truncation error proposed by Kavetski et al. [Kavetski D, Binning P, Sloan SW. Non-iterative time-stepping schemes with adaptive truncation error control for the solution of Richards equation. Water Resour Res 2002;38(10):1211, http://dx.doi.org/10.1029/2001WR000720.] is implemented into the code of the HydroGeoSphere model. This time-stepping scheme is applied for the first time to the low-Rayleigh-number thermal Elder problem of free convection in porous media [van Reeuwijk M, Mathias SA, Simmons CT, Ward JD. Insights from a pseudospectral approach to the Elder problem. Water Resour Res 2009;45:W04416, http://dx.doi.org/10.1029/2008WR007421.], and to the solutal [Shikaze SG, Sudicky EA, Schwartz FW. Density-dependent solute transport in discretely-fractured geological media: is prediction possible? J Contam Hydrol 1998;34:273-91] problem of free convection in fractured-porous media. Numerical simulations demonstrate that the proposed scheme efficiently limits the temporal truncation error to a user-defined tolerance by controlling the time-step size. The non-iterative second-order time-stepping scheme can be applied to (i) thermal and solutal variable-density flow problems, (ii) linear and non-linear density functions, and (iii) problems including porous and fractured-porous media.
Development of real-time diagnostics and feedback algorithms for JET in view of the next step
NASA Astrophysics Data System (ADS)
Murari, A.; Joffrin, E.; Felton, R.; Mazon, D.; Zabeo, L.; Albanese, R.; Arena, P.; Ambrosino, G.; Ariola, M.; Barana, O.; Bruno, M.; Laborde, L.; Moreau, D.; Piccolo, F.; Sartori, F.; Crisanti, F.; de la Luna, E.; Sanchez, J.; Contributors, EFDA-JET
2005-03-01
Real-time control of many plasma parameters will be an essential aspect in the development of reliable high performance operation of next step tokamaks. The main prerequisites for any feedback scheme are the precise real-time determination of the quantities to be controlled, requiring top quality and highly reliable diagnostics, and the availability of robust control algorithms. A new set of real-time diagnostics was recently implemented on JET to prove the feasibility of determining, with high accuracy and time resolution, the most important plasma quantities. Some of the signals now routinely provided in real time at JET are: (i) the internal inductance and the main confinement quantities obtained by calculating the Shafranov integrals from the pick-up coils with 2 ms time resolution; (ii) the electron temperature profile, from electron cylotron emission every 10 ms; (iii) the ion temperature and plasma toroidal velocity profiles, from charge exchange recombination spectroscopy, provided every 50 ms; and (iv) the safety factor profile, derived from the inversion of the polarimetric line integrals every 2 ms. With regard to feedback algorithms, new model-based controllers were developed to allow a more robust control of several plasma parameters. With these new tools, several real-time schemes were implemented, among which the most significant is the simultaneous control of the safety factor and the plasma pressure profiles using the additional heating systems (LH, NBI, ICRH) as actuators. The control strategy adopted in this case consists of a multi-variable model-based technique, which was implemented as a truncated singular value decomposition of an integral operator. This approach is considered essential for systems like tokamak machines, characterized by a strong mutual dependence of the various parameters and the distributed nature of the quantities, the plasma profiles, to be controlled. First encouraging results were also obtained using non
The Influence of Time Spent in Outdoor Play on Daily and Aerobic Step Count in Costa Rican Children
ERIC Educational Resources Information Center
Morera Castro, Maria del Rocio
2011-01-01
The purpose of this study is to examine the influence of time spent in outdoor play (i.e., on weekday and weekend days) on daily (i.e., average step count) and aerobic step count (i.e., average moderate to vigorous physical activity [MVPA] during the weekdays and weekend days) in fifth grade Costa Rican children. It was hypothesized that: (a)…
Tremblay, Jean Christophe; Carrington, Tucker Jr.
2004-12-15
If the Hamiltonian is time dependent it is common to solve the time-dependent Schroedinger equation by dividing the propagation interval into slices and using an (e.g., split operator, Chebyshev, Lanczos) approximate matrix exponential within each slice. We show that a preconditioned adaptive step size Runge-Kutta method can be much more efficient. For a chirped laser pulse designed to favor the dissociation of HF the preconditioned adaptive step size Runge-Kutta method is about an order of magnitude more efficient than the time sliced method.
Potenza, M. A. C.; Sanvito, T.; Argentiere, S.; Cella, C.; Paroli, B.; Lenardi, C.; Milani, P.
2015-01-01
The behavior of nanoparticles in biological systems is determined by their dimensions, size distribution, shape, surface chemistry, density, drug loading and stability; the characterization of these parameters in realistic conditions and the possibility to follow their evolution in vitro and in vivo are, in most of the cases, far from the capabilities of the standard characterization technologies. Optical techniques such as dynamic light scattering (DLS) are, in principle, well suited for in line characterization of nanoparticle, however their fail in characterizing the evolution of nanoparticle in solution where change in particle dimension and density is present. Here we present an in-line optical technique based on single particle extinction and scattering (SPES) overcoming the limitations typical of DLS and allowing for the efficient characterization of nanoparticle polydispersity, index of refraction and degradation dynamics in solution. Using SPES, we characterized the evolution of PLGA nanoparticles with different structures and drug payloads in solution and we compared the results with DLS. Our results suggest that SPES could be used as a process analytical technology for pharmaceutical nanoparticle production. PMID:26667064
NASA Astrophysics Data System (ADS)
Potenza, M. A. C.; Sanvito, T.; Argentiere, S.; Cella, C.; Paroli, B.; Lenardi, C.; Milani, P.
2015-12-01
The behavior of nanoparticles in biological systems is determined by their dimensions, size distribution, shape, surface chemistry, density, drug loading and stability; the characterization of these parameters in realistic conditions and the possibility to follow their evolution in vitro and in vivo are, in most of the cases, far from the capabilities of the standard characterization technologies. Optical techniques such as dynamic light scattering (DLS) are, in principle, well suited for in line characterization of nanoparticle, however their fail in characterizing the evolution of nanoparticle in solution where change in particle dimension and density is present. Here we present an in-line optical technique based on single particle extinction and scattering (SPES) overcoming the limitations typical of DLS and allowing for the efficient characterization of nanoparticle polydispersity, index of refraction and degradation dynamics in solution. Using SPES, we characterized the evolution of PLGA nanoparticles with different structures and drug payloads in solution and we compared the results with DLS. Our results suggest that SPES could be used as a process analytical technology for pharmaceutical nanoparticle production.
ERIC Educational Resources Information Center
Mucientes, A. E.; de la Pena, M. A.
2009-01-01
The concentration-time integrals method has been used to solve kinetic equations of parallel-consecutive first-order reactions with a reversible step. This method involves the determination of the area under the curve for the concentration of a given species against time. Computer techniques are used to integrate experimental curves and the method…
A Step Response Based Mixed-Signal BIST Approach for Continuous-time Linear Circuits
NASA Technical Reports Server (NTRS)
Walker, Alvernon; Lala, P. K.
2001-01-01
A new Mixed-Signal Built-in self-test approach that is based upon the step response of a reconfigurable (or multifunction) analog block is presented in this paper. The technique requires the overlapping step response of the Circuit Under Test (CUT) for two circuit configurations. Each configuration can be realized by changing the topology of the CUT or by sampling two CUT nodes with differing step responses. The technique can effectively detect both soft and hard faults and does not require an analog-to-digital converter (ADC) and/or digital-to-analog converter(DAC). It also does not require any precision voltage sources or comparators. This approach does not require any additional analog circuits to realize the test signal generator and sample circuits. The paper is concluded with the application of the proposed approach to a circuit found in the work of Epstein et al and two ITC 97 analog benchmark circuits.
NASA Astrophysics Data System (ADS)
Alerskans, Emy; Kaas, Eigil
2016-04-01
In semi-Lagrangian models used for climate and NWP the trajectories are normally/often determined kinematically. Here we propose a new method for calculating trajectories in a more dynamically consistent way by pre-integrating the governing equations in a pseudo-Lagrangian manner using a short time step. Only non-advective adiabatic terms are included in this calculation, i.e., the Coriolis and pressure gradient force plus gravity in the momentum equations, and the divergence term in the continuity equation. This integration is performed with a forward-backward time step. Optionally, the tendencies are filtered with a local space filter, which reduces the phase speed of short wave gravity and sound waves. The filter relaxes the time step limitation related to high frequency oscillations without compromising locality of the solution. The filter can be considered as an alternative to less local or global semi-implicit solvers. Once trajectories are estimated over a complete long advective time step the full set of governing equations is stepped forward using these trajectories in combination with a flux form semi-Lagrangian formulation of the equations. The methodology is designed to improve consistency and scalability on massively parallel systems, although here it has only been verified that the technique produces realistic results in a shallow water model and a 2D model based on the full Euler equations.
Rushford, Michael C.
2002-01-01
An optical monitoring instrument monitors etch depth and etch rate for controlling a wet-etching process. The instrument provides means for viewing through the back side of a thick optic onto a nearly index-matched interface. Optical baffling and the application of a photoresist mask minimize spurious reflections to allow for monitoring with extremely weak signals. A Wollaston prism enables linear translation for phase stepping.
Corrales, Louis R.; Devanathan, Ram
2006-09-01
Non-equilibrium molecular dynamics simulation trajectories must in principle conserve energy along the entire path. Processes exist in high-energy primary knock-on atom cascades that can affect the energy conservation, specifically during the ballistic phase where collisions bring atoms into very close proximities. The solution, in general, is to reduce the time step size of the simulation. This work explores the effects of variable time step algorithms and the effects of specifying a maximum displacement. The period of the ballistic phase can be well characterized by methods developed in this work to monitor the kinetic energy dissipation during a high-energy cascade.
ERIC Educational Resources Information Center
Thompson, Ron; Robinson, Denise
2008-01-01
The unprecedented degree of attention given to the learning and skills sector in England by successive New Labour governments has led to a significant increase in what is expected of the teaching workforce. To help meet these expectations, a "step change" in the quality of initial teacher training for the sector is promised, alongside provisions…
Multi-Step Ahead Predictions for Critical Levels in Physiological Time Series.
ElMoaqet, Hisham; Tilbury, Dawn M; Ramachandran, Satya Krishna
2016-07-01
Standard modeling and evaluation methods have been classically used in analyzing engineering dynamical systems where the fundamental problem is to minimize the (mean) error between the real and predicted systems. Although these methods have been applied to multi-step ahead predictions of physiological signals, it is often more important to predict clinically relevant events than just to match these signals. Adverse clinical events, which occur after a physiological signal breaches a clinically defined critical threshold, are a popular class of such events. This paper presents a framework for multi-step ahead predictions of critical levels of abnormality in physiological signals. First, a performance metric is presented for evaluating multi-step ahead predictions. Then, this metric is used to identify personalized models optimized with respect to predictions of critical levels of abnormality. To address the paucity of adverse events, weighted support vector machines and cost-sensitive learning are used to optimize the proposed framework with respect to statistical metrics that can take into account the relative rarity of such events. PMID:27244754
Energy Science and Technology Software Center (ESTSC)
2014-06-01
ARKode is part of a software family called SUNDIALS: SUite of Nonlinear and Differential/ALgebraic equation Solvers [1]. The ARKode solver library provides an adaptive-step time integration package for stiff, nonstiff and multi-rate systems of ordinary differential equations (ODEs) using Runge Kutta methods [2].
Operational flood control of a low-lying delta system using large time step Model Predictive Control
NASA Astrophysics Data System (ADS)
Tian, Xin; van Overloop, Peter-Jules; Negenborn, Rudy R.; van de Giesen, Nick
2015-01-01
The safety of low-lying deltas is threatened not only by riverine flooding but by storm-induced coastal flooding as well. For the purpose of flood control, these deltas are mostly protected in a man-made environment, where dikes, dams and other adjustable infrastructures, such as gates, barriers and pumps are widely constructed. Instead of always reinforcing and heightening these structures, it is worth considering making the most of the existing infrastructure to reduce the damage and manage the delta in an operational and overall way. In this study, an advanced real-time control approach, Model Predictive Control, is proposed to operate these structures in the Dutch delta system (the Rhine-Meuse delta). The application covers non-linearity in the dynamic behavior of the water system and the structures. To deal with the non-linearity, a linearization scheme is applied which directly uses the gate height instead of the structure flow as the control variable. Given the fact that MPC needs to compute control actions in real-time, we address issues regarding computational time. A new large time step scheme is proposed in order to save computation time, in which different control variables can have different control time steps. Simulation experiments demonstrate that Model Predictive Control with the large time step setting is able to control a delta system better and much more efficiently than the conventional operational schemes.
Hue, Kien Duong Thi; Tuan, Trung Vu; Thi, Hanh Tien Nguyen; Bich, Chau Tran Nguyen; Anh, Huy Huynh Le; Wills, Bridget A; Simmons, Cameron P
2011-11-01
Dengue is mosquito-borne virus infection that annually causes ~50 million clinically apparent cases worldwide. An internally controlled one-step real-time multiplex RT-PCR assay was developed for detection and quantitation of DENV RNA in plasma sample by using specific primers and fluorogenic TaqMan probes. All primers and probes targeted sequences near the 3' end of the NS5 gene. The method comprised two multiplex assays and was validated for sensitivity, specificity, linearity, reproducibility and precision. An internal control template was spiked into each clinical specimen to provide quality assurance for each experimental step. The assay allowed for detection of between 0.5 and 3 infectious particles per mL, is rapid and has been operationally characterized in 287 Vietnamese dengue patients from two therapeutic intervention trials at the Hospital for Tropical Diseases, Ho Chi Minh City, Viet Nam. PMID:21843553
Meyer, Chad D.; Balsara, Dinshaw S.; Aslam, Tariq D.
2014-01-15
Parabolic partial differential equations appear in several physical problems, including problems that have a dominant hyperbolic part coupled to a sub-dominant parabolic component. Explicit methods for their solution are easy to implement but have very restrictive time step constraints. Implicit solution methods can be unconditionally stable but have the disadvantage of being computationally costly or difficult to implement. Super-time-stepping methods for treating parabolic terms in mixed type partial differential equations occupy an intermediate position. In such methods each superstep takes “s” explicit Runge–Kutta-like time-steps to advance the parabolic terms by a time-step that is s{sup 2} times larger than a single explicit time-step. The expanded stability is usually obtained by mapping the short recursion relation of the explicit Runge–Kutta scheme to the recursion relation of some well-known, stable polynomial. Prior work has built temporally first- and second-order accurate super-time-stepping methods around the recursion relation associated with Chebyshev polynomials. Since their stability is based on the boundedness of the Chebyshev polynomials, these methods have been called RKC1 and RKC2. In this work we build temporally first- and second-order accurate super-time-stepping methods around the recursion relation associated with Legendre polynomials. We call these methods RKL1 and RKL2. The RKL1 method is first-order accurate in time; the RKL2 method is second-order accurate in time. We verify that the newly-designed RKL1 and RKL2 schemes have a very desirable monotonicity preserving property for one-dimensional problems – a solution that is monotone at the beginning of a time step retains that property at the end of that time step. It is shown that RKL1 and RKL2 methods are stable for all values of the diffusion coefficient up to the maximum value. We call this a convex monotonicity preserving property and show by examples that it is very
NASA Astrophysics Data System (ADS)
Chen, Chun-Jung; Chang, Allen Y.; Tsai, Chang-Lung; Lee, Chih-Jen; Chou, Li-Ping; Shin, Tien-Hao
2012-04-01
A modified Waveform Relaxation algorithm with transmission line calculation ability is proposed to perform large-scale circuit simulation for MOSFET circuits with lossy coupled transmission lines. The adopted full time-domain transmission line calculation algorithm, based on the Method of Characteristic, has been equipped with a time step control scheme to improve the calculation efficiency. All proposed methods have been implemented in a simulation program to simulate several circuits. The simulation results well justify the success of proposed methods.
NASA Technical Reports Server (NTRS)
Rosenfeld, Moshe
1990-01-01
The development, validation and application of a fractional step solution method of the time-dependent incompressible Navier-Stokes equations in generalized coordinate systems are discussed. A solution method that combines a finite-volume discretization with a novel choice of the dependent variables and a fractional step splitting to obtain accurate solutions in arbitrary geometries was previously developed for fixed-grids. In the present research effort, this solution method is extended to include more general situations, including cases with moving grids. The numerical techniques are enhanced to gain efficiency and generality.
A dissociative fluorescence enhancement technique for one-step time-resolved immunoassays
Mukkala, Veli-Matti; Hakala, Harri H. O.; Mäkinen, Pauliina H.; Suonpää, Mikko U.; Hemmilä, Ilkka A.
2010-01-01
The limitation of current dissociative fluorescence enhancement techniques is that the lanthanide chelate structures used as molecular probes are not stable enough in one-step assays with high concentrations of complexones or metal ions in the reaction mixture since these substances interfere with lanthanide chelate conjugated to the detector molecule. Lanthanide chelates of diethylenetriaminepentaacetic acid (DTPA) are extremely stable, and we used EuDTPA derivatives conjugated to antibodies as tracers in one-step immunoassays containing high concentrations of complexones or metal ions. Enhancement solutions based on different β-diketones were developed and tested for their fluorescence-enhancing capability in immunoassays with EuDTPA-labelled antibodies. Characteristics tested were fluorescence intensity, analytical sensitivity, kinetics of complex formation and signal stability. Formation of fluorescent complexes is fast (5 min) in the presented enhancement solution with EuDTPA probes withstanding strong complexones (ethylenediaminetetra acetate (EDTA) up to 100 mM) or metal ions (up to 200 μM) in the reaction mixture, the signal is intensive, stable for 4 h and the analytical sensitivity with Eu is 40 fmol/L, Tb 130 fmol/L, Sm 2.1 pmol/L and Dy 8.5 pmol/L. With the improved fluorescence enhancement technique, EDTA and citrate plasma samples as well as samples containing relatively high concentrations of metal ions can be analysed using a one-step immunoassay format also at elevated temperatures. It facilitates four-plexing, is based on one chelate structure for detector molecule labelling and is suitable for immunoassays due to the wide dynamic range and the analytical sensitivity. Figure PMID:21161513
Kasa, Sawinee; Faksri, Kiatichai; Kaewkes, Wanlop; Lulitanond, Viraphong; Namwat, Wises
2015-01-01
Mycobacterium tuberculosis (M. tb) is a causative agent of tuberculosis, a worldwide public health problem. In recent years, the incidence of human mycobacterial infection due to species other than M. tb has increased. However, the lack of specific, rapid, and inexpensive methods for identification of mycobacterial species remains a pressing problem. A diagnostic test was developed for mycobacterial strain differentiation utilizing a double-step multiplex real time PCR together with melting curve analysis for identifying and distinguishing among M. tb, M. bovis BCG, other members of M. tb. complex, M. avium, and non-tuberculosis mycobacteria. The assay was tested using 167 clinical sputum samples in comparison with acid-fast staining and culturing. Using only the first step (step A) the assay achieved sensitivity and specificity of 81% and 95%, respectively. The detection limit was equivalent to 50 genome copies. PMID:26513906
Kuroda, Yoshihiro; Nisky, Ilana; Uranishi, Yuki; Imura, Masataka; Okamura, Allison M; Oshiro, Osamu
2013-01-01
We present a novel algorithm for real-time detection of the onset of surface electromyography signal in step-tracking wrist movements. The method identifies abrupt increase of the quasi-tension signal calculated from sEMG resulting from the step-by-step recruitment of activated motor units. We assessed the performance of our proposed algorithm using both simulated and real sEMG signals, and compared with two existing detection methods. Evaluation with simulated sEMG showed that the detection accuracy of our method is robust to different signal-to-noise ratios, and that it outperforms the existing methods in terms of bias when the noise is large (low SNR). Evaluation with real sEMG analysis also indicated better detection performance compared to existing methods. PMID:24110123
Ultra-fast consensus of discrete-time multi-agent systems with multi-step predictive output feedback
NASA Astrophysics Data System (ADS)
Zhang, Wenle; Liu, Jianchang
2016-04-01
This article addresses the ultra-fast consensus problem of high-order discrete-time multi-agent systems based on a unified consensus framework. A novel multi-step predictive output mechanism is proposed under a directed communication topology containing a spanning tree. By predicting the outputs of a network several steps ahead and adding this information into the consensus protocol, it is shown that the asymptotic convergence factor is improved by a power of q + 1 compared to the routine consensus. The difficult problem of selecting the optimal control gain is solved well by introducing a variable called convergence step. In addition, the ultra-fast formation achievement is studied on the basis of this new consensus protocol. Finally, the ultra-fast consensus with respect to a reference model and robust consensus is discussed. Some simulations are performed to illustrate the effectiveness of the theoretical results.
Finn, John M.
2015-03-01
Properties of integration schemes for solenoidal fields in three dimensions are studied, with a focus on integrating magnetic field lines in a plasma using adaptive time stepping. It is shown that implicit midpoint (IM) and a scheme we call three-dimensional leapfrog (LF) can do a good job (in the sense of preserving KAM tori) of integrating fields that are reversible, or (for LF) have a 'special divergence-free' property. We review the notion of a self-adjoint scheme, showing that such schemes are at least second order accurate and can always be formed by composing an arbitrary scheme with its adjoint. We also review the concept of reversibility, showing that a reversible but not exactly volume-preserving scheme can lead to a fractal invariant measure in a chaotic region, although this property may not often be observable. We also show numerical results indicating that the IM and LF schemes can fail to preserve KAM tori when the reversibility property (and the SDF property for LF) of the field is broken. We discuss extensions to measure preserving flows, the integration of magnetic field lines in a plasma and the integration of rays for several plasma waves. The main new result of this paper relates to non-uniform time stepping for volume-preserving flows. We investigate two potential schemes, both based on the general method of Ref. [11], in which the flow is integrated in split time steps, each Hamiltonian in two dimensions. The first scheme is an extension of the method of extended phase space, a well-proven method of symplectic integration with non-uniform time steps. This method is found not to work, and an explanation is given. The second method investigated is a method based on transformation to canonical variables for the two split-step Hamiltonian systems. This method, which is related to the method of non-canonical generating functions of Ref. [35], appears to work very well.
NASA Technical Reports Server (NTRS)
Janus, J. Mark; Whitfield, David L.
1990-01-01
Improvements are presented of a computer algorithm developed for the time-accurate flow analysis of rotating machines. The flow model is a finite volume method utilizing a high-resolution approximate Riemann solver for interface flux definitions. The numerical scheme is a block LU implicit iterative-refinement method which possesses apparent unconditional stability. Multiblock composite gridding is used to orderly partition the field into a specified arrangement of blocks exhibiting varying degrees of similarity. Block-block relative motion is achieved using local grid distortion to reduce grid skewness and accommodate arbitrary time step selection. A general high-order numerical scheme is applied to satisfy the geometric conservation law. An even-blade-count counterrotating unducted fan configuration is chosen for a computational study comparing solutions resulting from altering parameters such as time step size and iteration count. The solutions are compared with measured data.
NASA Astrophysics Data System (ADS)
Li, H.; Zhang, Z.; Chen, X.
2012-12-01
It is widely accepted that they are oversampled in spatial grid spacing and temporal time step in the high speed medium if uniform grids are used for the numerical simulation. This oversampled grid spacing and time step will lower the efficiency of the calculation, especially high velocity contrast exists. Based on the collocated-grid finite-difference method (FDM), we present an algorithm of spatial discontinuous grid, with localized grid blocks and locally varying time steps, which will increase the efficiency of simulation of seismic wave propagation and earthquake strong ground motion. According to the velocity structure, we discretize the model into discontinuous grid blocks, and the time step of each block is determined according to the local stability. The key problem of the discontinuous grid method is the connection between grid blocks with different grid spacing. We use a transitional area overlapped by both of the finer and the coarser grids to deal with the problem. In the transitional area, the values of finer ghost points are obtained by interpolation from the coarser grid in space and time domain, while the values of coarser ghost points are obtained by downsampling from the finer grid. How to deal with coarser ghost points can influent the stability of long time simulation. After testing different downsampling methods and finally we choose the Gaussian filtering. Basically, 4th order Rung-Kutta scheme will be used for the time integral for our numerical method. For our discontinuous grid FDM, discontinuous time steps for the coarser and the finer grids will be used to increase the simulation efficiency. Numerical tests indicate that our method can provide a stable solution even for the long time simulation without any additional filtration for grid spacing ratio n=2. And for larger grid spacing ratio, Gaussian filtration could be used to preserve the stability. With the collocated-grid FDM, which is flexible and accurate in implementation of free
ERIC Educational Resources Information Center
Nothard, Michaela; Irvine, Susan; Theobald, Maryanne; Staton, Sally; Pattinson, Cassandra; Thorpe, Karen
2015-01-01
Daytime sleep is a significant part of the daily routine for children attending early childhood education and care (ECEC) services in Australia and many other countries. The practice of sleep-time can account for a substantial portion of the day in ECEC and often involves a mandated sleep/rest period for all children, including older…
Belevich, Nikolai P; Bertsova, Yulia V; Verkhovskaya, Marina L; Baykov, Alexander A; Bogachev, Alexander V
2016-02-01
Bacterial Na(+)-translocating NADH:quinone oxidoreductase (Na(+)-NQR) uses a unique set of prosthetic redox groups-two covalently bound FMN residues, a [2Fe-2S] cluster, FAD, riboflavin and a Cys4[Fe] center-to catalyze electron transfer from NADH to ubiquinone in a reaction coupled with Na(+) translocation across the membrane. Here we used an ultra-fast microfluidic stopped-flow instrument to determine rate constants and the difference spectra for the six consecutive reaction steps of Vibrio harveyi Na(+)-NQR reduction by NADH. The instrument, with a dead time of 0.25 ms and optical path length of 1 cm allowed collection of visible spectra in 50-μs intervals. By comparing the spectra of reaction steps with the spectra of known redox transitions of individual enzyme cofactors, we were able to identify the chemical nature of most intermediates and the sequence of electron transfer events. A previously unknown spectral transition was detected and assigned to the Cys4[Fe] center reduction. Electron transfer from the [2Fe-2S] cluster to the Cys4[Fe] center and all subsequent steps were markedly accelerated when Na(+) concentration was increased from 20 μM to 25 mM, suggesting coupling of the former step with tight Na(+) binding to or occlusion by the enzyme. An alternating access mechanism was proposed to explain electron transfer between subunits NqrF and NqrC. According to the proposed mechanism, the Cys4[Fe] center is alternatively exposed to either side of the membrane, allowing the [2Fe-2S] cluster of NqrF and the FMN residue of NqrC to alternatively approach the Cys4[Fe] center from different sides of the membrane. PMID:26655930
NASA Astrophysics Data System (ADS)
Zasche, P.
2016-03-01
An easy step-by-step manual of PHOEBE is presented. It should serve as a starting point for the first time users of PHOEBE analyzing the eclipsing binary light curve. It is demonstrated on one particular detached system also with the downloadable data and the whole procedure is described easily till the final trustworthy fit is being reached.
Finite-difference modeling with variable grid-size and adaptive time-step in porous media
NASA Astrophysics Data System (ADS)
Liu, Xinxin; Yin, Xingyao; Wu, Guochen
2014-04-01
Forward modeling of elastic wave propagation in porous media has great importance for understanding and interpreting the influences of rock properties on characteristics of seismic wavefield. However, the finite-difference forward-modeling method is usually implemented with global spatial grid-size and time-step; it consumes large amounts of computational cost when small-scaled oil/gas-bearing structures or large velocity-contrast exist underground. To overcome this handicap, combined with variable grid-size and time-step, this paper developed a staggered-grid finite-difference scheme for elastic wave modeling in porous media. Variable finite-difference coefficients and wavefield interpolation were used to realize the transition of wave propagation between regions of different grid-size. The accuracy and efficiency of the algorithm were shown by numerical examples. The proposed method is advanced with low computational cost in elastic wave simulation for heterogeneous oil/gas reservoirs.
Gavrea, B. I.; Anitescu, M.; Potra, F. A.; Mathematics and Computer Science; Univ. of Pennsylvania; Univ. of Maryland
2008-01-01
In this work we present a framework for the convergence analysis in a measure differential inclusion sense of a class of time-stepping schemes for multibody dynamics with contacts, joints, and friction. This class of methods solves one linear complementarity problem per step and contains the semi-implicit Euler method, as well as trapezoidal-like methods for which second-order convergence was recently proved under certain conditions. By using the concept of a reduced friction cone, the analysis includes, for the first time, a convergence result for the case that includes joints. An unexpected intermediary result is that we are able to define a discrete velocity function of bounded variation, although the natural discrete velocity function produced by our algorithm may have unbounded variation.
Review of RT-LAB and Steps Involved for Implementation of a Simulink Model from MATLAB to REAL-TIME
NASA Astrophysics Data System (ADS)
Mikkili, Suresh; Panda, Anup Kumar
2013-11-01
In recent days, every researcher wants to develop his/her model in real-time. Simulation tools have been widely used for the design and improvement of electrical systems since the mid-twentieth century. The evolution of simulation tools has progressed in step with the evolution of computing technologies. Now a days, computing technologies have improved dramatically in performance and become widely available at a steadily decreasing cost. Consequently, simulation tools have also seen dramatic performance gains and steady cost decreases. Researchers and engineers now have access to affordable, high-performance simulation tools that were previously too cost prohibitive, except for the largest manufacturers and utilities.This article has introduced a specific class of digital simulator known as a real-time simulator by answering the questions "What is real-time simulation?" "Why is it needed" and "How it works". The latest trend in real-time simulation consists of exporting simulation models to FPGA. In this article, the steps involved for implementation of a model from MATLAB to REAL-TIME are provided in detail. The detailed real-time results are presented to support the feasibility of real-time digital simulator.
NASA Technical Reports Server (NTRS)
Ardalan, Sasan H.
1992-01-01
Two narrow-band radar systems are developed for high resolution target range estimation in inhomogeneous media. They are reformulations of two presently existing systems such that high resolution target range estimates may be achieved despite the use of narrow bandwidth radar pulses. A double sideband suppressed carrier radar technique originally derived in 1962, and later abandoned due to its inability to accurately measure target range in the presence of an interfering reflection, is rederived to incorporate the presence of an interfering reflection. The new derivation shows that the interfering reflection causes a period perturbation in the measured phase response. A high resolution spectral estimation technique is used to extract the period of this perturbation leading to accurate target range estimates independent of the signal-to-interference ratio. A non-linear optimal signal processing algorithm is derived for a frequency-stepped continuous wave radar system. The resolution enhancement offered by optimal signal processing of the data over the conventional Fourier Transform technique is clearly demonstrated using measured radar data. A method for modeling plane wave propagation in inhomogeneous media based on transmission line theory is derived and studied. Several simulation results including measurement of non-uniform electron plasma densities that develop near the heat tiles of a space re-entry vehicle are presented which verify the validity of the model.
Spectral (600-1050 nm) time exposures (99.6 μs) of a lightning stepped leader
NASA Astrophysics Data System (ADS)
Warner, Tom A.; Orville, Richard E.; Marshall, J. L.; Huggins, Kyle
2011-06-01
A cloud-to-ground lightning stepped leader has been recorded with a slitless spectrograph at a recording rate of 10,000 images per second at a distance of 0.6 km. Five sequential images of the leader spectra were recorded with an exposure (integration) time of 99.6 μs each over a spectral range from 600 to 1050 nm. These are the first stepped leader spectra covering the range 600 to 1050 nm. The last three spectra, obtained immediately before the return stroke, were analyzed at an altitude of between 108 and 122 m above a struck vehicle. The spectral emissions in the near infrared are dominated by neutral nitrogen and oxygen emissions, and Hα, with only a few emission lines from singly ionized nitrogen. A singly ionized nitrogen line at 661.1 nm is present in the first analyzed image, but not in the two subsequent images at the same height, which suggests a cooling of the channel. The emissions are integrated over a 99.6 μs exposure time and therefore show no evidence of stepping. The ensuing negative return stroke was detected by the National Lightning Detection Network and had an estimated peak current of -15.2 kA. One subsequent stroke was outside the field of view of the spectrograph. The flash occurred on 11 September 2009 near New Underwood, South Dakota, and the exact location of the first stroke is known because it struck a car traveling on Interstate 90. The stepped leader two-dimensional speed increased in the last four steps from 1.53 × 105 to 2.42 × 105 m/s with an average of 2.03 × 105 m/s.
Finite time step and spatial grid effects in δf simulation of warm plasmas
NASA Astrophysics Data System (ADS)
Sturdevant, Benjamin J.; Parker, Scott E.
2016-01-01
This paper introduces a technique for analyzing time integration methods used with the particle weight equations in δf method particle-in-cell (PIC) schemes. The analysis applies to the simulation of warm, uniform, periodic or infinite plasmas in the linear regime and considers the collective behavior similar to the analysis performed by Langdon for full-f PIC schemes [1,2]. We perform both a time integration analysis and spatial grid analysis for a kinetic ion, adiabatic electron model of ion acoustic waves. An implicit time integration scheme is studied in detail for δf simulations using our weight equation analysis and for full-f simulations using the method of Langdon. It is found that the δf method exhibits a CFL-like stability condition for low temperature ions, which is independent of the parameter characterizing the implicitness of the scheme. The accuracy of the real frequency and damping rate due to the discrete time and spatial schemes is also derived using a perturbative method. The theoretical analysis of numerical error presented here may be useful for the verification of simulations and for providing intuition for the design of new implicit time integration schemes for the δf method, as well as understanding differences between δf and full-f approaches to plasma simulation.
NASA Astrophysics Data System (ADS)
Guo, Peng; Cheng, Wenming; Wang, Yi
2015-11-01
This article considers the parallel machine scheduling problem with step-deteriorating jobs and sequence-dependent setup times. The objective is to minimize the total tardiness by determining the allocation and sequence of jobs on identical parallel machines. In this problem, the processing time of each job is a step function dependent upon its starting time. An individual extended time is penalized when the starting time of a job is later than a specific deterioration date. The possibility of deterioration of a job makes the parallel machine scheduling problem more challenging than ordinary ones. A mixed integer programming model for the optimal solution is derived. Due to its NP-hard nature, a hybrid discrete cuckoo search algorithm is proposed to solve this problem. In order to generate a good initial swarm, a modified Biskup-Hermann-Gupta (BHG) heuristic called MBHG is incorporated into the population initialization. Several discrete operators are proposed in the random walk of Lévy flights and the crossover search. Moreover, a local search procedure based on variable neighbourhood descent is integrated into the algorithm as a hybrid strategy in order to improve the quality of elite solutions. Computational experiments are executed on two sets of randomly generated test instances. The results show that the proposed hybrid algorithm can yield better solutions in comparison with the commercial solver CPLEX® with a one hour time limit, the discrete cuckoo search algorithm and the existing variable neighbourhood search algorithm.
Numerical simulation of diffusion MRI signals using an adaptive time-stepping method.
Li, Jing-Rebecca; Calhoun, Donna; Poupon, Cyril; Le Bihan, Denis
2014-01-20
The effect on the MRI signal of water diffusion in biological tissues in the presence of applied magnetic field gradient pulses can be modelled by a multiple compartment Bloch-Torrey partial differential equation. We present a method for the numerical solution of this equation by coupling a standard Cartesian spatial discretization with an adaptive time discretization. The time discretization is done using the explicit Runge-Kutta-Chebyshev method, which is more efficient than the forward Euler time discretization for diffusive-type problems. We use this approach to simulate the diffusion MRI signal from the extra-cylindrical compartment in a tissue model of the brain gray matter consisting of cylindrical and spherical cells and illustrate the effect of cell membrane permeability. PMID:24351275