scholarly journals Efficient implementation of advanced Richardson Extrapolation in an atmospheric chemical scheme

Author(s):  
Zahari Zlatev ◽  
Ivan Dimov ◽  
István Faragó ◽  
Krassimir Georgiev ◽  
Ágnes Havasi

AbstractThe numerical treatment of an atmospheric chemical scheme, which contains 56 species, is discussed in this paper. This scheme is often used in studies of air pollution levels in different domains, as, for example, in Europe, by large-scale environmental models containing additionally two other important physical processes—transport of pollutants in the atmosphere (advection) and diffusion phenomena. We shall concentrate our attention on the efficient numerical treatment of the chemical scheme by using Implicit Runge–Kutta Methods combined with accurate and efficient advanced versions of the Richardson Extrapolation. A Variable Stepsize Variable Formula Method is developed in order to achieve high accuracy of the calculated results within a reasonable computational time. Reliable estimations of the computational errors when the proposed numerical methods are used in the treatment of the chemical scheme will be demonstrated by presenting results from several representative runs and comparing these results with “exact” concentrations obtained by applying a very small stepsize during the computations. Results related to the diurnal variations of some of the chemical species will also be presented. The approach used in this paper does not depend on the particular chemical scheme and can easily be applied when other atmospheric chemical schemes are selected.

2019 ◽  
Vol 01 (02) ◽  
pp. 1950009 ◽  
Author(s):  
ZAHARI ZLATEV ◽  
IVAN DIMOV ◽  
ISTVÁN FARAGÓ ◽  
KRASSIMIR GEORGIEV ◽  
ÁGNES HAVASI

The treatment of large-scale air pollution models is not only important for modern society, but also an extremely difficult task. Five important physical and chemical processes: (1) horizontal advection, (2) horizontal diffusion, (3) vertical exchange, (4) emission of different pollutants and (5) dry and wet deposition have to be united and handled together. This leads to huge computational problems which have to be treated on modern high-speed parallel architectures by using advanced numerical methods. The computational difficulties are vastly increased when such models are used in the investigation of the impact of climatic changes on high pollution levels that are often exceeding certain critical concentrations and, therefore, are becoming harmful for plants, animals and human beings. There are two major reasons for the great increase in computational difficulties: (a) it is necessary to run the discretized data on fine spatial grid models over very long time-periods consisting of many consecutive years and (b) different scenarios are to be used in order to investigate the sensitivity of the pollution levels to systematic variations of several carefully selected key parameters. How the major difficulties can be resolved is explained here. Furthermore, many results are presented in order to demonstrate the ability of the model to handle successfully the huge computational tasks by using fine space and time discretization and by applying a chemical scheme with many chemical species participating in several hundred chemical reactions. Our results indicate that the climatic changes will often lead to some increases in the pollution levels.


2011 ◽  
Vol 2011 ◽  
pp. 1-16
Author(s):  
Zahari Zlatev ◽  
Ivan Dimov ◽  
István Faragó ◽  
Krassimir Georgiev ◽  
Ágnes Havasi ◽  
...  

Advection equations appear often in large-scale mathematical models arising in many fields of science and engineering. The Crank-Nicolson scheme can successfully be used in the numerical treatment of such equations. The accuracy of the numerical solution can sometimes be increased substantially by applying the Richardson Extrapolation. Two theorems related to the accuracy of the calculations will be formulated and proved in this paper. The usefulness of the combination consisting of the Crank-Nicolson scheme and the Richardson Extrapolation will be illustrated by numerical examples.


2019 ◽  
Author(s):  
Liqun Cao ◽  
Jinzhe Zeng ◽  
Mingyuan Xu ◽  
Chih-Hao Chin ◽  
Tong Zhu ◽  
...  

Combustion is a kind of important reaction that affects people's daily lives and the development of aerospace. Exploring the reaction mechanism contributes to the understanding of combustion and the more efficient use of fuels. Ab initio quantum mechanical (QM) calculation is precise but limited by its computational time for large-scale systems. In order to carry out reactive molecular dynamics (MD) simulation for combustion accurately and quickly, we develop the MFCC-combustion method in this study, which calculates the interaction between atoms using QM method at the level of MN15/6-31G(d). Each molecule in systems is treated as a fragment, and when the distance between any two atoms in different molecules is greater than 3.5 Å, a new fragment involved two molecules is produced in order to consider the two-body interaction. The deviations of MFCC-combustion from full system calculations are within a few kcal/mol, and the result clearly shows that the calculated energies of the different systems using MFCC-combustion are close to converging after the distance thresholds are larger than 3.5 Å for the two-body QM interactions. The methane combustion was studied with the MFCC-combustion method to explore the combustion mechanism of the methane-oxygen system.


2021 ◽  
Vol 12 (3) ◽  
pp. 212-231
Author(s):  
Issam El Hammouti ◽  
Azza Lajjam ◽  
Mohamed El Merouani

The berth allocation problem is one of the main concerns of port operators at a container terminal. In this paper, the authors study the berth allocation problem at the strategic level commonly known as the strategic berth template problem (SBTP). This problem aims to find the best berth template for a set of calling ships accepted to be served at the port. At strategic level, port operator can reject some ships to be served for avoid congestion. Since the computational complexity of the mathematical formulation proposed for SBTP, solution approaches presented so far for the problem are limited especially at level of large-scale instances. In order to find high quality solutions with a short computational time, this work proposes a population based memetic algorithm which combine a first-come-first-served (FCFS) technique, two genetics operators, and a simulating annealing algorithm. Different computational experiences and comparisons against the best known solutions so far have been presented to show the performance and effectiveness of the proposed method.


2016 ◽  
Author(s):  
Janek Meyer ◽  
Hannes Renzsch ◽  
Kai Graf ◽  
Thomas Slawig

While plain vanilla OpenFOAM has strong capabilities with regards to quite a few typical CFD-tasks, some problems actually require additional bespoke solvers and numerics for efficient computation of high-quality results. One of the fields requiring these additions is the computation of large-scale free-surface flows as found e.g. in naval architecture. This holds especially for the flow around typical modern yacht hulls, often planing, sometimes with surface-piercing appendages. Particular challenges include, but are not limited to, breaking waves, sharpness of interface, numerical ventilation (aka streaking) and a wide range of flow phenomenon scales. A new OF-based application including newly implemented discretization schemes, gradient computation and rigid body motion computation is described. In the following the new code will be validated against published experimental data; the effect on accuracy, computational time and solver stability will be shown by comparison to standard OF-solvers (interFoam / interDyMFoam) and Star CCM+. The code’s capabilities to simulate complex “real-world” flows are shown on a well-known racing yacht design.


2003 ◽  
Vol 125 (4) ◽  
pp. 234-241 ◽  
Author(s):  
Vincent Y. Blouin ◽  
Michael M. Bernitsas ◽  
Denby Morrison

In structural redesign (inverse design), selection of the number and type of performance constraints is a major challenge. This issue is directly related to the computational effort and, most importantly, to the success of the optimization solver in finding a solution. These issues are the focus of this paper, which provides and discusses techniques that can help designers formulate a well-posed integrated complex redesign problem. LargE Admissible Perturbations (LEAP) is a general methodology, which solves redesign problems of complex structures with, among others, free vibration, static deformation, and forced response amplitude constraints. The existing algorithm, referred to as the Incremental Method is improved in this paper for problems with static and forced response amplitude constraints. This new algorithm, referred to as the Direct Method, offers comparable level of accuracy for less computational time and provides robustness in solving large-scale redesign problems in the presence of damping, nonstructural mass, and fluid-structure interaction effects. Common redesign problems include several natural frequency constraints and forced response amplitude constraints at various frequencies of excitation. Several locations on the structure and degrees of freedom can be constrained simultaneously. The designer must exercise judgment and physical intuition to limit the number of constraints and consequently the computational time. Strategies and guidelines are discussed. Such techniques are presented and applied to a 2,694 degree of freedom offshore tower.


2010 ◽  
Vol 10 (14) ◽  
pp. 6645-6660 ◽  
Author(s):  
P. Huszar ◽  
D. Cariolle ◽  
R. Paoli ◽  
T. Halenka ◽  
M. Belda ◽  
...  

Abstract. In general, regional and global chemistry transport models apply instantaneous mixing of emissions into the model's finest resolved scale. In case of a concentrated source, this could result in erroneous calculation of the evolution of both primary and secondary chemical species. Several studies discussed this issue in connection with emissions from ships and aircraft. In this study, we present an approach to deal with the non-linear effects during dispersion of NOx emissions from ships. It represents an adaptation of the original approach developed for aircraft NOx emissions, which uses an exhaust tracer to trace the amount of the emitted species in the plume and applies an effective reaction rate for the ozone production/destruction during the plume's dilution into the background air. In accordance with previous studies examining the impact of international shipping on the composition of the troposphere, we found that the contribution of ship induced surface NOx to the total reaches 90% over remote ocean and makes 10–30% near coastal regions. Due to ship emissions, surface ozone increases by up to 4–6 ppbv making 10% contribution to the surface ozone budget. When applying the ship plume parameterization, we show that the large scale NOx decreases and the ship NOx contribution is reduced by up to 20–25%. A similar decrease was found in the case of O3. The plume parameterization suppressed the ship induced ozone production by 15–30% over large areas of the studied region. To evaluate the presented parameterization, nitrogen monoxide measurements over the English Channel were compared with modeled values and it was found that after activating the parameterization the model accuracy increases.


Author(s):  
Feng Jie Zheng ◽  
Fu Zheng Qu ◽  
Xue Guan Song

Reservoir-pipe-valve (RPV) systems are widely used in many industrial process. The pressure in an RPV system plays an important role in the safe operation of the system, especially during the sudden operation such as rapid valve opening/closing. To investigate the pressure especially the pressure fluctuation in an RPV system, a multidimensional and multiscale model combining the method of characteristics (MOC) and computational fluid dynamics (CFD) method is proposed. In the model, the reservoir is modeled by a zero-dimensional virtual point, the pipe is modeled by a one-dimensional MOC, and the valve is modeled by a three-dimensional CFD model. An interface model is used to connect the multidimensional and multiscale model. Based on the model, a transient simulation of the turbulent flow in an RPV system is conducted, in which not only the pressure fluctuation in the pipe but also the detailed pressure distribution in the valve are obtained. The results show that the proposed model is in good agreement with the full CFD model in both large-scale and small-scale spaces. Moreover, the proposed model is more computationally efficient than the CFD model, which provides a feasibility in the analysis of complex RPV system within an affordable computational time.


2021 ◽  
Author(s):  
Itsuki Sugita ◽  
Shohei Matsuyama ◽  
Hiroki Dobashi ◽  
Daisuke Komura ◽  
Shumpei Ishikawa

Here, we present Viola, a Python package that provides structural variant (SV; large scale genome DNA variations that can result in disease, e.g., cancer) signature analytical functions and utilities for custom SV classification, merging multi-SV-caller output files, and SV annotation. We demonstrate that Viola can extract biologically meaningful SV signatures from publicly available SV data for cancer and we evaluate the computational time necessary for annotation of the data.


Author(s):  
Juan Gea Bermúdez ◽  
Kaushik Das ◽  
Hardi Koduvere ◽  
Matti Juhani Koivisto

This paper proposes a mathematical model to simulate Day-ahead markets of large-scale multi-energy systems with high share of renewable energy. Furthermore, it analyses the importance of including unit commitment when performing such analysis. The results of the case study, which is performed for the North Sea region, show the influence of massive renewable penetration in the energy sector and increasing electrification of the district heating sector towards 2050, and how this impacts the role of other energy sources such as thermal and hydro. The penetration of wind and solar is likely to challenge the need for balancing in the system as well as the profitability of thermal units. The degree of influence of the unit commitment approach is found to be dependent on the configuration of the energy system. Overall, including unit commitment constraints with integer variables leads to more realistic behaviour of the units, at the cost of increasing considerably the computational time. Relaxing integer variables reduces significantly the computational time, without highly compromising the accuracy of the results. The proposed model, together with the insights from the study case, can be specially useful for system operators for optimal operational planning.


Sign in / Sign up

Export Citation Format

Share Document