Energy Technology Data Exchange (ETDEWEB)

A stable, accurate, and efficient implementation of MacCormack's explicit algorithm for the Parabolized Navier-Stokes equations is demonstrated. The familiar problem of decoding the conservative axial flux vector is solved, resulting in accurate, smooth dependent variable profiles through the viscous-layer sonic line. Source terms due to transformation of the parabolized governing equations into the computational plane and the equations into the computational plane and the resulting metric differencing have been identified and eliminated through inclusion of appropriate geometric conservation law terms. Test cases computed include two- and three-dimensional supersonic and hypersonic flow at laminar and turbulent Reynolds numbers. The computed results demonstrate very good agreement with experiment and with solutions of the full Navier-Stokes equations. Computational times required for the MacCormack explicit PNS code are approximately ...

Energy Technology Data Exchange (ETDEWEB)

Staggered arrays of dimples printed on opposite surfaces of a cooling channel is formulated numerically and optimized with hybrid multi-objective evolutionary algorithm and Pareto optimal front. As Pareto optimal front produces a set of optimal solutions, the trends of objective functions with design variables are predicted by hybrid multi-objective evolutionary algorithm. The problem is defined by three non-dimensional geometric design variables composed of dimpled channel height, dimple print diameter, dimple spacing, and dimple depth, to maximize heat transfer rate compromising with pressure drop. Twenty designs generated by Latin hypercube sampling were evaluated by Reynolds-averaged Navier-Stokes solver and the evaluated objectives were used to construct Pareto optimal front through hybrid multi-objective evolutionary algorithm. The optimum designs were grouped by k-means clustering technique and some of the clustered points were evaluated ...

Energy Technology Data Exchange (ETDEWEB)

Uranium enrichment by gas centrifugation relies on the pressure diffusion of the isotopes when their gaseous mixture is submitted to the pressure diffusion in the strong radial pressure gradient due to the centrifugal force in a very rapidly rotating cylinder. To compute the isotope separation in a gas centrifuge, the viscous compressible Navier-Stokes equations are solved for the mixture, then the diffusion equations are solved for the isotopes. As far as the separation performance is concerned, the models deal with the steady state of the axi-symmetric gas motion in the rotating cylinder. The numerical tools are briefly presented: first, the MOLINA code which numerically solves the axi-symmetric gas motion in the specific conditions of rapid rotation by use of a change of variables suited to the high density gradient and by use of a Marker And Cell method; then, the diffusion equations solver which uses a Finite Volume method with a mesh ...

CERN Document Server

A well-known diffuse interface model consists of the Navier-Stokes equations nonlinearly coupled with a convective Cahn-Hilliard type equation. This system describes the evolution of an incompressible isothermal mixture of binary-fluids and it has been investigated by many authors. Here we consider a variant of this model where the standard Cahn-Hilliard equation is replaced by its nonlocal version. More precisely, the gradient term in the free energy functional is replaced by a spatial con- volution operator acting on the order parameter phi. Therefore the coupling with the Navier-Stokes equations is difficult to handle even in two spatial dimensions because of the lack of regularity of phi. We establish the global existence of a weak solution.

CERN Document Server

The Cahn-Hilliard-Navier-Stokes system is based on a well-known diffuse interface model and describes the evolution of an incompressible isothermal mixture of binary fluids. A nonlocal variant consists of the Navier-Stokes equations suitably coupled with a nonlocal Cahn-Hilliard equation. The authors, jointly with P. Colli, have already proven the existence of a global weak solution to a nonlocal Cahn-Hilliard-Navier-Stokes system subject to no-slip and no-flux boundary conditions. Uniqueness is still an open issue even in dimension two. However, in this case, the energy identity holds. This property is exploited here to define, following J.M. Ball's approach, a generalized semiflow which has a global attractor. Through a similar argument, we can also show the existence of a (connected) global attractor for the convective nonlocal Cahn-Hilliard equation with a given velocity field, even in dimension three. Finally, we ...

Energy Technology Data Exchange (ETDEWEB)

The solution of the governing steady transport equations for momentum, heat and mass transfer in flowing fluids can be very difficult. These difficulties arise from the nonlinear, coupled, nonsymmetric nature of the system of algebraic equations that results from spatial discretization of the PDEs. In this manuscript the authors focus on evaluating a proposed nonlinear solution method based on an inexact Newton method with backtracking. In this context they use a particular spatial discretization based on a pressure stabilized Petrov-Galerkin finite element formulation of the low Mach number Navier-Stokes equations with heat and mass transport. The discussion considers computational efficiency, robustness and some implementation issues related to the proposed nonlinear solution scheme. Computational results are presented for several challenging CFD benchmark problems as well as two large scale 3D flow simulations.

British Library Electronic Table of Contents (United Kingdom)

A method for transforming the Euler and Navier-Stokes equations and a complete system of fluid dynamics equations in three dimensions to a closed system on any moving surface is proposed. As a result, for an arbitrary geometric configuration, the dimension of the equations is reduced by one, which makes them convenient for numerical simulation. The general principles of the method are described, and verifying examples are presented.

Energy Technology Data Exchange (ETDEWEB)

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 simulations from the geometric stage to complete ...

Energy Technology Data Exchange (ETDEWEB)

Different calculated models are presented for turbulent break-off flows and their classification, reflecting the sequence of historical development. The study was done based on equations of viscous liquid of the Navier-Stokes type with development of special phenomenological models of turbulence which take into consideration real properties of the break-off flows based on simpler models of flow presented in the work. In order to calculate two-dimensional turbulent flows, a method of viscous-nonviscous interaction is used. It employs numerical solutions for nonviscous flow and integrated methods of calculating the dissipative region. This method can be extended for calculating the transonic break-off flows, and also break-off in an incompressible fluid when there is cavitation.

Energy Technology Data Exchange (ETDEWEB)

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

Energy Technology Data Exchange (ETDEWEB)

This paper is an introduction course in modelling turbulent thermohydraulics, aimed at computational fluid dynamics users. No specific knowledge other than the Navier Stokes equations is required beforehand. Chapter I (which those who are not beginners can skip) provides basic ideas on turbulence physics and is taken up in a textbook prepared by the teaching team of the ENPC (Benque, Viollet). Chapter II describes turbulent viscosity type modelling and the 2k-{epsilon} two equations model. It provides details of the channel flow case and the boundary conditions. Chapter III describes the `standard` (R{sub ij}-{epsilon}) Reynolds tensions transport model and introduces more recent models called `feasible`. A second paper deals with heat transfer and the effects of gravity, and returns to the Reynolds stress transport model. (author). 37 refs.

British Library Electronic Table of Contents (United Kingdom)

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

International Nuclear Information System (INIS)

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

British Library Electronic Table of Contents (United Kingdom)

This paper describes a comprehensive model of wastewater treatment in secondary facultative ponds, which combines 3D hydrodynamics with a mechanistic water quality model. The hydrodynamics are based on the Navier-Stokes equation for incompressible fluids under shallow water and Boussinesq assumptions capturing the flow dynamics along length, breadth and depth of the pond. The water quality sub model is based on the Activated Sludge Model (ASM) concept, describing COD and nutrient removal as function of bacterial growth following Monod kinetics, except for Escherichia coli removal, which was modelled as first order decay. The model was implemented in the Delft3D software and was used to evaluate the effect of wind and the addition of baffles on the water flow pattern, temperature profiles i...

Energy Technology Data Exchange (ETDEWEB)

The analytical results of blowdown characteristics and thrust forces were compared with the experiments, which were performed as pipe whip and jet discharge tests under the PWR LOCA conditions. The blowdown thrust forces were obtained by Navier-Stokes momentum equation for a single-phase, homogeneous and separated two-phase flow, assuming critical pressure at the exit if a critical flow condition was satisfied. The following results are obtained: (1) The node-junction method is useful for both the analyses of the blowdown thrust force and of the water hammer phenomena. (2) The Henry-Fauske model for subcooled critical flow is effective for the analysis of the maximum thrust force under the PWR LOCA conditions. The jet thrust parameter of the analysis and experiment is equal to 1.08. (3) The thrust parameter of saturated blowdown has the same one with the value under pressurized condition when the stagnant pressure is chosen as the saturated one. (4) The dominant ...

Energy Technology Data Exchange (ETDEWEB)

A spectral element method embodying Large Eddy Simulation based on Re- Normalization Group theory for simulating Sub Grid Scale viscosity was chosen for this work. This method is embodied in a computer code called NEKTON. NEKTON solves the unsteady, 2D or 3D,incompressible Navier Stokes equations by a spectral element method. The code was later extended to include the variable density and multiple reactive species effects at low Mach numbers, and to compute transport of large particles governed by inertia. Transport of small particles is computed by treating them as trace species. Code computations were performed for a number of test conditions typical of flow past a deep tube bank in a boiler. Results indicate qualitatively correct behavior. Predictions of deposition rates and deposit shape evolution also show correct qualitative behavior. These simulations are the first attempts to compute flow field results at realistic flow Reynolds numbers ...

Science.gov (United States)

Shape optimization of heat transfer augmentation device employed in turbine blade internal cooling passage has been performed numerically using single as well as multi-objective optimization procedures. Polynomial response surface approximation method and multi-objective genetic algorithm are used for single and multi-objective optimizations, respectively. Problem to enhance heat transfer rate considering staggered dimples on single surface of cooling passage has been formulated, and Reynolds-averaged Navier-Stokes equations are solved to analyze the flow field and the heat transfer. Three design variables defining channel and dimple dimensions, and two objective functions related to Nusselt number and friction drag are employed. Latin hypercube sampling is used to generate sampling points in design space, and the evaluated objectives are used to generate a set of optimal designs. Optimal shapes show higher heat transfer rates in the case of lower channel height, ...

Energy Technology Data Exchange (ETDEWEB)

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 simulation results inside an ESP. The real velocity profile was obtained from the velocity ...

Energy Technology Data Exchange (ETDEWEB)

Full text of publication follows: Hydrogen release and distribution in nuclear power plant containment is an important safety issue. Selection of a proper turbulence model is important for accurate estimation of the mixing process. The selection of turbulence model is dictated by the best compromise between accuracy and computational efforts. For this, three different turbulence models, viz. Standard k-{epsilon}, RNG k-{epsilon} and Reynolds Stress Model, based on Reynolds averaged Navier Stokes equations (RANS) approach, were used. The computations were done using the CFD code FLUENT, which is based on the control volume methodology. The computational results were compared with the experimental results of HYMIS test facility, where helium was used to simulate hydrogen. The processes of helium plume rise, multiple plume merging, distribution and mixing were studied. Based on these computations, a simple analytical/empirical zone based model was ...

CERN Document Server

We show the existence of global-in-time weak solutions to a general class of coupled FENE-type bead-spring chain models that arise from the kinetic theory of dilute solutions of polymeric liquids with noninteracting polymer chains. The class of models involves the unsteady incompressible Navier-Stokes equations in a bounded domain in two or three space dimensions for the velocity and the pressure of the fluid, with an elastic extra-stress tensor appearing on the right-hand side in the momentum equation. The extra-stress tensor stems from the random movement of the polymer chains and is defined by the Kramers expression through the associated probability density function that satisfies a Fokker-Planck-type parabolic equation, a crucial feature of which is the presence of a center-of-mass diffusion term. We require no structural assumptions on the drag term in the Fokker-Planck equation; in particular, the drag term need not be corotational. With a square-integrable ...

Science.gov (United States)

The current status of a computational aeroacoustics (CAA) approach to simulate broadband noise is reviewed. The method rests on the use of steady Reynolds averaged Navier-Stokes (RANS) simulation to describe the time-averaged motion of turbulent flow. By means of synthetic turbulence the steady one-point statistics (e.g. turbulence kinetic energy) and turbulent length- and time-scales of RANS are translated into fluctuations having statistics that very accurately reproduce the initial RANS target-setting. The synthetic fluctuations are used to prescribe sound sources which drive linear perturbation equations. The whole approach represents a methodology to solve statistical noise theory with state-of-the-art CAA tools in the time-domain. A brief overview of the synthetic turbulence model and its numerical discretization in terms of the random particle-mesh (RPM) and fast random particle-mesh (FRPM) method is given. Results are presented for trailing-edge noise, slat ...

International Nuclear Information System (INIS)

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 an internally heated volume in ...

Energy Technology Data Exchange (ETDEWEB)

The objective of this manuscript is to fully derive a geophysical multiphase model able to ''accommodate'' different multiphase turbulence approaches; viz., the Reynolds Averaged Navier-Stokes (RANS), the Large Eddy Simulation (LES), or hybrid RANSLES. This manuscript is the first part of a larger geophysical multiphase project--lead by LANL--that aims to develop comprehensive modeling tools for large-scale, atmospheric, transient-buoyancy dusty jets and plume (e.g., plinian clouds, nuclear ''mushrooms'', ''supercell'' forest fire plumes) and for boundary-dominated geophysical multiphase gravity currents (e.g., dusty surges, diluted pyroclastic flows, dusty gravity currents in street canyons). LES is a partially deterministic approach constructed on either a spatial- or a temporal-separation between the large and small scales of the flow, whereas RANS is an ...

International Nuclear Information System (INIS)

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

Energy Technology Data Exchange (ETDEWEB)

The Hyperion project was developed to determine an algorithm for assessing the risk of hydrate plug formation in the pipeline transport oil-water-gas mixtures at low temperatures. The project is a collaboration between physicists, chemists and engineers within the petroleum industry. This paper provided an overview of the project and outlined results obtained as the project entered its third and final year. The main objective of the project has been to understand the inherent mitigation effects of some oils on gas hydrate formation as well as to develop methods of predicting the risk of hydrate plugging. To date, the project has extracted and studied natural inhibiting components (NICs) in oils. Molecular modelling techniques have been used to study hydrate and fluid interfaces in order to estimate the driving force of agglomeration and growth through mechanical surface stress and measurement of surface wave fluctuations. A scheme is also being developed measure the adhesiveness of the ...

Energy Technology Data Exchange (ETDEWEB)

Fully coupled, Newton-Krylov algorithms are investigated for solving strongly coupled, nonlinear systems of partial differential equations arising in the field of computational fluid dynamics. Primitive variable forms of the steady incompressible and compressible Navier-Stokes and energy equations that describe the flow of a laminar Newtonian fluid in two-dimensions are specifically considered. Numerical solutions are obtained by first integrating over discrete finite volumes that compose the computational mesh. The resulting system of nonlinear algebraic equations are linearized using Newton`s method. Preconditioned Krylov subspace based iterative algorithms then solve these linear systems on each Newton iteration. Selected Krylov algorithms include the Arnoldi-based Generalized Minimal RESidual (GMRES) algorithm, and the Lanczos-based Conjugate Gradients Squared (CGS), Bi-CGSTAB, and Transpose-Free Quasi-Minimal Residual (TFQMR) algorithms. Both Incomplete ...

Science.gov (United States)

Under support from the French Space Agency (CNES), a 3D+t dusty-gas model of Comet 67P/Churyumov-Gerasimenko is being developed, to compute, from the first 2014 Rosetta orbital data, the aerodynamic forces exerted on the Rosetta orbiter and on the descent lander. We report the recently developed dust dynamics part of the code. The multi-species (presently H2O and CO) gas code is optimized in terms of computational speed owing to the use of two complementary methods: (a) 3D+t Direct Simulation Monte Carlo (DSMC) runs in the non-equilibrium regions adjacent to the surface and very distant from it, and (b) solutions of the Navier-Stokes equations in-between. The model is used presently using Lamy et al. (Space Sci. Rev., 2007, 128, 23) coarse information on 67P nucleus shape and rotation, and a range of possible gas production rates Q for the early Rosetta observations at rh 3 AU (Q 1026 - 1027 s-1). In the interim version, simplifying assumptions are made with ...

International Nuclear Information System (INIS)

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

Science.gov (United States)

The unsteady aerodynamic forces of a model fruit fly wing in flapping motion were investigated by numerically solving the Navier-Stokes equations. The flapping motion consisted of translation and rotation [the translation velocity (u(t)) varied according to the simple harmonic function (SHF), and the rotation was confined to a short period around stroke reversal]. First, it was shown that for a wing of given geometry with u(t) varying as the SHF, the aerodynamic force coefficients depended only on five non-dimensional parameters, i.e. Reynolds number (Re), stroke amplitude (Phi), mid-stroke angle of attack (alpha(m)), non-dimensional duration of wing rotation (Delta tau(r)) and rotation timing [the mean translation velocity at radius of the second moment of wing area (U), the mean chord length (c) and c/U were used as reference velocity, length and time, respectively]. Next, the force coefficients were investigated for a case in which typical values of these ...

Science.gov (United States)

Dippold, V., S. Mohler , Validation of the Wind-US Unstructured Flow Solver .... Hamed, A. and A. Mohamed, Assessment of Shock Induced Flow Separation and ...... Nozzles for Hypersonic Propulsion, NASA CR 185197, AGARD Symposium ...

DEFF Research Database (Denmark)

Level set method based segmentation provides an efficient tool for topological and geometrical shape handling. Conventional level set surfaces are only $C^0$ continuous since the level set evolution involves linear interpolation to compute derivatives. Bajaj et al. present a higher order method to evaluate level set surfaces that are $C^2$ continuous, but are slow due to high computational burden. In this paper, we provide a higher order GPU based solver for fast and efficient segmentation of large volumetric images. We also extend the higher order method to multi-domain segmentation. Our streaming solver is efficient in memory usage.

Energy Technology Data Exchange (ETDEWEB)

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.

CERN Document Server

An update of the ODEtools Maple package, for the analytical solving of 1st and 2nd order ODEs using Lie group symmetry methods, is presented. The set of routines includes an ODE-solver and user-level commands realizing most of the relevant steps of the symmetry scheme. The package also includes commands for testing the returned results, and for classifying 1st and 2nd order ODEs.

Energy Technology Data Exchange (ETDEWEB)

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

Energy Technology Data Exchange (ETDEWEB)

Production decline-curve analysis is accepted for prediction of future performance of oil and gas wells and fields. One type of decline-curve analysis involves the hyperbolic-decline-curve equation. The challenge with this equation is to determine simultaneously three parameters. The new method creates a computer spreadsheet containing a set of production data. Then, the specific arrangement of cells and columns containing data and/or formulas allows for statistical fits of the data. Finally, setting-up the nonlinear solver allows for optimization of the three parameters by maximization of the square of the regression coefficient. Eight field cases demonstrated the repeatable and rapid curve-fitting of the method.

Energy Technology Data Exchange (ETDEWEB)

We show that the eigenvalue equations describing a cylindrical ideal magnetophydrodynamicsw (MHD) plasma interacting with a thin resistive wall can be put into the standard mathematical form: ??? = ??? ?. This is accomplished by using a finite element basis for the plasma, and by adding an extra degree of freedom corresponding to the electrical current in the thin wall. The standard form allows the use of linear eigenvalue solvers, without additional interations, to compute the complete spectrum of plasma modes in the presence of a surrounding restrictive wall at arbitrary separation. We show that our method recovers standard results in the limits of (1) an infinitely resistive wall (no wall), and (2) a zero resistance wall (ideal wall).

British Library Electronic Table of Contents (United Kingdom)

This paper presents an innovative way in which university education can help pre-service teachers become better problem-solvers. The central idea is to use the ?Technology Fair?? as a means for promoting pre-service teachers pedagogical content knowledge about technological problem solving skills. This innovation is supported with results from a study carried out in autumn 2004. The purpose of the study was to investigate the influence of a procedure of working with primary school children to complete and present a technology fair project, on the educational value and meanings attached to problem solving skills by pre-service teachers. Pre-tests, mid-test and post-tests were administered to the pre-service teachers before, during, and after the preparation of the technology fair, respectiv...

Energy Technology Data Exchange (ETDEWEB)

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 model is also applied to industrial combustor geometries. (Author)

British Library Electronic Table of Contents (United Kingdom)

Ramjet burners are known to produce highly unsteady operating conditions with strong couplings between combustion, acoustics and flow dynamics. Predicting such operating limit-cycles still remains a difficult task for Computational Fluid Dynamics (CFD) although recent use of Large Eddy Simulation (LES) clearly opens new possibilities. The main difficulties for LES are to properly address numerically specific flow features at the same time. For example, a proper representation of the acoustic ramjet eigenmodes necessitates for the solver to be able to treat shocks often present at the inflow conditions without interfering with the low Mach number flow in the region of combustion. Chemistry modelling is another difficulty and it is still not clear what level of description is sufficient to r...

Energy Technology Data Exchange (ETDEWEB)

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

International Nuclear Information System (INIS)

To investigate the 2D/1D fusion core analysis method, a 1D neutron transport problem solver, PEACH-ID, is developed. It is a code of method of characteristics (MOC), both the usual fiat-source step characteristics (SC) scheme and linear source (LS) approximation scheme are adopted for tracking calculation along the neutron flying trajectory. Exponential function interpolation table and fission source extrapolation are adopted as two major methods to accelerate the computational process. Numerical results demonstrate that PEACH-1D is accurate and efficient, and the proposed LS scheme is able to handle quite larger mesh division and deserves much more application in the MOC codes. (authors)

Energy Technology Data Exchange (ETDEWEB)

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

International Nuclear Information System (INIS)

British Nuclear Fuels plc owns and operates the near-surface Drigg disposal facility for low level radioactive waste. The long-term performance of the site is modelled by a suite of computer codes called MONDRIAN. One of the modules of MONDRIAN deals with the transport of radionuclides through the environment, and this paper reports on the current status of this module (BIOS). We derive the basic set of working equations from first principles and show clearly how the approximate nature of the final equations is arrived at. This is done by an averaging process leading to compartments, in and out of which radionuclides, solids and water can flow. The equations allow radioactive decay chains and an arbitrary number of compartments. There is also the facility to deal with changes in the rate coefficients, thereby simulating different environmental states. It is also possible to include the creation of new compartments arising as a consequence of climatic variations. In addition to ...

Science.gov (United States)

A series of studies conducted to identify the factors that block and unblock problem solving is described. Through the construction of an isomorph of the classic "water jar" problems developed by A. S. Luchins (1942) as a dynamic graphic micro-world, several factors involved in producing blocked states were identified. Subjects included 10 individuals and 10 pairs of subjects. By comparing the behavior of individuals tackling the "missionaries and cannibals" problem to that of pairs of subjects solving this problem, a study identified means by which problem solvers operating in a social context are able to overcome blocks that discourage individuals. These studies point to the importance of "reflection" (evaluation of problem-solving results) for flexible problem solving. This research suggests that teaching students to analyze what they have done will help them develop flexibility in using a new approach when blocked. The results may also account for the ...

Energy Technology Data Exchange (ETDEWEB)

This paper describes the development of a computational multiphase fluid dynamics (CMFD) model of the Fischer Tropsch (FT) process in a Slurry Bubble Column Reactor (SBCR). The CMFD model is fundamentally based which allows it to be applied to different industrial processes and reactor geometries. The NPHASE CMFD solver [1] is used as the robust computational platform. Results from the CMFD model include gas distribution, species concentration profiles, and local temperatures within the SBCR. This type of model can provide valuable information for process design, operations and troubleshooting of FT plants. An ensemble-averaged, turbulent, multi-fluid solution algorithm for the multiphase, reacting flow with heat transfer was employed. Mechanistic models applicable to churn turbulent flow have been developed to provide a fundamentally based closure set for the equations. In this four-field model formulation, two of the fields are used to track the gas phase (i.e., ...

CERN Document Server

An important issue towards a broader acceptance of answer-set programming (ASP) is the deployment of tools which support the programmer during the coding phase. In particular, methods for debugging an answer-set program are recognised as a crucial step in this regard. Initial work on debugging in ASP mainly focused on propositional programs, yet practical debuggers need to handle programs with variables as well. In this paper, we discuss a debugging technique that is directly geared towards non-ground programs. Following previous work, we address the central debugging question why some interpretation is not an answer set. The explanations provided by our method are computed by means of a meta-programming technique, using a uniform encoding of a debugging request in terms of ASP itself. Our method also permits programs containing comparison predicates and integer arithmetics, thus covering a relevant language class commonly supported by all state-of-the-art ASP ...

Science.gov (United States)

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

Science.gov (United States)

Novice problem solvers are rather sensitive to surface problem features, and they often resort to trial and error formula matching rather than identifying an appropriate solution approach. These observations have been interpreted to imply that novices structure their knowledge according to surface features rather than according to problem type categories. However, it may also be the case that novices do know problem types, but cannot map the problem at hand to a known type, because they fail to create a sufficiently well-elaborated problem representation. This study aims to distinguish between these explanations. In this study novice physics students at high and low levels of proficiency completed two problem-sorting tasks from the domain of electricity and magnetism, one with and one without elaboration support. Results confirm that these students do distinguish problem types in accordance with their required solution approaches, and that their problem-sorting ...