Construction of reachable sets of controlled systems with second order of accuracy with respect to time step

2013 ◽  
Vol 23 (02) ◽  
pp. 273-305 ◽  
Author(s):  
JIANYAO YAO ◽  
G. R. LIU ◽  
DONG QIAN ◽  
CHUNG-LUNG CHEN ◽  
GEORGE X. XU

A computational fluid dynamics (CFD) solver based on the gradient smoothing method (GSM) with moving mesh enabled is presented in this paper. The GSM uses unstructured meshes which could be generated and remeshed easily. The spatial derivatives of field variables at nodes and midpoints of cell edges are calculated using the gradient smoothing operations. The presented GSM codes use second-order Roes upwind flux difference splitting method and second-order 3-level backward differencing scheme for the compressible Navier–Stokes equations with moving mesh, and the second-order of accuracy for both the spatial and temporal discretization is ensured. The spatial discretization accuracy is verified using the method of manufactured solutions (MMS) on both structured and unstructured triangle meshes, and the results show that the observed order of accuracy achieves 2 even when highly distorted meshes are used. The temporal discretization accuracy is verified using the results with different time step lengths, and second-order accuracy is also obtained. Therefore, it is confirmed that the proposed GSM-CFD solver is a uniform second-order scheme.


2020 ◽  
Vol 20 (4) ◽  
pp. 769-782
Author(s):  
Amiya K. Pani ◽  
Vidar Thomée ◽  
A. S. Vasudeva Murthy

AbstractWe analyze a second-order in space, first-order in time accurate finite difference method for a spatially periodic convection-diffusion problem. This method is a time stepping method based on the first-order Lie splitting of the spatially semidiscrete solution. In each time step, on an interval of length k, of this solution, the method uses the backward Euler method for the diffusion part, and then applies a stabilized explicit forward Euler approximation on {m\geq 1} intervals of length {\frac{k}{m}} for the convection part. With h the mesh width in space, this results in an error bound of the form {C_{0}h^{2}+C_{m}k} for appropriately smooth solutions, where {C_{m}\leq C^{\prime}+\frac{C^{\prime\prime}}{m}}. This work complements the earlier study [V. Thomée and A. S. Vasudeva Murthy, An explicit-implicit splitting method for a convection-diffusion problem, Comput. Methods Appl. Math. 19 2019, 2, 283–293] based on the second-order Strang splitting.


2011 ◽  
Author(s):  
Allaberen Ashyralyev ◽  
Serhat Yilmaz ◽  
Theodore E. Simos ◽  
George Psihoyios ◽  
Ch. Tsitouras ◽  
...  

2019 ◽  
Vol 27 (4) ◽  
pp. 457-468 ◽  
Author(s):  
Allaberen Ashyralyev ◽  
Abdullah Said Erdogan ◽  
Ali Ugur Sazaklioglu

Abstract The present paper is devoted to the investigation of a source identification problem that describes the flow in capillaries in the case when an unknown pressure acts on the system. First and second order of accuracy difference schemes are presented for the numerical solution of this problem. Almost coercive stability estimates for these difference schemes are established. Additionally, some numerical results are provided by testing the proposed methods on an example.


2013 ◽  
Vol 80 (2) ◽  
Author(s):  
Ali Akbar Gholampour ◽  
Mehdi Ghassemieh ◽  
Mahdi Karimi-Rad

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.


2011 ◽  
Vol 1 (3) ◽  
pp. 215-234 ◽  
Author(s):  
Lizhen Chen ◽  
Jie Shen ◽  
Chuanju Xu

AbstractWe propose and analyze spectral direction splitting schemes for the incompressible Navier-Stokes equations. The schemes combine a Legendre-spectral method for the spatial discretization and a pressure-stabilization/direction splitting scheme for the temporal discretization, leading to a sequence of one-dimensional elliptic equations at each time step while preserving the same order of accuracy as the usual pressure-stabilization schemes. We prove that these schemes are unconditionally stable, and present numerical results which demonstrate the stability, accuracy, and efficiency of the proposed methods.


Author(s):  
A.M. Lipanov ◽  

In this work, an implicit method is proposed to numerically solve a system of the onedimensional nonstationary equations of gas dynamics transformed by the method of characteristics. Internal points of the channel for a solid-propellant charge are considered at a preignition period of the solid-propellant rocket engine operation. The use of the implicit method makes it possible to calculate the values of gas-dynamic parameters at nodal points of the regular coordinate grid. Calculations of the gas-dynamic parameters both when integrating over time and along the spatial coordinate are performed with the second order of accuracy. Both subsonic and supersonic flows are studied. It is shown that, when predicting the expected pressure value during the transition from one time layer to another with the second order of accuracy, the twenty-fold efficiency of the implicit method is achieved in comparison with the explicit difference method. The trial calculation is performed.


Author(s):  
Luca Casarsa ◽  
Pietro Giannattasio ◽  
Diego Micheli

A simple and efficient numerical model is presented for the simulation of pulse combustors. It is based on the numerical solution of the quasi-1D unsteady flow equations and on phenomenological sub-models of turbulence and combustion. The gas dynamics equations are solved by using the Flux Difference Splitting (FDS) technique, a finite-volume upwind numerical scheme, and ENO reconstructions to obtain second-order accurate non-oscillatory solutions. The numerical fluxes computed at the cell interfaces are used to transport also the reacting species, their formation energy and the turbulent kinetic energy. The combustion progress in each cell is evaluated explicitly at the end of each time step according to a second-order overall reaction kinetics. In this way, the computations of gas dynamic evolution and heat release are decoupled, which makes the model particularly simple and efficient. A comprehensive set of measurements has been performed on a small Helmholtz type pulse-jet in order to validate the model. Air and fuel consumptions, wall temperatures, pressure cycles in both combustion chamber and tail-pipe, and instantaneous thrust have been recorded in different operating conditions of the device. The comparison between numerical and experimental results turns out to be satisfactory in all the working conditions of the pulse-jet. In particular, accurate predictions are obtained of the device operating frequency and of shape, amplitude and phase of the pressure waves in both combustion chamber and tail-pipe.


Sign in / Sign up

Export Citation Format

Share Document