Meshing Method for Transient Computation With Complex Geometry in CFD

2002 ◽  
Author(s):  
J. L. Kueny ◽  
B. Dore´ ◽  
G. Coppens
Keyword(s):  
Complex Geometry ◽  
Mesh Adaptation ◽  
Moving Mesh ◽  
Computational Time ◽  
Tangential Displacement ◽  
Time Step ◽  
Local Observation ◽  
Complex Boundary ◽  
Accuracy Of Results ◽  
Main Ideas

The description of mesh evolution during a transient computation with moving walls and mesh adaptation has to respect many rules. Using a good mesh for each computational time step is important for accuracy of results. The complexity of geometry can make this objective more complex. A method has been developed to obtain a good moving mesh description with complex boundary geometry. It is based on a local observation of boundary movement and can be resumed by two main ideas: • Add cells where the volume of solution domain increases. • Slide the mesh where the boundary has a tangential displacement.

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.

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

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.

scholarly journals A Sequential Approach to Numerical Simulations of Solidification with Domain and Time Decomposition

2019 ◽  
Vol 9 (10) ◽  
pp. 1972 ◽  
Author(s):  
Elzbieta Gawronska
Keyword(s):  
Parallel Computing ◽  
Computational Cost ◽  
Fixed Time ◽  
Simulation Software ◽  
Computational Method ◽  
Computational Time ◽  
Time Step ◽  
Sequential Approach ◽  
Time Partitioning ◽  
Physical Phenomena

Progress in computational methods has been stimulated by the widespread availability of cheap computational power leading to the improved precision and efficiency of simulation software. Simulation tools become indispensable tools for engineers who are interested in attacking increasingly larger problems or are interested in searching larger phase space of process and system variables to find the optimal design. In this paper, we show and introduce a new approach to a computational method that involves mixed time stepping scheme and allows to decrease computational cost. Implementation of our algorithm does not require a parallel computing environment. Our strategy splits domains of a dynamically changing physical phenomena and allows to adjust the numerical model to various sub-domains. We are the first (to our best knowledge) to show that it is possible to use a mixed time partitioning method with various combination of schemes during binary alloys solidification. In particular, we use a fixed time step in one domain, and look for much larger time steps in other domains, while maintaining high accuracy. Our method is independent of a number of domains considered, comparing to traditional methods where only two domains were considered. Mixed time partitioning methods are of high importance here, because of natural separation of domain types. Typically all important physical phenomena occur in the casting and are of high computational cost, while in the mold domains less dynamic processes are observed and consequently larger time step can be chosen. Finally, we performed series of numerical experiments and demonstrate that our approach allows reducing computational time by more than three times without losing the significant precision of results and without parallel computing.

scholarly journals Improvement of anisotropic porosity models with a merging technique

2018 ◽  
Vol 40 ◽  
pp. 06023
Author(s):  
Martin Bruwier ◽  
Pierre Archambeau ◽  
Sébastien Erpicum ◽  
Michel Pirotton ◽  
Benjamin Dewals
Keyword(s):  
Shallow Water ◽  
Flow Characteristics ◽  
Threshold Value ◽  
Computational Grid ◽  
Computational Time ◽  
Time Step ◽  
Shallow Water Models ◽  
The Cost ◽  
The Impact ◽  
Porosity Model

Anisotropic porosity shallow-water models are used to take into account detailed topographic information through porosity parameters multiplying the various terms of the shallow-water equations. A storage porosity is assigned to each cell to reflect the void fraction in the cell and a conveyance porosity is used at each edge to reproduce the impact of subgrid obstacles on the flux terms. To guaranty the numerical stability, the time step depends on the value of the porosity parameters. This may hamper severely the computational efficiency in the presence of cells with low values of storage porosity. Cartesian grids are particularly sensitive to such a case since the meshing stems directly from the choice of the grid size. In this paper, this problem is addressed by using an original merging technique consisting in merging cells with a storage porosity lower than a threshold value with neighbouring cells. The model was tested for modelling a prismatic channel with different orientations between the Cartesian computational grid and the channel direction. The results show that the standard anisotropic porosity model (without merging) improves the reproduction of the flow characteristics; but at the cost of a significantly higher computational time. In contrast, the computational time is drastically reduced and the accuracy preserved when the merging technique is used with the porosity model.

Sign in / Sign up

Export Citation Format

Share Document