Search results

The purpose of this study is to investigate the errors arising from the numerical treatment of model processes, paying particular attention to the impact of key system features including widely variable dispersion coefficients, spatiotemporal velocities of algal cells, and the aggregation of algae from single cells to large colonies. An advection–dispersion model has been presented to describe the vertical transport of colonial and motile harmful algae in a lake environment.

Model performance is examined for two different numerical treatments of the advective term: first-order upwind and quadratic upwind with a stability-preserving flux limiter (SMART). To determine how these schemes impact predictions, comparisons are made across a sequence of models with increasing complexity.

Using first-order upwinding for advection–dispersion calculations with a time oscillating velocity field leads to oscillatory numerical dispersion. Subjecting an initially uniform distribution of large-sized algal colonies to a spatiotemporal velocity creates a concentration pulse, which reaches a steady-state width at high-grid Peclet numbers when using the SMART scheme; the pulse exhibits contraction–expansion behavior throughout a velocity cycle at all Peclet numbers when using first-order upwinding. When aggregation dynamics are included with advection-dominated spatiotemporal transport, results indicate the SMART scheme predicts larger peak concentration values than those predicted by first-order upwind, but peak location and the time to large colony appearance remain largely unchanged between the two advective schemes.

To the best of the authors’ knowledge, this study is the first numerical investigation of a novel advection–dispersion model of vertical algal transport. In addition, a generalized expression for the effective dispersion coefficient of temporally variable flow fields is presented.

The purpose of this paper is to investigate the ways to diminish or eliminate numerical diffusion and dispersion. Numerical dispersion and diffusion are present in the predicted macrosegregation profiles reported in the literature and they hinder the interpretation of the simulation results. With the motivation to eliminate these numerical problems by employing appropriate meshes, simulations of macrosegregation in a billet direct‐chill cast from a multi‐component aluminium alloy has been performed.

First the idea that numerical dispersion could be alleviated by refining the structured mesh size is tested and the extent of this mesh refining to overcome these numerical problems is discussed. Second the link of numerical dispersion and diffusion to the type of mesh used is investigated.

Unstructured mesh eliminates the numerical dispersion present in the structured mesh while it introduces the numerical diffusion. It is concluded by performing calculations with the same settings but different meshes that, although refining the structured mesh could alleviate the numerical oscillation, it increases the computation time dramatically. Therefore the best solution to overcome these numerical problems is the employment of a hybrid mesh consisting of both structured and unstructured mesh.

This work reveals the reasons behind the numerical dispersion and diffusion in macrosegregation modelling and gives a practical solution.

The classical advection-dispersion equation (ADE) model cannot accurately depict the gas transport process in natural geological formations. This paper aims to study the behavior of CO₂ transport in fractal porous media by using an effective Hausdorff fractal derivative advection-dispersion equation (HFDADE) model.

Anomalous dispersion behaviors of CO₂ transport are effectively characterized by the investigation of time and space Hausdorff derivatives on non-Euclidean fractal metrics. The numerical simulation has been performed with different Hausdorff fractal dimensions to reveal characteristics of the developed fractal ADE in fractal porous media. Numerical experiments focus on the influence of the time and space fractal dimensions on flow velocity and dispersion coefficient.

The physical mechanisms of parameters in the Hausdorff fractal derivative model are analyzed clearly. Numerical results demonstrate that the proposed model can well fit the history of gas production data and it can be a powerful technique for depicting the early arrival and long-tailed phenomenon by incorporating a fractal dimension.

To the best of the authors’ knowledge, first time these results are presented.

This paper aims to tackle the problem of thermo‐solutal convection and macrosegregation during ingot solidification of metal alloys. Complex flow structures associated with the development of channels segregate and sharp gradients in the solutal field call for the implementation of accurate methods for numerical modeling of alloy solidification. In particular, the solute transport equation is convection dominated and requires special non‐oscillarity type high‐order schemes to handle the regions of channels segregates.

In the present study, a time‐splitting approach has been adopted to separately handle solute advection and diffusion. This splitting technique allows the application of accurate total variation dimensioning (TVD) schemes for solution of solute advection. Applications of second‐order Lax‐Wendroff TVD SUPERBEE and fifth‐order weighted essentially non‐oscillatory (WENO) schemes are described in the present article. Classical numerical solution of solute transport using hybrid and central‐difference schemes are also employed for the purpose of comparisons. Numerical simulations for solidification of Pb‐18%Sn in a two‐dimensional rectangular cavity have been carried out using different numerical schemes.

Numerical results show the difficulty of obtaining grid‐independent solutions with respect to local details in the region of channels. Grid convergence patterns and numerical uncertainty are found to be dependent on the applied scheme. In general, the first‐order hybrid scheme is diffusive and under predicts the formation of channels. The second‐order central‐difference scheme brings about oscillations with possible non‐physical extremes of solute composition in the region of channel segregates due to sharp gradients in the solutal field. The results obtained using TVD and WENO schemes contain no oscillations and show an excellent capture of channels formation and resolution of the interface between solute‐rich and depleted bands. Different stages of channels formation are followed by analyzing thermo‐solutal convection and macrosegregation at different times during solidification.

Accurate prediction of local variation in the solutal and flow fields in the channels regions requires grid refinement up to scales in the order of microscopic dendrite arm spacing. This imposes limitations in terms of large computational time and applicability of available macroscopic models based on classical volume‐averaging techniques.

The present study is very useful for numerical simulation of macrosegregation during ingot casting of metal alloys.

The paper provides the methodology and application of TVD schemes to predict channel segregates during columnar solidification of metal alloys. It also demonstrates the limitations of classical schemes for simulation of alloy solidification.

A computer model based on the fractional step method is presented for modelling density coupled mass transport in groundwater. Although several models utilising the fractional step method had been developed previously, all were based on the Eulerian solution approach. The model developed by the authors uses the Langrangian approach which has some inherent advantages and disadvantages. The problems associated with the implementation of the fractional step method and techniques by which they were overcome are discussed. The performance of the model is examined and results obtained for standard problems are compared with those from other computer packages.

The NIRVANA project is concerned with the development of a nonoscillatory, integrally reconstructed, volume‐averaged numerical advection scheme. The conservative, flux‐based finite‐volume algorithm is built on an explicit, single‐step, forward‐in‐time update of the cell‐average variable, without restrictions on the size of the time‐step. There are similarities with semi‐Lagrangian schemes; a major difference is the introduction of a discrete integral variable, guaranteeing conservation. The crucial step is the interpolation of this variable, which is used in the calculation of the fluxes; the (analytic) derivative of the interpolant then gives sub‐cell behaviour of the advected variable. In this paper, basic principles are described, using the simplest possible conditions: pure one‐dimensional advection at constant velocity on a uniform grid. Piecewise Nth‐degree polynomial interpolation of the discrete integral variable leads to an Nth‐order advection scheme, in both space and time. Nonoscillatory results correspond to convexity preservation in the integrated variable, leading naturally to a large‐Δt generalisation of the universal limited. More restrictive TVD constraints are also extended to large Δt. Automatic compressive enhancement of step‐like profiles can be achieved without exciting “stair‐casing”. One‐dimensional simulations are shown for a number of different interpolations. In particular, convexity‐limited cubic‐spline and higher‐order polynomial schemes give very sharp, nonoscillatory results at any Courant number, without clipping of extrema. Some practical generalisations are briefly discussed.

The purpose of this paper is to conduct the accurate analysis and systematic characterisation of realistic generalised bi‐isotropic and lossy chiral metamaterial 3D applications at microwave frequencies.

An accuracy‐adjustable time‐domain methodology is developed. The technique uses a convex combination of optimal stencils along with an advanced wavefield decomposition to precisely model the highly dispersive, double negative nature of chirality and constitutive parameters. Furthermore, open‐region radiation or scattering problems are terminated through a pertinently modified perfectly matched layer (PML) of variable depth.

The paper reveals that the proposed algorithm is versatile in the generation of adaptive stencils that attain a very natural way of manipulating continuity conditions at material interfaces. Thus, when periodic structures with split‐ring resonators are to be modelled, the resulting schemes attain optimal precision and minimised dispersion errors. Numerical validation proves these merits via diverse demanding structures of curved shape and multiple layers.

The new technique introduces a family of piecewise polynomials and spatial discretization criteria which lead to additional degrees of freedom for the discrete vectors of the application. In this manner, grid dual is intrinsically embedded in the physical profile of the problem, without resorting to the simplified conventions of other approaches. Moreover, singularity points or demanding geometric discontinuities are properly manipulated, even via coarse lattice resolutions. Thus, the overall accuracy is significantly improved and the computational requirements remain in very logical and affordable levels.

A dust cloud is formed by a high-pressure air blast in dust explosion experiments in the spherical 20 L chamber. The state of the dust cloud has a significant impact on the dust explosion. However, it is difficult to observe the dust distribution in the chamber during the dust dispersion. Numerical simulation was used to examine the dust distribution in the chamber with the rebound nozzle in this work. The paper aims to discuss these issues.

Through a series numerical simulations, the influences of the dust particle size and the pressure for dust dispersion on the have been analyzed, and the results are discussed.

Dust in the spherical 20 L chamber is in the state of very intensifying motion within 30 ms from dispersion starting. Dust in the chamber reaches a uniform state beyond 200 ms. The pressure for dust dispersion should be higher than 0.5 MPa for the aluminum dusts of larger than 50. The higher blast pressure is not always applicable to achieve a uniform dispersion. There is a best blast pressure value for a given dust to achieve a uniform dispersion in the spherical 20 L chamber.

Dust cloud generation is essential for understanding dust explosions. Dust cloud deflagration parameters depend on the uniformity and concentration of dusts dispersed by a high-pressure air blast. Numerical simulation was used to examine the multiphase flow of the dust air mixture in this work. Through a series numerical simulations, the influences of the dust particle size and the pressure for dust dispersion on the have been analyzed, and the results are discussed. The data are useful for understanding the basics of dust cloud formation.

The data are useful for evaluating dust explosion experimental parameters.

Dispersible uniformity has a strong impact on measured parameters of dust explosion in a chamber. However, it is difficult to observe the dust particles distribution during the dispersion. Numerical simulation was used to examine the dust particles distribution and its influencing factors during the dispersion in this work. New finding is: the approach to examine the distribution of dust particles dispersed by a high-pressure blast in a chamber; the variation of dispersible uniformity and its influencing factors when dust is injected into the spherical 20 L chamber by high-pressure air blast.

The dispersive behaviour of waves in softening problems is analysed. Attention is focused on the influence of the numerical scheme on the dispersion characteristics in the process of localization of deformation. Distinction has been made between softening models defined in a standard plasticity framework and in a gradient‐dependent plasticity theory. Waves in a standard softening plasticity continuum do not disperse but due to spatial discretization dispersion is introduced which results in a mesh size dependent length scale effect. On the other hand, wave propagation in a gradient‐dependent softening plasticity continuum is dispersive. By carrying out the dispersion analysis on the discretized system the influence of numerical dispersion on material dispersion can be quantified which enables us to determine the accuracy for the solution of the localization zone. For a modelling with and without the inclusion of strain gradients accuracy considerations with respect to mass discretization, finite element size, time integration scheme and time step have been carried out.

In this work, an SFEM is proposed for solving acoustic problems by redistributing the entries in the mass matrix to “tune” the balance between “stiffness” and “mass” of discrete equation systems, aiming to minimize the dispersion error. The paper aims to discuss this issue.

This is done by simply shifting the four integration points’ locations when computing the entries of the mass matrix in the scheme of SFEM, while ensuring the mass conservation. The proposed method is devised for bilinear quadratic elements.

The balance between “stiffness” and “mass” of discrete equation systems is critically important in simulating wave propagation problems such as acoustics. A formula is also derived for possibly the best mass redistribution in terms of minimizing dispersion error reduction. Both theoretical and numerical examples demonstrate that the present method possesses distinct advantages compared with the conventional SFEM using the same quadrilateral mesh.

After introducing the mass-redistribution technique, the magnitude of the leading relative dispersion error (the quadratic term) of MR-SFEM is bounded by (5/8), which is much smaller than that of original SFEM models with traditional mass matrix (13/4) and consistence mass matrix (2). Owing to properly turning the balancing between stiffness and mass, the MR-SFEM achieves higher accuracy and much better natural eigenfrequencies prediction than the original SFEM does.