scholarly journals Thermal modeling of hydrogen storage by absorption in a magnesium hydrides tank

2014 ◽  
Vol 5 ◽  
pp. A21 ◽  
Author(s):  
Karim Lahmer ◽  
Rachid Bessaïh ◽  
Angel Scipioni ◽  
Mohammed El Ganaoui
Keyword(s):  
Hydrogen Absorption ◽  
Numerical Scheme ◽  
Computational Time ◽  
Governing Equations ◽  
Time Step ◽  
Kinetic Reaction ◽  
Fully Implicit ◽  
Temporal Profiles ◽  
Tank Geometry ◽  
Reaction Equations

This paper summarizes numerical results of hydrogen absorption simulated in an axisymmetric tank geometry containing magnesium hydride heated to 300 °C and at moderate storage pressure 1 MPa. The governing equations are solved with a fully implicit finite volume numerical scheme used by a commercial software FLUENT. The effect of the different kinetic reaction equations modeling hydrogen absorption was studied by the introduction of a specific subroutine at each time step in order to consider which one will provide results close to available experimental results. Spatial and temporal profiles of temperature and concentration in hydride bed are plotted. Results show that suitable method for our two-dimensional study is a CV-2D technique because it generates the smallest error especially during the beginning of the reaction. Also, its computational time is the shortest one compared to the other methods.

Is Hydraulics Over-Used or Under-Used in Storm Drainage?

1994 ◽  
Vol 29 (1-2) ◽  
pp. 53-61
Author(s):  
Ben Chie Yen
Keyword(s):  
Unsteady Flow ◽  
Computational Time ◽  
Time Step ◽  
Flow Routing ◽  
Kinematic Wave Equation ◽  
Storm Drainage ◽  
Saint Venant Equations ◽  
Unsteady Flow Routing ◽  
Selection Of

Urban drainage models utilize hydraulics of different levels. Developing or selecting a model appropriate to a particular project is not an easy task. Not knowing the hydraulic principles and numerical techniques used in an existing model, users often misuse and abuse the model. Hydraulically, the use of the Saint-Venant equations is not always necessary. In many cases the kinematic wave equation is inadequate because of the backwater effect, whereas in designing sewers, often Manning's formula is adequate. The flow travel time provides a guide in selecting the computational time step At, which in turn, together with flow unsteadiness, helps in the selection of steady or unsteady flow routing. Often the noninertia model is the appropriate model for unsteady flow routing, whereas delivery curves are very useful for stepwise steady nonuniform flow routing and for determination of channel capacity.

scholarly journals A Fully Implicit Finite Volume Scheme for a Seawater Intrusion Problem in Coastal Aquifers

Water ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 1639
Author(s):  
Abdelkrim Aharmouch ◽  
Brahim Amaziane ◽  
Mustapha El Ossmani ◽  
Khadija Talali
Keyword(s):  
Finite Volume ◽  
Coastal Aquifers ◽  
Nonlinear Partial Differential Equations ◽  
Time Integration ◽  
Mass Conservation ◽  
Coupled System ◽  
Finite Volume Scheme ◽  
Time Step ◽  
Volume Scheme ◽  
Fully Implicit

We present a numerical framework for efficiently simulating seawater flow in coastal aquifers using a finite volume method. The mathematical model consists of coupled and nonlinear partial differential equations. Difficulties arise from the nonlinear structure of the system and the complexity of natural fields, which results in complex aquifer geometries and heterogeneity in the hydraulic parameters. When numerically solving such a model, due to the mentioned feature, attempts to explicitly perform the time integration result in an excessively restricted stability condition on time step. An implicit method, which calculates the flow dynamics at each time step, is needed to overcome the stability problem of the time integration and mass conservation. A fully implicit finite volume scheme is developed to discretize the coupled system that allows the use of much longer time steps than explicit schemes. We have developed and implemented this scheme in a new module in the context of the open source platform DuMu X . The accuracy and effectiveness of this new module are demonstrated through numerical investigation for simulating the displacement of the sharp interface between saltwater and freshwater in groundwater flow. Lastly, numerical results of a realistic test case are presented to prove the efficiency and the performance of the method.

On the Friction Coefficient in the Numerical Simulation of Tidal Wave Propagation

2010 ◽  
Author(s):  
Z. Y. Song ◽  
C. Cheng ◽  
F. M. Xu ◽  
J. Kong
Keyword(s):  
Numerical Simulation ◽  
Wave Propagation ◽  
Friction Factor ◽  
Tidal Wave ◽  
Computational Time ◽  
Numerical Dissipation ◽  
Courant Number ◽  
Time Step ◽  
Large Time Step ◽  
The Relationship

Based on the analytical solution of one-dimensional simplified equation of damping tidal wave and Heuristic stability analysis, the precision of numerical solution, computational time and the relationship between the numerical dissipation and the friction dissipation are discussed with different numerical schemes in this paper. The results show that (1) when Courant number is less than unity, the explicit solution of tidal wave propagation has higher precision and requires less computational time than the implicit one; (2) large time step is allowed in the implicit scheme in order to reduce the computational time, but the precision of the solution also reduce and the calculation precision should be guaranteed by reducing the friction factor: (3) the friction factor in the implicit solution is related to Courant number, presented as the determined friction factor is smaller than the natural value when Courant number is larger than unity, and their relationship formula is given from the theoretical analysis and the numerical experiments. These results have important application value for the numerical simulation of the tidal wave.

A General Rezoning Technique for KIVA3V Internal Combustion Engines CFD Simulations

2010 ◽  
Author(s):  
Randy P. Hessel ◽  
Ettore Musu ◽  
Salvador M. Aceves ◽  
Daniel L. Flowers
Keyword(s):  
Combustion Chamber ◽  
Internal Combustion Engines ◽  
Ad Hoc ◽  
National Laboratory ◽  
Internal Combustion ◽  
Computational Time ◽  
Combustion Engines ◽  
Combustion Chambers ◽  
Time Step ◽  
Engine Cycle

A computational mesh is required when performing CFD-combustion modeling of internal combustion engines. For combustion chambers with moving pistons and valves, like those in typical cars and trucks, the combustion chamber shape changes continually in response to piston and valve motion. The combustion chamber mesh must then also change at each time step to reflect that change in geometry. The method of changing the mesh from one computational time step to the next is called rezoning. This paper introduces a new method of mesh rezoning for the KIVA3V CFD-combustion program. The standard KIVA3V code from Los Alamos National Laboratory comes with standard rezoners that very nicely handle mesh motion for combustion chambers whose mesh does not include valves and for those with flat heads employing vertical valves. For pent-roof and wedge-roof designs KIVA3V offers three rezoners to choose from, the choice depending on how similar a combustion chamber is to the sample combustion chambers that come with KIVA3V. Often, the rezoners must be modified for meshes of new combustion chamber geometries to allow the mesh to successfully capture change in geometry during the full engine cycle without errors. There is no formal way to approach these modifications; typically this requires a long trial and error process to get a mesh to work for a full engine cycle. The benefit of the new rezoner is that it replaces the three existing rezoners for canted valve configurations with a single rezoner and has much greater stability, so the need for ad hoc modifications of the rezoner is greatly reduced. This paper explains how the new rezoner works and gives examples of its use.

scholarly journals Comparison of numerical schemes of river flood routing with an inertial approximation of the Saint Venant equations

RBRH ◽  
2018 ◽  
Vol 23 (0) ◽  
Author(s):  
Alice César Fassoni-Andrade ◽  
Fernando Mainardi Fan ◽  
Walter Collischonn ◽  
Artur César Fassoni ◽  
Rodrigo Cauduro Dias de Paiva
Keyword(s):  
Numerical Stability ◽  
Flood Routing ◽  
Computational Time ◽  
Numerical Schemes ◽  
Time Step ◽  
Simple Formulation ◽  
One Dimensional ◽  
Saint Venant Equations ◽  
Original Scheme ◽  
The One

ABSTRACT The one-dimensional flow routing inertial model, formulated as an explicit solution, has advantages over other explicit models used in hydrological models that simplify the Saint-Venant equations. The main advantage is a simple formulation with good results. However, the inertial model is restricted to a small time step to avoid numerical instability. This paper proposes six numerical schemes that modify the one-dimensional inertial model in order to increase the numerical stability of the solution. The proposed numerical schemes were compared to the original scheme in four situations of river’s slope (normal, low, high and very high) and in two situations where the river is subject to downstream effects (dam backwater and tides). The results are discussed in terms of stability, peak flow, processing time, volume conservation error and RMSE (Root Mean Square Error). In general, the schemes showed improvement relative to each type of application. In particular, the numerical scheme here called Prog Q(k+1)xQ(k+1) stood out presenting advantages with greater numerical stability in relation to the original scheme. However, this scheme was not successful in the tide simulation situation. In addition, it was observed that the inclusion of the hydraulic radius calculation without simplification in the numerical schemes improved the results without increasing the computational time.

scholarly journals MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma

2018 ◽  
Vol 615 ◽  
pp. A67 ◽  
Author(s):  
P. A. González-Morales ◽  
E. Khomenko ◽  
T. P. Downes ◽  
A. de Vicente
Keyword(s):  
Numerical Scheme ◽  
Operator Splitting ◽  
Conservation Equation ◽  
Ambipolar Diffusion ◽  
Point Of View ◽  
Integration Time ◽  
Solar Photosphere ◽  
Time Step ◽  
Diffusion Term ◽  
Fluid Approximation

The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere and low chromosphere. Under the single-fluid formalism the main non-ideal effects are described by a series of extra terms in the generalised induction equation and in the energy conservation equation. These effects are: Ohmic diffusion, ambipolar diffusion, the Hall effect, and the Biermann battery effect. From the point of view of the numerical solution of the single-fluid equations, when ambipolar diffusion or Hall effects dominate can introduce severe restrictions on the integration time step and can compromise the stability of the numerical scheme. In this paper we introduce two numerical schemes to overcome those limitations. The first of them is known as super time-stepping (STS) and it is designed to overcome the limitations imposed when the ambipolar diffusion term is dominant. The second scheme is called the Hall diffusion scheme (HDS) and it is used when the Hall term becomes dominant. These two numerical techniques can be used together by applying Strang operator splitting. This paper describes the implementation of the STS and HDS schemes in the single-fluid code MANCHA3D. The validation for each of these schemes is provided by comparing the analytical solution with the numerical one for a suite of numerical tests.

scholarly journals IL-GLOBO (1.0) – integrated Lagrangian particle model and Eulerian general circulation model GLOBO: development of the vertical diffusion module

2014 ◽  
Vol 7 (5) ◽  
pp. 2181-2191 ◽  
Author(s):  
D. Rossi ◽  
A. Maurizi
Keyword(s):  
General Circulation Model ◽  
General Circulation ◽  
Spline Interpolation ◽  
Circulation Model ◽  
Linear Interpolation ◽  
Computational Time ◽  
Integration Time ◽  
Particle Model ◽  
Vertical Diffusion ◽  
Time Step

Abstract. The development and validation of the vertical diffusion module of IL-GLOBO, a Lagrangian transport model coupled online with the Eulerian general circulation model GLOBO, is described. The module simulates the effects of turbulence on particle motion by means of a Lagrangian stochastic model (LSM) consistently with the turbulent diffusion equation used in GLOBO. The implemented LSM integrates particle trajectories, using the native σ-hybrid coordinates of the Eulerian component, and fulfils the well-mixed condition (WMC) in the general case of a variable density profile. The module is validated through a series of 1-D offline numerical experiments by assessing its accuracy in maintaining an initially well-mixed distribution in the vertical. A dynamical time-step selection algorithm with constraints related to the shape of the diffusion coefficient profile is developed and discussed. Finally, the skills of a linear interpolation and a modified Akima spline interpolation method are compared, showing that both satisfy the WMC with significant differences in computational time. A preliminary run of the fully integrated 3-D model confirms the result only for the Akima interpolation scheme while the linear interpolation does not satisfy the WMC with a reasonable choice of the minimum integration time step.

A New Unconditionally Stable Time Integration Method for Analysis of Nonlinear Structural Dynamics

2013 ◽  
Vol 80 (2) ◽  
Author(s):  
Ali Akbar Gholampour ◽  
Mehdi Ghassemieh ◽  
Mahdi Karimi-Rad
Keyword(s):  
Time Integration ◽  
Second Order ◽  
Numerical Dispersion ◽  
Computational Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Similar Time ◽  
Unconditionally Stable ◽  
Average Acceleration

A new time integration scheme is presented for solving the differential equation of motion with nonlinear stiffness. In this new implicit method, it is assumed that the acceleration varies quadratically within each time step. By increasing the order of acceleration, more terms of the Taylor series are used, which are expected to have responses with better accuracy than the classical methods. By considering this assumption and employing two parameters δ and α, a new family of unconditionally stable schemes is obtained. The order of accuracy, numerical dissipation, and numerical dispersion are used to measure the accuracy of the proposed method. Second order accuracy is achieved for all values of δ and α. The proposed method presents less dissipation at the lower modes in comparison with Newmark's average acceleration, Wilson-θ, and generalized-α methods. Moreover, this second order accurate method can control numerical damping in the higher modes. The numerical dispersion of the proposed method is compared with three unconditionally stable methods, namely, Newmark's average acceleration, Wilson-θ, and generalized-α methods. Furthermore, the overshooting effect of the proposed method is compared with these methods. By evaluating the computational time for analysis with similar time step duration, the proposed method is shown to be faster in comparison with the other methods.

A Modular Code for Real Time Dynamic Simulation of Gas Turbines in Simulink

2002 ◽  
Vol 128 (3) ◽  
pp. 506-517 ◽  
Author(s):  
S. M. Camporeale ◽  
B. Fortunato ◽  
M. Mastrovito
Keyword(s):  
Mathematical Model ◽  
Real Time ◽  
Gas Turbine ◽  
Dynamic Behavior ◽  
Gas Turbines ◽  
Simulation Software ◽  
Computational Time ◽  
High Fidelity ◽  
Time Step ◽  
The Mathematical Model

A high-fidelity real-time simulation code based on a lumped, nonlinear representation of gas turbine components is presented. The code is a general-purpose simulation software environment useful for setting up and testing control equipments. The mathematical model and the numerical procedure are specially developed in order to efficiently solve the set of algebraic and ordinary differential equations that describe the dynamic behavior of gas turbine engines. For high-fidelity purposes, the mathematical model takes into account the actual composition of the working gases and the variation of the specific heats with the temperature, including a stage-by-stage model of the air-cooled expansion. The paper presents the model and the adopted solver procedure. The code, developed in Matlab-Simulink using an object-oriented approach, is flexible and can be easily adapted to any kind of plant configuration. Simulation tests of the transients after load rejection have been carried out for a single-shaft heavy-duty gas turbine and a double-shaft aero-derivative industrial engine. Time plots of the main variables that describe the gas turbine dynamic behavior are shown and the results regarding the computational time per time step are discussed.

Sign in / Sign up

Export Citation Format

Share Document