Robin-type wall functions and their numerical implementation

The paper is devoted to numerical implementation of the wall functions of Robin-type for modeling near-wall turbulent flows. The wall functions are based on the transfer of a boundary condition from a wall to some intermediate boundary near the wall. The boundary conditions on the intermediate boundary are of Robin-type and represented in a differential form. The wall functions are formulated in an analytical easy-to-implement form, can take into account the source terms, and do not include free parameters. The relation between the wall functions of Robin type and the theory of Calderon–Ryaben'kii's potentials is demonstrated. A universal robust approach to the implementation of the Robin-type wall functions in finite-volume codes is provided. The example of an impinging jet is considered.     

On the boundary of numerical ranges of operators

In this note, a new characterization of an attained boundary point of the numerical range of an operator on a Hilbert space is given. As an application, we point out a gap in the proof of the main result in the paper [M.T. Chien, L. Yeh, On the boundary of the numerical range of a matrix, Appl. Math. Lett. 23 (2010) 725–727].     

Energetic BEM-FEM coupling for the numerical solution of the damped wave equation

Time-dependent problems modeled by hyperbolic partial differential equations can be reformulated in terms of boundary integral equations and solved via the boundary element method. In this context, the analysis of damping phenomena that occur in many physics and engineering problems is a novelty. Starting from a recently developed energetic space-time weak formulation for the coupling of boundary integral equations and hyperbolic partial differential equations related to wave propagation problems, we consider here an extension for the damped wave equation in layered media. A coupling algorithm is presented, which allows a flexible use of finite element method and boundary element method as local discretization techniques. Stability and convergence, proved by energy arguments, are crucial in guaranteeing accurate solutions for simulations on large time intervals. Several numerical benchmarks, whose numerical results confirm theoretical ones, are illustrated and discussed.     

A problem on heat conduction in an insulated wire solved by numerical inverse Laplace transform

A problem of transient heat conduction in an insulated wire is solved by use of Laplace transform and numerical inversion. The problem is solved for the radiation boundary condition and also for the boundary condition of no heat flux through the outer surface of the insulation. The results are presented both numerically with four significant figures and graphically. Asymptotic expansions are derived for small and large values of the time variable. The numerical inversion of the Laplace transform is checked by comparison with the asymptotic expansions and with the numerical results obtained by a numerical inversion formula utilizing one more abscissa than the previous one.     

A numerical method to compute the dissolution of second phases in ternary alloys

Dissolution of stoichiometric multi-component particles in ternary alloys is an important process occurring during the heat treatment of as-cast aluminium alloys prior to hot extrusion. A mathematical model is proposed to describe such a process. In this model an equation is given to determine the position of the particle interface in time, using two diffusion equations which are coupled by nonlinear boundary conditions at the interface. Some results concerning existence, uniqueness, and monotonicity are given. Furthermore, for an unbounded domain an analytical approximation is derived. The main part of this work is the development of a numerical solution method. Finite differences are used on a grid which changes in time. The discretization of the boundary conditions is important to obtain an accurate solution. The resulting nonlinear algebraic system is solved by the Newton-Raphson method. Numerical experiments illustrate the accuracy of the numerical method. The numerical solution is compared with the analytical approximation.     

Smale空间上的广群C~*-代数的迹

本文证明由拓扑混合的Smale空间上的渐进等价关系定义的广群C*-代数及其相应的Ruelle代数有唯一的迹态;在拓扑可迁的情形下,证明此C*-代数的迹态构成了一个单形,此单形顶点的个数等于“Smale谱分解”中基本空间的个数,单形的重心是该C*-代数的唯一的αa-不变迹态;此回答了I.Putnam的一个猜测.     

A numerical method for solving an inverse thermoacoustic problem

In this paper, we consider an inverse problem of determining the initial condition of an initial boundary value problem for the wave equation with some additional information about solving a direct initial boundary value problem. The information is obtained from measurements at the boundary of the solution domain. The purpose of our paper is to construct a numerical algorithm for solving the inverse problem by an iterative method called a method of simple iteration (MSI) and to study the resolution quality of the inverse problem as a function of the number and location of measurement points. Three two-dimensional inverse problem formulations are considered. The results of our numerical calculations are presented. It is shown that the MSI decreases the objective functional at each iteration step. However, due to the ill-posedness of the inverse problem the difference between the exact and approximate solutions decreases up to some fixed number k _min, and then monotonically increases. This shows the regularizing properties of the MSI, and the iteration number can be considered a regularization parameter.     

Self-adaptive FEM numerical modeling of the mild-slope equation

Based on the linear wave theory, the mild-slope equation (MSE) is a preferred mathematical model to simulate nearshore wave propagation. A numerical model to solve the MSE is developed here on the basis of a self-adaptive finite element model (FEM) combined with an iterative method to determine the wave direction angle to the boundary and thus to improve the treatment of the boundary conditions. The numerical resolution of the waves into ideal domains and multidirectional waves through a breakwater gap shows that the numerical model developed here is effective in representing wave absorption at the absorbing boundaries and can be used to simulate multidirectional wave propagation. Finally, the simulated wave distribution in a real harbor shows that the numerical model can be used for engineering practice.     

Comparison of several numerical optimization methods

Several numerical methods for solving the nonlinear two-point boundary value problem associated with an optimum spacecraft trajectory are considered. A comparative evaluation of the methods is made to determine the relative merits of each method. Particular attention is given to such characteristics as the simplicity of formulation and implementation, the convergence sensitivity, the computing time required, and the computer storage requirements. The methods considered are the perturbation method, the quasilinearization method, and the gradient method. The numerical comparison is made by considering a two-dimensional, low-thrust, minimum-time, Earth-Mars trajectory.The authors are greatly indebted to Mr. Robert D. Witty, Lockheed Electronics Company, for providing the excellent programming support.     

An accurate numerical solution for the transient heating of an evaporating spherical droplet

A recently derived numerical algorithm for one-dimensional time-dependent Stefan problems is extended for the purposes of solving a moving boundary problem for the transient heating of an evaporating spherical droplet. The Keller box finite-difference scheme is used, in tandem with the so-called boundary immobilization method. An important component of the work is the careful use of variable transformations that must be built into the numerical algorithm in order to preserve second-order accuracy in both time and space - an issue not previously discussed in relation to this widely-used scheme. In addition, we demonstrate that our solution is in close agreement with the solution obtained using an alternative numerical scheme that employs an analytic solution of the heat conduction equation inside the droplet, for which the droplet radius was assumed to be a piecewise linear function of time. The advantages of the new method are discussed.     

Approximate method for the numerical solution of singular perturbation problems

We present an approximate method for the numerical solution of linear singularly perturbed two point boundary value problems in ordinary differential equations with a boundary layer on the left end of the underlying interval. It is motivated by the asymptotic behavior of singular perturbation problems. The original problem is divided into inner and outer region problems. The reduced problem is solved to obtain the terminal boundary condition. Then, a new inner region problem is created and solved as a two point boundary value problem. In turn, the outer region problem is also modified and the resulting problem is efficiently treated by employing the trapezoidal formula coupled with discrete invariant imbedding algorithm. The proposed method is iterative on the terminal point. Some numerical experiments have been included to demonstrate its applicability.     

A numerical algorithm for nonlinear multi-point boundary value problems

In this paper, an algorithm is presented for solving second-order nonlinear multi-point boundary value problems (BVPs). The method is based on an iterative technique and the reproducing kernel method (RKM). Two numerical examples are provided to show the reliability and efficiency of the present method.     

On numerical blow-up sets

In this paper we study numerical blow-up sets for semidicrete approximations of the heat equation with nonlinear boundary conditions. We prove that the blow-up set either concentrates near the boundary or is the whole domain.

     

A globally convergent numerical method for an inverse elliptic problem of optical tomography

A new globally convergent numerical method is developed for an inverse problem for the elliptic equation with the unknown potential. The boundary data simulating measurements in optical tomography are generated by the running source. Global convergence analysis is presented along with numerical experiments.     

Crossed products by finite group actions with certain tracial Rokhlin property

We introduce a special tracial Rokhlin property for unital C~*-algebras. Let A be a unital tracial rank zero C~*-algebra(or tracial rank no more than one C~*-algebra). Suppose that α : G → Aut(A) is an action of a finite group G on A, which has this special tracial Rokhlin property, and suppose that A is a α-simple C~*-algebra. Then, the crossed product C~*-algebra C~*(G, A, α) has tracia rank zero(or has tracial rank no more than one). In fact,we get a more general results.     

Convergence and stability of a numerical method for micromagnetics

The convergence and stability of a numerical method, which applies a nonconforming finite element method and an artificial boundary method to a multi-atomic Young measure relaxation model, for micromagnetics are analyzed. By revealing some key properties of the solution sets of both the continuous and discrete problems, we show that our numerical method is stable, and the solution set of the continuous problem is well approximated by those of the discrete problems. The performance of our method is also illustrated by some numerical examples. The research was supported in part by the Major State Basic Research Projects (2005CB321701), NSFC projects (10431050, 10571006, 10528102 and 10871011) and RFDP of China.     

Stability under Small Hilbert–Schmidt Perturbations for <Emphasis Type="Italic">C</Emphasis><Superscript>*</Superscript>-Algebras

This paper studies the tracial stability of C^*-algebras, which is a general property of stability of relations in a Hilbert–Schmidt-type norm defined by a trace on a C^*-algebra. Precise definitions are formulated in terms of tracial ultraproducts. For nuclear C^*-algebras, a characterization of matricial tracial stability in terms of approximation of tracial states by traces of finite-dimensional representations is obtained. For the nonnuclear case, new obstructions and counterexamples are constructed in terms of free entropy theory.     

Anisotropic mesh adaptation for numerical solution of boundary value problems

We present an efficient mesh adaptation algorithm that can be successfully applied to numerical solutions of a wide range of 2D problems of physics and engineering described by partial differential equations. We are interested in the numerical solution of a general boundary value problem discretized on triangular grids. We formulate a necessary condition for properties of the triangulation on which the discretization error is below the prescribed tolerance and control this necessary condition by the interpolation error. For a sufficiently smooth function, we recall the strategy how to construct the mesh on which the interpolation error is below the prescribed tolerance. Solving the boundary value problem we apply this strategy to the smoothed approximate solution. The novelty of the method lies in the smoothing procedure that, followed by the anisotropic mesh adaptation (AMA) algorithm, leads to the significant improvement of numerical results. We apply AMA to the numerical solution of an elliptic equation where the exact solution is known and demonstrate practical aspects of the adaptation procedure: how to control the ratio between the longest and the shortest edge of the triangulation and how to control the transition of the coarsest part of the mesh to the finest one if the two length scales of all the triangles are clearly different. An example of the use of AMA for the physically relevant numerical simulation of a geometrically challenging industrial problem (inviscid transonic flow around NACA0012 profile) is presented. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004.     

Normality and the numerical range

It is well known that if A is an n by n normal matrix, then the numerical range of A is the convex hull of its spectrum. The converse is valid for n ? 4 but not for larger n. In this spirit a characterization of normal matrices is given only in terms of the numerical range. Also, a characterization is given of matrices for which the numerical range coincides with the convex hull of the spectrum. A key observation is that the eigenvectors corresponding to any eigenvalue occuring on the boundary of the numerical range must be orthogonal to eigenvectors corresponding to all other eigenvalues.