Energy Technology Data Exchange (ETDEWEB)

Piping in nuclear power plants is vital to the proper operation and safety of these facilities. To assure safety in the unlikely event of a pipe break, it is necessary to evaluate the consequences from the resulting whipping pipe on neighboring components and structures. Numerical simulations allow for rapid evaluation of these consequences. Before simulations can be accepted, however, the methodology and computer codes must be validated against experimental results. This paper uses a probabilistic approach to validate pipe whip simulations against limited experimental results. Probabilistic analysis software was developed and coupled to existing deterministic finite element software. An example of a whipping pipe impacting against a reinforced concrete slab was simulated. The described probabilistic approach was used to ...

International Nuclear Information System (INIS)

In this paper, numerical simulation of a concentric tube heat exchanger is presented to determine the convective heat transfer coefficient and friction factor in a smooth tube. Increasing the convective heat transfer coefficient can increase heat transfer rate in a concentric tube heat exchanger from a given tubular surface area. This can be achieved by using heat transfer augmentation devices. This work constitutes the initial phase of the numerical simulation of heat transfer from tubes employing augmentation devices, such as twisted tapes, wire-coil inserts, for heat transfer enhancement. A computational fluid dynamics (CFD) simulation tool was developed with CFX software and the results obtained from the simulations are validated with the empirical correlations for a smooth tube heat exchanger. The difficulties associated with the ...

International Nuclear Information System (INIS)

The validation of the forecast model for early emergency response to nuclear accidents is evaluated by trace tests in atmosphere in Daya Bay nuclear power site. The simulation experiment of the Daya Bay nuclear power site shows that the particle spreading image and the time-integrated concentration distribution given by plume concentration prediction model can perform the variation of pathway of the pollutant transport, as well as the effects of topography on transport and diffusion of pollutants. The simulation of five trace tests in field shows that 59.1% of ratios between predicted results and observed results are within the range of 10, and 41% of ratios are within the range of 5 approximately. (authors)

Energy Technology Data Exchange (ETDEWEB)

This document reproduces the final project of David San Fabian Ayuso, presented on May 26, 2005, for the obtention of the engineer degree of the Carlos III University of Madrid. A single-phase, isothermal model, including both electron and proton transport, is introduced for the simulation of polymer electrolyte fuel cells (PEM). The model is implemented in the commercial code Fluent 6.0, through the use of UDFs (User Defined Functions). In order to validate the model, a single canal of a PEM monocell is simulated in three dimensions. The obtained result are qualitatively satisfactory. It is observed that it is not essential to solve the current collectors when a monocell is considered (and not a stack). in the present study, the number of nodes is the computational grid appears to be too low in the membrane zone in order to make a complete validation of the model. (Author) 20 refs.

British Library Electronic Table of Contents (United Kingdom)

SummaryA meshless numerical model is proposed to investigate shallow-water dam break flows in 1D open channels. The numerical model is to solve the shallow water equations (SWE) based on smoothed particle hydrodynamics (SPH). The concept of slice water particles (SWP) is adopted in the SPH-SWE formulation. The numerical sensitivity analysis is first performed to study the appropriate SWP number and variable smoothing length through dam break flows in an idealized 1D channel with dry/wet beds. Extensive validation by comparison with laboratory and field data is next conducted for four benchmark problems, including dam break flows through a rough flat channel, a rough bumpy channel with various downstream boundary conditions, a nonprismatic channel, and a realistic scale model of the Toce ri...

CERN Document Server

We examine the effect of a threshold bias on the power spectrum and the bispectrum in an ensemble of numerical simulations (Gaussian initial perturbations with power law spectra P(k) \\sim k^n, n=+1, 0, -1, -2) and compare our results with theoretical predictions. Our simulations are evolved sufficiently that on the scale where we apply the threshold the rms fluctuation has developed significantly into the nonlinear regime. Thus, predictions based on perturbation theory do not necessarily apply. Nevertheless, we find our results for the power spectrum, biased power simply amplified by a numerical factor, follow predicted trends, far beyond the regime where perturbation theory is expected to be valid. We find that the biased bispectrum continues to follow the so-called hierarchical form, with reduced three-point amplitude Q \\approx 1 in the strongly nonlinear regime, independent of ...

Energy Technology Data Exchange (ETDEWEB)

Analytical and numerical calculations are presented for a reflexing electron beam type of collective ion accelerator. These results are then compared to those obtained through experiment. By constraining one free parameter to experimental conditions, the self-similar solution of the ion energy distribution agrees closely with the experimental distribution. Hence the reflexing beam model appears to be a valid model for explaining the experimental data. Simulation shows in addition to the agreement with the experimental ion distribution that synchronization between accelerated ions and electric field is phase unstable. This instability seems to further restrict the maximum ion energy to several times the electron energy.

Science.gov (United States)

This viewgraph presentation provides an overview of research designed to provide a numerical

Energy Technology Data Exchange (ETDEWEB)

Although incompressible fluid flows can be regarded as a particular case of a general problem, numerical methods and the mathematical formulation aimed to solve compressible and incompressible flows have their own peculiarities, in such a way, that it is generally not possible to attain both regimes with a single approach. In this work, we start from a typically compressible formulation, slightly modified to make use of pressure variables and, through augmenting the stabilising parameters, we end up with a simplified model which is able to deal with a wide range of flow regimes, from supersonic to low speed gas flows. The resulting methodology is flexible enough to allow for the simulation of liquid flows as well. Examples using conservative and pressure variables are shown and the results are compared to those published in the literature, in order to validate the method. (author)

British Library Electronic Table of Contents (United Kingdom)

The use of numerical modeling in oil spill incidents is a well established technique that has proven to provide cost-effective and reasonable estimates of oil surface drift. Good predictability of such models depends highly on the quality of the input data of the incident and on the model calibration effort. This paper presents the results of simulating oil spillage trajectory in the Arabian (Persian) Gulf. The study employed a 3-D rectilinear hydrodynamic model combined with oil spill model. Typical representative environmental conditions of the Arabian Gulf were first setup into a hydrodynamic circulation model using data from various sources. The performance of the hydrodynamic model was then tested against measurements of tidal fluctuation and sea currents at selected locations. The sp...

Energy Technology Data Exchange (ETDEWEB)

Squeeze film dampers (SFD) are devices utilized to control vibrations of the shafts of high-speed rotating machinery. The SFD - squirrel cage combination is probably the most used system for tuning the stiffness and damping of the supports for rotors installed on ball bearings. Squeeze film dampers are essentially hydrodynamic bearings which contain the ball bearings housings of ball-bearings supported shafts. Consequently, the oil film within the SFD are influenced only by the precession and nutation of the shaft, that is the flow of the oil within the damper is not directly influenced by the spin of the rotor. However, in the classical theory, the flow in the thin film is also governed by the Reynolds equation. In this paper, some of the limits of the classical theory of the SFD are discussed and theoretical and experimental studies, which illustrate the ideas presented herein, are presented as well. The orbits of an unbalanced rotor that is supported by a ...

Science.gov (United States)

Turbulent flow in a 2D channel with repeated rectangular rib roughness was numerically simulated using a low Reynolds number form of the k-epsilon turbulence model. Friction factors and average Stanton numbers were calculated for various pitch to rib height ratios and bulk Reynolds numbers. Comparisons with experiment were generally adequate, with the predictions of friction superior to those for heat transfer. The effect of variable properties for channel flow was investigated, and the results showed a greater effect for friction than for heat transfer. Comparison with experiment yielded no clear conclusions. The turbulence model was also validated for a related problem, that of flow downstream of an abrupt pipe expansion.

British Library Electronic Table of Contents (United Kingdom)

A constitutive model for an explosive consists of three parts: equation of state of reactants, equation of state of products, and reaction rate. In the absence of a good physical theory, models are empirical and based on heuristics. For the HMX based plastic-bonded explosive PBX 9501, three models have previously been calibrated and used in numerical simulations. Here, these three models are compared. They illustrate several generic difficulties with developing and validating models. The available data is insufficient to determine a constitutive model uniquely. Different models can be calibrated to fit the available data yet differ significantly in regions of phase space of importance to applications. Comparing models gives an indication of the uncertainty in model predictions outside the ...

Energy Technology Data Exchange (ETDEWEB)

A validated computational fluid dynamics (CFD) computer code, ICRKFLO, was used to investigate the scale-up effects on the coke yields of thermal cracking riser factors. Comparisons were made for calculated coke yields of pilot- and commercial-scales riser units. Computational results show that the riser aspect ratio, reaction temperature, particle residence time, and particle/oil ratio have major impacts on the coke yield. A computational experiment was conducted to determine optimal operating conditions for a conceptual design of a commercial-scale riser unit. This experiment showed that the performance loss in scale-up from pilot to commercial scale may be almost completely recovered through optimizing the operating conditions after scale-up using the CFD simulations as a guide.

International Nuclear Information System (INIS)

Radio frequency (RF) sheaths are suspected of limiting the performance of present-day ion cyclotron range of frequencies (ICRFs) antennas over long pulses and should be minimized in future fusion devices. Within the simplest models, RF-sheath effects are quantified by the integral VRF = ? E|| ? dl where the parallel RF field E|| is linked with the slow wave. On 'long open field lines' with large toroidal extension on both sides of the antenna it was shown that VRF is excited by parallel RF currents j|| flowing on the antenna structure. In this paper, the validity of this simple sheath theory is tested experimentally on the Tore Supra (TS) ITER-like antenna prototype (ILP), together with antenna simulation and post-processing codes developed to compute VRF. The predicted poloidal localization of high-|VRF| zones is confronted to that inferred from experimental data analysis. Surface temperature distribution on ILP front face, as well as ...

International Nuclear Information System (INIS)

The radiative transport problem is cast in integral form using a transport kernel. The transport kernel has an explicit representation in terms of a Feynman Path Integral over all paths between selected points in a volume. This representation is setup in detail. Numerical evaluation of this Path Integral is formulated with a Frenet-Serret based procedure for generating valid random paths, and with a numerical evaluation of the weight for each valid path. Very early sanity checks of a numerical implementation are reported. Approaches to optimization are identified. (authors)

Energy Technology Data Exchange (ETDEWEB)

2D simulation modeling of a river basin is based on notion of a River Basin as an opened non equilibrium nature system. All components of the system (elevation, water, soil properties and others) are linked to single multi - functional model. The objective of flood simulation is one the model functions. The objectives of the paper are as follows: a mean scale basin floods computing and the model validation; a large scale basin simulation and demands for its validation; results and demands for flood hazard assessment. (orig.)

Energy Technology Data Exchange (ETDEWEB)

A numerical investigation of a humidification dehumidification desalination (HDD) process using solar energy is presented. The HDD system consists mainly of a concentrating solar water heating collector, flat plate solar air heating collector, humidifying tower and dehumidifying exchanger. Two separate circulating loops constitute the HDD system, the first for heating the feed water and the second for heating air. A mathematical model is developed, simulating the HDD system, to study the influence of the different system configurations, weather and operating conditions on the system productivity. The model validity is examined by comparing the theoretical and experimental results of the same authors. It is found that the results of the developed mathematical model are in good agreement with the experimental results and other published works. The results show also that the productivity of the unit is strongly influenced by ...

Science.gov (United States)

The PLANCK SURVEYOR is a European Space Agency satellite mission to image the very faint anisotropies in the temperature of the Cosmic Microwave Background (CMB) radiation. Maynooth is actively participating in an international collaboration of scientists involved with the optical design of the High Frequency Instrument (HFI). This paper outlines research which has been undertaken in Maynooth concerned with numerical modelling of the optical characteristics of the multi-frequency array of detectors making up the HFI instrument. In the study the commercial software package ZEMAX was used to model the coupling of the focal plane HFI detectors to the PLANCK telescope. This package is particularly useful in the optical design of PLANCK because of the powerful optimisation features of the software. It is thus possible to readily determine the optimum positioning of the detectors in the focal plane of the telescope. Although the package is based on a ...

International Nuclear Information System (INIS)

As a three-year joint university-industry effort, development of a generalized boiling transition analysis method has been started in 2002 aiming at enhanced capabilities of subchannel analysis for a wide variety of BWR-type fuel bundle geometry from ordinary BWR to tight lattice fuel bundles. For this purpose, five dominant factors affecting boiling transition phenomena have been identified on which our efforts of experimentation and numerical analyses are focused. In this report, as the first-year achievement, we will describe a master plan of the development and contents for experimental approaches to construct thermal-hydraulic databases. The databases will be utilized for the developments of constitutive equations to describe the basic characteristics of the elementary processes. The planned experiments are divided into two groups. One is air-water experiments at atmospheric pressure, and the other is steam-water experiments up to 1 MPa. The former group of ...

Science.gov (United States)

The Predictive Capability Maturity Model (PCMM) is a new model that can be used to assess the level of maturity of computational modeling and simulation (M&S) efforts. The development of the model is based on both the authors experience and their analysis of similar investigations in the past. The perspective taken in this report is one of judging the usefulness of a predictive capability that relies on the numerical solution to partial differential equations to better inform and improve decision making. The review of past investigations, such as the Software Engineering Institute's Capability Maturity Model Integration and the National Aeronautics and Space Administration and Department of Defense Technology Readiness Levels, indicates that a more restricted, more interpretable method is needed to assess the maturity of an M&S effort. The PCMM addresses six contributing elements to M&S: (1) representation and geometric ...

British Library Electronic Table of Contents (United Kingdom)

A three-dimensional numerical simulation of mammatus-like clouds is presented. A portion of a cirrus outflow anvil cloud is simulated including cloud ice and snow microphysical representations. The simulated mammatus clouds appear in a cellular pattern and are compared with the few available previously published physical observations of mammatus. Copyright Copyright 2006 Royal Meteorological Society

International Nuclear Information System (INIS)

A physical model of free-electron laser (FEL) amplifier with variable-parameter wiggler magnets for one-dimensional numerical simulation is presented and a numerical example is given. The wiggler parameters, efficiency of energy conversion between electron beam and laser field, laser intensity, phase-space distributions and energy spectrum of electrons are computed. The period of synchronous oscillation and saturation value of laser intensity agree with estimated one.

Energy Technology Data Exchange (ETDEWEB)

Biodiesel made from vegetable oil is among the most desirable of renewable energy sources because it can be a substitute for diesel oil. However, biodiesel from soybean or corn can be confronted with a food crisis. Microalgae is a new biodiesel source which contains high oil lipids with a high growth rate, and which also offers value-added products from the residue, such as cosmetics, health functional food or pharmaceuticals. Microalgae are best cultivated in photo-bioreactors (PBRs) where light, nutrients, carbon dioxide and temperature can be controlled. Despite the current availability of PBRs, only a few can be practically used for mass production. Computational fluid dynamics (CFD) was used in this study to design an optimum bubble-column PBR for mass production of microalgae. Multi-phase models including bubble movement, meshes and time step independent tests were considered to develop the 3-dimensional CFD model. Particle Image Velocimetry (PIV) tests were used to enhance and ...

International Nuclear Information System (INIS)

Numerous mercury intrusion porosimetry (MIP) studies have been carried out to investigate the pore structure in cement-based materials. However, the standard MIP often results in an underestimation of large pores and an overestimation of small pores because of its intrinsic limitation. In this paper, an innovative MIP method is developed in order to provide a more accurate estimation of pore size distribution. The new MIP measurements are conducted following a unique mercury intrusion procedure, in which the applied pressure is increased from the minimum to the maximum by repeating pressurization-depressurization cycles instead of a continuous pressurization followed by a continuous depressurization. Accordingly, this method is called pressurization-depressurization cycling MIP (PDC-MIP). By following the PDC-MIP testing sequence, the volumes of the throat pores and the corresponding ink-bottle pores can be determined at every pore size. These values are used to ...

Energy Technology Data Exchange (ETDEWEB)

In the year 1996 the spheres of competence of several DLR-Institutes working in the areas of fluid dynamics, reaction kinetics, combustion, numerical methods and laser measuring techniques have been brought together while contributing to the internal DLR project 'combustion chamber modelling (BKM)', in order to proceed with the computational simulation of combustion processes in combustion chambers of gas turbines. The main issue was the development of a research code for numerical simulation of fluid flow in real combustion chambers. Here the development of computational models of physical and chemical processes was emphasized, among other processes the formation of soot was treated. Moreover, a worldwide outstanding database of measured data for the purpose of code validation has been created within the framework of the BKM project using the laboratory facilities ...

International Nuclear Information System (INIS)

The thermal environment of a hypothetical large indoor hydrocarbon pool fire is more complex compared to outdoor fires and can be more severe for engulfed objects. In order to analyze potential thermal environments for interim storage of spent fuel casks or low-level radioactive waste packages engulfed in pool fires numerical simulations with the CFD fire code FDS3 were carried out for different storage configurations. In addition, data of indoor pool fire experiments were used to validate the model for this type of application. A series of pool fire experiments under different ventilation conditions and varied pool surface (1 m"2 - 4 m"2) inside a compartment of 3.6 m x 3.6 m x 5.7 m was conducted at iBMB (Institut fuer Baustoffe, Massivbau und Brandschutz) of Braunschweig University of Technology, Germany. The instrumentation included thermocouples, heatflux and pressure gauges, bi-directional flow probes and gas ...

Science.gov (United States)

terms of 11) numerical models, 2D SPH simulations, and analytic models. The 1D model was useful for determining the evolution over very long timescales. ...

Science.gov (United States)

... being pulled out of water, the movement is restricted ... This heating depended on the particle size ... hydrometeor injection into the upper tropopause. ...

International Nuclear Information System (INIS)

This paper reports on Monte Carlo path integral simulations of a Bianchi Type IX cosmology which yield a wave function which reflects the known chaotic dynamics of the classical system and predicts a significant probability for a correlation between large universe volume and large anisotropy.

Science.gov (United States)

In this paper a new procedure, addressed as Interpolation Damage Detecting Method (IDDM), is investigated as a possible mean for early detection and location of light damage in a structure struck by an earthquake. Damage is defined in terms of the accuracy of a spline function in interpolating the operational mode shapes (ODS) of the structure. At a certain location a decrease (statistically meaningful) of accuracy, with respect to a reference configuration, points out a localized variation of the operational shapes thus revealing the existence of damage. In this paper, the proposed method is applied to a numerical model of a multistory frame, simulating a damaged condition through a reduction of the story stiffness. Several damage scenarios have been considered and the results indicate the effectiveness of the method to assess and localize damage for the case of concentrated damage and for low to medium levels of noise in the recorded signals. ...

Energy Technology Data Exchange (ETDEWEB)

Coupled modeling of subsurface multiphase fluid and heat flow, solute transport and chemical reactions can be used for the assessment of acid mine drainage remediation, mineral deposition, waste disposal sites, hydrothermal convection, contaminant transport, and groundwater quality. Here they present a numerical simulation model, TOUGHREACT, which considers non-isothermal multi-component chemical transport in both liquid and gas phases. A wide range of subsurface thermo-physical-chemical processes is considered. The model can be applied to one-, two- or three-dimensional porous and fractured media with physical and chemical heterogeneity. The model can accommodate any number of chemical species present in liquid, gas and solid phases. A variety of equilibrium chemical reactions is considered, such as aqueous complexation, gas dissolution/exsolution, cation exchange, and surface complexation. Mineral dissolution/precipitation can proceed either ...

Energy Technology Data Exchange (ETDEWEB)

Numerical simulations of compressible flows are commonly based on the Euler equations when effects of viscosity are thought to be negligible. These equations admit singular solutions, even in cases where the initial and boundary conditions are smooth. So-called ''Euler solvers'' rely on numerical dissipation, explicitly or implicitly present in the scheme, to regularize the problem, such that physical solutions are selected.

International Nuclear Information System (INIS)

According to the present concept, the low and intermediate level wastes generated during the Cernavoda NPP operation will be disposed in a near surface repository. The Saligny site, placed in the NPP protected area, has been proposed for their disposal. Geologically, the main components of this site are the quaternary loess, the Precambrian and Pre-quaternary clays, the Eocene and Barremian limestone. Hydrologically, the site can be divided into a vadose zone down to 45-50 m and three distinct aquifers, two of them in the limestone beds and the third in the lenses of sand and limestone existing in the pre-quaternary clay layer. A large research program for site characterization was initiated in 1996. At present, the site characteristics requested for safety analysis have been experimentally measured on soil samples or calculated by different computer programs. Hundreds of experimental values of the density, porosity, hydraulic conductivity, soil-water retention, moisture content or ...

Energy Technology Data Exchange (ETDEWEB)

In this study we tried to develop the long-range transport system and find the way to prevent from the radiological emergency risk. For the study, meteorological forecast system in Korea Meteorological Administration is investigated. Numerical simulation is also carried out by the long-range transport model and Vis-5D. We surveyed the emergency preparedness for nuclear accidents which were ARAC in USA, RODOS in Europe and WSPEED in Japan and then investigated the processing of medium- and long-range atmospheric diffusion modeling system. We also studied on the application of KMA/NWPD model which are GDAPS and RDAPS. In the future, it is necessary to produce to the high resolution meteorological data from KMA/NWPD for the development of medium- and long-range atmospheric diffusion modeling system and construct the integrated system for data processing in real time. It was simulated by using micro-scale meteorological field ...

Science.gov (United States)

Three-dimensional computer simulation and analyses of the horizontal steam generator thermal-hydraulics of the WWER 1000 nuclear power plant have been performed for 50% and 75% partial loads, 100% nominal load and 110% over-load. Presented results show water and steam mass flow rate vectors, steam void fraction spatial distribution, recirculation zones, swell level position, water mass inventory on the shell side, and other important thermal-hydraulic parameters. The simulations have been performed with the computer code 3D ANA, based on the 'two-fluid' model approach. Steam-water interface transport processes, as well as tube bundle flow resistance, energy transfer, and steam generation within tube bundles are modelled with {sup c}losure laws{sup .} Applied approach implies non-equilibrium thermal and flow conditions. The model is solved by the control volume procedure, which has been extended in order to take into account ...

International Nuclear Information System (INIS)

A numerical analysis method for melting/solidification phenomena has been developed to evaluate a feasibility of several candidate techniques in the nuclear fuel cycle. Our method is based on the eXtended Finite Element Method (X-FEM) which has been used for moving boundary problems. Key technique of the X-FEM is to incorporate signed distance function into finite element interpolation to represent a discontinuous gradient of the temperature at a moving solid-liquid interface. Construction of the finite element equation, the technique of quadrature and the method to solve the equation are reported here. The numerical solutions of the one-dimensional Stefan problem, solidification in a two-dimensional square corner and melting of pure gallium are compared to the exact solutions or to the experimental data. Through these analyses, validity of the newly developed numerical analysis method has been ...

Energy Technology Data Exchange (ETDEWEB)

Three-dimensional magnetohydrodynamic simulations are executed in a full toroidal geometry to clarify the physical mechanisms of the Internal Reconnection Event (IRE), which is observed in the spherical tokamak experiments. The simulation results reproduce several main properties of IRE. Comparison between the numerical results and experimental observation indicates fairly good agreements regarding nonlinear behavior, such as appearance of localized helical distortion, appearance of characteristic conical shape in the pressure profile during thermal quench, and subsequent appearance of the m=2/n=1 type helical distortion of the torus. (author)

Energy Technology Data Exchange (ETDEWEB)

A rigorous mathematical model is developed to describe chemical absorption of carbon dioxide in a falling-film micro-contactor. Contrary to previous models described in literature, hydrodynamics in both phases is accurately captured by the Navier-Stokes equations. By this means, backmixing effects can be accounted for. Furthermore, the interface position is determined directly using the level-set method. The developed model is validated with experimental data obtained from literature and a good agreement is established. With the validated model, detailed studies are performed to investigate the impact of process parameters.

Energy Technology Data Exchange (ETDEWEB)

Forests on the Kola Peninsula in Northern Russia grow close to the northern tree line. They are subjected to both natural and anthropogenic stress factors. The Cu-Ni smelter 'Severonikel' (Lat. 67 deg 55'N; Long. 32 deg 57'E) near Monchegorsk is one of the two major sources of sulphur dioxide and heavy metals emissions on the Kola Peninsula. These emissions have caused significant deterioration of the surrounding vegetation. The thesis demonstrates how methods of Remote sensing, ground survey and mathematical modelling can be integrated for monitoring of the smelter's environmental impact on the surrounding vegetation: ground truth data are used for calibration of remote-sensed data, which further serve to verify mathematical models. The study aims were: * to estimate the scale of airborne sulphur pollution from the smelting industry on the Kola Peninsula and its effect on vegetation; * to assess spatial extent of the forest decline in the ...

Energy Technology Data Exchange (ETDEWEB)

The mainly nonstationary operation of a solar-heated receiver can be simulated with sufficient accuracy only if data about the dynamic behavior are available. For this reason, the dynamic behavior of a solar cavity receiver with parabolic dish collector is investigated. The development of a mathematical simulation considering heat transfer and storage processes is presented and the procedure for a numerical solution is illustrated. The performance of the calculation method is finally demonstrated by simulating the passage of a cloud.

Science.gov (United States)

This paper focuses on a theoretical model that predicts the temperature increase of Magneto-rheological (MR) fluid dampers experiencing a sinusoidal input motion. A theoretical model is developed to estimate the temperature rise based on the non-linear behavior of the MR fluid damper. This model is solved numerically, and the numerical solution is compared with a known linear solution and experimental result in order to validate the accuracy of the model. Also, a non-dimensional form of the governing equations are developed to examine the key parameters. The non-dimensional terms show the effect of external and internal parameters on the trends of heat dissipation as well as heat generation within the MR fluid damper.

Energy Technology Data Exchange (ETDEWEB)

The present paper is concerned with development and application of a so-called Effective Convection Model (ECM), which aims to provide a detailed, mechanistic description of heat transfer processes in a BWR lower plenum. The ECM is a Computational Fluid Dynamics (CFD)-like tool which employs a simpler and more effective approach to compute heat transfer by solving only energy conservation equation instead of solving the full set of Navier-Stokes and energy equations by a CFD code. We implement the ECM in a CFD code (Fluent), with detailed description of the ECM development, implementation and validation. A dual approach is used to validate the ECM, namely validation against experimental data and against heat transfer results obtained by CFD predictions in the same geometries and conditions. Insights gained from CFD simulations are also used to improve ECM. The ECM capability as an effective tool to ...

Energy Technology Data Exchange (ETDEWEB)

X-ray ablation of material from the first wall and other components of an ICF (Inertial Confinement Fusion) chamber is a major threat to the laser final optics. Material condensing on these optics after a shot may cause damage with subsequent laser shots. To ensure the successful operation of the ICF facility, removal rates must be predicted accurately. The goal for this dissertation is to develop an experimentally validated x-ray response model, with particular application to the National Ignition Facility (NIF). Accurate knowledge of the x-ray and debris emissions from ICF targets is a critical first step in the process of predicting the performance of the target chamber system. A number of 1-D numerical simulations of NIF targets have been run to characterize target output in terms of energy, angular distribution, spectrum, and pulse shape. Scaling of output characteristics with variations of both target yield and ...

Science.gov (United States)

In the past few decades the need for improved nuclear reactor safety analyses has led to a rapid development of advanced methods for multidimensional thermal-hydraulic analyses. These methods have become progressively more complex in order to account for the many physical phenomena anticipated during steady state and transient Light Water Reactor (LWR) conditions. The advanced thermal-hydraulic subchannel code COBRA-TF (Thurgood, M. J. et al., 1983) is used worldwide for best-estimate evaluations of the nuclear reactor safety margins. In the framework of a joint research project between the Pennsylvania State University (PSU) and AREVA NP GmbH, the theoretical models and numerics of COBRA-TF have been improved. Under the name F-COBRA-TF, the code has been subjected to an extensive verification and validation program and has been applied to variety of LWR steady state and transient simulations. To enable F-COBRA-TF for ...

CERN Document Server

Heterogeneous anisotropic diffusion problems arise in the various areas of science and engineering including plasma physics, petroleum engineering, and image processing. Standard numerical methods can produce spurious oscillations when they are used to solve those problems. A common approach to avoid this difficulty is to design a proper numerical scheme and/or a proper mesh so that the numerical solution validates the discrete counterpart (DMP) of the maximum principle satisfied by the continuous solution. A well known mesh condition for the DMP satisfaction by the linear finite element solution of isotropic diffusion problems is the non-obtuse angle condition that requires the dihedral angles of mesh elements to be non-obtuse. In this paper, a generalization of the condition, the so-called anisotropic non-obtuse angle condition, is developed for the finite element solution of heterogeneous anisotropic ...

Energy Technology Data Exchange (ETDEWEB)

The interaction between a pair of counter-rotating vortices and a lean premixed CH{sub 4}/O{sub 2}/N{sub 2} flame ({Phi} = + 0.55) has been studied by direct numerical simulations using detailed and reduced chemical reaction schemes. Results from the complex chemistry simulation are discussed with respect to earlier experiments and differences in the simulations using detailed and reduces chemistry are investigated. Transient evolutions of the flame surface and the total heat release rate are compared and modifications in the evolution of the local flame structure are displayed. (authors) 22 refs.

Science.gov (United States)

This paper reports that, to obtain better simulation results for a Canada deuterium uranium (CANDU) reactor operation, a new simulation method is developed that uses actual detector readings as a correction factor. Detector readings from a CANDU reactor are used to correct the calculated flux distribution during core calculation iterations. A suitable function is found to describe the relationship between the detector flux and the fluxes of mesh points around the detector. The new simulation method is tested by performing numerical calculations for the Wolsung reactor (a CANDU-600). The results show that the new method predicts the core state more accurately with fewer iterations.

Energy Technology Data Exchange (ETDEWEB)

An important component of the overall program to validate WIMS-AECL for use with RFSP in the analysis of CANDU-6 reactors for design and safety analysis calculations is the validation of calculations of incremental cross sections used to represent reactivity devices. A method has been developed for the calculation of the three-dimensional neutron flux distribution in and around CANDU reactor fuel channels and reactivity control devices. The methods is based on one- and two dimensional transport calculations with the WIMS-AECL lattice cell code, SPH homogenization, and three-dimensional flux calculations with finite-difference diffusion theory using the MULTICELL code. Simulations of Wolsung 1 Phase-B commissioning measurements and Point Lepreau restart tests have been performed, as a part of the program to validate WIMS-AECL lattice cell calculations for application to CANDU reactor ...

British Library Electronic Table of Contents (United Kingdom)

A mechanical model with local resonators is proposed as an acoustic metamaterial that exhibits an unusual frequency-dependent effective stiffness. If treated as an equivalent elastic solid, its effective Young's modulus can become unbounded or vanishingly small at two respective frequencies. Moreover, in a certain frequency range, the effective Young's modulus would become negative, resulting in a band gap that coincides with this frequency range. The wave attenuation behavior and mechanism are studied through numerical simulations on the acoustic metamaterial model. The capability of the metamaterial to selectively block or filter unwanted waves is demonstrated by a numerical simulation example.

Energy Technology Data Exchange (ETDEWEB)

The formation of a quasiequilibrium beam distribution matched to an alternating-gradient quadrupole focusing lattice by means of the adiabatic turn-on of the oscillating focusing field is studied numerically using particle-in-cell simulations. Quiescent beam propagation over several hundred lattice periods is demonstrated for a broad range of beam intensities and vacuum phase advances describing the strength of the oscillating focusing field. Properties of the matched-beam distribution are investigated. In particular, self-similar evolution of the beam density profile is observed over a wide range of system parameters. The numerical simulations are performed using the WARP particle-in-cell code.

British Library Electronic Table of Contents (United Kingdom)

Results are given of numerical simulation of three-dimensional pattern of flow of a two-phase steam-water mixture in the house of a PGV-1000 horizontal steam generator obtained using the BAGIRA best-estimate thermohydrodynamic computer codes. The space distributions of velocities and local void fractions in the steam generator housing for different modes of operation of power-generating unit are calculated and compared with available experimental data.

British Library Electronic Table of Contents (United Kingdom)

The aim of this work is the development and experimental validation of a detailed dynamic fuel cell model using the gPROMS modeling environment. The model is oriented towards optimization and control and it relies on material and energy balances as well as electrochemical equations including semi-empirical equations. For the experimental validation of the model a fully automated and integrated hydrogen fuel cell testing unit was used. The predictive power of the model has been compared with the data obtained during load change experiments. A sensitivity analysis has been employed to reveal the most critical empirical model parameters that should be estimated using a systematic estimation procedure. Model predictions are in good agreement with experimental data under a wide range of operati...

CERN Document Server

We have carried out numerical simulations of strongly gravitating systems based on the Einstein equations coupled to the relativistic hydrodynamic equations using adaptive mesh refinement (AMR) techniques. AMR coalescences of neutron stars can now be simulated with sufficient resolution covering the neutron stars while having the computational domain extend to the local wave zone. We show an AMR simulation carried out with a workstation having an accuracy equivalent to that of a 1025^3 regular uni-grid simulation, which is, to the best of our knowledge, larger than all previous simulations of similar NS systems on supercomputers. We believe the capability opens new possibilities in general relativistic simulations.

CERN Document Server

We study analytically and numerically the stability of the pressure-less, viscously spreading accretion ring. We show that the ring is unstable to small non-axisymmetric perturbations. To perform the perturbation analysis of the ring we use a stretching transformation of the time coordinate. We find that to 1st order, one-armed spiral structures, and to 2nd order additionally two-armed spiral features may appear. Furthermore, we identify a dispersion relation determining the instability of the ring. The theoretical results are confirmed in several simulations, using two different numerical methods. These computations prove independently the existence of a secular spiral instability driven by viscosity, which evolves into persisting leading and trailing spiral waves. Our results settle the question whether the spiral structures found in earlier simulations of the spreading ring are ...

CERN Document Server

We present a new set of 3.5 Post-Newtonian equations in which Newtonian hydrodynamics is coupled to the nonconservative effects of gravitational radiation emission. Our formalism differs in two significant ways from a similar 3.5 Post-Newtonian approach proposed by Blanchet (1993, 1997). Firstly we concentrate only on the radiation-reaction effects produced by a time-varying mass-current quadrupole $S_{ij}$. Secondly, we adopt a gauge in which the radiation-reaction force densities depend on the fourth time derivative of $S_{ij}$, rather than on the fifth, as in Blanchet's approach. This difference makes our formalism particularly well-suited to numerical implementation and could prove useful in performing fully numerical simulations of the recently discovered $r$-mode instability for rotating neutron stars subject to axial perturbations.

International Nuclear Information System (INIS)

A three-dimensional (3D) thermal-hydraulic model and a numerical procedure for the simulation and analysis of a steady-state, as well as transient operation of nuclear power plant components are presented. A two-fluid approach is applied to modelling of two-phase flow. Thermal-hydraulics of a horizontal steam generator in the WWER 1000 nuclear power plant has been simulated at the full load, steady-state operation. A comparison of the numerical results with data measured at the NPP Novovoronjezh shows good agreement. 3D numerical results can be used in plant design or retrofitting, in nuclear power plant operation and safety analysis and as improvement of existing one-dimensional thermal-hydraulics models of the horizontal steam generator which are assessed by system codes used for the nuclear power plant safety analyses. (author)

Energy Technology Data Exchange (ETDEWEB)

We demonstrate a two-dimensional device simulator for MOSFET structures that incorporates models for defect distributions and show predicted effects on device switching performance for various spatial distributions of defects in amorphous and polycrystalline silicon.

International Nuclear Information System (INIS)

The model of transient enhanced diffusion of ion-implanted As is formulated and the finite-difference method for numerical solution of the system of equations obtained is developed. The nonuniform distribution of point defects near the interface and more accurate description of arsenic clustering are simultaneously taken into account. Simulation of As diffusion during rapid annealing gives a reasonable agreement with the experimental data. (authors)

Science.gov (United States)

A three-dimensional numerical simulation of the loss-of-feed water transient at the horizontal steam generator of the Kola nuclear power plant is performed. Presented numerical results show transient change of integral steam generator parameters, such as steam generation rate, water mass inventory, outlet reactor coolant temperature, as well as detailed distribution of shell side thermal-hydraulic parameters: swell and collapsed levels, void fraction distributions, mass flux vectors, etc. Numerical results are compared with measurements at the Kola NPP. The agreement is satisfactory, while differences are close to or below the measurement uncertainties. Obtained numerical results are the first ones that give complete insight into the three-dimensional and transient horizontal steam generator thermal-hydraulics. Also, the presented results serve as benchmark tests for the assessment ...

CERN Document Server

In this paper we present a model and simulator for many clusters of heterogeneous PCs belonging to a local network. These clusters are assumed to be connected to each other through a global network and each cluster is managed via a local scheduler which is shared by many users. We validate our simulator by comparing the experimental and analytical results of a M/M/4 queuing system. These studies indicate that the simulator is consistent. Next, we do the comparison with a real batch system and we obtain an average error of 10.5\\% for the response time and 12\\% for the makespan. We conclude that the simulator is realistic and well describes the behaviour of a large-scale system. Thus we can study the scheduling of our system called \\dirac in a high throughput context. We justify our decentralized, adaptive and opportunistic approach in comparison to a centralized approach in such a ...

Science.gov (United States)

The understanding of transport processes that occur in naturally fractured geothermal systems is far from being complete. Often, evaluation and numerical simulations of fractured geothermal reservoirs, are carried out by assuming equivalent porous media and homogeneous petrophysical properties within big matrix blocks. The purpose of this paper, is to present a comparison between results obtained from numerical studies of a naturally fractured reservoir treated as a simple porous medium and the simulation of some real aspects of the fractured reservoir. A general conclusion outlines the great practical importance of considering even approximately, the true nature of such systems. Our results show that the homogeneous simplified evaluation of the energy resource in a fractured system, could result in unrealistic estimates of the reservoir capacity to generate electricity.

Energy Technology Data Exchange (ETDEWEB)

The present paper describes the application of the commercial CFD-code FLUENT 6.1 to the numerical simulation of the flow field in a high specific speed radial fan. Two-dimensional as well as three-dimensional computations are performed using the steady 'Frozen-Rotor method' as well as the unsteady 'Sliding-Mesh method'. The predicted performance curves of the radial fan are compared with the results obtained from measurements on a test stand. Due to the relatively low computation times, the 'Frozen-Rotor method' can be used routinely for design purposes. The results of this approach can be improved if the leakage flow through the gap between the inlet nozzle and the impeller is taken into account. (orig.)

International Nuclear Information System (INIS)

Designs of nuclear reactor plants aim for high performance under safety consideration. Because of large scale and high pressure/temperature conditions, data from costly mockup tests have been required to verify simulation codes of systems and components. Establishment of design by analysis (DBA) in nuclear engineering is required for development of next generation nuclear reactors. Recent powerful computers and simulation technique enable numerical analyses to predict realistic behaviors of thermo-fluid flow, structure and do on. The present report describes resent simulation results of complex gas-liquid two-phase flow, large scale structure dynamics and fluid-structure interaction. (author)

Energy Technology Data Exchange (ETDEWEB)

There is an ever increasing demand to perform flow simulations that incorporate the complete details of geometry as well as sophisticated flow physics. This has led to the development of numerical algorithms that can simulate the actual flow phenomena with greater fidelity. However, the success of these algorithms hinges on the grid that models the geometry. Grid generation methods for 2-D models have long existed and the general lack of complexity of the simpler 2-D models has not quite challenged the efforts in this area. However, demands for generating better 3-D geometric models for flow simulations involving complex geometries have completely changed the perspective of grid generation strategies. As a consequence, grid generation efforts have earned equal significance as that of numerical solver efforts.

Energy Technology Data Exchange (ETDEWEB)

Hyogo-ken Nanbu earthquake with a focus in the Akashi straits has given huge earthquake damages in and around Awaji Island and Kobe City in 1995. It is clear that the basement structure, which is steeply deepened at Kobe City from Rokko Mountains towards the coast, and the focus under this related closely to the local generation of strong ground motion. Generation process of the strong ground motion was discussed using 2D and 3D numerical simulation methods. The 3D pseudospectral method was used for the calculation. Space of 51.2km{times}25.6km{times}25.6km was selected for the calculation. This space was discretized with the lattice interval of 200m. Consequently, it was found that the basement structure with a steeply deepened basement, soft and weak geological structure thickly deposited on the basement, and earthquake faults running under the boundary of base rock and sediments related greatly to the generation of strong ground motion. ...

CERN Document Server

We present a dynamical spectral model for Large Eddy Simulation of the incompressible magnetohydrodynamic (MHD) equations based on the Eddy Damped Quasi Normal Markovian approximation. This model extends classical spectral Large Eddy Simulations for the Navier-Stokes equations to incorporate general (non Kolmogorovian) spectra as well as eddy noise. We derive the model for MHD and show that introducing a new eddy-damping time for the dynamics of spectral tensors in the absence of equipartition between the velocity and magnetic fields leads to better agreement with direct numerical simulations, an important point for dynamo computations.

International Nuclear Information System (INIS)

The paper is concerned with the simulation of practical testing situations which are too difficult and/or expensive to replicate in a laboratory environment. Numerical experiments are described which simulate the differential eddy current probe response to the build-up and chemical flushing of magnetite in the crevice gap of a PWR steam generator unit. The simulation results agree well with the only experimental data available to the authors and lead to the conclusion that conventional differential eddy current probes should be capable of characterizing crevice gap conditions with respect to the presence of magnetite. (author).

Energy Technology Data Exchange (ETDEWEB)

The simulation of processes in engineering and the physical sciences has progressed rapidly over the last several years. With rapid developments in supercomputers, parallel processing, numerical algorithms and software, scientists and engineers are now positioned to quantitatively simulate systems requiring many billions of arithmetic operations. The need to understand and assimilate such massive amounts of data has been a driving force in the development of both hardware and software to create visual representations of the underling physical systems. In this paper, and the accompanying videotape, the evolution and development of the visualization process in scientific computing will be reviewed. Specific applications and associated imaging hardware and software technology illustrate both the computational needs and the evolving trends. 6 refs.

International Nuclear Information System (INIS)

Maps of wind resources were used to develop a method of evaluating climatic changes and their potential impacts on wind energy resources. Global IPCC climate change scenarios were used to predict climatic conditions for the future, while past wind resource availability was simulated and validated using NCEP and NCAR reanalysis data as well as observed meteorological data from Environment Canada. The simulations were used to compare each 5-year period with a 50-year reference period. Regional scale climate change impacts were evaluated using a statistical dynamic down-scaling method. Advanced meteorological models were used to predict wind flow patterns across specific landscapes. The evolution of past wind resource availability was then simulated. Five-year wind resource simulations for a 50-year period were simulated at 25 km"2 wind speeds at 80 m above the ...

International Nuclear Information System (INIS)

A Monte Carlo method called transient curve fitting method was developed to calculate a eigenvalues by first simulating the existing neutrons and precursors in the system, then calculating the Eigen-distribution of neutron flux and calculating the a eigenvalues using the transient results based on the Eigen-distribution by the code TMCC. The results of this method are tested by calculating Godiva Benchmark problems and they agree well with the benchmark results. Then the reasonable results of Subcritical Facility in Tsinghua University are given by TMCC, and the results are compared with the experimental results measured by Rossi-a method. Even in the deep subcritical cases, the method can give results consistent with experimental results. (authors)

International Nuclear Information System (INIS)

The applicability of conventional strength hypotheses in the long-term creep range was examined by means of multiaxial loaded small and large-scale specimens of 1% CrMoNiV-steel at 550deg C. The experimentally determined deformation behaviour in the primary and secondary creep range was simulated by means of a finite element analysis. The creep deformation capacity is influenced by the degree of multiaxiality. The validity of a strength hypothesis depends on the type of load and on deformation resp. failure criteria. (orig.).

Energy Technology Data Exchange (ETDEWEB)

This paper describes TENT, a system designed for integrating the various tools and parts of a CFD simulation system (CFD - Computational Fluid Dynamics) and the simulation code itself in one interactive simulation environment. Initially the system was intended to be a TEstbed for Numerical Turbine simulation but it turned out that the approach was general enough to be suitable for more complex simulation tasks in the field of aircraft simulation. The paper will give a brief description of a CFD simulation system and its component interactions, which are the basis for describing the Component Architecture based on CORBA. This includes the control architecture using an event model and data flow considerations. Because most CFD codes use parallelization techniques based on Message Passing, a so called MIMD Object ...

Science.gov (United States)

We have developed a novel laser ray-tracing method to measure aberrations in optical systems. It consists of delivering narrow laser pencils (by a laser scanner), recording the spots that are formed on the image plane (with a CCD camera), and computing the position of each centroid. This approach could be considered an experimental (approximate) implementation of standard numerical ray tracing. Several tests and experiments, including a direct comparison with a Hartmann-Shack wave-front sensor, provided highly satisfactory results that confirmed the validity of the method and revealed potential advantages. PMID:18073906

DEFF Research Database (Denmark)

A simple closed-form equation to calculate the envelope correlation between any two receiver or transmitter antennas in a multi-input multi-output (MIMO) system of an arbitrary number of elements is derived. The equation uses the scattering parameters obtained at the antenna feed point to calculate the envelope correlation coefficient. This approach has the advantage that it does not require knowledge of the antenna radiation pattern. Numerical data that include conductor and permittivity loss are shown to validate the approach. Using the scattering parameters for calculating the envelope correlation is less laborious and relatively cheaper, as compared to the radiation pattern approach.

Energy Technology Data Exchange (ETDEWEB)

The plasma plume generated by a new type of high energy Janitor known as the railplug, is examined. The railplug is a miniaturized railgun that has the potential for improving ignition characteristics of combustible mixtures in engines. The objective of the study is to gain an uderstanding of the characteristics of the plasma created by a transparent railplug, and to validate a multidimensional computer simulation of the plasma and shock fronts. The nature of the plume emitted by the railplug was examined for three levels of electrical energy while firing into air at a pressure of 1 atm. The computer model is to be used to predict trends in railplug performance for various railplug designs, energies, and ambient conditions. The velocity of the plasma movement inside a transparent railplug was measured, as well as the velocity of the plume ejected from the cavity. A shock is produced at the initiation point of the arc and propagates down the ...

Science.gov (United States)

This investigation concerns the nonlinear dynamics of heat transfer from a surface using an upstream eddy promoter. A numerical model is developed for the coupled fluid flow and heat transfer process based on a modified MacCormack scheme. Numerical simulations are carried out to determine the response and heat transfer enhancement due to the promoter. The average heat transfer from a cavity floor is seen to be increased by a factor of approximately five over the unpromoted'' flow. Another interesting feature of the study is the nonlinear viscous flow dynamics from the cylinder-wall interaction which differ significantly from the familiar cylinder-free stream patterns.

International Nuclear Information System (INIS)

The objective of this study is to establish evaluation and verification guideline for the APR 1400 and to investigate the thermal-hydraulic characteristics for fluidic device is analyzed using FLUENT. The scope and major results of research are flow characteristics for fluidic device. In this study, three-dimensional numerical model for fluidic device is developed adequately for, and results are compared with experimental data performed by VAPER (VAlve Performance Evaluation test Rig) in KAERI with an aim to verify numerical simulation. In addition, the parametric study has also carried out to investigate the effect of major parameters such as velocity and pressure inside FD chamber.

British Library Electronic Table of Contents (United Kingdom)

A new aiming point guidance (APG) law in 3-dimensional space is derived analytically. The guidance law is capable of guiding the missile (interceptor) to a straight collision course with respect to the target. The missile turning rate command proportional to the error angle between the current missile heading and the predicted aim point is calculated based on the predicted target position. Performance of the APG law is numerically evaluated and compared with the proportional navigation guidance (PNG) law. The result of numerical simulation shows that the APG law is superior to the PNG law in terms of miss distance, time of interception and final stage missile turning rate demand.

International Nuclear Information System (INIS)

Energy savings from the use of daylighting in commercial buildings are realized through implementation of photoelectric lighting controls that dim electric lights when sufficient daylight is available to provide adequate workplane illumination. The dimming level of electric lighting is based on the signal of a photosensor. Current simulation approaches for such systems are based on the questionable assumption that the signal of the photosensor is proportional to the task illuminance. This paper presents a method that simulates the performance of photosensor controls considering the acceptance angle, angular sensitivity, placement of the photosensor within a space, and color correction filter. The method is based on the multiplication of two fisheye images: one generated from the angular sensitivity of the photosensor and the other from a 180- or 360-degree fisheye image of the space as ''seen'' by the photosensor. The paper includes a detailed ...

DEFF Research Database (Denmark)

The determination of conformational preferences in unfolded and disordered proteins is an important challenge in structural biology. We here describe an algorithm to optimize energy functions for the simulation of unfolded proteins. The procedure is based on the maximum likelihood principle and employs a fast and efficient gradient descent method to find the set of parameters of the energy function that best explain the experimental data. We first validate the method by using synthetic reference data, and subsequently apply the algorithms to data from nuclear magnetic resonance spin-labeling experiments on the Delta 131 Delta fragment of Staphylococcal nuclease. A significant strength of the procedure that we present is that it directly uses experimental data to optimize the energy parameters, without relying on the availability of high resolution structures. The procedure is fully general and can be applied to a range of experimental data and ...

Energy Technology Data Exchange (ETDEWEB)

The impact of the azimuthal wave refraction in the middle atmosphere on the distribution of gravity wave amplitudes, propagation azimuths, and other wave parameters is investigated using a numerical ray-tracing model of gravity wave propagation through a representative zonal mean reference model of geostrophic winds and temperature in the middle atmosphere. Simulations are first performed with only a single type of gravity wave to help explain some important effects which occur during the refraction process. Then, a multiray simulation is performed which traces a crude spectrum of waves from different altitudes through the atmosphere for every month of a climatological year. The simulated wave climatologies are compared with observations. 108 refs.

CERN Document Server

We examine the observable properties of simulated barred galaxies including radial mass profiles, edge-on structure and kinematics, bar lengths and pattern speed evolution for detailed comparison to real systems. We have run several simulations in which bars are created through inherent instabilities in self-consistent simulations of a realistic disc+halo galaxy model with a disc-dominated, flat rotation curve. These simulations were run at high (N=20M particles) and low (N=500K) resolution to test numerical convergence. We determine the pattern speeds in simulations directly from the phase angle of the bar versus time and the Tremaine-Weinberg method. Fundamental dynamics do not change between the high and low resolution, suggesting that convergence has been reached in this case. We find the higher resolution is needed to simulate ...

Energy Technology Data Exchange (ETDEWEB)

Turbulent premixed methane-air combustion in a cylindrical chamber is numerically simulated considering radiation effect. Reaction rates are considered as minimum rates between Arrhenius rates and eddy break up rates. A five step reduced mechanism is used. Turbulent modeling is done via standard k-{epsilon} model imposed by empirical inlet boundary conditions. Source terms of energy equation consist of reaction rates and radiation effects. The discrete ordinate method (DOM) is employed to solve the radiative transfer equation (RTE) and the weighted sum of gray gas model (WSGGM) is imposed to consider radiation effect of non-gray gases. The results indicate that in the case of turbulent combusting flows, the effect of radiation of gases can affect the temperature and species concentrations. The numerical results obtained considering radiation effect are closer to the experimental data than that of the case without radiation ...

International Nuclear Information System (INIS)

Results of three-dimensional laminar and standard K-#epsilon# turbulent numerical simulations of natural convection cooling of ten cubic aluminum blocks mounted on an insulated plate, facing a shrouding wall, are presented. This geometry is chosen so that comparison with experimental results is possible. The considered problem is of great practical importance because it simulates the case of heated electronic chips, mounted on printed board assemblies, which are frequently encountered in electronic industry applications. The problem is mathematically modeled by the three-dimensional conservation differential equations of mass, momentum, energy and turbulent kinetic energy and dissipation (for the turbulent flow model). IN this paper, these equations are numerically solved by a finite volume method and the laminar and turbulent results are compared to the experimental results obtained with similar ...

CERN Document Server

Dark matter particles form halos that contribute the major part of the mass of galaxy clusters. The formation of these cosmological structures have been investigated both observationally and in numerical simulations, which have confirmed the existence of a universal mass profile. However, the dynamic behaviour of dark matter in halos is not as well understood. We have used observations of 16 equilibrated galaxy clusters to show that the random velocities of dark matter particles are larger on average along the radial direction than along the tangential, and that the magnitude of this velocity anisotropy is radially varying. Our measurement implies that the collective behaviour of dark matter particles is fundamentally different from that of normal particles and the radial variation of the anisotropy velocity agrees with the predictions of numerical simulation.

British Library Electronic Table of Contents (United Kingdom)

In this research, the interactive effect of grain and specimen sizes on the flow stress of sheet metal in microforming is investigated via the tensile test of pure copper and numerical modeling. Models based on different assumptions are proposed to analyze the size effect phenomenon. It is found that the flow stress decreases linearly with the decrease of the ratio of specimen to grain sizes. The grain boundary thickness decreases and its volume fraction increases with the decrease of grain size. The variation of grain boundary thickness is not proportional to the variation of grain size. Furthermore, the fraction of grain boundary increases with the strain and the ratio of specimen to grain sizes. Based on the FE simulation, it is found that the simulated flow stress, which is modeled bas...

British Library Electronic Table of Contents (United Kingdom)

The FVM-LES-acoustic analogy method (FVM-LES-AAM), which is a hybrid prediction technique for the acoustical property computation, is presented and performed in this paper. The FVM-LES-AAM was developed by combining the finite volume method (FVM), the large eddy simulation (LES), and the Ffowcs Williams-Hawkings analogy algorithm (FWH-AA). To predict the acoustical properties of induction cookers, the FVM is used for discretizing the calculation field and building numerical equations, and the LES and FWH-AA are performed for computing the sound sources and predicting the far-field sound, respectively. Using the FVM with the unstructured grids method to discretize the control equation of Navier-Stokes was introduced for illuminating the above numerical simulation procedure. To prove the FVM...

International Nuclear Information System (INIS)

This study deals with testing the capability of the code PARET to simulate natural convection cooling phenomena under different boundary conditions. In addition to applying and testing some new options related to simulation of the control rod movement and studying the reactivity effect of thermal expansion fuel elements. The experiments of the simple thermal hydraulic loop of Missouri university about natural cooling phenomena in two narrow paralled channels were used to validate the code. The study indicate good results regarding the distribution of coolant flux velocity and clad temperature. In particular the heat transfer coefficient of natural convection has been calculated in good agreement with the experiment. On the other hand, the core of MNSR reactor has been modelled to simulate the reactor dynamic behaviour under natural convection cooling conditions for different initial power level. The ...

Science.gov (United States)

A novel approach to fault diagnosis is proposed using multiscale morphology analysis to extract impulsive features from the signals with strong background noise. Multiscale morphology is applied to one-dimensional signal by defining both the length and height scales of structuring elements (SEs). A local-peak-value based adaptive algorithm is also introduced. The new approach makes the selection of SEs more transparent and is independent of empirical rules. Both simulated impulsive and vibration signals of two defective roller bearings are employed to validate the proposed algorithm. The roller bearing faults presented in the validation include both inner and outer race faults. The test results show that the multiscale morphology analysis is effective and robust to extract morphological features.

British Library Electronic Table of Contents (United Kingdom)

Mechanobiological theories have been introduced to illustrate the interaction between biology and the local mechanical environment during bone healing. Although several theories have been proposed, a quantitative validation using histomorphometric data is still missing. In this study, in vivo histological data based on an ovine animal experiment was quantified and used to validate bone healing simulations focussing on the endochondral ossification process. The bone formation at different callus regions (periosteal and endosteal bone at the medial and lateral side) was analyzed from in vivo data and quantitatively compared with in silico results. A histomorphometric difference was found in medial and lateral hard callus formation 3 weeks after osteotomy. However, the same amount of new bone...

International Nuclear Information System (INIS)

Fluid dynamic analysis of a commercial, counter-flow Ranque-Hilsch Vortex Tube (Rh vt), Ex air 25 s cfm, has been performed in this work both experimentally and numerically; in particular Rh vt cooling power and temperature separation performances have been tested in both direct cooling employment (jet impingement) and indirect cooling employment (supplying cold plates). Experimental techniques, used in this work, revealed several difficulties to produce detailed information about velocity and temperature fields inside the tube and at both the exits. Hence numerical simulation of the flow inside the tube has been conducted using the commercial Cfd code Fluent 6.3.26. Compressible, turbulent, high swirling flow inside Rh vt has been simulated by using both Rans and Les approaches. In particular several turbulence closures have been used in the Rans simulations and results have been ...

Science.gov (United States)

We develop a 3-D model to simulate the synthetic aperture radar (SAR) image formation process of an undulated vegetation canopy such as corn grown in fields with large periodic drainage reliefs. We explain how the simulated SAR image of undulated vegetation medium is obtained by the convolution of a 2-D slice of the 3-D simulated SAR system point spread function [(PSF), emulating the SAR beam modeled by a cosine modulated Gaussian], with the 2-D projection of the observed undulated vegetation canopy (modeled with scatterers randomly distributed in 3-D undulated space) followed by the extraction of each look envelope, the summation of looks, and sampling in azimuthal and range directions. Our model is useful to study the parameters involved in the formation and the analysis of SAR images of undulated vegetation medium. Validation of simulations made with actual SAR images shows that ...

International Nuclear Information System (INIS)

Three dimensional temperature field and helium flow field of TBM are simulated using the general purpose computational fluid dynamics (CFD) code FLUENT. The temperature distribution of Be Armor, Be Pebble Bed, Li_4SiO_4 Pebble Bed, Structure Material of TBM, and helium flow field in the cooling pipe are presented. The research indicates that the work temperature of each material is under the material temperature allowed except some places where high temperature should be excluded in the design. The results will provide references for further optimized thermal hydraulic design of ITER China TBM. (authors)

British Library Electronic Table of Contents (United Kingdom)

Due to the low thermal conductivity of ceramics large temperature gradients are present through the powder particles during plasma spray deposition. As a result the particles often impinge at the substrate in a semi-molten form; which in turn substantially affects the final characteristics of the coating being formed. This study is dedicated to the novel modelling development and simulation of a semi-molten droplet impingement. The study examines the impingement process during impact, spreading and solidification of semi-molten zirconia. The simulation provides an insight to the heat transfer process during impact and solidification of a semi-molten powder particle and illustrates the freezing-induced break-up mechanism at the splat periphery.

International Nuclear Information System (INIS)

The mechanism proposed by Kormendy (1984) for the formation of counterrotating cores in elliptical galaxies is investigated using self-consistent numerical simulations of mergers between a high- and a low-luminosity elliptical galaxies. The conditions for a counterrotation to appear are determined, observational properties of the remnants are described, and the evolution of the structural and kinematic parameters of the larger galaxy is analyzed. It is shown that a counterrotation results only when the merging orbits are retrograde, due to a large change in the secondary spin during the merger. 36 refs.

Energy Technology Data Exchange (ETDEWEB)

To simulate the effects of the breakup channel on rearrangement amplitudes, the conventional coupled-reaction-channel (CRC) expansion is augmented by pseudoreaction channels. The construction of the projector for the extended CRC space is discussed, and transition-operator equations on this space are given. By solving the full and post-approximation forms of the CRC equations for a model three-particle problem, the crucial role played by the nonorthogonality terms is demonstrated.

Energy Technology Data Exchange (ETDEWEB)

This paper suggests how nonlinear periodic optimal control of a pumped storage plant can be realized. The control problem consists in maximization of the plant benefits over an operational period. An optimal control law is proposed, yielding a bang-off-bang mode of operation. An algorithm for numerical solution of the problem was developed, and its effectiveness was demonstrated by simulation experiments.

Energy Technology Data Exchange (ETDEWEB)

In this letter we report a clear and unambiguous observation of the out-of-plane quadrupole magnetic field suggested by numerical simulations in the reconnecting current sheet in the Magnetic Reconnection Experiment (MRX). Measurements show that the Hall effect is large in collisionless regime and becomes small as the collisionality increases, indicating that the Hall effect plays an important role in collisionless reconnection.

International Nuclear Information System (INIS)

A distinguishing design feature of CANDU nuclear reactors is the use of horizontal fuel channels housed in a horizontal vessel called a calandria, which is made of stainless steel 304L. Each channel consists of a Zr-2.5%Nb alloy pressure tube and an externally concentric Zr-2 calandria tube. The calandria tubes are joined to the end plates (tubesheets) of the calandria vessel by joints formed by roller expansion. The bores in the tubesheets are grooved. Roller expanded joints provide a cost effective means of joining dissimilar materials, require minimal space and no maintenance. The quality of these roller expanded joints is important from a sealing, strength and stress corrosion point of view. The roller expansion process consists of expanding the calandria tubes to deform them plastically against the bores and into the grooves of the tubesheets. Therefore, understanding the effect of the number, geometry and the pitch of the grooves on the quality of a roller expanded joint is very ...

Energy Technology Data Exchange (ETDEWEB)

There is a crescent interest in the scientific community in the study of natural circulation phenomenon. New generation of compact nuclear reactors uses the natural circulation of the fluid as a system of cooling and of residual heat removal in case of accident or shutdown. The objective of this paper is to present a study through the comparison of experimental data and numerical simulation for the natural circulation phenomenon in one and two-phase flow regime. An experimental circuit built with glass tubes is used for the experiments. Thus, it allows the thermal hydraulic phenomena visualization. There is an electric heater as the heat source, a heat exchanger as the heat sink and an expansion tank to accommodate fluid density excursions. The circuit instrumentation consists of thermocouples and pressure meters to better keep track of the flow and heat transfer phenomena. Instrumentation data acquisition is performed through a computer ...

International Nuclear Information System (INIS)

The influence that electromagnetic fields have on mechanical and electronic equipment was discussed. A three dimensional numerical simulation model of electric fields around high voltage power transmission lines was described. 4 refs., 1 tab., 5 figs.

International Nuclear Information System (INIS)

The Accelerator Transmutation Waste System is under development at the Los Alamos National Laboratory. The goal is to perform an independent verification of the feasibility of actinide and long-lived fission product burning in this system. The authors' work is divided into five tasks: (a) production of an actinide and long-lived fission product cross section library from JEF 2.2; (b) simulation using MCNP and KENO IV Monte Carlo codes, of the Accelerator Transmutation Waste System configurations existing in literature; (c) validation of HETC Monte Carlo code (production of spallation source); (d) validation of the cross sections by comparison of Keff and reaction rate results, calculated with MNCP and KENO IV, with experimental benchmarks and intercomparison between the authors' calculations of a PWR unit cell and the computations carried out with various codes and cross section libraries (NEACRP critically group data); and ...

Energy Technology Data Exchange (ETDEWEB)

A common practice of electric utilities experiencing transmission and distribution (T and D) system overloads is to expand the substation, add lines, or upgrade equipment, all of which are capital intensive options. In 1988, it was hypothesized that strategically sited photovoltaics (PV) could benefit parts of T and D systems near or at overloaded conditions. An evaluation methodology was developed and applied to a test case (Kerman Substation near Fresno, California). Analytical results suggested that the value of PV to the T and D system could substantially exceed its energy and generation capacity value. The importance of this finding indicated the need for empirical validation. This led to the construction of a 0.5 MW PV demonstration plant by Pacific Gas and Electric Company (PG and E) at Kerman, California as part of the PVUSA (PV for Utility Scale Applications) project. PVUSA is a national public-private partnership that is assessing and demonstrating the ...

CERN Document Server

An effective computer program for three dimensional relativistic hydrodynamical model has been developed. It implements a new approach to the early hot phase of relativistic heavy-ion collisions. The computer program simulates time-space evolution of nuclear matter in terms of ideal-fluid dynamics. Equations of motions of hydrodynamics are solved making use of finite difference methods. Commonly-used algorithms of numerical relativistic hydrodynamics RHLLE and MUSTA-FORCE have been applied in simulations. To speed-up calculations, parallel processing has been made available for solving hydrodynamical equations. The test results of simulations for 3D, 2D and Bjorken expansion are reported in this paper. As a next step we plan to implement the hadronization algorithm by implementing the continuous particle emission for freeze-out and comparing it with Cooper-Frye formula.

Energy Technology Data Exchange (ETDEWEB)

This contribution is an informal essay based on a talk delivered at the Institute for Mathematics and its Applications (IMA) in Minneapolis, under the summer program in molecular biology, July 18-22, 1994. I exclude many technical details, which can be found elsewhere, and instead focus on the basic ideas of molecular dynamics simulations, with the goal of conveying to students and non-specialists the key concepts of the theory and practice of large-scale simulations. Following a description of the basic idea in molecular dynamics, I discuss some of the practical details involved in simulations of large biological molecules, the numerical timestep problem, and approaches to this problem based on implicit-integration techniques. I end with a perspective of open challenges in the field and directions for future research. 79 refs., 6 figs., 1 tab.

International Nuclear Information System (INIS)

The energy deposited by the passage of a single, energetic, heavy-ion through a semiconductor produces dense electron-hole (eh) pair concentrations near the ion trajectory. The size, shape, and charge density of an ion track represent critical parameters for many models of single event phenomena. The authors describe the design and uses of possible semiconductor test structures for measuring the initial radial distribution of charge and subsequent charge transport in a high energy, heavy-ion track. Numerical simulations show how the test structure can resolve different radial distributions of charge within an ion track. The test structure simulations also show the importance of accurately representing ion track structure in single event effects simulations.

Energy Technology Data Exchange (ETDEWEB)

Current projections of the response of the biosphere to global climatic change indicate as much as 50% to 90% spatial displacement of extratropical biomes. The mechanism of spatial shift could be dominated by either (1) competitive displacement of northern biomes by southern biomes, or (2) drought-induced dieback of areas susceptible to change. The current suite of global biosphere models cannot distinguish between these two processes, thus determining the need for a mechanistically based biome model. The model is in an early stage of development and will require several enhancements, including explicit simulation of potential evapotranspiration, extension to boreal and tropical biomes, a shift from steady-state to transient dynamics, and validation on other continents.

Energy Technology Data Exchange (ETDEWEB)

The following paper is focused on the numerical modelling of the transient radon diffusion through radon-proof membranes during the measurement of their radon diffusion coefficient. The major aim of such numerical modelling is to increase the accuracy of radon diffusion coefficients derived from the measured data sets. The developed complex ''transient'' numerical model is able to calculate the radon diffusion coefficient with sufficient accuracy from almost any data set - even from a short-time measurement with a non-linear course of results. This numerical model can also be used for various analyses of transient radon transfer processes (e.g. for the calculation of radon distribution curves within the membrane). The following paper presents governing equations for the simulation model, together with a brief description of algorithms incorporated ...

British Library Electronic Table of Contents (United Kingdom)

Eight selected sulfonamide drugs were investigated as inhibitors of heat shock protein 90 (Hsp90). The investigation included simulated docking experiments to fit the selected compounds within the binding pocket of Hsp90. The selected molecules were found to readily fit within the ATP-binding pocket of Hsp90 in low-energy poses. The sulfonamides torsemide, sulfathiazole, and sulfadiazine were found to inhibit the ATPase activity of Hsp90 with IC50 values of 1.0, 2.6, and 1.5 mM, respectively. Our results suggest that these well-established sulfonamides can be good leads for subsequent optimization into potent Hsp90 inhibitors.

British Library Electronic Table of Contents (United Kingdom)

To form a licensing basis for a new methodology for a fuel channel safety analysis code for CANDU-6 nuclear reactor, a CATHENA model for a post-blowdown fuel channel analysis has been developed, and tested for a high temperature thermal-chemical experiment CS28-1 [Lei, Q.M., 1993. Post-test analysis of the 28-element high-temperature thermal-chemical experiment CS28-1. In: 4th International Conference on Simulation Methods in Nuclear Engineering, Montreal, PQ, 1993]. Pursuant to the objective of this investigation, the current study has focused on understanding the involved phenomena, their interrelations, and how to maintain a good accuracy of the temperature and H2 generation rate prediction without losing the important physics of the involved phenomena. The transient simulation results ...

International Nuclear Information System (INIS)

This study deals with testing the capacity of the code PARET to simulate natural circulation phenomena under different boundary conditions in addition to assessment of some new options related to simulation of control rod movement and the reactivity effect of thermal expansion fuel elements. the experiments of the simple thermal hydraulic loop of Missouri University about natural circulation phenomena in narrow parallel channel were used to validate the code. The results indicate good agreements regarding the evolution of coolant velocity and clad temperature. In particular the heat transfer coefficient of natural convection has been calculated in good agreement with the experiment. On the other hand, the core of MNSR reactor has been modelled to stimulate the reactor dynamic behaviour under natural circulation condition for different initial power level. The observed oscillations during the initial phase vanish gradually ...

Science.gov (United States)

A computational fluid dynamics (CFD) is coupled with a computational structural dynamics (CSD) to simulate the unsteady rotor flow with aeroelasticity effects. An unstructured upwind Navier-Stokes solver was developed for this simulation, with 2nd order time-accurate dual-time stepping method for temporal discretization and low Mach number preconditioning method. For turbulent flows, both the Spalart-Allmaras and Menter's SST model are available. Mesh deformation is achieved through a fast dynamic grid method called Delaunay graph map method for unsteady flow simulation. The rotor blades are modeled as Hodges & Dowell's nonlinear beams coupled flap-lag-torsion. The rotorcraft computational structural dynamics code employs the 15-dof beam finite element formulation for modeling. The structure code was validated by comparing the natural frequencies of a rotor model with UMARC. The flow and structure ...

Science.gov (United States)

In this dissertation, a two-phase, air-droplets, dilute, turbulent, and steady state flow in a horizontal rectangular channel, is modeled and numerically simulated using a modified KIVA-3V code. The deposition of different sizes of droplets on the walls of the channel is also studied. In this model, the interaction effects between the phases (two-way coupling) are considered by source terms in the momentum and energy equations for the continuous phase and by the instantaneous local velocity of the air in the droplet equation of motion, which includes the aerodynamic and gravitational forces. The turbulence is modeled by a k-? model. The interaction effects between the turbulence and the dispersed droplets are also taken into account. The effects of the turbulence on the droplets are modeled by a fluctuating component added to the local air velocity in the droplet equation of motion. The effects of the droplets on the turbulence are modeled by ...

Energy Technology Data Exchange (ETDEWEB)

Precondition for the low-NO operation of gas turbine burners is the generation of as homogeneous as possible a gaseous fuel/air mixture. Such a burner behavior can be achieved by selective adaptation of the flow and injection. A computing model is introduced that calculates the dynamically balanced GASs/droplet flow through numerical solution of the transport equations in a curved orthogonal coordinate system. The feedback reaction of the fuel droplets to the gas phase is considered by source terms and/or sink terms. Test computations were carried out for comparison with analytical solutions from the flow mechanics and validated by comparison with measurement results. Droplet motion and evaporation were checked by means of published experimental results on single droplets. For the validation of the developed computing technique, model tests were carried out with water in place of fuel. It is pointed out that the choice of ...

CERN Document Server

We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the IGM, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, evolving from a sub-Mpc scale at the beginning of reionization to R>10 Mpc towards the end. We present a detailed statistical description of our results, and ...

CERN Document Server

(Abridged) We use 3D SPH calculations with higher resolution, as well as with more realistic viscosity and sound-speed prescriptions than previous work to examine the eccentric instability which underlies the superhump phenomenon in semi-detached binaries. We illustrate the importance of the two-armed spiral mode in the generation of superhumps. Differential motions in the fluid disc cause converging flows which lead to strong spiral shocks once each superhump cycle. The dissipation associated with these shocks powers the superhump. We compare 2D and 3D results, and conclude that 3D simulations are necessary to faithfully simulate the disc dynamics. We ran our simulations for unprecedented durations, so that an eccentric equilibrium is established except at high mass ratios where the growth rate of the instability is very low. Our improved simulations give a closer match to the observed relationship ...

Energy Technology Data Exchange (ETDEWEB)

The transport of radon in concrete takes place through the complicated network of interconnected pores that is, at any time, the result of the process of hydration of cement and of moisture distribution and transport. Initially the microstructure of concrete depends on the mix proportions and curing conditions, its time-evolution being conditioned by its surrounding environment. Radon transport will be consequently a function of time, as it is influenced by the changing microstructure (total porosity and its distribution) and by the amount and distribution of the moisture contained in the pore system. A selection of information from the large amount of literature available on concrete is presented in chapter 2. A model that describes the process of hydration, of microstructure development and of moisture transport is presented in chapter 3. The physics of radon diffusion in homogeneous porous materials is outlined in chapter 4. The coupling of the numerical ...

International Nuclear Information System (INIS)

A numerical and experimental investigation is carried out in a solar thermochemical reactor for the thermal dissociation of ZnO at 2000 K using concentrated solar energy. The reactor consists of a cavity-receiver lined with ZnO particles and directly exposed to high-flux irradiation. A transient heat transfer model is formulated to link the rate of radiation, convection, and conduction heat transfer to the reaction kinetics. The radiosity and Monte Carlo methods are applied to obtain the distribution of net radiative fluxes at the internal surfaces of the reactor cavity and at the surface of the ZnO bed. Validation is accomplished in terms of the calculated and measured transient temperature profiles and chemical reaction rates.

British Library Electronic Table of Contents (United Kingdom)

A non-similar boundary layer analysis is presented to study the flow, heat and mass transfer characteristics of non-Darcian mixed convection of a non-Newtonian fluid from a vertical isothermal plate embedded in a homogeneous porous medium with the effect of Soret and Dufour and in the presence of either surface injection or suction. The value of the mixed-convection parameter lies between 0 and 1. In addition, the power-law model is used for non-Newtonian fluids with exponent n1 for dilatant fluids. Furthermore, the coordinates and dependent variables are transformed to yield computationally efficient numerical solutions that are valid over the entire range of mixed convection, from the pure forced-convection limit to the pure free-conve...

British Library Electronic Table of Contents (United Kingdom)

The offshore and inshore wave energy resource in Asturias (N Spain) is studied using wave buoy data and a hindcast dataset spanning 44 years (1958-2001). Offshore average wave power and annual wave energy values are found to exceed 30 kW/m and 250 MWh/m, respectively, at 7 of the 11 study sites. This substantial resource is characterised in terms of the sea states involved. Most of the energy is provided by IV quadrant waves with significant wave heights between 2 m and 5 m and energy periods between 11 s and 13 s. After analysing the offshore resource, numerical modelling is used to investigate the inshore wave patterns. A coastal wave model is validated with wave buoy data and applied to three case studies representative of storm, winter and summer conditions. Inshore wave energy concent...

British Library Electronic Table of Contents (United Kingdom)

The heating of water layer using microwave oven with a rectangular waveguide has been studied both numerically and experimentally. The mathematical model is validated with the experimental data. The transient Maxwell's equations are solved by using the Finite Difference Time Domain (FDTD) method to describe the electromagnetic field inside the waveguide and sample. The temperature profile and velocity field within sample are determined by the solutions of the momentum, energy and Maxwell's equations. In this study, the effects of physical parameters, e.g. microwave power level, placement of sample inside the waveguide, volume of sample, are studied. The distribution of electric field, temperature profile and velocity field are presented in details. The results show good agreement between s...

CERN Document Server

We obtain an elegant and useful description of the dynamics of Szekeres dust models (in their full generality) by means of "quasi--local" scalar variables constructed by suitable integral distributions that can be interpreted as weighed proper volume averages of the local covariant scalars. In terms of these variables, the field equations and basic physical and geometric quantities are formally identical to their corresponding expressions in the spherically symmetric LTB dust models. Since we can map every Szekeres model to a unique LTB model, rigorous results valid for the latter models can be readily generalized to a non--spherical Szekeres geometry. The new variables lead naturally to an initial value formulation in which all scalars are expressed as scaling laws in terms of their values at an arbitrary initial space slice. These variables also yield a significant simplification of numerical work, since the fluid flow evolution equations ...

Energy Technology Data Exchange (ETDEWEB)

We intend to solve equations governing turbulent plane-vertical isotherm and non isotherm jets by taking into account inflow conditions at the exit of the nozzle. The analysis is focused on the influence of these conditions on this type of flow. Two cases are considered (uniform and parabolic velocity and temperature profiles). A finite difference scheme is developed to solve the governing equations. This numeric model allows us to show that the region of fully developed regime begins much nearer the nozzle for the turbulent case than for the laminar flow case. Indeed, the turbulence increases the mixing between the incoming gas from the nozzle and the ambient fluid, and consequently the size of the potential core zone decreases. The results are compared to other works introducing mathematical variables based on the energy conservation for the case of the mixed convection and the momentum conservation for the forced convection, which allows the ...

International Nuclear Information System (INIS)

A detailed treatment is introduced to measure the dynamic stability of the relativistic electrons in a self-amplified spontaneous emission free-electron laser (FEL) system, which includes the numerical approach of the Kolmogorov entropy (entropy-like quantity), the general equations of motion for a charged particle and the method of monitoring the simulation accuracy. Numerical experiments reveal a new phenomenon that there exists the possibility of the transition from chaotic to non-chaotic phase-space trajectories of the strongly relativistic electrons due to the effect of their self-fields. The adiabatic magnetic field of a one-dimensional wiggler may have a slight influence on the electron transportation in the absence of the FEL fields, but substantially affects the dynamic stability of the electrons in the process of the FEL interaction. Moreover, the laser fields diminish the dynamic stability of the electrons as the ...

CERN Document Server

The Fourier transform of cosmological density perturbations can be represented in terms of amplitudes and phases for each Fourier mode. We investigate the phase evolution of these modes using a mixture of analytical and numerical techniques. Using a toy model of one-dimensional perturbations evolving under the Zel'dovich approximation as an initial motivation, we develop a statistic that quantifies the information content of the distribution of phases. Using numerical simulations beginning with more realistic Gaussian random-phase initial conditions, we show that the information content of the phases grows from zero in the initial conditions, first slowly and then rapidly when structures become non-linear. This growth of phase information can be expressed in terms of an effective entropy: Gaussian initial conditions are a maximum entropy realisation of the initial power spectrum, gravitational evolution decreases the phase ...

Energy Technology Data Exchange (ETDEWEB)

Combined gas and steam turbine processes based on direct coal firing show a high thermal efficiency. At RWTH Aachen, University of Technology, an experimental test furnace has been built to investigate the pressurized pulverized coal combustion (PPCC). The PPCC-facility has been constructed as a slag tap furnace. Particles hitting the walls at temperatures above the melting point cause slagging depositions and create a film flowing down the reactor walls. As a part of the PPCC-program different mathematical models have been developed and implemented into the CFD-code FLUENT to predict the behavior of slag films at the furnace walls. Numerical strategies and the mathematical models used are described in detail. 12 refs., 9 figs.

Energy Technology Data Exchange (ETDEWEB)

The scope of the present article is two-fold. Firstly, to conduct an experiment to provide the temperature-time history of the cooling of a hot ball bearing in quiescent ambient air. Secondly, to predict the temporal variation of the bearing under the hypothesis of natural convection, radiation or natural convection coexists with radiation for a non-vanishing total hemispherical emissivity of the surface of the bearing. Numerical solutions of the three governing nonlinear lumped heat equations were carried out with a Runge-Kutta-Fehlberg (RKF45) algorithm accounting for automatic step size control. The experimental data was obtained with chrome steel ball bearings of diameter 0.953 cm (7/16 in) heated in an electric oven to a pre-set temperature. The heated bearing was exposed later to ambient air at atmospheric temperature and pressure. (orig.)

British Library Electronic Table of Contents (United Kingdom)

A test bench was designed and assembled to carry out impact tests on samples and components. The system allows simple and rapid adjustment of the test parameters, such as the shot size and air pressure, with good repeatability of the results. Tests on steel and light alloys were carried out under both as-produced condition and on thin hard-coated samples. Significant reductions in dimple dimensions were seen after coating. FE models simulating the experiments overestimated the dimple depths, although the parameter trend was satisfactorily captured. The residual stresses from coating and shot peening determined numerically are believed to have been proven effective against fatigue.

DEFF Research Database (Denmark)

We present a theoretical study of the resonant interaction between dynamical localized states (discrete breathers) and linear electromagnetic excitations (EE's) in Josephson junction ladders. By making use of direct numerical simulations we find that such an interaction manifests itself by resonant steps and various sharp switchings (voltage jumps) in the current-voltage characteristics. Moreover, the power of ac oscillations away from the breather center (the breather tail) displays singularities as the externally applied dc bias decreases. All these features may be mapped to the spectrum of EE's that has been derived analytically and numerically. Using an improved analysis of the breather tail, a spectroscopy of the EE's is developed, The nature of breather instability driven by localized EE's is established.

Energy Technology Data Exchange (ETDEWEB)

The bulk removal of CO/sub 2/ from industrial gases is a frequently applied process in the field of gas-treating. Usually chemical solvents are used in order to increase both the solvent-loading capacity and the absorption rate. These solvents often are aqueous solutions of alkanolamines. The loaded-liquid is regenerated at high temperatures with the aid steam. In the present study absorption rates are experimentally determined in various mixtures of amines over a wide range of concentrations. Simultaneously, a model was developed for the description of mass transfer accompanied by parallel reversible chemical reactions. The equations of this mass transfer model could not be solved analytically and therefore numerical techniques were used to obtain an exact description of this process. This model also includes the asymptotic situations of homogeneous catalysts and shuttle mechanism for which in the literature several approximate solutions have been presented. The ...

Energy Technology Data Exchange (ETDEWEB)

The Stockpile Stewardship Program (SSP) is a single, highly integrated technical program for maintaining the surety and reliability of the U.S. nuclear stockpile. The SSP uses past nuclear test data along with current and future non-nuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. It includes stockpile surveillance, experimental research, development and engineering (D&E) programs, and an appropriately scaled production capability to support stockpile requirements. This integrated national program requires the continued use of current facilities and programs along with new experimental facilities and computational enhancements to support these programs. The Advanced Simulation and Computing Program (ASC) is a cornerstone of the SSP, providing simulation capabilities and computational resources to support the annual stockpile assessment ...

Energy Technology Data Exchange (ETDEWEB)

During the last years one focus of German PWR safety analysis was boron dilution events with the potential of reactivity transients. Coolant with a low boron concentration could be collected in localized areas of the reactor coolant system e.g. by separation of borated reactor coolant into highly concentrated and diluted fractions (inherent dilution) which can occur during reflux- condenser heat transfer after a small break loss of coolant accident with a limited availability of the emergency core cooling systems. The TUeV NORD SysTec was charged by German supervisory authorities with the assessment of the safety analyses of boron dilution events presented by the utilities. These analyses are based on the simulation of boron dilution and transport processes in conjunction with a number of dedicated experiments. The analyses shall demonstrate that boron dilution events cannot lead to recriticality of the core. Hence the boron concentration at the core inlet has to ...

Energy Technology Data Exchange (ETDEWEB)

We have developed a kinetic Monte Carlo (kMC) simulator that links atomic migration and binding energies determined primarily from first principles calculations to macroscopic phenomena and laboratory time scales. Input for the kMC simulation is obtained from a combination of ab initio planewave pseudopotential calculations, molecular dynamics simulations, and experimental data. The simulator is validated against an extensive series of experimental studies of the diffusion of B spikes in self-implanted Si. The implant energy, dose, and dose rate, as well as the detailed thermal history of the sample, are included. Good agreement is obtained with the experimental data for temperatures between 750 and 950 C and times from 15 to 255 s. At 1050o C we predict too little diffusion after 105 s compared to experiment: apparently, some mechanism which is not adequately represented by our ...

International Nuclear Information System (INIS)

The applicability of the thermal response of an electrically heated simulated rod mostly used in loss-of-coolant-accident (LOCA) experiments to that of a nuclear fuel rod is a concern for the safety evaluation of a reactor. The present analysis describes the characteristics of the thermal response for both electrically heated and nuclear fuel rods during typical reflood conditions for a PWR-LOCA. A model describing the radial temperature field in the rod is developed based on the scheme in HETRAP code by Malang and incorporated into a reflood analysis code, REFLA for that purpose. The calculations applied to the existing reflood tests gave good agreement with experiments, showing the validity of the present model. The analysis has shown that the nuclear fuel rod tends to give a lower clad temperature and a sooner quench time than the electrically heated rod in a typical reflood condition, due to the smaller gap heat transfer and smaller heat ...

CERN Document Server

Several observational studies suggest that solar wind dynamic pressure fluctuations can drive magnetospheric ultra-low frequency (ULF) waves on the dayside. To investigate this causal relationship, we present results from Lyon-Fedder-Mobarry (LFM) global, three-dimensional magnetohydrodynamic (MHD) simulations of the solar wind-magnetosphere interaction. These simulations are driven with synthetic solar wind input conditions, where idealized ULF dynamic pressure fluctuations are embedded in the upstream solar wind. In three of the simulations, a monochromatic, sinusoidal ULF oscillation is introduced into the solar wind dynamic pressure time series. In the fourth simulation, a continuum of ULF fluctuations over the 0-50 mHz frequency band is introduced into the solar wind dynamic pressure time series. In this numerical experiment, the idealized solar wind input conditions allow us ...

International Nuclear Information System (INIS)

The current trend towards the increased use of risk assessment in the regulation of nuclear power plants will inevitably result in changes in the analysis of fire protection systems and the methods of analysis. Before fire protection can be regulated on a risk basis, a consensus must be reached on a number of issues. One key issue is what types of computational tools will be allowable for analyzing fire events, and what types of scenarios those tools will be approved for use. Reaching this consensus will require an understanding of the types of computational tools available and their inherent advantages and disadvantages. To aid with this understanding, three different methods of fire simulation are applied to an oil pool fire test in the HDR (Heiss Dampf Reaktor) containment test facility. These methods are a hand calculation, the zone model code CFAST (Consolidated Model of Fire Growth and Smoke Transport), and the computational fluid dynamics code FDS (Fire ...

Energy Technology Data Exchange (ETDEWEB)

E-35 reported the specific air flow analysis program using a CFSV model (including a human body model) to develop a high-accuracy evaluation simulator for heterogeneous thermal environment based on human thermal physiological conditions. E-36 reported development of the radiation heat transfer analysis module by using a numerical integral method to calculate shape factors. Radiation heat transfer characteristics every body part were clarified quantitatively by this module. E-37 reported the results on continuous measurement of physiological conditions and thermal environment factors every monitor under controlled indoor thermal environment, to collect measurement data necessary for improving the accuracy of a prediction program for human body temperatures. E-38 reported the study result on physiological reaction and subjective evaluation through an examinee experiment in a blood controlled region, to collect basic data for the above ...

International Nuclear Information System (INIS)

In this paper, we calculated the radioactive concentration distribution of radioactive waste water, the temperature distribution of drained cooling water and the effect of implement from the Daya Bay Nuclear Power Plant on nearby waters range, discussed and analysed some problems of computational results and computation with Alternating Direction Implicit Method (ADI). The contents of the article included: the establishment of two-dimension tidal current equation, radioactive waste water pollutant dispersion equation and cooling water heat convection diffusion equation, the numerical difference calculation model of tidal current field, concentration field as well as temperature field, effect impingement with ADI method, numerical calculation results. The result of research showed that: when the Daya Bay Nuclear Power Plant is on normal operation and after the low level radioactive waste water and low temperature cooling water have been drained ...

CERN Document Server

A Langmuir film is a molecularly thin film on the surface of a fluid; we study the evolution of a Langmuir film with two co-existing fluid phases driven by an inter-phase line tension and damped by the viscous drag of the underlying subfluid. Experimentally, we study an 8CB Langmuir film via digitally-imaged Brewster Angle Microscopy (BAM) in a four-roll mill setup which applies a transient strain and images the response. When a compact domain is stretched by the imposed strain, it first assumes a bola shape with two tear-drop shaped reservoirs connected by a thin tether which then slowly relaxes to a circular domain which minimizes the interfacial energy of the system. We process the digital images of the experiment to extract the domain shapes. We then use one of these shapes as an initial condition for the numerical solution of a boundary-integral model of the underlying hydrodynamics and compare the subsequent images of the experiment to the ...

International Nuclear Information System (INIS)

The RELAP5 code has been developed for best-estimate transient simulation of light water reactor coolant systems during postulated accidents. The code models the coupled behavior of the reactor coolant system and the core for loss-of-coolant accidents and operational transients such as anticipated transient without scram, loss of offsite power, loss of feedwater, and loss of flow. A generic modeling approach is used that permits simulating a variety of thermal hydraulic systems. Control system and secondary system components are included to permit modeling of plant controls, turbines, condensers, and secondary feedwater systems. RELAP5/MOD3 code documentation is divided into seven volumes: Volume I presents modeling theory and associated numerical schemes; Volume II details instructions for code application and input data preparation; Volume III presents the results of developmental assessment cases that demonstrate and ...

Energy Technology Data Exchange (ETDEWEB)

A numerical model was developed to simulate the combined effects of heat transfer, magnetite oxidation, and NO{sub x} formation in a grate-kiln furnace for taconite pellet induration. Heat transfer from the flame in the kiln was described by the net radiation method. The shrinking core model was used to account for magnetite oxidation on the grate. A novel approach to oxidation of tumbling pellets in a kiln was derived. The Zeldovich mechanism was used to predict thermal NO generation. Temperature fluctuations in the gas streams were estimated with a clipped Gaussian probability density function. The thermal energy balances and mass balances resulted in coupled systems of first-order differential equations, which were solved numerically. The model is capable of predicting NO production and magnetite oxidation in agreement with observation of plant performance. Although the design of the grate-kiln system is for efficient ...

British Library Electronic Table of Contents (United Kingdom)

Abstract We present the first numerical,-N-body, hydrodynamical, chemical simulations of cosmic structure formation in the framework of non-Gaussian models. We study the impact of primordial non-Gaussianities on early chemistry (e-, H, H+, H-, He, He+, He++, H2, H+2, D, D+, HD, HeH+), molecular and atomic gas cooling, star formation, metal (C, O, Si, Fe, Mg, S) enrichment, Population-III (popIII) and Population-II-I (popII) transition and on the evolution of -visible- objects. We find that non-Gaussianities can have some consequences on baryonic structure formation at very early epochs, but the subsequent evolution at later times washes out any difference among the various models. When assuming reasonable values for primordial non-Gaussian perturbations, it turns out that they are responsi...

International Nuclear Information System (INIS)

A force measurement technique has been developed for large-scale aerodynamic models with a short test time. The technique is based on direct acceleration measurements, with miniature accelerometers mounted on a test model suspended by wires. Measuring acceleration at two different locations, the technique can eliminate oscillations from natural vibration of the model. The technique was used for drag force measurements on a 3 m long supersonic combustor model in the HIEST free-piston driven shock tunnel. A time resolution of 350 #mu#s is guaranteed during measurements, whose resolution is enough for ms order test time in HIEST. To evaluate measurement reliability and accuracy, measured values were compared with results from a three-dimensional Navier-Stokes numerical simulation. The difference between measured values and numerical simulation values was less than 5%. We conclude that this measurement ...

Energy Technology Data Exchange (ETDEWEB)

This report documents a one-dimensional numerical model, SAMFT1D, developed to simulate single-phase and multiphase fluid flow and solute transport in variably saturated porous media. The formulation of the governing equations and the numerical procedures used in the code for single-phase and multiphase flow and transport are presented. The code is constructed to handle single-phase as well as two or three-phase flow conditions using two integrated sets of computational modules. The fully implicit scheme is used in the code for both single-phase and multiphase flow simulations. Either the Crank-Nicholson scheme or the fully implicit scheme may be used in the transport simulation. The single-phase modules employ the Galerkin and upstream weighted residual finite element techniques to model flow and transport of water (aqueous phase) containing dissolved single-species contaminants ...

British Library Electronic Table of Contents (United Kingdom)

This work aims to contribute to the understanding of the influence of the ionospheric layer height (ILH) on the thin layer ionospheric model (TLIM) used to retrieve ionospheric information from the GNSS observations. Particular attention is paid to the errors caused on the estimation of the vertical total electron content (vTEC) and the GNSS satellites and receivers inter-frequency biases (IFB), by the use of an inappropriate ILH. The work relies upon numerical simulations performed with an empirical model of the Earth?s ionosphere: the model is used to create realistic but controlled ionospheric scenarios and the errors are evaluated after recovering those scenarios with the TLIM. The error assessment is performed in the Central and the northern part of the South American continents, a re...

CERN Document Server

After two decades of direct dynamical simulation of large-scale structure in the universe, it is safe to say the subject is now mature. Still, there are parts of the problem that are less well developed than others. In general, the collisionless dynamics of the dark matter component is better understood than the collisional gas dynamics of the baryonic component. In situations where the gas dynamics is relatively simple, such as the Lyman-$\\alpha$ forest and the intracluster medium in X-ray clusters, our ability to reproduce observational data has evolved rapidly, and the interpretive and predictive power of such experiments should now be taken seriously. A comparison of twelve gas dynamic codes to the problem of forming a single X-ray cluster shows that numerical inaccuracies are modest (typically below ten percent), leaving missing physics as the main source for large systematic differences between theory and observation. Galaxy formation, ...

Science.gov (United States)

Recently the author showed that the widely used simulation code TDA3D, even though a single frequency code, can be used to determine the power spectrum in the SASE process with excellent approximation in the exponential growth regime. In this paper, the author applies this method to the BNL Cornell Wiggler A SASE experiment as an example. When the gain is not very high, there are many modes in the radiation, which seems to make the analytical calculation very difficult. However, he shows that the increment of the radiation due to SASE over the spontaneous radiation can be expanded in terms of guided modes with rapid convergence. Thus when the spontaneous radiation is substracted from the SASE power during the calculation, there is a good agreement between the analytical theory and the numerical simulation.

CERN Document Server

We investigate the influence of the turbulence forcing on the mass distributions of gravitationally unstable cores by postprocessing data from simulations of non-selfgravitating isothermal supersonic turbulence with varying resolution. In one set of simulations solenoidal forcing is applied, while the second set uses purely compressive forcing to excite turbulent motions. From the resulting density field, we compute the mass distribution of gravitationally unstable cores by means of a clump-finding algorithm. Using the time-averaged probability density functions of the mass density, semi-analytic mass distributions are calculated from analytical theories. We apply stability criteria that are based on the Bonnor-Ebert mass resulting from the thermal pressure and from the sum of thermal and turbulent pressure. Although there are uncertainties in the application of the clump-finding algorithm, we find systematic differences in the mass ...

British Library Electronic Table of Contents (United Kingdom)

The effects of thermal radiation and viscous dissipation on magneto-hydrodynamic (MHD) unsteady free-convection flow over a semi-infinite vertical porous plate are analysed. The fluid considered is non-gray (absorption coefficient dependent on wave length). The Network Simulation Method is used to solve the boundary-layer equations based on the finite-difference formulation; only discretization of the spatial co-ordinates is necessary, while time remains as a real continuous variable. This method provides a solution for both transient and steady-state problems at the same time, and programming does not require manipulation of the sophisticated mathematical software that is inherent in other numerical methods. The velocity, temperature, local skin-friction and local Nusselt number are studi...

British Library Electronic Table of Contents (United Kingdom)

A suite of numerical experiments is conducted to extend our previous studies to explore the effects of linear, ambient wind shear on simulated mammatus-like clouds. Increasing values of unidirectional shear result in banded cloud bases, and for the strongest shears roll-type quasi-two-dimensional convection results rather than three-dimensional lobe morphology. Two observed soundings were used, and the sounding with a drier, shallower sub-cloud layer appeared to be most strongly affected by the presence of the ambient shear. The drier sounding also had less snow in the mammatus lobes owing to sublimation. Copyright Copyright 2009 Royal Meteorological Society

International Nuclear Information System (INIS)

The 540 MWe PHWR has two fully independent shutdown systems. The first shutdown system (SDS-1) comprises of 28 spring assisted, vertical gravity drop shut-off rods, each consisting of a cadmium absorber sandwiched between stainless steel tubes. The second shutdown system (SDS-2) constitutes six poison tanks connected to respective zircaloy injection tubes. This system is capable of high speed injection of gadolinium nitrate solution (in D2O) into the moderator through these tubes. Theoretical estimation was carried out at different injection pressures and different concentration of gadolinium nitrate solution to arrive at the limiting value of these parameters from reactivity consideration point of view. The plant measurements of SDS-2 actuations at 60 and 80 Kg/cm2 pressure of helium was used to validate and upgrade the estimation model. The paper gives the details of the validation details of SDS-2 actuation. (author)

International Nuclear Information System (INIS)

The Melt Vessel Interaction (MVI) project is concerned with the consequences of the interactions that a core melt, generated during a postulated severe accident in a light water reactor, may have with the pressure vessel. In particular, the issues concerned with the failure of the vessel bottom head are the focus of the research. The specific objectives of the project are to obtain data and develop validated models, which could be applied to prototypic plants, and accident conditions, for resolution of issues related to the melt vessel interactions. The project work has been performed by nine partners having varied responsibility. The work included a large number of experiments, with simulant materials, whose observations and results are employed, respectively, to understand the physical mechanisms and to develop validated models. Applications to the prototypic geometry and conditions have also been performed. This report ...

Energy Technology Data Exchange (ETDEWEB)

Some diagnostics at the National Ignition Facility (NIF), including the Gamma Reaction History (GRH) diagnostic, require multiple channels of data to achieve the required dynamic range. These channels need to be stitched together into a single time series, and they may have non-uniform and redundant time samples. We chose to apply the popular cubic smoothing spline technique to our stitching problem because we needed a general non-parametric method. We adapted one of the algorithms in the literature, by Hutchinson and deHoog, to our needs. The modified algorithm and the resulting code perform a cubic smoothing spline fit to multiple data channels with redundant time samples and missing data points. The data channels can have different, time-varying, zero-mean white noise characteristics. The method we employ automatically determines an optimal smoothing level by minimizing the Generalized Cross Validation (GCV) score. In order to automatically ...

British Library Electronic Table of Contents (United Kingdom)

The load-carrying capacity and ultimate deformation capacity of timber-concrete composite systems can be significantly influenced by the ductility of the connection between the two materials. The use of more ductile connections can increase the load-carrying capacity of the composite system as well as its ultimate deformation capacity. In this work, the potential increase that might be expected for these two parameters due to the use of ductile connections will be assessed through numerical simulations, taking the non-linear behaviour of the connections into account. Furthermore, the connection ductility required to achieve the maximum load capacity depends on the mechanical properties of the connection as well as on the geometric and mechanical properties of the composite system. There ar...

Science.gov (United States)

This paper presents a systematic study on the stability of a two-dimensional vibration isolation floating raft system with a time-delayed feedback control. Based on the generalized Sturm criterion, the critical control gain for the delay-independent stability region and critical time delays for the stability switches are derived. The critical conditions can provide a theoretical guidance of chaotification design for line spectra reduction. Numerical simulations verify the correctness of the approach. Bifurcation analyses reveal that chaotification is more likely to occur in unstable region defined by these critical conditions, and the stiffness of the floating raft and mass ratio are the sensitive parameters to reduce critical control gain.

International Nuclear Information System (INIS)

An extension of the Hodgkin-Huxley mathematical model for the propagation of nerve signal which takes into account dynamical heat transfer in biological tissue is derived and fine tuned with existing experimental data. The medium is heated by Joule's effect associated with action potential propagation, leading to characteristic thermal patterns in association with spiral and scroll waves. The introduction of heat transfer-necessary on physical grounds-provides a novel way to directly observe the movement, regular or chaotic, of the tip of spiral waves in numerical simulations and possibly in experiments regarding different biological excitable media.

British Library Electronic Table of Contents (United Kingdom)

The Monte Carlo ray-tracing method is applied and coupled with optical properties to predict the radiation performance of solar concentrator/cavity receiver systems. Several different cavity geometries are compared on the radiation performance. A flux density distribution measurement system for dish parabolic concentrators is developed. The contours of the flux distribution for target placements at different distances from the dish vertex of a solar concentrator are taken by using an indirect method with a Lambert and a charge coupled device (CCD) camera. Further, the measured flux distributions are compared with a Monte Carlo-predicted distribution. The results can be a valuable reference for the design and assemblage of the solar collector system.

British Library Electronic Table of Contents (United Kingdom)

A method to achieve NMR of dilute samples in the earth's magnetic field by applying para-hydrogen induced polarization is presented. Maximum achievable polarization enhancements were calculated by numerically simulating the experiment and compared to the experimental results and to the thermal equilibrium in the earth's magnetic field. Simultaneous 19F and 1H NMR detection on a sub-milliliter sample of a fluorinated alkyne at millimolar concentration (1018 nuclear spins) was realized with just one single scan. A highly resolved spectrum with a signal/noise ratio higher than 50:1 was obtained without using an auxiliary magnet or any form of radio frequency shielding.

Energy Technology Data Exchange (ETDEWEB)

Thermal striping is identified as one of the causes of thermal fatigue failure in nuclear power plants. Numerical studies of thermal striping require three-dimensional, unsteady turbulent modeling that resolves both large and small-scale turbulent motions. Benchmark studies were carried out using the LES turbulence model solved by the commercial CFD code FLUENT. Two types of mixing tee configurations were modeled to evaluate the performance of the CFD code. The simulation results presented in normalized average temperature and normalized fluctuating temperatures are in good agreement with measurements.

CERN Document Server

We study the production of gravitational waves from cosmic domain walls created during phase transition in the early universe. We investigate the process of formation and evolution of domain walls by running three dimensional lattice simulations. If we introduce an approximate discrete symmetry, walls become metastable and finally disappear. We calculate the spectrum of gravitational waves produced by collapsing metastable domain walls. Extrapolating the numerical results, we find the signal of gravitational waves produced by domain walls whose energy scale is around 10^10-10^12GeV will be observable in the next generation gravitational wave interferometers.

CERN Document Server

We consider the adiabatic demagnetization in the rotating reference frame (ADRF) of a system of dipolar coupled nuclear spins $s=1/2$ in the external magnetic field. The demagnetization starts with the offset of the external magnetic field (in frequency units) from the Larmor frequency being several times greater than the local dipolar field. For different subsystem sizes, we have found from numerical simulations the temperatures at which subsystems of a one-dimensional nine-spin chain and a plane nine-spin cluster become entangled. These temperatures are of the order of microkelvins and are almost independent of the subsystem size. There is a weak dependence of the temperature on the space dimension of the system.

Energy Technology Data Exchange (ETDEWEB)

This book of proceedings reports on 7 papers on combustion chemistry and formation of pollutants presented during the workshop organized by the `Combustion and Flames` section of the French society of thermal engineers. The chemistry of combustion is analyzed in various situations such as: turbojet engines, spark ignition engines, industrial burners, gas turbines etc... Numerical simulation is used to understand the physico-chemical processes involved in combustion, to describe the kinetics of oxidation, combustion and flame propagation, and to predict the formation of pollutants. (J.S.)

Energy Technology Data Exchange (ETDEWEB)

This conference is composed of 21 communications and 21 posters in the domain of building ventilation and indoor air quality; the main themes are: indoor air quality assessment and optimization; performance enhancement and optimization of ventilation systems and equipment; ventilation systems for renovated and rehabilitated buildings; French and European regulations, standardizations and certifications; experimental and numerical simulation studies concerning ventilation systems, air flow, temperature distribution, air quality, thermal comfort and ventilation-induced acoustic levels

International Nuclear Information System (INIS)

The asymmetric fingered structure for polysilicon thin-film transistors (AF-TFTs) is analysed in detail by combining experimental characteristics and two-dimensional numerical simulations. This structure allows an effective reduction of the kink effect and off-current, without introducing any additional series resistance. In addition, a substantial improvement in the device stability is also observed when compared to conventional TFT. The AF-TFT characteristics have been explained by considering a two transistor model.

Energy Technology Data Exchange (ETDEWEB)

Diffusion of radon in a bed of activated charcoal is described by diffusion equations. An analytical solution of these equations is presented for the case of constant radon concentration in the atmosphere. The solutions are given separately for short term and long term exposure. An analytical form of the calibration constant f for long term exposure and constant radon concentration in air, was found. A numerical method and computer code based on the method of finite elements is developed for the case of variable radon concentration in air. This program simulates radon adsorption by the activated charcoal bed, enabling determination of sensitivity. The dependence of sensitivity on different parameters, such as temperature, thickness of the charcoal, etc. was studied using this program. (orig.) 10 refs.

International Nuclear Information System (INIS)

This paper describes the use of numerical analysis for studying natural convection cooling systems for long term storage of heat producing radioactive materials, including special nuclear materials and nuclear waste. The paper explains the major design philosophy, and shares the experiences of numerical modeling. The strategy of storing radioactive material is to immobilize nuclear high-level waste by a vitrification process, convertion it into borosilicate glass, and cast the glass into stainless steel canisters. These canisters are seal welded, decontaminated, inspected, and temporarily stored in an underground vault until they can be sent to a geologic repository for permanent storage. These canisters generate heat by nuclear decay of radioactive isotopes. The function of the storage facility ventilation system is to ensure that the glass centerline temperature does not exceed the glass transition temperature during storage and the vault ...

Energy Technology Data Exchange (ETDEWEB)

A new air injection process called toe-to-heel air injection (THAI) has been developed to achieve very stable combustion performance in heavy oil wells by controlling gas override. In the THAI process, injected air migrates preferentially to the combustion front into a horizontal producer well by the oxygen diffusional gradient created under steady state conditions. This is possible because of the balance between the stoichiometric oxygen flux needed to burn the fuel (coke) laid down ahead of the combustion front and the removal of downstream combustion gases and mobilised fluids. This study involved a numerical simulation of the THAI process to increase the understanding of the processes involving fluid flow, heat transfer and chemical reactions taking place during in situ combustion. Two tests were conducted involving one dry and one wet combustion using heavy Wolf Lake Oil. The objective was to determine if the behaviour of the oil could be ...

Energy Technology Data Exchange (ETDEWEB)

In radiation therapy with hadron beams, conformal irradiation to a tumour can be achieved by using the properties of incident ions such as the high dose concentration around the Bragg peak. For the effective utilization of such properties, it is necessary to evaluate the volume irradiated with hadron beams and the deposited dose distribution in a patient's body. Several methods have been proposed for this purpose, one of which uses the positron emitters generated through fragmentation reactions between incident ions and target nuclei. In the previous paper, we showed that the maximum likelihood estimation (MLE) method could be applicable to the estimation of beam end-point from the measured positron emitting activity distribution for mono-energetic beam irradiations. In a practical treatment, a spread-out Bragg peak (SOBP) beam is used to achieve a uniform biological dose distribution in the whole target volume. Therefore, in the present paper, we proposed to extend the MLE ...

Energy Technology Data Exchange (ETDEWEB)

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). Note that this includes support for most popular parallel and serial computers; (2) Improved performance for all numerical kernels (e.g., time integrator, nonlinear and linear solvers) through state-of-the-art algorithms and novel techniques. (3) Device models which are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users ...

Energy Technology Data Exchange (ETDEWEB)

The groundwater flow and transport model of the Faultless underground nuclear test conducted at the Central Nevada Test Area (CNTA) was accepted by the state regulator and the environmental remediation efforts at the site have progressed to the stages of model validation and long-term monitoring design. This report discusses the long-term monitoring strategy developed for CNTA. Subsurface monitoring is an expensive and time-consuming process, and the design approach should be based on a solid foundation. As such, a thorough literature review of monitoring network design is first presented. Monitoring well networks can be designed for a number of objectives including aquifer characterization, parameter estimation, compliance monitoring, detection monitoring, ambient monitoring, and research monitoring, to name a few. Design methodologies also range from simple hydrogeologic intuition-based tools to sophisticated statistical- and optimization-based tools. When ...

Energy Technology Data Exchange (ETDEWEB)

Fugitive dust emissions from stockpiles in the open storage yards of industrial zones and the subsequent atmospheric dust dispersion have brought about many ecological and economical problems. This paper introduces a new approach to reduce pollutant emissions by simply changing pile configuration and without affecting the operational activity on these parks. Flow around piles of different geometries and for various wind conditions was studied using previously validated Computational Fluid Dynamic (CFD) simulations. Different pile height scenarios were investigated corresponding to a constant material volume and a fixed angle of repose under various wind magnitudes. The results obtained were integrated in order to evaluate the rate of dust emissions for the various configurations considered. It was found that, for the range of wind conditions and pile dimensions tested, an optimal geometry exists which corresponds to the lower emissions rate. 13 ...

British Library Electronic Table of Contents (United Kingdom)

In this study, the thermodynamic stability of the grain boundaries and the grain growth of nanocrystalline Palladium (Pd) at various temperatures were investigated. For this purpose, the Gibbs free energy curves of grain boundaries were plotted in terms of the excess volume by the use of the equation of state (EOS) and Song's thermodynamic models. The results showed that, according to the prediction of these models, the nanocrystalline growth in metals was stopped at the grain sizes less than the critical grain size. Also, the results of the temperature variations and its effect on the Gibbs free energy curves showed that by the increase of the temperature, the possibility for the stoppage of grain growth is facilitated and the critical grain size is increased. To investigate the validity ...

Energy Technology Data Exchange (ETDEWEB)

A method has been developed to optimize the energy resolution of a horizontally focusing monochromator. The method consists of determining the optimum radius of curvature of the cylindrically bent monochromator crystal by minimizing the total diffracted X-ray flux measured through an X-ray absorption foil. When measured at an absorption edge a global minimum can be identified, which corresponds to the minimum energy band accepted along the entire length of the crystal. Experimental verification of this method has been validated by comparing X-ray fluorescent scans taken for a series of crystal curvatures and by directly measuring the X-ray beam profiles at the corresponding points. The actual optical configuration and asymmetric-cut parameter of the focusing monochromator on beamline X4C at the National Synchrotron Light Source are modeled and ray tracing simulations of the optical system are compared with direct beam profile measurements.

International Nuclear Information System (INIS)

Dispersed flow film boiling is the heat transfer regime that occurs at high void fractions in a heated channel. The way this transfer mode is modelled in the NRC computer codes (RELAP5 and TRAC) and the validity of the assumption and empirical correlations used is discussed. An extensive review of the theoretical and experimental work related with heat transfer to highly dispersed mixtures reveals the basic deficiencies of these models: the investigation refers mostly to the typical conditions of low rate bottom reflooding, since the simulation of this physical situation by the computer codes has often showed poor results. The alternative models that are available in the literature are reviewed, and their merits and limits are highlighted. The modification that could improve the physics of the models implemented in the codes are identified. (author) 13 figs., 123 refs.

Energy Technology Data Exchange (ETDEWEB)

This paper presents an investigation into the feasibility and potential advantages of a small scale Miller cycle natural gas engine for applications such as domestic combined heat and power systems. The Miller cycle engine is compared to a standard Otto cycle engine using cycle analyses and multidimensional simulation, and basic engine design implications are discussed. It is found that the Miller cycle engine has a potential for improved fuel efficiency, but at the cost of a reduced power to weight ratio. A fuel efficiency advantage of 5{yields}10% compared to a standard Otto cycle engine appears possible, however it is stated that further investigations, in particular into the topic of engine friction, are required in order to validate the findings. (author)

Energy Technology Data Exchange (ETDEWEB)

Porosity evolution at reactive interfaces is a key process that governs the evolution and performances of many engineered systems that have important applications in earth and environmental sciences. This is the case, for example, at the interface between cement structures and clays in deep geological nuclear waste disposals. Although in a different transport regime, similar questions arise for permeable reactive barriers used for biogeochemical remediation in surface environments. The COMEDIE project aims at investigating the coupling between transport, hydrodynamics and chemistry when significant variations of porosity occur. The present work focuses on a numerical benchmark used as a design exercise for the future COMEDIE-2D experiment. The use of reactive transport simulation tools like Hytec and Crunch provides predictions of the physico-chemical evolutions that are expected during the future experiments in laboratory. Focus is given in ...

International Nuclear Information System (INIS)

In summary, a scaling analysis of a water-cooled Reactor Cavity Cooling System (RCCS) system was performed based on generic information on the RCCS design of PBMR. The analysis demonstrates that the water-cooled RCCS can be simulated at the ANL NSTF facility at a prototypic scale in the lateral direction and about half scale in the vertical direction. Because, by necessity, the scaling is based on a number of approximations, and because no analytical information is available on the performance of a reference water-cooled RCCS, the scaling analysis presented here needs to be 'validated' by analysis of the steady state and transient performance of a reference water-cooled RCCS design. The analysis of the RCCS performance by CFD and system codes presents a number of challenges including: strong 3-D effects in the cavity and the RCCS tubes; simulation of turbulence in flows characterized by natural circulation, high Rayleigh ...

Energy Technology Data Exchange (ETDEWEB)

An axisymmetric plasma screw pinch is an axisymmetric column of ionized gaseous plasma radially confined by forces from axial and azimuthal currents driven in the plasma and its surroundings. This dissertation is a contribution to detailed, high resolution computer simulation of dynamic plasma screw pinches in 2-d {ital rz}-coordinates. The simulation algorithm combines electron fluid and particle-in-cell (PIC) ion models to represent the plasma in a hybrid fashion. The plasma is assumed to be quasineutral; along with the Darwin approximation to the Maxwell equations, this implies application of Ampere`s law without displacement current. Electron inertia is assumed negligible so that advective terms in the electron momentum equation are ignored. Electrons and ions have separate scalar temperatures, and a scalar plasma electrical resistivity is assumed. Altemating-direction-implicit (ADI) methods are used to advance the electron fluid drift ...

Energy Technology Data Exchange (ETDEWEB)

An experimental and analytical study was performed to improve understanding of the dynamic impact behavior of carbon steel pipes. The test program addressed two types of pipe impact scenarios using both 2- and 4-in. Sch-80 pipes and elbows. Projectile-on-pipe tests simulated the behavior of a stationary target pipe which is impacted at its center by a larger, more rigid whipping pipe. These target pipes, which contained non-flowing water at about 290{degree}C temperature and ca 8.5 megapascals pressure, exhibited a peak deformation of up to 45% reduction in their diameter. For each test condition, the local deformation at the impact zone is a function of the peak impact force and impact velocity. Pipe-on-wall tests simulated the impact of an elbow at the free end of a cantilevered whipping pipe with a rigid barrier. The peak crush of the impact zone of the elbows was a function of peak impact force. Numerical analyses of ...