Energy Technology Data Exchange (ETDEWEB)

A new experimental technique to measure material shear strength at high pressures has been developed for use on magneto-hydrodynamic (MHD) drive pulsed power platforms. By applying an external static magnetic field to the sample region, the MHD drive directly induces a shear stress wave in addition to the usual longitudinal stress wave. Strength is probed by passing this shear wave through a sample material where the transmissible shear stress is limited to the sample strength. The magnitude of the transmitted shear wave is measured via a transverse VISAR system from which the sample strength is determined.

UK PubMed Central (United Kingdom)

The impacts of climate change on crop productivity are often assessed using simulations from a numerical climate model as an input to a crop simulation model. The precision of these predictions reflects...Full Text Available

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 ...

Energy Technology Data Exchange (ETDEWEB)

The present paper is concerned with the prediction method for the motions and tensions of mooring lines on slack-moored floating oceanic architectural building in coastal zone. The motions and the tensions of mooring lines were estimated by the time domain numerical simulation in the present method. The present method includes the nonlinearity of mooring system, slow varying wave drift force. The authors performed a time domain numerical simulation for motions and the tensions of mooring lines on floating barge in two component waves. Also, the results obtained by numerical simulation were compared with the experimental results. The validity of the present method was confirmed.

Science.gov (United States)

... that range from relatively low frequency strong ground motion modeling [Vidale et al., 1985] to high frequency reverse time migration of seismic ...

UK PubMed Central (United Kingdom)

Many simulation studies in biomedicine are based on a similar sequence of processing steps, starting from images and running through geometric model generation, assignment of tissue properties, numerical...Full Text Available

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 validate the numerical ...

International Nuclear Information System (INIS)

In this paper, finite-difference time-domain method is used to research the transient electromagnetic pulse (EMP) coupling to the cylinder object which is shielded in another shielding cavity with a hole numerically. The different coupling effectiveness of three transient EMP is simulated. The distribution of the electric field in the shielding cavity and the current on the cavity are researched. The shielding effectiveness is evaluated. (authors)

Energy Technology Data Exchange (ETDEWEB)

Numerical simulation of dynamic fracture of concrete slabs, impacted by steel projectiles, was carried out in this study. The concrete response was described by a continuum damage model. This continuum damage model was originally developed to study rock fragmentation and was modified in the present study with an emphasis on the post-limit structural response. The model was implemented into a transient dynamic explicit finite element code LS-DYNA2D and the code was then used for the numerical simulations. The specific impact configuration of this study follows the experiment series conducted by Hanchak et al. Comparisons between calculated results and measured data were made. Good agreements were found.

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)

A previous paper introduced the use of wave digital filters as a basic building block for power system simulation, particularly suitable for real-time applications. This paper stresses the simulation of non-linear and switching elements, emphasizing the advantages of the wave filters implementation. The digital structure is maintained even when non-linear components change their characteristics or power electronic devices switch their states. As a very important by-product, the suppression of numerical oscillations related to the trapezoidal rule is achieved in a rather simple way, with no effects on simulation results.

Energy Technology Data Exchange (ETDEWEB)

This paper describes a signal simulator which is designed to calibrate eddy current probes. By using simultaneously an excitation coil and an electrically conducting plane, the simulator can emulate signals which correspond to spot welds with different flaws. The presence of a conducting plane allows a wide variation of the amplitude and phase of the excitation coil voltage during the calibration of a probe. A mathematical model for the simulator is derived and used to study the variation of the impedance change upon the parameters of the probe. Numerical computation shows that the impedance change depends in an important way on the frequency of the excitation coil current.

Energy Technology Data Exchange (ETDEWEB)

RODbaffle heat exchanger is a kind of shell-and-tube heat exchanger with eminent performance. Because of the characteristics of tube bundle support structure, fluid flow on shellside is longitudinal and periodical, and fluid flow is of symmetry in landscape orientation. According to the fluid flow characteristics on shellside and ignoring the impact of inlet, outlet and shell wall, a periodic flow unit duct was taken as the simplified model of the shellside to perform numerical simulation by using CFD code, FLUENT. It was found that the both errors in magnitude of the main fluid velocities and heat transfer coefficients between results of simulation and that of experiment or correlations are less than 10%, and the errors of pressure drop between simulation and correlation are no more than 20%. The detailed characteristics and relation of fluid flow and heat transfer on shellside of the RODbaffle heat ...

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.

Energy Technology Data Exchange (ETDEWEB)

The improvement in the functions of the viscous flow calculation method VEGA-SHIP around a ship and the expansion of application range were described as the numerical flow simulation in ship and ocean engineering and at the same time application examples to the ocean engineering by the general-purpose flow simulation code FLOW-3D handling the non-steady flow with a free surface were introduced as the numerical simulation regarding such products as a water gate and a dam. In the VEGA-SHIP, water surface was handled as a fixed wall so that wave could not be calculated. Therefore, an algorithm for calculating wave on the water surface was added to the VEGA-SHIP and a calculation method simultaneously considering the creation of wave around the ship and viscosity was developed. The FLOW-3D was used to calculate the phenomenon where inside liquid moved greatly due to the oscillation of a ...

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 ...

British Library Electronic Table of Contents (United Kingdom)

This paper presents a numerical study on the formation history of coalbed methane (CBM) reservoir in the southeast edge of Ordos Basin, China. The coal seams studied belong to the Late Palaeozoic coal-bearing series. These coal seams have a burial history and experienced the process of subsidence, rapid subsidence alternated with uplift and then uplift, sequentially, and underwent the geothermal actions at normal, extremely high, and then normal temperatures, respectively. Coal organic matter of the coal seams matured in the Triassic Period and in the Late Jurassic to Early Cretaceous Period. The results from numerical simulation reveal that CBM reservoir evolution history can be classified into five stages, namely primary, initial, stagnant, active and dissipative stages. In the first (pr...

Energy Technology Data Exchange (ETDEWEB)

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)

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...

International Nuclear Information System (INIS)

Full text of publication follows: As the coolant experiences no phase change in the core, SCWRs, unlike LWRs, cannot use design criteria based on the critical heat flux concept. The commonly accepted practice in SCWRs is to specify cladding temperature limits that must be met during transient and accident events. Therefore for the design of the SCWR, it is very important to predict the heat transfer coefficient to the supercritical water coolant with great accuracy. Our recent study focuses on the critical issue of measuring heat transfer to supercritical water at prototypical SCWR conditions and to develop the tools to predict the SCWR thermal behavior. A heat transfer test loop using a surrogate fluids, CO_2, is under construction. The reason of using CO_2 instead of water is that (i) valuable insight of the physical phenomena can be obtained with this fluid, and (ii) some existing facilities already used surrogate fluids, which in general have lower critical pressure and ...

Energy Technology Data Exchange (ETDEWEB)

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.

CERN Document Server

Numerical results on the translocation of long biopolymers through mid-sized and wide pores are presented. The simulations are based on a novel methodology which couples molecular motion to a mesoscopic fluid solvent. Thousands of events of long polymers (up to 8000 monomers) are monitored as they pass through nanopores. Comparison between the different pore sizes shows that wide pores can host a larger number of multiple biopolymer segments, as compared to smaller pores. The simulations provide clear evidence of folding quantization in the translocation process as the biopolymers undertake multi-folded configurations, characterized by a well-defined integer number of folds. Accordingly, the translocation time is no longer represented by a single-exponent power law dependence on the length, as it is the case for single-file translocation through narrow pores. The folding quantization increases with the biopolymer length, ...

Energy Technology Data Exchange (ETDEWEB)

To represent a heterogeneous unsaturated fractured rock by its homogeneous equivalent, Monte Carlo simulations are used to obtain upscaled (effective) flow properties. In this study, we present a numerical procedure for upscaling the van Genuchten parameters of unsaturated fractured rocks by conducting Monte Carlo simulations of the unsaturated flow in a domain under gravity-dominated regime. The simulation domain can be chosen as the scale of block size in the field-scale modeling. The effective conductivity is computed from the steady-state flux at the lower boundary and plotted as a function of the averaging pressure head or saturation over the domain. The scatter plot is then fitted using van Genuchten model and three parameters, i.e., the saturated conductivity K{sub s}, the air-entry parameter {alpha}, the pore-size distribution parameter n, corresponding to this model are considered as the ...

International Nuclear Information System (INIS)

This paper describes the study of particles' dispersion in an isotropic turbulent flow. The particle's motion and the turbulent flow characteristics are calculated independently. While the particles' displacement is computed by the author's code, the flow is simulated with a commercial code : PowerFLOW. The particles and the flow are coupled through the relative velocity component of the aerodynamic force. When the simulated flow is turbulent, a turbulence regeneration model is used in order to get the flow instantaneous velocity. Validation of the method is done by comparing the particles' dispersion obtained with experimental results from literature and with the results calculated by FLUENT. Good accordance is found between numerical studies and experimental results. However, comparison between results of PowerFLOW coupled to the author's code and results from FLUENT shows differences when the particle's path goes through ...

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.

Energy Technology Data Exchange (ETDEWEB)

A new scheme to accurately retain kinetic electron effects in particle-in-cell (PIC) simulations for the case of electrostatic drift waves is presented. The splitting scheme, which is based on exact separation between adiabatic and on adiabatic electron responses, is shown to yield more accurate linear growth rates than the standard df scheme. The linear and nonlinear elliptic problems that arise in the splitting scheme are solved using a multi-grid solver. The multi-grid particle-in-cell approach offers an attractive path, both from the physics and numerical points of view, to simulate kinetic electron dynamics in global toroidal plasmas.

Energy Technology Data Exchange (ETDEWEB)

The results of this study show that Monte Carlo simulations of diffusion in homogeneous turbulence can be formulated in terms of the Eulerian space-time velocity autocorrelation function. Numerical results obtained using this approach agree with results obtained by Taylor (1921) using the Lagrangian autocorrelation function. Use was made of the equivalence of the Lagrangian and Eulerian Monte Carlo approaches to derive analytic relations between the Lagrangian integral time scale and the Eulerian integral space and time scales. These analytic results were verified by comparison to Monte Carlo simulations and to other theoretical results. They are in general agreement with many existing theories and semi-empirical relations.

British Library Electronic Table of Contents (United Kingdom)

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 pr...

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 ...

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 initial spectrum. In ...

Energy Technology Data Exchange (ETDEWEB)

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.

Energy Technology Data Exchange (ETDEWEB)

Stabilized formation temperatures were estimated at different depths in 40 wells from the Los Humeros geothermal field, Mexico, using the Horner and the spherical radial flow (SRF) methods. The results showed that the Horner method underestimates formation temperatures, while the SRF method gives temperatures that are closer to the true formation temperatures. This was supported by numerical simulation of a combined circulation and shut-in period in several wells, and results for well H-26 are presented. Numerical reproduction of logged temperature is more feasible if an initial temperature profile based on the SRF method is employed instead of using an initial temperature profile based on the Horner method. (Author)

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...

Energy Technology Data Exchange (ETDEWEB)

The microbial transport simulator (MTS) is a three-dimensional, three-phase, multiple-component numerical model that permits the study of the transport of microorganisms and nutrients in porous media. Microbial parameters incorporated into MTS include: Microbial growth and decay, microbial deposition, chemotaxis, diffusion, convective dispersion, tumbling, and nutrient consumption. Governing equations for microbial and nutrient transport are coupled with continuity and flow equations under conditions appropriate for a black oil reservoir. The model`s mathematical formulations and preparation procedures of data files for conducting simulations using MTS are described. A general background of microbial transport simulation is given in Section I and the governing equations, mechanisms, and numerical solutions of MTS are given in Section II. Explanations for preparing an input data file ...

Energy Technology Data Exchange (ETDEWEB)

The purpose of this paper is to present an analysis of lightning fault detection, location and protection using numeric distance relays applied in high voltage transmission lines, more specifically in the 500 kV transmission lines of CEMIG (Brazilian Energy Utility) between the Vespasiano 2 - Neves 1 (short line - 23.9 km) and Vespasiano 2 - Mesquita (long line - 148.6 km) substations. The analysis was based on the simulations results of numeric distance protective relays on power transmission lines, realized in September 02 to 06, 2002, at Siemens AG's facilities (Erlangen - Germany), using Real Time Digital Simulator (RTDS{sup TM}). Several lightning faults simulations were accomplished, in several conditions of the electrical power system where the protective relays would be installed. The results are presented not only with the times of lightning faults elimination, but ...

Energy Technology Data Exchange (ETDEWEB)

The Idaho National Engineering Laboratory (INEL) has been conducting geothermal reservoir research and testing sponsored by the US Department of Energy (DOE) since 1983. The INEL research program is primarily aimed at the development of reservoir engineering techniques for fractured geothermal reservoirs. Numerical methods have been developed which allow the simulation of fluid flow and heat transfer in complex fractured reservoirs. Sensitivity studies have illustrated the importance of incorporating the influence of fractures in reservoir simulations. Related efforts include fracture characterization, geochemical reaction kinetics and field testing.

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.

Science.gov (United States)

This dissertation reports on the application of numerical optimization techniques as applied to fuel cell simulation and design. Due to the "multi-physics" inherent in a fuel cell, which results in a highly coupled and non-linear behavior, an experimental program to analyze and improve the performance of fuel cells is extremely difficult. This program applies new optimization techniques with computational methods from the field of aerospace engineering to the fuel cell design problem. After an overview of fuel cell history, importance, and classification, a mathematical model of solid oxide fuel cells (SOFC) is presented. The governing equations are discretized and solved with computational fluid dynamics (CFD) techniques including unstructured meshes, non-linear solution methods, numerical derivatives with complex variables, and sensitivity analysis with adjoint methods. Following the validation of the fuel cell model in ...

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.

Energy Technology Data Exchange (ETDEWEB)

Incomplete convergence in numerical simulation such as computational physics simulations and/or Monte Carlo simulations can enter into the calculation of the objective function in an optimization problem, producing noise, bias, and topo- graphical inaccuracy in the objective function. These affect accuracy and convergence rate in the optimization problem. This paper is concerned with global searching of a diverse parameter space, graduating to accelerated local convergence to a (hopefully) global optimum, in a framework that acknowledges convergence uncertainty and manages model resolu- tion to efficiently reduce uncertainty in the final optimum. In its own right, the global-to-local optimization engine employed here (devised for noise tolerance) performs better than other classical and contemporary optimization approaches tried individually and in combination on the "industrial" test problem to be ...

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)

The aim of this study concerns the use of numerical methods for the resolution of the Reynolds Averaged Navier Stokes equations adapted to the simulation of the cooling of the trailing edge of a stator in a high pressure turbine. These methods, based on the elsA solver developed at ONERA, use a four steps Runge Kutta time discretization scheme and a Jameson centered space discretization scheme. The scheme is applied through a finite volume approach on control volume centered on the cells of a multi-block structured mesh. Turbulence is simulated either through the algebraic Michel model, or through the one-transport-equation Spalart-Allmaras model, or through the two-transport-equations k 1, k {omega} and k {epsilon} models, and through ASM model. A simulation of the flow in a bidimensional stator, without cooling, is carried out. The cooling, which is realized with trailing edge slots, is then ...

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 ...

International Nuclear Information System (INIS)

This thesis was carried out at Forsmark Nuclear Power Plant. The power plant in Forsmark consists of three boiling water reactors (BWR) which produce about 17% of Swedish electricity. In a BWR the nuclear reactions are used to boil water inside the reactor vessel. The water works both as a coolant and as a moderator and the resulting steam is used directly to run the turbines. A problem when running a BWR at low flow conditions is the density wave oscillations that might occur to the water flow inside the fuel assemblies. These oscillations arise due to the connection between power and flow rate in a heated channel with two-phase flow. In order to improve the stability performance of the channel an orifice plate is placed at the inlet of each fuel assembly. Today these orifice plates have sharp edges and a constant resistance coefficient. Experimental work has been done with progressive orifices, the edge of which is half-oval in shape. The advantage of progressive orifices is the ...

CERN Document Server

Numerical simulations predict that a large number of baryons reside in intergalactic space at temperatures between 10^5-10^7 K. Highly-ionized metals, such as O VII and O VIII, are good tracers of this ``warm-hot intergalactic medium'', or WHIM. For collisionally-ionized gas, the ionization fraction of each ion peaks at some particular temperature (``peak temperatures''), so different ions can therefore trace the IGM at different temperatures. We performed a hydrodynamic simulation to study the metal distributions in the IGM. We then draw random lines-of-sight across the simulated region and synthesize resonance absorption line spectra in a similar way to simulating the Ly-alpha forest. By studying the distribution functions of H- and He-like O, Si and Fe in a collisionally-ionized IGM and comparing with semi-analytic results based on the Press-Schechter formalism, we find: (1) ions ...

Energy Technology Data Exchange (ETDEWEB)

A free-electron laser (FEL) two-beam accelerator (TBA) is proposed, in which the FEL interaction takes place in a series of drive cavities, rather than in a waveguide. Each drive cavity is 'beat-coupled' to a section of the accelerating structure. This standing-wave TBA is investigated theoretically and numerically, with analyses included of microwave extraction, growth of the FEL signal through saturation, equilibrium longitudinal beam dynamics following saturation, and sensitivity of the microwave amplitude and phase to errors in current and energy. It is found that phase errors due to current jitter are substantially reduced from previous versions of the TBA. Analytic scalings and numerical simulations are used to obtain an illustrative TBA parameter set.

International Nuclear Information System (INIS)

The design and construction of a Browne-Buechner type round pole face magnetic electron spectrograph is described. The design is based on the first order transfer properties of the bending magnet including the focal surface, dispersion, magnification and resolution. Second order aberrations are compared wit numerical ray tracing results using measured field data to simulate the magnet. The numerical results are also compared with floating wire measurements on the system. Agreement is within one part in 10"4. Data extraction from the instrument is done with an array of copper detectors placed along the focal surface. The spectrograph has been successfully employed in free-electron laser (FEL) experiments for over a year.

CERN Document Server

The effect of the wind on the sustain of extreme water waves is investigated experimentally and numerically. A series of experiments conducted in the Large Air-Sea Interactions Facility (LASIF) showed that a wind blowing over a strongly nonlinear short wave group due to the linear focusing of a modulated wave train may increase the life time of the extreme wave event. The expriments suggested that the air flow separation that occurs on the leeward side of the steep crests may sustain longer the maximum of modulation of the focusing-defocusing cycle. Based on a Boundary-Integral Equation Method and a pressure distribution over the steep crests given by the Jeffreys'sheltering theory, similar numerical simulations have confirmed the experimental results

Energy Technology Data Exchange (ETDEWEB)

Theoretical investigations of fatty acids as a phase change material (PCM) for energy storage system have been conducted in this study. The selected fatty acids were capric acid, lauric acid, myristic acid, palmitic acid and stearic acid. For the two-dimensional simulation model based on the enthalpy approach, calculations have been made for the melt fraction with conduction only. Glass, stainless steel, tin, aluminium mixed, aluminium and copper were used as heat exchanger materials in the numerical calculations. Theoretical results show that capric acid was found good compatibility with latent heat storage system. The large value of thermal conductivity of heat exchanger materials did not make significant contribution on the melt fraction. (author)

CERN Document Server

The abstract is additional with repect to the paper published in JGCD. Ordinary Earth satellites are usually placed into five categories of special orbits: sun-synchronous orbits, orbits at the critical inclination, frozen orbits, repeating ground track orbits, and geostationary orbits. This paper investigates their counterparts around Mars and examines the basic nature of these orbits, which are of special interest for missions conducted around Mars, including Mars reconnaissance. Mars' gravity field is much more complicated, with relatively smaller J2, compared to Earth's, which makes the behaviors of these Martian orbits different from those of Earth. Analytical formulations and numerical simulations are used to analyze these Martian orbits and compare them with their Earth counterparts. First, mean element theory is employed to describe variations of orbital elements and give the constraint conditions for achieving these special orbits. ...

British Library Electronic Table of Contents (United Kingdom)

In energy industries, it is always of difficulty to produce high heat value(HHV) gas continuously using pulverized coal. In this paper, a new type furnace for partitioned alternative gasification using pulverized coal is developed, in which the oxidation and reduction reaction occur alternatively with the dropping of pulverized coal and finally HHV gas could be continuously obtained at the reduction zone exit and low heat value(LHV) gas at the oxidation zone exit. Furthermore, the gasification characteristics and their factors in furnace are numerically simulated under two dimensional model with a self-coded program, based on heat, mass and energy transfer as well as reaction principles. It is found that the producing rate of HHV gas is 1.10Nm3/kg with heat value of 11.72MJ/Nm3, however, t...

British Library Electronic Table of Contents (United Kingdom)

In order to reduce the size and cost of heat exchangers, an air-side wavy fin-and-tube heat transfer surface with three-row tubes needs to be replaced by two-row tubes with some appropriate enhancing techniques. The major purpose of the present paper is to search for such new structure by numerical simulation. First, longitudinal vortex generators of Delta-winglet type are tried. The influence of number and of arrangement of the winglets on the performance of the heat transfer surface is studied in detail. The numerical results show that the fin with two winglets aligned spanwise in the front and rear of each tube (Fin W6) has higher heat transfer capability than other enhanced structures with vortex generators, but it still unable to meet the heat transfer requirement. Then a combination ...

British Library Electronic Table of Contents (United Kingdom)

Biostimulation by nutrient application is a viable technology for restoring oil-contaminated beaches. Maximizing the nutrient residence time is key for achieving a rapid cost-effective cleanup. We considered the nutrient injection strategy through a perforated pipe at the high tide line and we simulated numerically beach hydraulics, which allowed us to estimate the optimal injection flow rate of nutrient solution. Our results indicate that the optimal application is one that starts following the falling high tide and lasts for half tidal cycle. The saturated wet-front of the nutrient solution on the beach surface would move seaward with the same speed of the falling tide keeping a constant distance with the tide line. The numerical results were generalized to beaches of wide ranges of hydr...

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)

Turbulent heat transfer performance of a fuel rod with three-dimensional trapezoidal spacer ribs for high temperature gas-cooled reactors was studied for various Reynolds numbers using an annular channel at the same coolant condition as the reactor operation, maximum outlet temperature of 1000 C and pressure of 4 MPa, and analytically by a numerical simulation using the k-#epsilon# turbulence model. The turbulent heat transfer coefficients of the fuel rod were 18-80% higher than those of a concentric smooth annulus at a region of Reynolds number exceeding 2000. On the other hand, the predicted average Nusselt number of the fuel rod agreed well with the empirical correlation obtained from the experimental data within a relative error of 10% with Reynolds number of more than 5000. It was verified that the numerical analysis results had sufficient accuracy. Furthermore, the numerical prediction could ...

Energy Technology Data Exchange (ETDEWEB)

Discussed in this report is a wavefield simulation in the 3-dimensional seismic survey. With the level of the object of exploration growing deeper and the object more complicated in structure, the survey method is now turning 3-dimensional. There are several modelling methods for numerical calculation of 3-dimensional wavefields, such as the difference method, pseudospectral method, and the like, all of which demand an exorbitantly large memory and long calculation time, and are costly. Such methods have of late become feasible, however, thanks to the advent of the parallel computer. As compared with the difference method, the pseudospectral method requires a smaller computer memory and shorter computation time, and is more flexible in accepting models. It outputs the result in fullwave just like the difference method, and does not cause wavefield numerical variance. As the computation platform, the parallel computer ...

Energy Technology Data Exchange (ETDEWEB)

The scalar mixing time scale, a key quantity in many turbulent combustion models, is investigated for reactive scalars in premixed combustion. Direct numerical simulations (DNS) of three-dimensional, turbulent Bunsen flames with reduced methane-air chemistry have been analyzed in the thin reaction zones regime. Previous conclusions from single step chemistry DNS studies are confirmed regarding the role of dilatation and turbulence-chemistry interactions on the progress variable dissipation rate. Compared to the progress variable, the mixing rates of intermediate species is found to be several times greater. The variation of species mixing rates are explained with reference to the structure of one-dimensional premixed laminar flames. According to this analysis, mixing rates are governed by the strong gradients which are imposed by flamelet structures at high Damkoehler numbers. This suggests a modeling approach to estimate the mixing rate of ...

Science.gov (United States)

The behavior of a small power system consisting of two interconnected generators is simulated in real-time by a prototype Digital Transient Network Analyzer (TNA). The prototype Digital TNA consists of two Computational Modules and one I/O Module. The Modules communicate with each other through ribbon cables. Each Computational Module simulates one Turbo-Generator, its Transformer, its Governor, Exciter, and Power System Stabilizer Systems. The numerical integration is shared by two TMS320C30 DSPs at a step-size of 100 microseconds in real-time. The I/O module post-processes the state variables and presents selected information for analog display. The paper presents oscillograms from a test program which includes symmetry checks and behavioral checks against well known waveforms of hunting oscillations, synchronization out-of-phase torques, and subsynchronous resonance phenomena. The success of the Digital TNA depends on: ...

Energy Technology Data Exchange (ETDEWEB)

The results of an experimental investigation of the performance of a small bore engine with a swirl chamber when operating as a dual fuel engine with commercial methane as the gaseous fuel were presented in this paper. The experiment involved using a 3-dimensional computational fluid dynamics model to predict the performance of the engine. A detailed chemical kinetics for the gaseous fuel component, consisting primarily of methane and a reduced detailed chemical kinetics for the diesel fuel while considering the turbulent combustion processes an associated performance of a dual fuel engine with a swirl chamber were incorporated in the simulation. The study experimentally and numerically investigated the effects of changes in the quantities of the liquid fuel pilot and gaseous fuels on the combustion processes, engine performance, cyclic variations, and emissions. The paper discussed the experimental approach and results. It also discussed the ...

CERN Document Server

This work investigates the impingement of a liquid microdroplet onto a glass substrate at different temperatures. A finite-element model is applied to simulate the transient fluid dynamics and heat transfer during the process. Results for impingement under both isothermal and non-isothermal conditions are presented for four liquids: isopropanol, water, dielectric fluid (FC-72) and eutectic tin-lead solder (63Sn-37Pb). The objective of the work is to select liquids for a combined numerical and experimental study involving a high resolution, laser-based interfacial temperature measurement to measure interfacial heat transfer during microdroplet deposition. Applications include spray cooling, micro-manufacturing and coating processes, and electronics packaging. The initial droplet diameter and impact velocity are 80 {\\mu}m and 5 m/s, respectively. For isothermal impact, our simulations with water and isopropanol show very ...

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)

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, ...

International Nuclear Information System (INIS)

In this paper, a time-domain equivalent circuit method is applied to solve dispersion of coupled-cavity travelling-wave tube (CCTWT). First, the time-domain circuit equations of CCTWT coupled-cavity chain are deduced from the equivalent circuit model. Then, the equations are solved numerically by fourth-order Runge-Kutta method and a program CTTDCP is developed using MATLAB. Last, a L-band CCTWT is calculated using CTTDCP and the cavity pass-band of this tube is computed to be 1.08-1.48 GHz, which is consistent with the experimental results and the simulation results of electromagnetic code and demonstrates the validity of the time-domain equivalent circuit method. In addition, a new design method which uses the equivalent circuit method and electromagnetic simulation together to optimize the cold cavity characteristics of CCTWT is proposed. (authors)

Energy Technology Data Exchange (ETDEWEB)

Recently, we showed (Yu, Phys. Rev. E 58 (1998) 4991) 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, we apply 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, we show 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 subtracted 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)

In this paper linear and nonlinear models of spacecraft attitude dynamics equations and gravity gradient moments are investigated. In addition, effects of gravity gradient moments on attitude dynamics of the satellite are studied. The purpose of this paper is to present a comparison between nonlinear and linear models of spacecraft attitude dynamics and gravity gradient moments in order to determine divergence of linear approximation from the nonlinear model. Simulation results indicate that designer of spacecraft attitude control subsystem should be meticulous in applying linear approximation of equations especially in low earth orbits. Consequently, finding an upper bound for small angle to keep the linear model valid and precise enough would be a vital part of using linear approximation...

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)

There has been increased theoretical and experimental research interest in autonomous mobile robots exhibiting cooperative behaviour. This paper provides consistent quantitative measures of organizational degree of a two-dimensional environment. We proved, by the way of numerical simulations, that the theoretically derived values of the feature are reliable measures of aggregation degree. The slope of the feature's dependence on memory radius leads to an optimization criterion for stochastic functional self-organization. We also described the intellectual heritages that have guided our research, as well as possible future developments. (author)

International Nuclear Information System (INIS)

Axial offset anomaly (AOA) is caused by the deposition of crud on the fuel cladding of a PWR. When significant levels of crud build up on the cladding, boron can accumulate in the pores of the crud as a concentrated solution or solid phase, and cause the flux depression. Numerous studies have been conducted on the primary water chemistry to reduce the amount of crud in the primary circuit to avoid radioactivity buildup and unexpected power transition in the plant. However, experiments on the crud are restricted in the laboratory because the crud is a highly radioactive material. The objective of this study is to develop a test method for simulating the deposition of crud in a nuclear power plant

Energy Technology Data Exchange (ETDEWEB)

This report details an investigation into the efficacy of two approaches to solving the radiation diffusion equation within a radiation hydrodynamic simulation. Because leading-edge scientific computing platforms have evolved from large single-node vector processors to parallel aggregates containing tens to thousands of individual CPU's, the ability of an algorithm to maintain high compute efficiency when distributed over a large array of nodes is critically important. The viability of an algorithm thus hinges upon the tripartite question of numerical accuracy, total time to solution, and parallel efficiency.

International Nuclear Information System (INIS)

Hung and Quyen's model (analysis of strategic interactions between players in the game theoretic framework) is first recapitulated. Solutions to the embargo game with the seller acting as the Stackelberg leader in both periods. Then the timing coordination issue is discussed and perfect equilibrium outcomes under different structures of leader-follower in the game are compared. Numerical simulations show that the structure corresponding to alternated leadership from one period to the other yields the perfect equilibrium outcome that is Pareto improving with respect to the structure where the buyer is first mover in both periods and Pareto dominant for some specific values of the parameters embedded in the game.

DEFF Research Database (Denmark)

Nonlinear distortion added by the loudspeaker in a hearing aid lowers the signal-to-noise ratio and may degrade the hearing aid user's ability to understand speech. The balancedarmature- type loudspeakers, predominantly used in hearing aids, are inherently nonlinear devices, as any displacement of the loudspeaker diaphragm inevitably changes the magnetic and electrical characteristics of the loudspeaker. A numerical time-domain model capable of describing these nonlinearities is presented. By simulation it is demonstrated how the output distortion could potentially be reduced significantly through careful design of the mechanical properties of the armature.

Energy Technology Data Exchange (ETDEWEB)

A numerical model has been developed to simulate the charging and discharge behaviour of an electrically rechargeable alkaline zinc-air battery. Further a galvanostatic experiment including three charge/discharge cycles has been performed. The cell voltages, the Zn electrode potentials versus a Zn reference, and the O{sub 2} electrode potentials versus a Zn reference calculated with the model are in fairly good agreement with the corresponding experimental data. The model is expected to be useful for zinc-air battery design and for analysis of experimental data. (author)

Science.gov (United States)

The electronic transitions and photodissociation of the bromine molecule were studied in the visible-near UV continuum using dynamic simulation. The molar extinction coefficients in this study were obtained in numerical calculations. The quantum yields of the spin-orbit Br*(2 P 1/2) product at different photon frequencies were determined. Time-dependent density functional theory was used to analyze the highest five occupied and lowest five unoccupied Br2 orbitals. The transition to the 1? u state was found to be most probable in the visible-near UV absorption range.

Energy Technology Data Exchange (ETDEWEB)

The primary objective is to perform a fundamental study of filtration of nanoparticles, and to obtain filtration knowledge necessary to design particle collection devices/systems for nanoparticle processing and for preventing nanoparticle emissions into the environment. The research covered a wide area relevant to nanoparticle filtration, under these main topics: (1) nanoparticle filtration and molecular dynamics simulation, (2) nanoparticle virtual impactor, (3) particle transport under low pressure, and (4) development of a high-throughput nanoparticle generator. A number of novel tools and numerical models have been developed under the DOE support.

International Nuclear Information System (INIS)

Environmental computer science is a new partial discipline of applied computer science, which makes use of methods and techniques of information processing in environmental protection. Thanks to the inter-disciplinary nature of environmental problems, computer science acts as a mediator between numerous disciplines and institutions in this sector. The handbook reflects the broad spectrum of state-of-the art environmental computer science. The following important subjects are dealt with: Environmental databases and information systems, environmental monitoring, modelling and simulation, visualization of environmental data and knowledge-based systems in the environmental sector. (orig.).

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.

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)

Longitudinal space charge (LSC) force can be a main effect driving the microbunching instability in the linac for an x-ray free-electron laser (FEL). In this paper, the LSC-induced beam modulation is studied using an integral equation approach that takes into account the transverse (radial) variation of LSC field for both the coasting beam limit and bunched beam. Changes of beam energy and the transverse beam size can be also incorporated. We discuss the validity of this approach and compare it with other analyses as well as numerical simulations.

CERN Document Server

An algorithm for pose and motion estimation using corresponding features in images and a digital terrain map is proposed. Using a Digital Terrain (or Digital Elevation) Map (DTM/DEM) as a global reference enables recovering the absolute position and orientation of the camera. In order to do this, the DTM is used to formulate a constraint between corresponding features in two consecutive frames. The utilization of data is shown to improve the robustness and accuracy of the inertial navigation algorithm. Extended Kalman filter was used to combine results of inertial navigation algorithm and proposed vision-based navigation algorithm. The feasibility of this algorithms is established through numerical simulations.

Energy Technology Data Exchange (ETDEWEB)

This paper deals with the efficient simulation of the dynamical behaviour of molten carbonate fuel cells (MCFCs). MCFCs allow an efficient and environmentally friendly energy production via electrochemical reactions. Their dynamics can be described by large scale systems of up to currently 22 nonlinear partial differential algebraic equations (PDAE). The paper also serves as a basis for later parameter identification and optimal control purposes. Therefore, the numerical simulations are particularly based on hierarchically embedded systems of PDAE, first of all in one space dimension. The PDAE are of mixed parabolic-hyperbolic type and are completed by nonlinear initial and boundary conditions of mixed type. For a series of embedded models in one space dimension, the vertical method of lines (MOL) is used throughout this paper. For the semi-discretization in space appropriate difference schemes are applied depending on the ...

Energy Technology Data Exchange (ETDEWEB)

This paper is dedicated to the numerical study of natural ventilation in a room through a large external opening using the CFD code Fluent. Bidimensional numerical simulations are performed for wind speeds up to 5.55 m-s{sup -1} (20 km-h{sup -1}). We propose a two-step solution procedure with grid refinement. When combined with appropriate thermal boundary conditions, this technique appears very efficient at limiting local convergence problems. Through an analysis of the flow pattern in the cavity, we qualitatively explain those results and define a critical Archimedes number. Calculated air change rates are compared to values deduced from empirical correlations of the literature. Agreement is in general poor, which can be explained by the differences between the conditions of our simulations and those of the experiments that lead to those correlations. However, from our ...

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)

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 this paper, molecular quantum computation is numerically studied with the quantum search algorithm (Grover's algorithm) by means of optimal control simulation. Qubits are implemented in the vibronic states of I_2, while gate operations are realized by optimally designed laser pulses. The methodological aspects of the simulation are discussed in detail. We show that the algorithm for solving a gate pulse-design problem has the same mathematical form as a state-to-state control problem in the density matrix formalism, which provides monotonically convergent algorithms as an alternative to the Krotov method. The sequential irradiation of separately designed gate pulses leads to the population distribution predicted by Grover's algorithm. The computational accuracy is reduced by the imperfect quality of the pulse design and by the electronic decoherence processes that are modeled by the non-Markovian master equation. ...

Energy Technology Data Exchange (ETDEWEB)

The generation of a defined swivel momentum in car door hinges depends on numerous constructional and technical manufacturing parameters. If these parameters and their influence are to be investigated, then in addition to detailed experiments with variations in the parameters, methods are also required which enable the measuring data produced to be assessed in such a way that, in general, the non-linear relationships between initial and target size can be described sufficiently accurately. This paper explains the parameter reduction necessary in the experimental investigation, gives the results of the data assessment with conventional statistical methods and describes in particular the use of artificial neural networks (ANN) to construct so-called 'neuro hinge models' on the basis of the data resulting from the experiments. Parameter variations can be simulated with the hinge models and in this way optimal constructional and ...

Energy Technology Data Exchange (ETDEWEB)

To significantly improve the simulation of climate by general circulation models (GCMs), systematic errors in representations of relevant processes must first be identified, and then reduced. This endeavor demands that the GCM parameterizations of unresolved processes, in particular, should be tested over a wide range of time scales, not just in climate simulations. Thus, a numerical weather prediction (NWP) methodology for evaluating model parameterizations and gaining insights into their behavior may prove useful, provided that suitable adaptations are made for implementation in climate GCMs. This method entails the generation of short-range weather forecasts by a realistically initialized climate GCM, and the application of six-hourly NWP analyses and observations of parameterized variables to evaluate these forecasts. The behavior of the parameterizations in such a weather-forecasting framework can provide insights on ...

Energy Technology Data Exchange (ETDEWEB)

In technical combustors, flow and combustion processes inside the combustion chamber determine the design and also the safety of the plant. A modelling technique was developed at TU Dresden which comprises both numerical simulations and experimental data and thus achieves a high degree of transparency in its visualisation of the processes inside the combustion chambers under investigation. With the experience gained, numerical models of commercial simulation software (FLUENT, CFX) can be validated by measured values, and potential applications of the programs can be investigated. The first stage of experimental modelling consisted in isothermal (cold) modelling in consideration of flow-relevant similarity characteristics. In the reactive, non-isothermal models, the fuel/air ratio and temperature conditions inside the original combustion chamber were modelled as well. In order to obtain detailed data on ...

Energy Technology Data Exchange (ETDEWEB)

A three-dimensional, three-phase, multiple-component numerical simulator was developed to investigate transport and growth of microorganisms in porous media and the impacts of microbial activities on oil recovery. The microbial activities modeled in this study included: (1) growth, retention, chemotaxis, and end product inhibition of growth, (2) the formation of metabolic products, and (3) the consumption of nutrients. Major mechanisms for microbial enhanced oil recovery (MEOR) processes were modeled as follows: (1) improvement in sweep efficiency of a displacement process due to in situ plugging of highly-permeable production zones by cell mass or due to improved mobility control achieved by increasing the viscosity of the displacing fluid with a biopolymer, and (2) solubilization and mobilization of residual oil in porous media due to the reduction of the interfacial tension between oleic and aqueous phases by the production of a ...

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 ILP-induced modifications of ...

Energy Technology Data Exchange (ETDEWEB)

A numerical study has been performed to investigate the hydrodynamic aspects of the pool boiling on horizontal-, vertical- and downward-facing surfaces. The FlowLab code, which is based on a Lattice-Boltzmann (LB) model of two-phase flows, is employed. Macroscopic properties, such as surface tension ({sigma}) and contact angle ({beta}), are implemented through the fluid-fluid (G{sub {sigma}}) and fluid-solid (G{sub t}) interaction potentials. The model is found to express a linear relation between the macroscopic properties ({sigma}, {beta}) and microscopic parameters (G{sub {sigma}}, G{sub t}). The simulation results on bubble departure diameter appear to have the same parametric dependence as the empirical correlation. Hydrodynamic aspects of two-phase flow regime transition mechanism are investigated for different surface-coolant configurations. Results of the LB simulation clearly demonstrate that not only the bubble ...

Energy Technology Data Exchange (ETDEWEB)

Simulation and experimental results on the transport of microbes and nutrients in one-dimensional cores are presented, and the development of a three-dimensional, three-phase, multiple-component numerical model to describe the microbial transport and oil recovery in porous media is described. The change of rock`s wettability and associated relative permeability values after microbial treatments were accounted for in the model for additional oil recovery. Porosity and permeability reductions due to cell clogging have been considered and the production of gas by microbial metabolism has been incorporated. Governing equations for microbial and nutrient transport are coupled with continuity and flow equations under conditions appropriate for a black oil reservoir. The computer simulator has been used to determine the effects of various transport parameters on microbial transport phenomena. The model can accurately describe the ...

CERN Document Server

In the present study, applications of the SPH method to industrial related issues are considered by starting from an existing open source 2D SPH code, namely the SPHYSICS code, which offers an effective ground for numerical developments, which are performed in order to bring an answer to industrial problems, such as simulations of solid/fluid coupling in a free surface flow context. The purpose of the present paper is therefore to expose the numerical developments which yield an enhanced version (referred to as "SPHYSIC2") of the initial code. Firstly, the different features added to obtain the operational code needed for engineering applications are described, and so are the problems raised on this way, offering a kind of review of SPH methods for engineers. Secondly, the validation of the proposed code is partially presented with two well known but difficult test cases, namely the classical "dam break" and "wedge entry" ...

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 the air flow rate, ...

Energy Technology Data Exchange (ETDEWEB)

Increase usage of magnesium components in the portable electronics industries has propelled considerable research interest on the development of magnesium forging. There have been some studies done in the area of sheet metal forming and forging of axis-symmetrical parts of magnesium alloys but there are very little work done on the plane strain backward extrusion of magnesium alloy. This paper deals with the analysis of plane strain extrusion of magnesium alloy, AZ31, using numerical modeling at different extrusion conditions of temperatures and friction coefficients. A commercially available finite element package, ANSYS/LSDYNA, was used in the finite element analysis. FE simulated results, material deformation, punch force and stress/strain distributions, were compared with actual experimental data. From the numerical analysis, a decrease in forging temperature from 300 C to 250 C will lead to a moderate increase in ...

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 ...

Energy Technology Data Exchange (ETDEWEB)

Moving out of the scientific community research laboratories, computational fluid dynamics (CFD) software packages are now allowing industrials to analyse and optimize industrial processes involving the use of gases, liquids and even some two-phase fluids. Their attractiveness and their impact stems out from the opportunity they offer to bring insight into an existing unit, or even at the design stage, by displaying the spatial distribution of process relevant variables such as temperature, concentration. The filling of the spacing in between a two-layer window is a simple example. This new opportunity of visualisation is at times an unique way, when the process environment is an opaque one, such as liquid metal flowing into a tundish or when measurements of flows may be a long and tedious work, such as flows within water treatment basins. This environment we are to investigate in order to optimize can also be a harsh one, due to its high temperature level for example. Such are ...

Energy Technology Data Exchange (ETDEWEB)

Excessive erosion of the labyrinth seal of a 100 MW geothermal turbine has been investigated. This study used computational fluid dynamics (CFD) and aims to identify one cause of erosion and a possible solution for substantially reducing it. The predictions were based upon a numerical calculation using a CFD model of the labyrinth seal with a water/steam flow containing hard solid particles and solved with a commercial CFD code: Fluent V5.0. The results confirmed the existence of flow conditions that play a major role in the rotor labyrinth seal erosion. Afterwards, the flow path was simulated with changes of rotor labyrinth seal geometry, which are indeed feasible of being implemented. The results confirmed that it is possible to reduce the erosion process by approximately 80% by incorporating a steam flow deflector in the fourth stage diaphragm, which changes the steam flow direction in the inlet zone to the rotor labyrinth seal channel, ...

Energy Technology Data Exchange (ETDEWEB)

This is a numerical verification of two groups of pipe whip tests sponsored or cosponsored by EPRI. Experimental data of the two pipe whip tests, one by Tennessee Valley Authority (TVA) and by FRAMATOME/CEA, were provided by EPRI. A nonlinear finite element code, ABAQUS-EPGEN, developed under partial sponsorship by EPRI was used for modeling the pipe whip tests. Beam elements together with an equivalent nonlinear spring element or a partial shell mesh were used to model pipes and elbow in the pipe whip tests. Material nonlinearity due to plasticity, strain rate effects, and temperature, as well as geometric nonlinearity due to large rotation and boundary conditions were included in the study. Effects of strain rate and modeling techniques were assessed. Results by current industry approach were also included as a reference solution. This report can be used as a guideline for numerical simulation of pipe whip phenomena. 74 ...

Energy Technology Data Exchange (ETDEWEB)

This paper describes the comprehensive modeling method of fluid and heat flows in the hydrological system. Based upon the concept that the hydrological system is composed of the air-water 2-phase fluids, and the rock phase, various natural processes are modeled including surface/subsurface 2-phase mass transfer, heat exchanges between fluid phases, between fluid and solid phases, and sensible/latent heat exchanges on the ground surface. A field-oriented numerical simulator is developed, in which a set of governing equations is solved for different variables on the surface and in subsurface regions, respectively. Results of an experimental study are presented, in which transient formulation of heat exchange between fluid/solid phases is examined through matching of the observation and calculated performances. (author)

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 ...

International Nuclear Information System (INIS)

DF (decontamination factor) of two types of air filter was investigated against "1"9"8Au aerosol particles with AMAD (activity median aerodynamic diameter) = 1.2 #mu#m and #sigma#g (geometric standard deviation) = 2.0. One of the filters, SO filter, is HEPA (high efficiency particulate air) filter media with 46 mm#phi#, and the other is HE-40T (cellulose-glass fiber filter for dust sampling) filter media with 40 mm#phi#. DF of the SO filter was found to be over 5 x 10"4 at air flow rate of 1 l/min. And DF of HE-40 T filter was determined to be 94#+-#4 at the same airflow rate. In the multistage HE-40 T filter system, DF decreased at each succeeding stage. It was confirmed by numerical simulation that the decrease of DF was not inconsistent with the decrease of collection performance in each successive filter. It was also confirmed that DF could be predicted by the numerical simulaiton, if particle penetration function of ...

Energy Technology Data Exchange (ETDEWEB)

The convection heat loss from cavity receiver in parabolic dish solar thermal power system can significantly reduce the efficiency and consequently the cost effectiveness of the system. It is important to assess this heat loss and subsequently improve the thermal performance of the receiver. This paper aims to present a comprehensive review and systematic summarization of the state of the art in the research and progress in this area. The efforts include the convection heat loss mechanism, experimental and numerical investigations on the cavity receivers with varied shapes that have been considered up to date, and the Nusselt number correlations developed for convection heat loss prediction as well as the wind effect. One of the most important features of this paper is that it has covered numerous cavity literatures encountered in various other engineering systems, such as those in electronic cooling devices and buildings. The studies related ...

Science.gov (United States)

The main purpose of this paper is to explore a numerical algorithm for determining the contact stress when a circular crowned roller is compressed between two plates. To start with, the deformation curve on a plate surface will be derived by using the contact mechanical model. Then, the contact stress distribution along the roller which occurs on the plate surface is divided into three parts: from the center of contact to the edge, the edge and apart from the contact line. The first part is calculated by the elastic contact theorem for the contact subjected to nominal stress between non-crowned parts of roller and plates, the second part is obtained by the classical Hertzian contact solution for the contact between crowned parts of roller and plates, and the third part is simulated as exponential decay. In order to overcome the defect of the half space theorem, in which a plate with infinite thickness is assumed initially, a weighting method is ...

International Nuclear Information System (INIS)

It is well known that one of the main advantages of the high speed forming (HSF) processes is the improvement in the forming limits of the used materials.Using the Electromagnetic Forming (EMF) technology two materials have been tested with different mechanical and physical properties: the AA5754 aluminium and the AZ31B magnesium alloys.The EMF process principle can be described as follows: A significant amount of electrical energy is stored in a bank of capacitors which are suddenly discharged releasing all the stored energy. This electric discharge runs through a coil which generates an intense transient magnetic field. At the same time transient Eddy currents are induced in the electrically conductive part placed some millimetres far from the coil. Another intense magnetic field is generated due to those Eddy currents but on the opposite direction as the one generated by the coil. A big magnetic repulsion force is created between the part and the coil. This magnetic repulsion ...

Energy Technology Data Exchange (ETDEWEB)

This report provides a user's manual for PROTOCOL, a comprehensive coupled kinetic/equilibrium computer program for analyzing the dissolution reactions of solids with aqueous solutions, specifically applied to the potential corrosion of vitrified nuclear waste by groundwater. The capabilities and available options are summarized as well as instructions for setting up and running problems. Also described in this report and included in the PROTOCOL software package are MASTER, a master file of species thermodynamic data, MANEQL, a preprocessor program and POSTP, a postprocessor. POSTP provides offline plotting using the CRAY-1 DISSPLA 9.0 graphics library. PROTOCOL is operational on the CDC-7600 and CRAY-1 computers at the Lawrence Livermore National Laboratory. 7 references, 10 figures, 2 tables.

Science.gov (United States)

Self-amplifying spontaneous radiation free-electron lasers, such as the LCLS or the European X-FEL, rely on the incoherent, spontaneous radiation as the seed for the amplifying process. Though this method overcomes the need for an external seed source one drawback is the incoherence of the effective seed signal. The FEL process allows for a natural growth of the coherence because the radiation phase information is spread out within the bunch due to slippage and diffraction of the radiation field. However, at short wavelengths this spreading is not sufficient to achieve complete coherence. In this presentation we report on the results of numerical simulations of the LCLS X-ray FEL. From the obtained radiation field distribution the coherence properties are extracted to help to characterize the FEL as a light source.

British Library Electronic Table of Contents (United Kingdom)

An investigation into a three-dimensional, curved shock wave interacting with a three-dimensional, curved boundary layer on a slender body is presented. Three different nose profiles mounted on a cylindrical body were tested in a supersonic wind tunnel and numerically simulated by solving the Navier?Stokes equations. The conical and hemispherical nose profiles tested were found to generate shock waves of sufficient strength to separate the boundary layer on the cylinder, while the shock wave generated by the ogival profile did not separate the boundary layer. For the separated flow, separation was found to occur predominantly on the windward side of the cylinder with the lee-side remaining shielded from the direct impact of the incident shock wave. A thickening of the boundary layer on the...

CERN Document Server

We formulate a complete theory of Edge Radiation based on a novel method relying on Fourier Optics techniques. Similar types of radiation like Transition Undulator Radiation are addressed in the framework of the same formalism. Special attention is payed in discussing the validity of approximations upon which the theory is built. Our study makes consistent use of both similarity techniques and comparisons with numerical results from simulation. We discuss both near and far zone. Physical understanding of many asymptotes is discussed. Based on the solution of the field equation with a tensor Green's function technique, we also discuss an analytical model to describe the presence of a vacuum chamber. In particular, explicit calculations for a circular vacuum chamber are reported. Finally, we consider the use of Edge Radiation as a tool for electron beam diagnostics. We discuss Coherent Edge Radiation, Extraction of Edge Radiation by a mirror, and ...

British Library Electronic Table of Contents (United Kingdom)

We present structure calculations of neutral and singly ionized Mg clusters of up to 30 atoms, as well as Na clusters of up to 10 atoms. The calculations have been performed using density functional theory (DFT) within the local (spin-)density approximation, ion cores are described by pseudopotentials. We have utilized a new algorithm for solving the Kohn-Sham equations that is formulated entirely in coordinate space and, thus, permits straightforward control of the spatial resolution. Our numerical method is particularly suitable for modern parallel computer architectures; we have thus been able to combine an unrestricted simulated annealing procedure with electronic structure calculations of high spatial resolution, corresponding to a plane-wave cutoff of 954 eV for Mg. We report the geo...

Energy Technology Data Exchange (ETDEWEB)

Polyethylene (UHMWPE) wear in current knee prosthesis causes prosthesis loosening after no more than 15 years. In this work, a steady state one-dimensional lubrication model with non- Newtonian fluid, porous elastic layer on tibial component, ultra-filtration mechanism of fluid and some features of the surface roughness is studied through a numerical technique based on the Finite Element Method. The results show that the UHMWPE stiffness makes difficult the lubrication mechanism of the artificial joint and promotes abrasive and fatigue wear. Nevertheless, the use of compliant porous materials on the tibial component could reduce friction and wear. Moreover, the ultra-filtration mechanism promotes efficiency on the joint.

Science.gov (United States)

We present a numerical model of vapor-dominated reservoirs which is based on the well-known conceptual model of White, Muffler, and Truesdell. Computer simulations show that upon heat recharge at the base, a single phase liquid-dominated geothermal reservoir in fractured rock with low matrix permeability will evolve into a two-phase reservoir with B.P.D. (boiling point-for-depth) pressure and temperature profiles. A rather limited discharge event through cracks in the caprock, involving loss of only a few percent of fluids in place, is sufficient to set the system off to evolve a vapor-dominated state. The attributes of this state are discussed, and some features requiring further clarification are identified. 26 refs., 5 figs.

International Nuclear Information System (INIS)

We present the design, analysis, and results of the high brightness electron beam experiments currently under investigation at Sandia National Laboratories. The anticipated beam parameters are the following: energy 12 MeV, current 35 endash 40 kA, rms radius 0.5 mm, and pulse duration 40 ns full width at half-maximum. The accelerator is SABRE, a pulsed linear inductive voltage adder modified to higher impedance, and the electron source is a magnetically immersed foilless electron diode. 20 endash 30 T solenoidal magnets are required to insulate the diode and contain the beam to its extremely small-sized (1 mm) envelope. These experiments are designed to push the technology to produce the highest possible electron current in a submillimeter radius beam. Design, numerical simulations, and experimental results are presented. copyright 1997 American Institute of Physics.

Energy Technology Data Exchange (ETDEWEB)

Water vapor is well known to be a critical component in many aspects of atmospheric research, such as radiative transfer and cloud and aerosol processes. This requires both improved measurements of the columnar water vapor and its profiles in the atmosphere in a wide range of conditions, and adjustment of water vapor parameterizations in radiation codes including the perfection of spectroscopic parameters. In this paper we will present the results of comparison of our calculations and downward solar fluxes measured with Rotating Shadowband Spectroradiometer under conditions of horizontally homogeneous clouds. We also will discuss the sensitivity of atmospheric radiation characteristics to variations of water vapor in the band 940 nm: these results may be useful for development of new methods of retrieval of the total column water vapor content (WVC) in the atmosphere from data of radiation observations.

British Library Electronic Table of Contents (United Kingdom)

The present investigation focuses on the study of design parameters for a generic fan-in-wing configuration. The objective is to evaluate the performance of a fan-in-wing concept in a Short-Take-Off and Landing situation. The results presented are also of interest for transition between hover and a wing-sustained flight for a Vertical Take Off and Landing concept. Steady Reynolds Averaged Navier-Stokes computations using an actuator disk technique are performed. The results show a good agreement with experimental data and simulation results obtained in previous investigations for a generic wind tunnel model. A parametric study is conducted for a more realistic wing planform. For a single fan installed in the wing, the streamwise and spanwise influence on the wing aerodynamics is studied. I...

Energy Technology Data Exchange (ETDEWEB)

This paper describes a nonlinear dynamic analysis of TVA high energy line pipe whip tests using the ABAQUS-EPGEN code. The analysis considers the effects of large deformation and strain rate on resisting moment and energy absorption capability. The numerical results of impact forces, impact velocities, pipe strains, and reaction forces at pipe supports are compared to the TVA test data. The calculated pipe whip impact time and forces are also compared with those predicted using current industry practice. The calculated pipe support reaction forces are found to be in good agreement with the TVA test data except for some peak values at the very beginning of the pipe break. These peaks are believed to be due to stress wave propagation which cannot be addressed by the ABAQUS code. Both elbow crushing and strain rate have been approximately simulated. The effects are found to be important for pipe whip impact evaluation.

International Nuclear Information System (INIS)

The recognition that natural convection offers the prospect of an important inherent safety feature for liquid metal cooled reactor systems has provided the impetus for a world-wide research effort over the past decade. Whilst this research has been based on experiment, both plant experiments and out-of-pile experiments, the enormous advances in the development of computing power in recent years have enabled complementary programmes of mathematical modelling through numerical simulation of the transport equations in three spatial dimensions. These not only offer considerable promise for the designer in projecting the behaviour from experiments and prototype plant to full scale plant, they have also proved to be of considerable value in helping us to interpret and understand the results of the experiments themselves. This paper attempts to review the progress made with the emphasis on decay heat removal by natural convection in the pool-type ...

CERN Document Server

The finite element method and dimensional analysis have been applied in the present paper to study a hydraulic impact, which is utilized in a non-explosive rock breaking technology in mining industry. The impact process of a high speed piston on liquid water, previously introduced in a borehole drilled in rock, is numerically simulated. The research is focused on the influences of all the parameters involved in the technology on the largest principal stress in the rock, which is considered as one of the key factors to break the rock. Our detailed parametric investigation reveals that the variation of the isotropic rock material properties, especially its density, has no significant influence on the largest principal stress. The influences of the depth of the hole and the depth of the water column are also very small. On the other hand, increasing the initial kinetic energy of the piston can dramatically increase the largest principal stress and ...

British Library Electronic Table of Contents (United Kingdom)

Regional Ocean Modeling System (ROMS v 3.0), a three-dimensional numerical ocean model, was previously enhanced for shallow water applications by including wave-induced radiation stress forcing provided through coupling to wave propagation models (SWAN, REF/DIF). This enhancement made it suitable for surf zone applications as demonstrated using examples of obliquely incident waves on a planar beach and rip current formation in longshore bar trough morphology (Haas and Warner, 2009). In this contribution, we present an update to the coupled model which implements a wave roller model and also a modified method of the radiation stress term based on Mellor (2008, 2011a,b,in press) that includes a vertical distribution which better simulates non-conservative (i.e., wave breaking) processes and ...

Energy Technology Data Exchange (ETDEWEB)

The Energy Laboratory of the University of Houston has developed a computer simulation program called CREAM (Cavity Radiation Exchange Analysis Model) for application to the solar central receiver system. CREAM contains a geometric configuration factor generator based on Nusselt's method. A formulation of Nusselt's method provides the basis for a FORTRAN subroutine NUSSELT. Numerical results from NUSSELT are compared to analytic values and values from Sparrow's method. Sparrow's method is based on a double contour integral and its reduction to a single integral which is a approximated by Gaussian methods. Nusselt's method is adequate for the intended engineering applications, but Sparrow's method is an order of magnitude more efficient in many situations.

CERN Document Server

The association of at least some long gamma-ray bursts with type Ic supernova explosions has been established beyond reasonable doubt. Theoretically, the challenge is to explain the presence of a light hyper-relativistic flow propagating through a massive stellar core without losing those properties. We discuss the role of the jet-star interaction in shaping the properties of the outflow emerging on the surface of the star. We show that the nature of the inner engine is hidden from the observer for most of the evolution, well beyond the time of the jet breakout on the stellar surface. The discussion is based on analytical considerations as well as high resolution numerical simulations. Finally, the observational consequences of the scenario are addressed in light of the present capabilities.

Science.gov (United States)

A full-scale model of a rigid sectorial ring and a set of bellows which was covered with a temperature control layer and a coolant supply equipment was constructed to verify the adequacy of the heating and cooling system designed for JT-60. To cool and heat it effectively and to decrease the temperature differences among the various parts in the vacuum vessel, heater units and cooling pipes were located on the surface of the ring. The temperature control layer is to heat the vessel to 500/degree/sub //C within 70 hours and to maintain it at the temperature for 48 hours for simulating the state of bakeout. Subsequently the vessel is cooled down within about the same time as in heating. Prior to the series of tests, numerical analyses were performed to predict the cooling and heating efficiencies on the model and to examine the method of the temperature regulation.

British Library Electronic Table of Contents (United Kingdom)

Intrinsic spherical aberrations of electron lenses have been the major resolution limiting factor in electron microscopes for several decades. While effective correctors have recently been implemented, an alternative to correct these aberrations is to circumvent them by scaling down lens dimensions by several orders of magnitude. We have fabricated electrostatic lenses exhibiting one micrometer diameter apertures and evaluated their beam forming properties against predictions from numerical ray tracing simulations. It turns out that it is routinely possible to shape a paraxial low-energy electron beam by such micron-sized lenses. Beam profiles have been measured both at a distant detector as well as in a plane close to the lens. It is shown that the lens can form a parallel beam extending ...

International Nuclear Information System (INIS)

An estimation method of plasma density based on surface plasmons theory for surface-wave plasmas is proposed. The number of standing-wave is obtained directly from the discharge image, and the propagation constant is calculated with the trim size of the apparatus in this methods, then plasma density can be determined with the value of 9.1 x 1017 m-3. Plasma density is measured using a Langmuir probe, the value is 8.1 x 1017 m-3 which is very close to the predicted value of surface plasmons theory. Numerical simulation is used to check the number of standing-wave by the finite-difference time-domain (FDTD) method also. All results are compatible both of theoretical analysis and experimental measurement. (authors)

International Nuclear Information System (INIS)

The effect of Kerr nonlinearity on an Airy beam is investigated by using the nonlinear Schroedinger equation. Based on the moments method, the evolution of the Airy beam width in the rms sense is analytically described. Numerical simulations indicate that the central parts of the major lobe of the Airy beam initially give rise to radial compression during propagation in a focusing medium, even though the rms beam width broadens. The partial collapse of the center parts of the major lobe of the beam appear below the threshold for a global collapse. The evolutions of the field distributions of the Airy beams are different during propagation in different Kerr media while the beams still travel along the parabolic trajectory just as the beam propagates in free space.

British Library Electronic Table of Contents (United Kingdom)

Phreatic groundwater pumping is affecting water availability for crops in areas with a shallow water table. This can reduce crop growth and so affect farm income. There is a need for a generic and transparent method to assess the agricultural damage caused by water table drawdown. This paper proposes such a method that consists of 'damage tables' relating agricultural production losses to the groundwater regime for different soil/crop combinations found in Northern Belgium. The damage tables are constructed based on numerous simulations with the agrohydrological model SWAP, in which the bottom boundary conditions are gradually changed to reflect different groundwater regimes. The credibility of the resulting metamodel is assessed in three ways: using (1) field data, (2) an existing local e...

Energy Technology Data Exchange (ETDEWEB)

The experimental foam treatments conducted on two production wells, each in a different miscible flood pool, were described, including aspects of the laboratory work, numerical coning simulation, well selection process and foam placement strategy. One of the wells was injected with a foam-forming surfactant, as well as with a surface-generated preformed foam. Neither of these treatments had any effect on the performance of this well. In the other well, preformed foam treatment was successful in reducing gas-oil ratio, resulting in improved oil production over a 14-month period. Future foam injection projects are being planned, based on the experiences gained in this field experiment. 13 refs., 21 figs.

British Library Electronic Table of Contents (United Kingdom)

A model is developed to represent mechanical strain, stress-enhanced dissolution, and shear dilation as innately hysteretic and interlinked processes in rough contacting fractures. The model is incorporated into a numerical simulator designed to examine permeability change and thermal exchange in chemically active and deformable fractured reservoirs. A candidate engineered geothermal reservoir system (EGS) is targeted. The mechanistic model is able to distinguish differences between the evolution of fluid transmission characteristics of (1) small scale, closely spaced fractures, and (2) large-scale, more widely spaced fractures. Alternate realizations of fracture frequency and scale, exhibiting identical initial bulk permeability, lead to significantly different conclusions regarding perme...

International Nuclear Information System (INIS)

This paper develops a fluid dynamic model using finite difference methods characterizing flow in phantoms simulating in vivo conditions and to compare those results with velocity encoded MR images. The phantom consisted of 1-inch (2.5-cm) tubing with semicircular insert and fluid with viscosity, T1, and T2 comparable to blood. Numeric solutions to Navier-Stokes equations for this system were obtained using finite difference methods, with velocity input function of zero at walls and parabolic at both ends. In resulting color raster (CR) images, color temperature represented velocity value. In velocity-encoded MR images acquired under the same flow conditions, phase is proportional to average velocity during application of flow-encoding gradients. Because these gradients are applied along one direction per acquisition, magnitude and direction of velocity are obtained.

Energy Technology Data Exchange (ETDEWEB)

By means of a numerical simulation model (PHOENICS) the air flow and the distribution of the air concentration in an operating room of an hospital have been calculated. Comparisons were made for 11 different conditions (m[sup 3]/h air flow, speed of the air supply in the centre, the middle and outside of the room), of which six cases are discussed in this article, in order to optimize the air flow pattern in the room. The calculations were conducted for the so-called Howorth system. It is shown that the inlet velocity profile with a higher inlet velocity in the center zone of the hood gives the best results. The calculated results were validated by some measurements and experiments with smoke. 7 figs., 2 tabs., 8 refs.

Energy Technology Data Exchange (ETDEWEB)

In this work, we investigated self-heating related instability in polysilicon thin film transistors (poly-Si TFTs) fabricated on polyimide (PI) substrates. Indeed, when Joule heating becomes relevant, the temperature of the active layer can substantially rise, since the devices are fabricated on thermally insulating substrates. As a result, electrical instability is triggered and attributed to the generation of interface states, due to the Si-H bond breaking, and charge trapping into the gate insulator. In addition, by using 3-dimensional numerical simulations, coupling the thermodynamic and transport models, we analyzed the temperature distribution of the device under operating conditions and found that self-heating is more severe for devices fabricated on plastic substrates.

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.

Energy Technology Data Exchange (ETDEWEB)

Three-dimensional dynamics of a compact toroid (CT) plasmoid, which is injected into a magnetized target plasma region is investigated by using magnetohydrodynamic (MHD) numerical simulations. It is found that the process of the CT penetration into this region is much more complicated than what has been analyzed so far by using a conducting sphere (CS) model. The injected CT suffers from a tilting instability, which grows with the similar time scale as the CT penetration. The instability is accompanied by magnetic reconnection between the CT magnetic field and the target magnetic field, which disrupts the magnetic configuration of the CT. Magnetic reconnection plays a role to supply the high density plasma initially confined in the CT magnetic field into the target region. Also, the penetration depth of the CT high density plasma is examined. It is shown to be shorter than that estimated from the CS model. The CT high density plasma is ...

CERN Document Server

Lyman-limit absorption systems can play many important roles during and after cosmological reionization. Unfortunately, due to the prohibitively large dynamic range required, it is impossible to self-consistently include these systems in cosmological simulations. Using fast and versatile semi-numeric simulations, we systematically explore the spatial distribution of absorption systems during and following reionization. We self-calibrate the resulting number of absorbers to the mean free path (mfp) of the ionizing ultraviolet background (UVB), and present results at a given mfp and neutral hydrogen fraction. We use a simple optical depth criterion to identify the locations of absorbers. Our approach is fairly robust to uncertainties such as missing subgrid structure. Unlike at lower redshifts where the UVB is relatively uniform, at higher redshifts the fluctuations in the UVB and the HII morphology of reionization can drive ...

Energy Technology Data Exchange (ETDEWEB)

In the development of the man-made island siting technology of nuclear power plants, assessing the stability of the seawall against large ocean waves and earthquakes is indispensable. Concerning with the seismic stability of the seawall, prediction of the deformation like sliding and settlement of the seawall during earthquake including the armour units in front of the caisson becomes important factor. For this purpose, the authors have developed the two-dimensional DEM-FEM coupled analysis method (SEAWALL-2D) to predict the deformation of the seawall covered with the armour units during earthquake. In this method, movements of the armour units are calculated in DEM analysis part and deformation of the caisson, rubble moundsand seabed and back fill are calculated in FEM analysis part taking the nonlinearity of the soil materials based on the effective stress into account. Numerical simulations of dynamic centrifuge model tests of the seawall ...

International Nuclear Information System (INIS)

In order to study the boundary conditions for the occurrence of ratcheting, a testing rig was erected. In this rig, small-scale mock-ups of the first wall were tested under mechanical conditions similar to those expected during disruptions. These mock-ups were made of stainless steel AISI 316. They were heated on the surface and water cooled from inside. Disruption forces were simulated by a high-speed hydraulic bending system by impact loads. The mock-ups were heated by a high-frequency generator and cooled by water through two cooling channels.In addition to the experimental work, simplified analytical analysis and numerical calculations were carried out using the finite element program ABAQUS. In these calculations, the temperature and stress fields in the mock-ups were investigated and the influence of the different constitutive equations was studied.The experimental results confirm the ratcheting behaviour predicted by the calculations. It ...

International Nuclear Information System (INIS)

This conference and trade show provided a forum to present new, emerging and existing technologies that will play a key role in ensuring that hydrocarbon exploitation occurs in a sustainable manner, where costs and resource recovery are optimized while minimizing environmental impacts. The challenge of developing maturing basins, frontier lands and unconventional resources was discussed. In particular, the application of thermal recovery methods such as steam assisted gravity drainage (SAGD) for the in-situ recovery of heavy oils and bitumens was reviewed and compared with the performance of alternative enhanced energy recovery technologies such as vapor extraction (VAPEX). The conference was attended by local and international participants in the petroleum sector, including managers, engineers, technologists, geologists, and other petroleum industry specialists. The conference featured 132 presentations addressing topics such as improved recovery of conventional oil and gas; heavy oil ...

Science.gov (United States)

A butterfly valve is a type of flow control device, typically used to regulate a fluid flowing. Currently, FEA is often used to predict the safety in the design of valve disc. Also, the study about the affection of butterfly valve's disc to the valve flow characteristics by using CFD has been done by many researchers. Along with the development of computer technique, design and analysis of computer experiments has becoming more and more important in engineering design and optimization. Hereinto Kriging model is one popular analysis approach for the purpose of creating a cheap ``meta-model'' as a surrogate to a computationally expensive simulation model. In this paper, the numerical analysis considered the strength, pressure loss coefficient and weight of valve disc simultaneously is investigated to improve the shape of a traditional butterfly valve disc. Firstly, an initial model of butterfly valve is made to evaluate the performance of the ...

International Nuclear Information System (INIS)

We simulate mountain-scale coupled heat transfer and gas flow at Yucca Mountain. A coupled rock-gas flow and heat transfer model, TGIF2, is used to simulate mountain-scale two-dimensional transient heat transfer and gas flow. The model is first verified against an analytical solution for the problem of an infinite horizontal layer of fluid heated from below. Our numerical results match very well with the analytical solution. Then, we obtain transient temperature and gas flow distributions inside the mountain. These distributions are used by a transient semianalytical particle tracker to obtain carbon-14 travel times for particles starting at different locations within the repository. Assuming that the repository is filled with 30-year-old waste at an initial areal power density of 57 kw/acre, we find that repository temperatures remain above 60 degrees C for more than 10,000 years. Carbon-14 travel times to the surface are ...

International Nuclear Information System (INIS)

The paper reports on the results of hydrodynamical simulations of the late phase of the common envelope stage of a binary consisting of a 2-solar mass red giant and a 1-solar mass main-sequence companion. The numerical results demonstrate that sufficient energy is released from the orbit to eject the mass within the common envelope without requiring the main-sequence companion to spiral into the white dwarf core of the red giant star. At the end of the simulation the orbital decay time scale increases rapidly to more than 160 yr. The long decay time scale reflects the removal of mass from the common envelope and its subsequent spin-up to near corotation. The ratio of the orbital decay time scale to the mass-loss time scale from the common envelope increases to more than 700, and the mass contained within the common envelope decreases to about 0.01 solar mass or less. It is argued that further orbital decay will be small and ...

CERN Document Server

We compare the predictions of four different algorithms for the distribution of ionized gas during the Epoch of Reionization. These algorithms are all used to run a 100 Mpc/h simulation of reionization with the same initial conditions. Two of the algorithms are state-of-the-art ray-tracing radiative transfer codes that use disparate methods to calculate the ionization history. The other two algorithms are fast but more approximate schemes based on iterative application of a smoothing filter to the underlying source and density fields. We compare these algorithms' resulting ionization and 21 cm fields using several different statistical measures. The two radiative transfer schemes are in excellent agreement with each other (with the cross-correlation coefficient of the ionization fields >0.8 for k 0.6 for k < 1 h/Mpc). When used to predict the 21cm power spectrum at different times during reionization, all ionization algorithms agree with one another at the ...

International Nuclear Information System (INIS)

In addition to the usual circularly symmetric TM/sub 010/ mode used to accelerate particles in an rf linac, there is a large number of modes with cos phi or sin phi dependence, for example the TM/sub 1xx/ modes. These latter modes possess a uniform magnetic (dipole) field near the axis of symmetry and therefore can deflect the beam away from the axis. Any portion of an accelerated beam that is off-axis will drive these modes, so that subsequent portions of the beam will be deflected. This deflected beam will then resonantly drive the same modes in downstream cavities, so that still later portions of the beam will be more severely deflected, and so on. In this paper are reported the results of numerical simulations of this so-called cumulative beam-breakup instability. The simulation assumes that only the TM/sub 110/ mode acts to deflect the beam, and further assumes that this mode does not couple from one accelerating ...

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 of the DLR, which ...

Energy Technology Data Exchange (ETDEWEB)

This second report covers the period August 1 to September 30, 1991. During this period we modified our approach to Task 1 (Hybrid Armature study) with the addition of a numerical simulation of the steady state plasma brush length with a 1-D, four-phase equation of state code. A test fixture was designed and fabrication begun on a plasma brush test fixture for hybrid studies. It will be powered by a 140 kJ, 1 MA capacitor bank on loan to the project. Diagnostic data collection with 20 channels of digitizers and data reduction with a microvax were also arranged. Sliding friction of a solid armature was simulated and found to form a vapor layer at velocities below 2 km/s. The simulations indicated that higher interface pressure results in vapor boundary layer formation at lower velocity. If there is also a minimum interface pressure required to commutate high currents, then an optimum interface pressure ...

Energy Technology Data Exchange (ETDEWEB)

One of the Eddy Current Testing issues in aeronautics is the inspection of fastened structures to detect flaws nearby rivets which can grow because of mechanical stress. EADS and the CEA LIST have started a collaborative work with the support of the Ile-de-France Region to develop a simulation tool of EC fastened structures testing, integrated to the CIVA platform, aimed at conceiving testing methods, optimizing and qualifying it. The volume integral method using the Green dyadics formalism has been chosen in order to get a fast resolution of Maxwell equations. A first milestone was to build a simulation model of multilayer structures testing, thanks to the use of the multilayer Green dyads. Because of the rivet volume, 60 times bigger than the one of a typical flaw, a large number of discretization cells are needed. Therefore an iterative method has been developed in order to numerically solve large calculation zones. ...

Science.gov (United States)

Not Available

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 concentration measurements. A ...

International Nuclear Information System (INIS)

In this study, a high-resolution characteristic-based finite-volume (FV) method on unstructured grids [Int. J. Numer. Method Eng. 50 (2001) 11; Int. J. Heat Fluid Flow 21 (2000) 432] is extended by a matrix-free implicit dual-time stepping scheme for the numerical simulation of steady and unsteady flow and heat transfer with porous media. The method has been used to study the characteristics of a complex problem: flow and heat transfer in a channel with multiple discrete porous blocks, which was originally proposed by Huang and Vafai [J. Thermophys. Heat Transfer 8 (3) (1994) 563]. In addition, flow and heat transfer in a channel partially or fully filled with porous layers and containing solid protruding blocks with constant heat flux on its lower surface are also investigated in details. Hydrodynamic and heat transfer results are reported for both steady and transient flow cases. In particular, the effects of Darcy and ...

Science.gov (United States)

In this dissertation, we design and analyze efficient numerical methods for obtaining accurate solutions to model problems arising in fuel cells. A basic fuel cell model consists of five principles of conservation, namely, mass, momentum, species, charges (electrons and ions), and thermal energy. Overall, transport equations couple with electrochemical processes through source terms to describe reaction kinetics and electro-osmotic drag in the polymer electrolyte. To model multiphase species transport in the porous media and the gas channel of fuel cells, we consider a multiphase mixture model framework. The diffusivity of the two-phase mixture water conservation equation in this model is nonlinear, discontinuous, and degenerate. To handle this difficulty, we developed efficient and fast nonlinear iterative solvers based on the Kirchhoff transformation and nonlinear Dirichlet-Neumann domain decomposition methods. To model the coupling between the multiphase flow in ...

Energy Technology Data Exchange (ETDEWEB)

Numerical simulations are used to predict the migration of radionuclides from the disposal units at Material Disposal Area G through the vadose zone and into the main aquifer in support of a radiological performance assessment and composite analysis for the site. The calculations are performed with the finite element code, FEHM. The transport of nuclides through the vadose zone is computed using a three-dimensional model that describes the complex mesa top geology of the site. The model incorporates the positions and inventories of thirty-four disposal pits and four shaft fields located at Area G as well as those of proposed future pits and shafts. Only three nuclides, C-14, Tc-99, and I-129, proved to be of concern for the groundwater pathway over a 10,000-year period. The spatial and temporal flux of these three nuclides from the vadose zone is applied as a source term for the three-dimensional saturated zone model of the main aquifer that ...

Energy Technology Data Exchange (ETDEWEB)

Calibration-Free Laser-Induced Breakdown Spectroscopy (CF-LIBS) has been proposed several years ago as an approach for quantitative analysis of Laser-Induced Breakdown Spectroscopy spectra. Recently developed refinement of the spectral processing method is described in the present work. Accurate quantitative results have been demonstrated for several metallic alloys. However, the degree of accuracy that can be achieved with Calibration-Free Laser-Induced Breakdown Spectroscopy analysis of generic samples still needs to be thoroughly investigated. The authors have undertaken a systematic study of errors and biasing factors affecting the calculation in the Calibration-Free Laser-Induced Breakdown Spectroscopy spectra processing. These factors may be classified in three main groups: 1) experimental aberrations (intensity fluctuations and inaccuracy in the correction for spectral efficiency of a detection system), 2) inaccuracy in theoretical parameters used for calculations (Stark ...

Energy Technology Data Exchange (ETDEWEB)

An Eulerian/Langrangian model was developed to study the two-phase flow in a two-dimensional, axisymmetric, direct-injection diesel engine, under motored and firing conditions. The liquid phase was tracked in a Lagrangian sense using discrete droplet packets (DDP). The gaseous phase was solved using an implicit, iterative, finite-difference, scheme of the control volume variety. Gaseous-phase turbulence was modeled using a compressible two-equation k/epsilon model. Turbulence effects on the liquid phase were modeled by superimposing a randomly oriented turbulent eddy velocity proportional to the root of the turbulent kinetic energy, on each packet's mean relative velocity. Droplet-droplet interactions were neglected. The general knowledge obtained during validation studies, was used to simulate the direct-injection two-stroke diesel engine under motored conditions. The simulation showed that by beginning of the injection, (15/sup 0/ ...

International Nuclear Information System (INIS)

The study of tropospheric kinetics underlies global change because key greenhouse gases are photochemically active. Modeling of tropospheric chemistry on a global scale is essential because some indirect greenhouse gases are short-lived and interact in a non-linear fashion. It is also extremely challenging, however; the global change grid is extensive in both the physical and temporal domains, and critical lower atmospheric species include the organics and their oxidized derivatives, which are numerous. Several types of optimization may be incorporated into kinetics modules to enhance their ability to simulate the complete lower atmospheric gas phase chemical system. (1) The photochemical integrator can be accelerated by avoiding matrix and iterative solutions and by establishing families. Accuracy and mass conservation are sacrificed in the absence of iteration, but atom balancing is restorable post hoc. (2) Chemistry can be arranged upon the ...

International Nuclear Information System (INIS)

In order to predict or to evaluate the structural reliability of a variety of nuclear components under actual changing operating conditions, it is necessary to take the detailed geometrical features of components into consideration on one hand, and thermo-mechanical loading conditions with space and time variations during the lifetime, on the other hand. This consideration leads the authors to develop a super large-scale structural analysis program based on a parallel computational FEM with a variable scale of PC cluster architecture, whose prototype system is presented in a companion paper in this Symposium. Development of the PC cluster technology is one of the most important subjects in large-scale structural simulation, reducing computer costs and thus making simulation more acceptable to the engineering community. This paper presents the developmental status of the parallel FEM solver, demonstrating computational capabilities with three ...

Science.gov (United States)

Not Available

Energy Technology Data Exchange (ETDEWEB)

No abstract prepared.

Energy Technology Data Exchange (ETDEWEB)

The control problem associated to a class of horizontal-axis fixed-pitch variable-speed low-power wind energy conversion systems, working in the partial load region, consisting in the energy conversion maximization, is approached here under the assumption that the wind turbine model and its parameters are poorly known. Using a new approach derived from the optimum seeking methods category, generically called Maximum Power Point Tracking (MPPT), the proposed control solution aims at driving the average position of the operating point near to optimality. Instead of inducing sinusoidal search signals, the wind turbulence is here used as search disturbance. The high-speed shaft's average rotational speed is slowly adjusted using the Fast Fourier Transform processing of some available measures from the system as an estimate of the operating point's position/distance to optimality. Numerical simulations are used for preliminary ...

International Nuclear Information System (INIS)

Natural convection flow is established in KMRR (Korea Multi-Purpose Research Reactor) reflector tank at the loss of reflector circulator. To simulate the reflector tank natural convection flow with high temperatures at the inner shell and bottom plate due to nuclear heating, experimental and numerical studies in an open cavity with 'L' type heated wall made by the combination of a vertical and horizontal plate were performed. It was confirmed through these studies that the heat transfer rates were highest at the lower region of the vertical plate and the inlet region of horizontal plate and comparatively high at the middle portion of both plates. The heat transfer rate distribution of this trend shows a desirable trend for the effective natural convection cooling of KMRR reflector tank. It was also confirmed that the average Nusselts numbers at the 'L' type heated wall were lower than those obtained from the existing natural convection heat ...

CERN Document Server

This paper addresses the fine-scale axisymmetric structure exhibited in Saturn's A and B-rings. We aim to explain both the periodic microstructure on 150-220m, revealed by the Cassini UVIS and RSS instruments, and the irregular variations in brightness on 1-10km, reported by the Cassini ISS. We propose that the former structures correspond to the peaks and troughs of the nonlinear wavetrains that form naturally in a viscously overstable disk. The latter variations on longer scales may correspond to modulations and defects in the wavetrains' amplitudes and wavelength. We explore these ideas using a simple hydrodynamical model which captures the correct qualitative behaviour of a disk of inelastically colliding particles, while also permitting us to make progress with analytic and semi-analytic techniques. Specifically, we calculate a family of travelling nonlinear density waves and determine their stability properties. Detailed numerical ...

Energy Technology Data Exchange (ETDEWEB)

One of the main objectives of the Finnish Research Programme on Climate Change (SILMU) has been to assess the possible impacts of future changes in climate due to the enhanced greenhouse effect on natural systems and human activities in Finland. In order to address this objective, it was first necessary to specify the types of climate changes to be expected in the Finnish region. Estimates of future climate are conventionally obtained using numerical models, which simulate the evolution of the future climate in response to radiative forcing due to changes in the composition of the atmosphere (i.e. of greenhouse gases and aerosols). However, there are large uncertainties in the model estimates because current knowledge and understanding of atmospheric processes remains incomplete. Since accurate predictions of climate change are not available, an alternative approach is to develop scenarios. These are plausible projections which reflect the best ...