STAB: A kinetic, three-dimensional, one-group digital computer program

International Nuclear Information System (INIS)

Curtis, A.R.; Tyror, J.G.; Wrigley, H.E.

1961-10-01

A computer program for solving the one-group, time dependent, three-dimensional diffusion equation together with auxiliary equations representing heat transfer, xenon production and control rod movements, is presented. The equations and the methods of solution are discussed. (author)

VNAP2: a computer program for computation of two-dimensional, time-dependent, compressible, turbulent flow

Energy Technology Data Exchange (ETDEWEB)

Cline, M.C.

1981-08-01

VNAP2 is a computer program for calculating turbulent (as well as laminar and inviscid), steady, and unsteady flow. VNAP2 solves the two-dimensional, time-dependent, compressible Navier-Stokes equations. The turbulence is modeled with either an algebraic mixing-length model, a one-equation model, or the Jones-Launder two-equation model. The geometry may be a single- or a dual-flowing stream. The interior grid points are computed using the unsplit MacCormack scheme. Two options to speed up the calculations for high Reynolds number flows are included. The boundary grid points are computed using a reference-plane-characteristic scheme with the viscous terms treated as source functions. An explicit artificial viscosity is included for shock computations. The fluid is assumed to be a perfect gas. The flow boundaries may be arbitrary curved solid walls, inflow/outflow boundaries, or free-jet envelopes. Typical problems that can be solved concern nozzles, inlets, jet-powered afterbodies, airfoils, and free-jet expansions. The accuracy and efficiency of the program are shown by calculations of several inviscid and turbulent flows. The program and its use are described completely, and six sample cases and a code listing are included.

Transpiration and film cooling boundary layer computer program. Volume 1: Numerical solutions of the turbulent boundary layer equations with equilibrium chemistry

Science.gov (United States)

Levine, J. N.

1971-01-01

A finite difference turbulent boundary layer computer program has been developed. The program is primarily oriented towards the calculation of boundary layer performance losses in rocket engines; however, the solution is general, and has much broader applicability. The effects of transpiration and film cooling as well as the effect of equilibrium chemical reactions (currently restricted to the H2-O2 system) can be calculated. The turbulent transport terms are evaluated using the phenomenological mixing length - eddy viscosity concept. The equations of motion are solved using the Crank-Nicolson implicit finite difference technique. The analysis and computer program have been checked out by solving a series of both laminar and turbulent test cases and comparing the results to data or other solutions. These comparisons have shown that the program is capable of producing very satisfactory results for a wide range of flows. Further refinements to the analysis and program, especially as applied to film cooling solutions, would be aided by the acquisition of a firm data base.

The Development Of A Computer Program For Thermohydraulic Analysis Of Kartini Reactor

International Nuclear Information System (INIS)

Ischaq, Ma'sum; Syarip; BS, Edi Trijono; Suyamto

1996-01-01

A computer programming that could be used to calculate the fuel temperature and the reactor core, has been developed. The inside fuel temperature was calculated by explicit method. First, the differential partial equations were arranged the temperature as a function of the fuel radius and lime. T = f(r.t), and then the equations were transformed into the finite difference equations that could be solved numerically by the computer. The convection heal transfer coefficient between the fuel and the coolant was calculated basically by the free convection phenomena that followed the equation Nu = f (Gr, Pr). By this computer programming, the fuel and the core temperature in a certain condition of the reactor power and the fluid could be predicted

Four-Cylinder Stirling-Engine Computer Program

Science.gov (United States)

Daniele, C. J.; Lorenzo, C. F.

1986-01-01

Computer program developed for simulating steady-state and transient performance of four-cylinder Stirling engine. In model, four cylinders interconnected by four working spaces. Each working space contains seven volumes: one for expansion space, heater, cooler, and compression space and three for regenerator. Thermal time constant for regenerator mass associated with each regenator gas volume. Former code generates results very quickly, since it has only 14 state variables with no energy equation. Current code then used to study various aspects of Stirling engine in much more detail. Program written in FORTRAN IV for use on IBM 370 computer.

Computations of Wall Distances Based on Differential Equations

Science.gov (United States)

Tucker, Paul G.; Rumsey, Chris L.; Spalart, Philippe R.; Bartels, Robert E.; Biedron, Robert T.

2004-01-01

The use of differential equations such as Eikonal, Hamilton-Jacobi and Poisson for the economical calculation of the nearest wall distance d, which is needed by some turbulence models, is explored. Modifications that could palliate some turbulence-modeling anomalies are also discussed. Economy is of especial value for deforming/adaptive grid problems. For these, ideally, d is repeatedly computed. It is shown that the Eikonal and Hamilton-Jacobi equations can be easy to implement when written in implicit (or iterated) advection and advection-diffusion equation analogous forms, respectively. These, like the Poisson Laplacian term, are commonly occurring in CFD solvers, allowing the re-use of efficient algorithms and code components. The use of the NASA CFL3D CFD program to solve the implicit Eikonal and Hamilton-Jacobi equations is explored. The re-formulated d equations are easy to implement, and are found to have robust convergence. For accurate Eikonal solutions, upwind metric differences are required. The Poisson approach is also found effective, and easiest to implement. Modified distances are not found to affect global outputs such as lift and drag significantly, at least in common situations such as airfoil flows.

Newnes circuit calculations pocket book with computer programs

CERN Document Server

Davies, Thomas J

2013-01-01

Newnes Circuit Calculations Pocket Book: With Computer Programs presents equations, examples, and problems in circuit calculations. The text includes 300 computer programs that help solve the problems presented. The book is comprised of 20 chapters that tackle different aspects of circuit calculation. The coverage of the text includes dc voltage, dc circuits, and network theorems. The book also covers oscillators, phasors, and transformers. The text will be useful to electrical engineers and other professionals whose work involves electronic circuitry.

The SIMRAND 1 computer program: Simulation of research and development projects

Science.gov (United States)

Miles, R. F., Jr.

1986-01-01

The SIMRAND I Computer Program (Version 5.0 x 0.3) written in Microsoft FORTRAN for the IBM PC microcomputer and its compatibles is described. The SIMRAND I Computer Program comprises eleven modules-a main routine and ten subroutines. Two additional files are used at compile time; one inserts the system or task equations into the source code, while the other inserts the dimension statements and common blocks. The SIMRAND I Computer Program can be run on most microcomputers or mainframe computers with only minor modifications to the computer code.

HAMOC: a computer program for fluid hammer analysis

International Nuclear Information System (INIS)

Johnson, H.G.

1975-12-01

A computer program has been developed for fluid hammer analysis of piping systems attached to a vessel which has undergone a known rapid pressure transient. The program is based on the characteristics method for solution of the partial differential equations of motion and continuity. Column separation logic is included for situations in which pressures fall to saturation values

BOOK REVIEW: Advanced Topics in Computational Partial Differential Equations: Numerical Methods and Diffpack Programming

Science.gov (United States)

Katsaounis, T. D.

2005-02-01

equations in Diffpack can be used to derive fully implicit solvers for systems. The proposed techniques are illustrated in terms of two applications, namely a system of PDEs modelling pipeflow and a two-phase porous media flow. Stochastic PDEs is the topic of chapter 7. The first part of the chapter is a simple introduction to stochastic PDEs; basic analytical properties are presented for simple models like transport phenomena and viscous drag forces. The second part considers the numerical solution of stochastic PDEs. Two basic techniques are presented, namely Monte Carlo and perturbation methods. The last part explains how to implement and incorporate these solvers into Diffpack. Chapter 8 describes how to operate Diffpack from Python scripts. The main goal here is to provide all the programming and technical details in order to glue the programming environment of Diffpack with visualization packages through Python and in general take advantage of the Python interfaces. Chapter 9 attempts to show how to use numerical experiments to measure the performance of various PDE solvers. The authors gathered a rather impressive list, a total of 14 PDE solvers. Solvers for problems like Poisson, Navier--Stokes, elasticity, two-phase flows and methods such as finite difference, finite element, multigrid, and gradient type methods are presented. The authors provide a series of numerical results combining various solvers with various methods in order to gain insight into their computational performance and efficiency. In Chapter 10 the authors consider a computationally challenging problem, namely the computation of the electrical activity of the human heart. After a brief introduction on the biology of the problem the authors present the mathematical models involved and a numerical method for solving them within the framework of Diffpack. Chapter 11 and 12 are closely related; actually they could have been combined in a single chapter. Chapter 11 introduces several mathematical

Symbolic computation of analytic approximate solutions for nonlinear fractional differential equations

Science.gov (United States)

Lin, Yezhi; Liu, Yinping; Li, Zhibin

2013-01-01

The Adomian decomposition method (ADM) is one of the most effective methods to construct analytic approximate solutions for nonlinear differential equations. In this paper, based on the new definition of the Adomian polynomials, Rach (2008) [22], the Adomian decomposition method and the Padé approximants technique, a new algorithm is proposed to construct analytic approximate solutions for nonlinear fractional differential equations with initial or boundary conditions. Furthermore, a MAPLE software package is developed to implement this new algorithm, which is user-friendly and efficient. One only needs to input the system equation, initial or boundary conditions and several necessary parameters, then our package will automatically deliver the analytic approximate solutions within a few seconds. Several different types of examples are given to illustrate the scope and demonstrate the validity of our package, especially for non-smooth initial value problems. Our package provides a helpful and easy-to-use tool in science and engineering simulations. Program summaryProgram title: ADMP Catalogue identifier: AENE_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AENE_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 12011 No. of bytes in distributed program, including test data, etc.: 575551 Distribution format: tar.gz Programming language: MAPLE R15. Computer: PCs. Operating system: Windows XP/7. RAM: 2 Gbytes Classification: 4.3. Nature of problem: Constructing analytic approximate solutions of nonlinear fractional differential equations with initial or boundary conditions. Non-smooth initial value problems can be solved by this program. Solution method: Based on the new definition of the Adomian polynomials [1], the Adomian decomposition method and the Pad

Program package for the computation of lenses and deflectors

International Nuclear Information System (INIS)

Lencova, B.; Wisselink, G.

1990-01-01

In this paper a set of computer programs for the design of electrostatic and magnetic electron lenses and for the design of multipoles for electron microscopy and lithography is described. The two-dimensional field computation is performed by the finite-element method. In order to meet the high demands on accuracy, the programs include the use of a variable step in the fine mesh made with an automeshing procedure, improved methods for coefficient evaluation, a fast solution procedure for the linear equations, and modified algorithms for computation of multipoles and electrostatic lenses. They allow for a fast and accurate computation of electron optical elements. For the input and modification of data, and for presentation of results, graphical menu driven programs written for personal computers are used. For the computation of electron optical properties axial fields are used. (orig.)

Development of C++ Application Program for Solving Quadratic Equation in Elementary School in Nigeria

Science.gov (United States)

Bandele, Samuel Oye; Adekunle, Adeyemi Suraju

2015-01-01

The study was conducted to design, develop and test a c++ application program CAP-QUAD for solving quadratic equation in elementary school in Nigeria. The package was developed in c++ using object-oriented programming language, other computer program that were also utilized during the development process is DevC++ compiler, it was used for…

Viscous wing theory development. Volume 2: GRUMWING computer program user's manual

Science.gov (United States)

Chow, R. R.; Ogilvie, P. L.

1986-01-01

This report is a user's manual which describes the operation of the computer program, GRUMWING. The program computes the viscous transonic flow over three-dimensional wings using a boundary layer type viscid-inviscid interaction approach. The inviscid solution is obtained by an approximate factorization (AFZ)method for the full potential equation. The boundary layer solution is based on integral entrainment methods.

Computer programs for optical dendrometer measurements of standing tree profiles

Science.gov (United States)

Jacob R. Beard; Thomas G. Matney; Emily B. Schultz

2015-01-01

Tree profile equations are effective volume predictors. Diameter data for building these equations are collected from felled trees using diameter tapes and calipers or from standing trees using optical dendrometers. Developing and implementing a profile function from the collected data is a tedious and error prone task. This study created a computer program, Profile...

Automation of reliability evaluation procedures through CARE - The computer-aided reliability estimation program.

Science.gov (United States)

Mathur, F. P.

1972-01-01

Description of an on-line interactive computer program called CARE (Computer-Aided Reliability Estimation) which can model self-repair and fault-tolerant organizations and perform certain other functions. Essentially CARE consists of a repository of mathematical equations defining the various basic redundancy schemes. These equations, under program control, are then interrelated to generate the desired mathematical model to fit the architecture of the system under evaluation. The mathematical model is then supplied with ground instances of its variables and is then evaluated to generate values for the reliability-theoretic functions applied to the model.

Parallel computation for solving the tridiagonal linear system of equations

International Nuclear Information System (INIS)

Ishiguro, Misako; Harada, Hiroo; Fujii, Minoru; Fujimura, Toichiro; Nakamura, Yasuhiro; Nanba, Katsumi.

1981-09-01

Recently, applications of parallel computation for scientific calculations have increased from the need of the high speed calculation of large scale programs. At the JAERI computing center, an array processor FACOM 230-75 APU has installed to study the applicability of parallel computation for nuclear codes. We made some numerical experiments by using the APU on the methods of solution of tridiagonal linear equation which is an important problem in scientific calculations. Referring to the recent papers with parallel methods, we investigate eight ones. These are Gauss elimination method, Parallel Gauss method, Accelerated parallel Gauss method, Jacobi method, Recursive doubling method, Cyclic reduction method, Chebyshev iteration method, and Conjugate gradient method. The computing time and accuracy were compared among the methods on the basis of the numerical experiments. As the result, it is found that the Cyclic reduction method is best both in computing time and accuracy and the Gauss elimination method is the second one. (author)

COLUMN2 - A computer program for simulating migration

International Nuclear Information System (INIS)

Nielsen, O.J.; Bo, P.; Carlsen, L.

1985-10-01

COLUMN2 is a 1D FORTRAN77 computer program designed for studies of the effects of various physicochemical processes on migration. It solves the solute transport equation and cant take into account dispersion, sorption, ion exchange, first and second order homogeneous chemical reactions. Spacial variations of input pulses and retention factors are possible. The method of solution is based on a finite difference discretion followed by the application of the method of characteristics and two separate grid systems. This report explains the mathematical and numerical methods used, describes the necessary input, contains a number of test examples, provides a listing of the program and explains how to acquire the program, adapt it to other computers and run it. This report serves as a manual for the program. (author)

Activity computer program for calculating ion irradiation activation

Science.gov (United States)

Palmer, Ben; Connolly, Brian; Read, Mark

2017-07-01

A computer program, Activity, was developed to predict the activity and gamma lines of materials irradiated with an ion beam. It uses the TENDL (Koning and Rochman, 2012) [1] proton reaction cross section database, the Stopping and Range of Ions in Matter (SRIM) (Biersack et al., 2010) code, a Nuclear Data Services (NDS) radioactive decay database (Sonzogni, 2006) [2] and an ENDF gamma decay database (Herman and Chadwick, 2006) [3]. An extended version of Bateman's equation is used to calculate the activity at time t, and this equation is solved analytically, with the option to also solve by numeric inverse Laplace Transform as a failsafe. The program outputs the expected activity and gamma lines of the activated material.

MINI-TRAC code: a driver program for assessment of constitutive equations of two-fluid model

International Nuclear Information System (INIS)

Akimoto, Hajime; Abe, Yutaka; Ohnuki, Akira; Murao, Yoshio

1991-05-01

MINI-TRAC code, a driver program for assessment of constitutive equations of two-fluid model, has been developed to perform assessment and improvement of constitutive equations of two-fluid model widely and efficiently. The MINI-TRAC code uses one-dimensional conservation equations for mass, momentum and energy based on the two-fluid model. The code can work on a personal computer because it can be operated with a core memory size less than 640 KB. The MINI-TRAC code includes constitutive equations of TRAC-PF1/MOD1 code, TRAC-BF1 code and RELAP5/MOD2 code. The code is modulated so that one can easily change constitutive equations to perform a test calculation. This report is a manual of the MINI-TRAC code. The basic equations, numerics, constitutive, equations included in the MINI-TRAC code will be described. The user's manual such as input description will be presented. The program structure and contents of main variables will also be mentioned in this report. (author)

Nonlinear Boltzmann equation for the homogeneous isotropic case: Minimal deterministic Matlab program

Science.gov (United States)

Asinari, Pietro

2010-10-01

The homogeneous isotropic Boltzmann equation (HIBE) is a fundamental dynamic model for many applications in thermodynamics, econophysics and sociodynamics. Despite recent hardware improvements, the solution of the Boltzmann equation remains extremely challenging from the computational point of view, in particular by deterministic methods (free of stochastic noise). This work aims to improve a deterministic direct method recently proposed [V.V. Aristov, Kluwer Academic Publishers, 2001] for solving the HIBE with a generic collisional kernel and, in particular, for taking care of the late dynamics of the relaxation towards the equilibrium. Essentially (a) the original problem is reformulated in terms of particle kinetic energy (exact particle number and energy conservation during microscopic collisions) and (b) the computation of the relaxation rates is improved by the DVM-like correction, where DVM stands for Discrete Velocity Model (ensuring that the macroscopic conservation laws are exactly satisfied). Both these corrections make possible to derive very accurate reference solutions for this test case. Moreover this work aims to distribute an open-source program (called HOMISBOLTZ), which can be redistributed and/or modified for dealing with different applications, under the terms of the GNU General Public License. The program has been purposely designed in order to be minimal, not only with regards to the reduced number of lines (less than 1000), but also with regards to the coding style (as simple as possible). Program summaryProgram title: HOMISBOLTZ Catalogue identifier: AEGN_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGN_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License No. of lines in distributed program, including test data, etc.: 23 340 No. of bytes in distributed program, including test data, etc.: 7 635 236 Distribution format: tar

Nonlinear streak computation using boundary region equations

Energy Technology Data Exchange (ETDEWEB)

Martin, J A; Martel, C, E-mail: juanangel.martin@upm.es, E-mail: carlos.martel@upm.es [Depto. de Fundamentos Matematicos, E.T.S.I Aeronauticos, Universidad Politecnica de Madrid, Plaza Cardenal Cisneros 3, 28040 Madrid (Spain)

2012-08-01

The boundary region equations (BREs) are applied for the simulation of the nonlinear evolution of a spanwise periodic array of streaks in a flat plate boundary layer. The well-known BRE formulation is obtained from the complete Navier-Stokes equations in the high Reynolds number limit, and provides the correct asymptotic description of three-dimensional boundary layer streaks. In this paper, a fast and robust streamwise marching scheme is introduced to perform their numerical integration. Typical streak computations present in the literature correspond to linear streaks or to small-amplitude nonlinear streaks computed using direct numerical simulation (DNS) or the nonlinear parabolized stability equations (PSEs). We use the BREs to numerically compute high-amplitude streaks, a method which requires much lower computational effort than DNS and does not have the consistency and convergence problems of the PSE. It is found that the flow configuration changes substantially as the amplitude of the streaks grows and the nonlinear effects come into play. The transversal motion (in the wall normal-streamwise plane) becomes more important and strongly distorts the streamwise velocity profiles, which end up being quite different from those of the linear case. We analyze in detail the resulting flow patterns for the nonlinearly saturated streaks and compare them with available experimental results. (paper)

Kmonodium, a Program for the Numerical Solution of the One-Dimensional Schrodinger Equation

Science.gov (United States)

Angeli, Celestino; Borini, Stefano; Cimiraglia, Renzo

2005-01-01

A very simple strategy for the solution of the Schrodinger equation of a particle moving in one dimension subjected to a generic potential is presented. This strategy is implemented in a computer program called Kmonodium, which is free and distributed under the General Public License (GPL).

TRUMP3-JR: a finite difference computer program for nonlinear heat conduction problems

International Nuclear Information System (INIS)

Ikushima, Takeshi

1984-02-01

Computer program TRUMP3-JR is a revised version of TRUMP3 which is a finite difference computer program used for the solution of multi-dimensional nonlinear heat conduction problems. Pre- and post-processings for input data generation and graphical representations of calculation results of TRUMP3 are avaiable in TRUMP3-JR. The calculation equations, program descriptions and user's instruction are presented. A sample problem is described to demonstrate the use of the program. (author)

An object-oriented programming paradigm for parallelization of computational fluid dynamics

International Nuclear Information System (INIS)

Ohta, Takashi.

1997-03-01

We propose an object-oriented programming paradigm for parallelization of scientific computing programs, and show that the approach can be a very useful strategy. Generally, parallelization of scientific programs tends to be complicated and unportable due to the specific requirements of each parallel computer or compiler. In this paper, we show that the object-oriented programming design, which separates the parallel processing parts from the solver of the applications, can achieve the large improvement in the maintenance of the codes, as well as the high portability. We design the program for the two-dimensional Euler equations according to the paradigm, and evaluate the parallel performance on IBM SP2. (author)

Computer programming and computer systems

CERN Document Server

Hassitt, Anthony

1966-01-01

Computer Programming and Computer Systems imparts a "reading knowledge? of computer systems.This book describes the aspects of machine-language programming, monitor systems, computer hardware, and advanced programming that every thorough programmer should be acquainted with. This text discusses the automatic electronic digital computers, symbolic language, Reverse Polish Notation, and Fortran into assembly language. The routine for reading blocked tapes, dimension statements in subroutines, general-purpose input routine, and efficient use of memory are also elaborated.This publication is inten

A new kind of science - or a not-so-new kind of computer program?

CERN Multimedia

Naiditch, D

2003-01-01

Critical discussion of Stephen Wolfram's theory that "scientists give up their unwieldy equations and instead employ the types of computational rules used in cellular automata (CA) and related computer programs" (2 pages).

HYDRA-II: A hydrothermal analysis computer code: Volume 1, Equations and numerics

International Nuclear Information System (INIS)

McCann, R.A.

1987-04-01

HYDRA-II is a hydrothermal computer code capable of three-dimensional analysis of coupled conduction, convection, and thermal radiation problems. This code is especially appropriate for simulating the steady-state performance of spent fuel storage systems. The code has been evaluated for this application for the US Department of Energy's Commercial Spent Fuel Management Program. HYDRA-II provides a finite difference solution in Cartesian coordinates to the equations governing the conservation of mass, momentum, and energy. A cylindrical coordinate system may also be used to enclose the Cartesian coordinate system. This exterior coordinate system is useful for modeling cylindrical cask bodies. The difference equations for conservation of momentum are enhanced by the incorporation of directional porosities and permeabilities that aid in modeling solid structures whose dimensions may be smaller than the computational mesh. The equation for conservation of energy permits of modeling of orthotropic physical properties and film resistances. Several automated procedures are available to model radiation transfer within enclosures and from fuel rod to fuel rod. The documentation of HYDRA-II is presented in three separate volumes. This volume, Volume I - Equations and Numerics, describes the basic differential equations, illustrates how the difference equations are formulated, and gives the solution procedures employed. Volume II - User's Manual contains code flow charts, discusses the code structure, provides detailed instructions for preparing an input file, and illustrates the operation of the code by means of a model problem. The final volume, Volume III - Verification/Validation Assessments, presents results of numerical simulations of single- and multiassembly storage systems and comparisons with experimental data. 4 refs

Symbolic computation of analytic approximate solutions for nonlinear differential equations with initial conditions

Science.gov (United States)

Lin, Yezhi; Liu, Yinping; Li, Zhibin

2012-01-01

The Adomian decomposition method (ADM) is one of the most effective methods for constructing analytic approximate solutions of nonlinear differential equations. In this paper, based on the new definition of the Adomian polynomials, and the two-step Adomian decomposition method (TSADM) combined with the Padé technique, a new algorithm is proposed to construct accurate analytic approximations of nonlinear differential equations with initial conditions. Furthermore, a MAPLE package is developed, which is user-friendly and efficient. One only needs to input a system, initial conditions and several necessary parameters, then our package will automatically deliver analytic approximate solutions within a few seconds. Several different types of examples are given to illustrate the validity of the package. Our program provides a helpful and easy-to-use tool in science and engineering to deal with initial value problems. Program summaryProgram title: NAPA Catalogue identifier: AEJZ_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEJZ_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 4060 No. of bytes in distributed program, including test data, etc.: 113 498 Distribution format: tar.gz Programming language: MAPLE R13 Computer: PC Operating system: Windows XP/7 RAM: 2 Gbytes Classification: 4.3 Nature of problem: Solve nonlinear differential equations with initial conditions. Solution method: Adomian decomposition method and Padé technique. Running time: Seconds at most in routine uses of the program. Special tasks may take up to some minutes.

Program system for computation of the terrestrial gamma-radiation field

International Nuclear Information System (INIS)

Kirkegaard, P.; Loevborg, L.

1979-02-01

A system of computer programs intended for solution of the plane one-dimensional photon transport equation in the case of two adjacent media is described, and user's guides for the programs are given. One medium represents a natural ground with uniformly distributed potassium, uranium, and thorium gamma-ray emitters. The other medium is usually air with no radioactive contaminants. The solution method is the double-P 1 approximation with logarithmic energy spacing. The complete data-processing system GB contains the transport-theory code GAMP1, the code GFX for computation of scalar flux and dose rate, and a number of auxiliary programs and data files. (author)

Development of computer program for safety of nuclear power plant against tsunami

International Nuclear Information System (INIS)

Jin, S. B.; Choi, K. R.; Lee, S. K.; Cho, Y. S.

2001-01-01

The main objective of this study is the development of a computer program to check the safety of nuclear power plants along the coastline of the Korean Peninsula. The computer program describes the propagation and associated run-up process of tsunamis by solving linear and nonlinear shallow-water equations with finite difference methods. The computer program has been applied to several ideal and simplified problems. Obtained numerical solutions are compared to existing and available solutions and measurements. A very good agreement between numerical solutions and existing measurement is observed. The computer program developed in this study can be to check the safety analysis of nuclear power plants against tsunamis. The program can also be used to study the propagation of tsunamis for a long distance, and associated run-up and run-down process along a shoreline. Furthermore, the computer program can be used to provide the proper design criteria of coastal facilities and structures

A computer program for fitting smooth surfaces to an aircraft configuration and other three dimensional geometries

Science.gov (United States)

Craidon, C. B.

1975-01-01

A computer program that uses a three-dimensional geometric technique for fitting a smooth surface to the component parts of an aircraft configuration is presented. The resulting surface equations are useful in performing various kinds of calculations in which a three-dimensional mathematical description is necessary. Programs options may be used to compute information for three-view and orthographic projections of the configuration as well as cross-section plots at any orientation through the configuration. The aircraft geometry input section of the program may be easily replaced with a surface point description in a different form so that the program could be of use for any three-dimensional surface equations.

Solving large sets of coupled equations iteratively by vector processing on the CYBER 205 computer

International Nuclear Information System (INIS)

Tolsma, L.D.

1985-01-01

The set of coupled linear second-order differential equations which has to be solved for the quantum-mechanical description of inelastic scattering of atomic and nuclear particles can be rewritten as an equivalent set of coupled integral equations. When some type of functions is used as piecewise analytic reference solutions, the integrals that arise in this set can be evaluated analytically. The set of integral equations can be solved iteratively. For the results mentioned an inward-outward iteration scheme has been applied. A concept of vectorization of coupled-channel Fortran programs, based on this integral method, is presented for the use on the Cyber 205 computer. It turns out that, for two heavy ion nuclear scattering test cases, this vector algorithm gives an overall speed-up of about a factor of 2 to 3 compared to a highly optimized scalar algorithm for a one vector pipeline computer

PRECONDITIONED CONJUGATE-GRADIENT 2 (PCG2), a computer program for solving ground-water flow equations

Science.gov (United States)

Hill, Mary C.

1990-01-01

This report documents PCG2 : a numerical code to be used with the U.S. Geological Survey modular three-dimensional, finite-difference, ground-water flow model . PCG2 uses the preconditioned conjugate-gradient method to solve the equations produced by the model for hydraulic head. Linear or nonlinear flow conditions may be simulated. PCG2 includes two reconditioning options : modified incomplete Cholesky preconditioning, which is efficient on scalar computers; and polynomial preconditioning, which requires less computer storage and, with modifications that depend on the computer used, is most efficient on vector computers . Convergence of the solver is determined using both head-change and residual criteria. Nonlinear problems are solved using Picard iterations. This documentation provides a description of the preconditioned conjugate gradient method and the two preconditioners, detailed instructions for linking PCG2 to the modular model, sample data inputs, a brief description of PCG2, and a FORTRAN listing.

Parallel Computation of the Jacobian Matrix for Nonlinear Equation Solvers Using MATLAB

Science.gov (United States)

Rose, Geoffrey K.; Nguyen, Duc T.; Newman, Brett A.

2017-01-01

Demonstrating speedup for parallel code on a multicore shared memory PC can be challenging in MATLAB due to underlying parallel operations that are often opaque to the user. This can limit potential for improvement of serial code even for the so-called embarrassingly parallel applications. One such application is the computation of the Jacobian matrix inherent to most nonlinear equation solvers. Computation of this matrix represents the primary bottleneck in nonlinear solver speed such that commercial finite element (FE) and multi-body-dynamic (MBD) codes attempt to minimize computations. A timing study using MATLAB's Parallel Computing Toolbox was performed for numerical computation of the Jacobian. Several approaches for implementing parallel code were investigated while only the single program multiple data (spmd) method using composite objects provided positive results. Parallel code speedup is demonstrated but the goal of linear speedup through the addition of processors was not achieved due to PC architecture.

Experimental quantum computing to solve systems of linear equations.

Science.gov (United States)

Cai, X-D; Weedbrook, C; Su, Z-E; Chen, M-C; Gu, Mile; Zhu, M-J; Li, Li; Liu, Nai-Le; Lu, Chao-Yang; Pan, Jian-Wei

2013-06-07

Solving linear systems of equations is ubiquitous in all areas of science and engineering. With rapidly growing data sets, such a task can be intractable for classical computers, as the best known classical algorithms require a time proportional to the number of variables N. A recently proposed quantum algorithm shows that quantum computers could solve linear systems in a time scale of order log(N), giving an exponential speedup over classical computers. Here we realize the simplest instance of this algorithm, solving 2×2 linear equations for various input vectors on a quantum computer. We use four quantum bits and four controlled logic gates to implement every subroutine required, demonstrating the working principle of this algorithm.

EVOLUT - a computer program for fast burnup evaluation

International Nuclear Information System (INIS)

Craciunescu, T.; Dobrin, R.; Stamatescu, L.; Alexa, A.

1999-01-01

EVOLUT is a computer program for burnup evaluation. The input data consist on the one hand of axial and radial gamma-scanning profiles (for the experimental evaluation of the number of nuclei of a fission product - the burnup monitor - at the end of irradiation) and on the other hand of the history of irradiation (the time length and values proportional to the neutron flux for each step of irradiation). Using the equation of evolution of the burnup monitor the flux values are iteratively adjusted, by a multiplier factor, until the calculated number of nuclei is equal to the experimental one. The flux values are used in the equation of evolution of the fissile and fertile nuclei to determine the fission number and consequently the burnup. EVOLUT was successfully used in the analysis of several hundreds of CANDU and TRIGA-type fuel rods. We appreciate that EVOLUT is a useful tool in the burnup evaluation based on gamma spectrometry measurements. EVOLUT can be used on an usual AT computer and in this case the results are obtained in a few minutes. It has an original and user-friendly graphical interface and it provides also output in script MATLAB files for graphical representation and further numerical analysis. The computer program needs simple data and it is valuable especially when a large number of burnup analyses are required quickly. (authors)

Computations of concentration of radon and its decay products against time. Computer program; Obliczanie koncentracji radonu i jego produktow rozpadu w funkcji czasu. Program komputerowy

Energy Technology Data Exchange (ETDEWEB)

Machaj, B. [Institute of Nuclear Chemistry and Technology, Warsaw (Poland)

1996-12-31

This research is aimed to develop a device for continuous monitoring of radon in the air, by measuring alpha activity of radon and its short lived decay products. The influence of alpha activity variation of radon and its daughters on the measured results is of importance and requires a knowledge of this variation with time. Employing the measurement of alpha radiation of radon and of its short lived decay products, require knowledge of radon concentration variation and its decay products against the time. A computer program in Turbo Pascal language was therefore developed performing the computations employing the known relations involved, the program being adapted for IBM PC computers. The presented program enables computation of activity of {sup 222}Rn and its daughter products: {sup 218}Po, {sup 214}Pb, {sup 214}Bi and {sup 214}Po every 1 min within the period of 0-255 min for any state of radiation equilibrium between the radon and its daughter products. The program permits also to compute alpha activity of {sup 222}Rn + {sup 218}Po + {sup 214}Po against time and the total alpha activity at selected interval of time. The results of computations are stored on the computer hard disk in ASCII format and are used a graphic program e.g. by DrawPerfect program to make diagrams. Equations employed for computation of the alpha activity of radon and its decay products as well as the description of program functions are given. (author). 2 refs, 4 figs.

Combustion of hydrogen-air jets in local chemical equilibrium: A guide to the CHARNAL computer program

Science.gov (United States)

Spalding, D. B.; Launder, B. E.; Morse, A. P.; Maples, G.

1974-01-01

A guide to a computer program, written in FORTRAN 4, for predicting the flow properties of turbulent mixing with combustion of a circular jet of hydrogen into a co-flowing stream of air is presented. The program, which is based upon the Imperial College group's PASSA series, solves differential equations for diffusion and dissipation of turbulent kinetic energy and also of the R.M.S. fluctuation of hydrogen concentration. The effective turbulent viscosity for use in the shear stress equation is computed. Chemical equilibrium is assumed throughout the flow.

Fluid/Structure Interaction Studies of Aircraft Using High Fidelity Equations on Parallel Computers

Science.gov (United States)

Guruswamy, Guru; VanDalsem, William (Technical Monitor)

1994-01-01

Abstract Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft such as High Speed Civil Transport (HSCT) and Advanced Subsonic Transport (AST) which can experience complex flow/structure interactions. HSCT can experience vortex induced aeroelastic oscillations whereas AST can experience transonic buffet associated structural oscillations. Both aircraft may experience a dip in the flutter speed at the transonic regime. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional super computers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overcome the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENSAERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations.

A stable computational scheme for stiff time-dependent constitutive equations

International Nuclear Information System (INIS)

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

1977-01-01

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

Expressions of manipulator kinematic equations via symbolic computation

International Nuclear Information System (INIS)

Sasaki, Shinobu

1993-09-01

While it is simple in principle to determine the position and orientation of the manipulator hand, its computational process has been regarded as extremely laborious since trigonometric functions must be calculated many times in operations of revolute or rotation. Due to development of a general class of kinematic algorithm based on iterative methods, however, we have come to a satisfactory settlement of this problem. In the present article, we consider to construct symbolic kinematic equations in an automatic fashion making use of the algorithm. To this end, recursive expressions are applied to a symbolic computation system REDUCE. As a concrete result, a complete kinematic model for a six-jointed arm having all kinematic attributes is provided. Together with work space analysis, the computer-aided generation of kinematic equations in symbolic form will serve to liberate us from their cumbersome derivations. (author)

COMPUTER PROGRAM FOR CALCULATION MICROCHANNEL HEAT EXCHANGERS FOR AIR CONDITIONING SYSTEMS

Directory of Open Access Journals (Sweden)

Olga V. Olshevska

2016-08-01

Full Text Available Creating a computer program to calculate microchannel air condensers to reduce design time and carrying out variant calculations. Software packages for thermophysical properties of the working substance and the coolant, the correlation equation for calculating heat transfer, aerodynamics and hydrodynamics, the thermodynamic equations for the irreversible losses and their minimization in the heat exchanger were used in the process of creating. Borland Delphi 7 is used for creating software package.

Abstract of programs for nuclear reactor calculation and kinetic equations solution

International Nuclear Information System (INIS)

Marakazov, A.A.

1977-01-01

The collection includes about 50 annotations of programmes,developed in the Kurchatov Atomic Energy Institute in 1971-1976. The programmes are intended for calculating the neutron flux, for solving systems of multigroup equations in P 3 approximation, for calculating the reactor cell, for analysing the system stability, breeding ratio etc. The programme annotations are compiled according to the following diagram: 1.Programme title. 2.Computer type. 3.Physical problem. 4.Solution method. 5.Calculation limitations. 6.Characteristic computer time. 7.Programme characteristic features. 8.Bound programmes. 9.Programme state. 10.Literature allusions in the programme. 11.Required memory resourses. 12.Programming language. 13.Operation system. 14.Names of authors and place of programme adjusting

TURBO: a computer program for two-dimensional incompressible fluid flow analysis using a two-equations turbulence model

International Nuclear Information System (INIS)

Botelho, D.A.; Moreira, M.L.

1991-06-01

The Reynolds turbulent transport equations for an incompressible fluid are integrated on a bi-dimensional staggered grid, for velocity and pressure, using the SIMPLER method. With the resulting algebraic relations it was developed the TURBO program, which final objectives are the thermal stratification and natural convection analysis of nuclear reactor pools. This program was tested in problems applications with analytic or experimental solutions previously known. (author)

Hypersonic Shock Wave Computations Using the Generalized Boltzmann Equation

Science.gov (United States)

Agarwal, Ramesh; Chen, Rui; Cheremisin, Felix G.

2006-11-01

Hypersonic shock structure in diatomic gases is computed by solving the Generalized Boltzmann Equation (GBE), where the internal and translational degrees of freedom are considered in the framework of quantum and classical mechanics respectively [1]. The computational framework available for the standard Boltzmann equation [2] is extended by including both the rotational and vibrational degrees of freedom in the GBE. There are two main difficulties encountered in computation of high Mach number flows of diatomic gases with internal degrees of freedom: (1) a large velocity domain is needed for accurate numerical description of the distribution function resulting in enormous computational effort in calculation of the collision integral, and (2) about 50 energy levels are needed for accurate representation of the rotational spectrum of the gas. Our methodology addresses these problems, and as a result the efficiency of calculations has increased by several orders of magnitude. The code has been validated by computing the shock structure in Nitrogen for Mach numbers up to 25 including the translational and rotational degrees of freedom. [1] Beylich, A., ``An Interlaced System for Nitrogen Gas,'' Proc. of CECAM Workshop, ENS de Lyon, France, 2000. [2] Cheremisin, F., ``Solution of the Boltzmann Kinetic Equation for High Speed Flows of a Rarefied Gas,'' Proc. of the 24th Int. Symp. on Rarefied Gas Dynamics, Bari, Italy, 2004.

Programming in biomolecular computation

DEFF Research Database (Denmark)

Hartmann, Lars Røeboe; Jones, Neil; Simonsen, Jakob Grue

2011-01-01

Our goal is to provide a top-down approach to biomolecular computation. In spite of widespread discussion about connections between biology and computation, one question seems notable by its absence: Where are the programs? We identify a number of common features in programming that seem...... conspicuously absent from the literature on biomolecular computing; to partially redress this absence, we introduce a model of computation that is evidently programmable, by programs reminiscent of low-level computer machine code; and at the same time biologically plausible: its functioning is defined...... by a single and relatively small set of chemical-like reaction rules. Further properties: the model is stored-program: programs are the same as data, so programs are not only executable, but are also compilable and interpretable. It is universal: all computable functions can be computed (in natural ways...

Determining Balıkesir’s Energy Potential Using a Regression Analysis Computer Program

Directory of Open Access Journals (Sweden)

Bedri Yüksel

2014-01-01

Full Text Available Solar power and wind energy are used concurrently during specific periods, while at other times only the more efficient is used, and hybrid systems make this possible. When establishing a hybrid system, the extent to which these two energy sources support each other needs to be taken into account. This paper is a study of the effects of wind speed, insolation levels, and the meteorological parameters of temperature and humidity on the energy potential in Balıkesir, in the Marmara region of Turkey. The relationship between the parameters was studied using a multiple linear regression method. Using a designed-for-purpose computer program, two different regression equations were derived, with wind speed being the dependent variable in the first and insolation levels in the second. The regression equations yielded accurate results. The computer program allowed for the rapid calculation of different acceptance rates. The results of the statistical analysis proved the reliability of the equations. An estimate of identified meteorological parameters and unknown parameters could be produced with a specified precision by using the regression analysis method. The regression equations also worked for the evaluation of energy potential.

RICE: a computer program for multicomponent chemically reactive flows at all speeds

International Nuclear Information System (INIS)

Rivard, W.C.; Farmer, O.A.; Butler, T.D.

1974-11-01

The fluid dynamics of chemically reactive mixtures are calculated at arbitrary flow speeds with the RICE program. The dynamics are governed by the two-dimensional, time-dependent Navier-Stokes equations together with the species transport equations and the mass-action rate equations for the chemical reactions. The mass and momentum equations for the mixture are solved implicitly by the ICE technique. The equations for total energy and species transport are solved explicitly while the chemical rate equations are solved implicitly with a time step that may be a submultiple of the hydrodynamic time step. Application is made to continuous wave HF chemical lasers to compute the supersonic mixing and chemical reactions that take place in the lasing cavity. (U.S.)

Geochemical modelling. Column 2: a computer program for simulation of migration

International Nuclear Information System (INIS)

Nielsen, O.J.; Carlsen, L.; Bo, P.

1985-01-01

COLUMN2 is a 1D FORTRAN77 computer program designed for studies of the effects of various physicochemical processes on migration. It solves the solute transport equation and can take into account dispersion, sorption, ion exchange, first and second order homogeneous chemical reactions. Spatial variations of input pulses and retention factors are possible. The method of solution is based on a finite difference discretion followed by the application of the method of characteristics and two separate grid systems. This report explains the mathematical and numerical methods used, describes the necessary input, contains a number of test examples, provides a listing of the program and explains how to acquire the program, adapt it to other computers and run it. This report serves as a manual for the program

Statistically and Computationally Efficient Estimating Equations for Large Spatial Datasets

KAUST Repository

Sun, Ying; Stein, Michael L.

2014-01-01

For Gaussian process models, likelihood based methods are often difficult to use with large irregularly spaced spatial datasets, because exact calculations of the likelihood for n observations require O(n3) operations and O(n2) memory. Various approximation methods have been developed to address the computational difficulties. In this paper, we propose new unbiased estimating equations based on score equation approximations that are both computationally and statistically efficient. We replace the inverse covariance matrix that appears in the score equations by a sparse matrix to approximate the quadratic forms, then set the resulting quadratic forms equal to their expected values to obtain unbiased estimating equations. The sparse matrix is constructed by a sparse inverse Cholesky approach to approximate the inverse covariance matrix. The statistical efficiency of the resulting unbiased estimating equations are evaluated both in theory and by numerical studies. Our methods are applied to nearly 90,000 satellite-based measurements of water vapor levels over a region in the Southeast Pacific Ocean.

Statistically and Computationally Efficient Estimating Equations for Large Spatial Datasets

KAUST Repository

Sun, Ying

2014-11-07

For Gaussian process models, likelihood based methods are often difficult to use with large irregularly spaced spatial datasets, because exact calculations of the likelihood for n observations require O(n3) operations and O(n2) memory. Various approximation methods have been developed to address the computational difficulties. In this paper, we propose new unbiased estimating equations based on score equation approximations that are both computationally and statistically efficient. We replace the inverse covariance matrix that appears in the score equations by a sparse matrix to approximate the quadratic forms, then set the resulting quadratic forms equal to their expected values to obtain unbiased estimating equations. The sparse matrix is constructed by a sparse inverse Cholesky approach to approximate the inverse covariance matrix. The statistical efficiency of the resulting unbiased estimating equations are evaluated both in theory and by numerical studies. Our methods are applied to nearly 90,000 satellite-based measurements of water vapor levels over a region in the Southeast Pacific Ocean.

Systematic Equation Formulation

DEFF Research Database (Denmark)

Lindberg, Erik

2007-01-01

A tutorial giving a very simple introduction to the set-up of the equations used as a model for an electrical/electronic circuit. The aim is to find a method which is as simple and general as possible with respect to implementation in a computer program. The “Modified Nodal Approach”, MNA, and th......, and the “Controlled Source Approach”, CSA, for systematic equation formulation are investigated. It is suggested that the kernel of the P Spice program based on MNA is reprogrammed....

Computational Challenge of Fractional Differential Equations and the Potential Solutions: A Survey

Directory of Open Access Journals (Sweden)

Chunye Gong

2015-01-01

Full Text Available We present a survey of fractional differential equations and in particular of the computational cost for their numerical solutions from the view of computer science. The computational complexities of time fractional, space fractional, and space-time fractional equations are O(N2M, O(NM2, and O(NM(M + N compared with O(MN for the classical partial differential equations with finite difference methods, where M, N are the number of space grid points and time steps. The potential solutions for this challenge include, but are not limited to, parallel computing, memory access optimization (fractional precomputing operator, short memory principle, fast Fourier transform (FFT based solutions, alternating direction implicit method, multigrid method, and preconditioner technology. The relationships of these solutions for both space fractional derivative and time fractional derivative are discussed. The authors pointed out that the technologies of parallel computing should be regarded as a basic method to overcome this challenge, and some attention should be paid to the fractional killer applications, high performance iteration methods, high order schemes, and Monte Carlo methods. Since the computation of fractional equations with high dimension and variable order is even heavier, the researchers from the area of mathematics and computer science have opportunity to invent cornerstones in the area of fractional calculus.

Development of a computer program for solving the neutronics equations of a multidimensional HTR core model

International Nuclear Information System (INIS)

Schaefer, A.

1979-02-01

A new code for efficient solution of the multidimensional stationary multi-group, diffusion equation, to be used within a HTGR-code model, is presented. The approximation and iteration methods are described. Spacial approximation is based on the QUABOX-coarse-mesh method, but iteration methods are different from QUABOX to give linear dependence of computation time on the number of energy groups. Results for various multidimensional multi-group problems, among them the THTR pebble bed reactor are analyzed. It is shown, that computational labor for a 3D-case is reduced by about a factor 30 in comparison with conventional finite-difference-methods. Thus 3D-full-core calculations appear to be feasible for large HTGR's. (orig.) [de

A computer program to evaluate optical systems

Science.gov (United States)

Innes, D.

1972-01-01

A computer program is used to evaluate a 25.4 cm X-ray telescope at a field angle of 20 minutes of arc by geometrical analysis. The object is regarded as a point source of electromagnetic radiation, and the optical surfaces are treated as boundary conditions in the solution of the electromagnetic wave propagation equation. The electric field distribution is then determined in the region of the image and the intensity distribution inferred. A comparison of wave analysis results and photographs taken through the telescope shows excellent agreement.

Solving the Coupled System Improves Computational Efficiency of the Bidomain Equations

KAUST Repository

Southern, J.A.

2009-10-01

The bidomain equations are frequently used to model the propagation of cardiac action potentials across cardiac tissue. At the whole organ level, the size of the computational mesh required makes their solution a significant computational challenge. As the accuracy of the numerical solution cannot be compromised, efficiency of the solution technique is important to ensure that the results of the simulation can be obtained in a reasonable time while still encapsulating the complexities of the system. In an attempt to increase efficiency of the solver, the bidomain equations are often decoupled into one parabolic equation that is computationally very cheap to solve and an elliptic equation that is much more expensive to solve. In this study, the performance of this uncoupled solution method is compared with an alternative strategy in which the bidomain equations are solved as a coupled system. This seems counterintuitive as the alternative method requires the solution of a much larger linear system at each time step. However, in tests on two 3-D rabbit ventricle benchmarks, it is shown that the coupled method is up to 80% faster than the conventional uncoupled method-and that parallel performance is better for the larger coupled problem.

Solving the Coupled System Improves Computational Efficiency of the Bidomain Equations

KAUST Repository

Southern, J.A.; Plank, G.; Vigmond, E.J.; Whiteley, J.P.

2009-01-01

The bidomain equations are frequently used to model the propagation of cardiac action potentials across cardiac tissue. At the whole organ level, the size of the computational mesh required makes their solution a significant computational challenge. As the accuracy of the numerical solution cannot be compromised, efficiency of the solution technique is important to ensure that the results of the simulation can be obtained in a reasonable time while still encapsulating the complexities of the system. In an attempt to increase efficiency of the solver, the bidomain equations are often decoupled into one parabolic equation that is computationally very cheap to solve and an elliptic equation that is much more expensive to solve. In this study, the performance of this uncoupled solution method is compared with an alternative strategy in which the bidomain equations are solved as a coupled system. This seems counterintuitive as the alternative method requires the solution of a much larger linear system at each time step. However, in tests on two 3-D rabbit ventricle benchmarks, it is shown that the coupled method is up to 80% faster than the conventional uncoupled method-and that parallel performance is better for the larger coupled problem.

Topics in numerical partial differential equations and scientific computing

CERN Document Server

2016-01-01

Numerical partial differential equations (PDEs) are an important part of numerical simulation, the third component of the modern methodology for science and engineering, besides the traditional theory and experiment. This volume contains papers that originated with the collaborative research of the teams that participated in the IMA Workshop for Women in Applied Mathematics: Numerical Partial Differential Equations and Scientific Computing in August 2014.

Error characterization for asynchronous computations: Proxy equation approach

Science.gov (United States)

Sallai, Gabriella; Mittal, Ankita; Girimaji, Sharath

2017-11-01

Numerical techniques for asynchronous fluid flow simulations are currently under development to enable efficient utilization of massively parallel computers. These numerical approaches attempt to accurately solve time evolution of transport equations using spatial information at different time levels. The truncation error of asynchronous methods can be divided into two parts: delay dependent (EA) or asynchronous error and delay independent (ES) or synchronous error. The focus of this study is a specific asynchronous error mitigation technique called proxy-equation approach. The aim of this study is to examine these errors as a function of the characteristic wavelength of the solution. Mitigation of asynchronous effects requires that the asynchronous error be smaller than synchronous truncation error. For a simple convection-diffusion equation, proxy-equation error analysis identifies critical initial wave-number, λc. At smaller wave numbers, synchronous error are larger than asynchronous errors. We examine various approaches to increase the value of λc in order to improve the range of applicability of proxy-equation approach.

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

Science.gov (United States)

Wang, Chi R.

2005-01-01

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation

Multithreaded transactions in scientific computing. The Growth06_v2 program

Science.gov (United States)

Daniluk, Andrzej

2009-07-01

Writing a concurrent program can be more difficult than writing a sequential program. Programmer needs to think about synchronization, race conditions and shared variables. Transactions help reduce the inconvenience of using threads. A transaction is an abstraction, which allows programmers to group a sequence of actions on the program into a logical, higher-level computation unit. This paper presents a new version of the GROWTHGr and GROWTH06 programs. New version program summaryProgram title: GROWTH06_v2 Catalogue identifier: ADVL_v2_1 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVL_v2_1.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 65 255 No. of bytes in distributed program, including test data, etc.: 865 985 Distribution format: tar.gz Programming language: Object Pascal Computer: Pentium-based PC Operating system: Windows 9x, XP, NT, Vista RAM: more than 1 MB Classification: 4.3, 7.2, 6.2, 8, 14 Catalogue identifier of previous version: ADVL_v2_0 Journal reference of previous version: Comput. Phys. Comm. 175 (2006) 678 Does the new version supersede the previous version?: Yes Nature of problem: The programs compute the RHEED intensities during the growth of thin epitaxial structures prepared using the molecular beam epitaxy (MBE). The computations are based on the use of kinematical diffraction theory. Solution method: Epitaxial growth of thin films is modelled by a set of non-linear differential equations [1]. The Runge-Kutta method with adaptive stepsize control was used for solving initial value problem for non-linear differential equations [2]. Reasons for new version: According to the users' suggestions functionality of the program has been improved. Moreover, new use cases have been added which make the handling of the program easier and more

Use of the computer program in a cloud computing

Directory of Open Access Journals (Sweden)

Radovanović Sanja

2013-01-01

Full Text Available Cloud computing represents a specific networking, in which a computer program simulates the operation of one or more server computers. In terms of copyright, all technological processes that take place within the cloud computing are covered by the notion of copying computer programs, and exclusive right of reproduction. However, this right suffers some limitations in order to allow normal use of computer program by users. Based on the fact that the cloud computing is virtualized network, the issue of normal use of the computer program requires to put all aspects of the permitted copying into the context of a specific computing environment and specific processes within the cloud. In this sense, the paper pointed out that the user of a computer program in cloud computing, needs to obtain the consent of the right holder for any act which he undertakes using the program. In other words, the copyright in the cloud computing is a full scale, and thus the freedom of contract (in the case of this particular restriction as well.

Computer models for kinetic equations of magnetically confined plasmas

International Nuclear Information System (INIS)

Killeen, J.; Kerbel, G.D.; McCoy, M.G.; Mirin, A.A.; Horowitz, E.J.; Shumaker, D.E.

1987-01-01

This paper presents four working computer models developed by the computational physics group of the National Magnetic Fusion Energy Computer Center. All of the models employ a kinetic description of plasma species. Three of the models are collisional, i.e., they include the solution of the Fokker-Planck equation in velocity space. The fourth model is collisionless and treats the plasma ions by a fully three-dimensional particle-in-cell method

Sequent Calculus and Equational Programming

Directory of Open Access Journals (Sweden)

Nicolas Guenot

2015-07-01

Full Text Available Proof assistants and programming languages based on type theories usually come in two flavours: one is based on the standard natural deduction presentation of type theory and involves eliminators, while the other provides a syntax in equational style. We show here that the equational approach corresponds to the use of a focused presentation of a type theory expressed as a sequent calculus. A typed functional language is presented, based on a sequent calculus, that we relate to the syntax and internal language of Agda. In particular, we discuss the use of patterns and case splittings, as well as rules implementing inductive reasoning and dependent products and sums.

SNOW: a digital computer program for the simulation of ion beam devices

International Nuclear Information System (INIS)

Boers, J.E.

1980-08-01

A digital computer program, SNOW, has been developed for the simulation of dense ion beams. The program simulates the plasma expansion cup (but not the plasma source itself), the acceleration region, and a drift space with neutralization if desired. The ion beam is simulated by computing representative trajectories through the device. The potentials are simulated on a large rectangular matrix array which is solved by iterative techniques. Poisson's equation is solved at each point within the configuration using space-charge densities computed from the ion trajectories combined with background electron and/or ion distributions. The simulation methods are described in some detail along with examples of both axially-symmetric and rectangular beams. A detailed description of the input data is presented

A PC-based computer program for simulation of containment pressurization

International Nuclear Information System (INIS)

Seifaee, F.

1990-01-01

This paper reports that a PC-based computer program has been developed to simulate a pressurized water reactor (PWR) containment during various transients. This containment model is capable of determining pressure and temperature history of a PWR containment in the event of a loss of coolant accident, as well as main steam line breaks inside the containment. Conservation of mass and energy equations are applied to the containment model. Development of the program is based on minimization of input specified information and user friendliness. Maximization of calculation efficiency is obtained by superseding the traditional trial and error procedure for determination of the state variables and implementation of an explicit solution for pressure. The program includes simplified models for active heat removal systems. The results are in close agreement between the present model and CONTEMPT-MOD5 computer code for pressure and temperature inside the containment

Computational issues of solving the 1D steady gradually varied flow equation

Directory of Open Access Journals (Sweden)

Artichowicz Wojciech

2014-09-01

Full Text Available In this paper a problem of multiple solutions of steady gradually varied flow equation in the form of the ordinary differential energy equation is discussed from the viewpoint of its numerical solution. Using the Lipschitz theorem dealing with the uniqueness of solution of an initial value problem for the ordinary differential equation it was shown that the steady gradually varied flow equation can have more than one solution. This fact implies that the nonlinear algebraic equation approximating the ordinary differential energy equation, which additionally coincides with the wellknown standard step method usually applied for computing of the flow profile, can have variable number of roots. Consequently, more than one alternative solution corresponding to the same initial condition can be provided. Using this property it is possible to compute the water flow profile passing through the critical stage.

Progress report of a research program in computational physics

International Nuclear Information System (INIS)

Guralnik, G.S.

1990-01-01

Task D's research is focused on the understanding of elementary particle physics through the techniques of quantum field theory. We make intensive use of computers to aid our research. During the last year we have made significant progress in understanding the weak interactions through the use of Monte Carlo methods as applied to the equations of quenched lattice QCD. We have launched a program to understand full (not quenched) lattice QCD on relatively large lattices using massively parallel computers. Because of our awareness that Monte Carlo methods might not be able to give a good solution to field theories with the computer power likely to be available to us for the forseeable future we have launched an entirely different numerical approach to study these problems. This ''Source Galerkin'' method is based on an algebraic approach to the field theoretic equations of motion and is (somewhat) related to variational and finite element techniques applied to a source rather than a coordinate space. The results for relatively simple problems are sensationally good. In particular, fermions can be treated in a way which allows them to retain their status as independent dynamical entities in the theory. 8 refs

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

Science.gov (United States)

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

2017-11-01

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

Scilab software as an alternative low-cost computing in solving the linear equations problem

Science.gov (United States)

Agus, Fahrul; Haviluddin

2017-02-01

Numerical computation packages are widely used both in teaching and research. These packages consist of license (proprietary) and open source software (non-proprietary). One of the reasons to use the package is a complexity of mathematics function (i.e., linear problems). Also, number of variables in a linear or non-linear function has been increased. The aim of this paper was to reflect on key aspects related to the method, didactics and creative praxis in the teaching of linear equations in higher education. If implemented, it could be contribute to a better learning in mathematics area (i.e., solving simultaneous linear equations) that essential for future engineers. The focus of this study was to introduce an additional numerical computation package of Scilab as an alternative low-cost computing programming. In this paper, Scilab software was proposed some activities that related to the mathematical models. In this experiment, four numerical methods such as Gaussian Elimination, Gauss-Jordan, Inverse Matrix, and Lower-Upper Decomposition (LU) have been implemented. The results of this study showed that a routine or procedure in numerical methods have been created and explored by using Scilab procedures. Then, the routine of numerical method that could be as a teaching material course has exploited.

Programming in biomolecular computation

DEFF Research Database (Denmark)

Hartmann, Lars Røeboe; Jones, Neil; Simonsen, Jakob Grue

2010-01-01

in a strong sense: a universal algorithm exists, that is able to execute any program, and is not asymptotically inefficient. A prototype model has been implemented (for now in silico on a conventional computer). This work opens new perspectives on just how computation may be specified at the biological level......., by programs reminiscent of low-level computer machine code; and at the same time biologically plausible: its functioning is defined by a single and relatively small set of chemical-like reaction rules. Further properties: the model is stored-program: programs are the same as data, so programs are not only...... executable, but are also compilable and interpretable. It is universal: all computable functions can be computed (in natural ways and without arcane encodings of data and algorithm); it is also uniform: new “hardware” is not needed to solve new problems; and (last but not least) it is Turing complete...

Introduction of computing in physics learning visual programing

International Nuclear Information System (INIS)

Kim, Cheung Seop

1999-12-01

This book introduces physics and programing, foundation of visual basic, grammar of visual basic, visual programing, solution of equation, calculation of matrix, solution of simultaneous equation, differentiation, differential equation, simultaneous differential equation and second-order differential equation, integration and solution of partial differential equation. It also covers basic language, terms of visual basic, usage of method, graphic method, step by step method, fails-position method, Gauss elimination method, difference method and Euler method.

SHOCK-JR: a computer program to analyze impact response of shipping container

International Nuclear Information System (INIS)

Ikushima, Takeshi; Nakazato, Chikara; Shimoda, Osamu; Uchino, Mamoru.

1983-02-01

The report is provided for using a computer program, SHOCK-JR, which is used to analyze the impact response of shipping containers. Descriptions are the mathematical model, method of analysis, structures of the program and the input and output variables. The program solves the equations of motion for a one-dimensional, lumped mass and nonlinear spring model. The solution procedure uses Runge-Kutta-Gill and Newmark-β methods. SHOCK-JR is a revised version of SHOCK, which was developed by ORNL. In SHOCK-JR, SI dimension is used and graphical output is available. (author)

PRETTA：A COMPUTER PROGRAM FOR PWR PRESSURIZER’S TRANSIENT THERMODYNAMICS

Institute of Scientific and Technical Information of China (English)

阿谢德; 徐济鋆

2001-01-01

A computer program PRETTA “Pressurizer Transient Thermodynamics Analysis” was developed for the prediction of pressurizer under transient conditions. It is based on the solution of the conservation laws of heat and mass applied to the three separate and non equilibrium thermodynamic regions. In the program all of the important thermal-hydraulics phenomena occurring in the pressurizer: stratification of the hot water and incoming cold water, bulk flashing and condensation, wall condensation, and interfacial heat and mass transfer have been considered. The bubble rising and rain-out models are developed to describe bulk flashing and condensation, respectively. To obtain the wall condensation rate, a one-dimensional heat conduction equation is solved by the pivoting method. The presented computer program will predict the pressure-time behavior of a PWR pressurizer during a variety of transients. The results obtained from the proposed mathematical model are in good agreement with available data on the CHASHMA nuclear power plant's pressurizer performance.

Relativistic Photoionization Computations with the Time Dependent Dirac Equation

Science.gov (United States)

2016-10-12

Naval Research Laboratory Washington, DC 20375-5320 NRL/MR/6795--16-9698 Relativistic Photoionization Computations with the Time Dependent Dirac... Photoionization Computations with the Time Dependent Dirac Equation Daniel F. Gordon and Bahman Hafizi Naval Research Laboratory 4555 Overlook Avenue, SW...Unclassified Unlimited Unclassified Unlimited 22 Daniel Gordon (202) 767-5036 Tunneling Photoionization Ionization of inner shell electrons by laser

ASYS: a computer algebra package for analysis of nonlinear algebraic equations systems

International Nuclear Information System (INIS)

Gerdt, V.P.; Khutornoj, N.V.

1992-01-01

A program package ASYS for analysis of nonlinear algebraic equations based on the Groebner basis technique is described. The package is written in REDUCE computer algebra language. It has special facilities to treat polynomial ideals of positive dimension, corresponding to algebraic systems with infinitely many solutions. Such systems can be transformed to an equivalent set of subsystems with reduced number of variables in completely automatic way. It often allows to construct the explicit form of a solution set in many problems of practical importance. Some examples and results of comparison with the standard Reduce package GROEBNER and special-purpose systems FELIX and A1PI are given. 21 refs.; 2 tabs

Nonlinear evolution equations and solving algebraic systems: the importance of computer algebra

International Nuclear Information System (INIS)

Gerdt, V.P.; Kostov, N.A.

1989-01-01

In the present paper we study the application of computer algebra to solve the nonlinear polynomial systems which arise in investigation of nonlinear evolution equations. We consider several systems which are obtained in classification of integrable nonlinear evolution equations with uniform rank. Other polynomial systems are related with the finding of algebraic curves for finite-gap elliptic potentials of Lame type and generalizations. All systems under consideration are solved using the method based on construction of the Groebner basis for corresponding polynomial ideals. The computations have been carried out using computer algebra systems. 20 refs

Decomposition and Cross-Product-Based Method for Computing the Dynamic Equation of Robots

Directory of Open Access Journals (Sweden)

Ching-Long Shih

2012-08-01

Full Text Available This paper aims to demonstrate a clear relationship between Lagrange equations and Newton-Euler equations regarding computational methods for robot dynamics, from which we derive a systematic method for using either symbolic or on-line numerical computations. Based on the decomposition approach and cross-product operation, a computing method for robot dynamics can be easily developed. The advantages of this computing framework are that: it can be used for both symbolic and on-line numeric computation purposes, and it can also be applied to biped systems, as well as some simple closed-chain robot systems.

Proxy-equation paradigm: A strategy for massively parallel asynchronous computations

Science.gov (United States)

Mittal, Ankita; Girimaji, Sharath

2017-09-01

Massively parallel simulations of transport equation systems call for a paradigm change in algorithm development to achieve efficient scalability. Traditional approaches require time synchronization of processing elements (PEs), which severely restricts scalability. Relaxing synchronization requirement introduces error and slows down convergence. In this paper, we propose and develop a novel "proxy equation" concept for a general transport equation that (i) tolerates asynchrony with minimal added error, (ii) preserves convergence order and thus, (iii) expected to scale efficiently on massively parallel machines. The central idea is to modify a priori the transport equation at the PE boundaries to offset asynchrony errors. Proof-of-concept computations are performed using a one-dimensional advection (convection) diffusion equation. The results demonstrate the promise and advantages of the present strategy.

ICASE Computer Science Program

Science.gov (United States)

1985-01-01

The Institute for Computer Applications in Science and Engineering computer science program is discussed in outline form. Information is given on such topics as problem decomposition, algorithm development, programming languages, and parallel architectures.

Inverse kinetics equations for on line measurement of reactivity using personal computer

International Nuclear Information System (INIS)

Ratemi, Wajdi; El Gadamsi, Walied; Beleid, Abdul Kariem

1993-01-01

Computer with their astonishing speed of calculations along with their easy connection to real systems, are very appropriate for digital measurements of real system variables. In the nuclear industry, such computer application will produce compact control rooms of real power plants, where information and results display can be obtained through push button concept. In our study, we use two personal computers for the purpose of simulation and measurement. One of them is used as a digital simulator to a real reactor, where we effectively simulate the reactor power through a cross talk network. The computed power is passed at certain chosen sampling time to the other computer. The purpose of the other computer is to use the inverse kinetics equations to calculate the reactivity parameter based on the received power and then it performs on line display of the power curve and the reactivity curve using color graphics. In this study, we use the one group version of the inverse kinetics algorithm which can easily be extended to larger group version. The language of programming used in Turbo BASIC, which is very comparable, in terms of efficiency, to FORTRAN language, besides its effective graphics routines. With the use of the extended version of the Inverse Kinetics algorithm, we can effectively apply this techniques of measurement for the purpose of on line display of the reactivity of the Tajoura Research Reactor. (author)

REFLECT: a program to integrate the wave equation through a plane stratified plasma

International Nuclear Information System (INIS)

Greene, J.W.

1975-01-01

A program was developed to integrate the wave equation through a plane stratified plasma with a general density distribution. The reflection and transmission of a plane wave are computed as a function of the angle of incidence. The polarization of the electric vector is assumed to be perpendicular to the plane of incidence. The model for absorption by classical inverse bremsstrahlung avoids the improper extrapolation of underdense formulae that are singular at the plasma critical surface. Surprisingly good agreement with the geometric-optics analysis of a linear layer was found. The system of ordinary differential equations is integrated by the variable-step, variable-order Adams method in the Lawrence Livermore Laboratory Gear package. Parametric studies of the absorption are summarized, and some possibilities for further development of the code are discussed. (auth)

FISPRO: a simplified computer program for general fission product formation and decay calculations

International Nuclear Information System (INIS)

Jiacoletti, R.J.; Bailey, P.G.

1979-08-01

This report describes a computer program that solves a general form of the fission product formation and decay equations over given time steps for arbitrary decay chains composed of up to three nuclides. All fission product data and operational history data are input through user-defined input files. The program is very useful in the calculation of fission product activities of specific nuclides for various reactor operational histories and accident consequence calculations

The digital computer

CERN Document Server

Parton, K C

2014-01-01

The Digital Computer focuses on the principles, methodologies, and applications of the digital computer. The publication takes a look at the basic concepts involved in using a digital computer, simple autocode examples, and examples of working advanced design programs. Discussions focus on transformer design synthesis program, machine design analysis program, solution of standard quadratic equations, harmonic analysis, elementary wage calculation, and scientific calculations. The manuscript then examines commercial and automatic programming, how computers work, and the components of a computer

Periodicity computation of generalized mathematical biology problems involving delay differential equations.

Science.gov (United States)

Jasim Mohammed, M; Ibrahim, Rabha W; Ahmad, M Z

2017-03-01

In this paper, we consider a low initial population model. Our aim is to study the periodicity computation of this model by using neutral differential equations, which are recognized in various studies including biology. We generalize the neutral Rayleigh equation for the third-order by exploiting the model of fractional calculus, in particular the Riemann-Liouville differential operator. We establish the existence and uniqueness of a periodic computational outcome. The technique depends on the continuation theorem of the coincidence degree theory. Besides, an example is presented to demonstrate the finding.

Markovian Monte Carlo program EvolFMC v.2 for solving QCD evolution equations

Science.gov (United States)

Jadach, S.; Płaczek, W.; Skrzypek, M.; Stokłosa, P.

2010-02-01

We present the program EvolFMC v.2 that solves the evolution equations in QCD for the parton momentum distributions by means of the Monte Carlo technique based on the Markovian process. The program solves the DGLAP-type evolution as well as modified-DGLAP ones. In both cases the evolution can be performed in the LO or NLO approximation. The quarks are treated as massless. The overall technical precision of the code has been established at 5×10. This way, for the first time ever, we demonstrate that with the Monte Carlo method one can solve the evolution equations with precision comparable to the other numerical methods. New version program summaryProgram title: EvolFMC v.2 Catalogue identifier: AEFN_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFN_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including binary test data, etc.: 66 456 (7407 lines of C++ code) No. of bytes in distributed program, including test data, etc.: 412 752 Distribution format: tar.gz Programming language: C++ Computer: PC, Mac Operating system: Linux, Mac OS X RAM: Less than 256 MB Classification: 11.5 External routines: ROOT ( http://root.cern.ch/drupal/) Nature of problem: Solution of the QCD evolution equations for the parton momentum distributions of the DGLAP- and modified-DGLAP-type in the LO and NLO approximations. Solution method: Monte Carlo simulation of the Markovian process of a multiple emission of partons. Restrictions:Limited to the case of massless partons. Implemented in the LO and NLO approximations only. Weighted events only. Unusual features: Modified-DGLAP evolutions included up to the NLO level. Additional comments: Technical precision established at 5×10. Running time: For the 10 6 events at 100 GeV: DGLAP NLO: 27s; C-type modified DGLAP NLO: 150s (MacBook Pro with Mac OS X v.10

ALMOD-JRC computer program

International Nuclear Information System (INIS)

Cacciabue, P.C.; Lisanti, B.; Tozzi, A.

1984-01-01

This paper discusses the details concerning the newly developed or modified models of the computer program ALMOD-JRC, originating from ALMOD 3/Rel 4. The most important argument for the implementation of the new models was the need to enlarge the spectrum of the simulated phenomena, and to improve the simulation of experimental facilities such as LOFT or LOBI. This has led to a better formulation of the heat transfer and pressure drops correlations and to the implementation of the treatment of the heat losses to structural materials. In particular a series of test cases on real power plants, a pre-test examination of a LOBI station blackout ATWS experiment and the post test analysis of the L9-3 experiment, show the ability of ALMOD-JRC to correctly simulate PWR incident sequences. Although in ALMOD-JRC the code capabilities have been expanded, the limitations of the original version of the program still hold for what concerns the treatment of the coolant thermohydraulics as homogeneous flow for the two phase conditions in the primary coolant circuit. The other interesting feature of the new code is the remarkably shorter running times obtained with the introduction of simplified numerical treatments for the solving equations, without significant loss of accuracy of results

Analyses of computer programs for the probabilistic estimation of design earthquake and seismological characteristics of the Korean Peninsula

International Nuclear Information System (INIS)

Lee, Gi Hwa

1997-11-01

The purpose of the present study is to develop predictive equations from simulated motions which are adequate for the Korean Peninsula and analyze and utilize the computer programs for the probabilistic estimation of design earthquakes. In part I of the report, computer programs for the probabilistic estimation of design earthquake are analyzed and applied to the seismic hazard characterizations in the Korean Peninsula. In part II of the report, available instrumental earthquake records are analyzed to estimate earthquake source characteristics and medium properties, which are incorporated into simulation process. And earthquake records are simulated by using the estimated parameters. Finally, predictive equations constructed from the simulation are given in terms of magnitude and hypocentral distances

MAGNETOHYDRODYNAMIC EQUATIONS (MHD GENERATION CODE

Directory of Open Access Journals (Sweden)

Francisco Frutos Alfaro

2017-04-01

Full Text Available A program to generate codes in Fortran and C of the full magnetohydrodynamic equations is shown. The program uses the free computer algebra system software REDUCE. This software has a package called EXCALC, which is an exterior calculus program. The advantage of this program is that it can be modified to include another complex metric or spacetime. The output of this program is modified by means of a LINUX script which creates a new REDUCE program to manipulate the magnetohydrodynamic equations to obtain a code that can be used as a seed for a magnetohydrodynamic code for numerical applications. As an example, we present part of the output of our programs for Cartesian coordinates and how to do the discretization.

Adolescents' Chunking of Computer Programs.

Science.gov (United States)

Magliaro, Susan; Burton, John K.

To investigate what children learn during computer programming instruction, students attending a summer computer camp were asked to recall either single lines or chunks of computer programs from either coherent or scrambled programs. The 16 subjects, ages 12 to 17, were divided into three instructional groups: (1) beginners, who were taught to…

Computer programs for the solution of systems of linear algebraic equations

Science.gov (United States)

Sequi, W. T.

1973-01-01

FORTRAN subprograms for the solution of systems of linear algebraic equations are described, listed, and evaluated in this report. Procedures considered are direct solution, iteration, and matrix inversion. Both incore methods and those which utilize auxiliary data storage devices are considered. Some of the subroutines evaluated require the entire coefficient matrix to be in core, whereas others account for banding or sparceness of the system. General recommendations relative to equation solving are made, and on the basis of tests, specific subprograms are recommended.

Computer programs as accounting object

Directory of Open Access Journals (Sweden)

I.V. Perviy

2015-03-01

Full Text Available Existing approaches to the regulation of accounting software as one of the types of intangible assets have been considered. The features and current state of the legal protection of computer programs have been analyzed. The reasons for the need to use patent law as a means of legal protection of individual elements of computer programs have been discovered. The influence of the legal aspects of the use of computer programs for national legislation to their accounting reflection has been analyzed. The possible options for the transfer of rights from computer programs copyright owners have been analyzed that should be considered during creation of software accounting system at the enterprise. Identified and analyzed the characteristics of computer software as an intangible asset under the current law. General economic characteristics of computer programs as one of the types of intangible assets have been grounded. The main distinguishing features of software compared to other types of intellectual property have been all ocated

NASA's computer science research program

Science.gov (United States)

Larsen, R. L.

1983-01-01

Following a major assessment of NASA's computing technology needs, a new program of computer science research has been initiated by the Agency. The program includes work in concurrent processing, management of large scale scientific databases, software engineering, reliable computing, and artificial intelligence. The program is driven by applications requirements in computational fluid dynamics, image processing, sensor data management, real-time mission control and autonomous systems. It consists of university research, in-house NASA research, and NASA's Research Institute for Advanced Computer Science (RIACS) and Institute for Computer Applications in Science and Engineering (ICASE). The overall goal is to provide the technical foundation within NASA to exploit advancing computing technology in aerospace applications.

Computer Program Newsletter No. 7

International Nuclear Information System (INIS)

Magnuson, W.G. Jr.

1982-09-01

This issue of the Computer Program Newsletter updates an earlier newsletter (Number 2, September 1979) and focuses on electrical network analysis computer programs. In particular, five network analysis programs (SCEPTRE, SPICE2, NET2, CALAHAN, and EMTP) will be described. The objective of this newsletter will be to provide a very brief description of the input syntax and semantics for each program, highlight their strong and weak points, illustrate how the programs are run at Lawrence Livermore National Laboratory using the Octopus computer network, and present examples of input for each of the programs to illustrate some of the features of each program. In a sense, this newsletter can be used as a quick reference guide to the programs

CACHE: an extended BASIC program which computes the performance of shell and tube heat exchangers

International Nuclear Information System (INIS)

Tallackson, J.R.

1976-03-01

An extended BASIC program, CACHE, has been written to calculate steady state heat exchange rates in the core auxiliary heat exchangers, (CAHE), designed to remove afterheat from High-Temperature Gas-Cooled Reactors (HTGR). Computationally, these are unbaffled counterflow shell and tube heat exchangers. The computational method is straightforward. The exchanger is subdivided into a user-selected number of lengthwise segments; heat exchange in each segment is calculated in sequence and summed. The program takes the temperature dependencies of all thermal conductivities, viscosities and heat capacities into account providing these are expressed algebraically. CACHE is easily adapted to compute steady state heat exchange rates in any unbaffled counterflow exchanger. As now used, CACHE calculates heat removal by liquid weight from high-temperature helium and helium mixed with nitrogen, oxygen and carbon monoxide. A second program, FULTN, is described. FULTN computes the geometrical parameters required as input to CACHE. As reported herein, FULTN computes the internal dimensions of the Fulton Station CAHE. The two programs are chained to operate as one. Complete user information is supplied. The basic equations, variable lists, annotated program lists, and sample outputs with explanatory notes are included

Kinetic equation solution by inverse kinetic method

International Nuclear Information System (INIS)

Salas, G.

1983-01-01

We propose a computer program (CAMU) which permits to solve the inverse kinetic equation. The CAMU code is written in HPL language for a HP 982 A microcomputer with a peripheral interface HP 9876 A ''thermal graphic printer''. The CAMU code solves the inverse kinetic equation by taking as data entry the output of the ionization chambers and integrating the equation with the help of the Simpson method. With this program we calculate the evolution of the reactivity in time for a given disturbance

CINESP - computational program of spatial kinetics for nuclear reactors in the one-two dimension multigroup diffusion theory

International Nuclear Information System (INIS)

Santos, R.S. dos

1993-01-01

This paper presents a computational program to solve numerically the reactor kinetics equations in the multigroup diffusion theory. One or two-dimensional problems in cylindrical or Cartesian geometries, with any number of energy and delayed-neutron precursors groups are dealt with. The main input and output of the program are briefly discussed. Various results demonstrate the accuracy and versatility of the program, when compared with other kinetics programs. (author)

Computer technology and computer programming research and strategies

CERN Document Server

Antonakos, James L

2011-01-01

Covering a broad range of new topics in computer technology and programming, this volume discusses encryption techniques, SQL generation, Web 2.0 technologies, and visual sensor networks. It also examines reconfigurable computing, video streaming, animation techniques, and more. Readers will learn about an educational tool and game to help students learn computer programming. The book also explores a new medical technology paradigm centered on wireless technology and cloud computing designed to overcome the problems of increasing health technology costs.

Program MASTERCALC: an interactive computer program for radioanalytical computations. Description and operating instructions

International Nuclear Information System (INIS)

Goode, W.

1980-10-01

MASTERCALC is a computer program written to support radioanalytical computations in the Los Alamos Scientific Laboratory (LASL) Environmental Surveillance Group. Included in the program are routines for gross alpha and beta, 3 H, gross gamma, 90 Sr and alpha spectroscopic determinations. A description of MASTERCALC is presented and its source listing is included. Operating instructions and example computing sessions are given for each type of analysis

A programming approach to computability

CERN Document Server

Kfoury, A J; Arbib, Michael A

1982-01-01

Computability theory is at the heart of theoretical computer science. Yet, ironically, many of its basic results were discovered by mathematical logicians prior to the development of the first stored-program computer. As a result, many texts on computability theory strike today's computer science students as far removed from their concerns. To remedy this, we base our approach to computability on the language of while-programs, a lean subset of PASCAL, and postpone consideration of such classic models as Turing machines, string-rewriting systems, and p. -recursive functions till the final chapter. Moreover, we balance the presentation of un solvability results such as the unsolvability of the Halting Problem with a presentation of the positive results of modern programming methodology, including the use of proof rules, and the denotational semantics of programs. Computer science seeks to provide a scientific basis for the study of information processing, the solution of problems by algorithms, and the design ...

Program to solve the multigroup discrete ordinates transport equation in (x,y,z) geometry

International Nuclear Information System (INIS)

Lathrop, K.D.

1976-04-01

Numerical formulations and programming algorithms are given for the THREETRAN computer program which solves the discrete ordinates, multigroup transport equation in (x,y,z) geometry. An efficient, flexible, and general data-handling strategy is derived to make use of three hierarchies of storage: small core memory, large core memory, and disk file. Data management, input instructions, and sample problem output are described. A six-group, S 4 , 18 502 mesh point, 2 800 zone, k/sub eff/ calculation of the ZPPR-4 critical assembly required 144 min of CDC-7600 time to execute to a convergence tolerance of 5 x 10 -4 and gave results in good qualitative agreement with experiment and other calculations. 6 references

Solutions to estimation problems for scalar hamilton-jacobi equations using linear programming

KAUST Repository

Claudel, Christian G.; Chamoin, Timothee; Bayen, Alexandre M.

2014-01-01

This brief presents new convex formulations for solving estimation problems in systems modeled by scalar Hamilton-Jacobi (HJ) equations. Using a semi-analytic formula, we show that the constraints resulting from a HJ equation are convex, and can be written as a set of linear inequalities. We use this fact to pose various (and seemingly unrelated) estimation problems related to traffic flow-engineering as a set of linear programs. In particular, we solve data assimilation and data reconciliation problems for estimating the state of a system when the model and measurement constraints are incompatible. We also solve traffic estimation problems, such as travel time estimation or density estimation. For all these problems, a numerical implementation is performed using experimental data from the Mobile Century experiment. In the context of reproducible research, the code and data used to compute the results presented in this brief have been posted online and are accessible to regenerate the results. © 2013 IEEE.

An introduction to differential equations using MATLAB

CERN Document Server

Butt, Rizwan

2016-01-01

An Introduction to Differential Equations using MATLAB exploits the symbolic, numerical, and graphical capabilitiesof MATLAB to develop a thorough understanding of differential equations algorithms. This book provides the readerwith numerous applications, m-files, and practical examples to problems. Balancing theoretical concepts withcomputational speed and accuracy, the book includes numerous short programs in MATLAB that can be used to solveproblems involving first-and higher-order differential equations, Laplace transforms, linear systems of differentialequations, numerical solutions of differential equations, computer graphics, and more. The author emphasizes thebasic ideas of analytical and numerical techniques and the uses of modern mathematical software (MATLAB) ratherthan relying only on complex mathematical derivations to engineers, mathematician, computer scientists, andphysicists or for use as a textbook in applied or computational courses.A CD-ROM with all the figures, codes, solutions, appendices...

Numerical computation of soliton dynamics for NLS equations in a driving potential

Directory of Open Access Journals (Sweden)

Marco Caliari

2010-06-01

Full Text Available We provide numerical computations for the soliton dynamics of the nonlinear Schrodinger equation with an external potential. After computing the ground state solution r of a related elliptic equation we show that, in the semi-classical regime, the center of mass of the solution with initial datum built upon r is driven by the solution to $ddot x=- abla V(x$. Finally, we provide examples and analyze the numerical errors in the two dimensional case when V is a harmonic potential.

Stochastic equations for complex systems theoretical and computational topics

CERN Document Server

Bessaih, Hakima

2015-01-01

Mathematical analyses and computational predictions of the behavior of complex systems are needed to effectively deal with weather and climate predictions, for example, and the optimal design of technical processes. Given the random nature of such systems and the recognized relevance of randomness, the equations used to describe such systems usually need to involve stochastics.  The basic goal of this book is to introduce the mathematics and application of stochastic equations used for the modeling of complex systems. A first focus is on the introduction to different topics in mathematical analysis. A second focus is on the application of mathematical tools to the analysis of stochastic equations. A third focus is on the development and application of stochastic methods to simulate turbulent flows as seen in reality.  This book is primarily oriented towards mathematics and engineering PhD students, young and experienced researchers, and professionals working in the area of stochastic differential equations ...

Computer program CDCID: an automated quality control program using CDC update

International Nuclear Information System (INIS)

Singer, G.L.; Aguilar, F.

1984-04-01

A computer program, CDCID, has been developed in coordination with a quality control program to provide a highly automated method of documenting changes to computer codes at EG and G Idaho, Inc. The method uses the standard CDC UPDATE program in such a manner that updates and their associated documentation are easily made and retrieved in various formats. The method allows each card image of a source program to point to the document which describes it, who created the card, and when it was created. The method described is applicable to the quality control of computer programs in general. The computer program described is executable only on CDC computing systems, but the program could be modified and applied to any computing system with an adequate updating program

Some Matrix Iterations for Computing Generalized Inverses and Balancing Chemical Equations

Directory of Open Access Journals (Sweden)

Farahnaz Soleimani

2015-11-01

Full Text Available An application of iterative methods for computing the Moore–Penrose inverse in balancing chemical equations is considered. With the aim to illustrate proposed algorithms, an improved high order hyper-power matrix iterative method for computing generalized inverses is introduced and applied. The improvements of the hyper-power iterative scheme are based on its proper factorization, as well as on the possibility to accelerate the iterations in the initial phase of the convergence. Although the effectiveness of our approach is confirmed on the basis of the theoretical point of view, some numerical comparisons in balancing chemical equations, as well as on randomly-generated matrices are furnished.

Computational programs for shielding calculation with transport of one dimensional and monoenergetic SN

International Nuclear Information System (INIS)

Nunes, Carlos Eduardo A.; Barros, Ricardo C.

2009-01-01

This paper describes a computational program for result simulation of neutron transport problems at one velocity with isotropic scattering in Cartesian onedimensional geometry. Describing the physical modelling, the next phase is a mathematical modelling of the physical problem for simulation of the neutron distribution. The mathematical modelling uses the linearized Boltzmann equation which represents a balance among the production and loss of particles. The formulation of the discrete ordinates S N consists of discretization of angular variables at N directions (discrete ordinates), and using a set of angular quadratures for the approximation of integral terms of scattering sources. The S N equations are numerically solved. This work describes three numerical methods: diamond difference, step and characteristic step. The paper also presents numerical results for illustration of the efficiency of the developed program

Proposed equation of state experimental program at NOVA/NIF

International Nuclear Information System (INIS)

Holmes, N.C.; Cauble, R.; Celliers, P.; Collins, G.; DaSilva, L.; Hammel, B.; Stewart, R.; Strand, O.; Sullivan, A.

1996-03-01

The high pressure equations of state of materials must be accurately known for confident design, simulation, and interpretation of ICF target experiments and for weapons. Here, the authors sketch out a program to perform precise and accurate equation of state (EOS) experiments using large, high-power lasers. The program is divided into three phases: (1) a driver qualification program which will determine the necessary drive and target design parameters; (2) characterization of low and high-Z standard materials; (3) accurate impedance-match experiments to determine equations of state of materials of interest relative to the qualified standards. Novel methods are proposed for the first phase which are an order of magnitude more sensitive than those applied previously. They identify and describe the choices for standards and the samples slated for initial studies. In all cases, their goal is to obtain data that is of sufficiently high quality to serve for design purposes and to validate theories of material response at high pressures of the order of 1--10 TPa. This is particularly important right now, for materials such as CH plastics and hydrogen (DT). In addition, they will suggest a program for EOS measurements for P > 10 TPa, the region of maximum theoretical EOS uncertainty for many materials

Computer Programming Education with Miranda

NARCIS (Netherlands)

Joosten, S.M.M.; van den Berg, Klaas

During the past four years, an experiment has been carried out with an introductory course in computer programming, based on functional programming. This article describes the background of this approach, the aim of the computer programming course, the outline and subject matter of the course parts

SOLVEX: a computer program for simulation of solvent extraction processes

International Nuclear Information System (INIS)

Scotten, W.C.

1975-09-01

SOLVEX is a FORTRAN IV computer program that simulates the dynamic behavior of solvent extraction processes conducted in mixer-settlers and centrifugal contactors. Two options permit terminating dynamic phases by time or by achieving steady state, and a third option permits artificial rapid close to steady state. Thus the program is well suited to multiple phases of dynamic problems and multiple input of steady state problems. Changes from the previous problem are the only inputs required for each succeeding problem. Distribution data can be supplied by two-variable third-power polynomial equations or by three-variable tables in any one of 16 different combinations involving phase concentrations or distribution coefficients (ratio of phase concentrations) or their logarithms

Mississippi Curriculum Framework for Computer Information Systems Technology. Computer Information Systems Technology (Program CIP: 52.1201--Management Information Systems & Business Data). Computer Programming (Program CIP: 52.1201). Network Support (Program CIP: 52.1290--Computer Network Support Technology). Postsecondary Programs.

Science.gov (United States)

Mississippi Research and Curriculum Unit for Vocational and Technical Education, State College.

This document, which is intended for use by community and junior colleges throughout Mississippi, contains curriculum frameworks for two programs in the state's postsecondary-level computer information systems technology cluster: computer programming and network support. Presented in the introduction are program descriptions and suggested course…

Diffusive Wave Approximation to the Shallow Water Equations: Computational Approach

KAUST Repository

Collier, Nathan; Radwan, Hany; Dalcin, Lisandro; Calo, Victor M.

2011-01-01

We discuss the use of time adaptivity applied to the one dimensional diffusive wave approximation to the shallow water equations. A simple and computationally economical error estimator is discussed which enables time-step size adaptivity

AV Programs for Computer Know-How.

Science.gov (United States)

Mandell, Phyllis Levy

1985-01-01

Lists 44 audiovisual programs (most released between 1983 and 1984) grouped in seven categories: computers in society, introduction to computers, computer operations, languages and programing, computer graphics, robotics, computer careers. Excerpts from "School Library Journal" reviews, price, and intended grade level are included. Names…

Programmable calculator programs to solve softwood volume and value equations.

Science.gov (United States)

Janet K. Ayer. Sachet

1982-01-01

This paper presents product value and product volume equations as programs for handheld calculators. These tree equations are for inland Douglas-fir, young-growth Douglas-fir, western white pine, ponderosa pine, and western larch. Operating instructions and an example are included.

Computer Networking Laboratory for Undergraduate Computer Technology Program

National Research Council Canada - National Science Library

Naghedolfeizi, Masoud

2000-01-01

...) To improve the quality of education in the existing courses related to computer networks and data communications as well as other computer science courses such programming languages and computer...

RECON: a computer program for analyzing repository economics. Documentation and user's manual

International Nuclear Information System (INIS)

Clark, L.L.; Cole, B.M.; McNair, G.W.; Schutz, M.E.

1983-05-01

From 1981 through 1983 the Pacific Northwest Laboratory has been developing a computer model named RECON to calculate repository costs from parametric data input. The objective of the program has been to develop the capability to evalute the effect on costs of changes in repository design parameters and operating scenario assumptions. This report documents the development of the model through March of 1983. Included in the report are: (1) descriptions of model development and the underlying equations, assumptions and definitions; (2) descriptions of data input either using card images or an interactive data input program; and (3) detailed listings of the program and definitions of program variables. Cost estimates generated using the model have been verified against independent estimates and good agreement has been obtained

Computing generalized Langevin equations and generalized Fokker-Planck equations.

Science.gov (United States)

Darve, Eric; Solomon, Jose; Kia, Amirali

2009-07-07

The Mori-Zwanzig formalism is an effective tool to derive differential equations describing the evolution of a small number of resolved variables. In this paper we present its application to the derivation of generalized Langevin equations and generalized non-Markovian Fokker-Planck equations. We show how long time scales rates and metastable basins can be extracted from these equations. Numerical algorithms are proposed to discretize these equations. An important aspect is the numerical solution of the orthogonal dynamics equation which is a partial differential equation in a high dimensional space. We propose efficient numerical methods to solve this orthogonal dynamics equation. In addition, we present a projection formalism of the Mori-Zwanzig type that is applicable to discrete maps. Numerical applications are presented from the field of Hamiltonian systems.

The computer code EURDYN - 1 M (release 1) for transient dynamic fluid-structure interaction. Pt.1: governing equations and finite element modelling

International Nuclear Information System (INIS)

Donea, J.; Fasoli-Stella, P.; Giuliani, S.; Halleux, J.P.; Jones, A.V.

1980-01-01

This report describes the governing equations and the finite element modelling used in the computer code EURDYN - 1 M. The code is a non-linear transient dynamic program for the analysis of coupled fluid-structure systems; It is designed for safety studies on LMFBR components (primary containment and fuel subassemblies)

An appraisal of computational techniques for transient heat conduction equation

International Nuclear Information System (INIS)

Kant, T.

1983-01-01

A semi-discretization procedure in which the ''space'' dimension is discretized by the finite element method is emphasized for transient problems. This standard methodology transforms the space-time partial differential equation (PDE) system into a set of ordinary differential equations (ODE) in time. Existing methods for transient heat conduction calculations are then reviewed. Existence of two general classes of time integration schemes- implicit and explicit is noted. Numerical stability characteristics of these two methods are elucidated. Implicit methods are noted to be numerically stable, permitting large time steps, but the cost per step is high. On the otherhand, explicit schemes are noted to be inexpensive per step, but small step size is required. Low computational cost of the explicit schemes make it very attractive for nonlinear problems. However, numerical stability considerations requiring use of very small time steps come in the way of its general adoption. Effectiveness of the fourth-order Runge-Kutta-Gill explicit integrator is then numerically evaluated. Finally we discuss some very recent works on development of computational algorithms which not only achieve unconditional stability, high accuracy and convergence but involve computations on matrix equations of elements only. This development is considered to be very significant in the light of our experience gained for simple heat conduction calculations. We conclude that such algorithms have the potential for further developments leading to development of economical methods for general transient analysis of complex physical systems. (orig.)

AGRIS: Description of computer programs

International Nuclear Information System (INIS)

Schmid, H.; Schallaboeck, G.

1976-01-01

The set of computer programs used at the AGRIS (Agricultural Information System) Input Unit at the IAEA, Vienna, Austria to process the AGRIS computer-readable data is described. The processing flow is illustrated. The configuration of the IAEA's computer, a list of error messages generated by the computer, the EBCDIC code table extended for AGRIS and INIS, the AGRIS-6 bit code, the work sheet format, and job control listings are included as appendixes. The programs are written for an IBM 370, model 145, operating system OS or VS, and require a 130K partition. The programming languages are PL/1 (F-compiler) and Assembler

Recommended programming practices to facilitate the portability of science computer programs

International Nuclear Information System (INIS)

Anon.

1983-01-01

This standard recommends programming practices to facilitate the portability of computer programs prepared for scientific and engineering computations. These practices are intended to simplify implementation, conversion, and modification of computer programs

Integer programming theory, applications, and computations

CERN Document Server

Taha, Hamdy A

1975-01-01

Integer Programming: Theory, Applications, and Computations provides information pertinent to the theory, applications, and computations of integer programming. This book presents the computational advantages of the various techniques of integer programming.Organized into eight chapters, this book begins with an overview of the general categorization of integer applications and explains the three fundamental techniques of integer programming. This text then explores the concept of implicit enumeration, which is general in a sense that it is applicable to any well-defined binary program. Other

Development of a Computer Program for the Analysis Logistics of PWR Spent Fuels

International Nuclear Information System (INIS)

Choi, Heui Joo; Choi, Jong Won; Cha, Jeong Hun

2008-01-01

It is expected that the temporary storage facilities at the nuclear power plants will be full of the spent fuels within 10 years. Provided that a centralized interim storage facility is constructed along the coast of the Korean peninsula to solve this problem, a substantial amount of spent fuels should be transported by sea or by land every year. In this paper we developed a computer program for the analysis of transportation logistics of the spent fuels from 4 different nuclear power plant sites to the hypothetical centralized interim storage facility and the final repository. Mass balance equations were used to analyze the logistics between the nuclear power plants and the interim storage facility. To this end a computer program, CASK, was developed by using the VISUAL BASIC language. The annual transportation rates of spent fuels from the four nuclear power plant sites were determined by using the CASK program. The parameter study with the program illustrated the easiness of logistics analysis. The program could be used for the cost analysis of the spent fuel transportation as well.

What do reversible programs compute?

DEFF Research Database (Denmark)

Axelsen, Holger Bock; Glück, Robert

2011-01-01

Reversible computing is the study of computation models that exhibit both forward and backward determinism. Understanding the fundamental properties of such models is not only relevant for reversible programming, but has also been found important in other fields, e.g., bidirectional model...... transformation, program transformations such as inversion, and general static prediction of program properties. Historically, work on reversible computing has focussed on reversible simulations of irreversible computations. Here, we take the viewpoint that the property of reversibility itself should...... are not strictly classically universal, but that they support another notion of universality; we call this RTM-universality. Thus, even though the RTMs are sub-universal in the classical sense, they are powerful enough as to include a self-interpreter. Lifting this to other computation models, we propose r...

Mixtures Equation Pilot Program to Reduce Animal Testing

Science.gov (United States)

EPA is announcing the start of a pilot program to evaluate the usefulness and acceptability of a mathematical tool (the GHS Mixtures Equation), which is used in the Globally Harmonized System of Classification and Labeling of Chemicals (GHS).

Enhancing instruction in Fuels and Combustion Laboratory via a developed computer-assisted program for establishing efficient coal-diesel oil mixture (CDOM) fuel proportions

Energy Technology Data Exchange (ETDEWEB)

Maglaya, A.B. [La Salle University, Manila (Philippines). Dept. of Mechanical Engineering

2004-07-01

This paper discusses the relevance of digital computation in Fuels and Combustion Laboratory experiments used by the senior students of the Department of Mechanical Engineering, De La Salle University-Manila, Philippines. One of the students' experiments involved the determination of the most efficient CDOM fuel proportion as alternative fuel to diesel oil for steam generators and other industrial applications. Theoretical calculations show that it requires tedious and repetitive computations. A computer-assisted program was developed to lessen the time-consuming activities. The formulation of algorithms were based on the system of equations of the heat interaction between the CDOM fuel, combustion air and products of combustion and by applying the principles of mass and energy equations (or the First Law of Thermodynamics) for reacting systems were utilized. The developed computer-assisted program output verified alternative fuel selected through actual experimentation.

An algorithm of computing inhomogeneous differential equations for definite integrals

OpenAIRE

Nakayama, Hiromasa; Nishiyama, Kenta

2010-01-01

We give an algorithm to compute inhomogeneous differential equations for definite integrals with parameters. The algorithm is based on the integration algorithm for $D$-modules by Oaku. Main tool in the algorithm is the Gr\\"obner basis method in the ring of differential operators.

Functional programming for computer vision

Science.gov (United States)

Breuel, Thomas M.

1992-04-01

Functional programming is a style of programming that avoids the use of side effects (like assignment) and uses functions as first class data objects. Compared with imperative programs, functional programs can be parallelized better, and provide better encapsulation, type checking, and abstractions. This is important for building and integrating large vision software systems. In the past, efficiency has been an obstacle to the application of functional programming techniques in computationally intensive areas such as computer vision. We discuss and evaluate several 'functional' data structures for representing efficiently data structures and objects common in computer vision. In particular, we will address: automatic storage allocation and reclamation issues; abstraction of control structures; efficient sequential update of large data structures; representing images as functions; and object-oriented programming. Our experience suggests that functional techniques are feasible for high- performance vision systems, and that a functional approach simplifies the implementation and integration of vision systems greatly. Examples in C++ and SML are given.

Analysis and computation of the elastic wave equation with random coefficients

KAUST Repository

Motamed, Mohammad; Nobile, Fabio; Tempone, Raul

2015-01-01

We consider the stochastic initial-boundary value problem for the elastic wave equation with random coefficients and deterministic data. We propose a stochastic collocation method for computing statistical moments of the solution or statistics

A computer program incorporating Pitzer's equations for calculation of geochemical reactions in brines

Science.gov (United States)

Plummer, Niel; Parkhurst, D.L.; Fleming, G.W.; Dunkle, S.A.

1988-01-01

The program named PHRQPITZ is a computer code capable of making geochemical calculations in brines and other electrolyte solutions to high concentrations using the Pitzer virial-coefficient approach for activity-coefficient corrections. Reaction-modeling capabilities include calculation of (1) aqueous speciation and mineral-saturation index, (2) mineral solubility, (3) mixing and titration of aqueous solutions, (4) irreversible reactions and mineral water mass transfer, and (5) reaction path. The computed results for each aqueous solution include the osmotic coefficient, water activity , mineral saturation indices, mean activity coefficients, total activity coefficients, and scale-dependent values of pH, individual-ion activities and individual-ion activity coeffients , and scale-dependent values of pH, individual-ion activities and individual-ion activity coefficients. A data base of Pitzer interaction parameters is provided at 25 C for the system: Na-K-Mg-Ca-H-Cl-SO4-OH-HCO3-CO3-CO2-H2O, and extended to include largely untested literature data for Fe(II), Mn(II), Sr, Ba, Li, and Br with provision for calculations at temperatures other than 25C. An extensive literature review of published Pitzer interaction parameters for many inorganic salts is given. Also described is an interactive input code for PHRQPITZ called PITZINPT. (USGS)

Solving Coupled Gross--Pitaevskii Equations on a Cluster of PlayStation 3 Computers

Science.gov (United States)

Edwards, Mark; Heward, Jeffrey; Clark, C. W.

2009-05-01

At Georgia Southern University we have constructed an 8+1--node cluster of Sony PlayStation 3 (PS3) computers with the intention of using this computing resource to solve problems related to the behavior of ultra--cold atoms in general with a particular emphasis on studying bose--bose and bose--fermi mixtures confined in optical lattices. As a first project that uses this computing resource, we have implemented a parallel solver of the coupled time--dependent, one--dimensional Gross--Pitaevskii (TDGP) equations. These equations govern the behavior of dual-- species bosonic mixtures. We chose the split--operator/FFT to solve the coupled 1D TDGP equations. The fast Fourier transform component of this solver can be readily parallelized on the PS3 cpu known as the Cell Broadband Engine (CellBE). Each CellBE chip contains a single 64--bit PowerPC Processor Element known as the PPE and eight ``Synergistic Processor Element'' identified as the SPE's. We report on this algorithm and compare its performance to a non--parallel solver as applied to modeling evaporative cooling in dual--species bosonic mixtures.

Computer modeling of the dynamic processes in the Maryland University Training Reactor - (MUTR)

International Nuclear Information System (INIS)

White, Bernard H. IV; Ebert, David

1988-01-01

The simulator described in this paper models the behaviour of the Maryland University Training Reactor (MUTR). The reactor is a 250 kW, TRIGA reactor. The computer model is based on a system of five primary equations and eight auxiliary equations. The primary equations consist of the prompt jump approximation, a heat balance equation for the fuel and the moderator, and iodine and xenon buildup equations. For the comparison with the computer program, data from the reactor was acquired by using a personal computer (pc) which contained a Strawberry Tree data acquisition Card, connected to the reactor. The systems monitored by the pc were: two neutron detectors, fuel temperature, water temperature, three control rod positions and the period meter. The time differenced equations were programmed in the basic language. It has been shown by this paper, that the MUTR power rise from low power critical to high power, can be modelled by a relatively simple computer program. The program yields accurate agreement considering the simplicity of the program. The steady state error between the reactor and computer power is 4.4%. The difference in steady state temperatures, 112 deg. C and 117 deg. C, of the reactor and computer program, respectively, also yields a 4.5% error. Further fine tuning of the coefficients will yield higher accuracies

Computational time analysis of the numerical solution of 3D electrostatic Poisson's equation

Science.gov (United States)

Kamboh, Shakeel Ahmed; Labadin, Jane; Rigit, Andrew Ragai Henri; Ling, Tech Chaw; Amur, Khuda Bux; Chaudhary, Muhammad Tayyab

2015-05-01

3D Poisson's equation is solved numerically to simulate the electric potential in a prototype design of electrohydrodynamic (EHD) ion-drag micropump. Finite difference method (FDM) is employed to discretize the governing equation. The system of linear equations resulting from FDM is solved iteratively by using the sequential Jacobi (SJ) and sequential Gauss-Seidel (SGS) methods, simulation results are also compared to examine the difference between the results. The main objective was to analyze the computational time required by both the methods with respect to different grid sizes and parallelize the Jacobi method to reduce the computational time. In common, the SGS method is faster than the SJ method but the data parallelism of Jacobi method may produce good speedup over SGS method. In this study, the feasibility of using parallel Jacobi (PJ) method is attempted in relation to SGS method. MATLAB Parallel/Distributed computing environment is used and a parallel code for SJ method is implemented. It was found that for small grid size the SGS method remains dominant over SJ method and PJ method while for large grid size both the sequential methods may take nearly too much processing time to converge. Yet, the PJ method reduces computational time to some extent for large grid sizes.

Symbolic computation and solitons of the nonlinear Schroedinger equation in inhomogeneous optical fiber media

International Nuclear Information System (INIS)

Li Biao; Chen Yong

2007-01-01

In this paper, the inhomogeneous nonlinear Schroedinger equation with the loss/gain and the frequency chirping is investigated. With the help of symbolic computation, three families of exact analytical solutions are presented by employing the extended projective Riccati equation method. From our results, many previous known results of nonlinear Schroedinger equation obtained by some authors can be recovered by means of some suitable selections of the arbitrary functions and arbitrary constants. Of optical and physical interests, soliton propagation and soliton interaction are discussed and simulated by computer, which include snake-soliton propagation and snake-solitons interaction, boomerang-like soliton propagation and boomerang-like solitons interaction, dispersion managed (DM) bright (dark) soliton propagation and DM solitons interaction

SCOTCH: a program for solution of the one-dimensional, two-group, space-time neutron diffusion equations with temperature feedback of multi-channel fluid dynamics for HTGR cores

International Nuclear Information System (INIS)

Ezaki, Masahiro; Mitake, Susumu; Ozawa, Tamotsu

1979-06-01

The SCOTCH program solves the one-dimensional (R or Z), two-group reactor kinetics equations with multi-channel temperature transients and fluid dynamics. Sub-program SCOTCH-RX simulates the space-time neutron diffusion in radial direction, and sub-program SCOTCH-AX simulates the same in axial direction. The program has about 8,000 steps of FORTRAN statement and requires about 102 kilo-words of computer memory. (author)

A homotopy method for solving Riccati equations on a shared memory parallel computer

International Nuclear Information System (INIS)

Zigic, D.; Watson, L.T.; Collins, E.G. Jr.; Davis, L.D.

1993-01-01

Although there are numerous algorithms for solving Riccati equations, there still remains a need for algorithms which can operate efficiently on large problems and on parallel machines. This paper gives a new homotopy-based algorithm for solving Riccati equations on a shared memory parallel computer. The central part of the algorithm is the computation of the kernel of the Jacobian matrix, which is essential for the corrector iterations along the homotopy zero curve. Using a Schur decomposition the tensor product structure of various matrices can be efficiently exploited. The algorithm allows for efficient parallelization on shared memory machines

Debugging a high performance computing program

Science.gov (United States)

Gooding, Thomas M.

2013-08-20

Methods, apparatus, and computer program products are disclosed for debugging a high performance computing program by gathering lists of addresses of calling instructions for a plurality of threads of execution of the program, assigning the threads to groups in dependence upon the addresses, and displaying the groups to identify defective threads.

Computational partial differential equations using Matlab

CERN Document Server

Li, Jichun

2008-01-01

Brief Overview of Partial Differential Equations The parabolic equations The wave equations The elliptic equations Differential equations in broader areasA quick review of numerical methods for PDEsFinite Difference Methods for Parabolic Equations Introduction Theoretical issues: stability, consistence, and convergence 1-D parabolic equations2-D and 3-D parabolic equationsNumerical examples with MATLAB codesFinite Difference Methods for Hyperbolic Equations IntroductionSome basic difference schemes Dissipation and dispersion errors Extensions to conservation lawsThe second-order hyperbolic PDE

TRIGASIM: A computer program to simulate a TRIGA Mark I Reactor

International Nuclear Information System (INIS)

Ruby, Lawrence

1992-01-01

A Fortran-77 computer program has been written which simulates the operation of a TRIGA Mark I Reactor. The 'operator' has options at 1-second intervals, of raising rods, lowering rods, maintaining rods steady, dropping a rod, or scramming the reactor. Results are printed to the screen, and to 2 output files - a tabular record and a logarithmic plot of the power. The Point Kinetic Equations are programmed with 6 delayed groups, quasi-static power feedback, and forward differencing. A pulsing option is available, with simulation which employs the Fuchs Model. A pulse-tail model has been devised to simulate behavior for a few minutes following a pulse. Both graphic and tabular output are also available for the pulses. (author)

Some Matrix Iterations for Computing Generalized Inverses and Balancing Chemical Equations

OpenAIRE

Soleimani, Farahnaz; Stanimirovi´c, Predrag; Soleymani, Fazlollah

2015-01-01

An application of iterative methods for computing the Moore–Penrose inverse in balancing chemical equations is considered. With the aim to illustrate proposed algorithms, an improved high order hyper-power matrix iterative method for computing generalized inverses is introduced and applied. The improvements of the hyper-power iterative scheme are based on its proper factorization, as well as on the possibility to accelerate the iterations in the initial phase of the convergence. Although the ...

Solution of Fractional Order System of Bagley-Torvik Equation Using Evolutionary Computational Intelligence

Directory of Open Access Journals (Sweden)

Muhammad Asif Zahoor Raja

2011-01-01

Full Text Available A stochastic technique has been developed for the solution of fractional order system represented by Bagley-Torvik equation. The mathematical model of the equation was developed with the help of feed-forward artificial neural networks. The training of the networks was made with evolutionary computational intelligence based on genetic algorithm hybrid with pattern search technique. Designed scheme was successfully applied to different forms of the equation. Results are compared with standard approximate analytic, stochastic numerical solvers and exact solutions.

Effective computation of stochastic protein kinetic equation by reducing stiffness via variable transformation

Energy Technology Data Exchange (ETDEWEB)

Wang, Lijin, E-mail: ljwang@ucas.ac.cn [School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing 100049 (China)

2016-06-08

The stochastic protein kinetic equations can be stiff for certain parameters, which makes their numerical simulation rely on very small time step sizes, resulting in large computational cost and accumulated round-off errors. For such situation, we provide a method of reducing stiffness of the stochastic protein kinetic equation by means of a kind of variable transformation. Theoretical and numerical analysis show effectiveness of this method. Its generalization to a more general class of stochastic differential equation models is also discussed.

TEMPEST: A three-dimensional time-dependent computer program for hydrothermal analysis: Volume 1, Numerical methods and input instructions

International Nuclear Information System (INIS)

Trent, D.S.; Eyler, L.L.; Budden, M.J.

1983-09-01

This document describes the numerical methods, current capabilities, and the use of the TEMPEST (Version L, MOD 2) computer program. TEMPEST is a transient, three-dimensional, hydrothermal computer program that is designed to analyze a broad range of coupled fluid dynamic and heat transfer systems of particular interest to the Fast Breeder Reactor thermal-hydraulic design community. The full three-dimensional, time-dependent equations of motion, continuity, and heat transport are solved for either laminar or turbulent fluid flow, including heat diffusion and generation in both solid and liquid materials. 10 refs., 22 figs., 2 tabs

SciCADE 95: International conference on scientific computation and differential equations

Energy Technology Data Exchange (ETDEWEB)

NONE

1995-12-31

This report consists of abstracts from the conference. Topics include algorithms, computer codes, and numerical solutions for differential equations. Linear and nonlinear as well as boundary-value and initial-value problems are covered. Various applications of these problems are also included.

Solving the Fokker-Planck equation on a massively parallel computer

International Nuclear Information System (INIS)

Mirin, A.A.

1990-01-01

The Fokker-Planck package FPPAC had been converted to the Connection Machine 2 (CM2). For fine mesh cases the CM2 outperforms the Cray-2 when it comes to time-integrating the difference equations. For long Legendre expansions the CM2 is also faster at computing the Fokker-Planck coefficients. 3 refs

Finite Volume Based Computer Program for Ground Source Heat Pump System

Energy Technology Data Exchange (ETDEWEB)

Menart, James A. [Wright State University

2013-02-22

This report is a compilation of the work that has been done on the grant DE-EE0002805 entitled ?Finite Volume Based Computer Program for Ground Source Heat Pump Systems.? The goal of this project was to develop a detailed computer simulation tool for GSHP (ground source heat pump) heating and cooling systems. Two such tools were developed as part of this DOE (Department of Energy) grant; the first is a two-dimensional computer program called GEO2D and the second is a three-dimensional computer program called GEO3D. Both of these simulation tools provide an extensive array of results to the user. A unique aspect of both these simulation tools is the complete temperature profile information calculated and presented. Complete temperature profiles throughout the ground, casing, tube wall, and fluid are provided as a function of time. The fluid temperatures from and to the heat pump, as a function of time, are also provided. In addition to temperature information, detailed heat rate information at several locations as a function of time is determined. Heat rates between the heat pump and the building indoor environment, between the working fluid and the heat pump, and between the working fluid and the ground are computed. The heat rates between the ground and the working fluid are calculated as a function time and position along the ground loop. The heating and cooling loads of the building being fitted with a GSHP are determined with the computer program developed by DOE called ENERGYPLUS. Lastly COP (coefficient of performance) results as a function of time are provided. Both the two-dimensional and three-dimensional computer programs developed as part of this work are based upon a detailed finite volume solution of the energy equation for the ground and ground loop. Real heat pump characteristics are entered into the program and used to model the heat pump performance. Thus these computer tools simulate the coupled performance of the ground loop and the heat pump

Application of viscoplastic constitutive equations in finite element programs

International Nuclear Information System (INIS)

Hornberger, K.; Stamm, H.

1987-04-01

The general mathematical formulation of frequently used viscoplastic constitutive equations is explained and Robinson's model is discussed in more detail. The implementation of viscoplastic constitutive equations into Finite Element programs (such as ABAQUS) is described using Robinson's model as an example. For the numerical integration both an explicit (explicit Euler) and an implicit (generalized midpoint rule) integration scheme is utilized in combination with a time step control strategy. In the implicit integration scheme, convergence in solving a system of nonlinear algebraic equation is improved introducing a projection method. The efficiency of the implemented procedures is demonstrated for different homogeneous load cases as well as for creep loading and strain controlled cyclic loading of a perforated plate. (orig./HP) [de

Computational Acoustics: Computational PDEs, Pseudodifferential Equations, Path Integrals, and All That Jazz

Science.gov (United States)

Fishman, Louis

2000-11-01

The role of mathematical modeling in the physical sciences will be briefly addressed. Examples will focus on computational acoustics, with applications to underwater sound propagation, electromagnetic modeling, optics, and seismic inversion. Direct and inverse wave propagation problems in both the time and frequency domains will be considered. Focusing on fixed-frequency (elliptic) wave propagation problems, the usual, two-way, partial differential equation formulation will be exactly reformulated, in a well-posed manner, as a one-way (marching) problem. This is advantageous for both direct and inverse considerations, as well as stochastic modeling problems. The reformulation will require the introduction of pseudodifferential operators and their accompanying phase space analysis (calculus), in addition to path integral representations for the fundamental solutions and their subsequent computational algorithms. Unlike the more traditional, purely numerical applications of, for example, finite-difference and finite-element methods, this approach, in effect, writes the exact, or, more generally, the asymptotically correct, answer as a functional integral and, subsequently, computes it directly. The overall computational philosophy is to combine analysis, asymptotics, and numerical methods to attack complicated, real-world problems. Exact and asymptotic analysis will stress the complementary nature of the direct and inverse formulations, as well as indicating the explicit structural connections between the time- and frequency-domain solutions.

The psychology of computer programming

CERN Document Server

Weinberg, Gerald Marvin

1998-01-01

This landmark 1971 classic is reprinted with a new preface, chapter-by-chapter commentary, and straight-from-the-heart observations on topics that affect the professional life of programmers. Long regarded as one of the first books to pioneer a people-oriented approach to computing, The Psychology of Computer Programming endures as a penetrating analysis of the intelligence, skill, teamwork, and problem-solving power of the computer programmer. Finding the chapters strikingly relevant to today's issues in programming, Gerald M. Weinberg adds new insights and highlights the similarities and differences between now and then. Using a conversational style that invites the reader to join him, Weinberg reunites with some of his most insightful writings on the human side of software engineering. Topics include egoless programming, intelligence, psychological measurement, personality factors, motivation, training, social problems on large projects, problem-solving ability, programming language design, team formati...

The Julia programming language: the future of scientific computing

Science.gov (United States)

Gibson, John

2017-11-01

Julia is an innovative new open-source programming language for high-level, high-performance numerical computing. Julia combines the general-purpose breadth and extensibility of Python, the ease-of-use and numeric focus of Matlab, the speed of C and Fortran, and the metaprogramming power of Lisp. Julia uses type inference and just-in-time compilation to compile high-level user code to machine code on the fly. A rich set of numeric types and extensive numerical libraries are built-in. As a result, Julia is competitive with Matlab for interactive graphical exploration and with C and Fortran for high-performance computing. This talk interactively demonstrates Julia's numerical features and benchmarks Julia against C, C++, Fortran, Matlab, and Python on a spectral time-stepping algorithm for a 1d nonlinear partial differential equation. The Julia code is nearly as compact as Matlab and nearly as fast as Fortran. This material is based upon work supported by the National Science Foundation under Grant No. 1554149.

FLANGE: a computer program for the analysis of flanged joints with ring-type gaskets

International Nuclear Information System (INIS)

Rodabaugh, E.C.; O'Hara, F.M. Jr.; Moore, S.E.

1976-01-01

The computer program FLANGE was written to calculate not only the stresses due to moment loads on the flange ring but also stresses due to internal pressure; stresses due to a temperature difference between the hub and ring; and stresses due to the variations in bolt load that result from pressure, hub-ring temperature gradient, and/or bolt-ring temperature difference. The program FLANGE is applicable to tapered-hub, straight, and blind flanges. The analysis method is based on the differential equations for thin plates and shells. The stresses due to moment loading calculated by the two methods are essentially identical for identical boundary conditions. A description of the general model of flanges used in the theoretical development of the computer code is provided. The actual mathematical expressions for calculating stresses and displacements due to moment and pressure loads are derived. Example calculations, listings, and flowcharts of the program and its subroutines are included as appendices

Exploring Poetry through Interactive Computer Programs.

Science.gov (United States)

Nimchinsky, Howard; Camp, Jocelyn

The goal of a project was to design, test, and evaluate several computer programs that allow students in introductory literature and poetry courses to explore a poem in detail and, through a dialogue with the program, to develop their own interpretation of it. Computer programs were completed on poems by Robert Frost and W.H. Auden. Both programs…

Wing-Body Aeroelasticity Using Finite-Difference Fluid/Finite-Element Structural Equations on Parallel Computers

Science.gov (United States)

Byun, Chansup; Guruswamy, Guru P.; Kutler, Paul (Technical Monitor)

1994-01-01

In recent years significant advances have been made for parallel computers in both hardware and software. Now parallel computers have become viable tools in computational mechanics. Many application codes developed on conventional computers have been modified to benefit from parallel computers. Significant speedups in some areas have been achieved by parallel computations. For single-discipline use of both fluid dynamics and structural dynamics, computations have been made on wing-body configurations using parallel computers. However, only a limited amount of work has been completed in combining these two disciplines for multidisciplinary applications. The prime reason is the increased level of complication associated with a multidisciplinary approach. In this work, procedures to compute aeroelasticity on parallel computers using direct coupling of fluid and structural equations will be investigated for wing-body configurations. The parallel computer selected for computations is an Intel iPSC/860 computer which is a distributed-memory, multiple-instruction, multiple data (MIMD) computer with 128 processors. In this study, the computational efficiency issues of parallel integration of both fluid and structural equations will be investigated in detail. The fluid and structural domains will be modeled using finite-difference and finite-element approaches, respectively. Results from the parallel computer will be compared with those from the conventional computers using a single processor. This study will provide an efficient computational tool for the aeroelastic analysis of wing-body structures on MIMD type parallel computers.

Generalized fish life-cycle poplulation model and computer program

International Nuclear Information System (INIS)

DeAngelis, D.L.; Van Winkle, W.; Christensen, S.W.; Blum, S.R.; Kirk, B.L.; Rust, B.W.; Ross, C.

1978-03-01

A generalized fish life-cycle population model and computer program have been prepared to evaluate the long-term effect of changes in mortality in age class 0. The general question concerns what happens to a fishery when density-independent sources of mortality are introduced that act on age class 0, particularly entrainment and impingement at power plants. This paper discusses the model formulation and computer program, including sample results. The population model consists of a system of difference equations involving age-dependent fecundity and survival. The fecundity for each age class is assumed to be a function of both the fraction of females sexually mature and the weight of females as they enter each age class. Natural mortality for age classes 1 and older is assumed to be independent of population size. Fishing mortality is assumed to vary with the number and weight of fish available to the fishery. Age class 0 is divided into six life stages. The probability of survival for age class 0 is estimated considering both density-independent mortality (natural and power plant) and density-dependent mortality for each life stage. Two types of density-dependent mortality are included. These are cannibalism of each life stage by older age classes and intra-life-stage competition

RECON: a computer program for analyzing repository economics. Documentation and user's manual. Revision 1

International Nuclear Information System (INIS)

Clark, L.L.; Schutz, M.E.; Luksic, A.T.

1985-07-01

From 1981 through 1984 the Pacific Northwest Laboratory has been developing a computer model named RECON to calculate repository costs from parametric data input. The objective of the program has been to develop the capability to evaluate the effect on costs of changes in repository design parameters and operating scenario assumptions. This report documents the development of the model through September of 1984. Included in the report are: (1) descriptions of model development and the underlying equations, assumptions and definitions; (2) descriptions of data input using either card images or an interactive data input program; and (3) detailed listings of the program and definitions of program variables. Cost estimates generated using the model have been verified against independent estimates and good agreement has been obtained. 2 refs

The Y2K program for scientific-analysis computer programs at AECL

International Nuclear Information System (INIS)

Popovic, J.; Gaver, C.; Chapman, D.

1999-01-01

The evaluation of scientific-analysis computer programs for year-2000 compliance is part of AECL' s year-2000 (Y2K) initiative, which addresses both the infrastructure systems at AECL and AECL's products and services. This paper describes the Y2K-compliance program for scientific-analysis computer codes. This program involves the integrated evaluation of the computer hardware, middleware, and third-party software in addition to the scientific codes developed in-house. The project involves several steps: the assessment of the scientific computer programs for Y2K compliance, performing any required corrective actions, porting the programs to Y2K-compliant platforms, and verification of the programs after porting. Some programs or program versions, deemed no longer required in the year 2000 and beyond, will be retired and archived. (author)

The Y2K program for scientific-analysis computer programs at AECL

International Nuclear Information System (INIS)

Popovic, J.; Gaver, C.; Chapman, D.

1999-01-01

The evaluation of scientific analysis computer programs for year-2000 compliance is part of AECL's year-2000 (Y2K) initiative, which addresses both the infrastructure systems at AECL and AECL's products and services. This paper describes the Y2K-compliance program for scientific-analysis computer codes. This program involves the integrated evaluation of the computer hardware, middleware, and third-party software in addition to the scientific codes developed in-house. The project involves several steps: the assessment of the scientific computer programs for Y2K compliance, performing any required corrective actions, porting the programs to Y2K-compliant platforms, and verification of the programs after porting. Some programs or program versions, deemed no longer required in the year 2000 and beyond, will be retired and archived. (author)

Numerical methods for differential equations and applications

International Nuclear Information System (INIS)

Ixaru, L.G.

1984-01-01

This book is addressed to persons who, without being professionals in applied mathematics, are often faced with the problem of numerically solving differential equations. In each of the first three chapters a definite class of methods is discussed for the solution of the initial value problem for ordinary differential equations: multistep methods; one-step methods; and piecewise perturbation methods. The fourth chapter is mainly focussed on the boundary value problems for linear second-order equations, with a section devoted to the Schroedinger equation. In the fifth chapter the eigenvalue problem for the radial Schroedinger equation is solved in several ways, with computer programs included. (Auth.)

Parallel computing solution of Boltzmann neutron transport equation

International Nuclear Information System (INIS)

Ansah-Narh, T.

2010-01-01

The focus of the research was on developing parallel computing algorithm for solving Eigen-values of the Boltzmam Neutron Transport Equation (BNTE) in a slab geometry using multi-grid approach. In response to the problem of slow execution of serial computing when solving large problems, such as BNTE, the study was focused on the design of parallel computing systems which was an evolution of serial computing that used multiple processing elements simultaneously to solve complex physical and mathematical problems. Finite element method (FEM) was used for the spatial discretization scheme, while angular discretization was accomplished by expanding the angular dependence in terms of Legendre polynomials. The eigenvalues representing the multiplication factors in the BNTE were determined by the power method. MATLAB Compiler Version 4.1 (R2009a) was used to compile the MATLAB codes of BNTE. The implemented parallel algorithms were enabled with matlabpool, a Parallel Computing Toolbox function. The option UseParallel was set to 'always' and the default value of the option was 'never'. When those conditions held, the solvers computed estimated gradients in parallel. The parallel computing system was used to handle all the bottlenecks in the matrix generated from the finite element scheme and each domain of the power method generated. The parallel algorithm was implemented on a Symmetric Multi Processor (SMP) cluster machine, which had Intel 32 bit quad-core x 86 processors. Convergence rates and timings for the algorithm on the SMP cluster machine were obtained. Numerical experiments indicated the designed parallel algorithm could reach perfect speedup and had good stability and scalability. (au)

NORIA-SP: A finite element computer program for analyzing liquid water transport in porous media

International Nuclear Information System (INIS)

Hopkins, P.L.; Eaton, R.R.; Bixler, N.E.

1991-12-01

A family of finite element computer programs has been developed at Sandia National Laboratories (SNL) most recently, NORIA-SP. The original NORIA code solves a total of four transport equations simultaneously: liquid water, water vapor, air, and energy. Consequently, use of NORIA is computer-intensive. Since many of the applications for which NORIA is used are isothermal, we decided to ''strip'' the original four-equation version, leaving only the liquid water equation. This single-phase version is NORIA-SP. The primary intent of this document is to provide the user of NORIA-SP an accurate user's manual. Consequently, the reader should refer to the NORIA manual if additional detail is required regarding the equation development and finite element methods used. The single-equation version of the NORIA code (NORIA-SP) has been used most frequently for analyzing various hydrological scenarios for the potential underground nuclear waste repository at Yucca Mountain in western Nevada. These analyses are generally performed assuming a composite model to represent the fractured geologic media. In this model the material characteristics of the matrix and the fractures are area weighted to obtain equivalent material properties. Pressure equilibrium between the matrix and fractures is assumed so a single conservation equation can be solved. NORIA-SP is structured to accommodate the composite model. The equations for water velocities in both the rock matrix and the fractures are presented. To use the code for problems involving a single, nonfractured porous material, the user can simply set the area of the fractures to zero

Application of the reactor kinetics equations to the reactor safety analysis

International Nuclear Information System (INIS)

Sdouz, G.

1976-01-01

The reactor kinetics equations which can be solved by the computer program AIREK-III are used to describe the behavior of fast reactivity transients. By supplementing this computer program it was possible to solve additional safety problems, e.g. the course of reactor excursions induced by any form of reactivity input, the control of reactivity input as a function of a threshold-energy and the computation of produced energy. (author)

A computational procedure for finding multiple solutions of convective heat transfer equations

International Nuclear Information System (INIS)

Mishra, S; DebRoy, T

2005-01-01

In recent years numerical solutions of the convective heat transfer equations have provided significant insight into the complex materials processing operations. However, these computational methods suffer from two major shortcomings. First, these procedures are designed to calculate temperature fields and cooling rates as output and the unidirectional structure of these solutions preclude specification of these variables as input even when their desired values are known. Second, and more important, these procedures cannot determine multiple pathways or multiple sets of input variables to achieve a particular output from the convective heat transfer equations. Here we propose a new method that overcomes the aforementioned shortcomings of the commonly used solutions of the convective heat transfer equations. The procedure combines the conventional numerical solution methods with a real number based genetic algorithm (GA) to achieve bi-directionality, i.e. the ability to calculate the required input variables to achieve a specific output such as temperature field or cooling rate. More important, the ability of the GA to find a population of solutions enables this procedure to search for and find multiple sets of input variables, all of which can lead to the desired specific output. The proposed computational procedure has been applied to convective heat transfer in a liquid layer locally heated on its free surface by an electric arc, where various sets of input variables are computed to achieve a specific fusion zone geometry defined by an equilibrium temperature. Good agreement is achieved between the model predictions and the independent experimental results, indicating significant promise for the application of this procedure in finding multiple solutions of convective heat transfer equations

Solving the neutron diffusion equation on combinatorial geometry computational cells for reactor physics calculations

International Nuclear Information System (INIS)

Azmy, Y. Y.

2004-01-01

An approach is developed for solving the neutron diffusion equation on combinatorial geometry computational cells, that is computational cells composed by combinatorial operations involving simple-shaped component cells. The only constraint on the component cells from which the combinatorial cells are assembled is that they possess a legitimate discretization of the underlying diffusion equation. We use the Finite Difference (FD) approximation of the x, y-geometry diffusion equation in this work. Performing the same combinatorial operations involved in composing the combinatorial cell on these discrete-variable equations yields equations that employ new discrete variables defined only on the combinatorial cell's volume and faces. The only approximation involved in this process, beyond the truncation error committed in discretizing the diffusion equation over each component cell, is a consistent-order Legendre series expansion. Preliminary results for simple configurations establish the accuracy of the solution to the combinatorial geometry solution compared to straight FD as the system dimensions decrease. Furthermore numerical results validate the consistent Legendre-series expansion order by illustrating the second order accuracy of the combinatorial geometry solution, the same as standard FD. Nevertheless the magnitude of the error for the new approach is larger than FD's since it incorporates the additional truncated series approximation. (authors)

GEOSURF: a computer program for modeling adsorption on mineral surfaces from aqueous solution

Science.gov (United States)

Sahai, Nita; Sverjensky, Dimitri A.

1998-11-01

A new program, GEOSURF, has been developed for calculating aqueous and surface speciation consistent with the triple-layer model of surface complexation. GEOSURF is an extension of the original programs MINEQL, MICROQL and HYDRAQL. We present, here, the basic algorithm of GEOSURF along with a description of the new features implemented. GEOSURF is linked to internally consistent data bases for surface species (SURFK.DAT) and for aqueous species (AQSOL.DAT). SURFK.DAT contains properties of minerals such as site densities, and equilibrium constants for adsorption of aqueous protons and electrolyte ions on a variety of oxides and hydroxides. The Helgeson, Kirkham and Flowers version of the extended Debye-Huckel Equation for 1:1 electrolytes is implemented for calculating aqueous activity coefficients. This permits the calculation of speciation at ionic strengths greater than 0.5 M. The activity of water is computed explicitly from the osmotic coefficient of the solution, and the total amount of electrolyte cation (or anion) is adjusted to satisfy the electroneutrality condition. Finally, the use of standard symbols for chemical species rather than species identification numbers is included to facilitate use of the program. One of the main limitations of GEOSURF is that aqueous and surface speciation can only be calculated at fixed pH and at fixed concentration of total adsorbate. Thus, the program cannot perform reaction-path calculations: it cannot determine whether or not a solution is over- or under-saturated with respect to one or more solid phases. To check the proper running of GEOSURF, we have compared results generated by GEOSURF with those from two other programs, HYDRAQL and EQ3. The Davies equation and the "bdot" equation, respectively, are used in the latter two programs for calculating aqueous activity coefficients. An example of the model fit to experimental data for rutile in 0.001 M-2.0 M NaNO 3 is included.

Computer Program for Calculation of Complex Chemical Equilibrium Compositions, Rocket Performance, Incident and Reflected Shocks, and Chapman-Jouguet Detonations. Interim Revision, March 1976

Science.gov (United States)

Gordon, S.; Mcbride, B. J.

1976-01-01

A detailed description of the equations and computer program for computations involving chemical equilibria in complex systems is given. A free-energy minimization technique is used. The program permits calculations such as (1) chemical equilibrium for assigned thermodynamic states (T,P), (H,P), (S,P), (T,V), (U,V), or (S,V), (2) theoretical rocket performance for both equilibrium and frozen compositions during expansion, (3) incident and reflected shock properties, and (4) Chapman-Jouguet detonation properties. The program considers condensed species as well as gaseous species.

Computationally efficient statistical differential equation modeling using homogenization

Science.gov (United States)

Hooten, Mevin B.; Garlick, Martha J.; Powell, James A.

2013-01-01

Statistical models using partial differential equations (PDEs) to describe dynamically evolving natural systems are appearing in the scientific literature with some regularity in recent years. Often such studies seek to characterize the dynamics of temporal or spatio-temporal phenomena such as invasive species, consumer-resource interactions, community evolution, and resource selection. Specifically, in the spatial setting, data are often available at varying spatial and temporal scales. Additionally, the necessary numerical integration of a PDE may be computationally infeasible over the spatial support of interest. We present an approach to impose computationally advantageous changes of support in statistical implementations of PDE models and demonstrate its utility through simulation using a form of PDE known as “ecological diffusion.” We also apply a statistical ecological diffusion model to a data set involving the spread of mountain pine beetle (Dendroctonus ponderosae) in Idaho, USA.

An Evaluation of Kernel Equating: Parallel Equating with Classical Methods in the SAT Subject Tests[TM] Program. Research Report. ETS RR-09-06

Science.gov (United States)

Grant, Mary C.; Zhang, Lilly; Damiano, Michele

2009-01-01

This study investigated kernel equating methods by comparing these methods to operational equatings for two tests in the SAT Subject Tests[TM] program. GENASYS (ETS, 2007) was used for all equating methods and scaled score kernel equating results were compared to Tucker, Levine observed score, chained linear, and chained equipercentile equating…

A Model of Computation for Bit-Level Concurrent Computing and Programming: APEC

Science.gov (United States)

Ajiro, Takashi; Tsuchida, Kensei

A concurrent model of computation and a language based on the model for bit-level operation are useful for developing asynchronous and concurrent programs compositionally, which frequently use bit-level operations. Some examples are programs for video games, hardware emulation (including virtual machines), and signal processing. However, few models and languages are optimized and oriented to bit-level concurrent computation. We previously developed a visual programming language called A-BITS for bit-level concurrent programming. The language is based on a dataflow-like model that computes using processes that provide serial bit-level operations and FIFO buffers connected to them. It can express bit-level computation naturally and develop compositionally. We then devised a concurrent computation model called APEC (Asynchronous Program Elements Connection) for bit-level concurrent computation. This model enables precise and formal expression of the process of computation, and a notion of primitive program elements for controlling and operating can be expressed synthetically. Specifically, the model is based on a notion of uniform primitive processes, called primitives, that have three terminals and four ordered rules at most, as well as on bidirectional communication using vehicles called carriers. A new notion is that a carrier moving between two terminals can briefly express some kinds of computation such as synchronization and bidirectional communication. The model's properties make it most applicable to bit-level computation compositionally, since the uniform computation elements are enough to develop components that have practical functionality. Through future application of the model, our research may enable further research on a base model of fine-grain parallel computer architecture, since the model is suitable for expressing massive concurrency by a network of primitives.

Computational Fluid Dynamics

International Nuclear Information System (INIS)

Myeong, Hyeon Guk

1999-06-01

This book deals with computational fluid dynamics with basic and history of numerical fluid dynamics, introduction of finite volume method using one-dimensional heat conduction equation, solution of two-dimensional heat conduction equation, solution of Navier-Stokes equation, fluid with heat transport, turbulent flow and turbulent model, Navier-Stokes solution by generalized coordinate system such as coordinate conversion, conversion of basic equation, program and example of calculation, application of abnormal problem and high speed solution of numerical fluid dynamics.

A solution of the thermal neutron diffusion equation for a two-region cyclindrical system program for ODRA-1305 computer

International Nuclear Information System (INIS)

Drozdowicz, K.; Woznicka, U.

1982-01-01

The program in FORTRAN for the ODRA-1305 computer is described. The dependence of the decay constant of the thermal neutron flux upon the dimensions of the two-region concentric cylindrical system is the result of the program. The solution (with a constant neutron flux in the inner medium assumed) is generally obtained in the one-group diffusion approximation by the method of the perturbation calculation. However, the energy distribution of the thermal neutron flux and the diffusion cooling are taken into account. The program is written for the case when the outer medium is hydrogenous. The listing of the program and an example of calculation results are included. (author)

Symbolic Computations and Exact and Explicit Solutions of Some Nonlinear Evolution Equations in Mathematical Physics

International Nuclear Information System (INIS)

Oezis, Turgut; Aslan, Imail

2009-01-01

With the aid of symbolic computation system Mathematica, several explicit solutions for Fisher's equation and CKdV equation are constructed by utilizing an auxiliary equation method, the so called G'/G-expansion method, where the new and more general forms of solutions are also constructed. When the parameters are taken as special values, the previously known solutions are recovered. (general)

USERDA computer program summaries. Numbers 177--239

International Nuclear Information System (INIS)

1975-10-01

Since 1960 the Argonne Code Center has served as a U. S. Atomic Energy Commission information center for computer programs developed and used primarily for the solution of problems in nuclear physics, reactor design, reactor engineering and operation. The Center, through a network of registered installations, collects, validates, maintains, and distributes a library of these computer programs and publishes a compilation of abstracts describing them. In 1972 the scope of the Center's activities was officially expanded to include computer programs developed in all of the U. S. Atomic Energy Commission program areas and the compilation and publication of this report. The Computer Program Summary report contains summaries of computer programs at the specification stage, under development, being checked out, in use, or available at ERDA offices, laboratories, and contractor installations. Programs are divided into the following categories: cross section and resonance integral calculations; spectrum calculations, generation of group constants, lattice and cell problems; static design studies; depletion, fuel management, cost analysis, and reactor economics; space-independent kinetics; space--time kinetics, coupled neutronics--hydrodynamics--thermodynamics and excursion simulations; radiological safety, hazard and accident analysis; heat transfer and fluid flow; deformation and stress distribution computations, structural analysis and engineering design studies; gamma heating and shield design programs; reactor systems analysis; data preparation; data management; subsidiary calculations; experimental data processing; general mathematical and computing system routines; materials; environmental and earth sciences; space sciences; electronics and engineering equipment; chemistry; particle accelerators and high-voltage machines; physics; controlled thermonuclear research; biology and medicine; and data

One-dimensional unstable eigenfunction and manifold computations in delay differential equations

International Nuclear Information System (INIS)

Green, Kirk; Krauskopf, Bernd; Engelborghs, Koen

2004-01-01

In this paper we present a new numerical technique for computing the unstable eigenfunctions of a saddle periodic orbit in a delay differential equation. This is used to obtain the necessary starting data for an established algorithm for computing one-dimensional (1D) unstable manifolds of an associated saddle fixed point of a suitable Poincare map. To illustrate our method, we investigate an intermittent transition to chaos in a delay system describing a semiconductor laser subject to phase-conjugate feedback

Multi-level programming paradigm for extreme computing

International Nuclear Information System (INIS)

Petiton, S.; Sato, M.; Emad, N.; Calvin, C.; Tsuji, M.; Dandouna, M.

2013-01-01

In order to propose a framework and programming paradigms for post peta-scale computing, on the road to exa-scale computing and beyond, we introduced new languages, associated with a hierarchical multi-level programming paradigm, allowing scientific end-users and developers to program highly hierarchical architectures designed for extreme computing. In this paper, we explain the interest of such hierarchical multi-level programming paradigm for extreme computing and its well adaptation to several large computational science applications, such as for linear algebra solvers used for reactor core physic. We describe the YML language and framework allowing describing graphs of parallel components, which may be developed using PGAS-like language such as XMP, scheduled and computed on supercomputers. Then, we propose experimentations on supercomputers (such as the 'K' and 'Hooper' ones) of the hybrid method MERAM (Multiple Explicitly Restarted Arnoldi Method) as a case study for iterative methods manipulating sparse matrices, and the block Gauss-Jordan method as a case study for direct method manipulating dense matrices. We conclude proposing evolutions for this programming paradigm. (authors)

HAPS, a Handy Analog Programming System

DEFF Research Database (Denmark)

Højberg, Kristian Søe

1975-01-01

HAPS (Hybrid Analog Programming System) is an analog compiler that can be run on a minicomputer in an interactive mode. Essentially HAPS is written in FORTRAN. The equations to be programmed for an ana log computer are read in by using a FORTRAN-like notation. The input must contain maximum...... and emphasizes the limitations HAPS puts on equation structure, types of computing circuit, scaling, and static testing....

SARTEMP2 - A computer program to calculate power and temperatures in a transport flask during a criticality accident

International Nuclear Information System (INIS)

Shaw, P.M.

1983-04-01

The computer code SARTEMP2, an extended version of the original SARTEMP program, which calculates the power and temperatures in a transport flask during a hypothetical criticality accident is described. The accident arises, it is assumed, during the refilling of the flask with water, bringing the system to delayed critical. As the water level continues to rise, reactivity is added causing the power to rise, and thus temperatures in the fuel, clad and water to increase. The point kinetics equations are coupled to the one-dimensional heat conduction equation. The model used, the method of solution of the equations and the input data required are given. (author)

Development of a computer program for drop time and impact velocity of the rod cluster control assembly

International Nuclear Information System (INIS)

Choi, K.-S.; Yim, J.-S.; Kim, I.-K.; Kim, K.-T.

1993-01-01

In PWR the rod cluster control assembly (RCCA) for shutdown is released upon the action of the control drive mechanism and falls down through the guide thimble by its weight. Drop time and impact velocity of the RCCA are two key parameters with respect to reactivity insertion time and the mechanical integrity of fuel assembly. Therefore, the precise control of the drop time and impact velocity is prerequisite to modifying the existing design features of the RCCA and guide thimble or newly designing them. During its falling down into the core, the RCCA is retarded by various forces acting on it such as flow resistance and friction caused by the RCCA movement, buoyancy mechanical friction caused by contacting inner surface of the guide thimble, etc. However, complicated coupling of the various forces makes it difficult to derive an analytical dynamic equation for the drop time and impact velocity. This paper deals with the development of a computer program containing an analytical dynamic equation applicable to the Korean Fuel Assembly (KOFA) loaded in the Korean nuclear power plants. The computer program is benchmarked with an available single control rod drop tests. Since the predicted values are in good agreements with the test results, the computer program developed in this paper can be employed to modify the existing design features of the RCCA and guide thimble and to develop their new design features for advanced nuclear reactors. (author)

Stochastic Computational Approach for Complex Nonlinear Ordinary Differential Equations

International Nuclear Information System (INIS)

Khan, Junaid Ali; Raja, Muhammad Asif Zahoor; Qureshi, Ijaz Mansoor

2011-01-01

We present an evolutionary computational approach for the solution of nonlinear ordinary differential equations (NLODEs). The mathematical modeling is performed by a feed-forward artificial neural network that defines an unsupervised error. The training of these networks is achieved by a hybrid intelligent algorithm, a combination of global search with genetic algorithm and local search by pattern search technique. The applicability of this approach ranges from single order NLODEs, to systems of coupled differential equations. We illustrate the method by solving a variety of model problems and present comparisons with solutions obtained by exact methods and classical numerical methods. The solution is provided on a continuous finite time interval unlike the other numerical techniques with comparable accuracy. With the advent of neuroprocessors and digital signal processors the method becomes particularly interesting due to the expected essential gains in the execution speed. (general)

Efficient Numerical Methods for Stochastic Differential Equations in Computational Finance

KAUST Repository

Happola, Juho

2017-09-19

Stochastic Differential Equations (SDE) offer a rich framework to model the probabilistic evolution of the state of a system. Numerical approximation methods are typically needed in evaluating relevant Quantities of Interest arising from such models. In this dissertation, we present novel effective methods for evaluating Quantities of Interest relevant to computational finance when the state of the system is described by an SDE.

Efficient Numerical Methods for Stochastic Differential Equations in Computational Finance

KAUST Repository

Happola, Juho

2017-01-01

Stochastic Differential Equations (SDE) offer a rich framework to model the probabilistic evolution of the state of a system. Numerical approximation methods are typically needed in evaluating relevant Quantities of Interest arising from such models. In this dissertation, we present novel effective methods for evaluating Quantities of Interest relevant to computational finance when the state of the system is described by an SDE.

Analysis of computer programming languages

International Nuclear Information System (INIS)

Risset, Claude Alain

1967-01-01

This research thesis aims at trying to identify some methods of syntax analysis which can be used for computer programming languages while putting aside computer devices which influence the choice of the programming language and methods of analysis and compilation. In a first part, the author proposes attempts of formalization of Chomsky grammar languages. In a second part, he studies analytical grammars, and then studies a compiler or analytic grammar for the Fortran language

Programming Unconventional Computers: Dynamics, Development, Self-Reference

Directory of Open Access Journals (Sweden)

Susan Stepney

2012-10-01

Full Text Available Classical computing has well-established formalisms for specifying, refining, composing, proving, and otherwise reasoning about computations. These formalisms have matured over the past 70 years or so. Unconventional Computing includes the use of novel kinds of substrates–from black holes and quantum effects, through to chemicals, biomolecules, even slime moulds–to perform computations that do not conform to the classical model. Although many of these unconventional substrates can be coerced into performing classical computation, this is not how they “naturally” compute. Our ability to exploit unconventional computing is partly hampered by a lack of corresponding programming formalisms: we need models for building, composing, and reasoning about programs that execute in these substrates. What might, say, a slime mould programming language look like? Here I outline some of the issues and properties of these unconventional substrates that need to be addressed to find “natural” approaches to programming them. Important concepts include embodied real values, processes and dynamical systems, generative systems and their meta-dynamics, and embodied self-reference.

SIMULTANEOUS DIFFERENTIAL EQUATION COMPUTER

Science.gov (United States)

Collier, D.M.; Meeks, L.A.; Palmer, J.P.

1960-05-10

A description is given for an electronic simulator for a system of simultaneous differential equations, including nonlinear equations. As a specific example, a homogeneous nuclear reactor system including a reactor fluid, heat exchanger, and a steam boiler may be simulated, with the nonlinearity resulting from a consideration of temperature effects taken into account. The simulator includes three operational amplifiers, a multiplier, appropriate potential sources, and interconnecting R-C networks.

Computational study on full-wave inversion based on the elastic wave-equation; Dansei hado hoteishiki full wave inversion no model keisan ni yoru kento

Energy Technology Data Exchange (ETDEWEB)

Uesaka, S [Kyoto University, Kyoto (Japan). Faculty of Engineering; Watanabe, T; Sassa, K [Kyoto University, Kyoto (Japan)

1997-05-27

Algorithm is constructed and a program developed for a full-wave inversion (FWI) method utilizing the elastic wave equation in seismic exploration. The FWI method is a method for obtaining a physical property distribution using the whole observed waveforms as the data. It is capable of high resolution which is several times smaller than the wavelength since it can handle such phenomena as wave reflection and dispersion. The method for determining the P-wave velocity structure by use of the acoustic wave equation does not provide information about the S-wave velocity since it does not consider S-waves or converted waves. In an analysis using the elastic wave equation, on the other hand, not only P-wave data but also S-wave data can be utilized. In this report, under such circumstances, an inverse analysis algorithm is constructed on the basis of the elastic wave equation, and a basic program is developed. On the basis of the methods of Mora and of Luo and Schuster, the correction factors for P-wave and S-wave velocities are formulated directly from the elastic wave equation. Computations are performed and the effects of the hypocenter frequency and vibration transmission direction are examined. 6 refs., 8 figs.

Preschool Cookbook of Computer Programming Topics

Science.gov (United States)

Morgado, Leonel; Cruz, Maria; Kahn, Ken

2010-01-01

A common problem in computer programming use for education in general, not simply as a technical skill, is that children and teachers find themselves constrained by what is possible through limited expertise in computer programming techniques. This is particularly noticeable at the preliterate level, where constructs tend to be limited to…

Simple, accurate equations for human blood O2 dissociation computations.

Science.gov (United States)

Severinghaus, J W

1979-03-01

Hill's equation can be slightly modified to fit the standard human blood O2 dissociation curve to within plus or minus 0.0055 fractional saturation (S) from O less than S less than 1. Other modifications of Hill's equation may be used to compute Po2 (Torr) from S (Eq. 2), and the temperature coefficient of Po2 (Eq. 3). Variations of the Bohr coefficient with Po2 are given by Eq. 4. S = (((Po2(3) + 150 Po2)(-1) x 23,400) + 1)(-1) (1) In Po2 = 0.385 In (S-1 - 1)(-1) + 3.32 - (72 S)(-1) - 0.17(S6) (2) DELTA In Po2/delta T = 0.058 ((0.243 X Po2/100)(3.88) + 1)(-1) + 0.013 (3) delta In Po2/delta pH = (Po2/26.6)(0.184) - 2.2 (4) Procedures are described to determine Po2 and S of blood iteratively after extraction or addition of a defined amount of O2 and to compute P50 of blood from a single sample after measuring Po2, pH, and S.

FastChem: A computer program for efficient complex chemical equilibrium calculations in the neutral/ionized gas phase with applications to stellar and planetary atmospheres

Science.gov (United States)

Stock, Joachim W.; Kitzmann, Daniel; Patzer, A. Beate C.; Sedlmayr, Erwin

2018-06-01

For the calculation of complex neutral/ionized gas phase chemical equilibria, we present a semi-analytical versatile and efficient computer program, called FastChem. The applied method is based on the solution of a system of coupled nonlinear (and linear) algebraic equations, namely the law of mass action and the element conservation equations including charge balance, in many variables. Specifically, the system of equations is decomposed into a set of coupled nonlinear equations in one variable each, which are solved analytically whenever feasible to reduce computation time. Notably, the electron density is determined by using the method of Nelder and Mead at low temperatures. The program is written in object-oriented C++ which makes it easy to couple the code with other programs, although a stand-alone version is provided. FastChem can be used in parallel or sequentially and is available under the GNU General Public License version 3 at https://github.com/exoclime/FastChem together with several sample applications. The code has been successfully validated against previous studies and its convergence behavior has been tested even for extreme physical parameter ranges down to 100 K and up to 1000 bar. FastChem converges stable and robust in even most demanding chemical situations, which posed sometimes extreme challenges for previous algorithms.

On the kinetics of the aluminum-water reaction during exposure in high-heat flux test loops: 1, A computer program for oxidation calculations

International Nuclear Information System (INIS)

Pawel, R.E.

1988-01-01

The ''Griess Correlation,'' in which the thickness of the corrosion product on aluminum alloy surfaces is expressed as a function of time and temperature for high-flux-reactor conditions, was rewritten in the form of a simple, general rate equation. Based on this equation, a computer program that calculates oxide-layer thickness for any given time-temperature transient was written. 4 refs

Computation of high Reynolds number internal/external flows

Science.gov (United States)

Cline, M. C.; Wilmoth, R. G.

1981-01-01

A general, user oriented computer program, called VNAP2, was developed to calculate high Reynolds number, internal/ external flows. The VNAP2 program solves the two dimensional, time dependent Navier-Stokes equations. The turbulence is modeled with either a mixing-length, a one transport equation, or a two transport equation model. Interior grid points are computed using the explicit MacCormack Scheme with special procedures to speed up the calculation in the fine grid. All boundary conditions are calculated using a reference plane characteristic scheme with the viscous terms treated as source terms. Several internal, external, and internal/external flow calculations are presented.

A Fortran program for the numerical integration of the one-dimensional Schroedinger equation using exponential and Bessel fitting methods

International Nuclear Information System (INIS)

Cash, J.R.; Raptis, A.D.; Simos, T.E.

1990-01-01

An efficient algorithm is described for the accurate numerical integration of the one-dimensional Schroedinger equation. This algorithm uses a high-order, variable step Runge-Kutta like method in the region where the potential term dominates, and an exponential or Bessel fitted method in the asymptotic region. This approach can be used to compute scattering phase shifts in an efficient and reliable manner. A Fortran program which implements this algorithm is provided and some test results are given. (orig.)

Cable Connected Spinning Spacecraft, 1. the Canonical Equations, 2. Urban Mass Transportation, 3

Science.gov (United States)

Sitchin, A.

1972-01-01

Work on the dynamics of cable-connected spinning spacecraft was completed by formulating the equations of motion by both the canonical equations and Lagrange's equations and programming them for numerical solution on a digital computer. These energy-based formulations will permit future addition of the effect of cable mass. Comparative runs indicate that the canonical formulation requires less computer time. Available literature on urban mass transportation was surveyed. Areas of the private rapid transit concept of urban transportation are also studied.

Direct Computation of Sound Radiation by Jet Flow Using Large-scale Equations

Science.gov (United States)

Mankbadi, R. R.; Shih, S. H.; Hixon, D. R.; Povinelli, L. A.

1995-01-01

Jet noise is directly predicted using large-scale equations. The computational domain is extended in order to directly capture the radiated field. As in conventional large-eddy-simulations, the effect of the unresolved scales on the resolved ones is accounted for. Special attention is given to boundary treatment to avoid spurious modes that can render the computed fluctuations totally unacceptable. Results are presented for a supersonic jet at Mach number 2.1.

Computer Tutorial Programs in Physics.

Science.gov (United States)

Faughn, Jerry; Kuhn, Karl

1979-01-01

Describes a series of computer tutorial programs which are intended to help college students in introductory physics courses. Information about these programs, which are either calculus or algebra-trig based, is presented. (HM)

A computer program for activation analysis

International Nuclear Information System (INIS)

Rantanen, J.; Rosenberg, R.J.

1983-01-01

A computer program for calculating the results of activation analysis is described. The program comprises two gamma spectrum analysis programs, STOAV and SAMPO and one program for calculating elemental concentrations, KVANT. STOAV is based on a simple summation of channels and SAMPO is based on fitting of mathematical functions. The programs are tested by analyzing the IAEA G-1 test spectra. In the determination of peak location SAMPO is somewhat better than STOAV and in the determination of peak area SAMPO is more than twice as accurate as STOAV. On the other hand, SAMPO is three times as expensive as STOAV with the use of a Cyber 170 computer. (author)

32 CFR 701.125 - Computer matching program.

Science.gov (United States)

2010-07-01

... counterintelligence matches done in the course of performing a background check for security clearances of Federal... 32 National Defense 5 2010-07-01 2010-07-01 false Computer matching program. 701.125 Section 701... OF THE NAVY DOCUMENTS AFFECTING THE PUBLIC DON Privacy Program § 701.125 Computer matching program...

Macroscopic balance equations for two-phase flow models

International Nuclear Information System (INIS)

Hughes, E.D.

1979-01-01

The macroscopic, or overall, balance equations of mass, momentum, and energy are derived for a two-fluid model of two-phase flows in complex geometries. These equations provide a base for investigating methods of incorporating improved analysis methods into computer programs, such as RETRAN, which are used for transient and steady-state thermal-hydraulic analyses of nuclear steam supply systems. The equations are derived in a very general manner so that three-dimensional, compressible flows can be analysed. The equations obtained supplement the various partial differential equation two-fluid models of two-phase flow which have recently appeared in the literature. The primary objective of the investigation is the macroscopic balance equations. (Auth.)

Computer program to solve two-dimensional shock-wave interference problems with an equilibrium chemically reacting air model

Science.gov (United States)

Glass, Christopher E.

1990-08-01

The computer program EASI, an acronym for Equilibrium Air Shock Interference, was developed to calculate the inviscid flowfield, the maximum surface pressure, and the maximum heat flux produced by six shock wave interference patterns on a 2-D, cylindrical configuration. Thermodynamic properties of the inviscid flowfield are determined using either an 11-specie, 7-reaction equilibrium chemically reacting air model or a calorically perfect air model. The inviscid flowfield is solved using the integral form of the conservation equations. Surface heating calculations at the impingement point for the equilibrium chemically reacting air model use variable transport properties and specific heat. However, for the calorically perfect air model, heating rate calculations use a constant Prandtl number. Sample calculations of the six shock wave interference patterns, a listing of the computer program, and flowcharts of the programming logic are included.

Computational chemistry

Science.gov (United States)

Arnold, J. O.

1987-01-01

With the advent of supercomputers, modern computational chemistry algorithms and codes, a powerful tool was created to help fill NASA's continuing need for information on the properties of matter in hostile or unusual environments. Computational resources provided under the National Aerodynamics Simulator (NAS) program were a cornerstone for recent advancements in this field. Properties of gases, materials, and their interactions can be determined from solutions of the governing equations. In the case of gases, for example, radiative transition probabilites per particle, bond-dissociation energies, and rates of simple chemical reactions can be determined computationally as reliably as from experiment. The data are proving to be quite valuable in providing inputs to real-gas flow simulation codes used to compute aerothermodynamic loads on NASA's aeroassist orbital transfer vehicles and a host of problems related to the National Aerospace Plane Program. Although more approximate, similar solutions can be obtained for ensembles of atoms simulating small particles of materials with and without the presence of gases. Computational chemistry has application in studying catalysis, properties of polymers, all of interest to various NASA missions, including those previously mentioned. In addition to discussing these applications of computational chemistry within NASA, the governing equations and the need for supercomputers for their solution is outlined.

Program packages for dynamics systems analysis and design

International Nuclear Information System (INIS)

Athani, V.V.

1976-01-01

The development of computer program packages for dynamic system analysis and design are reported. The purpose of developing these program packages is to take the burden of writing computer programs off the mind of the system engineer and to enable him to concentrate on his main system analysis and design work. Towards this end, four standard computer program packages have been prepared : (1) TFANA - starting from system transfer function this program computes transient response, frequency response, root locus and stability by Routh Hurwitz criterion, (2) TFSYN - classical synthesis using algebraic method of Shipley, (3) MODANA - starting from state equations of the system this program computes solution of state equations, controllability, observability and stability, (4) OPTCON - This program obtains solutions of (i) linear regulator problem, (ii) servomechanism problems and (iii) problem of pole placement. The paper describes these program packages with the help of flowcharts and illustrates their use with the help of examples. (author)

SIMMER-II: A computer program for LMFBR disrupted core analysis

Energy Technology Data Exchange (ETDEWEB)

Bohl, W.R.; Luck, L.B.

1990-06-01

SIMMER-2 (Version 12) is a computer program to predict the coupled neutronic and fluid-dynamics behavior of liquid-metal fast reactors during core-disruptive accident transients. The modeling philosophy is based on the use of general, but approximate, physics to represent interactions of accident phenomena and regimes rather than a detailed representation of specialized situations. Reactor neutronic behavior is predicted by solving space (r,z), energy, and time-dependent neutron conservation equations (discrete ordinates transport or diffusion). The neutronics and the fluid dynamics are coupled via temperature- and background-dependent cross sections and the reactor power distribution. The fluid-dynamics calculation solves multicomponent, multiphase, multifield equations for mass, momentum, and energy conservation in (r,z) or (x,y) geometry. A structure field with nine density and five energy components; a liquid field with eight density and six energy components; and a vapor field with six density and on energy component are coupled by exchange functions representing a modified-dispersed flow regime with a zero-dimensional intra-cell structure model.

SIMMER-II: A computer program for LMFBR disrupted core analysis

International Nuclear Information System (INIS)

Bohl, W.R.; Luck, L.B.

1990-06-01

SIMMER-2 (Version 12) is a computer program to predict the coupled neutronic and fluid-dynamics behavior of liquid-metal fast reactors during core-disruptive accident transients. The modeling philosophy is based on the use of general, but approximate, physics to represent interactions of accident phenomena and regimes rather than a detailed representation of specialized situations. Reactor neutronic behavior is predicted by solving space (r,z), energy, and time-dependent neutron conservation equations (discrete ordinates transport or diffusion). The neutronics and the fluid dynamics are coupled via temperature- and background-dependent cross sections and the reactor power distribution. The fluid-dynamics calculation solves multicomponent, multiphase, multifield equations for mass, momentum, and energy conservation in (r,z) or (x,y) geometry. A structure field with nine density and five energy components; a liquid field with eight density and six energy components; and a vapor field with six density and on energy component are coupled by exchange functions representing a modified-dispersed flow regime with a zero-dimensional intra-cell structure model

Thermal transport properties of helium, helium--air mixtures, water, and tubing steel used in the CACHE program to compute HTGR auxiliary heat exchanger performance

International Nuclear Information System (INIS)

Tallackson, J.R.

1976-02-01

A description is presented of the thermal transport properties of the materials involved in digital computer calculations of heat transfer rates by the core auxiliary heat exchangers in large HTGR nuclear steam supply systems. These materials are pure helium, mixtures of helium with common gases having molecular weights in the range of 28 to 32, alloy steel tubing, and water. For use in programmed computations the viscosity, thermal conductivity, and specific heat are represented primarily by equations augmented by curves and tabulations. Materials supporting the development and selection of the property equations are included

Computer local construction of a general solution for the Chew-Low equations

International Nuclear Information System (INIS)

Gerdt, V.P.

1980-01-01

General solution of the dynamic form of the Chew-Low equations in the vicinity of the restpoint is considered. A method for calculating coefficients of series being members of such solution is suggested. The results of calculations, coefficients of power series and expansions carried out by means of the SCHOONSCHIP and SYMBAL systems are given. It is noted that the suggested procedure of the Chew-Low equation solutions basing on using an electronic computer as an instrument for analytical calculations permits to obtain detail information on the local structure of general solution

The Multigroup Neutron Diffusion Equations/1 Space Dimension

Energy Technology Data Exchange (ETDEWEB)

Linde, Sven

1960-06-15

A description is given of a program for the Ferranti Mercury computer which solves the one-dimensional multigroup diffusion equations in plane, cylindrical or spherical geometry, and also approximates automatically a two-dimensional solution by separating the space variables. In section A the method of calculation is outlined and the preparation of data for two group problems is described. The spatial separation of two-dimensional equations is considered in section B. Section C covers the multigroup equations. These parts are self contained and include all information required for the use of the program. Details of the numerical methods are given in section D. Three sample problems are solved in section E. Punching and operating instructions are given in an appendix.

The Multigroup Neutron Diffusion Equations/1 Space Dimension

International Nuclear Information System (INIS)

Linde, Sven

1960-06-01

A description is given of a program for the Ferranti Mercury computer which solves the one-dimensional multigroup diffusion equations in plane, cylindrical or spherical geometry, and also approximates automatically a two-dimensional solution by separating the space variables. In section A the method of calculation is outlined and the preparation of data for two group problems is described. The spatial separation of two-dimensional equations is considered in section B. Section C covers the multigroup equations. These parts are self contained and include all information required for the use of the program. Details of the numerical methods are given in section D. Three sample problems are solved in section E. Punching and operating instructions are given in an appendix

Characteristics and possibilities of computer program for fast assessment of primary frequency control of electric power interconnections

Directory of Open Access Journals (Sweden)

Ivanović Milan

2011-01-01

Full Text Available This paper presents the possibilities and practical features of a computer program for fast assessment of the effects of primary frequency regulation of electric power interconnections. It is based on two methods. The first one is the analytical method, which applies analytical expressions for the non-zero initial conditions, with a range of benefits provided by the analytical form, allowing consideration of possible structural changes in the power system during the analysis process. The second is a simulation method, with recurrent application of suitable drafted, fully decoupled difference equations. Capabilities and features of this computer program have been identified in case of isolated power system of Serbia, and then for the case of a widespread appreciation of its surrounding.

Simulation of the Schroedinger equation on SHAC

International Nuclear Information System (INIS)

Stewart, A.

1976-01-01

A simulation of the Schroedinger wave equation for the hydrogen atom, on SHAC, a simple homogeneous analogue computer primarily intended for use in schools, is described. Due to the incorporation of FET switches very high speed switching from initial conditions to compute modes is possible. The techniques employed in the multiplier and divider are discussed and the flow diagram for the Schroedinger program shown. Results and photographs are discussed. (U.K.)

A first course in structural equation modeling

CERN Document Server

Raykov, Tenko

2012-01-01

In this book, authors Tenko Raykov and George A. Marcoulides introduce students to the basics of structural equation modeling (SEM) through a conceptual, nonmathematical approach. For ease of understanding, the few mathematical formulas presented are used in a conceptual or illustrative nature, rather than a computational one.Featuring examples from EQS, LISREL, and Mplus, A First Course in Structural Equation Modeling is an excellent beginner's guide to learning how to set up input files to fit the most commonly used types of structural equation models with these programs. The basic ideas and methods for conducting SEM are independent of any particular software.Highlights of the Second Edition include: Review of latent change (growth) analysis models at an introductory level Coverage of the popular Mplus program Updated examples of LISREL and EQS A CD that contains all of the text's LISREL, EQS, and Mplus examples.A First Course in Structural Equation Modeling is intended as an introductory book for students...

Functional Programming in Computer Science

Energy Technology Data Exchange (ETDEWEB)

Anderson, Loren James [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Davis, Marion Kei [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2016-01-19

We explore functional programming through a 16-week internship at Los Alamos National Laboratory. Functional programming is a branch of computer science that has exploded in popularity over the past decade due to its high-level syntax, ease of parallelization, and abundant applications. First, we summarize functional programming by listing the advantages of functional programming languages over the usual imperative languages, and we introduce the concept of parsing. Second, we discuss the importance of lambda calculus in the theory of functional programming. Lambda calculus was invented by Alonzo Church in the 1930s to formalize the concept of effective computability, and every functional language is essentially some implementation of lambda calculus. Finally, we display the lasting products of the internship: additions to a compiler and runtime system for the pure functional language STG, including both a set of tests that indicate the validity of updates to the compiler and a compiler pass that checks for illegal instances of duplicate names.

An iterative algorithm for solving the multidimensional neutron diffusion nodal method equations on parallel computers

International Nuclear Information System (INIS)

Kirk, B.L.; Azmy, Y.Y.

1992-01-01

In this paper the one-group, steady-state neutron diffusion equation in two-dimensional Cartesian geometry is solved using the nodal integral method. The discrete variable equations comprise loosely coupled sets of equations representing the nodal balance of neutrons, as well as neutron current continuity along rows or columns of computational cells. An iterative algorithm that is more suitable for solving large problems concurrently is derived based on the decomposition of the spatial domain and is accelerated using successive overrelaxation. This algorithm is very well suited for parallel computers, especially since the spatial domain decomposition occurs naturally, so that the number of iterations required for convergence does not depend on the number of processors participating in the calculation. Implementation of the authors' algorithm on the Intel iPSC/2 hypercube and Sequent Balance 8000 parallel computer is presented, and measured speedup and efficiency for test problems are reported. The results suggest that the efficiency of the hypercube quickly deteriorates when many processors are used, while the Sequent Balance retains very high efficiency for a comparable number of participating processors. This leads to the conjecture that message-passing parallel computers are not as well suited for this algorithm as shared-memory machines

Computer Series, 99: Bits and Pieces, 39.

Science.gov (United States)

Moore, John W., Ed.

1989-01-01

Presents five computer programs: (1) Accurate Numerical Solutions of the One-Dimensional Schrodinger Equation; (2) NMR Simulation and Interactive Drill/Interpretation; (3) A Simple Computer Program for the Calculation of 13C-NMR Chemical Shifts; (4) Constants of 1:1 Complexes from NMR or Spectrophotometric Measurements; and (5) Saturation…

Computer programs simplify optical system analysis

Science.gov (United States)

1965-01-01

The optical ray-trace computer program performs geometrical ray tracing. The energy-trace program calculates the relative monochromatic flux density on a specific target area. This program uses the ray-trace program as a subroutine to generate a representation of the optical system.

Computer Programming Languages for Health Care

Science.gov (United States)

O'Neill, Joseph T.

1979-01-01

This paper advocates the use of standard high level programming languages for medical computing. It recommends that U.S. Government agencies having health care missions implement coordinated policies that encourage the use of existing standard languages and the development of new ones, thereby enabling them and the medical computing community at large to share state-of-the-art application programs. Examples are based on a model that characterizes language and language translator influence upon the specification, development, test, evaluation, and transfer of application programs.

HTA educational outreach program and change the equation participation

Science.gov (United States)

Gordon, Robert

2013-05-01

In this presentation, Hitachi High Technologies America (HTA) introduces its Educational Outreach Program and explains it's involvement with Change The Equation (CTEq), a nonprofit, nonpartisan, CEO-led initiative that is mobilizing the business community to improve the quality of science, technology, engineering and mathematics (STEM) learning in the United States.

MAKSIMA-CHEMIST: a program for Mass Action Kinetics Simulation by Automatic Chemical Equation Manipulation and Integration using Stiff Techniques

International Nuclear Information System (INIS)

Carver, M.B.; Hanley, D.V.; Chaplin, K.R.

1979-02-01

MAKSIMA-CHEMIST was written to compute the kinetics of simultaneous chemical reactions. The ordinary differential equations, which are automatically derived from the stated chemical equations, are difficult to integrate, as they are coupled in a highly nonlinear manner and frequently involve a large range in the magnitude of the reaction rates. They form a classic 'stiff' differential equaton set which can be integrated efficiently only by recently developed advanced techniques. The new program also contains provision for higher order chemical reactions, and has a dynamic storage and decision feature. This permits it to accept any number of chemical reactions and species, and choose an integraton scheme which will perform most efficiently within the available memory. Sparse matrix techniques are used when the size and structure of the equation set is suitable. Finally, a number of post-analysis options are available, including printer and Calcomp plots of transient response of selected species, and graphical representation of the reaction matrix. (auth)

Development and validation of the computer program TNHXY

International Nuclear Information System (INIS)

Xolocostli M, V.; Valle G, E. del; Alonso V, G.

2003-01-01

This work describes the development and validation of the computer program TNHXY (Neutron Transport with Nodal Hybrid schemes in X Y geometry), which solves the discrete-ordinates neutron transport equations using a discontinuous Bi-Linear (DBiL) nodal hybrid method. One of the immediate applications of TNHXY is in the analysis of nuclear fuel assemblies, in particular those of BWRs. Its validation was carried out by reproducing some results for test or benchmark problems that some authors have solved using other numerical techniques. This allows to ensure that the program will provide results with similar accuracy for other problems of the same type. To accomplish this two benchmark problems have been solved. The first problem consists in a BWR fuel assembly in a 7x7 array without and with control rod. The results obtained with TNHXY are consistent with those reported for the TWOTRAN code. The second benchmark problem is a Mixed Oxide (MOX) fuel assembly in a 10x10 array. This last problem is known as the WPPR benchmark problem of the NEA Data Bank and the results are compared with those obtained with commercial codes like HELIOS, MCNP-4B and CPM-3. (Author)

Development and adjustment of programs for solving systems of linear equations

International Nuclear Information System (INIS)

Fujimura, Toichiro

1978-03-01

Programs for solving the systems of linear equations have been adjusted and developed in expanding the scientific subroutine library SSL. The principal programs adjusted are based on the congruent method, method of product form of the inverse, orthogonal method, Crout's method for sparse system, and acceleration of iterative methods. The programs developed are based on the escalator method, direct parallel residue method and block tridiagonal method for band system. Described are usage of the programs developed and their future improvement. FORTRAN lists with simple examples in tests of the programs are also given. (auth.)

The RANDOM computer program: A linear congruential random number generator

Science.gov (United States)

Miles, R. F., Jr.

1986-01-01

The RANDOM Computer Program is a FORTRAN program for generating random number sequences and testing linear congruential random number generators (LCGs). The linear congruential form of random number generator is discussed, and the selection of parameters of an LCG for a microcomputer described. This document describes the following: (1) The RANDOM Computer Program; (2) RANDOM.MOD, the computer code needed to implement an LCG in a FORTRAN program; and (3) The RANCYCLE and the ARITH Computer Programs that provide computational assistance in the selection of parameters for an LCG. The RANDOM, RANCYCLE, and ARITH Computer Programs are written in Microsoft FORTRAN for the IBM PC microcomputer and its compatibles. With only minor modifications, the RANDOM Computer Program and its LCG can be run on most micromputers or mainframe computers.

Gender Differences in the Use of Computers, Programming, and Peer Interactions in Computer Science Classrooms

Science.gov (United States)

Stoilescu, Dorian; Egodawatte, Gunawardena

2010-01-01

Research shows that female and male students in undergraduate computer science programs view computer culture differently. Female students are interested more in the use of computers than in doing programming, whereas male students see computer science mainly as a programming activity. The overall purpose of our research was not to find new…

Gender differences in the use of computers, programming, and peer interactions in computer science classrooms

Science.gov (United States)

Stoilescu, Dorian; Egodawatte, Gunawardena

2010-12-01

Research shows that female and male students in undergraduate computer science programs view computer culture differently. Female students are interested more in the use of computers than in doing programming, whereas male students see computer science mainly as a programming activity. The overall purpose of our research was not to find new definitions for computer science culture but to see how male and female students see themselves involved in computer science practices, how they see computer science as a successful career, and what they like and dislike about current computer science practices. The study took place in a mid-sized university in Ontario. Sixteen students and two instructors were interviewed to get their views. We found that male and female views are different on computer use, programming, and the pattern of student interactions. Female and male students did not have any major issues in using computers. In computing programming, female students were not so involved in computing activities whereas male students were heavily involved. As for the opinions about successful computer science professionals, both female and male students emphasized hard working, detailed oriented approaches, and enjoying playing with computers. The myth of the geek as a typical profile of successful computer science students was not found to be true.

Converting differential-equation models of biological systems to membrane computing.

Science.gov (United States)

Muniyandi, Ravie Chandren; Zin, Abdullah Mohd; Sanders, J W

2013-12-01

This paper presents a method to convert the deterministic, continuous representation of a biological system by ordinary differential equations into a non-deterministic, discrete membrane computation. The dynamics of the membrane computation is governed by rewrite rules operating at certain rates. That has the advantage of applying accurately to small systems, and to expressing rates of change that are determined locally, by region, but not necessary globally. Such spatial information augments the standard differentiable approach to provide a more realistic model. A biological case study of the ligand-receptor network of protein TGF-β is used to validate the effectiveness of the conversion method. It demonstrates the sense in which the behaviours and properties of the system are better preserved in the membrane computing model, suggesting that the proposed conversion method may prove useful for biological systems in particular. Copyright © 2013 Elsevier Ireland Ltd. All rights reserved.

Computer technique for evaluating collimator performance

International Nuclear Information System (INIS)

Rollo, F.D.

1975-01-01

A computer program has been developed to theoretically evaluate the overall performance of collimators used with radioisotope scanners and γ cameras. The first step of the program involves the determination of the line spread function (LSF) and geometrical efficiency from the fundamental parameters of the collimator being evaluated. The working equations can be applied to any plane of interest. The resulting LSF is applied to subroutine computer programs which compute corresponding modulation transfer function and contrast efficiency functions. The latter function is then combined with appropriate geometrical efficiency data to determine the performance index function. The overall computer program allows one to predict from the physical parameters of the collimator alone how well the collimator will reproduce various sized spherical voids of activity in the image plane. The collimator performance program can be used to compare the performance of various collimator types, to study the effects of source depth on collimator performance, and to assist in the design of collimators. The theory of the collimator performance equation is discussed, a comparison between the experimental and theoretical LSF values is made, and examples of the application of the technique are presented

A CAD (Classroom Assessment Design) of a Computer Programming Course

Science.gov (United States)

Hawi, Nazir S.

2012-01-01

This paper presents a CAD (classroom assessment design) of an entry-level undergraduate computer programming course "Computer Programming I". CAD has been the product of a long experience in teaching computer programming courses including teaching "Computer Programming I" 22 times. Each semester, CAD is evaluated and modified…

Recovery Act: Finite Volume Based Computer Program for Ground Source Heat Pump Systems

Energy Technology Data Exchange (ETDEWEB)

James A Menart, Professor

2013-02-22

This report is a compilation of the work that has been done on the grant DE-EE0002805 entitled Finite Volume Based Computer Program for Ground Source Heat Pump Systems. The goal of this project was to develop a detailed computer simulation tool for GSHP (ground source heat pump) heating and cooling systems. Two such tools were developed as part of this DOE (Department of Energy) grant; the first is a two-dimensional computer program called GEO2D and the second is a three-dimensional computer program called GEO3D. Both of these simulation tools provide an extensive array of results to the user. A unique aspect of both these simulation tools is the complete temperature profile information calculated and presented. Complete temperature profiles throughout the ground, casing, tube wall, and fluid are provided as a function of time. The fluid temperatures from and to the heat pump, as a function of time, are also provided. In addition to temperature information, detailed heat rate information at several locations as a function of time is determined. Heat rates between the heat pump and the building indoor environment, between the working fluid and the heat pump, and between the working fluid and the ground are computed. The heat rates between the ground and the working fluid are calculated as a function time and position along the ground loop. The heating and cooling loads of the building being fitted with a GSHP are determined with the computer program developed by DOE called ENERGYPLUS. Lastly COP (coefficient of performance) results as a function of time are provided. Both the two-dimensional and three-dimensional computer programs developed as part of this work are based upon a detailed finite volume solution of the energy equation for the ground and ground loop. Real heat pump characteristics are entered into the program and used to model the heat pump performance. Thus these computer tools simulate the coupled performance of the ground loop and the heat pump. The

Solution of the transient Fourier heat conduction equation in r,phi geometry

International Nuclear Information System (INIS)

Kowa, E.; Ehnis, L.

1978-11-01

The two-dimensional transient Fourier heat conduction equation is solved in r,phi geometry for anisotropic materials with the computer program TERFI. The Alternating-Direction-Implicit method is used for the solution of this equation with specified start- and boundary conditions, temperature dependent material properties and space dependent heat sources. The solution area is devided in a mesh grid by the finite difference method. Slidely non-orthogonaly geometry (displacement of mesh grid) can be regarded. There were some difficulties in the treatment of the boundary conditions for the circularly-closed solution area because of the continuity of temperature and heat flux on the 0 0 /360 0 -line. This problem can be solved by an iterativ method with different starting points for the solution scheme. Emphasis was put on reaching reasonable computer time for the iteration. The computer code TERFI, programed in FORTRAN IV, is a modul of the program system RSYST. As an example the temperature distribution of a PWR fuel rod is calculated. (orig.) [de

Parallel computing: numerics, applications, and trends

National Research Council Canada - National Science Library

Trobec, Roman; Vajteršic, Marián; Zinterhof, Peter

2009-01-01

... and/or distributed systems. The contributions to this book are focused on topics most concerned in the trends of today's parallel computing. These range from parallel algorithmics, programming, tools, network computing to future parallel computing. Particular attention is paid to parallel numerics: linear algebra, differential equations, numerica...

Programming Paradigms in Computer Science Education

OpenAIRE

Bolshakova, Elena

2005-01-01

Main styles, or paradigms of programming – imperative, functional, logic, and object-oriented – are shortly described and compared, and corresponding programming techniques are outlined. Programming languages are classified in accordance with the main style and techniques supported. It is argued that profound education in computer science should include learning base programming techniques of all main programming paradigms.

MP Salsa: a finite element computer program for reacting flow problems. Part 1--theoretical development

Energy Technology Data Exchange (ETDEWEB)

Shadid, J.N.; Moffat, H.K.; Hutchinson, S.A.; Hennigan, G.L.; Devine, K.D.; Salinger, A.G.

1996-05-01

The theoretical background for the finite element computer program, MPSalsa, is presented in detail. MPSalsa is designed to solve laminar, low Mach number, two- or three-dimensional incompressible and variable density reacting fluid flows on massively parallel computers, using a Petrov-Galerkin finite element formulation. The code has the capability to solve coupled fluid flow, heat transport, multicomponent species transport, and finite-rate chemical reactions, and to solver coupled multiple Poisson or advection-diffusion- reaction equations. The program employs the CHEMKIN library to provide a rigorous treatment of multicomponent ideal gas kinetics and transport. Chemical reactions occurring in the gas phase and on surfaces are treated by calls to CHEMKIN and SURFACE CHEMKIN, respectively. The code employs unstructured meshes, using the EXODUS II finite element data base suite of programs for its input and output files. MPSalsa solves both transient and steady flows by using fully implicit time integration, an inexact Newton method and iterative solvers based on preconditioned Krylov methods as implemented in the Aztec solver library.

COMMIX-PPC: A three-dimensional transient multicomponent computer program for analyzing performance of power plant condensers

International Nuclear Information System (INIS)

Chien, T.H.; Domanus, H.M.; Sha, W.T.

1993-02-01

The COMMIX-PPC computer pregrain is an extended and improved version of earlier COMMIX codes and is specifically designed for evaluating the thermal performance of power plant condensers. The COMMIX codes are general-purpose computer programs for the analysis of fluid flow and heat transfer in complex Industrial systems. In COMMIX-PPC, two major features have been added to previously published COMMIX codes. One feature is the incorporation of one-dimensional equations of conservation of mass, momentum, and energy on the tube stile and the proper accounting for the thermal interaction between shell and tube side through the porous-medium approach. The other added feature is the extension of the three-dimensional conservation equations for shell-side flow to treat the flow of a multicomponent medium. COMMIX-PPC is designed to perform steady-state and transient. Three-dimensional analysis of fluid flow with heat transfer tn a power plant condenser. However, the code is designed in a generalized fashion so that, with some modification, it can be used to analyze processes in any heat exchanger or other single-phase engineering applications. Volume I (Equations and Numerics) of this report describes in detail the basic equations, formulation, solution procedures, and models for a phenomena. Volume II (User's Guide and Manual) contains the input instruction, flow charts, sample problems, and descriptions of available options and boundary conditions

Application of computers in a radiological survey program

International Nuclear Information System (INIS)

Berven, B.A.; Blair, M.S.; Doane, R.W.; Little, C.A.; Perdue, P.T.

1984-01-01

Computers have become increasingly important in data analysis and data management as well as assisting in report preparation in the Oak Ridge National Laboratory (ORNL) Radiological Survey Activities (RASA) Program. The primary function of the RASA program is to collect, analyze, report, and manage data collected to characterize the radiological condition of potentially contaminated sites identified in the Department of Energy's (DOE) remedial action programs. Three different computer systems are routinely utilized in ORNL/RASA operations. Two of these systems are employed in specific functions. A Nuclear Data (ND) 682 is used to perform isotopic analysis of gamma spectroscopic data generated by high-purity germanium detectors for air, water and soil samples. The ND682 employs a 16,000-channel analyzer that is routinely used with four germanium spectrometers. Word processing and data management are accomplished using the INtext system implemented on a DEC PDP-11 computer. A group of personal computers are used to perform a diverse number of functions. These computer systems are Commodore Business Machines (CBM) Model 8032 with a dual floppy disk storage medium and line printers (with optional X-Y plotters). The CBM's are utilized for: (1) data analysis -- raw data from radiation detection instrumentation are stored and manipulated with customized computer programs; (2) data reduction -- raw data are converted into report-ready tables using customized programs; (3) data management -- radionuclide data on each air, water and soil sample are stored on diskettes along with location of archived samples; and (4) program management -- site surveys and report status are tracked by computer files as well as program budget information to provide contemporary information of program status

Computation of reactor control rod drop time under accident conditions

International Nuclear Information System (INIS)

Dou Yikang; Yao Weida; Yang Renan; Jiang Nanyan

1998-01-01

The computational method of reactor control rod drop time under accident conditions lies mainly in establishing forced vibration equations for the components under action of outside forces on control rod driven line and motion equation for the control rod moving in vertical direction. The above two kinds of equations are connected by considering the impact effects between control rod and its outside components. Finite difference method is adopted to make discretization of the vibration equations and Wilson-θ method is applied to deal with the time history problem. The non-linearity caused by impact is iteratively treated with modified Newton method. Some experimental results are used to validate the validity and reliability of the computational method. Theoretical and experimental testing problems show that the computer program based on the computational method is applicable and reliable. The program can act as an effective tool of design by analysis and safety analysis for the relevant components

Computational fluid dynamics on a massively parallel computer

Science.gov (United States)

Jespersen, Dennis C.; Levit, Creon

1989-01-01

A finite difference code was implemented for the compressible Navier-Stokes equations on the Connection Machine, a massively parallel computer. The code is based on the ARC2D/ARC3D program and uses the implicit factored algorithm of Beam and Warming. The codes uses odd-even elimination to solve linear systems. Timings and computation rates are given for the code, and a comparison is made with a Cray XMP.

On the Dynamic Programming Approach for the 3D Navier-Stokes Equations

International Nuclear Information System (INIS)

Manca, Luigi

2008-01-01

The dynamic programming approach for the control of a 3D flow governed by the stochastic Navier-Stokes equations for incompressible fluid in a bounded domain is studied. By a compactness argument, existence of solutions for the associated Hamilton-Jacobi-Bellman equation is proved. Finally, existence of an optimal control through the feedback formula and of an optimal state is discussed

Exploring Students' Understanding of Ordinary Differential Equations Using Computer Algebraic System (CAS)

Science.gov (United States)

Maat, Siti Mistima; Zakaria, Effandi

2011-01-01

Ordinary differential equations (ODEs) are one of the important topics in engineering mathematics that lead to the understanding of technical concepts among students. This study was conducted to explore the students' understanding of ODEs when they solve ODE questions using a traditional method as well as a computer algebraic system, particularly…

Computing with linear equations and matrices

International Nuclear Information System (INIS)

Churchhouse, R.F.

1983-01-01

Systems of linear equations and matrices arise in many disciplines. The equations may accurately represent conditions satisfied by a system or, more likely, provide an approximation to a more complex system of non-linear or differential equations. The system may involve a few or many thousand unknowns and each individual equation may involve few or many of them. Over the past 50 years a vast literature on methods for solving systems of linear equations and the associated problems of finding the inverse or eigenvalues of a matrix has been produced. These lectures cover those methods which have been found to be most useful for dealing with such types of problem. References are given where appropriate and attention is drawn to the possibility of improved methods for use on vector and parallel processors. (orig.)

DCHAIN: A user-friendly computer program for radioactive decay and reaction chain calculations

International Nuclear Information System (INIS)

East, L.V.

1994-05-01

A computer program for calculating the time-dependent daughter populations in radioactive decay and nuclear reaction chains is described. Chain members can have non-zero initial populations and be produced from the preceding chain member as the result of radioactive decay, a nuclear reaction, or both. As presently implemented, chains can contain up to 15 members. Program input can be supplied interactively or read from ASCII data files. Time units for half-lives, etc. can be specified during data entry. Input values are verified and can be modified if necessary, before used in calculations. Output results can be saved in ASCII files in a format suitable for including in reports or other documents. The calculational method, described in some detail, utilizes a generalized form of the Bateman equations. The program is written in the C language in conformance with current ANSI standards and can be used on multiple hardware platforms

Introduction to the Atari Computer. A Program Written in the Pilot Programming Language.

Science.gov (United States)

Schlenker, Richard M.

Designed to be an introduction to the Atari microcomputers for beginners, the interactive computer program listed in this document is written in the Pilot programing language. Instructions are given for entering and storing the program in the computer memory for use by students. (MES)

Solution of the kinetic equation in the P3-approximation in a plane geometry

International Nuclear Information System (INIS)

Vlasov, Yu.A.

1975-01-01

A method and a program are described for solving single-velocity kinetic equations of neutron transfer for the plane geometry in the finite-difference approximation. A difference high-accuracy scheme and a matrix factorization method are used for the differential-difference equation systems. The program is written in the ALGOL-60 language and is adapted for M-20, M-220, M-222 and BESM-4 computers

Development of a computer program for the simulation of one-dimensional fixed- and moving-bed reactors

International Nuclear Information System (INIS)

Hartner, P.

1996-11-01

Chemical reactors with a flow through a bed of solid particles are of great importance in the processing industry. Modern computational tools allow for an improved characterization of the complex facts in such reactors leading to new opportunities of optimizing the reactor operation and environmental effects. This thesis is concerned with the development of the one-dimensional simulation software REASIM. The program covers the effects within a reacting bed and is designed for fixed and moving beds. To describe the reactor the balances for energy, momentum and mass are solved. The drying of the particles, pyrolysis and chemical gas-solid and gas-gas reactions are considered. For the description of the chemical gas-solid reactions a particle model for porous solids is developed. The calculation of mass transfer and of chemical reactions is strictly separated. All parameters necessary for the model can be measured in the laboratory. The model equations form a system of partial differential equations. This system is transformed to a set of ordinary differential equations. It is found that the best discretization method is the method of finite differences with the upwind-scheme for situations where convection is strong. The program has a modular structure making it is easy to replace parts of the program by new, improved modules if they become available. (author)

Solution of three-dimensional energy equation using finite element method

International Nuclear Information System (INIS)

Bhasin, V.; Singh, R.K.; Dutta, B.K.; Kushwaha, H.S.

1993-01-01

In the present work an attempt has been made to formulate an efficient 3-D finite element program for solving coupled momentum-energy equation with unsymmetric frontal solver and a suitable upwinding scheme. Based on the above solution technique of energy equation it can be concluded that upwinding scheme can lead to fairly accurate and smooth results even with coarse mesh. Otherwise the mesh size requirement will be extremely stringent for most of the practical problems. With upwinding the additional computer time required is marginally more. This effort has resulted in getting practical solution for large size real life problems in nuclear industry. The program was used for computation of temperature field in heavy water moderator of Madras Atomic Power Station (MAPS) reactor, in new mode of operation. (author). 9 refs., 7 figs

A Computer Program for Modeling the Conversion of Organic Waste to Energy

Directory of Open Access Journals (Sweden)

Pragasen Pillay

2011-11-01

Full Text Available This paper presents a tool for the analysis of conversion of organic waste into energy. The tool is a program that uses waste characterization parameters and mass flow rates at each stage of the waste treatment process to predict the given products. The specific waste treatment process analysed in this paper is anaerobic digestion. The different waste treatment stages of the anaerobic digestion process are: conditioning of input waste, secondary treatment, drying of sludge, conditioning of digestate, treatment of digestate, storage of liquid and solid effluent, disposal of liquid and solid effluents, purification, utilization and storage of combustible gas. The program uses mass balance equations to compute the amount of CH4, NH3, CO2 and H2S produced from anaerobic digestion of organic waste, and hence the energy available. Case studies are also presented.

Computer programming and architecture the VAX

CERN Document Server

Levy, Henry

2014-01-01

Takes a unique systems approach to programming and architecture of the VAXUsing the VAX as a detailed example, the first half of this book offers a complete course in assembly language programming. The second describes higher-level systems issues in computer architecture. Highlights include the VAX assembler and debugger, other modern architectures such as RISCs, multiprocessing and parallel computing, microprogramming, caches and translation buffers, and an appendix on the Berkeley UNIX assembler.

Human Memory Organization for Computer Programs.

Science.gov (United States)

Norcio, A. F.; Kerst, Stephen M.

1983-01-01

Results of study investigating human memory organization in processing of computer programming languages indicate that algorithmic logic segments form a cognitive organizational structure in memory for programs. Statement indentation and internal program documentation did not enhance organizational process of recall of statements in five Fortran…

Structured Design Language for Computer Programs

Science.gov (United States)

Pace, Walter H., Jr.

1986-01-01

Box language used at all stages of program development. Developed to provide improved productivity in designing, coding, and maintaining computer programs. BOX system written in FORTRAN 77 for batch execution.

Computational Materials Program for Alloy Design

Science.gov (United States)

Bozzolo, Guillermo

2005-01-01

The research program sponsored by this grant, "Computational Materials Program for Alloy Design", covers a period of time of enormous change in the emerging field of computational materials science. The computational materials program started with the development of the BFS method for alloys, a quantum approximate method for atomistic analysis of alloys specifically tailored to effectively deal with the current challenges in the area of atomistic modeling and to support modern experimental programs. During the grant period, the program benefited from steady growth which, as detailed below, far exceeds its original set of goals and objectives. Not surprisingly, by the end of this grant, the methodology and the computational materials program became an established force in the materials communitiy, with substantial impact in several areas. Major achievements during the duration of the grant include the completion of a Level 1 Milestone for the HITEMP program at NASA Glenn, consisting of the planning, development and organization of an international conference held at the Ohio Aerospace Institute in August of 2002, finalizing a period of rapid insertion of the methodology in the research community worlwide. The conference, attended by citizens of 17 countries representing various fields of the research community, resulted in a special issue of the leading journal in the area of applied surface science. Another element of the Level 1 Milestone was the presentation of the first version of the Alloy Design Workbench software package, currently known as "adwTools". This software package constitutes the first PC-based piece of software for atomistic simulations for both solid alloys and surfaces in the market.Dissemination of results and insertion in the materials community worldwide was a primary focus during this period. As a result, the P.I. was responsible for presenting 37 contributed talks, 19 invited talks, and publishing 71 articles in peer-reviewed journals, as

Computation of rational solutions for a first-order nonlinear differential equation

Directory of Open Access Journals (Sweden)

Djilali Behloul

2011-09-01

Full Text Available In this article, we study differential equations of the form $y'=sum A_i(xy^i/sum B_i(xy^i$ which can be elliptic, hyperbolic, parabolic, Riccati, or quasi-linear. We show how rational solutions can be computed in a systematic manner. Such results are most likely to find applications in the theory of limit cycles as indicated by Gine et al [4].

A highly efficient parallel algorithm for solving the neutron diffusion nodal equations on shared-memory computers

International Nuclear Information System (INIS)

Azmy, Y.Y.; Kirk, B.L.

1990-01-01

Modern parallel computer architectures offer an enormous potential for reducing CPU and wall-clock execution times of large-scale computations commonly performed in various applications in science and engineering. Recently, several authors have reported their efforts in developing and implementing parallel algorithms for solving the neutron diffusion equation on a variety of shared- and distributed-memory parallel computers. Testing of these algorithms for a variety of two- and three-dimensional meshes showed significant speedup of the computation. Even for very large problems (i.e., three-dimensional fine meshes) executed concurrently on a few nodes in serial (nonvector) mode, however, the measured computational efficiency is very low (40 to 86%). In this paper, the authors present a highly efficient (∼85 to 99.9%) algorithm for solving the two-dimensional nodal diffusion equations on the Sequent Balance 8000 parallel computer. Also presented is a model for the performance, represented by the efficiency, as a function of problem size and the number of participating processors. The model is validated through several tests and then extrapolated to larger problems and more processors to predict the performance of the algorithm in more computationally demanding situations

Computational Physics Program of the National MFE Computer Center

International Nuclear Information System (INIS)

Mirin, A.A.

1984-12-01

The principal objective of the computational physics group is to develop advanced numerical models for the investigation of plasma phenomena and the simulation of present and future magnetic confinement devices. A summary of the groups activities is presented, including computational studies in MHD equilibria and stability, plasma transport, Fokker-Planck, and efficient numerical and programming algorithms. References are included

ARDISC (Argonne Dispersion Code): computer programs to calculate the distribution of trace element migration in partially equilibrating media

International Nuclear Information System (INIS)

Strickert, R.; Friedman, A.M.; Fried, S.

1979-04-01

A computer program (ARDISC, the Argonne Dispersion Code) is described which simulates the migration of nuclides in porous media and includes first order kinetic effects on the retention constants. The code allows for different absorption and desorption rates and solves the coupled migration equations by arithmetic reiterations. Input data needed are the absorption and desorption rates, equilibrium surface absorption coefficients, flow rates and volumes, and media porosities

XMDS2: Fast, scalable simulation of coupled stochastic partial differential equations

Science.gov (United States)

Dennis, Graham R.; Hope, Joseph J.; Johnsson, Mattias T.

2013-01-01

XMDS2 is a cross-platform, GPL-licensed, open source package for numerically integrating initial value problems that range from a single ordinary differential equation up to systems of coupled stochastic partial differential equations. The equations are described in a high-level XML-based script, and the package generates low-level optionally parallelised C++ code for the efficient solution of those equations. It combines the advantages of high-level simulations, namely fast and low-error development, with the speed, portability and scalability of hand-written code. XMDS2 is a complete redesign of the XMDS package, and features support for a much wider problem space while also producing faster code. Program summaryProgram title: XMDS2 Catalogue identifier: AENK_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AENK_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 2 No. of lines in distributed program, including test data, etc.: 872490 No. of bytes in distributed program, including test data, etc.: 45522370 Distribution format: tar.gz Programming language: Python and C++. Computer: Any computer with a Unix-like system, a C++ compiler and Python. Operating system: Any Unix-like system; developed under Mac OS X and GNU/Linux. RAM: Problem dependent (roughly 50 bytes per grid point) Classification: 4.3, 6.5. External routines: The external libraries required are problem-dependent. Uses FFTW3 Fourier transforms (used only for FFT-based spectral methods), dSFMT random number generation (used only for stochastic problems), MPI message-passing interface (used only for distributed problems), HDF5, GNU Scientific Library (used only for Bessel-based spectral methods) and a BLAS implementation (used only for non-FFT-based spectral methods). Nature of problem: General coupled initial-value stochastic partial differential equations. Solution method: Spectral method

Application of computers in a Radiological Survey Program

International Nuclear Information System (INIS)

Berven, B.A.; Blair, M.S.; Doane, R.W.; Little, C.A.; Perdue, P.T.

1984-01-01

A brief description of some of the applications of computers in a radiological survey program is presented. It has been our experience that computers and computer software have allowed our staff personnel to more productively use their time by using computers to perform the mechanical acquisition, analyses, and storage of data. It is hoped that other organizations may similarly profit from this experience. This effort will ultimately minimize errors and reduce program costs

Translator program converts computer printout into braille language

Science.gov (United States)

Powell, R. A.

1967-01-01

Computer program converts print image tape files into six dot Braille cells, enabling a blind computer programmer to monitor and evaluate data generated by his own programs. The Braille output is printed 8 lines per inch.

Taylor's series method for solving the nonlinear point kinetics equations

International Nuclear Information System (INIS)

Nahla, Abdallah A.

2011-01-01

Highlights: → Taylor's series method for nonlinear point kinetics equations is applied. → The general order of derivatives are derived for this system. → Stability of Taylor's series method is studied. → Taylor's series method is A-stable for negative reactivity. → Taylor's series method is an accurate computational technique. - Abstract: Taylor's series method for solving the point reactor kinetics equations with multi-group of delayed neutrons in the presence of Newtonian temperature feedback reactivity is applied and programmed by FORTRAN. This system is the couples of the stiff nonlinear ordinary differential equations. This numerical method is based on the different order derivatives of the neutron density, the precursor concentrations of i-group of delayed neutrons and the reactivity. The r th order of derivatives are derived. The stability of Taylor's series method is discussed. Three sets of applications: step, ramp and temperature feedback reactivities are computed. Taylor's series method is an accurate computational technique and stable for negative step, negative ramp and temperature feedback reactivities. This method is useful than the traditional methods for solving the nonlinear point kinetics equations.

Introduction of the ASP3D Computer Program for Unsteady Aerodynamic and Aeroelastic Analyses

Science.gov (United States)

Batina, John T.

2005-01-01

A new computer program has been developed called ASP3D (Advanced Small Perturbation 3D), which solves the small perturbation potential flow equation in an advanced form including mass-consistent surface and trailing wake boundary conditions, and entropy, vorticity, and viscous effects. The purpose of the program is for unsteady aerodynamic and aeroelastic analyses, especially in the nonlinear transonic flight regime. The program exploits the simplicity of stationary Cartesian meshes with the movement or deformation of the configuration under consideration incorporated into the solution algorithm through a planar surface boundary condition. The new ASP3D code is the result of a decade of developmental work on improvements to the small perturbation formulation, performed while the author was employed as a Senior Research Scientist in the Configuration Aerodynamics Branch at the NASA Langley Research Center. The ASP3D code is a significant improvement to the state-of-the-art for transonic aeroelastic analyses over the CAP-TSD code (Computational Aeroelasticity Program Transonic Small Disturbance), which was developed principally by the author in the mid-1980s. The author is in a unique position as the developer of both computer programs to compare, contrast, and ultimately make conclusions regarding the underlying formulations and utility of each code. The paper describes the salient features of the ASP3D code including the rationale for improvements in comparison with CAP-TSD. Numerous results are presented to demonstrate the ASP3D capability. The general conclusion is that the new ASP3D capability is superior to the older CAP-TSD code because of the myriad improvements developed and incorporated.

General, Interactive Computer Program for the Solution of the Schrodinger Equation

Science.gov (United States)

Griffin, Donald C.; McGhie, James B.

1973-01-01

Discusses an interactive computer algorithm which allows beginning students to solve one- and three-dimensional quantum problems. Included is an example of the Thomas-Fermi-Dirac central field approximation. (CC)

Computer-Aided Corrosion Program Management

Science.gov (United States)

MacDowell, Louis

2010-01-01

This viewgraph presentation reviews Computer-Aided Corrosion Program Management at John F. Kennedy Space Center. The contents include: 1) Corrosion at the Kennedy Space Center (KSC); 2) Requirements and Objectives; 3) Program Description, Background and History; 4) Approach and Implementation; 5) Challenges; 6) Lessons Learned; 7) Successes and Benefits; and 8) Summary and Conclusions.

Numerical solution of conservation equations in the transient model for the system thermal - hydraulics in the Korsar computer code

International Nuclear Information System (INIS)

Yudov, Y.V.

2001-01-01

The functional part of the KORSAR computer code is based on the computational unit for the reactor system thermal-hydraulics and other thermal power systems with water cooling. The two-phase flow dynamics of the thermal-hydraulic network is modelled by KORSAR in one-dimensional two-fluid (non-equilibrium and nonhomogeneous) approximation with the same pressure of both phases. Each phase is characterized by parameters averaged over the channel sections, and described by the conservation equations for mass, energy and momentum. The KORSAR computer code relies upon a novel approach to mathematical modelling of two-phase dispersed-annular flows. This approach allows a two-fluid model to differentiate the effects of the liquid film and droplets in the gas core on the flow characteristics. A semi-implicit numerical scheme has been chosen for deriving discrete analogs the conservation equations in KORSAR. In the semi-implicit numerical scheme, solution of finite-difference equations is reduced to the problem of determining the pressure field at a new time level. For the one-channel case, the pressure field is found from the solution of a system of linear algebraic equations by using the tri-diagonal matrix method. In the branched network calculation, the matrix of coefficients in the equations describing the pressure field is no longer tri-diagonal but has a sparseness structure. In this case, the system of linear equations for the pressure field can be solved with any of the known classical methods. Such an approach is implemented in the existing best-estimate thermal-hydraulic computer codes (TRAC, RELAP5, etc.) For the KORSAR computer code, we have developed a new non-iterative method for calculating the pressure field in the network of any topology. This method is based on the tri-diagonal matrix method and performs well when solving the thermal-hydraulic network problems. (author)

An efficient computer based wavelets approximation method to solve Fuzzy boundary value differential equations

Science.gov (United States)

Alam Khan, Najeeb; Razzaq, Oyoon Abdul

2016-03-01

In the present work a wavelets approximation method is employed to solve fuzzy boundary value differential equations (FBVDEs). Essentially, a truncated Legendre wavelets series together with the Legendre wavelets operational matrix of derivative are utilized to convert FB- VDE into a simple computational problem by reducing it into a system of fuzzy algebraic linear equations. The capability of scheme is investigated on second order FB- VDE considered under generalized H-differentiability. Solutions are represented graphically showing competency and accuracy of this method.

New exact travelling wave solutions for two potential coupled KdV equations with symbolic computation

International Nuclear Information System (INIS)

Yang Zonghang

2007-01-01

We find new exact travelling wave solutions for two potential KdV equations which are presented by Foursov [Foursov MV. J Math Phys 2000;41:6173-85]. Compared with the extended tanh-function method, the algorithm used in our paper can obtain some new kinds of exact travelling wave solutions. With the aid of symbolic computation, some novel exact travelling wave solutions of the potential KdV equations are constructed

EDEF: a program for solving the neutron diffusion equation using microcomputers

International Nuclear Information System (INIS)

Fernandes, A.; Maiorino, J.R.

1990-01-01

This work presents the development of a program to solve the two-group two-dimensional diffusion equation (with a buckling option to simulate axial leakage) applying the finite element method. It has been developed to microcomputers compatibles to the IBM-PC. Among the facilities of the program, we can mention the simplicity to represent two-dimensional complex domains, the input through a pre-processor and the output in which the fluxes are presented graphically. The program also calculates the multiplication factor, the peaking factor and the power distribution. (author) [pt

Solution of the two-dimensional space-time reactor kinetics equation by a locally one-dimensional method

International Nuclear Information System (INIS)

Chen, G.S.; Christenson, J.M.

1985-01-01

In this paper, the authors present some initial results from an investigation of the application of a locally one-dimensional (LOD) finite difference method to the solution of the two-dimensional, two-group reactor kinetics equations. Although the LOD method is relatively well known, it apparently has not been previously applied to the space-time kinetics equations. In this investigation, the LOD results were benchmarked against similar computational results (using the same computing environment, the same programming structure, and the same sample problems) obtained by the TWIGL program. For all of the problems considered, the LOD method provided accurate results in one-half to one-eight of the time required by the TWIGL program

Differential equations, mechanics, and computation

CERN Document Server

Palais, Richard S

2009-01-01

This book provides a conceptual introduction to the theory of ordinary differential equations, concentrating on the initial value problem for equations of evolution and with applications to the calculus of variations and classical mechanics, along with a discussion of chaos theory and ecological models. It has a unified and visual introduction to the theory of numerical methods and a novel approach to the analysis of errors and stability of various numerical solution algorithms based on carefully chosen model problems. While the book would be suitable as a textbook for an undergraduate or elementary graduate course in ordinary differential equations, the authors have designed the text also to be useful for motivated students wishing to learn the material on their own or desiring to supplement an ODE textbook being used in a course they are taking with a text offering a more conceptual approach to the subject.

PyR@TE. Renormalization group equations for general gauge theories

Science.gov (United States)

Lyonnet, F.; Schienbein, I.; Staub, F.; Wingerter, A.

2014-03-01

Although the two-loop renormalization group equations for a general gauge field theory have been known for quite some time, deriving them for specific models has often been difficult in practice. This is mainly due to the fact that, albeit straightforward, the involved calculations are quite long, tedious and prone to error. The present work is an attempt to facilitate the practical use of the renormalization group equations in model building. To that end, we have developed two completely independent sets of programs written in Python and Mathematica, respectively. The Mathematica scripts will be part of an upcoming release of SARAH 4. The present article describes the collection of Python routines that we dubbed PyR@TE which is an acronym for “Python Renormalization group equations At Two-loop for Everyone”. In PyR@TE, once the user specifies the gauge group and the particle content of the model, the routines automatically generate the full two-loop renormalization group equations for all (dimensionless and dimensionful) parameters. The results can optionally be exported to LaTeX and Mathematica, or stored in a Python data structure for further processing by other programs. For ease of use, we have implemented an interactive mode for PyR@TE in form of an IPython Notebook. As a first application, we have generated with PyR@TE the renormalization group equations for several non-supersymmetric extensions of the Standard Model and found some discrepancies with the existing literature. Catalogue identifier: AERV_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AERV_v1_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 924959 No. of bytes in distributed program, including test data, etc.: 495197 Distribution format: tar.gz Programming language: Python. Computer

Choosing of mode and calculation of multiple regression equation parameters in X-ray radiometric analysis

International Nuclear Information System (INIS)

Mamikonyan, S.V.; Berezkin, V.V.; Lyubimova, S.V.; Svetajlo, Yu.N.; Shchekin, K.I.

1978-01-01

A method to derive multiple regression equations for X-ray radiometric analysis is described. Te method is realized in the form of the REGRA program in an algorithmic language. The subprograms included in the program are describe. In analyzing cement for Mg, Al, Si, Ca and Fe contents as an example, the obtainment of working equations in the course of calculations by the program is shown to simpliy the realization of computing devices in instruments for X-ray radiometric analysis

Symbolic computation of exact solutions expressible in rational formal hyperbolic and elliptic functions for nonlinear partial differential equations

International Nuclear Information System (INIS)

Wang Qi; Chen Yong

2007-01-01

With the aid of symbolic computation, some algorithms are presented for the rational expansion methods, which lead to closed-form solutions of nonlinear partial differential equations (PDEs). The new algorithms are given to find exact rational formal polynomial solutions of PDEs in terms of Jacobi elliptic functions, solutions of the Riccati equation and solutions of the generalized Riccati equation. They can be implemented in symbolic computation system Maple. As applications of the methods, we choose some nonlinear PDEs to illustrate the methods. As a result, we not only can successfully obtain the solutions found by most existing Jacobi elliptic function methods and Tanh-methods, but also find other new and more general solutions at the same time

PABLM: a computer program to calculate accumulated radiation doses from radionuclides in the environment

Energy Technology Data Exchange (ETDEWEB)

Napier, B.A.; Kennedy, W.E. Jr.; Soldat, J.K.

1980-03-01

A computer program, PABLM, was written to facilitate the calculation of internal radiation doses to man from radionuclides in food products and external radiation doses from radionuclides in the environment. This report contains details of mathematical models used and calculational procedures required to run the computer program. Radiation doses from radionuclides in the environment may be calculated from deposition on the soil or plants during an atmospheric or liquid release, or from exposure to residual radionuclides in the environment after the releases have ended. Radioactive decay is considered during the release of radionuclides, after they are deposited on the plants or ground, and during holdup of food after harvest. The radiation dose models consider several exposure pathways. Doses may be calculated for either a maximum-exposed individual or for a population group. The doses calculated are accumulated doses from continuous chronic exposure. A first-year committed dose is calculated as well as an integrated dose for a selected number of years. The equations for calculating internal radiation doses are derived from those given by the International Commission on Radiological Protection (ICRP) for body burdens and MPC's of each radionuclide. The radiation doses from external exposure to contaminated water and soil are calculated using the basic assumption that the contaminated medium is large enough to be considered an infinite volume or plane relative to the range of the emitted radiations. The equations for calculations of the radiation dose from external exposure to shoreline sediments include a correction for the finite width of the contaminated beach.

PABLM: a computer program to calculate accumulated radiation doses from radionuclides in the environment

International Nuclear Information System (INIS)

Napier, B.A.; Kennedy, W.E. Jr.; Soldat, J.K.

1980-03-01

A computer program, PABLM, was written to facilitate the calculation of internal radiation doses to man from radionuclides in food products and external radiation doses from radionuclides in the environment. This report contains details of mathematical models used and calculational procedures required to run the computer program. Radiation doses from radionuclides in the environment may be calculated from deposition on the soil or plants during an atmospheric or liquid release, or from exposure to residual radionuclides in the environment after the releases have ended. Radioactive decay is considered during the release of radionuclides, after they are deposited on the plants or ground, and during holdup of food after harvest. The radiation dose models consider several exposure pathways. Doses may be calculated for either a maximum-exposed individual or for a population group. The doses calculated are accumulated doses from continuous chronic exposure. A first-year committed dose is calculated as well as an integrated dose for a selected number of years. The equations for calculating internal radiation doses are derived from those given by the International Commission on Radiological Protection (ICRP) for body burdens and MPC's of each radionuclide. The radiation doses from external exposure to contaminated water and soil are calculated using the basic assumption that the contaminated medium is large enough to be considered an infinite volume or plane relative to the range of the emitted radiations. The equations for calculations of the radiation dose from external exposure to shoreline sediments include a correction for the finite width of the contaminated beach

a Recursive Approach to Compute Normal Forms

Science.gov (United States)

HSU, L.; MIN, L. J.; FAVRETTO, L.

2001-06-01

Normal forms are instrumental in the analysis of dynamical systems described by ordinary differential equations, particularly when singularities close to a bifurcation are to be characterized. However, the computation of a normal form up to an arbitrary order is numerically hard. This paper focuses on the computer programming of some recursive formulas developed earlier to compute higher order normal forms. A computer program to reduce the system to its normal form on a center manifold is developed using the Maple symbolic language. However, it should be stressed that the program relies essentially on recursive numerical computations, while symbolic calculations are used only for minor tasks. Some strategies are proposed to save computation time. Examples are presented to illustrate the application of the program to obtain high order normalization or to handle systems with large dimension.

Questioning the quantity equation using an agent-based computational model

DEFF Research Database (Denmark)

Bruun, Charlotte

2000-01-01

by Stutzel (1954), argues that the functional relationship may as well be negative. Even focusing the money needed to carry out transactions, there is no immediate answer to the question of the functional relationship between trade turnover and money demand. An agent-based computational model is used......In the literature we find two opposing hypotheses relating the volume of money to the volume of transactions or national income. The classic hypothesis, implicitly entailed in the quantity equation, argues that this relation must be positive, while an opposing hypothesis, most strongly presented...

Diffusive Wave Approximation to the Shallow Water Equations: Computational Approach

KAUST Repository

Collier, Nathan

2011-05-14

We discuss the use of time adaptivity applied to the one dimensional diffusive wave approximation to the shallow water equations. A simple and computationally economical error estimator is discussed which enables time-step size adaptivity. This robust adaptive time discretization corrects the initial time step size to achieve a user specified bound on the discretization error and allows time step size variations of several orders of magnitude. In particular, in the one dimensional results presented in this work feature a change of four orders of magnitudes for the time step over the entire simulation.

Description of mathematical models and computer programs

International Nuclear Information System (INIS)

1977-01-01

The paper gives a description of mathematical models and computer programs for analysing possible strategies for spent fuel management, with emphasis on economic analysis. The computer programs developed, describe the material flows, facility construction schedules, capital investment schedules and operating costs for the facilities used in managing the spent fuel. The computer programs use a combination of simulation and optimization procedures for the economic analyses. Many of the fuel cycle steps (such as spent fuel discharges, storage at the reactor, and transport to the RFCC) are described in physical and economic terms through simulation modeling, while others (such as reprocessing plant size and commissioning schedules, interim storage facility commissioning schedules etc.) are subjected to economic optimization procedures to determine the approximate lowest-cost plans from among the available feasible alternatives

Fluid dynamics computer programs for NERVA turbopump

Science.gov (United States)

Brunner, J. J.

1972-01-01

During the design of the NERVA turbopump, numerous computer programs were developed for the analyses of fluid dynamic problems within the machine. Program descriptions, example cases, users instructions, and listings for the majority of these programs are presented.

Advanced Certification Program for Computer Graphic Specialists. Final Performance Report.

Science.gov (United States)

Parkland Coll., Champaign, IL.

A pioneer program in computer graphics was implemented at Parkland College (Illinois) to meet the demand for specialized technicians to visualize data generated on high performance computers. In summer 1989, 23 students were accepted into the pilot program. Courses included C programming, calculus and analytic geometry, computer graphics, and…

DUKSUP: A Computer Program for High Thrust Launch Vehicle Trajectory Design and Optimization

Science.gov (United States)

Spurlock, O. Frank; Williams, Craig H.

2015-01-01

From the late 1960s through 1997, the leadership of NASAs Intermediate and Large class unmanned expendable launch vehicle projects resided at the NASA Lewis (now Glenn) Research Center (LeRC). One of LeRCs primary responsibilities --- trajectory design and performance analysis --- was accomplished by an internally-developed analytic three dimensional computer program called DUKSUP. Because of its Calculus of Variations-based optimization routine, this code was generally more capable of finding optimal solutions than its contemporaries. A derivation of optimal control using the Calculus of Variations is summarized including transversality, intermediate, and final conditions. The two point boundary value problem is explained. A brief summary of the codes operation is provided, including iteration via the Newton-Raphson scheme and integration of variational and motion equations via a 4th order Runge-Kutta scheme. Main subroutines are discussed. The history of the LeRC trajectory design efforts in the early 1960s is explained within the context of supporting the Centaur upper stage program. How the code was constructed based on the operation of the AtlasCentaur launch vehicle, the limits of the computers of that era, the limits of the computer programming languages, and the missions it supported are discussed. The vehicles DUKSUP supported (AtlasCentaur, TitanCentaur, and ShuttleCentaur) are briefly described. The types of missions, including Earth orbital and interplanetary, are described. The roles of flight constraints and their impact on launch operations are detailed (such as jettisoning hardware on heating, Range Safety, ground station tracking, and elliptical parking orbits). The computer main frames on which the code was hosted are described. The applications of the code are detailed, including independent check of contractor analysis, benchmarking, leading edge analysis, and vehicle performance improvement assessments. Several of DUKSUPs many major impacts on

A simplified model for computing equation of state of argon plasma

International Nuclear Information System (INIS)

Wang Caixia; Tian Yangmeng

2006-01-01

The paper present a simplified new model of computing equation of state and ionization degree of Argon plasma, which based on Thomas-Fermi (TF) statistical model: the authors fitted the numerical results of the ionization potential calculated by Thomas-Fermi statistical model and gained the analytical function of the potential versus the degree of ionization, then calculated the ionization potential and the average degree of ionization for Argon versus temperature and density in local thermal equilibrium case at 10-1000 eV. The results calculated of this simplified model are basically in agreement with several sets of theory data and experimental data. This simplified model can be used to calculation of the equation of state of plasmas mixture and is expected to have a more wide use in the field of EML technology involving the strongly ionized plasmas. (authors)

MOOSE: A parallel computational framework for coupled systems of nonlinear equations

International Nuclear Information System (INIS)

Gaston, Derek; Newman, Chris; Hansen, Glen; Lebrun-Grandie, Damien

2009-01-01

Systems of coupled, nonlinear partial differential equations (PDEs) often arise in simulation of nuclear processes. MOOSE: Multiphysics Object Oriented Simulation Environment, a parallel computational framework targeted at the solution of such systems, is presented. As opposed to traditional data-flow oriented computational frameworks, MOOSE is instead founded on the mathematical principle of Jacobian-free Newton-Krylov (JFNK). Utilizing the mathematical structure present in JFNK, physics expressions are modularized into 'Kernels,' allowing for rapid production of new simulation tools. In addition, systems are solved implicitly and fully coupled, employing physics-based preconditioning, which provides great flexibility even with large variance in time scales. A summary of the mathematics, an overview of the structure of MOOSE, and several representative solutions from applications built on the framework are presented.

Computer code qualification program for the Advanced CANDU Reactor

International Nuclear Information System (INIS)

Popov, N.K.; Wren, D.J.; Snell, V.G.; White, A.J.; Boczar, P.G.

2003-01-01

Atomic Energy of Canada Ltd (AECL) has developed and implemented a Software Quality Assurance program (SQA) to ensure that its analytical, scientific and design computer codes meet the required standards for software used in safety analyses. This paper provides an overview of the computer programs used in Advanced CANDU Reactor (ACR) safety analysis, and assessment of their applicability in the safety analyses of the ACR design. An outline of the incremental validation program, and an overview of the experimental program in support of the code validation are also presented. An outline of the SQA program used to qualify these computer codes is also briefly presented. To provide context to the differences in the SQA with respect to current CANDUs, the paper also provides an overview of the ACR design features that have an impact on the computer code qualification. (author)

Fermilab advanced computer program multi-microprocessor project

International Nuclear Information System (INIS)

Nash, T.; Areti, H.; Biel, J.

1985-06-01

Fermilab's Advanced Computer Program is constructing a powerful 128 node multi-microprocessor system for data analysis in high-energy physics. The system will use commercial 32-bit microprocessors programmed in Fortran-77. Extensive software supports easy migration of user applications from a uniprocessor environment to the multiprocessor and provides sophisticated program development, debugging, and error handling and recovery tools. This system is designed to be readily copied, providing computing cost effectiveness of below $2200 per VAX 11/780 equivalent. The low cost, commercial availability, compatibility with off-line analysis programs, and high data bandwidths (up to 160 MByte/sec) make the system an ideal choice for applications to on-line triggers as well as an offline data processor

The transition equation of the state intensities for exciton model and the calculation program

International Nuclear Information System (INIS)

Yu Xian; Zheng Jiwen; Liu Guoxing; Chen Keliang

1995-01-01

An equation set of the exciton model is given and calculation program is developed. The process of approaching to equilibrium state has been investigated with the program for 12 C + 64 Ni reaction at energy 72 MeV

The FORCE: A portable parallel programming language supporting computational structural mechanics

Science.gov (United States)

Jordan, Harry F.; Benten, Muhammad S.; Brehm, Juergen; Ramanan, Aruna

1989-01-01

This project supports the conversion of codes in Computational Structural Mechanics (CSM) to a parallel form which will efficiently exploit the computational power available from multiprocessors. The work is a part of a comprehensive, FORTRAN-based system to form a basis for a parallel version of the NICE/SPAR combination which will form the CSM Testbed. The software is macro-based and rests on the force methodology developed by the principal investigator in connection with an early scientific multiprocessor. Machine independence is an important characteristic of the system so that retargeting it to the Flex/32, or any other multiprocessor on which NICE/SPAR might be imnplemented, is well supported. The principal investigator has experience in producing parallel software for both full and sparse systems of linear equations using the force macros. Other researchers have used the Force in finite element programs. It has been possible to rapidly develop software which performs at maximum efficiency on a multiprocessor. The inherent machine independence of the system also means that the parallelization will not be limited to a specific multiprocessor.

A computer program for the estimation of time of death

DEFF Research Database (Denmark)

Lynnerup, N

1993-01-01

In the 1960s Marshall and Hoare presented a "Standard Cooling Curve" based on their mathematical analyses on the postmortem cooling of bodies. Although fairly accurate under standard conditions, the "curve" or formula is based on the assumption that the ambience temperature is constant and that t......In the 1960s Marshall and Hoare presented a "Standard Cooling Curve" based on their mathematical analyses on the postmortem cooling of bodies. Although fairly accurate under standard conditions, the "curve" or formula is based on the assumption that the ambience temperature is constant...... cooling of bodies is presented. It is proposed that by having a computer program that solves the equation, giving the length of the cooling period in response to a certain rectal temperature, and which allows easy comparison of multiple solutions, the uncertainties related to ambience temperature...

Finding all solutions of nonlinear equations using the dual simplex method

Science.gov (United States)

Yamamura, Kiyotaka; Fujioka, Tsuyoshi

2003-03-01

Recently, an efficient algorithm has been proposed for finding all solutions of systems of nonlinear equations using linear programming. This algorithm is based on a simple test (termed the LP test) for nonexistence of a solution to a system of nonlinear equations using the dual simplex method. In this letter, an improved version of the LP test algorithm is proposed. By numerical examples, it is shown that the proposed algorithm could find all solutions of a system of 300 nonlinear equations in practical computation time.

3D Relativistic Hydrodynamic Computations Using Lattice-QCD-Inspired Equations of State

International Nuclear Information System (INIS)

Hama, Yogiro; Andrade, Rone P.G.; Grassi, Frederique; Socolowski, Otavio; Kodama, Takeshi; Tavares, Bernardo; Padula, Sandra S.

2006-01-01

In this communication, we report results of three-dimensional hydrodynamic computations, by using equations of state with a critical end point as suggested by the lattice QCD. Some of the results are an increase of the multiplicity in the mid-rapidity region and a larger elliptic-flow parameter v 2 . We discuss also the effcts of the initial-condition fluctuations and the continuous emission

3D Relativistic Hydrodynamic Computations Using Lattice-QCD-Inspired Equations of State

Energy Technology Data Exchange (ETDEWEB)

Hama, Yogiro [Instituto de Fisica, Universidade de Sao Paulo (Brazil); Andrade, Rone P.G. [Instituto de Fisica, Universidade de Sao Paulo (Brazil); Grassi, Frederique [Instituto de Fisica, Universidade de Sao Paulo (Brazil); Socolowski, Otavio [Instituto Tecnologico da Aeronautica (Brazil); Kodama, Takeshi [Instituto de Fisica, Universidade Federal do Rio de Janeiro (Brazil); Tavares, Bernardo [Instituto de Fisica, Universidade Federal do Rio de Janeiro (Brazil); Padula, Sandra S. [Instituto de Fisica Teorica, Universidade Estadual Paulista (Brazil)

2006-08-07

In this communication, we report results of three-dimensional hydrodynamic computations, by using equations of state with a critical end point as suggested by the lattice QCD. Some of the results are an increase of the multiplicity in the mid-rapidity region and a larger elliptic-flow parameter v{sub 2}. We discuss also the effcts of the initial-condition fluctuations and the continuous emission.

The NASA computer science research program plan

Science.gov (United States)

1983-01-01

A taxonomy of computer science is included, one state of the art of each of the major computer science categories is summarized. A functional breakdown of NASA programs under Aeronautics R and D, space R and T, and institutional support is also included. These areas were assessed against the computer science categories. Concurrent processing, highly reliable computing, and information management are identified.

Case Studies of Liberal Arts Computer Science Programs

Science.gov (United States)

Baldwin, D.; Brady, A.; Danyluk, A.; Adams, J.; Lawrence, A.

2010-01-01

Many undergraduate liberal arts institutions offer computer science majors. This article illustrates how quality computer science programs can be realized in a wide variety of liberal arts settings by describing and contrasting the actual programs at five liberal arts colleges: Williams College, Kalamazoo College, the State University of New York…

Computer programming in the UK undergraduate mathematics curriculum

Science.gov (United States)

Sangwin, Christopher J.; O'Toole, Claire

2017-11-01

This paper reports a study which investigated the extent to which undergraduate mathematics students in the United Kingdom are currently taught to programme a computer as a core part of their mathematics degree programme. We undertook an online survey, with significant follow-up correspondence, to gather data on current curricula and received replies from 46 (63%) of the departments who teach a BSc mathematics degree. We found that 78% of BSc degree courses in mathematics included computer programming in a compulsory module but 11% of mathematics degree programmes do not teach programming to all their undergraduate mathematics students. In 2016, programming is most commonly taught to undergraduate mathematics students through imperative languages, notably MATLAB, using numerical analysis as the underlying (or parallel) mathematical subject matter. Statistics is a very popular choice in optional courses, using the package R. Computer algebra systems appear to be significantly less popular for compulsory first-year courses than a decade ago, and there was no mention of logic programming, functional programming or automatic theorem proving software. The modal form of assessment of computing modules is entirely by coursework (i.e. no examination).

Introduction of handheld computing to a family practice residency program.

Science.gov (United States)

Rao, Goutham

2002-01-01

Handheld computers are valuable practice tools. It is important for residency programs to introduce their trainees and faculty to this technology. This article describes a formal strategy to introduce handheld computing to a family practice residency program. Objectives were selected for the handheld computer training program that reflected skills physicians would find useful in practice. TRGpro handheld computers preloaded with a suite of medical reference programs, a medical calculator, and a database program were supplied to participants. Training consisted of four 1-hour modules each with a written evaluation quiz. Participants completed a self-assessment questionnaire after the program to determine their ability to meet each objective. Sixty of the 62 participants successfully completed the training program. The mean composite score on quizzes was 36 of 40 (90%), with no significant differences by level of residency training. The mean self-ratings of participants across all objectives was 3.31 of 4.00. Third-year residents had higher mean self-ratings than others (mean of group, 3.62). Participants were very comfortable with practical skills, such as using drug reference software, and less comfortable with theory, such as knowing the different types of handheld computers available. Structured training is a successful strategy for introducing handheld computing to a residency program.

Neuro-evolutionary computing paradigm for Painlevé equation-II in nonlinear optics

Science.gov (United States)

Ahmad, Iftikhar; Ahmad, Sufyan; Awais, Muhammad; Ul Islam Ahmad, Siraj; Asif Zahoor Raja, Muhammad

2018-05-01

The aim of this study is to investigate the numerical treatment of the Painlevé equation-II arising in physical models of nonlinear optics through artificial intelligence procedures by incorporating a single layer structure of neural networks optimized with genetic algorithms, sequential quadratic programming and active set techniques. We constructed a mathematical model for the nonlinear Painlevé equation-II with the help of networks by defining an error-based cost function in mean square sense. The performance of the proposed technique is validated through statistical analyses by means of the one-way ANOVA test conducted on a dataset generated by a large number of independent runs.

Systematic control of large computer programs

International Nuclear Information System (INIS)

Goedbloed, J.P.; Klieb, L.

1986-07-01

A package of CCL, UPDATE, and FORTRAN procedures is described which facilitates the systematic control and development of large scientific computer programs. The package provides a general tool box for this purpose which contains many conveniences for the systematic administration of files, editing, reformating of line printer output files, etc. In addition, a small number of procedures is devoted to the problem of structured development of a large computer program which is used by a group of scientists. The essence of the method is contained in three procedures N, R, and X for the creation of a new UPDATE program library, its revision, and execution, resp., and a procedure REVISE which provides a joint editor - UPDATE session which combines the advantages of the two systems, viz. speed and rigor. (Auth.)

Computer Aided Design System for Developing Musical Fountain Programs

Institute of Scientific and Technical Information of China (English)

刘丹; 张乃尧; 朱汉城

2003-01-01

A computer aided design system for developing musical fountain programs was developed with multiple functions such as intelligent design, 3-D animation, manual modification and synchronized motion to make the development process more efficient. The system first analyzed the music form and sentiment using many basic features of the music to select a basic fountain program. Then, this program is simulated with 3-D animation and modified manually to achieve the desired results. Finally, the program is transformed to a computer control program to control the musical fountain in time with the music. A prototype system for the musical fountain was also developed. It was tested with many styles of music and users were quite satisfied with its performance. By integrating various functions, the proposed computer aided design system for developing musical fountain programs greatly simplified the design of the musical fountain programs.

Subchannel analysis program for boiling water reactor fuel bundles based on five conservation equations of two-phase flow

International Nuclear Information System (INIS)

Bessho, Y.; Uchikawa, S.

1985-01-01

A subchannel analysis program, MENUETT, is developed for evaluation of thermal-hydraulic characteristics in boiling water reactor fuel bundles. This program is based on five conservation equations of two-phase flow with the drift-flux correlation. The cross flows are calculated separately for liquid and vapor phases from the lateral momentum conservation equation. The effects of turbulent mixing and void drift are accounted for in the program. The conservation equations are implicitly differentiated with the convective terms by the donor-cell method, and are solved iteratively in the axial and lateral directions. Data of the 3 X 3 rod bundle experiments are used for program verification. The lateral distributions of equilibrium quality and mass flow rate at the bundle exit calculated by the program compare satisfactorily with the experimental results

TRAFIC, a computer program for calculating the release of metallic fission products from an HTGR core

International Nuclear Information System (INIS)

Smith, P.D.

1978-02-01

A special purpose computer program, TRAFIC, is presented for calculating the release of metallic fission products from an HTGR core. The program is based upon Fick's law of diffusion for radioactive species. One-dimensional transient diffusion calculations are performed for the coated fuel particles and for the structural graphite web. A quasi steady-state calculation is performed for the fuel rod matrix material. The model accounts for nonlinear adsorption behavior in the fuel rod gap and on the coolant hole boundary. The TRAFIC program is designed to operate in a core survey mode; that is, it performs many repetitive calculations for a large number of spatial locations in the core. This is necessary in order to obtain an accurate volume integrated release. For this reason the program has been designed with calculational efficiency as one of its main objectives. A highly efficient numerical method is used in the solution. The method makes use of the Duhamel superposition principle to eliminate interior spatial solutions from consideration. Linear response functions relating the concentrations and mass fluxes on the boundaries of a homogeneous region are derived. Multiple regions are numerically coupled through interface conditions. Algebraic elimination is used to reduce the equations as far as possible. The problem reduces to two nonlinear equations in two unknowns, which are solved using a Newton Raphson technique

Computer Programs for Calculating and Plotting the Stability Characteristics of a Balloon Tethered in a Wind

Science.gov (United States)

Bennett, R. M.; Bland, S. R.; Redd, L. T.

1973-01-01

Computer programs for calculating the stability characteristics of a balloon tethered in a steady wind are presented. Equilibrium conditions, characteristic roots, and modal ratios are calculated for a range of discrete values of velocity for a fixed tether-line length. Separate programs are used: (1) to calculate longitudinal stability characteristics, (2) to calculate lateral stability characteristics, (3) to plot the characteristic roots versus velocity, (4) to plot the characteristic roots in root-locus form, (5) to plot the longitudinal modes of motion, and (6) to plot the lateral modes for motion. The basic equations, program listings, and the input and output data for sample cases are presented, with a brief discussion of the overall operation and limitations. The programs are based on a linearized, stability-derivative type of analysis, including balloon aerodynamics, apparent mass, buoyancy effects, and static forces which result from the tether line.

01010000 01001100 01000001 01011001: Play Elements in Computer Programming

Science.gov (United States)

Breslin, Samantha

2013-01-01

This article explores the role of play in human interaction with computers in the context of computer programming. The author considers many facets of programming including the literary practice of coding, the abstract design of programs, and more mundane activities such as testing, debugging, and hacking. She discusses how these incorporate the…

Computation of point reactor dynamics equations with thermal feedback via weighted residue method

International Nuclear Information System (INIS)

Suo Changan; Liu Xiaoming

1986-01-01

Point reactor dynamics equations with six groups of delayed neutrons have been computed via weighted-residual method in which the delta function was taken as a weighting function, and the parabolic with or without exponential factor as a trial function respectively for an insertion of large or smaller reactivity. The reactivity inserted into core can be varied with time, including insertion in forms of step function, polynomials up to second power and sine function. A thermal feedback of single flow channel model was added in. The thermal equations concerned were treated by use of a backward difference technique. A WRK code has been worked out, including implementation of an automatic selection of time span based on an input of error requirement and of an automatic change between computation with large reactivity and that with smaller one. On the condition of power varied slowly and without feedback, the results are not sensitive to the selection of values of time span. At last, the comparison of relevant results has shown that the agreement is quite well

Finite difference computing with exponential decay models

CERN Document Server

Langtangen, Hans Petter

2016-01-01

This text provides a very simple, initial introduction to the complete scientific computing pipeline: models, discretization, algorithms, programming, verification, and visualization. The pedagogical strategy is to use one case study – an ordinary differential equation describing exponential decay processes – to illustrate fundamental concepts in mathematics and computer science. The book is easy to read and only requires a command of one-variable calculus and some very basic knowledge about computer programming. Contrary to similar texts on numerical methods and programming, this text has a much stronger focus on implementation and teaches testing and software engineering in particular. .

An introduction to Python and computer programming

CERN Document Server

Zhang, Yue

2015-01-01

This book introduces Python programming language and fundamental concepts in algorithms and computing. Its target audience includes students and engineers with little or no background in programming, who need to master a practical programming language and learn the basic thinking in computer science/programming. The main contents come from lecture notes for engineering students from all disciplines, and has received high ratings. Its materials and ordering have been adjusted repeatedly according to classroom reception. Compared to alternative textbooks in the market, this book introduces the underlying Python implementation of number, string, list, tuple, dict, function, class, instance and module objects in a consistent and easy-to-understand way, making assignment, function definition, function call, mutability and binding environments understandable inside-out. By giving the abstraction of implementation mechanisms, this book builds a solid understanding of the Python programming language.

COBRA - 3C/KFKI: a digital computer program for steady and transient thermal-hydraulic analysis of rod bundle nuclear fuel elements

International Nuclear Information System (INIS)

Vigassy, J.; Kovacs, L.M.

1977-11-01

COBRA-3C/KFKI is a digital computer program for the CDC-3300 computer in FORTRAN language. The program is a revised version of the original COBRA-3C code. The code calculates steady-state and transient flow and enthalpy transport in rod-bundle nuclear fuel elements in both boiling and nonboiling conditions. The mathematical model is formulated by dividing the bundle flow area into flow subchannels that are assumed to contain one-dimensional flow and are coupled to each other by turbulent and diversion crossflow mixing. The program neglects sonic velocity propagation but allows for a temporal and spatial acceleration of the diversion crossflow in the transverse momentum equation. A semiexplicit finite-difference scheme is used to perform a boundary-value solution where the boundary conditions are the inlet enthalpy, inlet flow rate and exit pressure. (D.P.)

Parallel computing for homogeneous diffusion and transport equations in neutronics

International Nuclear Information System (INIS)

Pinchedez, K.

1999-06-01

Parallel computing meets the ever-increasing requirements for neutronic computer code speed and accuracy. In this work, two different approaches have been considered. We first parallelized the sequential algorithm used by the neutronics code CRONOS developed at the French Atomic Energy Commission. The algorithm computes the dominant eigenvalue associated with PN simplified transport equations by a mixed finite element method. Several parallel algorithms have been developed on distributed memory machines. The performances of the parallel algorithms have been studied experimentally by implementation on a T3D Cray and theoretically by complexity models. A comparison of various parallel algorithms has confirmed the chosen implementations. We next applied a domain sub-division technique to the two-group diffusion Eigen problem. In the modal synthesis-based method, the global spectrum is determined from the partial spectra associated with sub-domains. Then the Eigen problem is expanded on a family composed, on the one hand, from eigenfunctions associated with the sub-domains and, on the other hand, from functions corresponding to the contribution from the interface between the sub-domains. For a 2-D homogeneous core, this modal method has been validated and its accuracy has been measured. (author)

SAFE users manual. Volume 4. Computer programs

International Nuclear Information System (INIS)

Grady, L.M.

1983-06-01

Documentation for the Safeguards Automated Facility Evaluation (SAFE) computer programs is presented. The documentation is in the form of subprogram trees, program abstracts, flowcharts, and listings. Listings are provided on microfiche

Using trees to compute approximate solutions to ordinary differential equations exactly

Science.gov (United States)

Grossman, Robert

1991-01-01

Some recent work is reviewed which relates families of trees to symbolic algorithms for the exact computation of series which approximate solutions of ordinary differential equations. It turns out that the vector space whose basis is the set of finite, rooted trees carries a natural multiplication related to the composition of differential operators, making the space of trees an algebra. This algebraic structure can be exploited to yield a variety of algorithms for manipulating vector fields and the series and algebras they generate.

A Computer-Aided Writing Program for Learning Disabled Adolescents.

Science.gov (United States)

Fais, Laurie; Wanderman, Richard

The paper describes the application of a computer-assisted writing program in a special high school for learning disabled and dyslexic students and reports on a study of the program's effectiveness. Particular advantages of the Macintosh Computer for such a program are identified including use of the mouse pointing tool, graphic icons to identify…

Computer Assistance for Writing Interactive Programs: TICS.

Science.gov (United States)

Kaplow, Roy; And Others

1973-01-01

Investigators developed an on-line, interactive programing system--the Teacher-Interactive Computer System (TICS)--to provide assistance to those who were not programers, but nevertheless wished to write interactive instructional programs. TICS had two components: an author system and a delivery system. Underlying assumptions were that…

Scientific computing

CERN Document Server

Trangenstein, John A

2017-01-01

This is the third of three volumes providing a comprehensive presentation of the fundamentals of scientific computing. This volume discusses topics that depend more on calculus than linear algebra, in order to prepare the reader for solving differential equations. This book and its companions show how to determine the quality of computational results, and how to measure the relative efficiency of competing methods. Readers learn how to determine the maximum attainable accuracy of algorithms, and how to select the best method for computing problems. This book also discusses programming in several languages, including C++, Fortran and MATLAB. There are 90 examples, 200 exercises, 36 algorithms, 40 interactive JavaScript programs, 91 references to software programs and 1 case study. Topics are introduced with goals, literature references and links to public software. There are descriptions of the current algorithms in GSLIB and MATLAB. This book could be used for a second course in numerical methods, for either ...

Finite-element discretization of 3D energy-transport equations for semiconductors

Energy Technology Data Exchange (ETDEWEB)

Gadau, Stephan

2007-07-01

In this thesis a mathematical model was derived that describes the charge and energy transport in semiconductor devices like transistors. Moreover, numerical simulations of these physical processes are performed. In order to accomplish this, methods of theoretical physics, functional analysis, numerical mathematics and computer programming are applied. After an introduction to the status quo of semiconductor device simulation methods and a brief review of historical facts up to now, the attention is shifted to the construction of a model, which serves as the basis of the subsequent derivations in the thesis. Thereby the starting point is an important equation of the theory of dilute gases. From this equation the model equations are derived and specified by means of a series expansion method. This is done in a multi-stage derivation process, which is mainly taken from a scientific paper and which does not constitute the focus of this thesis. In the following phase we specify the mathematical setting and make precise the model assumptions. Thereby we make use of methods of functional analysis. Since the equations we deal with are coupled, we are concerned with a nonstandard problem. In contrary, the theory of scalar elliptic equations is established meanwhile. Subsequently, we are preoccupied with the numerical discretization of the equations. A special finite-element method is used for the discretization. This special approach has to be done in order to make the numerical results appropriate for practical application. By a series of transformations from the discrete model we derive a system of algebraic equations that are eligible for numerical evaluation. Using self-made computer programs we solve the equations to get approximate solutions. These programs are based on new and specialized iteration procedures that are developed and thoroughly tested within the frame of this research work. Due to their importance and their novel status, they are explained and

Mastering cloud computing foundations and applications programming

CERN Document Server

Buyya, Rajkumar; Selvi, SThamarai

2013-01-01

Mastering Cloud Computing is designed for undergraduate students learning to develop cloud computing applications. Tomorrow's applications won't live on a single computer but will be deployed from and reside on a virtual server, accessible anywhere, any time. Tomorrow's application developers need to understand the requirements of building apps for these virtual systems, including concurrent programming, high-performance computing, and data-intensive systems. The book introduces the principles of distributed and parallel computing underlying cloud architectures and specifical

Essential partial differential equations analytical and computational aspects

CERN Document Server

Griffiths, David F; Silvester, David J

2015-01-01

This volume provides an introduction to the analytical and numerical aspects of partial differential equations (PDEs). It unifies an analytical and computational approach for these; the qualitative behaviour of solutions being established using classical concepts: maximum principles and energy methods.   Notable inclusions are the treatment of irregularly shaped boundaries, polar coordinates and the use of flux-limiters when approximating hyperbolic conservation laws. The numerical analysis of difference schemes is rigorously developed using discrete maximum principles and discrete Fourier analysis. A novel feature is the inclusion of a chapter containing projects, intended for either individual or group study, that cover a range of topics such as parabolic smoothing, travelling waves, isospectral matrices, and the approximation of multidimensional advection–diffusion problems.   The underlying theory is illustrated by numerous examples and there are around 300 exercises, designed to promote and test unde...

The Dynamic Geometrisation of Computer Programming

Science.gov (United States)

Sinclair, Nathalie; Patterson, Margaret

2018-01-01

The goal of this paper is to explore dynamic geometry environments (DGE) as a type of computer programming language. Using projects created by secondary students in one particular DGE, we analyse the extent to which the various aspects of computational thinking--including both ways of doing things and particular concepts--were evident in their…

Computer Series, 38.

Science.gov (United States)

Moore, John W., Ed.

1983-01-01

Discusses numerical solution of the one-dimension Schrodinger equation. A PASCAL computer program for the Apple II which performs the calculations is available from the authors. Also discusses quantization and perturbation theory using microcomputers, indicating benefits of using the addition of a perturbation term to harmonic oscillator as an…

The Harwell TAILS computer program user's manual

International Nuclear Information System (INIS)

Rouse, K.D.; Cooper, M.J.

1980-11-01

The Harwell TAILS computer program is a versatile program for crystal structure refinement through the analysis of neutron or X-ray diffraction data from single crystals or powders. The main features of the program are described and details are given of the data input and output specifications. (author)

The Computational Physics Program of the national MFE Computer Center

International Nuclear Information System (INIS)

Mirin, A.A.

1989-01-01

Since June 1974, the MFE Computer Center has been engaged in a significant computational physics effort. The principal objective of the Computational Physics Group is to develop advanced numerical models for the investigation of plasma phenomena and the simulation of present and future magnetic confinement devices. Another major objective of the group is to develop efficient algorithms and programming techniques for current and future generations of supercomputers. The Computational Physics Group has been involved in several areas of fusion research. One main area is the application of Fokker-Planck/quasilinear codes to tokamaks. Another major area is the investigation of resistive magnetohydrodynamics in three dimensions, with applications to tokamaks and compact toroids. A third area is the investigation of kinetic instabilities using a 3-D particle code; this work is often coupled with the task of numerically generating equilibria which model experimental devices. Ways to apply statistical closure approximations to study tokamak-edge plasma turbulence have been under examination, with the hope of being able to explain anomalous transport. Also, we are collaborating in an international effort to evaluate fully three-dimensional linear stability of toroidal devices. In addition to these computational physics studies, the group has developed a number of linear systems solvers for general classes of physics problems and has been making a major effort at ascertaining how to efficiently utilize multiprocessor computers. A summary of these programs are included in this paper. 6 tabs

Numerical solution of neutron transport equations in discrete ordinates and slab geometry

International Nuclear Information System (INIS)

Serrano Pedraza, F.

1985-01-01

An unified formalism to solve numerically, between other equation, the neutron transport in discrete ordinates, slab geometry, several energy groups and independents of time, has been developed recently. Such a formalism cover some of the conventional schemes as diamond difference, (WDD) characteristic step (SC) lineal characteristic (LC), quadratic characteristic (QC) and lineal discontinuous. Unified formation gives before hand the convergence order of the previously selected scheme. In fact it allows besides to generate a big amount of numerical schemes, with which is also possible to solve numerical equations as soon as neutron transport. The essential purpose of this work was to solve the neutron transport equations in slab geometry and discrete ordinates considering several energy groups without to take under advisement time dependence based in the above mentioned unified formalism. To reach this purpose it was necesary to design a computer code with the name TNOD1 (Neutron transport in discrete ordinates and 1 dimension) which includes each one of the schemes already pointed out. there exist two numerical schemes, also recently developed, quadratic continuous (QC) and cubic continuous (CN), although covered by unified formalism, it has been possible to include them inside this computer code without make substantial changes in its structure. In chapter I, derivative of neutron transport equation independent of time is taken, for angular flux, including boundary conditions and discontinuity. In chapter II the neutron transport equations are obtained in multigroups, independents of time, for approximation of discrete ordinates. Description of theory related with unified formalism and its relationship with mentioned discretization schemes is presented in chapter III. Chapter IV describes the computer code developed and finally, in chapter V different numerical results obtained with TNOD1 program are shown. In Appendix A theorems and mathematical arguments used

Thinking processes used by high-performing students in a computer programming task

Directory of Open Access Journals (Sweden)

Marietjie Havenga

2011-07-01

Full Text Available Computer programmers must be able to understand programming source code and write programs that execute complex tasks to solve real-world problems. This article is a trans- disciplinary study at the intersection of computer programming, education and psychology. It outlines the role of mental processes in the process of programming and indicates how successful thinking processes can support computer science students in writing correct and well-defined programs. A mixed methods approach was used to better understand the thinking activities and programming processes of participating students. Data collection involved both computer programs and students’ reflective thinking processes recorded in their journals. This enabled analysis of psychological dimensions of participants’ thinking processes and their problem-solving activities as they considered a programming problem. Findings indicate that the cognitive, reflective and psychological processes used by high-performing programmers contributed to their success in solving a complex programming problem. Based on the thinking processes of high performers, we propose a model of integrated thinking processes, which can support computer programming students. Keywords: Computer programming, education, mixed methods research, thinking processes.  Disciplines: Computer programming, education, psychology

Computation of high Reynolds number internal/external flows

International Nuclear Information System (INIS)

Cline, M.C.; Wilmoth, R.G.

1981-01-01

A general, user oriented computer program, called VNAP2, has been developed to calculate high Reynolds number, internal/external flows. VNAP2 solves the two-dimensional, time-dependent Navier-Stokes equations. The turbulence is modeled with either a mixing-length, a one transport equation, or a two transport equation model. Interior grid points are computed using the explicit MacCormack scheme with special procedures to speed up the calculation in the fine grid. All boundary conditions are calculated using a reference plane characteristic scheme with the viscous terms treated as source terms. Several internal, external, and internal/external flow calculations are presented

Structured Parallel Programming Patterns for Efficient Computation

CERN Document Server

McCool, Michael; Robison, Arch

2012-01-01

Programming is now parallel programming. Much as structured programming revolutionized traditional serial programming decades ago, a new kind of structured programming, based on patterns, is relevant to parallel programming today. Parallel computing experts and industry insiders Michael McCool, Arch Robison, and James Reinders describe how to design and implement maintainable and efficient parallel algorithms using a pattern-based approach. They present both theory and practice, and give detailed concrete examples using multiple programming models. Examples are primarily given using two of th

Software survey: VOSviewer, a computer program for bibliometric mapping

NARCIS (Netherlands)

N.J.P. van Eck (Nees Jan); L. Waltman (Ludo)

2010-01-01

textabstractWe present VOSviewer, a freely available computer program that we have developed for constructing and viewing bibliometric maps. Unlike most computer programs that are used for bibliometric mapping, VOSviewer pays special attention to the graphical representation of bibliometric maps.

Method and computer program product for maintenance and modernization backlogging

Science.gov (United States)

Mattimore, Bernard G; Reynolds, Paul E; Farrell, Jill M

2013-02-19

According to one embodiment, a computer program product for determining future facility conditions includes a computer readable medium having computer readable program code stored therein. The computer readable program code includes computer readable program code for calculating a time period specific maintenance cost, for calculating a time period specific modernization factor, and for calculating a time period specific backlog factor. Future facility conditions equal the time period specific maintenance cost plus the time period specific modernization factor plus the time period specific backlog factor. In another embodiment, a computer-implemented method for calculating future facility conditions includes calculating a time period specific maintenance cost, calculating a time period specific modernization factor, and calculating a time period specific backlog factor. Future facility conditions equal the time period specific maintenance cost plus the time period specific modernization factor plus the time period specific backlog factor. Other embodiments are also presented.

hypoDD-A Program to Compute Double-Difference Hypocenter Locations

Science.gov (United States)

Waldhauser, Felix

2001-01-01

HypoDD is a Fortran computer program package for relocating earthquakes with the double-difference algorithm of Waldhauser and Ellsworth (2000). This document provides a brief introduction into how to run and use the programs ph2dt and hypoDD to compute double-difference (DD) hypocenter locations. It gives a short overview of the DD technique, discusses the data preprocessing using ph2dt, and leads through the earthquake relocation process using hypoDD. The appendices include the reference manuals for the two programs and a short description of auxiliary programs and example data. Some minor subroutines are presently in the c language, and future releases will be in c. Earthquake location algorithms are usually based on some form of Geiger’s method, the linearization of the travel time equation in a first order Taylor series that relates the difference between the observed and predicted travel time to unknown adjustments in the hypocentral coordinates through the partial derivatives of travel time with respect to the unknowns. Earthquakes can be located individually with this algorithm, or jointly when other unknowns link together the solutions to indivdual earthquakes, such as station corrections in the joint hypocenter determination (JHD) method, or the earth model in seismic tomography. The DD technique (described in detail in Waldhauser and Ellsworth, 2000) takes advantage of the fact that if the hypocentral separation between two earthquakes is small compared to the event-station distance and the scale length of velocity heterogeneity, then the ray paths between the source region and a common station are similar along almost the entire ray path (Fréchet, 1985; Got et al., 1994). In this case, the difference in travel times for two events observed at one station can be attributed to the spatial offset between the events with high accuracy. DD equations are built by differencing Geiger’s equation for earthquake location. In this way, the residual between

Bridging the Knowledge Gaps between Richards' Equation and Budyko Equation

Science.gov (United States)

Wang, D.

2017-12-01

The empirical Budyko equation represents the partitioning of mean annual precipitation into evaporation and runoff. Richards' equation, based on Darcy's law, represents the movement of water in unsaturated soils. The linkage between Richards' equation and Budyko equation is presented by invoking the empirical Soil Conservation Service curve number (SCS-CN) model for computing surface runoff at the event-scale. The basis of the SCS-CN method is the proportionality relationship, i.e., the ratio of continuing abstraction to its potential is equal to the ratio of surface runoff to its potential value. The proportionality relationship can be derived from the Richards' equation for computing infiltration excess and saturation excess models at the catchment scale. Meanwhile, the generalized proportionality relationship is demonstrated as the common basis of SCS-CN method, monthly "abcd" model, and Budyko equation. Therefore, the linkage between Darcy's law and the emergent pattern of mean annual water balance at the catchment scale is presented through the proportionality relationship.

Computer programs for display. [magnetic tapes - project planning/NASA programs

Science.gov (United States)

1975-01-01

The developments of an information storage and retrieval system are presented. Computer programs used in the system are described; the programs allow display messages to be placed on disks in an off-line environment permitting a more efficient use of memory. A time table that shows complete and scheduled developments of the system is given.

The NEA computer program library: a possible GDMS application

International Nuclear Information System (INIS)

Schuler, W.

1978-01-01

NEA Computer Program library maintains a series of eleven sequential computer files, used for linked applications in managing their stock of computer codes for nuclear reactor calculations, storing index and program abstract information, and administering their service to requesters. The high data redundancy beween the files suggests that a data base approach would be valid and this paper suggests a possible 'schema' for an CODASYL GDMS

78 FR 15734 - Privacy Act of 1974; Computer Matching Program

Science.gov (United States)

2013-03-12

... 1974; Computer Matching Program AGENCY: Department of Homeland Security/U.S. Citizenship and... computer matching program between the Department of Homeland Security/U.S. Citizenship and Immigration... Protection Act of 1988 (Pub. L. 100-503) and the Computer Matching and Privacy Protection Amendments of 1990...

78 FR 38724 - Privacy Act of 1974; Computer Matching Program

Science.gov (United States)

2013-06-27

... 1974; Computer Matching Program AGENCY: Department of Homeland Security/U.S. Citizenship and... Agreement that establishes a computer matching program between the Department of Homeland Security/U.S... and Privacy Protection Act of 1988 (Pub. L. 100-503) and the Computer Matching and Privacy Protection...

78 FR 15733 - Privacy Act of 1974; Computer Matching Program

Science.gov (United States)

2013-03-12

... 1974; Computer Matching Program AGENCY: Department of Homeland Security/U.S. Citizenship and... computer matching program between the Department of Homeland Security/U.S. Citizenship and Immigration... Protection Act of 1988 (Pub. L. 100-503) and the Computer Matching and Privacy Protection Amendments of 1990...

User's manual for computer program BASEPLOT

Science.gov (United States)

Sanders, Curtis L.

2002-01-01

The checking and reviewing of daily records of streamflow within the U.S. Geological Survey is traditionally accomplished by hand-plotting and mentally collating tables of data. The process is time consuming, difficult to standardize, and subject to errors in computation, data entry, and logic. In addition, the presentation of flow data on the internet requires more timely and accurate computation of daily flow records. BASEPLOT was developed for checking and review of primary streamflow records within the U.S. Geological Survey. Use of BASEPLOT enables users to (1) provide efficiencies during the record checking and review process, (2) improve quality control, (3) achieve uniformity of checking and review techniques of simple stage-discharge relations, and (4) provide a tool for teaching streamflow computation techniques. The BASEPLOT program produces tables of quality control checks and produces plots of rating curves and discharge measurements; variable shift (V-shift) diagrams; and V-shifts converted to stage-discharge plots, using data stored in the U.S. Geological Survey Automatic Data Processing System database. In addition, the program plots unit-value hydrographs that show unit-value stages, shifts, and datum corrections; input shifts, datum corrections, and effective dates; discharge measurements; effective dates for rating tables; and numeric quality control checks. Checklist/tutorial forms are provided for reviewers to ensure completeness of review and standardize the review process. The program was written for the U.S. Geological Survey SUN computer using the Statistical Analysis System (SAS) software produced by SAS Institute, Incorporated.

EZLP: An Interactive Computer Program for Solving Linear Programming Problems. Final Report.

Science.gov (United States)

Jarvis, John J.; And Others

Designed for student use in solving linear programming problems, the interactive computer program described (EZLP) permits the student to input the linear programming model in exactly the same manner in which it would be written on paper. This report includes a brief review of the development of EZLP; narrative descriptions of program features,…

HAL/SM language specification. [programming languages and computer programming for space shuttles

Science.gov (United States)

Williams, G. P. W., Jr.; Ross, C.

1975-01-01

A programming language is presented for the flight software of the NASA Space Shuttle program. It is intended to satisfy virtually all of the flight software requirements of the space shuttle. To achieve this, it incorporates a wide range of features, including applications-oriented data types and organizations, real time control mechanisms, and constructs for systems programming tasks. It is a higher order language designed to allow programmers, analysts, and engineers to communicate with the computer in a form approximating natural mathematical expression. Parts of the English language are combined with standard notation to provide a tool that readily encourages programming without demanding computer hardware expertise. Block diagrams and flow charts are included. The semantics of the language is discussed.

Direct integration of the inverse Radon equation for X-ray computed tomography.

Science.gov (United States)

Libin, E E; Chakhlov, S V; Trinca, D

2016-11-22

A new mathematical appoach using the inverse Radon equation for restoration of images in problems of linear two-dimensional x-ray tomography is formulated. In this approach, Fourier transformation is not used, and it gives the chance to create the practical computing algorithms having more reliable mathematical substantiation. Results of software implementation show that for especially for low number of projections, the described approach performs better than standard X-ray tomographic reconstruction algorithms.

Auxiliary equation method for solving nonlinear partial differential equations

International Nuclear Information System (INIS)

Sirendaoreji,; Jiong, Sun

2003-01-01

By using the solutions of an auxiliary ordinary differential equation, a direct algebraic method is described to construct several kinds of exact travelling wave solutions for some nonlinear partial differential equations. By this method some physically important nonlinear equations are investigated and new exact travelling wave solutions are explicitly obtained with the aid of symbolic computation

A new computational method for simulation of charge transport in semiconductor radiation detectors

International Nuclear Information System (INIS)

Holban, I.

1993-01-01

An effective computational method for simulation of charge transport in semiconductor radiation detectors is the purpose of the present work. Basic equations for analysis include (1) Poisson's equations, (2) continuity equation for electrons and holes, (3) rate equations for deep levels, (4) current equation for electrons and holes and (5) boundary conditions. The system of equations is discretized and equidistant space and time grids is brought. The nonlinearity of the problem is overcome by using Newton-Raphson iteration scheme. Instead of solving a nonlinear boundary problem we resolve a linear matrix equation. Our computation procedure becomes very efficient using a sparse matrix. The computed program allows to calculate the charge collection efficiency and transient response for arbitrary electric fields when trapping and detrapping effects are present. The earlier literature results are reproduced. (Author)

Evaluation of a Digital Game-Based Learning Program for Enhancing Youth Mental Health: A Structural Equation Modeling of the Program Effectiveness.

Science.gov (United States)

Huen, Jenny My; Lai, Eliza Sy; Shum, Angie Ky; So, Sam Wk; Chan, Melissa Ky; Wong, Paul Wc; Law, Y W; Yip, Paul Sf

2016-10-07

Digital game-based learning (DGBL) makes use of the entertaining power of digital games for educational purposes. Effectiveness assessment of DGBL programs has been underexplored and no attempt has been made to simultaneously model both important components of DGBL: learning attainment (ie, educational purposes of DGBL) and engagement of users (ie, entertaining power of DGBL) in evaluating program effectiveness. This study aimed to describe and evaluate an Internet-based DGBL program, Professor Gooley and the Flame of Mind, which promotes mental health to adolescents in a positive youth development approach. In particular, we investigated whether user engagement in the DGBL program could enhance their attainment on each of the learning constructs per DGBL module and subsequently enhance their mental health as measured by psychological well-being. Users were assessed on their attainment on each learning construct, psychological well-being, and engagement in each of the modules. One structural equation model was constructed for each DGBL module to model the effect of users' engagement and attainment on the learning construct on their psychological well-being. Of the 498 secondary school students that registered and participated from the first module of the DGBL program, 192 completed all 8 modules of the program. Results from structural equation modeling suggested that a higher extent of engagement in the program activities facilitated users' attainment on the learning constructs on most of the modules and in turn enhanced their psychological well-being after controlling for users' initial psychological well-being and initial attainment on the constructs. This study provided evidence that Internet intervention for mental health, implemented with the technologies and digital innovations of DGBL, could enhance youth mental health. Structural equation modeling is a promising approach in evaluating the effectiveness of DGBL programs.

Navier-Stokes equations by the finite element method

International Nuclear Information System (INIS)

Portella, P.E.

1984-01-01

A computer program to solve the Navier-Stokes equations by using the Finite Element Method is implemented. The solutions variables investigated are stream-function/vorticity in the steady case and velocity/pressure in the steady state and transient cases. For steady state flow the equations are solved simultaneously by the Newton-Raphson method. For the time dependent formulation, a fractional step method is employed to discretize in time and artificial viscosity is used to preclude spurious oscilations in the solution. The element used is the three node triangle. Some numerical examples are presented and comparisons are made with applications already existent. (Author) [pt

HEDPIN: a computer program to estimate pinwise power density

International Nuclear Information System (INIS)

Cappiello, M.W.

1976-05-01

A description is given of the digital computer program, HEDPIN. This program, modeled after a previously developed program, POWPIN, provides a means of estimating the pinwise power density distribution in fast reactor triangular pitched pin bundles. The capability also exists for computing any reaction rate of interest at the respective pin positions within an assembly. HEDPIN was developed in support of FTR fuel and test management as well as fast reactor core design and core characterization planning and analysis. The results of a test devised to check out HEDPIN's computational method are given, and the realm of application is discussed. Nearly all programming is in FORTRAN IV. Variable dimensioning is employed to make efficient use of core memory and maintain short running time for small problems. Input instructions, sample problem, and a program listing are also given

Computer programs in accelerator physics

International Nuclear Information System (INIS)

Keil, E.

1984-01-01

Three areas of accelerator physics are discussed in which computer programs have been applied with much success: i) single-particle beam dynamics in circular machines, i.e. the design and matching of machine lattices; ii) computations of electromagnetic fields in RF cavities and similar objects, useful for the design of RF cavities and for the calculation of wake fields; iii) simulation of betatron and synchrotron oscillations in a machine with non-linear elements, e.g. sextupoles, and of bunch lengthening due to longitudinal wake fields. (orig.)

On Solutions for Linear and Nonlinear Schrödinger Equations with Variable Coefficients: A Computational Approach

Directory of Open Access Journals (Sweden)

Gabriel Amador

2016-05-01

Full Text Available In this work, after reviewing two different ways to solve Riccati systems, we are able to present an extensive list of families of integrable nonlinear Schrödinger (NLS equations with variable coefficients. Using Riccati equations and similarity transformations, we are able to reduce them to the standard NLS models. Consequently, we can construct bright-, dark- and Peregrine-type soliton solutions for NLS with variable coefficients. As an important application of solutions for the Riccati equation with parameters, by means of computer algebra systems, it is shown that the parameters change the dynamics of the solutions. Finally, we test numerical approximations for the inhomogeneous paraxial wave equation by the Crank-Nicolson scheme with analytical solutions found using Riccati systems. These solutions include oscillating laser beams and Laguerre and Gaussian beams.

New Mexico district work-effort analysis computer program

Science.gov (United States)

Hiss, W.L.; Trantolo, A.P.; Sparks, J.L.

1972-01-01

The computer program (CAN 2) described in this report is one of several related programs used in the New Mexico District cost-analysis system. The work-effort information used in these programs is accumulated and entered to the nearest hour on forms completed by each employee. Tabulating cards are punched directly from these forms after visual examinations for errors are made. Reports containing detailed work-effort data itemized by employee within each project and account and by account and project for each employee are prepared for both current-month and year-to-date periods by the CAN 2 computer program. An option allowing preparation of reports for a specified 3-month period is provided. The total number of hours worked on each account and project and a grand total of hours worked in the New Mexico District is computed and presented in a summary report for each period. Work effort not chargeable directly to individual projects or accounts is considered as overhead and can be apportioned to the individual accounts and projects on the basis of the ratio of the total hours of work effort for the individual accounts or projects to the total New Mexico District work effort at the option of the user. The hours of work performed by a particular section, such as General Investigations or Surface Water, are prorated and charged to the projects or accounts within the particular section. A number of surveillance or buffer accounts are employed to account for the hours worked on special events or on those parts of large projects or accounts that require a more detailed analysis. Any part of the New Mexico District operation can be separated and analyzed in detail by establishing an appropriate buffer account. With the exception of statements associated with word size, the computer program is written in FORTRAN IV in a relatively low and standard language level to facilitate its use on different digital computers. The program has been run only on a Control Data Corporation

Injecting Artificial Memory Errors Into a Running Computer Program

Science.gov (United States)

Bornstein, Benjamin J.; Granat, Robert A.; Wagstaff, Kiri L.

2008-01-01

Single-event upsets (SEUs) or bitflips are computer memory errors caused by radiation. BITFLIPS (Basic Instrumentation Tool for Fault Localized Injection of Probabilistic SEUs) is a computer program that deliberately injects SEUs into another computer program, while the latter is running, for the purpose of evaluating the fault tolerance of that program. BITFLIPS was written as a plug-in extension of the open-source Valgrind debugging and profiling software. BITFLIPS can inject SEUs into any program that can be run on the Linux operating system, without needing to modify the program s source code. Further, if access to the original program source code is available, BITFLIPS offers fine-grained control over exactly when and which areas of memory (as specified via program variables) will be subjected to SEUs. The rate of injection of SEUs is controlled by specifying either a fault probability or a fault rate based on memory size and radiation exposure time, in units of SEUs per byte per second. BITFLIPS can also log each SEU that it injects and, if program source code is available, report the magnitude of effect of the SEU on a floating-point value or other program variable.

Solution of the Schrodinger Equation for One-Dimensional Anharmonic Potentials: An Undergraduate Computational Experiment

Science.gov (United States)

Beddard, Godfrey S.

2011-01-01

A method of solving the Schrodinger equation using a basis set expansion is described and used to calculate energy levels and wavefunctions of the hindered rotation of ethane and the ring puckering of cyclopentene. The calculations were performed using a computer algebra package and the calculations are straightforward enough for undergraduates to…

The Outlook for Computer Professions: 1985 Rewrites the Program.

Science.gov (United States)

Drake, Larry

1986-01-01

The author states that graduates of junior college programs who learn COBOL will continue to find jobs, but employers will increasingly seek college graduates when filling positions for computer programers and systems analysts. Areas of growth for computer applications (services, military, data communications, and artificial intelligence) are…

Near-Surface Seismic Velocity Data: A Computer Program For ...

African Journals Online (AJOL)

A computer program (NESURVELANA) has been developed in Visual Basic Computer programming language to carry out a near surface velocity analysis. The method of analysis used includes: Algorithms design and Visual Basic codes generation for plotting arrival time (ms) against geophone depth (m) employing the ...

Computational Cellular Dynamics Based on the Chemical Master Equation: A Challenge for Understanding Complexity.

Science.gov (United States)

Liang, Jie; Qian, Hong

2010-01-01

Modern molecular biology has always been a great source of inspiration for computational science. Half a century ago, the challenge from understanding macromolecular dynamics has led the way for computations to be part of the tool set to study molecular biology. Twenty-five years ago, the demand from genome science has inspired an entire generation of computer scientists with an interest in discrete mathematics to join the field that is now called bioinformatics. In this paper, we shall lay out a new mathematical theory for dynamics of biochemical reaction systems in a small volume (i.e., mesoscopic) in terms of a stochastic, discrete-state continuous-time formulation, called the chemical master equation (CME). Similar to the wavefunction in quantum mechanics, the dynamically changing probability landscape associated with the state space provides a fundamental characterization of the biochemical reaction system. The stochastic trajectories of the dynamics are best known through the simulations using the Gillespie algorithm. In contrast to the Metropolis algorithm, this Monte Carlo sampling technique does not follow a process with detailed balance. We shall show several examples how CMEs are used to model cellular biochemical systems. We shall also illustrate the computational challenges involved: multiscale phenomena, the interplay between stochasticity and nonlinearity, and how macroscopic determinism arises from mesoscopic dynamics. We point out recent advances in computing solutions to the CME, including exact solution of the steady state landscape and stochastic differential equations that offer alternatives to the Gilespie algorithm. We argue that the CME is an ideal system from which one can learn to understand "complex behavior" and complexity theory, and from which important biological insight can be gained.

Stochastic optimal control in infinite dimension dynamic programming and HJB equations

CERN Document Server

Fabbri, Giorgio; Święch, Andrzej

2017-01-01

Providing an introduction to stochastic optimal control in infinite dimension, this book gives a complete account of the theory of second-order HJB equations in infinite-dimensional Hilbert spaces, focusing on its applicability to associated stochastic optimal control problems. It features a general introduction to optimal stochastic control, including basic results (e.g. the dynamic programming principle) with proofs, and provides examples of applications. A complete and up-to-date exposition of the existing theory of viscosity solutions and regular solutions of second-order HJB equations in Hilbert spaces is given, together with an extensive survey of other methods, with a full bibliography. In particular, Chapter 6, written by M. Fuhrman and G. Tessitore, surveys the theory of regular solutions of HJB equations arising in infinite-dimensional stochastic control, via BSDEs. The book is of interest to both pure and applied researchers working in the control theory of stochastic PDEs, and in PDEs in infinite ...

A REDUCE program for symbolic computation of Puiseux expansions

International Nuclear Information System (INIS)

Gerdt, V.P.; Tiller, P.

1991-01-01

The program is described for computation of Puiseux expansions of algebraic functions. The Newton diagram method is used for construction of initial coefficients of all the Puiseux series at the given point. The program is written in computer algebra language Reduce. Some illustrative examples are given. 20 refs

Compiler Technology for Parallel Scientific Computation

Directory of Open Access Journals (Sweden)

Can Özturan

1994-01-01

Full Text Available There is a need for compiler technology that, given the source program, will generate efficient parallel codes for different architectures with minimal user involvement. Parallel computation is becoming indispensable in solving large-scale problems in science and engineering. Yet, the use of parallel computation is limited by the high costs of developing the needed software. To overcome this difficulty we advocate a comprehensive approach to the development of scalable architecture-independent software for scientific computation based on our experience with equational programming language (EPL. Our approach is based on a program decomposition, parallel code synthesis, and run-time support for parallel scientific computation. The program decomposition is guided by the source program annotations provided by the user. The synthesis of parallel code is based on configurations that describe the overall computation as a set of interacting components. Run-time support is provided by the compiler-generated code that redistributes computation and data during object program execution. The generated parallel code is optimized using techniques of data alignment, operator placement, wavefront determination, and memory optimization. In this article we discuss annotations, configurations, parallel code generation, and run-time support suitable for parallel programs written in the functional parallel programming language EPL and in Fortran.

RFQ modeling computer program

International Nuclear Information System (INIS)

Potter, J.M.

1985-01-01

The mathematical background for a multiport-network-solving program is described. A method for accurately numerically modeling an arbitrary, continuous, multiport transmission line is discussed. A modification to the transmission-line equations to accommodate multiple rf drives is presented. An improved model for the radio-frequency quadrupole (RFQ) accelerator that corrects previous errors is given. This model permits treating the RFQ as a true eight-port network for simplicity in interpreting the field distribution and ensures that all modes propagate at the same velocity in the high-frequency limit. The flexibility of the multiport model is illustrated by simple modifications to otherwise two-dimensional systems that permit modeling them as linear chains of multiport networks

Evaluating Equating Results: Percent Relative Error for Chained Kernel Equating

Science.gov (United States)

Jiang, Yanlin; von Davier, Alina A.; Chen, Haiwen

2012-01-01

This article presents a method for evaluating equating results. Within the kernel equating framework, the percent relative error (PRE) for chained equipercentile equating was computed under the nonequivalent groups with anchor test (NEAT) design. The method was applied to two data sets to obtain the PRE, which can be used to measure equating…

The Caltech Concurrent Computation Program - Project description

Science.gov (United States)

Fox, G.; Otto, S.; Lyzenga, G.; Rogstad, D.

1985-01-01

The Caltech Concurrent Computation Program wwhich studies basic issues in computational science is described. The research builds on initial work where novel concurrent hardware, the necessary systems software to use it and twenty significant scientific implementations running on the initial 32, 64, and 128 node hypercube machines have been constructed. A major goal of the program will be to extend this work into new disciplines and more complex algorithms including general packages that decompose arbitrary problems in major application areas. New high-performance concurrent processors with up to 1024-nodes, over a gigabyte of memory and multigigaflop performance are being constructed. The implementations cover a wide range of problems in areas such as high energy and astrophysics, condensed matter, chemical reactions, plasma physics, applied mathematics, geophysics, simulation, CAD for VLSI, graphics and image processing. The products of the research program include the concurrent algorithms, hardware, systems software, and complete program implementations.

Computational Experience with Globally Convergent Descent Methods for Large Sparse Systems of Nonlinear Equations

Czech Academy of Sciences Publication Activity Database

Lukšan, Ladislav; Vlček, Jan

1998-01-01

Roč. 8, č. 3-4 (1998), s. 201-223 ISSN 1055-6788 R&D Projects: GA ČR GA201/96/0918 Keywords : nonlinear equations * Armijo-type descent methods * Newton-like methods * truncated methods * global convergence * nonsymmetric linear systems * conjugate gradient -type methods * residual smoothing * computational experiments Subject RIV: BB - Applied Statistics, Operational Research

Heat exchanger performance analysis programs for the personal computer

International Nuclear Information System (INIS)

Putman, R.E.

1992-01-01

Numerous utility industry heat exchange calculations are repetitive and thus lend themselves to being performed on a Personal Computer. These programs may be regarded as engineering tools which, when put together, can form a Toolbox. However, the practicing Results Engineer in the utility industry desires not only programs that are robust as well as easy to use but can also be used both on desktop and laptop PC's. The latter also offer the opportunity to take the computer into the plant or control room, and use it there to process test or operating data right on the spot. Most programs evolve through the needs which arise in the course of day-to-day work. This paper describes several of the more useful programs of this type and outlines some of the guidelines to be followed when designing personal computer programs for use by the practicing Results Engineer

GAUSS VII: a computer program for the analysis of #betta#-ray spectra from Ge semiconductor spectrometers

International Nuclear Information System (INIS)

McCullagh, C.M.; Helmer, R.G.

1982-10-01

A description is given of a computer program, GAUSS VII, which has been written to determine #betta#-ray spectra from Ge semiconductor spectrometers. The preliminary portions of the program can determine the energy and width calibration equations, loacte individual peaks and define peak regions that are significantly above the local spectral background. The user may edit these lists of peaks and regions. Each peak region is fitted with one or more components in which the peaks are represented by a Gaussian function or a Gaussian with one or two additive exponential tails on the low-energy side and one on the high-energy side. A step-like background function can be used with each component. The program will automatically recycle to add one or more components to a region if needed to improve the fit. The #betta#-ray energies and intensities are computed from resulting Gaussian positions and peak areas. To allow the user to determine the best results, the results from the analyses for each region with different numbers of components can be printed and line-printer plots of the fits to the data can be made. The quality of the results depends primarily on the ability of the program to define a good spectral region for each analysis and the ability to recycle to determine the proper number of components

Computer-based Programs in Speech Therapy of Dyslalia and Dyslexia- Dysgraphia

Directory of Open Access Journals (Sweden)

Mirela Danubianu

2010-04-01

Full Text Available During the last years, the researchers and therapists in speech therapy have been more and more concerned with the elaboration and use of computer programs in speech disorders therapy. The main objective of this study was to evaluate the therapeutic effectiveness of computer-based programs for the Romanian language in speech therapy. Along the study, we will present the experimental research through assessing the effectiveness of computer programs in the speech therapy for speech disorders: dyslalia, dyslexia and dysgraphia. Methodologically, the use of the computer in the therapeutic phases was carried out with the help of some computer-based programs (Logomon, Dislex-Test etc. that we elaborated and we experimented with during several years of therapeutic activity. The sample used in our experiments was composed of 120 subjects; two groups of 60 children with speech disorders were selected for both speech disorders: 30 for the experimental ('computer-based' group and 30 for the control ('classical method' group. The study hypotheses verified whether the results, obtained by the subjects within the experimental group, improved significantly after using the computer-based program, compared to the subjects within the control group, who did not use this program but got a classical therapy. The hypotheses were confirmed for the speech disorders included in this research; the conclusions of the study confirm the advantages of using computer-based programs within speech therapy by correcting these disorders, as well as due to the positive influence these programs have on the development of children’s personality.

Computing Programs for Determining Traffic Flows from Roundabouts

Science.gov (United States)

Boroiu, A. A.; Tabacu, I.; Ene, A.; Neagu, E.; Boroiu, A.

2017-10-01

For modelling road traffic at the level of a road network it is necessary to specify the flows of all traffic currents at each intersection. These data can be obtained by direct measurements at the traffic light intersections, but in the case of a roundabout this is not possible directly and the literature as well as the traffic modelling software doesn’t offer ways to solve this issue. Two sets of formulas are proposed by which all traffic flows from the roundabouts with 3 or 4 arms are calculated based on the streams that can be measured. The objective of this paper is to develop computational programs to operate with these formulas. For each of the two sets of analytical relations, a computational program was developed in the Java operating language. The obtained results fully confirm the applicability of the calculation programs. The final stage for capitalizing these programs will be to make them web pages in HTML format, so that they can be accessed and used on the Internet. The achievements presented in this paper are an important step to provide a necessary tool for traffic modelling because these computational programs can be easily integrated into specialized software.

A linear programming manual

Science.gov (United States)

Tuey, R. C.

1972-01-01

Computer solutions of linear programming problems are outlined. Information covers vector spaces, convex sets, and matrix algebra elements for solving simultaneous linear equations. Dual problems, reduced cost analysis, ranges, and error analysis are illustrated.

OpenMx: An Open Source Extended Structural Equation Modeling Framework

Science.gov (United States)

Boker, Steven; Neale, Michael; Maes, Hermine; Wilde, Michael; Spiegel, Michael; Brick, Timothy; Spies, Jeffrey; Estabrook, Ryne; Kenny, Sarah; Bates, Timothy; Mehta, Paras; Fox, John

2011-01-01

OpenMx is free, full-featured, open source, structural equation modeling (SEM) software. OpenMx runs within the "R" statistical programming environment on Windows, Mac OS-X, and Linux computers. The rationale for developing OpenMx is discussed along with the philosophy behind the user interface. The OpenMx data structures are…

Review of Kaganove's solution for the reactor point kinetics equations

International Nuclear Information System (INIS)

Couto, R.T.; Santo, A.C.F. de.

1993-09-01

A review of Kaganove's method for the reactor point kinetics equations solution is performed. This was method chosen to calculate the power in ATR, a computer program for the analysis of reactivity transients. The reasons for this choice and the adaptation of the method to the purposes of ATR are presented. (author)

Computer program for computing the properties of seventeen fluids. [cryogenic liquids

Science.gov (United States)

Brennan, J. A.; Friend, D. G.; Arp, V. D.; Mccarty, R. D.

1992-01-01

The present study describes modifications and additions to the MIPROPS computer program for calculating the thermophysical properties of 17 fluids. These changes include adding new fluids, new properties, and a new interface to the program. The new program allows the user to select the input and output parameters and the units to be displayed for each parameter. Fluids added to the MIPROPS program are carbon dioxide, carbon monoxide, deuterium, helium, normal hydrogen, and xenon. The most recent modifications to the MIPROPS program are the addition of viscosity and thermal conductivity correlations for parahydrogen and the addition of the fluids normal hydrogen and xenon. The recently added interface considerably increases the program's utility.

Automatic control of commercial computer programs

International Nuclear Information System (INIS)

Rezvov, B.A.; Artem'ev, A.N.; Maevskij, A.G.; Demkiv, A.A.; Kirillov, B.F.; Belyaev, A.D.; Artem'ev, N.A.

2010-01-01

The way of automatic control of commercial computer programs is presented. The developed connection of the EXAFS spectrometer automatic system (which is managed by PC for DOS) is taken with the commercial program for the CCD detector control (which is managed by PC for Windows). The described complex system is used for the automation of intermediate amplitude spectra processing in EXAFS spectrum measurements at Kurchatov SR source

On a computer implementation of the block Gauss–Seidel method for normal systems of equations

Directory of Open Access Journals (Sweden)

Alexander I. Zhdanov

2016-12-01

Full Text Available This article focuses on the modification of the block option Gauss-Seidel method for normal systems of equations, which is a sufficiently effective method of solving generally overdetermined, systems of linear algebraic equations of high dimensionality. The main disadvantage of methods based on normal equations systems is the fact that the condition number of the normal system is equal to the square of the condition number of the original problem. This fact has a negative impact on the rate of convergence of iterative methods based on normal equations systems. To increase the speed of convergence of iterative methods based on normal equations systems, for solving ill-conditioned problems currently different preconditioners options are used that reduce the condition number of the original system of equations. However, universal preconditioner for all applications does not exist. One of the effective approaches that improve the speed of convergence of the iterative Gauss–Seidel method for normal systems of equations, is to use its version of the block. The disadvantage of the block Gauss–Seidel method for production systems is the fact that it is necessary to calculate the pseudoinverse matrix for each iteration. We know that finding the pseudoinverse is a difficult computational procedure. In this paper, we propose a procedure to replace the matrix pseudo-solutions to the problem of normal systems of equations by Cholesky. Normal equations arising at each iteration of Gauss–Seidel method, have a relatively low dimension compared to the original system. The results of numerical experimentation demonstrating the effectiveness of the proposed approach are given.

The computation of equating errors in international surveys in education.

Science.gov (United States)

Monseur, Christian; Berezner, Alla

2007-01-01

Since the IEA's Third International Mathematics and Science Study, one of the major objectives of international surveys in education has been to report trends in achievement. The names of the two current IEA surveys reflect this growing interest: Trends in International Mathematics and Science Study (TIMSS) and Progress in International Reading Literacy Study (PIRLS). Similarly a central concern of the OECD's PISA is with trends in outcomes over time. To facilitate trend analyses these studies link their tests using common item equating in conjunction with item response modelling methods. IEA and PISA policies differ in terms of reporting the error associated with trends. In IEA surveys, the standard errors of the trend estimates do not include the uncertainty associated with the linking step while PISA does include a linking error component in the standard errors of trend estimates. In other words, PISA implicitly acknowledges that trend estimates partly depend on the selected common items, while the IEA's surveys do not recognise this source of error. Failing to recognise the linking error leads to an underestimation of the standard errors and thus increases the Type I error rate, thereby resulting in reporting of significant changes in achievement when in fact these are not significant. The growing interest of policy makers in trend indicators and the impact of the evaluation of educational reforms appear to be incompatible with such underestimation. However, the procedure implemented by PISA raises a few issues about the underlying assumptions for the computation of the equating error. After a brief introduction, this paper will describe the procedure PISA implemented to compute the linking error. The underlying assumptions of this procedure will then be discussed. Finally an alternative method based on replication techniques will be presented, based on a simulation study and then applied to the PISA 2000 data.

Comparison of preconditioned generalized conjugate gradient methods to two-dimensional neutron and photon transport equation

International Nuclear Information System (INIS)

Chen, G.S.

1997-01-01

We apply and compare the preconditioned generalized conjugate gradient methods to solve the linear system equation that arises in the two-dimensional neutron and photon transport equation in this paper. Several subroutines are developed on the basis of preconditioned generalized conjugate gradient methods for time-independent, two-dimensional neutron and photon transport equation in the transport theory. These generalized conjugate gradient methods are used. TFQMR (transpose free quasi-minimal residual algorithm), CGS (conjuage gradient square algorithm), Bi-CGSTAB (bi-conjugate gradient stabilized algorithm) and QMRCGSTAB (quasi-minimal residual variant of bi-conjugate gradient stabilized algorithm). These sub-routines are connected to computer program DORT. Several problems are tested on a personal computer with Intel Pentium CPU. (author)

Radiation and pregnancy: a self-teaching computer program

International Nuclear Information System (INIS)

Kumar, Pratik; Rehani, M.M.

1994-01-01

Two self-interacting computer programs have been developed. The first program which consists of fifteen topics apprises the users with broad spectrum of radiation risks to the unborn during pregnancy and the status of various views in this regard. Another program estimates the dose to uterus in sixteen radiological examinations depending upon the radiographic parameters used. The dose to uterus and hence to the fetus calculated by computer program in different radiographic views have been found to be in agreement with that reported in NRPB-R200 survey report. The two programs combined provide a better understanding of the rather confusing situation regarding dilemma about termination of pregnancy following inadvertent radiation exposure, apprehension about radiation effect in the minds of prescribing doctor and patients, dose estimation and advice to pregnant workers and like. (author). 10 refs

Use of computer programs STLK1 and STWT1 for analysis of stream-aquifer hydraulic interaction

Science.gov (United States)

Desimone, Leslie A.; Barlow, Paul M.

1999-01-01

Quantifying the hydraulic interaction of aquifers and streams is important in the analysis of stream base fow, flood-wave effects, and contaminant transport between surface- and ground-water systems. This report describes the use of two computer programs, STLK1 and STWT1, to analyze the hydraulic interaction of streams with confined, leaky, and water-table aquifers during periods of stream-stage fuctuations and uniform, areal recharge. The computer programs are based on analytical solutions to the ground-water-flow equation in stream-aquifer settings and calculate ground-water levels, seepage rates across the stream-aquifer boundary, and bank storage that result from arbitrarily varying stream stage or recharge. Analysis of idealized, hypothetical stream-aquifer systems is used to show how aquifer type, aquifer boundaries, and aquifer and streambank hydraulic properties affect aquifer response to stresses. Published data from alluvial and stratifed-drift aquifers in Kentucky, Massachusetts, and Iowa are used to demonstrate application of the programs to field settings. Analytical models of these three stream-aquifer systems are developed on the basis of available hydrogeologic information. Stream-stage fluctuations and recharge are applied to the systems as hydraulic stresses. The models are calibrated by matching ground-water levels calculated with computer program STLK1 or STWT1 to measured ground-water levels. The analytical models are used to estimate hydraulic properties of the aquifer, aquitard, and streambank; to evaluate hydrologic conditions in the aquifer; and to estimate seepage rates and bank-storage volumes resulting from flood waves and recharge. Analysis of field examples demonstrates the accuracy and limitations of the analytical solutions and programs when applied to actual ground-water systems and the potential uses of the analytical methods as alternatives to numerical modeling for quantifying stream-aquifer interactions.

Web Program for Development of GUIs for Cluster Computers

Science.gov (United States)

Czikmantory, Akos; Cwik, Thomas; Klimeck, Gerhard; Hua, Hook; Oyafuso, Fabiano; Vinyard, Edward

2003-01-01

WIGLAF (a Web Interface Generator and Legacy Application Facade) is a computer program that provides a Web-based, distributed, graphical-user-interface (GUI) framework that can be adapted to any of a broad range of application programs, written in any programming language, that are executed remotely on any cluster computer system. WIGLAF enables the rapid development of a GUI for controlling and monitoring a specific application program running on the cluster and for transferring data to and from the application program. The only prerequisite for the execution of WIGLAF is a Web-browser program on a user's personal computer connected with the cluster via the Internet. WIGLAF has a client/server architecture: The server component is executed on the cluster system, where it controls the application program and serves data to the client component. The client component is an applet that runs in the Web browser. WIGLAF utilizes the Extensible Markup Language to hold all data associated with the application software, Java to enable platform-independent execution on the cluster system and the display of a GUI generator through the browser, and the Java Remote Method Invocation software package to provide simple, effective client/server networking.

Chemical Equation Balancing.

Science.gov (United States)

Blakley, G. R.

1982-01-01

Reviews mathematical techniques for solving systems of homogeneous linear equations and demonstrates that the algebraic method of balancing chemical equations is a matter of solving a system of homogeneous linear equations. FORTRAN programs using this matrix method to chemical equation balancing are available from the author. (JN)

Integrating computer programs for engineering analysis and design

Science.gov (United States)

Wilhite, A. W.; Crisp, V. K.; Johnson, S. C.

1983-01-01

The design of a third-generation system for integrating computer programs for engineering and design has been developed for the Aerospace Vehicle Interactive Design (AVID) system. This system consists of an engineering data management system, program interface software, a user interface, and a geometry system. A relational information system (ARIS) was developed specifically for the computer-aided engineering system. It is used for a repository of design data that are communicated between analysis programs, for a dictionary that describes these design data, for a directory that describes the analysis programs, and for other system functions. A method is described for interfacing independent analysis programs into a loosely-coupled design system. This method emphasizes an interactive extension of analysis techniques and manipulation of design data. Also, integrity mechanisms exist to maintain database correctness for multidisciplinary design tasks by an individual or a team of specialists. Finally, a prototype user interface program has been developed to aid in system utilization.

PARALLEL SOLUTION METHODS OF PARTIAL DIFFERENTIAL EQUATIONS

Directory of Open Access Journals (Sweden)

Korhan KARABULUT

1998-03-01

Full Text Available Partial differential equations arise in almost all fields of science and engineering. Computer time spent in solving partial differential equations is much more than that of in any other problem class. For this reason, partial differential equations are suitable to be solved on parallel computers that offer great computation power. In this study, parallel solution to partial differential equations with Jacobi, Gauss-Siedel, SOR (Succesive OverRelaxation and SSOR (Symmetric SOR algorithms is studied.

Numerical Solution of Heun Equation Via Linear Stochastic Differential Equation

Directory of Open Access Journals (Sweden)

Hamidreza Rezazadeh

2014-05-01

Full Text Available In this paper, we intend to solve special kind of ordinary differential equations which is called Heun equations, by converting to a corresponding stochastic differential equation(S.D.E.. So, we construct a stochastic linear equation system from this equation which its solution is based on computing fundamental matrix of this system and then, this S.D.E. is solved by numerically methods. Moreover, its asymptotic stability and statistical concepts like expectation and variance of solutions are discussed. Finally, the attained solutions of these S.D.E.s compared with exact solution of corresponding differential equations.

A Tangible Programming Tool for Children to Cultivate Computational Thinking

Science.gov (United States)

Wang, Danli; Liu, Zhen

2014-01-01

Game and creation are activities which have good potential for computational thinking skills. In this paper we present T-Maze, an economical tangible programming tool for children aged 5–9 to build computer programs in maze games by placing wooden blocks. Through the use of computer vision technology, T-Maze provides a live programming interface with real-time graphical and voice feedback. We conducted a user study with 7 children using T-Maze to play two levels of maze-escape games and create their own mazes. The results show that T-Maze is not only easy to use, but also has the potential to help children cultivate computational thinking like abstraction, problem decomposition, and creativity. PMID:24719575

A Tangible Programming Tool for Children to Cultivate Computational Thinking

Directory of Open Access Journals (Sweden)

Danli Wang

2014-01-01

Full Text Available Game and creation are activities which have good potential for computational thinking skills. In this paper we present T-Maze, an economical tangible programming tool for children aged 5–9 to build computer programs in maze games by placing wooden blocks. Through the use of computer vision technology, T-Maze provides a live programming interface with real-time graphical and voice feedback. We conducted a user study with 7 children using T-Maze to play two levels of maze-escape games and create their own mazes. The results show that T-Maze is not only easy to use, but also has the potential to help children cultivate computational thinking like abstraction, problem decomposition, and creativity.

Computation of a numerically satisfactory pair of solutions of the differential equation for conical functions of non-negative integer orders

NARCIS (Netherlands)

T.M. Dunster (Mark); A. Gil (Amparo); J. Segura (Javier); N.M. Temme (Nico)

2014-01-01

textabstractWe consider the problem of computing satisfactory pair of solutions of the differential equation for Legendre functions of non-negative integer order $\\mu$ and degree $-\\frac12+i\\tau$, where $\\tau$ is a non-negative real parameter. Solutions of this equation are the conical functions

The AAHA Computer Program. American Animal Hospital Association.

Science.gov (United States)

Albers, J W

1986-07-01

The American Animal Hospital Association Computer Program should benefit all small animal practitioners. Through the availability of well-researched and well-developed certified software, veterinarians will have increased confidence in their purchase decisions. With the expansion of computer applications to improve practice management efficiency, veterinary computer systems will further justify their initial expense. The development of the Association's veterinary computer network will provide a variety of important services to the profession.

Using the Computer in Special Vocational Programs. Inservice Activities.

Science.gov (United States)

Lane, Kenneth; Ward, Raymond

This inservice manual is intended to assist vocational education teachers in using the techniques of computer-assisted instruction in special vocational education programs. Addressed in the individual units are the following topics: the basic principles of computer-assisted instruction (TRS-80 computers and typing on a computer keyboard); money…

Programs for Testing Processor-in-Memory Computing Systems

Science.gov (United States)

Katz, Daniel S.

2006-01-01

The Multithreaded Microbenchmarks for Processor-In-Memory (PIM) Compilers, Simulators, and Hardware are computer programs arranged in a series for use in testing the performances of PIM computing systems, including compilers, simulators, and hardware. The programs at the beginning of the series test basic functionality; the programs at subsequent positions in the series test increasingly complex functionality. The programs are intended to be used while designing a PIM system, and can be used to verify that compilers, simulators, and hardware work correctly. The programs can also be used to enable designers of these system components to examine tradeoffs in implementation. Finally, these programs can be run on non-PIM hardware (either single-threaded or multithreaded) using the POSIX pthreads standard to verify that the benchmarks themselves operate correctly. [POSIX (Portable Operating System Interface for UNIX) is a set of standards that define how programs and operating systems interact with each other. pthreads is a library of pre-emptive thread routines that comply with one of the POSIX standards.

Computer program user's manual for FIREFINDER digital topographic data verification library dubbing system

Science.gov (United States)

Ceres, M.; Heselton, L. R., III

1981-11-01

This manual describes the computer programs for the FIREFINDER Digital Topographic Data Verification-Library-Dubbing System (FFDTDVLDS), and will assist in the maintenance of these programs. The manual contains detailed flow diagrams and associated descriptions for each computer program routine and subroutine. Complete computer program listings are also included. This information should be used when changes are made in the computer programs. The operating system has been designed to minimize operator intervention.

Scientific Computing in the CH Programming Language

Directory of Open Access Journals (Sweden)

Harry H. Cheng

1993-01-01

Full Text Available We have developed a general-purpose block-structured interpretive programming Ianguage. The syntax and semantics of this language called CH are similar to C. CH retains most features of C from the scientific computing point of view. In this paper, the extension of C to CH for numerical computation of real numbers will be described. Metanumbers of −0.0, 0.0, Inf, −Inf, and NaN are introduced in CH. Through these metanumbers, the power of the IEEE 754 arithmetic standard is easily available to the programmer. These metanumbers are extended to commonly used mathematical functions in the spirit of the IEEE 754 standard and ANSI C. The definitions for manipulation of these metanumbers in I/O; arithmetic, relational, and logic operations; and built-in polymorphic mathematical functions are defined. The capabilities of bitwise, assignment, address and indirection, increment and decrement, as well as type conversion operations in ANSI C are extended in CH. In this paper, mainly new linguistic features of CH in comparison to C will be described. Example programs programmed in CH with metanumbers and polymorphic mathematical functions will demonstrate capabilities of CH in scientific computing.

Energy consumption program: A computer model simulating energy loads in buildings

Science.gov (United States)

Stoller, F. W.; Lansing, F. L.; Chai, V. W.; Higgins, S.

1978-01-01

The JPL energy consumption computer program developed as a useful tool in the on-going building modification studies in the DSN energy conservation project is described. The program simulates building heating and cooling loads and computes thermal and electric energy consumption and cost. The accuracy of computations are not sacrificed, however, since the results lie within + or - 10 percent margin compared to those read from energy meters. The program is carefully structured to reduce both user's time and running cost by asking minimum information from the user and reducing many internal time-consuming computational loops. Many unique features were added to handle two-level electronics control rooms not found in any other program.

Analytical calculation of an invariant curve for Chew-Low equations

International Nuclear Information System (INIS)

Gerdt, V.P.

1978-01-01

The local structure of the one-parameter set of invariant curves for Chew-Low equations having the form of the convergent series is considered. Coefficients of this series βsub(i)(C) are polynomials in set parameter C. The transition to the general solution of Chew-Low equations is carried out by replacing the parameter C by arbitrary even, real, meromorphic function C(w) with the property C(w+1)=-C(w). The procedure for calculation of coefficients βsub(i)(C), which is based on the solution of nonlinear functional eqtions, following from Chew-Low equations, is developed. First twelve coefficients βsub(i)(C) are calculated analytically by computer, using program system SCHOONSCHIP

Computer program for calculating thermodynamic and transport properties of fluids

Science.gov (United States)

Hendricks, R. C.; Braon, A. K.; Peller, I. C.

1975-01-01

Computer code has been developed to provide thermodynamic and transport properties of liquid argon, carbon dioxide, carbon monoxide, fluorine, helium, methane, neon, nitrogen, oxygen, and parahydrogen. Equation of state and transport coefficients are updated and other fluids added as new material becomes available.

Basic design of parallel computational program for probabilistic structural analysis

International Nuclear Information System (INIS)

Kaji, Yoshiyuki; Arai, Taketoshi; Gu, Wenwei; Nakamura, Hitoshi

1999-06-01

In our laboratory, for 'development of damage evaluation method of structural brittle materials by microscopic fracture mechanics and probabilistic theory' (nuclear computational science cross-over research) we examine computational method related to super parallel computation system which is coupled with material strength theory based on microscopic fracture mechanics for latent cracks and continuum structural model to develop new structural reliability evaluation methods for ceramic structures. This technical report is the review results regarding probabilistic structural mechanics theory, basic terms of formula and program methods of parallel computation which are related to principal terms in basic design of computational mechanics program. (author)

Basic design of parallel computational program for probabilistic structural analysis

Energy Technology Data Exchange (ETDEWEB)

Kaji, Yoshiyuki; Arai, Taketoshi [Japan Atomic Energy Research Inst., Tokai, Ibaraki (Japan). Tokai Research Establishment; Gu, Wenwei; Nakamura, Hitoshi

1999-06-01

In our laboratory, for `development of damage evaluation method of structural brittle materials by microscopic fracture mechanics and probabilistic theory` (nuclear computational science cross-over research) we examine computational method related to super parallel computation system which is coupled with material strength theory based on microscopic fracture mechanics for latent cracks and continuum structural model to develop new structural reliability evaluation methods for ceramic structures. This technical report is the review results regarding probabilistic structural mechanics theory, basic terms of formula and program methods of parallel computation which are related to principal terms in basic design of computational mechanics program. (author)

Application research of computational mass-transfer differential equation in MBR concentration field simulation.

Science.gov (United States)

Li, Chunqing; Tie, Xiaobo; Liang, Kai; Ji, Chanjuan

2016-01-01

After conducting the intensive research on the distribution of fluid's velocity and biochemical reactions in the membrane bioreactor (MBR), this paper introduces the use of the mass-transfer differential equation to simulate the distribution of the chemical oxygen demand (COD) concentration in MBR membrane pool. The solutions are as follows: first, use computational fluid dynamics to establish a flow control equation model of the fluid in MBR membrane pool; second, calculate this model by adopting direct numerical simulation to get the velocity field of the fluid in membrane pool; third, combine the data of velocity field to establish mass-transfer differential equation model for the concentration field in MBR membrane pool, and use Seidel iteration method to solve the equation model; last but not least, substitute the real factory data into the velocity and concentration field model to calculate simulation results, and use visualization software Tecplot to display the results. Finally by analyzing the nephogram of COD concentration distribution, it can be found that the simulation result conforms the distribution rule of the COD's concentration in real membrane pool, and the mass-transfer phenomenon can be affected by the velocity field of the fluid in membrane pool. The simulation results of this paper have certain reference value for the design optimization of the real MBR system.

JAC2D: A two-dimensional finite element computer program for the nonlinear quasi-static response of solids with the conjugate gradient method

International Nuclear Information System (INIS)

Biffle, J.H.; Blanford, M.L.

1994-05-01

JAC2D is a two-dimensional finite element program designed to solve quasi-static nonlinear mechanics problems. A set of continuum equations describes the nonlinear mechanics involving large rotation and strain. A nonlinear conjugate gradient method is used to solve the equations. The method is implemented in a two-dimensional setting with various methods for accelerating convergence. Sliding interface logic is also implemented. A four-node Lagrangian uniform strain element is used with hourglass stiffness to control the zero-energy modes. This report documents the elastic and isothermal elastic/plastic material model. Other material models, documented elsewhere, are also available. The program is vectorized for efficient performance on Cray computers. Sample problems described are the bending of a thin beam, the rotation of a unit cube, and the pressurization and thermal loading of a hollow sphere

JAC3D -- A three-dimensional finite element computer program for the nonlinear quasi-static response of solids with the conjugate gradient method

International Nuclear Information System (INIS)

Biffle, J.H.

1993-02-01

JAC3D is a three-dimensional finite element program designed to solve quasi-static nonlinear mechanics problems. A set of continuum equations describes the nonlinear mechanics involving large rotation and strain. A nonlinear conjugate gradient method is used to solve the equation. The method is implemented in a three-dimensional setting with various methods for accelerating convergence. Sliding interface logic is also implemented. An eight-node Lagrangian uniform strain element is used with hourglass stiffness to control the zero-energy modes. This report documents the elastic and isothermal elastic-plastic material model. Other material models, documented elsewhere, are also available. The program is vectorized for efficient performance on Cray computers. Sample problems described are the bending of a thin beam, the rotation of a unit cube, and the pressurization and thermal loading of a hollow sphere

Overview of the history of neutronics computations

International Nuclear Information System (INIS)

Gelbard, E.M.

1992-01-01

This paper provides some background and perspective on the history of neutronics computations. The author started his carerr at Bettis Atomic Power Laboratory (BAPL) in 1954. At that time, computers were primitive and unfamiliar to most of us. Bettis did own an IBM Card Programmed Computer, which, as its name suggested, had a punched-card memory, and some of the industry's first cross-section preparation codes (the early MUFTs) were written for this machine by Ben Mount. There was also at BAPL another digital machine, which spewed output everywhere on punched tape. But the laboratory's most powerful machine was an analog computer that solved two-group diffusion equations in two dimensions. An excellent physicist of that period argued that analog machines were the wave of the future because digital computers could never solve the diffusion equation in two dimensions. At the time, this seemed like a very reasonable prediction; yet it was not very long before Bettis was renting time on UNIVACs all over the country precisely to solve the two-dimensional diffusion equation. Since then, we have learned that it is always risky to predict future trends in the computer industry

The integral equation method applied to eddy currents

International Nuclear Information System (INIS)

Biddlecombe, C.S.; Collie, C.J.; Simkin, J.; Trowbridge, C.W.

1976-04-01

An algorithm for the numerical solution of eddy current problems is described, based on the direct solution of the integral equation for the potentials. In this method only the conducting and iron regions need to be divided into elements, and there are no boundary conditions. Results from two computer programs using this method for iron free problems for various two-dimensional geometries are presented and compared with analytic solutions. (author)

Introduction to programming multiple-processor computers

International Nuclear Information System (INIS)

Hicks, H.R.; Lynch, V.E.

1985-04-01

FORTRAN applications programs can be executed on multiprocessor computers in either a unitasking (traditional) or multitasking form. The latter allows a single job to use more than one processor simultaneously, with a consequent reduction in wall-clock time and, perhaps, the cost of the calculation. An introduction to programming in this environment is presented. The concepts of synchronization and data sharing using EVENTS and LOCKS are illustrated with examples. The strategy of strong synchronization and the use of synchronization templates are proposed. We emphasize that incorrect multitasking programs can produce irreproducible results, which makes debugging more difficult

NEWBOX: A computer program for parameter estimation in diffusion problems

International Nuclear Information System (INIS)

Nestor, C.W. Jr.; Godbee, H.W.; Joy, D.S.

1989-01-01