Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

Eigenfunctions and fundamental solutions of the fractional Laplace and Dirac operators: the Riemann-Liouville case

In this paper we study eigenfunctions and fundamental solutions for the three parameter fractional Laplace operator $\Delta_+^{(\alpha,\beta,\gamma)}:= D_{x_0^+}^{1+\alpha} +D_{y_0^+}^{1+\beta} +D_{z_0^+}^{1+\gamma},$ where $(\alpha, \beta, \gamma) \in \,]0,1]^3$, and the fractional derivatives $D_{x_0^+}^{1+\alpha}$, $D_{y_0^+}^{1+\beta}$, $D_{z_0^+}^{1+\gamma}$ are in the Riemann-Liouville sense. Applying operational techniques via two-dimensional Laplace transform we describe a complete family of eigenfunctions and fundamental solutions of the operator $\Delta_+^{(\alpha,\beta,\gamma)}$ in classes of functions admitting a summable fractional derivative. Making use of the Mittag-Leffler function, a symbolic operational form of the solutions is presented. From the obtained family of fundamental solutions we deduce a family of fundamental solutions of the fractional Dirac operator, which factorizes the fractional Laplace operator. We apply also the method of separation of variables to obtain eigenfunctions and fundamental solutions.

Eigenfunctions and fundamental solutions of the Caputo fractional Laplace and Dirac operators

In this paper, by using the method of separation of variables, we obtain eigenfunctions and fundamental solutions for the three parameter fractional Laplace operator defined via fractional Caputo derivatives. The solutions are expressed using the Mittag-Leffler function and we show some graphical representations for some parameters. A family of fundamental solutions of the corresponding fractional Dirac operator is also obtained. Particular cases are considered in both cases.

A meshless method for an inverse two-phase one-dimensional nonlinear Stefan problem

We extend a meshless method of fundamental solutions recently proposed by the authors for the one-dimensional two-phase inverse linear Stefan problem, to the nonlinear case. In this latter situation the free surface is also considered unknown which is more realistic from the practical point of view. Building on the earlier work, the solution is approximated in each phase by a linear combination of fundamental solutions to the heat equation. The implementation and analysis are more complicated in the present situation since one needs to deal with a nonlinear minimization problem to identify the free surface. Furthermore, the inverse problem is ill-posed since small errors in the input measured data can cause large deviations in the desired solution. Therefore, regularization needs to be incorporated in the objective function which is minimized in order to obtain a stable solution. Numerical results are presented and discussed. © 2014 IMACS.

The development of a wax layer on the interior wall of a circular pipe transporting heated oil

In this article we develop a simple model to describe the evolution of a depositional wax layer on the inner surface of a circular pipe transporting heated oil, which contains dissolved wax. When the outer pipe surface is cooled sufficiently, the growth of a wax layer is initiated on the inner pipe wall, and this evolves to a saturated steady state thickness. The model proposed is based on fundamental balances of heat flow from the oil, into the wax layer, and across the pipe wall. We present an analysis of the model, examine a relevant asymptotic limit in which the full details of the solution to the model are available and develop an efficient numerical method (based on the method of fundamental solutions) for producing approximations of the model solution. The mathematical structure of the model is that of a free boundary evolution problem of generalised Stefan type. © The Author, 2014.

The development of a wax layer on the interior wall of a circular pipe transporting heated oil

In this article we develop a simple model to describe the evolution of a depositional wax layer on the inner surface of a circular pipe transporting heated oil, which contains dissolved wax. When the outer pipe surface is cooled sufficiently, the growth of a wax layer is initiated on the inner pipe wall, and this evolves to a saturated steady state thickness. The model proposed is based on fundamental balances of heat flow from the oil, into the wax layer, and across the pipe wall. We present an analysis of the model, examine a relevant asymptotic limit in which the full details of the solution to the model are available and develop an efficient numerical method (based on the method of fundamental solutions) for producing approximations of the model solution. The mathematical structure of the model is that of a free boundary evolution problem of generalised Stefan type. © The Author, 2014.

Identification of a time-dependent bio-heat blood perfusion coefficient

In the paper the identification of the time-dependent blood perfusion coefficient is formulated as an inverse problem. The bio-heat conduction problem is transformed into the classical heat conduction problem. Then the transformed inverse problem is solved using the method of fundamental solutions together with the Tikhonov regularization. Some numerical results are presented in order to demonstrate the accuracy and the stability of the proposed meshless numerical algorithm.

Stability of periodic solutions obtained via the averaging method for nonsmooth Lipschitz system

"Vegeu el resum a l'inici del document del fitxer adjunt."

Potentiometric titration as a method of determining the stoichiometric dissociation constants of weak acids from the thermodynamic dissociation constants of aqueous solutions at 298.15 K

The simple single-ion activity coefficient equation originating from the Debye-Hückel theory was used to determine the thermodynamic and stoichiometric dissociation constants of weak acids from data concerning galvanic cells. Electromotive force data from galvanic cells without liquid junctions, which was obtained from literature, was studied in conjuction with the potentiometric titration data relating to aqueous solutions at 298.15 K. The dissociation constants of weak acids could be determined by the presented techniques and almost all the experimental data studied could be interpreted within the range of experimental error. Potentiometric titration has been used here and the calculation methods were developed to obtain the thermodynamic and stoichiometric dissociation constants of some weak acids in aqueous solutions at 298.15 K. The ionic strength of titrated solutions were adjusted using an inert electrolyte, namely, sodium or potassium chloride. Salt content alonedetermines the ionic strength. The ionic strength of the solutions studied varied from 0.059 mol kg-1 to 0.37 mol kg-1, and in some cases up to 1.0 mol kg-1. The following substances were investigated using potentiometric titration: aceticacid, propionic acid, L-aspartic acid, L-glutamic acid and bis(2,2-dimethyl-3-oxopropanol) amine.

Chelating Adsorbents in Purification of Hydrometallurgical Solutions

In this thesis, equilibrium and dynamic sorption properties of weakly basic chelating adsorbents were studied to explain removal of copper, nickel from a concentrated zinc sulfate solution in a hydrometallurgical process. Silica-supported chelating composites containing either branched poly(ethyleneimine) (BPEI) or 2-(aminomethyl)pyridine (AMP) as a functional group were used. The adsorbents are commercially available from Purity Systems Inc, USA as WP-1® and CuWRAM®, respectively. The fundamental interactions between the adsorbents, sulfuric acid and metal sulfates were studied in detail and the results were used to find the best conditions for removal of copper and nickel from an authentic ZnSO4 process solution. In particular, the effect of acid concentration and temperature on the separation efficiency was considered. Both experimental and modeling aspectswere covered in all cases. Metal sorption is considerably affected by the chemical properties of the studied adsorbents and by the separation conditions. In the case of WP-1, acid affinity is so high that column separation of copper, nickel and zinc has to be done using the adsorbent in base-form. On the other hand, the basicity of CuWRAM is significantly lower and protonated adsorbent can be used. Increasing temperature decreases the basicity and the metals affinity of both adsorbents, but the uptake capacities remain practically unchanged. Moreover, increasing temperature substantially enhances intra-particle mass transport and decreases viscosities thus allowing significantly higher feed flow rates in the fixed-bed separation. The copper selectivity of both adsorbents is very high even in the presence of a 250-fold excess of zinc. However, because of the basicity of WP-1, metal precipitation is a serious problem and therefore only CuWRAM is suitable for the practical industrial application. The optimum temperature for copper removal appears to be around 60 oC and an alternative solution purification method is proposed. The Ni/Zn selectivity of both WP-1 and CuWRAM is insufficient for removal of the very small amounts of nickel present in the concentrated ZnSO4 solution.

Fast and accurate numerical solutions in some problems of particle and radiation transport: synthetic acceleration for the method of short characteristics, Doppler-broadened scattering kernel, remote sensing of the cryosphere

The aim of this work is to present various aspects of numerical simulation of particle and radiation transport for industrial and environmental protection applications, to enable the analysis of complex physical processes in a fast, reliable, and efficient way. In the first part we deal with speed-up of numerical simulation of neutron transport for nuclear reactor core analysis. The convergence properties of the source iteration scheme of the Method of Characteristics applied to be heterogeneous structured geometries has been enhanced by means of Boundary Projection Acceleration, enabling the study of 2D and 3D geometries with transport theory without spatial homogenization. The computational performances have been verified with the C5G7 2D and 3D benchmarks, showing a sensible reduction of iterations and CPU time. The second part is devoted to the study of temperature-dependent elastic scattering of neutrons for heavy isotopes near to the thermal zone. A numerical computation of the Doppler convolution of the elastic scattering kernel based on the gas model is presented, for a general energy dependent cross section and scattering law in the center of mass system. The range of integration has been optimized employing a numerical cutoff, allowing a faster numerical evaluation of the convolution integral. Legendre moments of the transfer kernel are subsequently obtained by direct quadrature and a numerical analysis of the convergence is presented. In the third part we focus our attention to remote sensing applications of radiative transfer employed to investigate the Earth's cryosphere. The photon transport equation is applied to simulate reflectivity of glaciers varying the age of the layer of snow or ice, its thickness, the presence or not other underlying layers, the degree of dust included in the snow, creating a framework able to decipher spectral signals collected by orbiting detectors.

Editorial reorganization of an alphabetically ordered collection of short stories: analysis of problems and hypothesis of possible solutions

This paper aims to show, analyze and solve the problems related to the translation of the book with meaning-bond alphabetically ordered chapter titles La vita non è in ordine alfabetico, by the Italian writer Andrea Bajani. The procedure is inevitable for a possible translation of the book, and it is necessary to have a preliminary pattern to follow. After translating the whole book, not only is a revision fundamental, but a restructure and reorganization of the collection may be required, hence, what this thesis offers is a scheme to start from, together with an analysis of the possible problems that may arise, and a useful method to find a solution.

A differential method for the determination of partial molal heat capacities and specific heats of dilute solutions /

Thesis (Ph.D.)--University of California, Berkeley, 1922.