On the modelling of solid state reactions.Synthesis of YAG

There is a model of yttrium aluminium garnet (YAG) synthesis presented in this article. The developed model is based on nonlinear reaction–diffusion partial differential equations. The solution was carried out numerically using finite difference techniques. We got dependability curves for diffusion and reaction rates and offered possible method to localize values of diffusion and reaction rate constants precisely enough.AMS subject classification: 35K57, 65M06     

On obtaining localized bonding schemes from delocalized molecular-orbital wave functions

A straightforward procedure is proposed for expanding a molecular orbital determinantal wave function into a set of determinantal wave functions composed of atomic orbitals localized at the atoms of a molecule. By employing this method, atomic orbital determinants and their weights can be derived for a molecule from the computed molecular-orbital wave function. The procedure permits the interpretation of a molecular orbital determinantal wave function in terms of bonding schemes related to the classic resonance structures used by organic chemists. By using the unrestricted molecular orbital determinant, bonding schemes and their weights are obtained for butadiene, the butadiene radical cation and the acrylonitrile radical anion. Their dominant bonding schemes are in accord with the relevant resonance structures for these molecules. For the butadiene radical cation and the acrylonitrile anion they are shown to be compatible with the accepted mechanisms of the electrochemical coupling reactions of butadiene and acrylonitrile. Received: 7 August 1996 / Accepted: 18 March 1997     

Use of Local Lagrange Interpolation for calculation of retention indices in linear temperature-programmed gas chromatrography

Summary The influence of the isothermal temperature, program rate, initial temperature and flow rate on retention indices was studied. The methods of Kováts, Van Den Dool and Local Lagrange Interpolation are compared. Ten experimental measurements were carried out on a capillary column coated with OV-101 stationary phase.     

Trivariate Polynomial Natural Spline for 3D Scattered Data Hermit Interpolation

Consider a kind of Hermit interpolation for scattered data of 3D by trivariate polynomial natural spline,such that the objective energy functional (with natural boundary conditions) is minimal.By the s...     

A BGK‐penalization‐based asymptotic‐preserving scheme for the multispecies Boltzmann equation

An asymptotic‐preserving (AP) scheme is efficient in solving multiscale problems where kinetic and hydrodynamic regimes coexist. In this article, we extend the BGK‐penalization‐based AP scheme, originally introduced by Filbet and Jin for the single species Boltzmann equation (Filbet and Jin, J Comput Phys 229 (2010) 7625–7648), to its multispecies counterpart. For the multispecies Boltzmann equation, the new difficulties arise due to: (1) the breaking down of the conservation laws for each species and (2) different convergence rates to equilibria for different species in disparate masses systems. To resolve these issues, we find a suitable penalty function—the local Maxwellian that is based on the mean velocity and mean temperature and justify various asymptotic properties of this method. This AP scheme does not contain any nonlinear nonlocal implicit solver, yet it can capture the fluid dynamic limit with time step and mesh size independent of the Knudsen number. Numerical examples demonstrate the correct asymptotic‐behavior of the scheme. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Linear stability of WENO schemes coupled with explicit Runge–Kutta schemes

This paper focuses on the results of the linear stability analysis of the finite‐difference weighted essentially non‐oscillatory (WENO) schemes with optimal weights. The standard WENO schemes between the third and 11th order, the order‐optimised WENO schemes of the sixth and eighth order and the bandwidth‐optimised WENO schemes of the third and fourth order are considered. Several explicit Runge–Kutta schemes including the recently published strong stability‐preserving explicit Runge–Kutta schemes are considered for time discretisation. The stability limits as well as dissipation and dispersion properties dependent on the Courant–Friedrichs–Lewy number are presented for a hyperbolic model equation. The different combinations of space and time discretisation schemes are compared in terms of their accuracy and efficiency. For a parabolic model equation, the viscous term is discretised with high‐order central differences. The stability limits for the parabolic problem are presented as well. Numerical results of linear test cases are shown. Copyright © 2011 John Wiley & Sons, Ltd.     

A thermodynamic and dynamic subgrid closure model for two‐material cells

A general and robust subgrid closure model for two‐material cells is proposed. The conservative quantities of the entire cell are apportioned between two materials, and then, pressure and velocity are fully or partially equilibrated by modeling subgrid wave interactions. An unconditionally stable and entropy‐satisfying solution of the processes has been successfully found. The solution is valid for arbitrary level of relaxation. The model is numerically designed with care for general materials and is computationally efficient without recourse to subgrid iterations or subcycling in time. The model is implemented and tested in the Lagrange‐remap framework. Two interesting results are observed in 1D tests. First, on the basis of the closure model without any pressure and velocity relaxation, a material interface can be resolved without creating numerical oscillations and/or large nonphysical jumps in the problem of the modified Sod shock tube. Second, the overheating problem seen near the wall surface can be solved by the present entropy‐satisfying closure model. The generality, robustness, and efficiency of the model make it useful in principle in algorithms, such as ALE methods, volume of fluid methods, and even some mixture models, for compressible two‐phase flow computations. Copyright © 2013 John Wiley & Sons, Ltd.     

On the signal degradation induced by the interpolation and the sampling rate reduction in aeroacoustics hybrid methods

Taking place within more extensive work that focuses on hybrid methods in aeroacoustics, the present study is devoted to the data transfer operations that are to be performed between two stages of a hybrid calculation. More precisely, the article focuses on two typical operations that usually accompany such data transfer, which are (i) the sampling rate reduction and (ii) the interpolation of the unsteady perturbations to be transmitted from one stage to the other. First part of the paper analyzes the two main issues of such operations, which are the spuriousing and the aliasing phenomena. For doing so, the usual notions of the interpolation theory are revisited before they are synthesized within an original approach. The here proposed formalism allows to understand better both the spuriousing and the aliasing phenomena, as well as to accurately predict the impact of the latter on the data to be transmitted in terms of signal degradation. Second part of the paper provides an illustration and a validation of these theoretical developments via a direct application to a typical aerodynamic noise problem (aeroacoustic emission by a 2D cylinder cross flow). There, it is further shown how the here proposed formalism can help in improving aeroacoustics hybrid calculations by predicting (and thus possibly minimizing) the bias to be induced on the acoustic extrapolation stage because of the aliasing and/or spuriousing effects inherited from the sampling rate reduction and/or interpolation of CFD data—which is likely to occur in any hybrid scenario. Copyright © 2012 John Wiley & Sons, Ltd.     

A numerical analysis of a class of generalized Boussinesq‐type equations using continuous/discontinuous FEM

An improved class of Boussinesq systems of an arbitrary order using a wave surface elevation and velocity potential formulation is derived. Dissipative effects and wave generation due to a time‐dependent varying seabed are included. Thus, high‐order source functions are considered. For the reduction of the system order and maintenance of some dispersive characteristics of the higher‐order models, an extra O(μ²ⁿ⁺²) term (n ∈ ?) is included in the velocity potential expansion. We introduce a nonlocal continuous/discontinuous Galerkin FEM with inner penalty terms to calculate the numerical solutions of the improved fourth‐order models. The discretization of the spatial variables is made using continuous P₂ Lagrange elements. A predictor‐corrector scheme with an initialization given by an explicit Runge–Kutta method is also used for the time‐variable integration. Moreover, a CFL‐type condition is deduced for the linear problem with a constant bathymetry. To demonstrate the applicability of the model, we considered several test cases. Improved stability is achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

Kinematics of layered reinforced-concrete planar beam finite elements with embedded transversal cracking

In this work crack formation and development is addressed and implemented in a planar layered reinforced-concrete beam element. The crack initiation and growth is described using the strength criterion in conjunction with exact kinematics of the interlayer connection. In this way a novel embedded-discontinuity beam finite element is derived in which the tensile stresses in concrete at the crack position reaching the tensile strength will trigger a crack to open. Since the element is multi-layered, in this way the crack is allowed to propagate through the depth of the beam. The cracked layer(s) will involve discontinuity in the cross-sectional rotation equal to the crack-profile angle, as well as a discontinuity in the position vector of the layer’s reference line. A bond–slip relationship is superimposed onto this model in a kinematically consistent manner with reinforcement being treated as an additional layer of zero thickness with its own material parameters and a constitutive law implemented in the multi-layered beam element.