Objectives(i) Investigate the correlation between Antarctic sea-ice and equatorial sea-surface temperature anomalies in a realistically forced oceanmodelsimulation of the last 50 years. (ii) Determine whether and how the enormous seasonal change in distribution of sea-ice modifies the seasonal cycle at the Equator. (iii) Determine the detailed pathways of wave propagation both in a historically-forced simulation and in response to realistic perturbations. (iv) Quantify the amplitude of the response i [continued...]DescriptionIt is well known that the equatorial ocean-atmosphere system plays a key role in global climate events such as the El Nino-Southern Oscillation (ENSO) phenomenon. There is now compelling evidence that changes in the Antarctic can strongly and quickly affect the equatorial ocean and the ENSO cycle. Observations demonstrate statistically ...
Regional OceanModeling System (ROMS v 3.0), a three-dimensional numerical oceanmodel, 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 ...
A parameterized surface wave-induced vertical mixing (B"v) is incorporated into the Modular OceanModel version 4 (MOM4). A comparison of the two numerical experiments with and without B"v shows that B"v can significantly improve the upper-ocean (20-100m) simulation in summer. The simulated upper-ocean temperature errors are reduced in summer due to the surface wave-induced vertical diffusive heat flux. The non-breaking-wave-induced vertical mixing can increase the probability of the simulated SST biases between -1^oC and 1^oC from 64% to 76% in the Southern Hemisphere (60^oS-10^oS) in January, and from 66% to 75% in the Northern Hemisphere (10^oN-60^oN) in July. The averaged mixed layer depth (MLD) simulated by the MOM4 without B"v is 14.4m shallower than the observations in 10^oS-60^oS i...
A study for Cs-137 radionuclide dispersion in the marine environment through of compartmental model (Box Model) is presented. The modelsimulates the surface water contamination caused by direct atmospheric deposition, surface wash off, desorption from sediments and transfer with the ground water of accidentally released radionuclides. For this study the model was applied to the North Sea, near to Sellafield, based on the transfer coefficients obtain at the literature. The results obtained are in good agreement with the literature, being that the model developed can be applied in to the brazilian coastal regions. (author). 7 refs, 7 figs.
The Subseabed Disposal Program has recently completed a 30-day, 0.286-scale laboratory simulation of an In Situ Heat Transfer Experiment (ISHTE), scheduled for eventual deployment in the deep ocean 600 nautical miles north of Hawaii. 10 figures.
The Subseabed Disposal Program has recently completed a 30-day, 0.286-scale laboratory simulation of an In Situ Heat Transfer Experiment (ISHTE), scheduled for eventual deployment in the deep ocean 600 nautical miles north of Hawaii. 10 figures.
The tidal motion of the ocean water through the ambient magnetic field, generates secondary electric field. This motionally induced electric field can be detected in the sea or inland and has a potential for electrical soundings of the Earth. A first goal of the paper is to gain an understanding of the global distribution of the electric signal due to tidal ocean flow. We simulate the electric signals for two tidal constituents - lunar semidiurnal (M2) and diurnal (O1) tides. We assume a realistic Earth's conductivity model with a surface thin shell and 1-D mantle underneath. Simulations demonstrate that in some coastal regions the amplitudes of the electric field can reach 100 mV/km and 10 mV/km for M2 and O1 tides respectively. The changes of lithosphere resistance produce detectable changes in the tidal electric signals. We show that our predictions are in a good agreement with ...
period T. Also shown is the returned chirp from a specular reflector at ...... algorithms is a Brown (1977) model of the return waveform, which assumes a ...... Townsend, W.F., 1980: An initial assessment of the performance achieved by ...
The climate modeling community has focused recently on improving our understanding of certain processes, such as cloud feedbacks and ocean circulation, that are deemed critical to climate-change prediction. Although attention to such processes is warranted, emphasis on these areas has diminished a general appreciation of the role played by the large-scale dynamics of the extratropical atmosphere. Lack of interest in extratropical dynamics may reflect the assumption that these dynamical processes are a non-problem as far as climate modeling is concerned, since general circulation models (GCMs) calculate motions on this scale from first principles. Nevertheless, serious shortcomings in our ability to understand and simulate large-scale dynamics exist. Partly due to a paucity of standard GCM diagnostic calculations of large-scale motions and their transports of heat, momentum, ...
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.
This paper presents results, numeric and experimental, due to installation operation of a Risers Uphold Sub-Surface Buoy, (BSR). This kind of installation developed by PETROBRAS is unique in the world. The work of BSR installation was based on a numeric pre analysis to verify the system and determine the main parameters to be experimentally verified. The second phase of the work was the experimental analysis in a deep water ocean basin. s. The work describes the BSR and their main accessories, the experimental environment and the model constructed in aluminum in a 1:12 scale and the main results. (author)
Understanding the processes responsible for coastal change is important for managing our coastal resources, both natural and economic. The current scientific understanding of coastal sediment transport and geology suggests that examining coastal processes at regional scales can lead to significant insight into how the coastal zone evolves. To better identify the significant processes affecting our coastlines and how those processes create coastal change we developed a Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) Modeling System, which is comprised of the Model Coupling Toolkit to exchange data fields between the oceanmodel ROMS, the atmosphere model WRF, the wave model SWAN, and the sediment capabilities of the Community Sediment Transport Model. This formulation builds upon ...
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 ...
Separate abstracts were prepared for 20 papers in this conference report and includes all papers in the ''engineering'' section (program E). Ocean thermal energy conversion (OTEC), moorings, cables, corrosion protection, corrosion testing, coatings, offshore platforms, pipelines, OTEC models and pilot plants, and energy transfer are topics discussed.
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 ...
BackgroundFood supply from the ocean is constrained by the shortage of domesticated and selected fish. Development of genomic models of economically important fishes should assist...Full Text Available
Ocean thermal energy conversion (OTEC) plants discharge large volumes of cold water into the upper ocean. A three-dimensional, limited-area model was developed to investigate the regional influence of the far-field effluent plume created by the negatively buoyant discharge. The model was applied to discharges from a 40-MW/sub e/ OTEC plant into coastal waters characterized by various ambient ocean conditions. A typical ambient temperature structure and nutrient distribution, as well as the behavior of the effluent plume itself, were strongly modified by the discharge-induced circulation. Although temperature perturbations in the plume were small, upward entrainment of nutrients from below the thermocline was significant. The regional influence of discharges from an 80-MW/sub e/ OTEC plant, the interactions between the discharges from two adjacent 40-MW/sub e/ OTEC plants, and the ...
Political leaders in numerous nations argue for an upper limit of the global average surface temperature of 2 K above the pre-industrial level, in order to attempt to avoid the most serious impacts of climate change. This paper analyzes what this limit implies in terms of radiative forcing, emissions pathways and abatement costs, for a range of assumptions on rate of ocean heat uptake and climate sensitivity. The primary aim is to analyze the importance of ocean heat uptake for radiative forcing pathways that temporarily overshoot the long-run stabilization forcing, yet keep the temperature increase at or below the 2 K limit. In order to generate such pathways, an integrated climate-economy model, MiMiC, is used, in which the emissions pathways generated represent the least-cost solution o...
This book, an updated version of the 1985 edition, contains thirteen chapters, beginning with a preface which provides the objective of the book. The primary objective is to offer a comprehensive survey of the biological, hydrological, mathematical, and biochemical aspects of stream, lake, estuary, and ocean pollution analysis. The book also contains ten appendices of useful tables and nomographs of pertinent data. This book provides a very good summary and review of stream, lake, estuary, and ocean pollution. This book is recommended for environmental engineering students, environmental consulting engineers, and regulatory personnel. It provides an excellent summary of the field of stream and lake analysis and modeling. Every chapter includes a significant number of questions and pertinent references.
The use of global three-dimensional (3-D) models with satellite observations of CO2 in inverse modeling studies is an area of growing importance for understanding Earth s carbon cycle. Here we use the GEOS-Chem model (version 8-02-01) CO2 mode with multiple modifications in order to assess their impact on CO2 forward simulations. Modifications include CO2 surface emissions from shipping (0.19 PgC yr 1), 3-D spatially-distributed emissions from aviation (0.16 PgC yr 1), and 3-D chemical production of CO2 (1.05 PgC yr 1). Although CO2 chemical production from the oxidation of CO, CH4 and other carbon gases is recognized as an important contribution to global CO2, it is typically accounted for by conversion from its precursors at the surface rather than in the free troposphere. We base our model 3-D spatial distribution of CO2 chemical production on monthly-averaged loss rates of CO (a ...
We explore the impact of obliquity variations on planetary habitability in hypothetical systems with high mutual inclination. For the hypothetical systems, we restrict our exploration to systems consisting of a solar-mass star, an Earth-mass planet at 1 AU, and 1 or 2 giant planets. We verify that these systems are stable for 108 years with N-body simulations. We then calculate the obliquity variations induced by the orbital architecture on the Earth-mass planets. We find that in some cases the spin axes can rotate through 360 degrees in as little as 10,000 years (John is that right? Can you look through the systems and find the most extreme case of obliquity variation?) Next, we run energy balance models (EBM) on the terrestrial planets to assess surface temperature and ice coverage on the planets' oceans. Finally, we explore differences in the outer edge of the habitable zone for planets with rapid obliquity variations. ...
It is often advantageous to generate power with combinations of wind and ocean waves. In fact ocean waves, their generation, propagation, dissipation are directly related to wind velocity and its duration oven the sea. In this paper an attempt has been made to demonstrate statistically to present some advantages with combined wind and ocean wave power generation. Even though many conceptual techniques and methods are possible to harness combined power generation, it is important to test feasibility of combined output as well as individual outputs mathematically. One of the major advantages of combined wind and wave power generation is to improve probability of continuous power supply (it minimises the interruptions and compensates power fluctuations with one another). Some of the major wave characteristics like wave Height (H), Time period (T), Wave length (L) significantly influence wave power generation. Interestingly, ...
An observational and modeling study was performed over the Arabian Gulf region to investigate the coastal circulations and aerosol transport in the area. Climatological data and observations from the United Arab Emirates' Unified Aerosol Experiment were used to develop a better understanding of the complex meteorological processes in the Arabian Gulf region. Climatological data suggests that sea breezes occur on more than 77 percent of days in all months of the year and land breezes occur on more than 70 percent of the days. The occurrence of the sea and land breeze circulations are higher (90-99 percent) during the summer months when large-scale weather patterns are quiescent. Measurements of aerosol concentrations taken during the Unified Aerosol Experiment(2004) experiment are used to investigate aerosols, namely dust, transported in the Arabian Gulf region. Vertical profiles of dust concentration along with vertical profiles of potential temperature and wind ...
This site explains how temperature, pressure, and salinity work together to determine the density of ocean water. The three density layers of the ocean are described by means of text description and a graphic illustration.
Remotely sensed ocean color data and numerical modeling have been used to study the phenology of both spring and fall phytoplankton blooms (FPBs) in the Nova Scotian Shelf (NSS)-Gulf of Maine (GoM) region. The ocean color data reveal a general pattern of westward progression of the spring phytoplankton bloom (SPB), and an eastward progression of the FPB in the NSS-GoM region. The spatial pattern of mean chlorophyll concentration in spring is similar to that in fall, with a lower concentration in the NSS and higher in the GoM. Interannually, there is a weak but significant tendency for years with earlier (delayed) SPBs to be followed by delayed (earlier) FPBs, but the mean chlorophyll concentrations during SPBs are not correlated with those during FPBs. The interannual variability of SPB ti...
Sep 1, 2011 ... 2005. Computer simulationmodeling of recreation use: current status, case studies, and future directions. Gen. Tech. Rep. RMRS-GTR-143. ...
Decadal-scale climate variations over the Pacific Ocean and its surroundings are strongly related to the so-called Pacific decadal oscillation (PDO) which is coherent with wintertime climate over North America and Asian monsoon, and have important impacts on marine ecosystems and fisheries. In a near-term climate prediction covering the period up to 2030, we require knowledge of the future state of internal variations in the climate system such as the PDO as well as the global warming signal. We perform sets of ensemble hindcast and forecast experiments using a coupled atmosphere-ocean climate model to examine the predictability of internal variations on decadal timescales, in addition to the response to external forcing due to changes in concentrations of greenhouse gases and aerosols, volcanic activity, and solar cycle variations. Our results highlight that an initialization of the upper-ocean state ...
A simple analytical/numerical model has been developed for computing the evolution, over periods of up to a few hours, of the current and temperature profile in the upper layer of the ocean. The model is based upon conservation laws for heat and momentum, and employs an eddy diffusion parameterisation which is dependent on both the wind speed and the wind stress applied at the sea surface. Other parameters such as the bulk-skin surface temperature difference and CO$_2$ flux are determined by application of the Molecular Oceanic Boundary Layer Model (MOBLAM) of Schluessel and Soloviev. A similar model, for the current profile only, predicts a temporary increase in wave breaking intensity and decrease in wave height under conditions where the wind speed increases suddenly, such as, for example, during gusts and squalls. The model results are ...
This series of animations depicts the results of a computer simulation of the effects of deforestation on the surface temperature and precipitation in the Amazon region. Six animations are shown - a control simulation, a deforestation simulation, and the difference between the two for both surface temperature and rainfall. The period simulated is the second half of 1987.
DescriptionRemote sensing of the sea surface from satellites in near-polar orbits has contributed greatly to our understanding of the links between physical and biological processes in marine systems. However most of this progress has been made in open oceanic waters or major upwelling areas, and many unresolved problems are encountered in coastal regions and shelf seas. In these optically complex waters, quantitative remote sensing requires a more sophisticated interpretation strategy than that implemente [continued...
In recent years the modelling of interannual climate variability has been studied, the atmospheric energy and water cycles, and climate simulations with the ECHAM3 model. In addition, the climate simulations of several models have been compared with special emphasis in the area of northern Europe
;Contents: Comparison of Simulation Modules of TRANSYT and INTEGRATION Models; Evaluation of SCATSIM-RTA Adaptive Traffic Network SimulationModel; Comparison NETSIM, NETFLO I, and NETFLO II Traffic SimulationModels for Fixed-Time Signal Control; Traffic Flow Simulation Through Parallel Processing; Cluster Analysis as Tool in Traffic Engineering; Traffic Platoon Dispersion Modeling on Arterial Streets; Hybrid Model for Estimating Permitted Left-Turn Saturations Flow Rate; and Passing Sight Distance and Overtaking Dilemma on Two-Lane Roads.
The impacts of climate change on crop productivity are often assessed using simulations from a numerical climate model as an input to a crop simulationmodel. The precision of these predictions reflects...Full Text Available
In connection with the ERP project 'Dynamic modelling of staged gasification processes' a gasification simulator has been constructed. The simulator consists of: a mathematical model of the gasification process developed at Technical University of Denmark, a user interface programme, IGSS, and a communication interface between the two programmes. (BA)
2D simulationmodeling of a river basin is based on notion of a River Basin as an opened non equilibrium nature system. All components of the system (elevation, water, soil properties and others) are linked to single multi - functional model. The objective of flood simulation is one the model functions. The objectives of the paper are as follows: a mean scale basin floods computing and the model validation; a large scale basin simulation and demands for its validation; results and demands for flood hazard assessment. (orig.)
This paper presents an approach to comparing computer run time of building simulation programs. The computing run time of a simulation program depends on several key factors, including the calculation algorithm and modeling capabilities of the program, the run period, the simulation time step, the complexity of the energy models, the run control settings, and the software and hardware configurations of the computer that is used to make the simulation runs. To demonstrate the approach, simulation runs are performed for several representative DOE-2.1E and EnergyPlus energy models. The computer run time of these energy models are then compared and analyzed.
terms of 11) numerical models, 2D SPH simulations, and analytic models. The 1D model was useful for determining the evolution over very long timescales. ...
Engine model: An engine and governor model adapted from a heavy lift helicopter simulation provides a restistic time delay between aerodynamic rotor ...
The metabolism of organic matter in the coastal regions of the ocean may significantly affect the oceanic carbon budget. This paper describes the high percentage (30%) of oceanic oxidation that occurs in coastal areas and the impact of this metabolism on the carbon cycle and climate of earth. Organic metabolism in the ocean appears to be a source of carbon dioxide release into the atmosphere. Oxidation in the coastal zone is of special interest, as it is likely influenced by anthropogenic activity. Recommendations for future research on this topic are proposed. 129 refs., 2 figs., 6 tabs.
Part I of this thesis is on long-term iceberg collision-risk assessment methods for fixed offshore structures. Estimates for the long-term probability that an iceberg will hit a fixed offshore structure are based on estimates of the total volume of produced iceberg, a model for the size distribution of icebergs and a description of how iceberg trajectories ``fill`` the plane (the ocean`s surface). Part II of the thesis discusses methods for analysis of low resolution data in the sense that the size of the area covered by each pixel is not small compared to the dominating floes. Still, this type of data can provide estimates of the ice extent (and hence the ice edge) and ice concentration. The estimates of the ice edge (at least) seem to be physically significant in the sense that ``in some way`` they can reflect ice drift in an image time sequence. A priori this is far from obvious. A central idea is that the ``ice edge`` is close to a ...
The Finnish simulation code APROS and especially its 5-equation model is applied to modelling the horizontal steam generator. Different nodalizations are used in the secondary side of different models. Simulation results of the stationary state run are compared with results of RELAP5/MOD2 calculations and with an experimental plant data. (2 refs., 3 figs., 4 tabs.).
The Finnish simulation code APROS and especially its 5-equation model is applied to modelling the horizontal steam generator. Different nodalizations are used in the secondary side of different models. Simulation results of the stationary state run are compared with results of RELAP5/MOD2 calculations and with an experimental plant data. (2 refs., 3 figs., 4 tabs.).
The authors review two types of simulators for the analysis of cosmic ray neutron-induced soft errors (SE's). One of them is the neutron-induced soft error simulator (NISES). A recently proposed nuclear reaction theory forms the foundation for the nuclear reaction database used in NISES. The other simulator, the simplified simulator MBGR, is based on a modified version of the burst generation rate (BGR) model. Both simulators accurately simulate neutron-induced SE rates (SER's). MBGR actually provides an easier and quicker estimation of neutron-induced SER's than NISES. On the other hand, NISES covers more applications; it simulates neutron-induced charge collection, multiple-bit SE, and [alpha]-induced SE analysis.
The authors review two types of simulators for the analysis of cosmic ray neutron-induced soft errors (SE's). One of them is the neutron-induced soft error simulator (NISES). A recently proposed nuclear reaction theory forms the foundation for the nuclear reaction database used in NISES. The other simulator, the simplified simulator MBGR, is based on a modified version of the burst generation rate (BGR) model. Both simulators accurately simulate neutron-induced SE rates (SER's). MBGR actually provides an easier and quicker estimation of neutron-induced SER's than NISES. On the other hand, NISES covers more applications; it simulates neutron-induced charge collection, multiple-bit SE, and #alpha#-induced SE analysis
The paper sets up a simulation of a geothermal heat pump with a non-azeotropic mixture. The model is modified and verified with experimental data. The results of the simulation show that the systematic model can predict the performance within {+-}12% of the experimental data. As a result, some improvements can be provided on the basis of the simulation platform. (Author)
In this project we determined primary production and optical variability in the shelf and slope waters off of Cape Hatteras, N.C. These processes were addressed in conjunction with other Ocean Margins Program investigators, during the Spring Transition period and during Summer. We found that there were significant differences in measured parameters between Spring and Summer, enabling us to develop seasonally specific carbon production and ecosystem models as well as seasonal and regional algorithm improvements for use in remote sensing applications.
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
... program incorporates the proven concepts of advanced distributed learning, modeling and simulation, and electronic digital libraries (Figure 1) to ...
12], the Predictive Capability Maturity Model (PCMM) [ref. 13], the Uncertainty Structure scale [ref. 11], and the Simulation Readiness Level (SRL) scale ...
Both the accuracy and real-time are quite strictly required for the RCS thermal-hydraulic model of the simulation of nuclear power station. So, a large amount of reasonable simplifications should be carried on. The author introduces the RCS thermal-hydraulic model of the full scope simulator of Daya Bay Nuclear Power Station, and emphasizes on the fundamental hypotheses and computation method of its kernel part--DEFI programme.
This is our response to a comment by Walter Eifler on our paper `A simple model for the short-time evolution of near-surface current and temperature profiles' (arXiv:physics/0503186, accepted for publication in Deep-Sea Research II). Although Eifler raises genuine issues regarding our model's validity and applicability, we are nevertheless of the opinion that it is of value for the short-term evolution of the upper-ocean profiles of current and temperature. The fact that the effective eddy viscosity tends to infinity for infinite time under a steady wind stress may not be surprising. It can be interpreted as a vertical shift of the eddy viscosity profile and an increase in the size of the dominant turbulent eddies under the assumed conditions of small stratification and infinite water depth.
A seasonal forecasting technique to produce probabilistic and deterministic streamflow forecasts for 23 basins in Norway and northern Sweden is developed in this work. Large scale circulation and moisture fields, forecasted by the ECHAM4.5 model 4 months in advance, are used to forecast spring flows. The technique includes model output statistics (MOS) based on a non-linear Neural Network (NN) approach. Results show that streamflow forecasts from Global Circulation Model (GCM) predictions, for the Scandinavia region are viable and highest skill values were found for basins located in south-western Norway. The physical interpretation of the forecasting skill is that stations close to the Norwegian coast are directly exposed to prevailing winds from the Atlantic ocean, which constitute the principal source of predictive information from the atmosphere on the seasonal timescale.
The Atmospheric Release Advisory Capability (ARAC) at the Lawrence Livermore National Laboratory uses the U.S. Navy's Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) to supply high-resolution wind data for use in its real-time dispersion modeling system. ARAC has used COAMPS products to support several events and exercises, and COAMPS forecasts appear accurate, based on qualitative examination. Recently ARAC has developed a quantitative verification system which calculates COAMPS error and bias statistics, comparing COAMPS forecasts of various lengths with observational data. This paper shows how this system has been used to guide ARAC operators, who need an estimate of the likely behavior of COAMPS forecasts of various lengths in different regions, seasons, and weather patterns.
This paper presents a new approach to improve the in silico modeling of ultrasound transducer arrays. While current simulation tools accurately predict the theoretical element...Full Text Available
The spectre of major climate change caused by the greenhouse effect has generated intensive research, heated scientific debate and a concerted international effort to draft agreements for the reduction of greenhouse gas emissions. This report of Scientific Perspectives on the greenhouse problem explains the technical issues in the debate in language readily understandable to the non-specialist. The inherent complexities of attempts to simulate the earth's climate are explained, particularly with regard to the effects of clouds and the circulation of the oceans, which together represent the largest factors of uncertainty in current global warming forecasts. Results of the search for the 'greenhouse signal' in existing climate records aredescribed in chapter 3 (part two). Chapter 5 (part two) develops a projection of 21st-century warming based on relatively firm evidence of the earth's actual response to known increases in greenhouse gas ...
A detailed knowledge about the dynamics of phytoplanktonic photosynthesis and respiration is crucial for the determination of primary productivity in open oceans as well as for biotechnological applications. The dynamics are best studied in photobioreactors that are able to simulate natural conditions in such, that light can be modulated not only diurnally but also mimicking effects of solar elevation angle from sunrise to sunset, variable cloudiness, light modulation in refractory sun flecks due to water waves, or light intermittence due to turbulent flow in dense suspensions. In addition, high performance photobioreactors ought to be able to monitor in real time photosynthetic and respiratory activities as well as culture growth. Here, we demonstrate performance of a newly designed bench...
Cardiovascular modelling has been a major research subject for the last decades. Different cardiac models have been developed at a cellular level as well as at the whole organ level. Most of...Full Text Available
This paper describes a simulationmodel for predicting end-use energy consumption in residential sectors of a city or region. In this model, the annual energy consumption of a dwelling is simulated from the occupants' schedule of living activities, weather data and energy efficiencies of appliances and dwellings. By summing up the simulation results for various household categories, total energy consumption for the residential sector in a region can be estimated. In this paper, energy consumption for Osaka City is simulated. The result is compared with statistical data. The effects of energy efficiency standards and urban heat island phenomena are examined. (author)
Computer simulations of agricultural finance and management may be helpful to students preparing for agricultural careers. Simulations would provide problem definition and problem solving in realistic interactions with complex and ambiguous situations similar to those encountered in farming. The paper describes a major financial problem for farmers, explores the implications of that problem for students, and outlines a computer-based simulationmodel which would give students experience in dealing with the problem. Data sources for simulations, potentials of interactive exercises among students, and advantages of using real data with existing ambiguities are described. The board game "AGRICULTURE" is discussed as a prototype for simulation. (LFL)
The Formerly Utilized Sites Remedial Action Program (FUSRAP) is to evaluate the radiological conditions at former MED-US AEC sites. Purpose of the Ocean FUSRAP program is to assess the feasibility of ocean disposal of FUSRAP waste which contains trace natural radioactive materials. This paper presents soil characterization information on the Middlesex, NJ, Sampling Plant site, and reports preliminary ocean dispersion calculations. (DLC)
A series of speeches is presented which cover political aspects of ocean energy, legislation, the status of the DOE Ocean Energy Program, the Maritime Administration OTEC Support Program. Six workshop reports are given covering ocean energy technology and systems concepts, power systems concepts, environmental developments, energy transfer and utilization, construction and test capabilities, and industry development. (LEW)
The magnicon is a highly attractive candidate to be the RF source for a future multi-Tev linear collider. Physical models and computer codes have been developed which can provide start-to-end self-consistent simulations of a magnicon, including precise simulations of the high-convergence electron gun, RF-system, magnetic system, and beam collector. The 3-D beam dynamics simulations include realistic fields, finite beam size and transverse space charge effects. The codes allow one to provide steady-state simulations of the entire tube, so as to evaluate transient process of magnicon excitation, parasitic mode self-excitation, stability analysis, and tolerance analysis. The results of the simulations are found to be in good agreement with magnicon experiments. A brief description of the physical models and simulation codes ...
The adjoint of an ocean general circulation model is used as a tool for investigating the causes of changes in ENSO SST indices. We identify adjoint Kelvin and Rossby waves in the sensitivities to sea level and wind stress at earlier times, which can be traced back for more than a year through western and weak eastern boundary reflections. Depending on the thermocline depth the first and second baroclinic modes are excited. The sensitivities to the heat flux and SST are local and decay in about a month. The sensitivities to the fluxes are converted into the influence of SST using the adjoint of a statistical atmosphere model. Focusing on SST perturbations in the index region itself, we recover, up to a scale factor, the delayed oscillator concept.
The purpose of the Ocean FUSRAP Program, which began in March 1981, is to assess the technical, environmental, and institutional feasibility of disposing, in the ocean and on the ocean floor, of FUSRAP waste which contains traces of natural radioactive materials. The initial planning has focused on the Middlesex, New Jersey, Sampling Plant site and surrounding properties, which contain on the order of 100,000 metric tons of material. The Belgian Congo uranium ore and other uranium ores used by the United States were handled at the sampling plant site and have since been removed.
A detailed component-based simulationmodel of a geothermal heat pump system has been calibrated to monitored data taken from a family housing unit located at Fort Polk, Louisiana. The simulationmodel represents the housing unit, geothermal heat pump, ground heat exchanger, thermostat, blower, and ground-loop pump. Each of these component models was tuned to better match the measured data from the site. These tuned models were then interconnected to form the system model. The system model was then exercised in order to demonstrate its capabilities.
A detailed component-based simulationmodel of a geothermal heat pump system has been calibrated to monitored data taken from a family housing unit located at Fort Polk, Louisiana. The simulationmodel represents the housing unit, geothermal heat pump, ground heat exchanger, thermostat, blower, and ground-loop pump. Each of these component models was 'tuned' to better match the measured data from the site. These tuned models were then interconnect to form the system model. The system model was then exercised in order to demonatrate its capabilities.
This paper introduces a new way of estimating parameters in a Brownian motion regime switching asset model to incorporate volatility clustering. The regime switching model is then applied to pricing of up-and-in barrier call options. We take the probability of crossing the barrier between simulation points into account, and we increase accuracy in simulations by importance sampling. The regime switching model is compared to the Normal Inverse Gaussian model and the traditional Black-Scholes model, and option prices from the regime switching model are compared to the closed form expression of up-and-in barrier calls in a Black-Scholes market.
In 2009 the European Space Agency satellite mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer) was launched. Its objectives are the precise and detailed determination of the Earth's gravity field and geoid. Its core instrument, a three axis gravitational gradiometer, measures the gravity gradient components V xx , V yy , V zz and V xz (second-order derivatives of the gravity potential V) with high precision and V xy , V yz with low precision, all in the instrument reference frame. The long wavelength gravity field is recovered from the orbit, measured by GPS (Global Positioning System). Characteristic elements of the mission are precise star tracking, a Sun-synchronous and very low (260 km) orbit, angular control by magnetic torquing and an extremely stiff and thermally stable instrument environment. GOCE is complementary to GRACE (Gravity Recovery and Climate Experiment), another satellite gravity mission, launched in 2002. While ...
A simple and continuous model for the on-state current of polysilicon thin-film transistors, suitable for implementation in circuit simulators, is presented. The model includes the potential barrier at the grain boundaries, the channel length modulation and the excess current due to impact ionization. Comparison between measured output characteristics and the model shows excellent agreement over wide range of bias voltages and for devices with different gate lengths.
This document describes the High Performance Electrical Modeling and Simulation (HPEMS) Global Verification Test Suite (VERTS). The VERTS is a regression test suite used for verification of the electrical circuit simulation codes currently being developed by the HPEMS code development team. This document contains descriptions of the Tier I test cases.
This document describes the High Performance Electrical Modeling and Simulation (HPEMS) Global Verification Test Suite (VERTS). The VERTS is a regression test suite used for verification of the electrical circuit simulation codes currently being developed by the HPEMS code development team. This document contains descriptions of the Tier I test cases.
This demonstration project uses modern simulation techniques to illustrate the important technologies and design variables that an auto-designer would consider in production a high efficiency, low emissions vehicle. Simulation and modeling techniques use the idea of capturing the relationships between real components of the systems with mathematical equations. These equations are then solved on a computer to simulate the behavior or performance of the system under various conditions. In the current demonstration project, we focus on many variations of a hydrogen-powered vehicle.
A distributed watershed hydrologic model soil and water assessment tool (SWAT, 2000 Version) was applied to simulate stream flow and nutrient loadings (including total nitrogen (TN) and total phosphorus (TP)) in Gucheng Lake Basin, which is located in the lower reach of Yangtze River Basin. Model parameters related to stream flow and nutrient loadings were calibrated and validated by the observed value, and the simulation showed that the simulated values were reasonably comparable to the observed data, suggesting the validity of SWAT model. The simulated surface runoff and nutrient yield between two stages, viz. 1981-2000 and 1951-1960, were compared, basically considering land reclamation in the mid-1970s. The result indicated that, the average annual runoff and output of TN and TP in 198...
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.
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.
The report provides and overview of calculation models for the simulation of airflows and deals comprehensively with field and multi-zone models as well as the coupling of individual zone and multi-zone models. Examples of calculations are given. figs., tabs., refs.
A new 3D model is developed to simulate the self-oscillation of the elongated vocal folds. This model allows for large deformation and longitudinal displacement. The displacement boundary condition...Full Text Available
This paper reports on Monte Carlo path integral simulations of a Bianchi Type IX cosmology which yield a wave function which reflects the known chaotic dynamics of the classical system and predicts a significant probability for a correlation between large universe volume and large anisotropy.
The percentages of total airflows over the nasal respiratory and olfactory epithelium of female rabbits were calculated from computational fluid dynamics (CFD) simulations of steady-state inhalation....Full Text Available
The author describes the purpose, classification and main process of acceptance test of full scope simulator of Daya Bay Nuclear Power Plant, including the correction of non-conformance items which are discovered during the performance of acceptance tests. The results of the acceptance tests show that the model accepted by the full scope simulator of Daya Bay Nuclear Power Plant is fully able to cope with the simulation of normal transients and incidental transients and the performance of the simulator indeed compiled with the technical specifications which are defined n the relevant contracts.
Details of the interaction of photons with tissue phantoms are elucidated using Monte Carlo simulations. In particular, photon sampling volumes and photon pathlengths are determined for a variety of scattering and absorption parameters. The Monte Carlo simulations are specifically designed to model light delivery and collection geometries relevant to clinical applications of optical biopsy techniques. The Monte Carlo simulations assume that light is delivered and collected by two, nearly-adjacent optical fibers and take into account the numerical aperture of the fibers as well as reflectance and refraction at interfaces between different media. To determine the validity of the Monte Carlo simulations for modeling the interactions between the photons and the tissue phantom in these geometries, the simulations were compared to measurements of ...
This paper describes a 1993 field experiment entitled Marine Aerosol Properties and Thermal Imager Performance Trial (MAPTIP) conducted by NATO AC/243 Panel 04/RSG.8 and 04/RSG.5 in the Dutch coastal waters. Objectives were: to improve and validate vertical marine aerosol models by providing an extensive set of aerosol and meteorological measurements, within a coastal environment at different altitudes and for a range of meteorological conditions; make aerosol and meteorological observations in the first 10 m of the ocean surface with a view to extending existing aerosol models to incorporate near-surface effects; and to assess marine boundary layer effects on thermal Imaging systems. Aerosol and meteorological instruments, as well as thermal imagers and calibrated targets, were utilized. This network of instrumentation has provided a comprehensive database of aerosol size distribution profiles and relevant meteorological ...
The point-defect-impurity pair diffusion model proposed recently by Mulvaney and Richardson is adopted and modified to simulate the coupled diffusion of phosphorus and self-interstitials in phosphorus-implanted silicon. The assumption of implantation-induced, but empirically determined initial interstitial distributions of Gaussian shape allows a simulation of the net effect of transient enhanced diffusion. As a result an improved modeling of phosphorus diffusion in silicon is achieved for a broad range of ion-implantation and annealing conditions. (author).
This work deals with a time-resolved optical study of the avalanche and streamer formation phases leading to breakdown in atmospheric nitrogen. The authors present the results obtained for nitrogen, from experiments and two-dimensional modelsimulations. This model is used to obtain a better insight in the relevant mechanisms and processes by a comparison of measurements and simulation data. The trends of externally measured quantities correspond with those predicted by the model.
The Taiwan Strait is an important channel in the west Pacific Ocean transporting water and chemical constituents between the East China Sea and the South China Sea. Due to its complex bottom topography, alternating monsoon forcing and conjunction of several current systems [such as the Zhejiang-Fujian (Zhe-Min) Coastal Current, the Kuroshio intrusion and the extension of the South China Sea Warm Current], the physical and biogeochemical processes and ecosystem dynamics in the Taiwan Strait vary significantly both in space and in time. Our recent interdisciplinary studies, combining in situ and remote sensing observations with numerical modeling, allow us to address several important issues concerning the Taiwan Strait. The temporal and spatial variation of circulation in the Taiwan Strait ...
This paper proposes a high accuracy numerical method to model oil spill trajectories using a particle-tracking algorithm. The Euler method, used to calculate oil trajectories, can give adequate solutions in most open ocean applications. However, this method may not predict accurate particle trajectories in certain highly non-uniform velocity fields near coastal zones or in river problems. Simple numerical experiments show that the Euler method may also introduce artificial numerical dispersion that could lead to overestimation of spill areas. This article proposes a fourth-order Runge-Kutta method with fourth-order velocity interpolation to calculate oil trajectories that minimise these problems. The algorithm is implemented in the OilTrack model to predict oil trajectories following the 'Nissos Amorgos' oil spill accident that occurred in the Gulf of Venezuela in 1997. Despite lack of adequate field ...
... UGV), in the air (UAV), or underwater (UUV) -across ... Figure 2 Tactical Mobile Robot MOUT concept: ... 10 robots to clear 20 rooms in 60 minutes ...
A catalytic converter of a ceramic monolith honeycomb substrate, coated with a washcoat of catalyst and attached to a natural gas/diesel dual fuel engine was simulated and studied experimentally. The paper describes the application of one-dimensional finite element model for the transient and steady state operation. Laminar flow was approximated using a dispersed plug flow model, and chemical kinetics were simulated using LHHW (Langmuir/ Hinshelwood/ Hougan/ Watson) type expressions. Simulation results were compared with experimental results for heating and cooling cycles which resulted from speed and load changes on the engine. The comparison showed a maximum difference between the two sets of emission levels of about 10 per cent, showing that the one-dimensional model is acceptable model for this dual fuel engine converter combination. 50 ...
... that range from relatively low frequency strong ground motion modeling [Vidale et al., 1985] to high frequency reverse time migration of seismic ...
Stochastic models for the simulation of global radiation are discussed. Thermal transients in the ground are analyzed. The performance of buried-pipe storage and a space heating system with long-term storage is described.
Aug 18, 2004 ... The Arabian Gulf region presents a challenge to meteorologists trying to simulate weather with computer models because sea-surface and land ...
... which may have embedded remote terminals, and (4 ... model the time-variant or stochastic aspects of ... on such system parameters as queue size, time ...
This contribution presents different approaches for the modeling of gas entrainment under water by a plunging jet. Since the generation of bubbles happens on a scale which is smaller than the bubbles, this process cannot be resolved in meso-scale simulations, which include the full length of the jet and its environment. This is why the gas entrainment has to be modeled in meso-scale simulations. In the frame of a Euler-Euler simulation, the local morphology of the phases has to be considered in the drag model. For example, the gas is a continuous phase above the water level but bubbly below the water level. Various drag models are tested and their influence on the gas void fraction below the water level is discussed. The algebraic interface area density (AIAD) model applies a drag coefficient for bubbles and a different ...
The world's oceans contain a complex mixture of micro-organisms that are for the most part, uncharacterized both genetically and biochemically. We report here a metagenomic study of the marine planktonic...Full Text Available
Only half of the carbon dioxide produced by burning fossil fuels stays in the atmosphere. It is thought that the oceans absorb the rest. To understand the world's climate it is necessary to know how the atmosphere, ocean and biosphere interact.
The report gives an initial assessment of the feasibility of disposing of flue gas desulfurization wastes in mines and in the ocean. The purpose of the assessment was to evaluate environmental, technical, regulatory, and economic aspects of the use of such disposal sites. Availab...
The characteristics of martensitic transformations in solids are largely determined by the elastic strain that develops as martensite particles grow and interact. To study the development of microstructure, a finite-element computer simulationmodel was constructed to mimic the transformation process. The transformation is athermal and simulated at each incremental step by transforming the cell which maximizes the decrease in the free energy. To determine the free energy change, the elastic energy developed during martensite growth is calculated from the theory of linear elasticity for elastically homogeneous media, and updated as the transformation proceeds.
This project is about the simulation design of an engine control unit (ECU) for an Otto cycle engine with electronic fuel injection (EFI). The simulation includes a model for the ECU as well as physical parameters of the engine, which allows closed-loop control and monitoring of various systems.This simulation has been realized using Simulink and Stateflow, which are components of Mathworks' MATLAB software.The program allows control of various parameters of the ECU, as well as the simulation of failures to verify that the designed ECU is fault-tolerant and can control the engine using an open loop control. The main function the ECU provides is fuel metering.Subsequently, this program could be used as a tool to quickly develop and test models of ECU in order to control an engine in laborat...
... Studies People Projects Opportunities Framework Critical Sectors Development Goals Geographic Region Geographic Scale Research Themes Printer-Friendly Center for Ocean Solutions Early Career Fellowship Program Location: Stanford, California Source: The Center for Ocean Solutions (“Ocean Solutions”) seeks one or more recent graduates who have received a JD, MBA or PhD in the natural, physical or social sciences in the last five years, and who ...
''Climates of the Oceans'' deals with the atmosphere over the world oceans and provides a treatment not only of the climatic elements such as temperature, pressure, wind, precipitation etc., but also of the circulation of the atmosphere and its changes throughout the year. The connection between sea and overlying air is examined through studies of the water and heat balance. Numerous maps and tables illustrate the elements and their fluctuations over the water and on islands and coasts.
As part of a larger elastic numerical modeling project, we have been investigating how energy reflected from steeply dipping interfaces is recorded using typical multicomponent acquisition geometries. Specifically, we have been interpreting how rcflection events from the flanks of salt dome structures are distributed on 3C and 4C phones for vertical seismic profiles (VSPs) and ocean bottom seismic (OBS) or land surface surveys. The ultimate goal of this investigation is to improve the structural imaging of steeply dipping interfaces and eventually to evaluate the usc of the recorded elastic wavefield for fluid description near these interfaces. In the current work, we focus on a common assumption used when processing converted wave reflection seismic data that most PP energy is recorded on the vertical geophone and/or the hydrophone and that most PS energy is recorded on the horizontal geophones. This is a useful assumption when it is valid, ...
This paper overviews the general features of computer program 'HEC-5, Simulation of Flood Control and Conservation Systems', with emphasis on the capabilities of the most recent release of HEC-5, Version 7.2, dated March 1991. HEC-5 can simulate the essential features and operation goals and constraints of simple or complex systems with simulation intervals ranging from minutes to one month. Single event flood analysis and period of record conservation analysis may be accomplished with the model. Flood control analysis includes balanced system operation for downstream damage centers with consideration of forecasted local flows and hydrologic routing. In addition, induced surcharge operation based on spillway gate regulation schedules can be simulated. Hydropower analysis may include run-of-river, peaking, and pumped storage plants as well as system power operation. ...
A parallel bundle of transmembrane (TM) alpha-helices surrounding a central pore is present in several classes of ion channel, including the nicotinic acetylcholine receptor (nAChR). We have modeled...Full Text Available
Most waste management activities are decided upon and carried out in a public or semi-public arena, typically involving the waste management organization, one or more regulators, and often other stakeholders and members of the public. In these environments, simulationmodeling can be a powerful tool in reaching a consensus on the best path forward, but only if the models that are developed are understood and accepted by all of the parties involved. These requirements for understanding and acceptance of the models constrain the appropriate software and model development procedures that are employed. This paper discusses requirements for both simulation software and for the models that are developed using the software. Requirements for the software include transparency, accessibility, flexibility, extensibility, quality assurance, ability to ...
We introduce in this paper a new method for reducing neurodynamical data to an effective diffusion equation, either experimentally or using simulations of biophysically detailed models. The dimensionality...Full Text Available
Calculations, drawing principally on developments at AERE Harwell, of the relaxation about lattice defects are reviewed with emphasis on the techniques required for such calculations. The principles of defect modelling are outlined and various programs developed for defect simulations are discussed. Particular calculations for metals, ionic crystals and oxides, are considered. (UK).
The efficacy of beta-lactam antibiotics and amikacin alone and in various combinations against Pseudomonas aeruginosa was studied in a rabbit modelsimulating a closed-space infection in a locally neutropenic...Full Text Available
This research aims to develop a simulation approach based on system dynamics modelling (SDM) and adaptive network based fuzzy inference system (ANFIS) for quantifying and reducing the bullwhip effect in a multi-product, multi-stage supply chain. The proposed model is comprised of three groups of variables influencing the bullwhip effect, namely the structure of a supply chain network, supply chain contributions (ordering process in regular situation or when a supplier has a promotion or shortage gaming) and supply chain performances (the number of defects and ordering lead time). As a result, a two layer simulationmodel is developed with three generic models. The flexibility of this proposed approach is its ability to model various types of ordering policies which are basic inventory poli...
We have conducted an agent-based simulation of chain bankruptcy. The propagation of credit risk on a network, i.e., chain bankruptcy, is the key to nderstanding largesized bankruptcies. In our model, decrease of revenue by the loss of accounts payable is modeled by an interaction term, and bankruptcy is defined as a capital deficit. Model parameters were estimated using financial data for 1,077 listed Japanese firms. Simulations of chain bankruptcy on the real transaction network consisting of those 1,077 firms were made with the estimated model parameters. Given an initial bankrupt firm, a list of chain bankrupt firms was obtained. This model can be used to detect high-risk links in a transaction network, for the management of chain bankruptcy.
Although computer simulation has played a central role in the study of nucleation and growth since the earliest molecular dynamics simulations almost 50 years ago, confusion surrounding the effect of finite size on such simulations have limited their applicability. Modeling solidification in molten tantalum on the BlueGene/L computer, we report here on the first atomistic simulation of solidification that verifies independence from finite size effects during the entire nucleation and growth process, up to the onset of coarsening. We show that finite size scaling theory explains the observed maximal grain sizes for systems up to about 8,000,000 atoms. For larger simulations, a cross-over from finite size scaling to more physical size-independent behavior is observed.
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.
In this paper we report on multi-physics, multi-billion macroparticle simulation of beam transport in a free electron laser (FEL) linac for future light source applications. The simulation includes a self-consistent calculation of 3D space-charge effects, short-range geometry wakefields, longitudinal coherent synchrotron radiation (CSR) wakefields, and detailed modeling of RF acceleration and focusing. We discuss the need for and the challenges associated with such large-scale simulation. Applications to the study of the microbunching instability in an FEL linac are also presented.
Abstract The electrochemical properties of a perfluorosulfonic acid (PFSA) membrane are estimated using a combination of molecular dynamics simulation and statistical thermodynamic model. We obtain all parameters in an ionic conductivity model from an atomistic simulation and remove all adjusted model parameters. From a microscopic point of view, the hydrated PFSA membrane shows micro-phase segregation which separated into hydrophilic and hydrophobic phases. Our present work originates with this phenomenon and we treat this phase segregation as if it is a continuous phase for each of which the proton (H+) is transported inside the PFSA membrane/solvent (water and alcohols) mixture. The chemical potential for a given system is estimated using a molecular simulation technique to predict the ...
We have simulated the transient enhanced diffusion (TED) of boron fluoride (BF_2"+) implanted in crystalline and germanium amorphized silicon. Based on recently published models, the effect of fluorine on boron diffusion in silicon has been introduced and a modelling has been suggested. In order to simulate the boron experimental profiles, we have assumed that fluorine forms clusters involving interstitial boron which reduces the junction depth. Experimental results indicate that fluorine behaviour depends on amorphization energy. Moreover, even no germanium preamorphization is performed, silicon is still amorphized by fluorine species. Hence, BF_2"+ implantation leads to an amorphous/crystalline (a/c) interface near the surface. An improvement of published models is suggested taking into account fluorine effects. The simulations satisfactory reproduce the SIMS ...
This report describes the test of a dynamic phase compensation unit for a wind turbine with directly connected induction generators. The compensation unit is based on thyristor switched capacitors, where conventional wind turbine compensations use mechanical contactors to switch the capacitors. The influence on power quality analysed, and influence on component lifetime is discussed. Besides, simulationmodels in Matlab/Simulink are presented, including a flicker meter model. (au)
Performance predictions have been made for attached-sunspace types of passively solar heated buildings. The predictions are based on hour-by-hour computer simulations using computer models developed in the framework of PASOLE, the Los Alamos Scientific Laboratory (LASL) passive solar energy simulation program. The models have been validated by detailed comparison with actual hourly temperature measurements taken in attached-sunspace test rooms at LASL.
Taking both white noise and colored environmental noise into account, a general stochastic logistic model under regime switching is proposed and studied. Sufficient conditions for extinction, nonpersistence in the mean, weak persistence, stochastic permanence and global attractivity are established. The critical number between weak persistence and extinction is obtained. Moreover, some simulation figures are introduced to illustrate the main results.
The validation of the forecast model for early emergency response to nuclear accidents is evaluated by trace tests in atmosphere in Daya Bay nuclear power site. The simulation experiment of the Daya Bay nuclear power site shows that the particle spreading image and the time-integrated concentration distribution given by plume concentration prediction model can perform the variation of pathway of the pollutant transport, as well as the effects of topography on transport and diffusion of pollutants. The simulation of five trace tests in field shows that 59.1% of ratios between predicted results and observed results are within the range of 10, and 41% of ratios are within the range of 5 approximately. (authors)
In the framework of the carbon storage two forestry methods are analyzed by the author: the standing capitalization (slowing of the crops) and the spices substitution (slow growing species by fast growing species). The construction of a model taking into account the all stage of the carbon implication (biomass, soil, forest products) offers simulation and quantitative results on these methods. (A.L.B.)
The use of numerical modeling in oil spill incidents is a well established technique that has proven to provide cost-effective and reasonable estimates of oil surface drift. Good predictability of such models depends highly on the quality of the input data of the incident and on the model calibration effort. This paper presents the results of simulating oil spillage trajectory in the Arabian (Persian) Gulf. The study employed a 3-D rectilinear hydrodynamic model combined with oil spill model. Typical representative environmental conditions of the Arabian Gulf were first setup into a hydrodynamic circulation model using data from various sources. The performance of the hydrodynamic model was then tested against measurements of tidal fluctuation and sea currents at selected locations. The sp...
This paper aims to provide a comprehensive review of the various models or simulations for predicting drug release from bulk-degrading systems. A brief description of bulk degradation processes and factors affecting the degradation rate, and consequently the release kinetics, is presented first. Next, several important classical models, often used as the basis for subsequent model development, are discussed. Both mathematical models and Monte-Carlo based simulations have been developed for controlled release from bulk-degrading systems. The mathematical models can be further subdivided into two categories. First, the diffusion-based models whose transport mechanism is mainly governed by diffusion, but with degradation-dependent diffusion coefficients. These are generally simpler and easier...
The aim of this work is the evaluation of the quality of a chemistry-transport model, not by a classical comparison with observations, but by the estimation of its uncertainties due to the input data, to the model formulation and to the numerical approximations. The study of these 3 sources of uncertainty is carried out with Monte Carlo simulations, with multi-modelsimulations and with comparisons between numerical schemes, respectively. A high uncertainty is shown for ozone concentrations. To overcome the uncertainty-related limitations, a strategy consists in using the overall forecasting. By combining several models (up to 48) on the basis of past observations, forecasts can be significantly improved. This work has been also the occasion of developing an innovative modeling system, named Polyphemus. (J.S.)
As the surroundings of objects of oil exploration grow more complicated, seismic survey methods have turned 3-dimensional and, in this report, several models are examined using the 3-dimensional simulation technology. The result obtained by the conventional wave tracking method is different from actual wavefields, and is unrealistic. The difference method among the fullwave modelling methods demands an exorbitantly long computation time and high cost. A pseudospectral method has been developed which is superior to the difference method, and has been put to practical use thanks to the advent of parallel computers. It is found that a 3-dimensional survey is mandatory in describing faults. After examining the SEG/EAGE Salt model, it is learned that the salt is well-developed and that 3-dimensional depth migration is required for sub-salt exploration. It is also found through simulation ...
Current advanced CMOS source/drain engineering involves the use of amorphizing implants with 3D geometry. Upon annealing, the induced transient enhanced diffusion (TED) can only be accurately predicted if the amorphized region is correctly modeled, as well as the formation and evolution of extended defects, particularly 3 1 1's and dislocation loops. In addition to the extended defects, already modeled in the atomistic kinetic Monte-Carlo simulator DADOS, we have developed a physically based modeling approach for the implant-induced damage build-up, amorphization and recrystallization, suitable to handle device-size process simulation. It is based on amorphous pockets (3D, irregular shape agglomerates of an arbitrary number of interstitials and vacancies, plus trapped impurities) with a size-dependent activation energy for recombination. The model is ...
Current advanced CMOS source/drain engineering involves the use of amorphizing implants with 3D geometry. Upon annealing, the induced transient enhanced diffusion (TED) can only be accurately predicted if the amorphized region is correctly modeled, as well as the formation and evolution of extended defects, particularly 3 1 1's and dislocation loops. In addition to the extended defects, already modeled in the atomistic kinetic Monte-Carlo simulator DADOS, we have developed a physically based modeling approach for the implant-induced damage build-up, amorphization and recrystallization, suitable to handle device-size process simulation. It is based on amorphous pockets (3D, irregular shape agglomerates of an arbitrary number of interstitials and vacancies, plus trapped impurities) with a size-dependent activation energy for recombination. The model is able to ...
CommAspen is a new agent-based model for simulating the interdependent effects of market decisions and disruptions in the telecommunications infrastructure on other critical infrastructures in the U.S. economy such as banking and finance, and electric power. CommAspen extends and modifies the capabilities of Aspen-EE, an agent-based model previously developed by Sandia National Laboratories to analyze the interdependencies between the electric power system and other critical infrastructures. CommAspen has been tested on a series of scenarios in which the communications network has been disrupted, due to congestion and outages. Analysis of the scenario results indicates that communications networks simulated by the model behave as their counterparts do in the real world. Results also show that the model could be used to analyze the economic impact of ...
This document reproduces the final project of David San Fabian Ayuso, presented on May 26, 2005, for the obtention of the engineer degree of the Carlos III University of Madrid. A single-phase, isothermal model, including both electron and proton transport, is introduced for the simulation of polymer electrolyte fuel cells (PEM). The model is implemented in the commercial code Fluent 6.0, through the use of UDFs (User Defined Functions). In order to validate the model, a single canal of a PEM monocell is simulated in three dimensions. The obtained result are qualitatively satisfactory. It is observed that it is not essential to solve the current collectors when a monocell is considered (and not a stack). in the present study, the number of nodes is the computational grid appears to be too low in the membrane zone in order to make a complete validation of the ...
This report summarizes some experiences gathered when using the semi-empirical temperature model that forms the heart of the steel casting advisory system `TUNDTEMP` in simulating the casting process. The model estimates the evolution of the liquid steel temperature in a tundish of a continuous casting machine through a temperature prediction at a couple of central time instants during the casting process. The results are compared with corresponding actual observations during the castings - measured steel temperatures and recorded casting difficulties and interruptions. Explanations for the possible discrepancies are suggested. Starting from a heuristically chosen set of model parameter values, some developing steps have been tested to reduce the simulation errors. The study concentrates mainly on temperature predictions for single heat castings. At the end of the report, some ...
There is an ever increasing demand to perform flow simulations that incorporate the complete details of geometry as well as sophisticated flow physics. This has led to the development of numerical algorithms that can simulate the actual flow phenomena with greater fidelity. However, the success of these algorithms hinges on the grid that models the geometry. Grid generation methods for 2-D models have long existed and the general lack of complexity of the simpler 2-D models has not quite challenged the efforts in this area. However, demands for generating better 3-D geometric models for flow simulations involving complex geometries have completely changed the perspective of grid generation strategies. As a consequence, grid generation efforts have earned equal significance as that of numerical solver efforts.
During the Eddy Experiment, two synchronous GPS receivers were flown at 1 km altitude to collect L1 signals and their reflections from the sea surface for assessment of altimetric precision and accuracy. Wind speed (U10) was around 10 m/s, and SWH up to 2 m. A geophysical parametric waveform model was used for retracking and estimation of the lapse between the direct and reflected signals with a 1-second precision of 3 m. The lapse was used to estimate the SSH along the track using a differential model. The RMS error of the 20 km averaged GNSS-R absolute altimetric solution with respect to Jason-1 SSH and a GPS buoy measurement was of 10 cm, with a 2 cm mean difference. Multipath and retracking parameter sensitivity due to the low altitude are suspected to have degraded accuracy. This result provides an important milestone on the road to a GNSS-R mesoscale altimetry space mission.
The precipitation of excess silicon interstitials into dislocation loops is modeled. This situation occurs when an amorphous layer is created at the surface in order to avoid boron channeling and form shallow p junctions. The modeling of the nucleation of these extended defects is included into the process simulator IMPACT-4. Their density and mean radius are calculated for several annealing times and temperatures and they are compared with experimental characterizations. This is the first step towards a full modeling of the complex processes involved in the transient enhanced diffusion of boron.
The plasma plume generated by a new type of high energy Janitor known as the railplug, is examined. The railplug is a miniaturized railgun that has the potential for improving ignition characteristics of combustible mixtures in engines. The objective of the study is to gain an uderstanding of the characteristics of the plasma created by a transparent railplug, and to validate a multidimensional computer simulation of the plasma and shock fronts. The nature of the plume emitted by the railplug was examined for three levels of electrical energy while firing into air at a pressure of 1 atm. The computer model is to be used to predict trends in railplug performance for various railplug designs, energies, and ambient conditions. The velocity of the plasma movement inside a transparent railplug was measured, as well as the velocity of the plume ejected from the cavity. A shock is produced at the initiation point of the arc and propagates down the ...
The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. In this paper, we test the present theoretical and observational models via a relativistic numerical simulation of the jets in the radio galaxy 3C 31. We use the parameters derived from the modelling presented by \\cite{lb02a,lb02b} as input parameters for the simulation of the evolution of the source, thus assuming that they have not varied over the lifetime of the source. We simulate about 10 % of the total lifetime of the jets in 3C 31. Realistic density and pressure gradients for the atmosphere are used. The simulation includes an equation of state for a two-component relativistic gas that allows a separate treatment of leptonic and baryonic matter. We compare our results with the modelling of the observational data of the source. Our results show ...
Accurate modeling of system response and scatter distribution is crucial for image reconstruction in emission tomography. Monte Carlo simulations are very well suited to calculate these quantities. However, Monte Carlo simulations are also slow and many simulated counts are needed to provide a sufficiently exact estimate of the detection probabilities. In order to overcome these problems, we propose to split the simulation into two parts, the detection system and the object to be imaged (the patient). A so-called 'virtual boundary' that separates these two parts is introduced. Within the patient, particles are simulated conventionally. Whenever a photon reaches the virtual boundary, its detection probability is calculated analytically by evaluating a multi-dimensional B-spline that depends on the photon position, direction and energy. The unknown B-spline knot ...
In accordance with the Nuclear Regulatory Commission regulation regarding groundwater travel times at geologic repositories, various models of unsaturated flow in fractured tuff have been developed and implemented to assess groundwater travel times at the potential repository at Yucca Mountain, Nevada. Kaplan used one-dimensional models to describe the uncertainty and sensitivity of travel times to various processes at Yucca Mountain. Robey and Arnold et al. used a two-dimensional equivalent continuum model (ECM) with inter- and intra-unit heterogeneity in an attempt to assess fast-flow paths through the unsaturated, fractured tuff at Yucca Mountain (GWTT-94). However, significant flow through the fractures in previous models was not simulated due to the characteristics of the ECM, which requires the matrix to be nearly saturated before flow through the fractures is initiated. In ...
The uptake of sup(99m)Tc-methylenediphosphonate (MDP) in different parts of rat femur was simulated using a local three-space model for tracer transfer. The model consisted of bone blood, bone ECF-space and space for tracer deposition. The measured sup(99m)Tc-MDP concentration in the systemic blood and the local bone blood flow measured by /sup 131/I-macroaggregated albumin microspheres were used as input parameters. The measured blood flow values were 6.3, 3.1 and 15.3 ml/100 g/min for proximal, middle and distal femur, respectively. the model parameters that gave the best fit to measured sup(99m)Tc-MDP uptake curves in computer simulation showed that bone blood flow, volume of ECF-space, permeability surface area product and accretion constant from ECF-space to space for tracer deposition were highest in distal and lowest in middle femur. The values corresponded to peak extraction ...
The uptake of sup(99m)Tc-methylenediphosphonate (MDP) in different parts of rat femur was simulated using a local three-space model for tracer transfer. The model consisted of bone blood, bone ECF-space and space for tracer deposition. The measured sup(99m)Tc-MDP concentration in the systemic blood and the local bone blood flow measured by "1"3"1I-macroaggregated albumin microspheres were used as input parameters. The measured blood flow values were 6.3, 3.1 and 15.3 ml/100 g/min for proximal, middle and distal femur, respectively. the model parameters that gave the best fit to measured sup(99m)Tc-MDP uptake curves in computer simulation showed that bone blood flow, volume of ECF-space, permeability surface area product and accretion constant from ECF-space to space for tracer deposition were highest in distal and lowest in middle femur. The values corresponded to peak extraction ...
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 ...
This paper proposes the use of particle swarm optimization method (PSO) for finite element (FE) model updating. The PSO method is compared to the existing methods that use simulated annealing (SA) or genetic algorithms (GA) for FE model for model updating. The proposed method is tested on an unsymmetrical H-shaped structure. It is observed that the proposed method gives updated natural frequencies the most accurate and followed by those given by an updated model that was obtained using the GA and a full FE model. It is also observed that the proposed method gives updated mode shapes that are best correlated to the measured ones, followed by those given by an updated model that was obtained using the SA and a full FE model. Furthermore, it is observed that the PSO achieves this accuracy at a computational speed that is ...
A constitutive model for an explosive consists of three parts: equation of state of reactants, equation of state of products, and reaction rate. In the absence of a good physical theory, models are empirical and based on heuristics. For the HMX based plastic-bonded explosive PBX 9501, three models have previously been calibrated and used in numerical simulations. Here, these three models are compared. They illustrate several generic difficulties with developing and validating models. The available data is insufficient to determine a constitutive model uniquely. Different models can be calibrated to fit the available data yet differ significantly in regions of phase space of importance to applications. Comparing models gives an indication of the uncertainty in model predictions ...
In this paper we describe a modelling project to improve a nuclear waste management program in charge of the creation of a new system for the permanent disposal of nuclear waste. SADT (Structural Analysis and Design Technique) is used in order to provide a work-flow description of the functions to be performed by the waste management program. This description is then translated into a number of Coloured Petri Nets (CPN or CP-nets) corresponding to different program functions where additional behavioural inscriptions provide basis for simulation. Each of these CP-nets is simulated to produce timed event charts that are useful for understanding the behaviour of the program functions under different scenarios. Then all the CPN models are linked together to form a single stand-alone application that is useful for validating the interaction and cooperation between the different program functions. A technique ...
As usual a shipping terminal contains a header and its diameter is smaller than the main pipeline diameter. This diameter reduction amplifies the effects caused by hydraulic transients. It was noticed during simulations that check valves without dynamic modeling may introduce some error in the maximum pressure results without a dynamic model approach. The current paper uses commercial pipeline simulation software to model the dynamic behavior of the check valves. It was studied the header diameter influence and how the check valve model type may change the maximum pressure in the pipeline. It was proved that even a pipeline that works with low pressure and small elevation drop needs a detailed valve modeling to prevent the calculation of unreal pressure values in the region with diameter reduction. (author)
Stochastic simulation has been employed in petroleum reservoir characterization as a modeling tool able to reconcile information from several different sources. It has the ability to preserve the variability of the modeled phenomena and permits transference of geological knowledge to numerical models of flux, whose predictions on reservoir constitute the main basis for reservoir management decisions. Several stochastic models have been used and/or suggested, depending on the nature of the phenomena to be described. Markov Random Fields (MRFs) appear as an alternative for the modeling of discrete variables, mainly reservoirs with mosaic architecture of facies. In this dissertation, the reader is introduced to the stochastic modeling by MRFs in a generic sense. The main aspects of the technique are reviewed. MRF Conceptual Background is ...
We develop a 3-D model to simulate the synthetic aperture radar (SAR) image formation process of an undulated vegetation canopy such as corn grown in fields with large periodic drainage reliefs. We explain how the simulated SAR image of undulated vegetation medium is obtained by the convolution of a 2-D slice of the 3-D simulated SAR system point spread function [(PSF), emulating the SAR beam modeled by a cosine modulated Gaussian], with the 2-D projection of the observed undulated vegetation canopy (modeled with scatterers randomly distributed in 3-D undulated space) followed by the extraction of each look envelope, the summation of looks, and sampling in azimuthal and range directions. Our model is useful to study the parameters involved in the formation and the analysis of SAR images of undulated vegetation medium. Validation of ...
In a Pb-Bi cooled direct contact steam generation fast reactor water is injected directly above the core, the produced steam is separated at the top and is send to the turbine. Neither the direct contact phenomenon nor the two-phase flow simulations in CFD have been thoroughly described yet. A first attempt in simulating such two-phase flow in 2D using the CFD code Fluent is presented in this paper. The volume of fluid explicit model was used. Other important simulation parameters were: pressure velocity relation PISO, discretization scheme body force weighted for pressure, second order upwind for momentum and CISCAM for void fraction. Boundary conditions were mass flow inlet (Pb-Bi 0 kg/s and steam 0.07 kg/s) and pressure outlet. The effect of mesh size (0.5 mm and 0.2 mm cells) was investigated as well as the effect of the turbulent model. It was found that using a fine mesh is ...
There is a model on the Oak Ridge Reservation which provides requirements for determining capacities and capabilities related to low-level, hazardous, and mixed wastes. In FY 1987, the model will be sufficiently advanced to provide various waste management scenarios. These scenarios will be compared technically, operationally, and financially by use of waste characterization data and process simulators that are currently under development. The results of the process simulations will be used to help identify waste treatment, storage, and disposal technologies that need to be demonstrated prior to full-scale development for DOE use. The information derived from this effort will be made available to all DOE facilities.
With a modelling approach based on cellular automata, five observed types of plant development can be simulated. In addition, the proposed model shows a strong tendency towards the formation of patches and a high degree of dynamical and structural instability leading to limits of predictability for the asymptotic solution chosen by the system among several possible metastable patterns (multistability). Further, external fluctuations can be shown to have advantages for certain plant types. The presented model unifies the fundamental dichotomy in vegetation dynamics between determinism (understood as predictability) and disorder (chance effects) by showing the outcome of both classical theories as special cases. (author) 2 figs., 4 refs.
Inductive electromagnetic launchers, or coilguns, use discrete solenoidal coils to accelerate a coaxial conductive armature. To date, Sandia has been using an internally developed code, SLINGSHOT, as a point-mass lumped circuit element simulation tool for modeling coilgun behavior for design and verification purposes. This code has shortcomings in terms of accurately modeling gun performance under stressful electromagnetic propulsion environments. To correct for these limitations, it was decided to attempt to closely couple two Sandia simulation codes, Xyce and ALEGRA, to develop a more rigorous simulation capability for demanding launch applications. This report summarizes the modifications made to each respective code and the path forward to completing interfacing between them.
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.
We review our recent work on an atomistic approach to the development of predictive process simulation tools. First principles methods, molecular dynamics simulations, and experimental results are used to construct a database of defect and dopant energetics in Si. This is used as input for kinetic Monte Carlo simulations. C and B trapping of the Si self- interstitial is shown to help explain the enormous disparity in its measured diffusivity. Excellent agreement is found between experiments and simulations of transient enhanced diffusion following 20-80 keV B implants into Si, and with those of 50 keV Si implants into complex B-doped structures. Our simulations predict novel behavior of the time evolution of the electrically active B fraction during annealing.
This paper presents a trial of probabilistic simulation for performance assessment of high-level waste (HLW) disposal using the Monte Carlo method. Based on the Reference Case conceptual model in H12, a new integrated simulation system which allowed us to rapidly evaluate the effect of data uncertainty was developed. The doses to hypothetical exposure group were compared with the results of H12 that were performed by using a point-wise approach, in terms of maximum total dose. This study showed that H12 results were consistent with results of probabilistic simulation and also showed that transmissivity had a strong influence for the uncertainty of the system performance in all simulation time. (author)
The traditional process for developing large-scale simulations is cumbersome, time consuming, costly, and in some cases, inadequate. The topics of software components and component-based software engineering are being explored by software professionals in academic and industrial settings. A component is a well-delineated, relatively independent, and replaceable part of a software system that performs a specific function. Many researchers have addressed the potential to derive a component-based approach to simulations in general, and a few have focused on military simulations in particular. In a component-based approach, functional or logical blocks of the simulation entities are represented as coherent collections of components satisfying explicitly defined interface requirements. A simulation is a top-level aggregate comprised of a collection of components that interact with each ...
The implementation of wind energy conversion systems in the built environment renewed the interest and the research on Vertical Axis Wind Turbines (VAWT), which in this application present several advantages over Horizontal Axis Wind Turbines (HAWT). The VAWT has an inherent unsteady aerodynamic behavior due to the variation of angle of attack with the angle of rotation, perceived velocity and consequentially Reynolds number. The phenomenon of dynamic stall is then an intrinsic effect of the operation of a Vertical Axis Wind Turbine at low tip speed ratios, having a significant impact in both loads and power. The complexity of the unsteady aerodynamics of the VAWT makes it extremely attractive to be analyzed using Computational Fluid Dynamics (CFD) models, where an approximation of the continuity and momentum equations of the Navier-Stokes equations set is solved. The complexity of the problem and the need for new design approaches for VAWT for the built ...
A computational fluid dynamics (CFD) is coupled with a computational structural dynamics (CSD) to simulate the unsteady rotor flow with aeroelasticity effects. An unstructured upwind Navier-Stokes solver was developed for this simulation, with 2nd order time-accurate dual-time stepping method for temporal discretization and low Mach number preconditioning method. For turbulent flows, both the Spalart-Allmaras and Menter's SST model are available. Mesh deformation is achieved through a fast dynamic grid method called Delaunay graph map method for unsteady flow simulation. The rotor blades are modeled as Hodges & Dowell's nonlinear beams coupled flap-lag-torsion. The rotorcraft computational structural dynamics code employs the 15-dof beam finite element formulation for modeling. The structure code was validated by comparing the natural frequencies of a rotor ...
A flowing bed kiln is a gas-solid reactor used in the civil nuclear fuel cycle for the successive conversion of uranium trioxide (UO{sub 3}) into uranium dioxide (UO{sub 2}) and then into uranium tetrafluoride (UF{sub 4}). A numerical model is developed which simulate the behaviour of this reactor in permanent regime. This model describes the physico-chemical phenomena involved, and combines a mechanistic approach in the vertical area of the kiln (resolution by the finite volumes method) and a systemic approach in the horizontal area, like in the model of cascade mixers. The first results have been obtained for reference operating conditions of the industrial kiln. Some possible improvements of the optimum temperature progression inside the kiln are evoked. (J.S.)
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 ...
The importance of temperature in the determination of the yield of an annual crop (groundnut; Arachis hypogaea L. in India) was assessed. Simulations from a regional climate model (PRECIS) were used with a crop model (GLAM) to examine crop growth under simulated current (1961-1990) and future (2071-2100) climates. Two processes were examined: the response of crop duration to mean temperature and the response of seed-set to extremes of temperature. The relative importance of, and interaction between, these two processes was examined for a number of genotypic characteristics, which were represented by using different values of crop model parameters derived from experiments.The impact of mean and extreme temperatures varied geographically, and depended upon the simulated genotypic properties....
A new numerical technique called the convolution-based particle tracking (CBPT) method is developed to simulate resident or flux-averaged solute concentrations in groundwater models. The method is valid for steady-state flow and linear transport processes such as sorption with a linear sorption isotherm and first-order decay. The CBPT method uses particle tracking to take advantage of the ability of particle-based approaches to maintain sharp fronts for advection-dominated transport problems common in groundwater modeling and because of the flexibility of the random walk method to simulate a wide range of possible forms of the dispersion tensor. Furthermore, the algorithm for carrying out the convolution and superposition calculation from particle tracking results is very efficient. We sho...
Ion implantation and subsequent annealing are essential stages in today's advanced CMOS processing. Although the dopant implanted profile can be accurately predicted by analytical fits calibrated with SIMS profiles, the damage has to be estimated with a binary collision approximation implant simulator. Some models have been proposed, like the '+n', in an attempt to simplify the anneal simulation. We have used the atomistic kinetic Monte Carlo dados to elucidate which are the implant modeling features most relevant in the simulation of transient enhanced diffusion (TED). For the experimental conditions studied we find that the spatial correlation of the I, V Frenkel pairs is not critical in order to yield the correct I supersaturation, that can be simulated just taking into account the net I-V excess distribution. In contrast to, simulate ...
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 ...
A unified physically based ion implantation damage model has been developed which successfully predicts both the impurity profiles and the damage profiles for a wide range of implant conditions for arsenic, phosphorus, BF{sub 2}, and boron implants into single-crystal silicon. In addition, the amorphous layer thicknesses predicted by this new damage model are also in excellent agreement with experimental measurements. This damage model is based on the physics of point defects in silicon, and explicitly simulates the defect production, diffusion, and their interactions which include interstitial-vacancy recombination, clustering of same type of defects, defect-impurity complex formation, emission of mobile defects from clusters, and surface effects for the first time. New computationally efficient algorithms have been developed to overcome the barrier of the excessive computational requirements. In ...
A unified physically based ion implantation damage model has been developed which successfully predicts both the impurity profiles and the damage profiles for a wide range of implant conditions for arsenic, phosphorus, BF_2, and boron implants into single-crystal silicon. In addition, the amorphous layer thicknesses predicted by this new damage model are also in excellent agreement with experimental measurements. This damage model is based on the physics of point defects in silicon, and explicitly simulates the defect production, diffusion, and their interactions which include interstitial-vacancy recombination, clustering of same type of defects, defect-impurity complex formation, emission of mobile defects from clusters, and surface effects for the first time. New computationally efficient algorithms have been developed to overcome the barrier of the excessive computational requirements. In addition, ...
The shock stand-off distance of a spherical model flying with transonic speeds is determined through numerical simulations. The model decelerates due to drag forces caused by the pressure and viscous shear stress at the model surface. Two-dimensional axisymmetric numerical codes with numerical grids fixed to the flying spherical model are used in the simulations. Numerically determined shock stand-off distances are compared with experimental data obtained in a previous study as well as with those obtained in our ballistic-range experiments. The numerical results and the experimental data are found to be in good agreement. In addition, the time-dependent shock stand-off distance of a decelerating model is investigated.
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...
An electrochemistry model was developed to analyse the J-V characteristics of a Proton Exchange Membrane (PEM) water electrolyzer for hydrogen production. The Butler-Volmer equation and water transport characteristics through electrolyte membrane were employed to simulate the electrode activation over-potential and membrane ohmic over-potential, respectively. The modeling results are found to agree reasonably well with experimental data published in the literature. The parametric simulations show that the ohmic over-potential is relatively small with typical water content in the membrane. Compared with the cathode over-potential, the anode over-potential is more significant and constitutes the major source of voltage loss. The high anode over-potential is due to the relatively slow oxidation kinetics, which is related to anode material property and microstructure. This model can be ...
The study abstracted aims at analyzing the economic effects of energy policies using simulations. Emphasis is on assessing the effects of decisions in favour of an expanding electric power industry based on nuclear energy and hard coal. The study discusses alternative assumptions concerning the cost involved, pollution abatement measures, shutdown and dismantling costs etc. Input-output analyses make the simulationmodel a consistent vehicle. Besides, the model in question is characterized by its explaining the investment demand of the different sectors by referring to the continuous annual development of innovative technologies.
We investigate the formation of nanostructures in 2D strained alloys on face centered cubic (111) surfaces by means of equilibrium Monte Carlo simulations. In the framework of an off-lattice model, we consider one monolayer of two bulk-immiscible adsorbates A and B with negative and positive misfit relative to the substrate, respectively. Simulations show that the adsorbates partly self-organize into island or stripe-like patterns. We show how these structures depend on the relative misfits, interaction, and concentration of components. The morphology is quite different for phase separation and intermixing regimes.
The solvation properties of a chlorine ion in small water clusters are investigated using state-of-the-art statistical mechanics. The simulations employ the polarizable water model developed recently by Dang [J. Chem. Phys. [bold 97], 2659 (1992)]. The ion--water interaction potentials are defined such that the successive binding energies for the ionic clusters, and the solvation enthalpy, bulk vertical binding energy, and structural properties of the aqueous solution agree with the best available results obtained from experiments. Simulated vertical electron binding energies of the ionic clusters Cl[sup [minus
The hot uniform aluminium plasma was produced by irradiating thin aluminium dotted foil smoothly with the 9th 0.53 ?m laser on Shenguang II laser facility. The emitted spectrum was measured from the front and tangential direction of the target with two crystal spectrometers, and the quantitative spectrum from the front of the target was obtained. The state of laser- produced plasma was simulated with the radiation hydrodynamics code MULTI-1D, and the emitted spectrum was calculated with the spectrum code of Collision-Radiation model under the simulated plasma state. The experimental spectrum accords with the simulated one. (authors)
The influence of the velocity profile at the inlet boundary on the simulation of air velocity distribution inside an electrostatic precipitator is presented in this study. Measurements and simulations were performed in a duct and an electrostatic precipitator (ESP). A four-hole cobra probe was used for the measurement of velocity distribution. The flow simulation was performed by using the computational fluid dynamics (CFD) code FLUENT. Numerical calculations for the air flow were carried out by solving the Reynolds-averaged Navier-Stokes equations coupled with the realizable k-{epsilon} turbulence model equations. Simulations were performed with two different velocity profiles at the inlet boundary - one with a uniform (ideal) velocity profile and the other with a non-uniform (real) velocity profile to demonstrate the effect of velocity inlet boundary condition on the flow ...
Maps of wind resources were used to develop a method of evaluating climatic changes and their potential impacts on wind energy resources. Global IPCC climate change scenarios were used to predict climatic conditions for the future, while past wind resource availability was simulated and validated using NCEP and NCAR reanalysis data as well as observed meteorological data from Environment Canada. The simulations were used to compare each 5-year period with a 50-year reference period. Regional scale climate change impacts were evaluated using a statistical dynamic down-scaling method. Advanced meteorological models were used to predict wind flow patterns across specific landscapes. The evolution of past wind resource availability was then simulated. Five-year wind resource simulations for a 50-year period were simulated at 25 km"2 wind speeds at 80 m above the ...
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 ...
We investigate the influence of lunar-like satellites on the infrared orbital light curves of Earth-analog extra-solar planets. Such light curves will be obtained by NASA's Terrestrial Planet Finder (TPF) and ESA's Darwin missions as a consequence of repeat observations to confirm the companion status of a putative planet. We use an energy balance model to calculate disk-averaged infrared (bolometric) fluxes from planet-satellite systems over a full orbital period (one year). The satellites are assumed to lack an atmosphere, have a low thermal inertia like that of the Moon and span a range of plausible radii. The planets are assumed to have thermal and orbital properties that mimic those of the Earth while their obliquities and orbital longitudes of inferior conjunction remain free parameters. Even if the gross thermal properties of the planet can be independently constrained (e.g. via spectroscopy or visible-wavelength detection of specular glint from a surface ...
The Atmospheric Radiation Monitor (ARM) uninhabited aerospace vehicle (UAV) program aims to develop measurement techniques and instruments suitable for a new class of high altitude, long endurance UAVs while supporting the climate community with valuable data sets. Using the Scaled Composites Proteus aircraft, ARM UAV participated in Mixed-Phase Arctic Cloud Experiment (M-PACE), obtaining unique data to help understand the interaction of clouds with solar and infrared radiation. Many measurements obtained using the Proteus were coincident with in-situ observations made by the UND Citation. Data from M-PACE are needed to understand interactions between clouds, the atmosphere and ocean in the Arctic, critical interactions given large-scale models suggest enhanced warming compared to lower latitudes is occurring.
The National Survey and Cadastre - Denmark (KMS) has for several years produced gravity anomaly maps over the oceans derived from satellite altimetry. During the last four years, KMS has also conducted airborne gravity surveys along the coast of Greenland dedicated to complement the existing onshore gravity coverage and fill in new data in the very-near coastal area, where altimetry data may contain gross errors. The airborne surveys extend from the coastline to approximately 100 km offshore, along 6000 km of coastline. An adequate merging of these different data sources is important for the use of gravity data especially, when computing geoid models in coastal regions. The presence of reliable marine gravity data for independent control offers an opportunity to study procedures for the merging of airborne and satellite data around Greenland. Two different merging techniques, both based on collocation, are investigated in this paper. ...
This study determined how well crude and fuel oils burn on water. Objectives were: (1) to measure the burning rates for several oils; (2) to determine whether adding heat improves the oils' combustibility; (3) to identify the conditions necessary to ignite fuels known to be difficult to ignite on ocean waters (e.g., diesel and Bunker C fuel oils); and (4) to evaluate the accuracy of an oil-burning model proposed by Thompson, Dawson, and Goodier (1979). Observations were made about how weathering and the thickness of the oil layer affect the combustion of crude and fuel oils. Nine oils commonly transported on the world's major waterways were tested. Burns were first conducted in Oklahoma under warm-weather conditions (approx. 30/sup 0/C) and later in Ohio under cold-weather conditions (approx. 0/sup 0/C to 10/sup 0/C).
A modified variable structure interacting multiple model (M-VSIMM) estimator for complex hybrid maneuver target tracking is presented. The M-VSIMM could potentially be applied to fire control systems (FCS) used on warships. Target model groups were designed using 3D dimensional dynamic target models. Optimal model group selection logic was proposed, contrary to the activation and termination logic in the original VSIMM. The system will respond faster with optimal model group selection logic. After performing simulations, the tracking performances of the Kalman, ?-?(-?), VDIE, IMM and M-VSIMM filters were compared under various maneuvering conditions.
This paper presents a numerical simulation with the EUROPLEXUS fast dynamics software of a pipe whipping phenomenon occurring in the thermal hydraulic conditions of a loss of coolant accident in a PWR primary circuit. Different physical phenomena take place simultaneously during the rupture and the whipping of the pipe such as plasticity, contact, large displacements, two-phase flow regime and fluid structure interaction. Two kinds of numerical models - a simplified pipeline model and a mixed 1D/3D model - are considered and compared throughout modelling and computation. Numerical results are compared with experimental data belonging to the Aquitaine II test campaign.
Recent developments in the physical model of 1 MeV to 100 TeV hadron and lepton interactions with nuclei and atoms are described. These include a new nuclear cross section library, a model for soft pion production, the cascade-exciton model, the dual parton model, deuteron-nucleus and neutrino-nucleus interaction models, detailed description of mu, pi and anti p absorption and a unified treatment of muon and charged hadron electromagnetic interactions with matter. New algorithms are implemented into the MARS13(98) Monte Carlo code and benchmarked against experimental data. The code capabilities to simulate cascades and generate a variety of results in complex media have been also enhanced.
This paper studies on the modelling of a butterfly-type control valve actuating by an on-off pneumatic solenoid valve. The mathematical model is composed of nonlinear differential equations three parts: (i) a solenoid valve, (ii) a pneumatic cylinder, (iii) a rotary-type butterfly valve. The flow characteristics of the butterfly control valve is analysed by a computer simulator, then its simple transfer function is identified from the step responses.
ObjectivesNot EnteredDescriptionMethane hydrate in an ice-like substance consisting of molecules of methane gas combined chemically with water. It is stable at high pressures and low temperatures. Since the ocean floor is normally cold, but temperatures increase with depth inside the Earth, such conditions normally prevail for a few hundred meters below the seafloor where the ocean depth is more than a few hundred metres. Methane beneath the ocean floor is formed by the decomposition of organic material. At the edges of the co [continued...
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.
Gas and stars in spiral galaxies are modelled with the DUAL code, using hydrodynamic and N-body techniques. The simulations reveal morphological differences mirroring the dual morphologies seen in B and K' band observations of many spiral galaxies. In particular, the gaseous images are more flocculent with lower pitch angles than the stellar images, and the stellar arm-interarm contrast correlates with the degree of morphological decoupling.
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)
The paper deals with a model in which the motion of a gas in plasma accerators and high-curent discharges in the present of a skin effect is treated as expulsion by a 'magnetic piston' under the action of surface current flow. Specifically examined is the situation where the initial gas pressure is negligible in comparison with the magnetic pressure, and the motion of the gas may be treated as a self-simulating one. A system of hydrodynamic equations of the problem is derived and integrated by Adam's method. Results are plotted and discussed.
A physical model of free-electron laser (FEL) amplifier with variable-parameter wiggler magnets for one-dimensional numerical simulation is presented and a numerical example is given. The wiggler parameters, efficiency of energy conversion between electron beam and laser field, laser intensity, phase-space distributions and energy spectrum of electrons are computed. The period of synchronous oscillation and saturation value of laser intensity agree with estimated one.
The Monte Carlo simulation method is applied to calculate the nucleon transfer part of the imaginary optical-model potential in colliding heavy ions. In the process of transfer of nucleons from one nucleus to the other, forward and backward neutrons and protons are treated separately by taking into consideration the effect of the Coulomb forces and the potential energy surfaces (driving forces). The Pauli exclusion principle is also incorporated properly in the calculation.
Power devices such as MOSFETSs and IGBTs, include parasitic structures that can give rise to destructive failures such as breakdown and latch-up. To determine a suitable strategy for device radiation hardening, simulation software like MEDICI-2D can be used to model the effects of technological modifications and device parameters that are difficult to measure experimentally. (authors).
Chemical composition of fine aerosol particles over the northern Chilean coastal waters was determined onboard the U.S. DOE G-1 aircraft during the VOCALS (VAMOS Ocean-Cloud-Atmosphere-Land Study) field campaign between October 16 and November 15, 2008. SO42-, NO3-, NH4+, and total organics (Org) were determined using an Aerodyne Aerosol Mass Spectrometer, and SO42-, NO3-, NH4+, Na+, Cl-, CH3SO3-, Mg2+, Ca2+, and K+ were determined using a particle-into-liquid sampler-ion chromatography technique. The results show the marine boundary layer (MBL) aerosol mass was dominated by non- sea-salt SO42- followed by Na+, Cl-, Org, NO3-, and NH4+, in decreasing importance; CH3SO3-, Ca2+, and K+ rarely exceeded their respective limits of detection. The SO42- aerosols were strongly acidic as the equivalent NH4+ to SO42- ratio was only ~0.25 on average. NaCl particles, presumably of sea-salt origin, showed chloride deficits but retained Cl- typically more than half the ...
The presented paper describes the formulation of the fractional mass balances which are used to calculate the concentration and size distribution of char in the balance cells of a circulating fluidized bed boiler. In case of using the shrinking core model additionally the distribution of conversion of the char particles is calculated. Basically two reaction models (a shrinking particle and a shrinking core model) are available to describe the combustion process. The model considers all essential physico/chemical phenomena (primary fragmentation, attrition, mass transfer to gas phase) which have to be taken into account for the description of the combustion process. The model is part of an overall simulation program for circulating fluidized bed boilers. This engineering simulator has been developed to return the main performance parameters ...
An essential part of fire risk assessment is the analysis of fire hazards and fire propagation. In this work, models and tools for two different aspects of numerical fire simulation have been developed. The primary objectives have been firstly to investigate the possibility of exploiting state-of-the-art fire models within probabilistic fire risk assessments and secondly to develop a computationally efficient solver of thermal radiation for the Fire Dynamics Simulator (FDS) code. In the first part of the work, an engineering tool for probabilistic fire risk assessment has been developed. The tool can be used to perform Monte Carlo simulations of fires and is called the Probabilistic Fire Simulator (PFS). In Monte Carlo simulation, the simulations are repeated multiple times, covering the whole range of variability of the ...
The basic relationships of the convolution/superposition dose calculation technique are reviewed, and a modelling technique that can be used for obtaining a satisfactory beam model for a commercially available convolution/superposition-based treatment planning system is described. A fluence energy spectrum for a Co-60 treatment machine obtained from a Monte Carlo simulation was used for modelling the fluence spectrum for a Rokus M machine. In order to achieve this model we measured the depth dose distribution and the dose profiles with a Wellhofer dosimetry system. The primary fluence was iteratively modelled by comparing the computed depth dose curves and beam profiles with the depth dose curves and crossbeam profiles measured in a water phantom. The objective of beam modelling is to build a model of the primary fluence ...
The aim of this work is the development and experimental validation of a detailed dynamic fuel cell model using the gPROMS modeling environment. The model is oriented towards optimization and control and it relies on material and energy balances as well as electrochemical equations including semi-empirical equations. For the experimental validation of the model a fully automated and integrated hydrogen fuel cell testing unit was used. The predictive power of the model has been compared with the data obtained during load change experiments. A sensitivity analysis has been employed to reveal the most critical empirical model parameters that should be estimated using a systematic estimation procedure. Model predictions are in good agreement with experimental data under a wide range of operati...
Federal regulators were set to vote on a plan to protect deep water corals and other sensitive fish habitats that will likely include a permanent ban on ...
The classification of petrels (Pterodroma spp.) from Round Island, near Mauritius in the Indian Ocean, has confounded researchers since their...Full Text Available
Deforestation farther inland causes delayed reactions. When roots no longer hold soil in place, it ends up muddying streams, rivers, and coastal waters. ...
A Self-Powered Underwater Robot for Ocean Exploration and Beyond ... first unmanned underwater vehicle (UUV) that is completely powered by renewable energy . ...
The company reports that with the help of outside companies, they completed all major fixes including high-voltage transmission lines to neighborhood ...
A strategy for the management of the disposal of low-level radioactive wastes into coastal zones and ocean waters has developed over the past three decades. While there has been a substantial increase in the number of international and required agreements there has also been a concomitant improvement in our understanding of the ocean and the processes at work within it. This has allowed more rigorous and reliable assessments to be made of the radiological consequences. With the continued development of basic guidelines of radiological protection by the International Commission on Radiological Protection (ICRP) methodologies have been formulated to derive the fundamental scientific requirement - the relationship between disposal or release rate to the ocean and the resultant radiation dose to exposed human populations.
Present and potential future changes to the global environment have important implications for marine pollution and for the air-sea exchange of both anthropogenic and natural substances. This report addresses three issues related to the potential impact of global change on the air-sea exchange of chemicals: Global change and the air-sea transfer of the nutrients nitrogen and iron. Global change and the air-sea exchange of gases. Oceanic responses to radiative and oxidative changes in the atmosphere. The deposition of atmospheric anthropogenic nitrogen has probably increased biological productivity in coastal regions along many continental margins. Atmospheric deposition of new nitrogen may also have increased productivity somewhat in mid-ocean regions. The projected future increases of nitrogen oxide emissions from Asia, Africa and South America will provide significant increases in the rate of deposition of oxidized nitrogen to the central ...
... 2008. Characterizing the relative contributions of large vessels to total ocean noise fields: a case study using the Gerry E. Studds Stellwagen Bank ...
Cheap air conditioning from oceanic temperature differences is a possible new competitor for coal in coastal regions. Coal suppliers planning expansion in growing economies need to be aware of the possibilities and limitations of this new competitor. Ocean thermal energy conversion (OTEC) used to be uneconomic but the National Energy Laboratory of Hawaii (NELH) has developed a very economic method of utilising the temperature difference between warm ocean surface and cold deep ocean to generate energy for air conditioning. The article outlines the process and mentions demonstrations of OTEC technology to produce electricity and air conditioning at a site in Kailua-Kona, Hawaii. Indonesia and India have viable sites for OTEC installations, and NELH is working with Madras State on a project to construct a commercial OTEC plant in India. 2 photos.
... Mineral and Organic Matter in the Coastal and Near-Coastal Ocean. ... Personal Author(s) : Stavn, Robert H. ; Spiering, Bruce A. ; Gould, Richard W ...
...00 Add to Cart Molten Carbonate Fuel Cells: Modeling , Analysis, Simulation, and Control by Kai Sundmacher (Editor), Achim Kienle (Editor), Hans Josef Pesch (Editor), Joachim F. Berndt (Editor), Gerhard Huppmann (Editor) ...
We analyzed the effect of short-term water deficits at different periods of sunflower (Helianthus annuus L.) leaf development on the spatial and temporal patterns of tissue expansion...Full Text Available
This paper describes a long term course of research into the fundamentals of two phase flow regime dynamics that had lead to a differential description of regime properties. The development has included experiment, microscale simulation, and mathematical averaging to produce a dynamic equation for interfacial area.
We use a novel normal mode analysis of an elastic network model drawn from configurations generated during microsecond all-atom molecular dynamics simulations to analyze the mechanism of auto-inhibition...Full Text Available
The Department of Defense (DoD) has long depended on military support functions enabled by space reconnaissance, surveillance, and target acquisition (RSTA) assets. Future generation satellite capabilities will further push technologies in space - if the right technologies are deployed in the right numbers and with the right payloads. Modeling and simulation play major parts in developing and deploying such assets: 1) system and operational requirements determination, and 2) assessment of military utility of such assets. Each area is critical in a system"s life cycle. Requirements determination cuts across the issues of doctrine, organizations, training, materiel, leader development, personnel, and facilities (DOTMLPF). Military utility assessments are necessary to explore and quantify the military worth/benefit of space RSTA assets to operational commanders. Each of these areas requires relevant modeling/simulation tools ...
A Monte Carlo simulation of the vacuum Bianchi type-IX (mixmaster) cosmology yields a significant correlation between large universe volume and high anisotropy. An analog of the model's chaotic classical behavior is seen in the break up of the universe wave function at large volume into fingers in the corners of the minisuperspace anisotropy potential.
Simulations with a regional chemical transport model show that anthropogenic emissions of volatile organic compounds and nitrogen oxides (NOx = NO + NO2) lead to a dramatic diurnal...Full Text Available
A new type of particle simulationmodel based on the gyrophase-averaged Vlasov and Poisson equations is presented. The reduced system, in which particle gyrations are removed from the equations of motion while the finite Larmor radius effects are still preserved, is most suitable for studying low frequency microinstabilities in magnetized plasmas. It is feasible to simulate an elongated system (L/sub parallel/ >> L/sub perpendicular/) with a three-dimensional grid using the present model without resorting to the usual mode expansion technique, since there is essentially no restriction on the size of ..delta..x/sub parallel/ in a gyrokinetic plasma. The new approach also enables us to further separate the time and spatial scales of the simulation from those associated with global transport through the use of multiple spatial scale expansion. Thus, the model ...
... Southall et al. 2007). However, based on preliminary simulationmodeling that attempted to allow for various uncertainties in ... 1986. Behavioral responses of gray whales to industrial noise:
BackgroundIn mammals, vasopressin (AVP) is released from magnocellular neurons of the hypothalamus when osmotic pressure exceeds a fixed set-point. AVP participates to the hydromineral...Full Text Available
The report describes a computer model developed by TVA to simulate the performance and determine the economics of coal cleaning, or coal cleaning combined with flue gas desulfurization (FGD), for power plant emission control processes over a wide range of user-specified condition...
We have used a sintered glass bead core to simulate the spaces and surfaces of reservoir rock in studies of the bacterial plugging phenomenon that affects waterflood oil recovery operations. The passage...Full Text Available
The Predictive Capability Maturity Model (PCMM) is a new model that can be used to assess the level of maturity of computational modeling and simulation (M&S) efforts. The development of the model is based on both the authors experience and their analysis of similar investigations in the past. The perspective taken in this report is one of judging the usefulness of a predictive capability that relies on the numerical solution to partial differential equations to better inform and improve decision making. The review of past investigations, such as the Software Engineering Institute's Capability Maturity Model Integration and the National Aeronautics and Space Administration and Department of Defense Technology Readiness Levels, indicates that a more restricted, more interpretable method is needed to assess the maturity of an M&S effort. The PCMM ...
Some authors have concluded that spiral structures and shocks do not develop if an adiabatic index gamma > 1.16 is adopted in accretion disc modelling, whilst others have claimed that they obtained well defined spirals and shocks adopting a gamma = 1.2 and a M_2/M_1 = 1 stellar mass ratio. In our opinion, it should be possible to develop spiral structures for low compressibility gas accretion discs if the primary component is a black hole. We considered a primary black hole of 8 solar mass and a small secondary component of 0.5 solar mass to favour spiral structures formations and possible spiral shocks via gas compression due to a strong gravitational attraction. We performed two 3D SPH simulations and two 2D SPH simulations and characterized a low compressibility model and a high compressibility model for each couple of simulations. 2D ...
For a direct injection diesel engine to exhibit its inherent low-emission characteristics, it is essential to improve control accuracy, in particular under transient operating conditions during which emissions are more noted. A model-based control system is established for exhaust gas recirculation (EGR) and air/fuel ratio, using a diesel-engine/vehicle model established on the CAD control system. Described herein are the results of simulation for optimizing control of transient conditions, and the emission-controlling effects confirmed by the field tests. The tests for verification of the established model have confirmed that there are good agreements between the simulated and observed results in terms of air/fuel ratio changing with EGR ratio, exhaust and intake gas and pressure levels, and vehicle speed. The model will be applicable to assessment of emission ...
The goal of this DOE Vision-21 project work scope is to develop an integrated suite of software tools that can be used to simulate and visualize advanced plant concepts. Existing process simulation software does not meet the DOE's objective of ''virtual simulation'' which is needed to evaluate complex cycles. The overall intent of the DOE is to improve predictive tools for cycle analysis, and to improve the component models that are used in turn to simulate the cycle. Advanced component models are available; however, a generic coupling capability that will link the advanced component models to the cycle simulation software remains to be developed. In the current project, the coupling of the cycle analysis and cycle component simulation software will be ...
In this report is described the work effort performed to provide NETL with VE-Suite based Virtual Engineering software and enhanced equipment models to support NETL's Advanced Process Engineering Co-simulation (APECS) framework for advanced power generation systems. Enhancements to the software framework facilitated an important link between APECS and the virtual engineering capabilities provided by VE-Suite (e.g., equipment and process visualization, information assimilation). Model enhancements focused on improving predictions for the performance of entrained flow coal gasifiers and important auxiliary equipment (e.g., Air Separation Units) used in coal gasification systems. In addition, a Reduced Order Model generation tool and software to provide a coupling between APECS/AspenPlus and the GE GateCycle simulation system were developed. CAPE-Open ...
Information about surface soil moisture conditions is of critical importance to real-world applications such as agricultural production, water resource management, flood prediction, fire prediction, water supply, military mobility, etc.. Near-surface soil moisture is currently available from non-ideal sensor configuration observations, and two missions targeted at measuring near-surface soil moisture with ideal sensor configuration are expected before the end of the decade (the European Space Agency (ESA), Soil Moisture and Ocean Salinity (SMOS) mission, and the National Aeronautics and Space Administration (NASA), Hydrospheric states "Hydros" mission). Though remote sensing can make spatially comprehensive measurements of surface soil moisture, it cannot provide information on the entire land surface hydrologic system, and the measurements represent only a snap shot in time. Alternatively, land surface hydrology process models may be used to ...
Traffic jams have become very serious at multiforked road intersections, and conventional pre-timed controls are less effective in such situations. In this article, a new traffic signal control system for multi-forked roads is proposed. First, the cellular automaton (CA) model is used to develop a traffic simulator for multi-forked roads. Next, a stochastic model of a traffic jam is built up. In addition, a new traffic signal control algorithm is designed using the optimization technique and a genetic algorithm (GA). Finally, the effectiveness of the proposed method is shown using actual traffic data with a traffic simulator.
We review the development and application of kinetic Monte Carlo simulations to investigate defect and dopant diffusion in ion implanted silicon. In these type of Monte Carlo models, defects and dopants are treated at the atomic scale, and move according to reaction rates given as input principles. These input parameters can be obtained from first principles calculations and/or empirical molecular dynamics simulations, or can be extracted from fits to experimental data. Time and length scales differing several orders of magnitude can be followed with this method, allowing for direct comparison with experiments. The different approaches are explained and some results presented.
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.
We analyse several amorphous Pdsub(1-x)Sisub(x) alloys obtained by simulation. The metalloid atoms environments can be studied in the same way as that described in the first paper of this series. The main part of this paper is devoted to the evolution of the density as a function of the concentration. In order to interpret the values obtained by the simulation, we develop two models: the first one assumes a complete disorder: the second one, which is shown to be better, assumes a tendency to a local ordering: Si atoms are assumed to have only Pd neighbours.
In FY 1984 three specific tasks which are all related to not-well-understood two-phase phenomena of importance to LWR accidents have been identified under the program. These three tasks are: (1) inverted annular flow experiments and modeling; (2) hot leg U-bend two-phase flow simulation study; and (3) development and evaluation of two-phase flow scaling criteria. Some of the important results obtained under Tasks (1) and (2) are reported in this paper.
This paper presents the results of a model for the impact of potential contamination produced by radionuclide leaking due to long-term physical, chemical and hydrogeological processes occurring in a LILW repository and in the natural environment. The analysis contains a deterministic and also a probabilistic approach for uncertainty assessment. The input function, i.e., the source term was evaluated by using AMBER code and the entire analysis was accomplished by using GOLDSIM, a powerful tool to support probabilistic simulation in management and decision-making in engineering and science. The results obtained were compared to previous simulations and uncertainty analyses (FEHM). (authors)
The underlying theory and model of capacitive voltage substations and its associated problems are described. This system is intended to supply electricity to rural communities directly from high-voltage transmission lines. This technique leads to important cost reductions when compared with conventional alternatives. However, in several transient cases, some problems such as ferroresonance and overvoltages could occur in these substations. Using a developed capacitive substation model these transient cases are simulated. A suitable filter is designed to damp overvoltages clue to ferroresonance. Analytical aspects of filter design are presented and some computer simulation tests are provided to highlight the usefulness of the proposed filter. (author)
The automatic control system for brain tissue temperature is studied theoretically for brain hypothermia treatment. In order to realize a human-friendly control mechanism, an automatic temperature regulation system is constructed to simulate brain hypothermia treatment by introducing a fuzzy algorithm for possible characteristic changes in patients. The brain temperature model is successfully realized to follow the desired temperature course automatically. The model reference fuzzy control of brain temperature based on water-cooling blankets is verified for clinical application to brain hypothermia treatments through various kinds of simulation experiment.
The SMART (System-integrated Modular Advanced ReacTor) which is a 330 MWt advanced integral PWR was developed by the KAERI (Korea Atomic Energy Institute) for electricity generation and seawater desalination. A thermal hydraulic evaluation and analysis of the SMART is performed by the TASS /SMR-S (Transient And Setpoint Simulation/System integrated Modular Reactor-Safety). The TASS/SMR-S code has various models reflecting the design features of the SMART such as the drift flux model, the core models (core power and core heat transfer model), the component models, and the specific models. One of the core models is the core heat transfer model. The role of this model is to calculate the heat flux and radial temperature profiles at a fuel rod surface using the relevant heat transfer ...
Abstract Suspension of solids in stirred reactor is widely used for catalytic reactions, dissolution, etc. Quality of solid suspension is an important parameter required for the reliable design, optimum performance, and scale up of the system. Quality of suspension depends on local characteristics of solid velocity and hold up profiles. The present work was focused on investigating quality of solid suspension using ultrasound velocity profiler (UVP) measurements and CFD simulations. The slip velocity measurements carried out with UVP were used to evaluate different drag correlations used in CFD simulations. Results discussed in this work would be useful for extending the applications of CFD models for simulating large stirred slurry reactors.
In this paper the simulation of hybrid processes containing membrane units will be discussed. For this purpose a user-defined module for simulation and design of membrane processes was implemented into the simulation program AspenPlus. The advantages can be summarized as follows: - any combination of membrane processes with all other units already implemented in AspenPlus is possible, including internal recycle streams, - utilization of the physical property models and data bases of AspenPlus is possible, - cost and sensitivity analysis can be performed. These benefits are demonstrated in detail for a membrane vapor recovery unit for the treatment of tank farm off-gas, for a two-stage reverse osmosis plant for organic/-organic separations and for a combination of distillation and pervaporation for the separation of a dimethylcarbonate/methanol mixture. (orig.)
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 ...
In this paper we present a model and simulator for many clusters of heterogeneous PCs belonging to a local network. These clusters are assumed to be connected to each other through a global network and each cluster is managed via a local scheduler which is shared by many users. We validate our simulator by comparing the experimental and analytical results of a M/M/4 queuing system. These studies indicate that the simulator is consistent. Next, we do the comparison with a real batch system and we obtain an average error of 10.5\\% for the response time and 12\\% for the makespan. We conclude that the simulator is realistic and well describes the behaviour of a large-scale system. Thus we can study the scheduling of our system called \\dirac in a high throughput context. We justify our decentralized, adaptive and opportunistic approach in comparison to a centralized approach in such a ...
We present results of multiple-time-scale simulations of 5, 10 and 15 keV low temperature ion implantation of arsenic on silicon (100), followed by high temperature anneals. The simulations start with a molecular dynamics (MD) calculation of the primary state of damage after 10ps. The results are then coupled to a kinetic Monte Carlo (MC) simulation of bulk defect diffusion and clustering. Dose accumulation is achieved considering that at low temperatures the damage produced in the lattice is stable. After the desired dose is accumulated, the system is annealed at 800 degrees C for several seconds. The results provide information on the evolution for the damage microstructure over macroscopic length and time scales and affords direct comparison to experimental results. We discuss the database of inputs to the MC model and how it affects the diffusion process.
This paper describes the modeling of horizontal steam generator with the TRACE code and calculation results of a loss-of-feedwater (LOF-10) experiment at the PACTEL facility. Parallel Channel Test Loop (PACTEL) is an integral test facility for a VVER-440 type nuclear reactor. The main objectives were to prepare a simulationmodel for its horizontal steam generator with the TRACE thermal hydraulic code and assess different modeling options of the code. PACTEL experiment LOF-10 was chosen for this assessment. The calculation results showed that TRACE is capable in simulating horizontal steam generator behavior both in steady state and during loss-of-feedwater transient. The phenomenon of heat transfer from primary to secondary side, steam superheating and flow reversal in the lowest heat exc...
In a shipping terminal, the pipeline in the pier usually has a diameter smaller than the main pipeline diameter. This diameter reduction amplifies the effects caused by hydraulic transients. It was noticed during simulations that check valves placed in the pier may introduce some error in the maximum pressure when it was modeled in a conventional fashion. The dynamic check valve model had been added to a commercial pipeline simulation software and the results was compared with a similar model using the method of characteristics. Considering a real state shipping terminal case, it was proved even a pipeline that works with low pressure and smooth elevation profile needs a detailed valve modeling to prevent the calculation of unreal pressure values in the region with diameter reduction. (author)
In this work, the pulse plasma nitriding process of pure iron is modelled using a kinetic model derived from the Fick's laws. This approach based on the error function model allows to simulate the growth kinetics of a biphase configuration consisting of {gamma}'-Fe{sub 4}N{sub 1-x} and {epsilon}-Fe{sub 2-3}N iron nitrides under specified conditions. By coupling the kinetics data about the diffusivities of nitrogen in {alpha}-Fe, {gamma}' and {epsilon} phases to the thermodynamic description of the Fe-N phase diagram, it is possible to evaluate the kinetic constant at each growing phase interface, the thickness of nitrided layer, and predict the nitrogen concentration profiles. The comparison between the experimental and theoretical results shows the reliability of the simulation results.
In this work, the pulse plasma nitriding process of pure iron is modelled using a kinetic model derived from the Fick's laws. This approach based on the error function model allows to simulate the growth kinetics of a biphase configuration consisting of #gamma#'-Fe_4N_1_-_x and #epsilon#-Fe_2_-_3N iron nitrides under specified conditions. By coupling the kinetics data about the diffusivities of nitrogen in #alpha#-Fe, #gamma#' and #epsilon# phases to the thermodynamic description of the Fe-N phase diagram, it is possible to evaluate the kinetic constant at each growing phase interface, the thickness of nitrided layer, and predict the nitrogen concentration profiles. The comparison between the experimental and theoretical results shows the reliability of the simulation results.
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...
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...
Various theories and examples of modelling, forecasting and optimization designing in the different parts of the petroleum and gas industries are presented, stochastic programming for long term planning in the refining industry, stochastic model for gasoline blending, feedstock optimization, location and sizing for offshore platforms, hydrocarbon exploration simulation rapid method, valuation of oil field development leases, economic models for petroleum allocation, models for oil supply market, trade embargo game theory, stochastic programming of gas contract portfolio management, scheduling transportation of oil and gas, strategic planning in an oil pipeline company, simulation of offshore oil terminal systems, hierarchical selection of oil and gas distribution systems.
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)
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)
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.
A model for the simulation of the time dependent behavior and the analysis of the equilibrium of the coupled system of storage ring and Free Electron Laser (FEL) is presented. The analysis comprises both amplifier and oscillator FEL. Bunch lengthening and energy widening due to wake forces are taken into account in a self-consistent way. The method is based on a mapping algorithm for means and correlations of the electron distribution function, pioneered by K. Hirata. The evolution of the laser field in the oscillator FEL is described by K. Hirata. The evolution of the laser field in the oscillator FEL is described by supermodes. The model is used to simulate an FEL in a small 500 MeV storage ring with 100 m circumference. Typical values for the output power, spatial, and spectral characteristics of the emitted radiation are presented.
The present work summarizes efforts on the simulation of natural convection cooling within the IFMIF target and test cell. The simulations have been performed with the STAR-CD code using the k-#omega# high-Reynolds number turbulence model. A dedicated thermohydraulic model has been devised including Lithium loop components. Nuclear heat production has been calculated by the Monte-Carlo code McDeLicious for different parts of the target and test cell walls and was used as input for the STAR-CD simulations. Helium atmospheres at several pressures from 0.1 to 10"-"5 MPa have been investigated. In order to limit the maximum temperature of the concrete walls to 80 deg. C it was necessary to add thermal insulation layers to the hot Lithium loop surfaces and a conceptual system of two cooling layers in different depths of the concrete walls.
This study describes the construction of a simple mathematical model its validation through the simulation in transient state of a real cycle performed by a solar refrigerator based on physical adsorption using an activated carbon/methanol pair. The deviation from experimental results was 4% for the cycled mass of methanol, 2.2 % for maximum collector average temperature, and 3 x 10{sup -3} for the theoretical cycle coefficient of performance. Additional simulations of the same cycle inputting values representing different types and larger amounts of activated carbon showed the possibility of increasing the cycled methanol mass up to about 150%. (author) 26 refs., 16 figs., 9 tabs.
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 ...
An important growth in the power of the photovoltaic systems connected to a grid has recently been observed. In spite of the advances in module technology, the problems in the system design increased, especially regarding the surface of the earth they occupy. In this work we propose a complete model for plant simulation with different wiring diagrams and under dynamic shading. Results obtained from simulations showed that the configuration with the lowest performance was that of only one serial-parallel group, whereas the highest efficiency corresponded to a design of groups of modules in parallel connected then in series. In general, a higher efficiency was obtained diminishing the quantity of modules in series and increasing their number in parallel. The simulationmodel proposed allows exploring different alternatives of wiring modules and finding the most efficient ...
We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the IGM, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, evolving from a sub-Mpc scale at the beginning of reionization to R>10 Mpc towards the end. We present a detailed statistical description of ...
The three-dimensional (3-D), multiphase, computational fluid dynamic (CFD) code FLUENT is used to simulated two-phase flow behaviour in a CANDU header manifold under low (natural circulation) flow conditions. This behaviour was previously inferred from experimental data. The CFD simulations reported here are being used to support these inferences and to obtain a better understanding of phase distribution in the header manifold. The simulations seem to show that the vapor-water mixture models in the FLUENT code do not capture properly phase separation in the header and proper phase branching at the header-feeder connections that have been observed in experiments at low flows. The simulations using discrete-phase model in FLUENT, which tracks the pathlines of the individual vapor bubbles in the water continuum phase, show interesting, complicated and, in some ...
In this article, we propose a model selection method, the Bayesian composite model space approach, to map quantitative trait loci (QTL) in a half-sib population for continuous and binary traits. In our method, the identity-by-descent-based variance component model is used. To demonstrate the performance of this model, the method was applied to map QTL underlying production traits on BTA6 in a Chinese half-sib dairy cattle population. A total of four QTLs were detected, whereas only one QTL was identified using the traditional least square (LS) method. We also conducted two simulation experiments to validate the efficiency of our method. The results suggest that the proposed method based on a multiple-QTL model is efficient in mapping multiple QTL for an outbred half-sib population and is more powerful than the LS method based on a single-QTL ...
Lean premixed combustion (LPC) of natural gas is of considerable interest in land-based gas turbines for power generation. However, modeling such combustors and adequately addressing the concerns of LPC, which include emissions of nitrogen oxides, carbon monoxide and unburned hydrocarbons, remains a significant challenge. In this paper, characteristics of published simulations of gas turbine combustion are summarised and methods of modeling turbulent combustion are reviewed. The velocity-composition PDF method is selected for implementation in a new comprehensive model that uses an unstructured-grid flow solver. Reduced mechanisms for methane combustion are evaluated in a partially stirred reactor model. Comprehensive model predictions of swirl-stabilised LPC of natural gas are compared with detailed measurements obtained in a laboratory-scale combustor. The ...
In this paper, a robust water level control system for the horizontal steam generator (SG) using the quantitative feedback theory (QFT) method is presented. To design a robust QFT controller for the nonlinear uncertain SG, control oriented linear models are identified. Then, the nonlinear system is modeled as an uncertain linear time invariant (LTI) system. The robust designed controller is applied to the nonlinear plant model. This nonlinear model is based on a locally linear neuro-fuzzy (LLNF) model. This model is trained using the locally linear model tree (LOLIMOT) algorithm. Finally, simulation results are employed to show the effectiveness of the designed QFT level controller. It is shown that it will ensure the entire designer's water level closed loop specifications.
Accurate modeling of the enhanced diffusion of boron during rapid thermal annealing has been accomplished by incorporating the effects of extended defect formation and annealing on enhanced diffusion into a multizone, semiempirical model. The multizone model divides the implant profile into three zones defining regions of different defects and diffusion enhancements. The model also contains the initial enhanced diffusion and the transient diffusion effects associated with the dissolution of defect clusters and the annealing of extended defects, respectively. The saturation time for transient-enhanced diffusion contains an exponential function of implant dose in order to model the increase in point defect generated with higher implant dose. As a result, the model accurately simulates the boron diffusion profile over a wide range of implant ...
An appropriate mathematical model can help researchers to simulate, evaluate, and control a proton exchange membrane fuel cell (PEMFC) stack system. Because a PEMFC is a nonlinear and strongly coupled system, many assumptions and approximations are considered during modeling. Therefore, some differences are found between model results and the real performance of PEMFCs. To increase the precision of the models so that they can describe better the actual performance, optimization of PEMFC model parameters is essential. In this paper, an artificial bee swarm optimization algorithm, called ABSO, is proposed for optimizing the parameters of a steady-state PEMFC stack model suitable for electrical engineering applications. For studying the usefulness of the proposed algorithm, ABSO-based results...
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 ...
Although experimental results so far agree with predictions of the standard model, it is widely felt to be incomplete. Many prospective theories beyond the standard model predict extra neutral gauge bosons, denoted by Z', which might be light enough to be accessible at the LHC. Observables sensitive to the properties of these extra gauge bosons might be used to discriminate between the different theories beyond the standard model. In the present work several of these observables (total decay width, leptonic cross-section and forward-backward asymmetries) are studied at generation level and with a full simulation in the ATLAS detector. The Z' {yields} e{sup +}e{sup -} decay channel was chosen and 2 values for the mass of Z': 1.5 TeV and 4 TeV. Background is studied as well and it is confirmed that a Z' boson could easily be discovered at the chosen masses. It is shown ...
A mathematical model of 137Cs migration in forest ecosystem is presented, which describes the behaviour of this radionuclide in the forest litter-soil system, trees, understory and forest animals. The model's parameters for different types of forest ecosystems are estimated and model's adequacy is tested through the use of independent experimental data. The sensitivity of the model's output variables is analyzed to variations in the most significant parameters. The differences in the seasonal and mean annual dynamics of 137Cs concentration in muscles of roe deers and mooses are shown to be defined by specific features of the diets of these animals and variations in 137Cs content in the main diet components. PMID:11402557
Current projections of the response of the biosphere to global climatic change indicate as much as 50% to 90% spatial displacement of extratropical biomes. The mechanism of spatial shift could be dominated by either (1) competitive displacement of northern biomes by southern biomes, or (2) drought-induced dieback of areas susceptible to change. The current suite of global biosphere models cannot distinguish between these two processes, thus determining the need for a mechanistically based biome model. The model is in an early stage of development and will require several enhancements, including explicit simulation of potential evapotranspiration, extension to boreal and tropical biomes, a shift from steady-state to transient dynamics, and validation on other continents.
The natural clayey soils are usually structural and unsaturated, which makes their mechanical properties quite different from the remolded saturated soils. A structural constitutive model is proposed to simulate the bonding-breakage micro-mechanism. In this model, the unsaturated soil element is divided into a cementation element and a friction element according to the binary medium theory, and the stress-strain coordination for these two elements is obtained. The cementation element is regarded as elastic, whereas the friction element is regarded as elastoplastic which can be described with the Gallipoli?s model. The theoretical formulation is verified with the comparative experiments of isotropic compressions on the saturated and unsaturated structural soils. Parametric analyses of the e...
Predicted by stochastic models and observed experimentally in a number of isomerization reactions, viscosity-induced solvent effects manifest themselves in a significant departure of the reaction rates from the values expected on the basis of transition state theory. These effects are well understood within the framework of stochastic models; however, the predictive power of such models is limited by the fact that their parameters are not readily available. Experiment and molecular dynamics (MD) simulations can provide such information and can serve as the testing grounds for various stochastic models. In real solvents, a change in viscosity is inevitably associated with variation of at least one of the three factors - temperature, pressure, or solvent identity, resulting in different solv...
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)
This paper introduces a new parsimonious structure for mixture of autoregressive models. the weighting coefficients are determined through latent random variables, following a hidden Markov model. We propose a dynamic programming algorithm for the application of forecasting. We also derive the limiting behavior of unconditional first moment of the process and an appropriate upper bound for the limiting value of the variance. This can be considered as long run behavior of the process. Finally we show convergence and stability of the second moment. Further, we illustrate the efficacy of the proposed model by simulation and forecasting.
CYRANO3 is the new EDF thermomechanical code developed to evaluate the overall fuel rod behavior under irradiation. In that context, this paper presents the phenomena to be simulated and the correlations adopted for modelling purposes. The empirical models presented are taken from the CYRANO2 code and a compilation of the relevant literature. The present revision corrects and supplements version B on the basis of its use during the software coding phase from January 1991 to May 1993. (authors). figs., tabs., 120 refs.
In modern quality engineering, dual response surface methodology is a powerful tool to model an industrial process by using both the mean and the standard deviation of the measurements as the responses. The least squares method in regression is often used to estimate the coefficients in the mean and standard deviation models, and various decision criteria are proposed by researchers to find the optimal conditions. Based on the inherent hierarchical structure of the dual response problems, we propose a Bayesian hierarchical approach to model dual response surfaces. Such an approach is compared with two frequentist least squares methods by using two real data sets and simulated data.
This thesis - Resource Management in Broadband Communication Networks - deals with different ways of optimizing the available resources of data- or telecommunication networks. Especially topics like optimal routing, load balancing and fast recovery of routes in case of link failures are covered. The first part gives a brief description of some of the existing protocols for routing and controlling resources, such as RSVP, OSPF, BGP, PNNI, etc. The remaining part concerns the following two fundamentally different approaches to resource management etc.: - Software agents - Simulated ants In the beginning of the part concerning software agents a description of what exactly constitutes software agents, according to the scientists and according to the organization, FIPA is given. After this, the main results from the IMPACT project are presented. The IMPACT project is an EU-project that aimed at developing a demonstration platform, where software agents handled virtually ...
The Stockpile Stewardship Program (SSP) is a single, highly integrated technical program for maintaining the surety and reliability of the U.S. nuclear stockpile. The SSP uses past nuclear test data along with current and future non-nuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. It includes stockpile surveillance, experimental research, development and engineering (D&E) programs, and an appropriately scaled production capability to support stockpile requirements. This integrated national program requires the continued use of current facilities and programs along with new experimental facilities and computational enhancements to support these programs. The Advanced Simulation and Computing Program (ASC) is a cornerstone of the SSP, providing simulation capabilities and computational resources to support the annual ...
The Stockpile Stewardship Program (SSP) is a single, highly integrated technical program for maintaining the safety and reliability of the U.S. nuclear stockpile. The SSP uses past nuclear test data along with current and future nonnuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. It includes stockpile surveillance, experimental research, development and engineering programs, and an appropriately scaled production capability to support stockpile requirements. This integrated national program requires the continued use of current facilities and programs along with new experimental facilities and computational enhancements to support these programs. The Advanced Simulation and Computing Program (ASC)1 is a cornerstone of the SSP, providing simulation capabilities and computational resources to support the annual stockpile ...
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 ...
We have simulated the transient enhanced diffusion (TED) of boron fluoride (BF{sub 2}{sup +}) implanted in crystalline and germanium amorphized silicon. Based on recently published models, the effect of fluorine on boron diffusion in silicon has been introduced and a modelling has been suggested. In order to simulate the boron experimental profiles, we have assumed that fluorine forms clusters involving interstitial boron which reduces the junction depth. Experimental results indicate that fluorine behaviour depends on amorphization energy. Moreover, even no germanium preamorphization is performed, silicon is still amorphized by fluorine species. Hence, BF{sub 2}{sup +} implantation leads to an amorphous/crystalline (a/c) interface near the surface. An improvement of published models is suggested taking into account fluorine effects. The simulations satisfactory ...
A nuclear power plant engineering (NPP) simulator may include many system functions such as thermal-hydraulic calculation, 3D reactor neutron kinetics model, control and protection system, display and operational human-machine interface, intelligent alarm system, etc. If all those functions are simulated using a single code, the size of the program will cause structural or managerial problems, and the hardware requirement could be tremendous. Thus, any local function errors or changes may affect the whole code, which will make development and maintenance extremely costly. So the nondistributed code is not considered to be flexible and feasible. Currently, distributed and parallel technologies have been applied to develop NPPs engineering simulator for safety analysis, verification of advanced main control room display, operation and intelligent alarm design. The distributed ...
Kelvin-Helmholtz instability of short gravity waves is examined in order to explain the recent findings of the decrease in momentum transfer from hurricane winds to sea waves. A three-fluid configuration of a foam layer between the atmosphere and the ocean is suggested to provide signifficant stabilization of the system and shifting the marginal critical wavelength to the shortwave part of the spectrum. It is conjectured that such stabilization leads to the observed drag reduction. The high contrasts in three fluid densities provide a universal mechanism for stabilizing surface perturbations.
... Rachel has also studied Geological Oceanography at the School of Ocean and Earth Sciences, University of Wales, Bangor, including a year spent at the Institute of Oceanographic Sciences Deacon Laboratory, Surrey, as a student researcher for the British Mid-Ocean Ridge Initiative (BRIDGE). Rachel has authored the following 1 Working Papers: 2002 ecm-2002-05 Towards Integrated Catchment/Coastal Management: Science, Policy ...
The aim of this paper is to provide a general view of wave energy resource assessment. First, a review of the origin of waves and the transformation they undergo as they propagate towards the coast through waters of decreasing depth is presented. Following this, the wave and wave-energy parameters and the statistics required for resource characterization are described. The various types of wave data and their usefulness for the present purposes are summarised. A common methodology for assessment of the wave energy resource is developed. Finally, a general description of the global open ocean resource is presented.
...ENSO Applications Center (PEAC) Institutional Affiliation: PEAC's core members are The National Oceanic and Atmospheric Administration - Office of Global Programs (NOAA/OGP) The NOAA National Weather Service - Pacific Region (NWS-PR), The University of Hawaii - School of Ocean and Earth Science and Technology UH/SOEST, The University of Guam - Water and Energy Research Institute (UOG/WERI), and a regional association of the USAPI Governments, the Pacific Basin ...
The oceanic phytoplancton productivity may essentially influence the total rate of the atmospheric CO{sub 2} absorption by the ocean - that is, a considerable amount of CO{sub 2} will be taken-up in the 50 micrometers thick layer near the air-sea interface. Even if phytoplancton production constitutes only 5% of the total oceanic biota production, this will increase the rate of CO{sub 2} absorption more than twice compared with the present estimates. The reason is that metabolic activity of phytoplancton leads to the emergence in a thin scin (50 micrometers, the average size of phytoplancton cells) layer near the water surface of an additional minimum in the CO{sub 2} partial pressure profile and of an additional maximum of {Delta} {sup 13}C in the same area. These two extremums cannot be detected if the corresponding characteristics are averaged over any microscopic area in the well mixing layer that is more than 1 meter ...
In this work, a comprehensive model including heat transfer, fluid flow and solidification was used to evaluate the performance of a novel delivery system for a single-belt steel casting process. This near-net-shape casting, although still in development, is one of the most promising routes for casting of low-carbon steel in large scale. This paper focuses on the performance of a low-Re {kappa}-{epsilon} model that was employed to simulate the three-dimensional turbulent flows fully-coupled with heat transfer and macroscopic solidification. Simulations were run for the intended conditions of caster operation, and the results obtained with the {kappa}-{epsilon} were compared with the results obtained with an ad hoc viscosity model, where the molecular viscosity was boosted 100 times, uniformly throughout the computational domain. A semi-analytical solution was employed for validating ...
A mathematical model was developed to simulate the combustion process of a direct injection diesel engine. The model considers the key features of such a combustion process; namely fuel droplet evaporation, spray formation, mixing, ignition of heterogeneous burning of the fuel. To model mixing, great emphasis was placed on treating both large and small scale aspects of mixing, the model divides the combustion chamber into three distinct regions of air, fuel spray, and a turbulent plume that is formed once ignition occurs. The fuel spray itself is considered to be composed of multiple zones, which in turn are made of cells of equal mass and of different thermodynamic states. Micromixing of the fuel and air cells, within each zone of the jet and the turbulent plume, is then achieved by utilizing the concept of Monte Carlo stochastic coalescence dispersion technique. The frequency of ...
The RELAP5 code has been developed for best-estimate transient simulation of light water reactor coolant systems during postulated accidents. The code models the coupled behavior of the reactor coolant system and the core for loss-of-coolant accidents and operational transients such as anticipated transient without scram, loss of offsite power, loss of feedwater, and loss of flow. A generic modeling approach is used that permits simulating a variety of thermal hydraulic systems. Control system and secondary system components are included to permit modeling of plant controls, turbines, condensers, and secondary feedwater systems. RELAP5/MOD3 code documentation is divided into seven volumes: Volume I presents modeling theory and associated numerical schemes; Volume II details instructions for code application and input data preparation; Volume III presents the ...
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 simulationmodel. 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 ...
The integrity of the cables located in the target room is very important in the Fire PSA, because the CDF and CCDP are changed according to the results of a cable integrity that depends on the surrounding gas temperature. The conservative assumptions used in the Fire PSA typically specify that all of the equipment and cables of a room would fail when a fire happens in the room. But the realistic assessment of a fire risk by using a fire simulation tool has become necessary in the Fire PSA as described in the ANS Fire PRA Standard. This paper evaluates the cable integrity of eight pump rooms in the nuclear power plant by using the CFAST zone fire model. The upper layer gas temperature of each room is estimated, and an analysis based on the results of modelsimulations is used to judge the cable integrity. According to the analysis results, the integrity of the cable located in the upper layer in the pump ...
The mechanism for afterglow suppression in codoped CsI:Tl,Eu reported in the preceding paper was investigated by combined radioluminescence and thermoluminescence experiments. Model rate equations informed by these experiments were employed to simulate afterglow. It was found that codoping with europium introduces deep electron traps, with room-temperature glow peaks, that effectively scavenge the electrons from shallow traps associated with thallium, thus suppressing afterglow in the time domain of tens of milliseconds.
Involving effects of media, opinion leader and other agents on the opinion of individuals of market society, a trader based model is developed and utilized to simulate price via supply and demand. Pronounced effects are considered with several weights and some personal differences between traders are taken into account. Resulting time series and probabilty distribution function involving a power law for price come out similar to the real ones.
Previous models and evaluation techniques based on the Monte Carlo simulation for the reliability assessment of mixed hydrothermal generation systems are extended. Models for reservoirs on multiple rivers and pumped storage plants are developed. Operating and water management policies are extended to include the relative economical effects of operating alternative sources of generation. These effects and the benefits are described and discussed using an extended IEEE reliability test system. (Author).
As a part of a phenomenological model, a method for simulating the wall/jet interaction in a direct injection diesel engine is proposed. The method is based on the application of the momentum conservation equation in the different directions in which the wall jet is spread, and takes into account both the interaction with the combustion chamber geometry and with swirl. It takes as initial conditions the results of calculating the free jet, which is divided into packages. The predictions provide good agreement with those by other researchers. (author).
In this study, we investigate an inter-industrial and inter-regional recycling system for industrial waste by the cement industry in Japan. We develop a linear programming model that represents cement production processes and waste transportation of all cement factories in Japan. We simulate cost and CO{sub 2}-minimizing systems. The result implies that making waste transportation more efficient in cost is an effective means for CO{sub 2} reduction. (author)
The Battery Life Estimator (BLE) Manual has been prepared to assist developers in their efforts to estimate the calendar life of advanced batteries for automotive applications. Testing requirements and procedures are defined by the various manuals previously published under the United States Advanced Battery Consortium (USABC). The purpose of this manual is to describe and standardize a method for estimating calendar life based on statistical models and degradation data acquired from typical USABC battery testing.
The Advanced Statistical Trajectory Regional Air Pollution (ASTRAP) modelsimulates long-term transport and deposition of oxides of and nitrogen. t is a potential screening tool for assessing long-term effects on regional visibility from sulfur emission sources. owever, a rigorou...
A simple model is described to simulate kinetic processes in dielectric barrier discharges for O_2/NO_x mixtures. A threshold of ozone production found experimentally is confirmed by the calculations of this modeling, and the underlying chemical reaction mechanisms are discussed. It is also found that the effects of diffusion processes in the period of the lifetime of O atoms are not important to micro-discharge channels with a large radius, i.e. larger than 150 #mu#m
On-orbit spectral calibration of hyperspectral imaging data is a key step for quantitatively analyzing them. Like the atmospheric correction, accurate spectral calibration is very necessary for improved studies of land or ocean surface properties. Based on the previous literatures, a new method which coupled an optimization algorithm was developed to simultaneously retrieve the central wavelength and the full width at half maximum (FWHM) of the hyperspectral sensor without needing the in situ reflectance spectra. Firstly, the Hyperion data set simulated using MODTRAN4 with the Hyperion spectral specification was used to test the new method, and the results indicated that the maximum error was less than 0.1 and 0.7 nm for central wavelength and FWHM respectively when the spectral shift is 5 nm. Then the algorithm was applied to the Hyperion data acquired on May 20, 2008 over Heihe River Basin and it was iteratively performed for each detector of ...
The difficulty and unstability of the parameter estimation of mass action model for radioimmunoassay are the major obstacles to its application in routine work. A unique robust curve fitting method for the mass action model was introduced. large numbers of simulated data were generated with Monte Carlo method to compare the accuracy and robustness of the robust fitting and the least squre fitting (Marquardt method) of the mass action model. The results demonstrated that the new robusst method significantly increased the accuracy and reliability of data processing with the mass action model. It also markedly reduced the influence of outliers in the fitting of RIA data, and it was always 'convergent'. This robust method could also be applied to other mathematical models which can be transformed into high degree equations and no initial estimates are needed for ...
The general objective of this workshop is to investigate and discuss methods by which uncertainties in mass balance models for toxics in the Great Lakes may be reduced. As described by the workshop prospectus, this paper is focused on problems of reducing (and quantifying) uncertainty as they relate to in situ field observations/system response measurements for the establishment of initial conditions, boundary conditions, calibration/confirmation data sets, and model post-audit data sets.'' I have taken this description to refer not only to the evaluation of uncertainty in the field observations themselves, but also to the uncertainty associated the analyses of in situ observations as they interact in the overall modeling process. Thus, I will be concerned here with quantification and reduction of uncertainty both (1) as they may be applied to descriptions of the system that is being ...
A mathematical model was developed to predict a moisture content profile during the thick layer re-wetting process of brown rice unpackaged and packaged with low density polyethylene (LDPE) and polybutylene terephthalate (PBT) films. Model validation was carried out by comparing predicted with measured moisture content derived from relative humidity data obtained from the brown rice re-wetting test for 6days at 25degreeC and 90% RH. The moisture standard errors of the model validation for brown rice packaged in LDPE and PBT films were 0.08% wet basis (wb) and 0.11%wb, respectively. It was concluded that the proposed re-wetting model could successfully describe the thick layer re-wetting of brown rice under the experimental conditions. Using this model, re-wetting simulations were carried o...
The effects of pore curvature and surface heterogeneity on the adsorption of water on a graphitic surface at 298 K were investigated using a Grand Canonical Monte Carlo (GCMC) simulation. Slit and cylindrical pores are used to study the curvature effects. To investigate the surface heterogeneity the functional group and the structural defect on the surface were specifically considered. The hydroxyl group (OH) is used as a model for the functional group and the water potential model proposed by Muller et al. is used to calculate the water interaction. For the homogeneous cylinder, the pore filling occurs at a pressure lower than the saturation pressure of the water model, while it is greater in the case of homogeneous slit pore. The size of hysteresis loop is more sensitive to the length of...
In this study, the shear viscosities of wood?plastic composites (WPC)/supercritical CO2 (sc-CO2) solutions at different die temperatures and sc-CO2 contents were measured by a slit die attached on a single-screw extrusion system. A theoretical model based on the Carreau-Yasuda model and Arrhenius equation was then developed to describe the shear-thinning behavior of the WPC/sc-CO2 solution. This model was incorporated into a finite element method software, POLYFLOW, to simulate the pressure and velocity field of WPC/sc-CO2 solution flowing through a foaming die. Based on the simulation results, the pressure drop rate and the location of the supersaturation point were estimated to investigate the cell nucleation. The saturation pressure of sc-CO2 in WPC was determined by a magnetic suspensi...
In order to evaluate postclosure off-site doses from low-level radioactive waste disposal facilities, a computer code was developed to simulate the radionuclide released from waste form, transport through vadose zone and transport in the saturated zone. This paper describes the methodology used to model these process. The radionuclide released from the waste is calculated using a model based on first order kinetics and the transport through porous media was determined using semi-analytical solution of the mass transport equation, considering the limiting case of unidirectional convective transport with three-dimensional dispersion in an isotropic medium. The results obtained in this work were compared with other codes, showing good agreement. (author)
Coumarin (1,2-benzopyrone) is occurring in food, and is also used in cosmetics. In order to perform a risk assessment for both oral and dermal exposure, we applied a physiologically based approach to model kinetics in humans by simulating both routes of exposure. The concentration-time profile in liver revealed a higher peak concentration (Cmax-hep) for the oral when compared to the dermal route. The area under the concentration-time curve in the liver (AUChep) was found the same for both routes if the same extent of absorption is assumed. Dose response information from published rat studies were used to identify the metric relevant for liver toxicity. Liver exposure levels resulting from doses and durations as outlined in the studies were simulated in a rat model. We obtained 31 data pair...
As lidar technology is able to provide fast data collection at a resolution of meters in an atmospheric volume, it is imperative to promote a modeling counterpart of the lidar capability. This paper describes an integrated capability based on data from a scanning water vapor lidar and a high-resolution hydrodynamic model (HIGRAD) equipped with a visualization routine (VIEWER) that simulates the lidar scanning. The purpose is to better understand the spatial and temporal representativeness of the lidar measurements and, in turn, to extend their utility in studying turbulence fields in the atmospheric boundary layer. Raman lidar water vapor data collected over the Pacific warm pool and the simulations with the HIGRAD code are used for identifying the underlying physics and potential aliasing effects of spatially resolved lidar measurements. This capability also helps improve the trade-off between ...
Variability in production is one of the largest factors that negatively impacts construction project performance. A common construction practice to protect production systems from variability is the use of buffers (Bf). Construction practitioners and researchers have proposed buffering approaches for different production situations, but these approaches have faced practical limitations in their application. A multiobjective analytic model (MAM) is proposed to develop a graphical solution for the design of Work-In-Process (WIP) Bf in order to overcome these practical limitations to Bf application, being demonstrated through the scheduling of repetitive building projects. Multiobjective analytic modeling is based on Simulation?Optimization (SO) modeling and Pareto Fronts concepts. Simulation...
This paper presents a mathematical model for simultaneous heat and moisture transfer under freezing-thawing in porous hygroscopic material. Even below O{sup o}C, unfrozen liquid water exists in the material which plays an important role in the processes. The model takes into account the existence of unfrozen liquid water. The system is treated as a three-(gas, liquid and solid) phase system of water. Moisture chemical potential is used as a moisture potential. Under freezing, moisture chemical potential is a unique function of local temperature. So, during freezing, two unknown variables to be solved are solid water content and moisture potential or temperature. A simple example of simulation for freezing processes is shown. (author)
A model for risk and reliability analysis of complex multifunctional production process systems is presented. The model employs Monte-Carlo and Markov Chain algorithms that uses a weighted index to train and simulate the fuzzy hazard data sets which represents failure outcomes of risk component transient and non-transient systems. Early simulation results shows that hazard rates and the risk of containment loss from typical floating production and storage offloading (FPSO)-Riser system for the risk components in parallel or series increases exponentially with time and decreases as safety ratings fraction increases. The reliability value decreases with time and safety fraction (SFAC) for all fuzzy hazard classifications. The results of the computed mean time before repair (MTBR) show that t...
The liquid/vapor coexistence density, the partial vapor pressure, and the heat of vaporization were calculated using Gibbs ensemble Monte Carlo simulation techniques. Long-range interactions such as charge-charge, charge-dipole, and dipole-dipole were evaluated using Ewald summation techniques. A polarizable potential model was used to describe the water-water interactions (Dang and Chang, J. Chem. Phys. 106, 8149, 1997). The model yields good agreement with the corresponding experimental data in the lower temperature region and moderate agreement in the higher temperature region. The critical temperature and density were estimated to be 565 K and 0.28 g/cm3.
We consider the exact path sampling of the squared Bessel process and some other continuous-time Markov processes, such as the CIR model, constant elasticity of variance diffusion model, and hypergeometric diffusions, which can all be obtained from a squared Bessel process by using a change of variable, time and scale transformation, and/or change of measure. All these diffusions are broadly used in mathematical finance for modelling asset prices, market indices, and interest rates. We show how the probability distributions of a squared Bessel bridge and a squared Bessel process with or without absorption at zero are reduced to randomized gamma distributions. Moreover, for absorbing stochastic processes, we develop a new bridge sampling technique based on conditioning on the first hitting time at zero. Such an approach allows us to simplify simulation schemes. New methods are illustrated with pricing ...
The integrity of the cables located in the target room is very important in the Fire PSA, because the CDF and CCDP are changed according to the results of a cable integrity that depends on the surrounding gas temperature. The conservative assumptions used in the Fire PSA typically specify that all of the equipment and cables of a room would fail when a fire happens in the room. But the realistic assessment of a fire risk by using a fire simulation tool has become necessary in the Fire PSA as described in the ANS Fire PRA Standard. This paper evaluates the cable integrity of eight pump rooms in the nuclear power plant by using the CFAST zone fire model. The upper layer gas temperature of each room is estimated, and an analysis based on the results of modelsimulations is used to judge the c...
The availability of hourly air relative humidity (HARH) data is a key requirement for the estimation of epidemic dynamics of plant fungal pathogens, in particular for the simulation of both the germination of the spores and the infection process. Most of the existing epidemic forecasting models require these data as input directly or indirectly, in the latter case for the estimation of leaf wetness duration. In many cases, HARH must be generated because it is not available in historical series and when there is the need to simulate epidemics either on a wide scale or with different climate scenarios. Thirteen modelling solutions (MS) for the generation of this variable were evaluated, with different input requirements and alternative approaches, on a large dataset including several sites a...
Studies on the multivariate autoregressive (MAR) analysis are carried out for the choice of the parameters for modelling the data obtained from various sensors optimally. Accordingly, the roles of the parameters on the analysis results are identified and the related ambiguities are reduced. Experimental investigations are carried out by means of synthesized reactor noise-like data obtained from a digital simulator providing simulated stochastic signals of an operating nuclear reactor so that the simulator constitutes a favourable tool for the present studies aimed. As the system is well defined with its known structure, precise comparison of the MAR analysis results with the true values is performed. With the help of the information gained through the studies carried out, conditions to be taken care of for optimal signal processing in MAR modelling are determined. Although the ...
An overview is provided of an evolving atomistic and kinetic model of semiconductor growth that unifies the main features of strain relaxation in low and high lattice misfit heteroepitaxy. The model reveals a kinetic pathway for dislocation formation during growth with little or no energy cost at low misfits, thus providing a way out of the longstanding dilemma of too high dislocation nucleation energies predicted by classical theories of the equilibrium behavior of a fixed number of particles at low misfits. The essential kinetic process underlying the model are identified on the basis of comparison of the predictions of kinetic Monte-Carlo simulations of growth with real-time or in-situ data obtained in such experiments as reflection high-energy electron diffraction (RHEED) and scanning probe microscopy (SPM). Relative significance of these atomistic kinetic processes is shown to naturally lead to ...
A transient stability program was developed for the Taiwan Power Company, which has a high proportion of nuclear generation in its power system. This program offers a new territory to investigate nuclear plant effects on the power system transient stability. This program also provides a high speed tool for the Taipower system operational planning. A generalized procedure of synchronous machine modeling for a large-scale stability study is presented. The merits and weaknesses of machine modeling can be comprehended through each item of this procedure. Three types of nonlinear synchronous machine modeling implemented into this stability program are derived by following this procedure. A robust subroutine was derived to perform the fourth order Runge-Kutta integration method, making the software programming neat and systematical. For simulating the nuclear plant influence on the system, this program ...
A transient stability program was developed for the Taiwan Power Company, which has a high proportion of nuclear generation in its power system. This program offers a new territory to investigate nuclear plant effects on the power system transient stability. This program also provides a high speed tool for the Taipower system operational planning. A generalized procedure of synchronous machine modeling for a large-scale stability study is presented. The merits and weaknesses of machine modeling can be comprehended through each item of this procedure. Three types of nonlinear synchronous machine modeling implemented into this stability program are derived by following this procedure. A robust subroutine was derived to perform the fourth order Runge-Kutta integration method, making the software programming neat and systematical. For simulating the nuclear plant influence on the system, this program ...
An analysis on landscape level was performed to investigate the growth of Scots pine stands (Pinus sylvestris L.) in the research area of 'Duebener Heide' under the impact of rapidly changing pollution loads and to evaluate their further development. Central to the work was a spatial database, installed using a geographic information system (GIS). This database enabled the statistical analysis of relationships and correlations between the growth of the pine stands, the influence of pollution impacts and the natural site characteristics. The results of emission and immission modelling gave detailed evaluation of the pollution load over the investigated area. The spatial information database was linked with a process-based growth model called FORSANA. The resulting regional model was used to simulate forest growth on stand level for variable time periods. The plausibility of the ...
Pre-stacking depth migration treatment is studied for the estimation of the fold configuration from seismic survey cross sections. The estimation of a velocity structure is necessary for the execution of such treatment, and the utilization of structural-geological knowledge is required for its interpretation. The concept of balanced cross section in relation to the fault-bend fold constructs a stratum structure model under conditions that the deformation during fold and fault formation is a planar strain, that there is no change in volume due to deformation, and that a fold is a parallel fold. In addition to the above geometric and kinetic approach, there is another fold formation process simulationmodel using a Newtonian fluid for study from the viewpoint of dynamics. This simulation stands on the presumption that the boundary contains a ramp that had been in presence before fold formation and that an ...
Concerns over the environment and energy security have recently prompted renewed interest in the U.S. in nuclear energy. Recognizing this, the U.S. Dept. of Energy has launched an initiative to revamp and modernize the role that modeling and simulation plays in the development and operation of nuclear facilities. This Nuclear Energy Advanced Modeling and Simulation (NEAMS) program represents a major investment in the development of new software, with one or more large multi-scale multi-physics capabilities in each of four technical areas associated with the nuclear fuel cycle, as well as additional supporting developments. In conjunction with this, we are designing a software architecture, computational environment, and component framework to integrate the NEAMS technical capabilities and make them more accessible to users. In this report of work very much in progress, we lay out the ...
A small break LOCA event simulation with no high pressure injection system available, known as International Atomic Energy Agency Standard Problem Exercise no. 4 (IAEA-SPE-4), was performed on the PMK-2 integral test facility in Budapest in 1993. This paper analyses the response of the PMK-2 facility, a model of VVER-440 nuclear power plant, using the latest released version MOD3.2.1.2 of the RELAP5 thermal-hydraulic code. After several years of the SPE-4 experiment analyses, many problems have emerged and been studied. Main goal of the present analyses was to study the main influencing parameters for adequate modelling of the hexagonal core channel with 19-rod bundle and phenomena during the core uncovery. Some influencing parameters have been identified, mostly on the primary side, but some also on the secondary side. This is exact simulation of main coolant pump coast down, hydro-accumulators water ...
A set of advanced models implemented into the simulator Sentaurus Process was applied to simulate ultra shallow junction formation by flash lamp annealing (FLA). The full path transient enhanced diffusion model includes equations for small interstitial clusters (I_2, I_3, I_4), #left brace#3 1 1#right brace# defects and dislocation loops. A dopant-point defect clustering model is used for dopant activation simulation. Several cluster types are considered: B_2, B_2I, B_2I_2, B_3I, B_3I_2, B_3I_3 for boron and As_2, As_2V, As_3, As_3V, As_4, As_4V for arsenic. Different point defect and dopant-point defect pair charge states are taken into account to obtain accurate results in the high doping level region. The flux expressions in the three-phase segregation model include a dependence on the doping level and point defect supersaturation. The ...
This guide describes the simulator`s governing equations, constitutive functions and numerical solution algorithms of the STOMP (Subsurface Transport Over Multiple Phases) simulator, a scientific tool for analyzing multiple phase subsurface flow and transport. The STOMP simulator`s fundamental purpose is to produce numerical predictions of thermal and hydrologic flow and transport phenomena in variably saturated subsurface environments, which are contaminated with volatile or nonvolatile organic compounds. Auxiliary applications include numerical predictions of solute transport processes including radioactive chain decay processes. In writing these guides for the STOMP simulator, the authors have assumed that the reader comprehends concepts and theories associated with multiple-phase hydrology, heat transfer, thermodynamics, radioactive chain decay, and nonhysteretic relative permeability, ...
A complete analysis is performed to simulate biomass production incorporated into a realistic whole farm situation, including or replacing a typical crop mix. Representative farms are constructed to accommodate such simulation. Four management systems are simulated for each firm, with each simulation depicting a different crop mix and/or use of different farming technologies and production methods. The first simulation was a base farm plan in which the operator would maintain the historical crop mix for the area, participate in all price support programs, and not participate in either a conservative reserve or a biomass production program. In the second simulation, the operator would again maintain the historical crop mix, would not participate in a conservation reserve or biomass production program, and would be ineligible to participate in any price support ...
A complete analysis is performed to simulate biomass production incorporated into a realistic whole farm situation, including or replacing a typical crop mix. Representative farms are constructed to accommodate such simulation. Four management systems are simulated for each firm, with each simulation depicting a different crop mix and/or use of different farming technologies and production methods. The first simulation was a base farm plan in which the operator would maintain the historical crop mix for the area, participate in all price support programs, and not participate in either a conservative reserve or a biomass production program. In the second simulation, the operator would again maintain the historical crop mix, would not participate in a conservation reserve or biomass production program, and would be ineligible to participate in any price support ...
In this article, we describe an approach to model the electromechanical behavior of the skeletal muscle based on the Huxley formulation. We propose a model that complies with a well established macroscopic behavior of striated muscles where force-length, force-velocity, and Mirsky-Parmley properties are taken into account. These properties are introduced at the microscopic scale and related to a tentative explanation of the phenomena. The method used integrates behavior ranging from the microscopic to the macroscopic scale, and allows the computation of the dynamics of the output force and stiffness controlled by EMG or stimulation parameters. The model can thus be used to simulate and carry out research to develop control strategies using electrical stimulation in the context of rehabilitation. Finally, through animal experiments, we estimated model parameters using a Sigma Point ...
The task of selecting reliable computer programs to simulate the performance of passive solar buildings has been assigned by the CEC to a specialized sub-group, the Modelling Sub-Group. This group has selected a detailed program for large computers, which has been commissioned and run by each participant. It has also given a reference for the analysis of a set of simple models, i.e. for mini/micro computers, pocket calculators and manual methods. These simulation techniques have been analyzed by answering an original questionnaire, making test runs of the models with building and climatic data, and comparing the results obtained with the large program. The best manual method has been included in the Passive Solar Handbook prepared by the Passive Solar Working Group. The MSG has also established a comparative table of all the selected simplified models, in order ...
The UNSAT-H model was developed at Pacific Northwest National Laboratory (PNNL) to assess the water dynamics of arid sites and, in particular, estimate recharge fluxes for scenarios pertinent to waste disposal facilities. During the last 4 years, the UNSAT-H model received support from the Immobilized Waste Program (IWP) of the Hanford Site's River Protection Project. This program is designing and assessing the performance of on-site disposal facilities to receive radioactive wastes that are currently stored in single- and double-shell tanks at the Hanford Site (LMHC 1999). The IWP is interested in estimates of recharge rates for current conditions and long-term scenarios involving the vadose zone disposal of tank wastes. Simulationmodeling with UNSAT-H is one of the methods being used to provide those estimates (e.g., Rockhold et al. 1995; Fayer et al. 1999). To achieve the above goals for ...
The topic of supervised learning within the conceptual framework of artificial neural network (ANN) models is addressed. An ANN is a parallel distributed processing system that consists of many computationally simple processing elements interconnected through uni-directional weighted connections. Such networks, which are roughly patterned after biological nervous systems, have been proposed for use in areas in which the traditional von Neumann computer architecture has been relatively unsuccessful. Learning in these networks is accomplished through the use of algorithms that adjust the values of the connection weights. The work presented here addresses the issue of improving the rate at which ANNs can learn to achieve the mapping of an input pattern to a desired output pattern. The most successful learning algorithms for accomplishing this task are based on gradient descent error minimization techniques. However, the large amount of training time that such ...
Two computational problems were worked on for this study. The first chapter examines the option of coal combustion using oxygen feed with carbon dioxide recycle to control the adiabatic flame temperature. Computer simulations using an existing state-of-the-art 3-dimensional computer code for turbulent reacting flows with reacting particles were employed to study the effects of increased carbon dioxide mole fraction on the char burnout, radiant heat transfer, metal partitioning, and NOx formation. The second chapter compares assumptions for the CO/CO{sub 2} ratio at the surface of mineral inclusions made in previous studies to predictions obtained from a pseudo-steady state kinetic model (SKIPPY) for a single porous particle. The detailed kinetic simulations from SKIPPY for varying particle sizes and bulk gas compositions were used to develop algebraic expressions for the CO/CO{sub 2} ratio that can be incorporated into ...
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 ...
An experiment was conducted on the structure of mesoscopic calculating method for cement-based composite materials. The purpose was utilization for the simulation of mesoscale failure processes and the development of new composite materials. Mesh split based on Voronoi polygons was applied as the discrete limit analysis of the cement-based composite materials containing aggregate particles taking the arrangement of random particles into consideration, and a two-dimensional rigid-body/spring model was employed as the model. The nonhomogeneity of concrete was represented by three kinds of different spring parameters which were aggregate, mortar matrix and their interfaces. The simulation results were summarized as follows: The stress-strain curve went into a non-linear state from the 30% of the maximum stress. The generation of interfacial crack started on the side of the aggregate, but the propagation of ...
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 ...
Several models have been proposed for the simulation of boron diffusion during annealing after implantation in silicon. It has been shown that transient enhanced diffusion (TED) tends to disappear at sub-keV implant energies. Under these conditions boron concentration is higher than the boron solubility limit value, precipitation phenomena occur. Extended defects ({l_brace}3 1 1{r_brace}, EOR) formation and boron precipitation affect both the redistribution during the annealing and activation of the boron. For these reasons, we have investigated the diffusion of low energy boron implanted in crystalline silicon and tested a complete simulation program which takes into account the effects of type I defects as a sink for self-interstitials and of boron precipitation. Experimental results have been simulated and consistent parameters have been found to fit the data. Boron enhanced diffusion effect has been ...
Several models have been proposed for the simulation of boron diffusion during annealing after implantation in silicon. It has been shown that transient enhanced diffusion (TED) tends to disappear at sub-keV implant energies. Under these conditions boron concentration is higher than the boron solubility limit value, precipitation phenomena occur. Extended defects (#left brace#3 1 1#right brace#, EOR) formation and boron precipitation affect both the redistribution during the annealing and activation of the boron. For these reasons, we have investigated the diffusion of low energy boron implanted in crystalline silicon and tested a complete simulation program which takes into account the effects of type I defects as a sink for self-interstitials and of boron precipitation. Experimental results have been simulated and consistent parameters have been found to fit the data. Boron enhanced diffusion effect ...
This report describes accomplishment of simulations of Pyrochemical Process Operation by using virtual engineering models. The pyrochemical process using molten salt electrorefining would introduce new technologies for new fuels of particle oxide, particle nitride and metallic fuels. This system is a batch treatment system of reprocessing and re-fabrication, which transports products of solid form from a process to next process. As a results, this system needs automated transport system for process operations by robotics. In this study, a simulation code system has been prepared, which provides virtual engineering environment to evaluate the pyrochemical process operation of a batch treatment system using handling robots. And the simulation study has been conducted to evaluate the required system functions, which are the function of handling robots, the interactions between robot and process equipment, ...
The present paper is concerned with development and application of a so-called Effective Convection Model (ECM), which aims to provide a detailed, mechanistic description of heat transfer processes in a BWR lower plenum. The ECM is a Computational Fluid Dynamics (CFD)-like tool which employs a simpler and more effective approach to compute heat transfer by solving only energy conservation equation instead of solving the full set of Navier-Stokes and energy equations by a CFD code. We implement the ECM in a CFD code (Fluent), with detailed description of the ECM development, implementation and validation. A dual approach is used to validate the ECM, namely validation against experimental data and against heat transfer results obtained by CFD predictions in the same geometries and conditions. Insights gained from CFD simulations are also used to improve ECM. The ECM capability as an effective tool to simulate heat transfer of ...
The present paper is concerned with development and application of a so-called Effective Convection Model (ECM), which aims to provide a detailed, mechanistic description of heat transfer processes in a BWR lower plenum. The ECM is a Computational Fluid Dynamics (CFD)-like tool which employs a simpler and more effective approach to compute heat transfer by solving only energy conservation equation instead of solving the full set of Navier-Stokes and energy equations by a CFD code. We implement the ECM in a CFD code (Fluent), with detailed description of the ECM development, implementation and validation. A dual approach is used to validate the ECM, namely validation against experimental data and against heat transfer results obtained by CFD predictions in the same geometries and conditions. Insights gained from CFD simulations are also used to improve ECM. The ECM capability as an effective tool to simulate heat transfer of ...
In this article two integral computational fluid dynamics methods for steady-state and transient vehicle aerodynamic simulations are described using a Chevrolet Corvette ZR-1 surface panel model. In the last decade, road-vehicle aerodynamics have become an important design consideration. Originally, the design of low-drag shapes was given high priority due to worldwide fuel shortages that occurred in the mid-seventies. More recently, there has been increased interest in the role aerodynamics play in vehicle stability and passenger safety. Consequently, transient aerodynamics and the aerodynamics of vehicle in yaw have become important issues at the design stage. While there has been tremendous progress in Navier-Stokes methodology in the last few years, the physics of bluff-body aerodynamics are still very difficult to model correctly. Moreover, the computational effort to perform Navier-Stokes ...
This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users Guide. The Xyce Parallel Electronic Simulator has been written to support, in a rigorous manner, the simulation needs of the Sandia National Laboratories electrical designers. It is targeted specifically to run on large-scale parallel computing platforms but also runs well on a variety of architectures including single processor workstations. It also aims to support a variety of devices and models specific to Sandia needs. This document is intended to complement the Xyce Users ...
The Stockpile Stewardship Program (SSP) is a single, highly integrated technical program for maintaining the safety and reliability of the U.S. nuclear stockpile. The SSP uses past nuclear test data along with future non-nuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. It includes stockpile surveillance, experimental research, development and engineering programs, and an appropriately scaled production capability to support stockpile requirements. This integrated national program will require the continued use of current facilities and programs along with new experimental facilities and computational enhancements to support these programs. The Advanced Simulation and Computing program (ASC) is a cornerstone of the SSP, providing simulation capabilities and computational resources to support the annual stockpile assessment ...
We study luminosities of millisecond pulsars in globular clusters by fitting the observed luminosity distribution with single and double power laws. We use simulations to model the observed distribution as the brighter part of some parent distribution for Terzan 5 and try to find a model which simultaneously agrees with the observed diffuse radio flux, total predicted number of pulsars and observed luminosity distribution. We find that wide ranges of parameters for log-normal and power-law distributions give such good models. No clear difference between the luminosity distributions of millisecond pulsars in globular clusters and normal disk pulsars was seen.
A non-linear mathematical model of dynamics of horizontal steam generator for nuclear power unit with WWER type reactor is presented. To realize this model the GEMMA-120 simulation language for computer Odra-1204 has been used. Necessity of taking into account disposited thermal storage capacities along tubulation of a primary cycle is demonstrated. A number of lumped elements of reactor division against a required static accuracy of calculations has been determined. (author).
Thermal striping is identified as one of the causes of thermal fatigue failure in nuclear power plants. Numerical studies of thermal striping require three-dimensional, unsteady turbulent modeling that resolves both large and small-scale turbulent motions. Benchmark studies were carried out using the LES turbulence model solved by the commercial CFD code FLUENT. Two types of mixing tee configurations were modeled to evaluate the performance of the CFD code. The simulation results presented in normalized average temperature and normalized fluctuating temperatures are in good agreement with measurements.
The production of many important commercial materials involves the evolution of a complex fluid through a cooling phase into a hardened product. Textile fibers, high-strength fibers(KEVLAR, VECTRAN), plastics, chopped-fiber compounds, and fiber optical cable are such materials. Industry desires to replace experiments with on-line, real time models of these processes. Solutions to the problems are not just a matter of technology transfer, but require a fundamental description and simulation of the processes. Goals of the project are to develop models that can be used to optimize macroscopic properties of the solid product, to identify sources of undesirable defects, and to seek boundary-temperature and flow-and-material controls to optimize desired properties.
The orbit response matrix (ORM) method [1] is applied to model the Fermilab Booster with parameters such as the BPM gains and rolls, and parameters in the lattice model, including the gradient errors and magnets rolls. We found that the gradients and rolls of the adjacent combined-function magnets were deeply correlated, preventing full determination of the model parameters. Suitable constraints of the parameters were introduced to guarantee an unique, equivalent solution. Simulations show that such solution preserves proper combinations of the adjacent parameters. The result shows that the gradient errors of combined-function magnets are within design limits.
In the study a preliminary mid-term analysis of the BARCOM test model is presented. The BARCOM test model is a 1:4 scale of an existing pressurized heavy water reactor (PHWR) pre-stressed concrete inner containment of 540 MW Tarapur Atomic Power Station 3 and 4 units in India. The goal of this midterm analysis is to illustrate the modeling approach and achieve a prediction of the failure mode. The analysis was carried out using ABAQUS/CAE and ABAQUS/EXPLICIT version 6.7-EF1 software
The RELAP5/MOD2 and MELCOR 1.8.1 codes have been used for simulate the PMK total loss of feedwater with secondary bleed and feed experiments done in a scale-model WWER-440 test facility. Nodalization studies and studies on several-core modelling options were also done. Good agreement was found between the calculations done by RELAP5/MOD2 and MELCOR 1.8.1 JY codes. (orig.) (5 refs., 31 figs.).
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.
In the present work we introduce a novel multi-agent model with the aim to reproduce the dynamics of a double auction market at microscopic time scale through a faithful simulation of the matching mechanics in the limit order book. The agents follow a noise decision making process where their actions are related to a stochastic variable, the market sentiment, which we define as a mixture of public and private information. The model, despite making just few basic assumptions over the trading strategies of the agents, is able to reproduce several empirical features of the high-frequency dynamics of the market microstructure not only related to the price movements but also to the deposition of the orders in the book.
In this study we tried to develop the long-range transport system and find the way to prevent from the radiological emergency risk. For the study, meteorological forecast system in Korea Meteorological Administration is investigated. Numerical simulation is also carried out by the long-range transport model and Vis-5D. We surveyed the emergency preparedness for nuclear accidents which were ARAC in USA, RODOS in Europe and WSPEED in Japan and then investigated the processing of medium- and long-range atmospheric diffusion modeling system. We also studied on the application of KMA/NWPD model which are GDAPS and RDAPS. In the future, it is necessary to produce to the high resolution meteorological data from KMA/NWPD for the development of medium- and long-range atmospheric diffusion modeling system and construct the integrated system for data processing in real time. It was ...
The summertime tropical cyclogenesis rate per unit area in the eastern Pacific ocean is arguably higher than in any other location in the world. Many if not most of these cyclones form from African easterly waves which cross Central America into the Pacific. Of order 25% of these waves intensify into cyclones. A significant fraction of east Pacific tropical cyclones undergoes landfall on the Mexican coast. Those which do not, generally dissipate over cold ocean waters north of the east Pacific intertropical convergence zone, often not far from land. The layer of warm ocean water which supports the development of east Pacific cyclones is unusually shallow and is structured by anticyclonic vortices which form by various processes and propagate slowly to the west. These vortices locally deepen the oceanic mixed layer and support stronger convection than their surroundings, possibly promoting ...
The Formerly Utilized Sites Remedial Action Program (FUSRAP) of the Department of Energy is designed to identify and evaluate the radiological conditions at sites formerly used by the Corps of Engineers Manhattan Engineer District and the US Atomic Energy Commission. Where required, remedial action will be instituted to remove potential restrictions on the use of the sites due to residual low-level radioactive contamination. A total of 31 sites that may require remedial action has been identified. The purpose of the Ocean FUSRAP Program, which began in March 1981, is to assess the technical, environmental, and institutional feasibility of disposing, in the ocean and on the ocean floor, of FUSRAP soil and rubble which contains traces of natural radioactive materials. The initial focus has been on the Middlesex, New Jersey, Sampling Plant site and surrounding properties, which contain on the order of 100,000 metric tons of ...
The Formerly Utilized Sites Remedial Action Program (FUSRAP) of the Department of Energy is designed to identify and evaluate the radiological conditions at sites formerly used by the Corps of Engineers Manhattan Engineer District and the US Atomic Energy Commission. Where required, remedial action will be instituted to remove potential restrictions on the use of the sites due to residual low-level radioactive contamination. A total of 31 sites that may require remedial action has been identified. The purpose of the Ocean FUSRAP Program, which began in March 1981, is to assess the technical, environmental, and institutional feasibility of disposing, in the ocean and on the ocean floor, of FUSRAP soil and rubble which contains traces of natural radioactive materials. The initial focus has been on the Middlesex, New Jersey, Sampling Plant site and surrounding properties, which contain on the order of 100,000 metric tons of ...
This paper reports on the first U.S. tanker converted to a floating production, storage, and offloading (FPSO) unit which takes up station in Gombe-Beta field off Gabon by Dec. 1. FPSO Ocean Producer will work under a 3 year, day rate contract let late in 1990 by Amoco-Gabon Bombe Marin co., a unit of Amoco Production Co. (OGJ, Dec. 24, 1990, p. 27). Gombe-Beta field is in the Atlantic Ocean about 70 miles south of Port Gentil, Gabon. Ocean Producer will be moored in 50 ft of water 3.7 miles off Gabon, with Bombe-Beta's unmanned production platform about 820 ft astern. The vessel will be held in position by a disconnectable, asymmetric, six point, spread mooring system, It is owned and operated by Oceaneering International Services Ltd. (OISL). Affiliate Oceaneering Production Systems (OPS) converted the 78,061 dwt oil tanker MT Baltimore Sea at a capital cost of $25 ...
The vortex-excited dynamics of a uniform pivoted cylinder in uniform and sheared flow was investigated experimentally. The experiments were numerically simulated using a diffusive Van der Pol oscillator model developed by Balasubramanian & Skop recently. Salient features of the experimental investigations and the numerical simulations are presented here. Comparisons between the experimentally recorded and numerically predicted structural response to vortex-excited vibrations, power spectral density measurements of near-wake velocity fluctuations and lock-in ranges are made. A comparison of the numerical predictions and the experimental data reveals good agreement.
Brachytherapy is a radiotherapeutic modality that makes use of radionuclides to deliver a high radiation dose to a well-defined volume while sparing surrounding healthy structures. At the National Institute for Cancer Research of Genova a High Dose Rate remote afterloading system provides Ir(192) endocavitary brachytherapy treatments. We studied the possibility to use the Geant4 Monte Carlo simulation toolkit in brachytherapy for calculation of complex physical parameters, not directly available by experiment al measurements, used in treatment planning dose deposition models.
A transient eddy current measurement method is presented to determine the thickness and conductivity of a conductive plate. The conductive plate is induced by an air-cored coil, the magnetic flux density along the axial is measured and the various signals corresponding to plates with different thickness and conductivity are calculated using a 3D transient eddy current simulator. Characteristic features are obtained from the transient response. A similarity-based modeling method is utilized in this study to estimate the thickness and conductivity of the conductive plate. (author)
Convective losses arising from buoyancy driven flow were calculated for a two-dimensional modelsimulating a solar cavity receiver. The TEMPEST code, capable of fully three-dimensional coupled thermal-hydraulic transient calculations, was used for the simulation. Predicted velocity and temperature results for a 2.59 m deep by 2.88 m high rectangular cavity with an aperture opening of 1.72 m were used to determine convective losses for prescribed interior wall temperatures and cavity orientation. Velocity vector and temperature isotherm plots were used to analyze flow characteristics.
Due to the low thermal conductivity of ceramics large temperature gradients are present through the powder particles during plasma spray deposition. As a result the particles often impinge at the substrate in a semi-molten form; which in turn substantially affects the final characteristics of the coating being formed. This study is dedicated to the novel modelling development and simulation of a semi-molten droplet impingement. The study examines the impingement process during impact, spreading and solidification of semi-molten zirconia. The simulation provides an insight to the heat transfer process during impact and solidification of a semi-molten powder particle and illustrates the freezing-induced break-up mechanism at the splat periphery.
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 ...
Within the APROS simulation environment we have built a thermo-hydraulic model of a conceptual fusion power plant which is water cooled and uses lithium-lead for tritium breeding. For the safety assessment of this design we have studied an accident sequence which starts from a loss or coolant flow then leads to first wall breach and pressurisation of the vacuum vessel. Simulations have revealed strong pressure transients which can be alleviated by design changes. One goal is to verify the adequacy of the containment design: it remains intact at least 14 h without any mitigating efforts. Estimates for radioactive releases are obtained. (author)
The performance of a solar chemical heat pipe was studied using CO{sub 2}reforming of methane as the endothermic reaction. A directly heated vertical reactor, packed with a rhodium catalyst was used. The solar tests were carried out in the Schaeffer solar furnace of the Weizmann Institute of Science. The power absorbed was up to 6.3 KW, the maximal flow rates of the gases reached 11,000 1/h, and the methane conversions reached 85%. A computer model was developed to simulate the process. Agreement of the calculations with the experimental results was quite satisfactory.
In 2009, the Atmospheric Radiation Measurement (ARM) Program and the Climate Change Prediction Program (CCPP) have been asked to produce joint science metrics. For CCPP, the second quarter metrics are reported in Evaluation of Simulated Precipitation in CCSM3: Annual Cycle Performance Metrics at Watershed Scales. For ARM, the metrics will produce and make available new continuous time series of radiative fluxes based on one year of observations from Barrow, Alaska, during the International Polar Year and report on comparisons of observations with baseline simulations of the Community Climate System Model (CCSM).
The overall objective of this project is to provide research guidance and quantification of research progress in the areas of direct and indirect coal liquefaction, coal/waste coprocessing, refining of coal-derived liquid fuels, and natural gas conversion. Specifically, the work is divided into two subtasks that relate to whether the technology application is direct or indirect. In subtask (a), Direct Coal Liquefaction technology is the subject of the analyses, and in subtask (b), Indirect Liquefaction, technologies will be evaluated in accordance with the priorities of the COR. Mitretek Systems has been developing detailed computer simulationmodels of direct and indirect coal and natural gas conversion systems for several years. These models are constantly being updated and improved as more data and better cost information becomes available. These models also include detailed refinery ...
Presented here is the calculation results of tidal field of Daya Bay by using pseudospectral method. According to the calculation results about the tidal field and measured value of wind in Daya Bay area, the diffusion of radionuclide from Daya Bay Nuclear Power Plant is simulated by using particle random walk model. The simulation results about the relative concentration distribution of the radionuclides and the trajectory of the distribution centrode are obtained, taking account of tidal current, wind and turbulent current. These results can provide information not only for the reasonable discharge of the radionuclides from Daya Bay Nuclear Power Plant, but also for understanding the potential effect on environment due to the normal operation of Daya Bay Nuclear Power Plant. (authors)
substrate). Within each group of simulations, three lubricant film thicknesses are studied to examine the effect of varying lubricant thickness. Statistical data are collected from each simulation and presented in this work. Via these data, together with the evolution, of atomic and molecular configurations, a very detailed picture of the properties of this thin film interface is presented. In particular, we conclude that perfluoropolyether lubricant forms distinct molecular layers when confined between two substrates, the rate of heat generation under shearing conditions typical of those in a head-disk interface is insufficient for thermal mechanisms to result directly in lubricant degradation, and mechanical stresses attained in the head-disk interface are unlikely to result in any significant degree of lubricant degradation. This thesis examines the tribology of a head-disk interface in an operating hard disk drive via non-equilibrium ...
The concepts of general flow and thermal boundary conditions are introduced to treat the interaction between indoor and outdoor thermal environments in computational fluid dynamics (CFD) simulations. The results from a multi-zone model are used to supply the general flow boundary conditions. The energy balance equation at wall-air interfaces is used to supply the general thermal boundary conditions. An example calculation in a 13-room building shows that infiltration influences indoor air flow patterns considerably. The air flow in a room ventilated by displacement is measured and simulated. Two surface coatings are considered, i.e. black walls and aluminum walls. Implementation of these two boundary conditions is essential in predicting air flow patterns, air quality, and thermal comfort in a real building. (author)
An innovative hybrid powder prepared using computer simulation allowed development of a new foundation having a fine, smooth texture that has never been achieved before. The optical structure/design of the powder was based on the results of measurements and analyses conducted on the optical characteristics of a baby's delicate skin, which is the envy of many women. To obtain the optimal optical characteristics, the finite differential time domain (FDTD) method for solving Maxwell's differential equation by difference and time domain was applied to the computer simulation method. For synthesis of the hybrid powder based on the optical model, a proprietary shape regulation coating technology was used in which flaky substrates were coated with microspherical forms of barium sulfate crystals. ...
The work presented in this paper is part of the cooperative research project AUTO-OPT carried out by twelve partners from the automotive industries. One major work package concerns the application of data mining methods in the area of automotive design. Suitable methods for data preparation and data analysis are developed. The objective of the work is the re-use of data stored in the crash-simulation department at BMW in order to gain deeper insight into the interrelations between the geometric variations of the car during its design and its performance in crash testing. In this paper a method for data analysis of finite element models and results from crash simulation is proposed and application to recent data from the industrial partner BMW is demonstrated. All necessary steps from data pre-processing to re-integration into the working environment of the engineer are covered.
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 ...
An integrated beam optics-nuclear processes framework is essential for accurate simulation of fragment separator beam dynamics. The code COSY INFINITY provides powerful differential algebraic methods for modeling and beam dynamics simulations in absence of beam-material interactions. However, these interactions are key for accurately simulating the dynamics of heavy ion fragmentation and fission. We have developed an extended version of the code that includes these interactions, and a set of new tools that allow efficient and accurate particle transport: by transfer map in vacuum and by Monte Carlo methods in materials. The new framework is presented, along with several examples from a preliminary layout of a fragment separator for a facility for rare isotope beams.
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.
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 ...
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 ...
The simulation of the Auxiliary Charcoal Bed (ACB) Vacuum System was performed to evaluate the original vacuum system design, detect and identify design deficiencies, investigate the effects of proposed corrections on system performance, and generally aid in refining the system design before construction and mockup testing. The simulation was performed by using the Advanced Continuous Simulation Language (ACSL). The vacuum system design goals are to provide approximately 20 SCFM of both booster gas and purge gas through the system and maintain a flow of approximately 40 SCFM with a velocity of 50 to 75 f/sec at the entrance to the cyclone separator. The model results showed that the original system design was incapable of meeting the system performance goals. Further simulations showed that the following modifications to the original vacuum system design were required to make the ...
This study deals with testing the capability of the code PARET to simulate natural convection cooling phenomena under different boundary conditions. In addition to applying and testing some new options related to simulation of the control rod movement and studying the reactivity effect of thermal expansion fuel elements. The experiments of the simple thermal hydraulic loop of Missouri university about natural cooling phenomena in two narrow paralled channels were used to validate the code. The study indicate good results regarding the distribution of coolant flux velocity and clad temperature. In particular the heat transfer coefficient of natural convection has been calculated in good agreement with the experiment. On the other hand, the core of MNSR reactor has been modelled to simulate the reactor dynamic behaviour under natural convection cooling conditions for different initial power level. The ...
Climates at regional scale are strongly dependent on the interaction between atmosphere and its lower boundary, the oceans and the land surface mosaic. Land surfaces influence climate through their albedo, and the aerodynamic roughness, the processes of the biosphere and many soil hydrological properties; all these factors vary considerably geographically. Land surfaces receive a certain portion of the solar irradiance depending on the cloudiness, atmospheric transparency and surface albedo. Short-wave solar irradiance is the source of the heat energy exchange at the earth`s surface and also regulates many biological processes, e.g. photosynthesis. Methods for estimating solar irradiance, atmospheric transparency and surface albedo were reviewed during the course of this project. The solar energy at earth`s surface is consumed for heating the soil and the lower atmosphere. Where moisture is available, evaporation is one of the key components of the surface energy ...
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 ...
Since its invention by Los Alamos during the Manhattan Project era the exploding bridgewire detonator (EBW) has seen tremendous use and study. Recent development of a laser-powered device with detonation properties similar to an EBW is reviving interest in the basic physics of the deflagration-to-detonation (DDT) process in both of these devices. Cutback experiments using both laser interferometry and streak camera observations are providing new insight into the initiation mechanism in EBWs. These measurements are being correlated to a DDT model of compaction to detonation and shock to detonation developed previously by Xu and Stewart. The DDT model is incorporated into a high-resolution, multi-material model code for simulating the complete process. Model formulation and the modeling issues required to describe the test data will be discussed.
Because the state of a free-floating space robot model is uncertain and sudden changes in the model parameters might undermine the stability of the system, this paper proposes a control strategy based on a variable structure neural integrated controller. This scheme does not need a precise space robot model, making use of the radial basis function neural network ability approach to learn about an uncertain model. The network weights are adjusted online in real-time. During the early period of the control phase and parameter changes, the variable structure controller compensates for the uncertain model which the neural network could not learn well. It also creates global asymptotic stability for the whole closed-loop system. Simulation results show that the controller can handle bad changea...
In general, macroscopic material failure is a manifestation of irreversible changes at the microscopic level. Many tissues, which may appear to be macroscopically homogeneous, are, at a fundamental microscopic level, a composite material. For example, cornea is composed of a hyaluronic acid matrix in which layers of collagen fibers are overlaid in a crossing pattern. The points where the collagen fibers intersect are potential nucleation sites for microscopic defects, which under the action of tensile stress, nucleate, grow and coalesce to form macroscopic failure planes, or spall planes. Using a model based on microstructural evolution, this paper examines the failure process during photoablation. Specifically, the paper describes a physically motivated, micromechanical model based on the nucleation and growth of spherical voids. This model is then used to simulate photoablation of cornea. Potential ...
As an integral part of the Fast Test Reactor Vibration Program for Reactor Internals, the flow-induced vibrational characteristics of scaled Fast Test Reactor core internal and peripheral components were assessed under scaled and simulated prototype flow conditions in the Hydraulic Core Mockup. The Hydraulic Core Mockup, a 0.285 geometric scale model, was designed to model the vibrational and hydraulic characteristics of the Fast Test Reactor. Model component vibrational characteristics were measured and determined over a range of 36 percent to 111 percent of the scaled prototype design flow. Selected model and prototype components were shaker tested to establish modal characteristics. The dynamic response of the Hydraulic Core Mockup components exhibited no anomalous flow-rate dependent or modal characteristics, and prototype response predictions were adjudged acceptable.
A model is described for simulating thermal hydraulic and chemical conditions within fuel crud deposits. Heat transfer takes place by wick boiling in which water flows through the porous deposit and evaporates into steam at the surface of chimneys. The transport and chemistry of dissolved species within the deposit is also modelled. This chemistry includes the equilibrium chemistry of Li/boric acid species, the equilibrium chemistry of Fe/Ni species and the radiolysis chemistry of water. The unique feature of this model is that the chemistry is coupled to the thermal hydraulics via the increase in the saturation temperature with the concentration of dissolved species. This has a profound effect on evaporative heat transfer within thick deposits, leading to conditions that explain the precipitation of LiBO{sub 2} and the possible formation of bonaccordite. The model helps understand ...
Electrochemical discharge machining (ECDM) is a promising hybrid process for high-performance machining of non-conductive glass. ECDM drilling has been found to have different characteristics and material removal mechanisms in discharge regime (less than 300??m in depth) and hydrodynamic regime (more than 300??m in depth); however, these regimes are never separately modeled in existing ECDM models, which leads to large prediction error, especially at low applied voltages and high machining depths. Until now, no model is particularly designed for discharge regime, where most material is removed. In this paper, a finite element based model for ECDM drilling in discharge regime is presented. Material removal subjected to a single spark was simulated using finite element method. The drilling d...
Spratly Islands, located in the southern part of the South China Sea (SCS), consist of more than 100 small islands, coral reefs and banks. Remote sensing is the only way to obtain a synoptic view of all of the islands in such a large area. It has been demonstrated that satellite synthetic aperture radar (SAR) imagery is a very powerful tool for monitoring meso-scale and small-scale ocean processes in a large area. In this study, satellite SAR images were used to study the ocean environment in the area of Spratly Islands. The aim was to understand the capability of satellite remote sensing to monitor ocean processes and provide information for future field studies. Two sets of high-resolution European Remote Sensing satellite (ERS)-2 SAR images over the entire Spratly Islands area were coll...
Various SST indices in the Indo-Pacific region have been proposed in the literature in light of a long-range seasonal forecasting of the Indian Summer Monsoon (ISM). However, the dynamics associated with these different indices have never been compared in detail. To this end, the present work re-examines the variabilities of ISM rainfall, onset and withdrawal dates at interannual timescales and explores their relationships with El Ni?o-Southern Oscillation (ENSO) and various modes of coupled variability in the Indian Ocean. Based on recent findings in the literature, five SST indices are considered here: Ni?o3.4 SST index in December?January both preceding [Nino(?1)] and following the ISM [Nino(0)], South East Indian Ocean (SEIO) SST in February?March, the Indian Ocean Basin (IOB) mode in ...
Man's wastes are polluting the oceans from various sources and delivered by various routes. The result has been the loss of habitat and the irreversible altering of marine ecosystems. Development in the coastal zones and offshore activities that disrupt biologically sensitive areas have led to international negotiations to regulate these impacts and develop a law of the sea. Basic to international cooperation, however, is the need to develop answers to questions about the ecological consequences of development programs and the carrying capacity of the oceans. Current information does not demonstrate long-term global impacts, making it difficult to detect and predict incremental changes and causal relationships. Society needs to determine how much environmental damage it will accept and then to develop appropriate technology, such as biomonitoring. Society also needs to reexamine its positions on technology-based regulatory controls to ...
The total amount of uranium dissolved in seawater at a uniform concentration of 3 mg U/m{sup 3} in the world's oceans is 4.5 billion tons. An adsorption method using polymeric adsorbents capable of specifically recovering uranium from seawater is reported to be economically feasible. A uranium-specific nonwoven fabric was used as the adsorbent packed in an adsorption cage 16 m{sup 2} in cross-sectional area and 16 cm in height. We submerged three adsorption cages in the Pacific Ocean at a depth of 20 m at 7 km offshore of Japan. The three adsorption cages consisted of stacks of 52 000 sheets of the uranium-specific non-woven fabric with a total mass of 350 kg. The total amount of uranium recovered by the nonwoven fabric was >1 kg in terms of yellow cake during a total submersion time of 240 days in the ocean.
The total amount of uranium dissolved in seawater at a uniform concentration of 3 mg U/m3 in the world's oceans is 4.5 billion tons. An adsorption method using polymeric adsorbents capable of specifically recovering uranium from seawater is reported to be economically feasible. A uranium-specific nonwoven fabric was used as the adsorbent packed in an adsorption cage 16 m2 in cross-sectional area and 16 cm in height. We submerged three adsorption cages in the Pacific Ocean at a depth of 20 m at 7 km offshore of Japan. The three adsorption cages consisted of stacks of 52 000 sheets of the uranium-specific non-woven fabric with a total mass of 350 kg. The total amount of uranium recovered by the nonwoven fabric was >1 kg in terms of yellow cake during a total submersion time of 240 days in the ocean.
This paper presents a new technique--Integrated Bayesian Uncertainty Estimator (IBUNE) to account for the major uncertainties of hydrologic rainfall-runoff predictions explicitly. The uncertainties from the input (forcing) data--mainly the precipitation observations and from the model parameters are reduced through a Monte Carlo Markov Chain (MCMC) scheme named Shuffled Complex Evolution Metropolis (SCEM) algorithm which has been extended to include a precipitation error model. Afterwards, the Bayesian Model Averaging (BMA) scheme is employed to further improve the prediction skill and uncertainty estimation using multiple model output. A series of case studies using three rainfall-runoff models to predict the streamflow in the Leaf River basin, Mississippi are used to examine the necessity and usefulness of this technique. The results suggests that ignoring either input forcings ...
This paper describes the derivation of a new equation that can be used to model the permeability behavior of a fractured, sorptive-elastic media, such as coal, under variable stress conditions commonly used during measurement of permeability data in the laboratory. The model is derived for cubic geometry under biaxial or hydrostatic confining pressures. The model is also designed to handle changes in permeability caused by adsorption and desorption of gases from the matrix blocks. The model equations can be used to calculate permeability changes caused by the production of methane from coal as well as the injection of gases, such as carbon dioxide, for sequestration in coal. Sensitivity analysis of the model found that each of the input variables can have a significant impact on the outcome of the permeability forecast as a function of changing pore pressure; thus, accurate input ...
A seaport simulationmodel, PORTSIM, has been developed for the Department of Defense (DOD) at Argonne National Laboratory. PORTSIM simulates the detailed processes of cargo loading and unloading in a seaport and provides throughput capability, resource utilization, and other important information on the bottlenecks in a seaport operation, which are crucial data in determining troop and equipment deployment capability. There are two key problems to solve in developing the HLA-compliant PORTSIM model. The first is the cargo object ownership transfer problem. In PORTSIM, cargo items, e.g. vehicles, containers, and pallets, are objects having asset attributes. Cargo comes to a seaport for loading or unloading. The ownership of a cargo object transfers from its carrier to the port and then from the port to a new carrier. Each owner of the cargo object is responsible for publishing and updating the ...
The papers collected in this volume were presented at the 4th Colloquium on Process Simulation held at Helsinki University of Technology, Espoo, Finland, June 11-13, 1997. In the more developed industrial nations, the processes for producing chemicals, energy, and materials encounter environmental concern and laws which challenge engineers to develop the processes towards more efficient, economical and safe operation. This necessitates more thorough understanding of the processes and phenomena involved. Formerly, the development of the processes was largely based on trial and error, whereas today, the development of computer performance together with the diversification of modelling software enables simulation of the processes. The increased capacity and possibilities for modelling the processes brought by the improved hardware and software, have generated a strong demand for more accurate mathematical ...
Fast Test Reactor core internal and peripheral components were assessed for flow-induced vibrational characteristics under scaled and simulated prototype flow conditions in the Hydraulic Core Mockup as an integral part of the Fast Test Reactor Vibration Program. The Hydraulic Core Mockup was an 0.285 geometric scale model of the Fast Test Reactor internals designed to simulate prototype vibrational and hydraulic characteristics. Using water to simulate sodium coolant, vibrational characteristics were measured and determined for selected model components over the scaled flow range of 36 to 110%. Additionally, in-situ shaker tests were conducted on selected Hydraulic Core Mockup outlet plenum components to establish modal characteristics. Most components exhibited resonant response at all test flow rates; however, the measured dynamic response was neither abnormal nor anomalously ...
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 ...
The CATHENA (formerly ATHENA) has been used to simulate the thermalhydraulic behaviour of the WOLSUNG-1 CANDU-600 reactor during the D_20 spill incident of 1984 November 25. A 4-inch (nominal) Liquid Relief Valve inadvertently opened in the reactor auxiliary system during normal reactor operation, resulting in a discharge of heavy water from the primary heat transport system. The valve remained open for approximately 29 minutes. CATHENA is an advanced thermalhydraulic computer code for analysis of postulated loss-of-coolant accidents (LOCA) and transient faults in CANDU nuclear reactors. A full two-fluid (six-equation) representation of the two-phase flow is used. Component models are used to represent pumps, valves, critical discharge, etc., which are necessary to describe the behaviour of the CANDU system under upset conditions. Heat transfer between the fluid and piping walls (or fuel) is modelled using applicable ...
A study was conducted to predict indoor ventilation, smoke movement and fire propagation in a new building currently under construction at the Ecole Polytechnique de Montreal. It was conducted in response to concerns regarding the impact that air quality, in normal operation conditions, and smoke concentration, in the event of a fire, may have on occupants of a building. A detailed three-dimensional model of the new building was constructed using the NURBS-based modeler, Rhino. Simulations of fire propagation and airflow ventilation were performed in different areas of the building. The flow pattern data was analyzed using the computational fluid dynamics (CFD) program called FLUENT. It predicted air flow conditions and estimated the mean age of air (MAA) in the room. Real-world geometries, such as diffuser inlets, have a significant impact on overall fluid flow behaviour and are necessary for this type of analysis. The ...
As part of a comprehensive simulation of a prototype locomotive propulsion system, a detailed model has been developed which predicts the dynamic response of an experimental 2-stroke, turbocharged and intercooled diesel engine. Engine fueling and brake torque are computed from regression equations derived from an extensive data base. Corrections are applied to the calculated steady-state torque to account for dynamic deviations of in-cylinder trapped air-fuel ratio from the steady-state value. The engine simulation accurately represents the operation of the turbo-charger, which is gear-driven at low turbocharger speeds, and freewheels through an overrunning clutch when exhaust energy accelerates the turbocharger beyond its geared speed. Engine fueling level, ie, rack, is determined from a dynamic simulation of an electrohydraulic governor which responds to the difference between the desired and the ...
With an increased use of fiber reinforced concrete in structural applications, proper characterization techniques and development of design guides are needed. This dissertation presents a multi-scale modeling approach for fiber and fabric reinforced cement-based composites. A micromechanics-based model of the yarn pullout mechanism due to the failure of the interfacial zone is presented. The effect of mechanical anchorage of transverse yarns is simulated using nonlinear spring elements. The yarn pullout mechanism was used in a meso-scale modeling approach to simulate the yarn bridging force in the crack evolution process. The tensile stress-strain response of a tension specimen that experiences distributed cracking can be simulated using a generalized finite difference approach. The stiffness degradation, tension stiffening, crack spacing evolution, and crack ...
We study the dynamical evolution of globular clusters containing primordial binaries, including full single and binary stellar evolution using our Monte Carlo cluster evolution code updated with an adaptation of the single and binary stellar evolution codes SSE and BSE from Hurley et al. We describe the modifications that we have made to the code. We present several test calculations and comparisons with existing studies to illustrate the validity of the code. We show that our code finds very good agreement with direct N-body simulations including primordial binaries and stellar evolution. We find significant differences in the evolution of the global properties of the simulated clusters using stellar evolution compared with simulations without any stellar evolution. In particular, we find that the mass loss from the stellar evolution acts as a significant energy production channel simply by reducing the total gravitational ...
Accurate crop growth monitoring and yield forecasting are significant to food security and sustainable development of agriculture. However, regional crop growth simulation faces the difficulties in determining the spatial distribution of some model parameters and initial conditions. In this study, regional biomasses at turn-green stage of winter wheat were re-estimated by linking WOFOST model and Soil Adjusted Vegetation Index (SAVI) synthesized from remote sensing data. Moreover, we proposed a way of combining evapotranspiration derived from satellite remote sensing data to crop grow simulationmodel. Thus, the regional initial available soil water and irrigation at earring stage were re-initialized and re-estimated by using remote sensing data. Those methods were well applied to simulate the growth and development for winter wheat at local site. After ...
This paper builds upon previous work [Sprigg and Ehlen, 2004] by introducing a bond market into a model of production and employment. The previous paper described an economy in which households choose whether to enter the labor and product markets based on wages and prices. Firms experiment with prices and employment levels to maximize their profits. We developed agent-based simulations using Aspen, a powerful economic modeling tool developed at Sandia, to demonstrate that multiple-firm economies converge toward the competitive equilibria typified by lower prices and higher output and employment, but also suffer from market noise stemming from consumer churn. In this paper we introduce a bond market as a mechanism for household savings. We simulate an economy of continuous overlapping generations in which each household grows older in the course of the simulation and continually ...
Glacier fluctuations exclusively due to internal variations in the climate system are simulated using downscaled integrations of the ECHAM4/OPYC coupled general circulation model (GCM). A process-based modeling approach using a mass balance model of intermediate complexity and a dynamic ice flow model considering simple shearing flow and sliding are applied. Multi-millennia records of glacier length fluctuations for Nigardsbreen (Norway) and Rhonegletscher (Switzerland) are simulated using auto-regressive processes determined by statistically downscaled GCM experiments. Return periods and probabilities of specific glacier length changes using GCM integrations excluding external forcings such as solar irradiation changes, volcanic or anthropogenic effects are analyzed and compared to historical glacier length records. Preindustrial fluctuations of the glaciers as ...
A force measurement technique has been developed for large-scale aerodynamic models with a short test time. The technique is based on direct acceleration measurements, with miniature accelerometers mounted on a test model suspended by wires. Measuring acceleration at two different locations, the technique can eliminate oscillations from natural vibration of the model. The technique was used for drag force measurements on a 3 m long supersonic combustor model in the HIEST free-piston driven shock tunnel. A time resolution of 350 #mu#s is guaranteed during measurements, whose resolution is enough for ms order test time in HIEST. To evaluate measurement reliability and accuracy, measured values were compared with results from a three-dimensional Navier-Stokes numerical simulation. The difference between measured values and numerical simulation values was less than ...
The design considerations and computational fluid dynamics (CFD) modeling of a high efficiency, low emissions, hydrogen-fueled engine for use as the prime mover of a series hybrid automobile is described. The series hybrid automobile uses the engine to generate electrical energy via a lightweight generator, the electrical energy is stored in a power peaking device (like a flywheel or ultracapacitor) and used as required to meet the tractive drive requirements (plus accessory loads) through an electrical motor. The engine/generator is stopped whenever the energy storage device is fully charged. Engine power output required was determined with a vehicle simulation code to be 15 to 20 kW steady state with peak output of 40 to 45 kW for hill climb. Combustion chamber and engine geometry were determined from a critical review of the hydrogen engine experiments in the literature combined with a simplified global engine model. Two ...
The aim of the work in this thesis is to investigate the convective and diffusive transport in the TM chemistry transport model, and to investigate some aspects of the consequences for NOx. The large inaccuracy and uncertainty in the description of processes like convection and turbulent diffusion, the strong dependence of the radiative forcing of ozone on its vertical distribution, and the strong dependence of the ozone production on the distribution of NOx, are the main motivation. The availability of the ERA-40 data, where convective data and vertical diffusion coefficients are archived, allows a study of the effect of different convective mass flux sets, and different vertical diffusion coefficients on the model-simulated distribution of tracers. In this thesis the following questions are addressed : (1) How large is the sensitivity of the (modelsimulated) distribution of ozone and nitrogen oxides ...
The objective of the current MITRE study is to estimate the cost of the coal derived liquid products that result from processing low rank coals in the Wilsonville and HRI CTSL configurations. The resultant cost of liquids is then compared to liquids produced from bituminous coals. The methodology used to accomplish the above objective was to utilize the MITRE coal liquefaction cost model to assess the overall technical and economic performance of low rank coals in both the Wilsonville and HRI configurations. Actual test data from Wilsonville and HRI was used as the basis for these analyses. This test data is modified to simulate complete resid recycle extinction and to standardize the operating performance of the deashing unit. This corrected yield data is then used in the simulationmodels to predict the overall performance and economics of a conceptual commercial coal liquefaction facility. 9 refs., 7 ...
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 ...
A dense Pd-Ag membrane reactor (MR) with 100% hydrogen selectivity packed with either Rh/La2O3 or Rh/La2O3-SiO2 as catalysts was used to carry out the dry reforming of methane. The membrane reactor simulation was performed using a well-known reactor model. For this purpose, we employed the equations derived from complete kinetic studies of the dry reforming of methane reaction in connection with both catalysts. In addition, we developed the kinetic equation for the reverse water gas shift reaction (RWGS). The combination of detailed kinetic studies with the measured permeation flux for the Pd-Ag membrane allowed a complete comparison between experimental and simulated operation variables. The variables studied for both catalysts were methane conversion and hydrogen permeation as a function...
In this paper we study theoretically a pulsed corona discharge plasma in air at atmospheric pressure. A self consistent one-dimensional model, using a flux-corrected transport numerical scheme, is presented to achieve it. The charged particle kinetics in this model are described by one-dimensional continuity equations coupled with Poisson's equations. The spatio-temporal local field, charge density variations and the velocity profile of an ionizing front are calculated to describe the dynamical behavior of corona discharge plasmas. The simulation results of a wire-in-cylinder corona discharge plasma explain the physical mechanisms of discharge processes. These results may also be apply to obtain the optimizing parameters for designing the plasma reactor.
Nano photonic materials are synthetically manufactured crystals at the nano scale with the target of creating a microstructure with a special electro-magnetic periodicity. Such nano photonic materials have the ability to control light propagation and thus are capable of creating photonic bandgaps in the frequency domain. We propose using nano photonic crystals as sensors to detect microdamage in composite materials. We demonstrate using a simulationmodel that a nano photonic sensor attached to a composite bar experiences a significant change in its bandgap profile when damage is induced in the composite bar. The model predicts the frequency response of the nano photonic sensor using the transfer matrix method. A damage metric to evaluate the change in the frequency response is developed. Successful developments of nano photonic sensors allow damage identification at scales not attainable using current sensing technologies.
Analytical and numerical methods have been used in this investigation to model the behavior of geothermal systems under exploitation. The work is divided into three parts: (1) development of a numerical code, (2) theoretical studies of geothermal systems, and (3) field applications. A new single-phase three-dimensional simulator, capable of solving heat and mass flow problems in a saturated, heterogeneous porous or fractured medium has been developed. The simulator uses the integrated finite difference method for formulating the governing equations and an efficient sparse solver for the solution of the linearized equations. In the theoretical studies, various reservoir engineering problems have been examined. These include (a) well-test analysis, (b) exploitation strategies, (c) injection into fractured rocks, and (d) fault-charged geothermal reservoirs.
The Linac Coherent Light Source (LCLS) will be the world's first x-ray free-electron laser (FEL). To ensure the vitality of FEL lasing, it is critical to preserve the high quality of the electron beam during acceleration and compression. The peak current and final energy are very sensitive to system jitter. To minimize this sensitivity, a longitudinal feedback system on the bunch length and energy is required, together with other diagnostics and feedback systems (e.g., on transverse phase space). Here, we describe a simulation framework, which includes a realistic jitter model for the LCLS accelerator system, the RF acceleration, structure wakefield, and second order optics. Simulation results show that to meet the tight requirements set by the FEL, such a longitudinal feedback system is mandatory.
Data from Monte Carlo transport codes are used to model radiobiological effects. We previously reported the Fourier analysis of ionization data generated by simulating a 500-keV proton traversing water. Here, we extend Fourier analysis to energy transfer data of another radiation type, a 14-MeV #alpha#-particle. A radiobiological model based on this frequency-domain analysis views cell as an information processing system . It lends itself naturally to traditional engineering analyses. One engineering principle-the output response of a linear system to random signal-is applied here to explain the fact that there is measurable difference in the magnitude of the biological effectiveness when a given biological system is irradiated with two different radiation types of the same Linear Energy Transfer (LET).
To form a licensing basis for a new methodology for a fuel channel safety analysis code for CANDU-6 nuclear reactor, a CATHENA model for a post-blowdown fuel channel analysis has been developed, and tested for a high temperature thermal-chemical experiment CS28-1 [Lei, Q.M., 1993. Post-test analysis of the 28-element high-temperature thermal-chemical experiment CS28-1. In: 4th International Conference on Simulation Methods in Nuclear Engineering, Montreal, PQ, 1993]. Pursuant to the objective of this investigation, the current study has focused on understanding the involved phenomena, their interrelations, and how to maintain a good accuracy of the temperature and H2 generation rate prediction without losing the important physics of the involved phenomena. The transient simulation results ...
An analytical model is proposed for the design and simulation of in-situ fracture of deep oil sand formations. This model is based on the finite element variational principle in conjunction with special empirical modules to characterize in-situ oil sands behavior. A computer code by the name of SANFRAC was developed to handle the dynamic fracture of formations induced by explosives. Simulation of hydraulic fracture processes can be treated by the same code as special cases using the quasi-static analysis. Numerical case studies by the SANFRAC code indicate that extensive horizontal fracture can be achieved by dynamic loads with proper fracture starters configured at the injection well. The unique advantage of the dynamic fracturing technique over the hydraulic fracture methods is also demonstrated by these studies.
When backward time travel through wormholes is taken into account, classical physics loses its determinism and allows simulation of some quantum behaviours. We show how it is possible to simulate a non-local wavefunction reduction-type effect, i.e. we present a mechanical analogy for the collapse of the wavefunction of an entangled state of two removed particles. This situation can be seen as the simplest EPR situation, i.e. the situation where there is just one direction to measure along the spin (or the correlated properties). We present no rigorous results here, just a different point of view about something that is generally thought to be impossible: modelling a quantum indeterministic and non-local behaviour with a mechanical system.
Several complex organic molecules are routinely detected in high abundances towards hot cores and hot corinos. For many of them, their paths of formation in space are uncertain, as gas phase reactions alone seem to be insufficient. In this paper, we investigate a possible solid-phase route of formation for methyl formate (HCOOCH3). We use a chemical model updated with recent results from an experiment where simulated grain surfaces were irradiated with 200 keV protons at 16 K, to simulate the effects of cosmic ray irradiation on grain surfaces. We find that this model may be sufficient to reproduce the observed methyl formate in dark clouds, but not that found in hot cores and corinos.
Motion Camouflage (MC) is illuminated as a novel strategy in counterattacking anti-satellite by way of stealth trajectory scheduling. The dynamics model of MC in space is developed and a quadratic function with three boundary constraints is employed for trajectory determination. Based on the model a scenario is set to run the simulation. The results indicate given the designed acceleration input, the predator will be moved following a prescribed route, which precisely locates the predator between two objects at each time instant. In the last approaching phase, the motion is achieved with a big bumping rate which guarantees the power of this striking. Methods for deriving minimum fuel cost in the fixed approaching duration and the minimum approaching duration in limited acceleration input are proposed and are verified in the simulation. At last, camouflage is recognized as a multi-faceted affair, in ...
Event generators that handle neutrino-nucleon interaction have been developed for the FLUKA code [1]. In earlier FLUKA versions only quasi-elastic (QEL) interactions were included, and the code relied on external event generators for the resonance (RES) and deep inelastic scattering (DIS). The new DIS+RES event generator is fully integrated in FLUKA and uses the same hadronization routines as those used for simulating hadron-nucleon interactions. Nuclear effects in neutrino-nucleus interactions are simulated within the same framework as in the FLUKA hadron-nucleus interaction model (PEANUT), thus profiting from its detailed physics modelling and longstanding benchmarking. The generators are available in the standard FLUKA distribution. They are presently under development and several improvements are planned to be implemented. The physics relevant to the neutrino-nucleon interactions and the results of ...
The research documented herein was funded by a research contract between the Research Reactors Division (RRD) of Oak Ridge National Laboratory (ORNL) and the University of Tennessee, Knoxville (UTK) Mechanical, Aerospace and Biomedical Engineering Department (MABE). The research was governed by a statement of work (SOW) which clearly defines nine specific tasks. This report is outlined to follow and document the results of each of these nine specific tasks. The primary goal of this phase of the research is to demonstrate, through verification and validation methods, that COMSOL is a viable simulation tool for thermal-hydraulic modeling of the High Flux Isotope Reactor (HFIR) core. A secondary goal of this two-dimensional phase of the research is to establish methodology and data base libraries that are also needed in the full three-dimensional COMSOL simulation to follow. COMSOL version 3.5a was used for all of the ...
High-energy ionizing radiation is a prominent modality for the treatment of many cancers. The approaches to electron dose calculation can be categorized into semi-empirical models (e.g. Fermi-Eyges, convolution-superposition) and probabilistic methods (e.g. Monte Carlo). A third approach to dose calculation has only recently attracted attention in the medical physics community. This approach is based on the deterministic kinetic equations of radiative transfer. We derive a macroscopic partial differential equation model for electron transport in tissue. This model involves an angular closure in the phase space. It is exact for the free streaming and the isotropic regime. We solve it numerically by a newly developed HLLC scheme based on Berthon et al (2007 J. Sci. Comput. 31 347-89) that exactly preserves the key properties of the analytical solution on the discrete level. We discuss several test cases taken from the medical ...
This investigation presents scale-up of a detailed heterogeneous geostatistical model to a full field reservoir simulationmodel, considering both single and two-phase flow properties. The model represents a typical low permeability Danish North Sea chalk reservoir and includes capillary pressure and saturation end-point variations. Two new up-scaling methods has been investigated, all based on fine scale simulation on a cross section of the geomodel. The first methods assumes piston style behaviour and a coupled viscosity is introduced into the basic Darcy`s equations. The second method is a modification of the JBN method traditionally applied in analysing results from core flooding experiments, which emerged as the most successful and therefore also the recommended method. 1. In addition to the up scaling work we review the Equivalent Radius Method for capillary pressure ...
We create realistic, full-sky, half-arcminute resolution simulations of the microwave sky matched to the most recent astrophysical observations. The primary purpose of these simulations is to test the data reduction pipeline for the Atacama Cosmology Telescope (ACT) experiment; however, we have widened the frequency coverage beyond the ACT bands and utilized the easily accessible HEALPix map format to make these simulations applicable to other current and near future microwave background experiments. Some of the novel features of these simulations are that the radio and infrared galaxy populations are correlated with the galaxy cluster and group populations, the primordial microwave background is lensed by the dark matter structure in the simulation via a ray-tracing code, the contribution to the thermal and kinetic Sunyaev-Zel'dovich (SZ) signals from galaxy clusters, ...
The SCDAP/RELAP5 code has been developed for best estimate transient simulation of light water reactor coolant systems during a severe accident. The code models the coupled behavior of the reactor coolant system, core, fission product released during a severe accident transient as well as large and small break loss of coolant accidents, operational transients such as anticipated transient without SCRAM, loss of off-site power, loss of feedwater, and loss of flow. A generic modeling approach is used that permits as much of a particular system to be modeled as necessary. Control system and secondary system components are included to permit modeling of plant controls, turbines, condensers, and secondary feedwater conditioning systems. This volume describes the organization and manner of the interface between severe accident models which are resident in the SCDAP ...
of satellite, aircraft and ground-based observations. In ..... swarm traps by Utah Department of Food and Agriculture (DAF) personnel. Preliminary Results ..... the 150th anniversary of the book's publication and the 200th anniversary of Darwin's birth. ...... cold war submarine missions to find that Arctic Ocean ...
BackgroundDinoflagellates are unicellular, often photosynthetic protists that play a major role in the dynamics of the Earth's oceans and climate. Sequencing of dinoflagellate nuclear...Full Text Available
Meso- and submesoscales (fronts, eddies, filaments) in surface ocean flow have a crucial influence on marine ecosystems. Their dynamics partly control the foraging behavior and the displacement of marine...Full Text Available
The ratio of particle velocity to translation speed exceeds 1 for 98% of .... The westward movement of the meandering jet and eddy system transports nutrient- rich ...... J. Hydrometeor., 4, 1044-1066. Gedney, N., P. M. Cox, H. Douville, ...
BackgroundIt has been reported that the El Niño–Southern Oscillation (ENSO) influences the interannual variation of endemic cholera in Bangladesh. There is increased...Full Text Available
Separate abstracts were prepared for 60 papers in this volume of conference proceedings. None other papers in the proceedings were previously abstracted for EDB. Also included in the proceedings but not indexed separately are six workshop summaries and a session introduction. (LEW)
Computer programs have been developed to define the temperature increase which would be needed to bring deep-ocean water into density equilibrium with surface water for locations where data are available. A series of continuous-flow studies on phytoplankton blooms resulting from mixtures of 80 percent deep and 20 percent surface water in 2000-liter concrete culturing vessels (''reactors'') has been completed. A quantitative determination of nutrient utilization and flow through a combined primary and secondary trophic level system has been completed. This study utilized the clam Tapes semidecussata, fed from phytoplankton grown in 80 percent deep and 20 percent surface water. An analysis of the fate of the deep water discharged from a floating OTEC plant indicates that horizontal containment of the resulting deep water: surface water mixture is necessary if conditions optimal for open-sea mariculture are to obtain. The design of ...
In the IEEE Marine Technology Society OCEANS 2009 Conference, Biloxi. October 26-29, 2009. Biloxi, MS. Al-Hamdan, M.; Estes, M.; Quattrochi, D.; Thom, R.; ...
provides global ocean-surface wind fields with more structures than numerical weather prediction (Geophys. Res. Lett., Vol 25, No. 6, 761-764). The NSCAT data sets was produced...
Records of past climates contained in ice cores, ocean sediments, and other archives show that large, abrupt, widespread climate changes have occurred repeatedly in the past. These changes were especially...Full Text Available
Foreward: Preface; Introduction; The natural science of global climate change; Land and water use; Coastal zones and oceans; Energy and industry; Energy and social systems; Technological change; and Sponsoring organizations, International Advisory Board, and project participants.
It is now widely accepted that novel infectious disease can be a leading cause of serious population decline and even outright extinction in some invertebrate and vertebrate groups (e.g., amphibians)....Full Text Available
Harmful algal blooms, caused by rapid growth and accumulation of certain microalgae in the ocean, pose considerable impacts on marine environments, aquatic industries and even public health. Here, we...Full Text Available
One tiny marine plant makes life on Earth possible: phytoplankton. These microscopic photosynthetic drifters form the basis of the marine food web, they regulate carbon in the atmosphere, and are responsible for half of the photosynthesis that takes place on this planet.
Adenosine triphosphate (ATP) is an indicator of living biomass in marine particulates. This report details the method used by Lawrence Berkeley Laboratory to analyze particulate ATP in samples taken from oligotrophic, tropical ocean waters. It represents a synthesis of previously published methods.
BackgroundOceans are iron-deficient and nutrient-poor environments. These conditions impart limitations on our understanding of and our ability to identify microorganisms from the...Full Text Available
Picocyanobacteria of the genus Synechococcus are important contributors to marine primary production and are ubiquitous in the world's oceans. This genus is genetically diverse, and...Full Text Available
In the groundwater flow modelling part of this work the effective dilution volume in the well scenario was estimated by means of transient simulations of groundwater flow and transport, which are coupled due to the varying salinity. Both deep, drilled wells and shallow surface wells in the vicinity of the repository were considered. The simulations covered the time period from the present to 1000 years after the present. Conceptually the fractured bedrock consists of planar fracture zones (with a high fracture density and a greater ability to conduct water) and the intact rock (in which the fracture density and the hydraulic conductivity are low). For them the equivalent-continuum model was applied separately. Thus, the fractured bedrock was considered as piecewise homogeneous (except for the depth dependence) and isotropic continuum with representative average characteristics. A generic ...
This paper describes the modeling of horizontal steam generator with the TRACE code and calculation results of a loss-of-feedwater (LOF-10) experiment at the PACTEL facility. Parallel Channel Test Loop (PACTEL) is an integral test facility for a VVER-440 type nuclear reactor. The main objectives were to prepare a simulationmodel for its horizontal steam generator with the TRACE thermal hydraulic code and assess different modeling options of the code. PACTEL experiment LOF-10 was chosen for this assessment. The calculation results showed that TRACE is capable in simulating horizontal steam generator behavior both in steady state and during loss-of-feedwater transient. The phenomenon of heat transfer from primary to secondary side, steam superheating and flow reversal in the lowest heat exchange tubes were studied in detail. Different nodalization options were introduced. In the ...
This paper describes the modeling of horizontal steam generator with the TRACE code and calculation results of a loss-of-feedwater (LOF-10) experiment at the PACTEL facility. Parallel Channel Test Loop (PACTEL) is an integral test facility for a VVER-440 type nuclear reactor. The main objectives were to prepare a simulationmodel for its horizontal steam generator with the TRACE thermal hydraulic code and assess different modeling options of the code. PACTEL experiment LOF-10 was chosen for this assessment. The calculation results showed that TRACE is capable in simulating horizontal steam generator behavior both in steady state and during loss-of-feedwater transient. The phenomenon of heat transfer from primary to secondary side, steam superheating and flow reversal in the lowest heat exchange tubes were studied in detail. Different nodalization options were introduced. In the ...
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 simulationmodel 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. ...
A scaling Monte Carlo method has been developed to calculate diffuse reflectance from multilayered media with a wide range of optical properties in the ultraviolet-visible wavelength range. This multilayered scaling method employs the photon trajectory information generated from a single baseline Monte Carlo simulation of a homogeneous medium to scale the exit distance and exit weight of photons for a new set of optical properties in the multilayered medium. The scaling method is particularly suited to simulating diffuse reflectance spectra or creating a Monte Carlo database to extract optical properties of layered media, both of which are demonstrated in this paper. Particularly, it was found that the root-mean-square error (RMSE) between scaled diffuse reflectance, for which the anisotropy factor and refractive index in the baseline simulation were, respectively, 0.9 and 1.338, and independently ...
Gamma scanning measurements are used to determine the radioactive nuclides from an irradiated fuel rod in a hot cell which is a routine task in the Irradiated Materials Examination Facility (IMEF) of the Korea Atomic Energy Research Institute (KAERI), Total gamma spectra of radioactive nuclides and the atomic ratio of {sup 134}Cs to {sup 137}Cs in a few segments of a single fuel rod are measured by the present gamma scanning system. For the purpose of determining the radioactive nuclide's distribution in detail within a fuel rod quantitatively, we attempt to upgrade the present gamma scanning system in this paper, therefore, an investigation of the feasibility of the tomographic gamma scanning (TGS) technique, which is one of the promising nondestructive assay (NDA) methods is proposed. The TGS technique is an efficient method and it has an accurate precision for the characterization of several fuel parameters. In this paper, a principle and an emission TGS and an iterative ...
For the performance analysis of Canadian deuterium uranium (CANDU) reactor shutdown system number 2 (SDS2), a computational fluid dynamics model of poison jet flow has been developed to estimate the flow field and poison concentration formed inside the CANDU reactor calandria. As the ratio of calandria shell radius over injection nozzle hole diameter is so large (1055), it is impractical to develop a full-size model encompassing the whole calandria shell. In order to reduce the model to a manageable size, a quarter of one-pitch length segment of the shell was modeled using symmetric nature of the jet; and the injected jet was treated as a source term to avoid the modeling difficulty caused by the big difference of the hole sizes. For the analysis of an actual CANDU-6 SDS2 poison injection, the grid structure was determined based on the results of two-dimensional real- and source-jet ...
Currently, x-ray mammography is the method of choice in breast cancer screening programmes. As the mammography technology moves from 2D imaging modalities to 3D, conventional computational phantoms do not have sufficient detail to support the studies of these advanced imaging systems. Studies of these 3D imaging systems call for a realistic and sophisticated computational model of the breast. DeBRa (Detailed Breast model for Radiological studies) is the most advanced, detailed, 3D computational model of the breast developed recently for breast imaging studies. A DeBRa phantom can be constructed to model a compressed breast, as in film/screen, digital mammography and digital breast tomosynthesis studies, or a non-compressed breast as in positron emission mammography and breast CT studies. Both the cranial-caudal and mediolateral oblique views can be modelled. The anatomical details ...
Recent efforts to identify methods of modeling pitting corrosion damage of high-level radioactive-waste containers are described. The need to develop models that can provide information useful to higher level system performance assessment models is emphasized, and examples of how this could be accomplished are described. Work to date has focused upon physically-based phenomenological stochastic models of pit initiation and growth. These models may provide a way to distill information from mechanistic theories in a way that provides the necessary information to the less detailed performance assessment models. Monte Carlo implementations of the stochastic theory have resulted in simulations that are, at least qualitatively, consistent with a wide variety of experimental data. The effects of environment on pitting corrosion have been included ...
In this paper, we evaluate the performance of the 1- and 5-site models of methane on the description of adsorption on graphite surfaces and in graphitic slit pores. These models have been known to perform well in the description of the fluid-phase behavior and vapor-liquid equilibria. Their performance in adsorption is evaluated in this work for nonporous graphitized thermal carbon black, and simulation results are compared with the experimental data of Avgul and Kiselev (Chemistry and Physics of Carbon; Dekker: New York, 1970; Vol. 6, p 1). On this nonporous surface, it is found that these models perform as well on isotherms at various temperatures as they do on the experimental isosteric heat for adsorption on a graphite surface. They are then tested for their performance in predicting the adsorption isotherms in graphitic slit pores, in which we would like to explore the effect of confinement on the ...
Customer Relationship Management (CRM) systems can help organizations manage customer interactions more effectively. Like many new technologies, CRM has been accompanied by vendor hype and stories of implementation failure. Work on critical success factors (CSFs) should encourage more appropriate implementation practice; however many CSF studies conclude with a list of factors but provide little further guidance. In particular, there is a need for stronger theoretical models of the entire CRM innovation process which can be used by managers to understand better the underlying causes of success and failure. This paper adopts a novel approach to this problem by firstly developing a conceptual model of CRM innovation and then converting this model into a dynamic simulationmodel. Some early s...
A new methodology for three-dimensional (3-D) numerical modeling of the continuous steel casting process is presented and evaluated. The results are compared with available experimental data obtained in a water-oil mold model for various flow rates and immersion depths of the submerged entry nozzle. The water-oil interface pattern at the top of the mold is also investigated. Wave formation and evolution in time is simulated with a volume tracking method. The predicted flow field inside the mold exhibits quite satisfactory agreement with the corresponding measurements, whereas some quantitative differences are observed referring to the final wave amplitude at low casting speeds. Conclusions of scientific and engineering importance concerning both the computer algorithm and the mold performance are also drawn.
In this study the results of mathematical simulation and experimental research of texture transition in stress assisted shear phase transformations are discussed. The orientation distribution function (ODF) calculation model to the any type crystalline lattice symmetry materials has been suggested. By crystalline lattice symmetry analysis this computer program determines necessary amount of the initial experimental data. The experimental verification supports this model well. The ODF calculation model after shear phase transformation with equal use of all possible variants of orientation relationship or with various degrees of variant selection worked out. In this study the experimental and calculating data of the direct and reversal shear transformations #alpha# <=> #gamma# in Fe-alloys and B2 <=> B19' in Ti-Ni alloys are discussed. The stress assisted shear phase transformation ...
Although managing fisheries is complex, policymakers must make decisions daily that affect the future of aquatic resources. To make these decisions, they often combine a mental model with lessons learned from formal models such as computer models. While such computer models are frequently interpreted as predictors of the future, they are often more valuable when used as tools for learning about fisheries and options for management. This approach, looking at computer simulations as learning laboratories, is what this study is all about since its objective is to investigate by means of a microworld if policymakers in Belgian fisheries have the policy instruments at hand to align the real fisheries world with their often completely different desired fisheries worlds. This study illustrates ba...
Objectives The optimal timing of pyeloplasty for children diagnosed with ureteropelvic junction obstruction (UPJO) after workup for antenatal hydronephrosis is disputed. We sought to examine the potential costs and clinical outcomes of treatment protocols featuring different indications for pediatric pyeloplasty using Markov models. Methods Cost and outcomes analysis using Markov modeling was performed for three treatment algorithms: medical management, immediate pyeloplasty (during the first year of life), and pyeloplasty after no improvement on imaging. The costs were determined from the perspective of the medical institution. The variables tracked during Markov modelsimulation included age at resolution of UPJO, the proportion of patients with worsened hydronephrosis, the number of pye...
Although state estimation is a mature technique which is widely available in the EMS industry, experience with practical implementation and use of state estimation and network security analysis functions at electric power utilities over the last two decades indicates recurring difficulties and problems attributable to inadequate external network modeling detail and data. This paper addresses the development of a set of guidelines for external network modeling and data exchange, based on the results of a recent project sponsored by EPRI. A general methodology is developed based on the results of a survey of a representative set of utilities and EMS suppliers, supplemented by subsequent analysis and simulation studies. Distinction is made between guidelines pertaining to external network topology, analog measurements, data exchange, and implementation procedures. A philosophy and approach for constructing and testing external ...
There are two difficulties encountered in modeling valveless micropumps using lumped-element methods. The pressure loss coefficient for fluidic diodes used in valveless pumps to rectify flow depends on the flow direction. A problem arises in choosing the proper loss correlation because the flow direction is not known a priori. Another problem is the quadratic form of the equation for the flow through the fluidic diodes, which brings about multiple solutions. The above problems become even more serious in multi-chamber cases. They are overcome in this study by suitably formulating the flow resistance. In addition, the flow inertia is accounted for in the unsteady model. The steady and unsteady models are evaluated by comparing with CFD simulations, which also serve to illustrate the flow fi...
A mathematical model for heat transfer during solidification in continuous casting of automobile steel, was established on researching under the influence of the solidifying process of bloom quality of CCM in the EAF steelmaking shop, at Shijiazhuang Iron and Steel Co. Ltd. Several steel grades were chosen to research, such as, 40Cr and 42CrMo. According to the results of the high temperature mechanical property tests of blooms, the respective temperature curves for controlling the solidification of different steels were acquired, and a simulating software was developed. The model was verified using two methods, which were bloom pinshooting and surface strand temperature measuring experiments. The model provided references for research on the solidifying process and optimization of a secon...
Understanding the impacts of climate change on viticulture is especially essential in those areas producing high-quality wines. In this work, we create an operational framework to investigate climate change impact on viticulture in the Tuscany region (central Italy) the viticulture industry of which relies on producing high-quality wines to compete in a global market. The framework includes (i) statistical downscaling of General Circulation Model (GCM) outputs for the period 1975?2099 to a local scale; (ii) the use of downscaling outputs as driving variables in specific simulationmodels; (iii) the spatial interpolation of model outputs to feed an economic and (iv) a quality model. The results show that as a consequence of a progressive increase in temperature and a decrease in rainfall, (...
Dispersed flow film boiling is the heat transfer regime that occurs at high void fractions in a heated channel. The way this transfer mode is modelled in the NRC computer codes (RELAP5 and TRAC) and the validity of the assumption and empirical correlations used is discussed. An extensive review of the theoretical and experimental work related with heat transfer to highly dispersed mixtures reveals the basic deficiencies of these models: the investigation refers mostly to the typical conditions of low rate bottom reflooding, since the simulation of this physical situation by the computer codes has often showed poor results. The alternative models that are available in the literature are reviewed, and their merits and limits are highlighted. The modification that could improve the physics of the models implemented in the codes are identified. (author) 13 figs., 123 refs.
The present paper concerns the development of a polymer electrolyte membrane (PEM) fuel cells model in order to analyse the behaviour of the cell varying its operating parameters such as temperature (70-120 ^oC), pressure (1-3 atm), relative humidity (35-100%) and CO content (0-200 ppm) in the feeding gas. The voltage output and the relative current of a generic PEM fuel cell has been predicted through a mathematical method based on semi-empirical correlations derived from the experimental data available in literature. The model has been implemented in a commercial process simulator (AspenPlus) and it has been validated by comparing the obtained results with further experimental data (not utilized for the model development). The comparison has underlined a good agreement with the experimen...
The bulk removal of CO/sub 2/ from industrial gases is a frequently applied process in the field of gas-treating. Usually chemical solvents are used in order to increase both the solvent-loading capacity and the absorption rate. These solvents often are aqueous solutions of alkanolamines. The loaded-liquid is regenerated at high temperatures with the aid steam. In the present study absorption rates are experimentally determined in various mixtures of amines over a wide range of concentrations. Simultaneously, a model was developed for the description of mass transfer accompanied by parallel reversible chemical reactions. The equations of this mass transfer model could not be solved analytically and therefore numerical techniques were used to obtain an exact description of this process. This model also includes the asymptotic situations of homogeneous catalysts and shuttle mechanism for which in the literature several ...