Most reservoir simulations use vertical flow performance (VFP) tables to represent flow in the tubing, which ignores the flow dynamics in wells and their downstream gathering and transportation networks. On the other hand, most dynamic well-pipeline flow models use pressure-rate equations to describe the inflow/injection from/to the reservoir, which ignores the flow and pressure transients in the near-wellbore regions. Obviously, neither of the two types of modeling can account for the transient reservoir-well/pipeline flow interactions that can be of great importance in many operation scenarios. To bridge this modeling gap, a joint industry program (JIP) was initiated based on a previous successful investigation [1] on the feasibility to implicitly couple a reservoir model to a dynamic well-pipeline multiphase flow model. The aim of the JIP was to deliver a simulation tool that can fulfill the industry's basic requirements on the modeling of transient flow interactions across the sandface, which should also lay the foundation for its future expansion in functionalities. This paper is intended to summarize the outcomes of the JIP. First of all, the paper discusses the need for integrated dynamic modeling and reviews earlier efforts on building integrated dynamic reservoir-well/pipeline systems. Secondly the paper describes the details of building the reservoir model and how to couple it to a well-pipeline flow model to assure numerical stability and simulation speed for cases of interest. Thirdly, the paper shows the PVT and fluid handling options the integrated simulation tool can provide. These are identical for the two models in order to keep the consistence of fluid properties particularly when the same fluid flows back-and-forth between the wellbore and the near-wellbore during transients. Two application cases based on the resulting model are presented in the paper. One case is about simulation of chemical squeezement hydraulics for refining the operational procedure, and the other case is about quantification of the pressure gradient in the near-wellbore formation for different well bean-up procedures in order to prevent sand production. The two cases demonstrate the advantages of using the coupled wellbore-reservoir transient modeling.

The paper describes an equation solver superstructure developed for a sequential modular dynamic process simulation system as part of a Eureka project with Norwegian and British participation. The purpose of the development was combining some of the advantages of equation based and purely sequential systems, enabling implicit treatment of key variables independent of module boundaries, and use of numerical integration techniques suitable for each individual type of variable. For training simulator applications the main advantages are gains in speed due to increased stability limits on time steps and improved consistency of simulation results. The system is split into an off-line analysis phase and an on-line equation solver. The off-line processing consists of automatic determination of the topological structure of the system connectivity from standard process description files and derivation of an optimized sparse matrix solution procedure for the resulting set of equations. The on-line routine collects equation coefficients from involved modules, solves the combined sets of structured equations, and stores the results appropriately. This method minimizes the processing cost during the actual simulation. The solver has been applied in the Veslefrikk training simulator project

Simulation of a planar Solid Oxide Fuel Cell stack requires the solution of the mass balances of the chemical species, the energy balances, the charge balance and the channel flow equations in order to compute the species concentrations, the temperature distributions, the current density and the channel flows. The unit cell geometry can be taken into account by combining detailed modeling of a unit cell with a homogenized model of a whole stack. In this study the effect of the asymmetric temperature distribution on the channel flows in a conventional cross-flow design has been investigated. The bidirectional cross-flow design is introduced, for which we can show more directional temperature and flow distributions