This paper presents two uniformly convergent numerical schemes for the two dimensional steady state discrete ordinates transport equation in the diffusive regime, which is valid up to the boundary and interface layers. A five-point node-centered and a four-point cell-centered tailored finite point schemes (TFPS) are introduced. The schemes first approximate the scattering coefficients and sources by piecewise constant functions and then use special solutions to the constant coefficient equation as local basis functions to formulate a discrete linear system. Numerically, both methods can not only capture the diffusion limit, but also exhibit uniform convergence in the diffusive regime, even with boundary layers. Numerical results show that the five-point scheme has first-order accuracy and the four-point scheme has second-order accuracy, uniformly with respect to the mean free path. Therefore a relatively coarse grid can be used to capture the two dimensional boundary and interface layers.

A direct-forcing immersed boundary method (DFIB) with both virtual force and heat source is developed here to solve Navier-Stokes and the associated energy transport equations to study some thermal flow problems caused by a moving rigid solid object within. The key point of this novel numerical method is that the solid object, stationary or moving, is first treated as fluid governed by Navier-Stokes equations for velocity and pressure, and by energy transport equation for temperature in every time step. An additional virtual force term is then introduced on the right hand side of momentum equations in the solid object region to make it act exactly as if it were a solid rigid body immersed in the fluid. Likewise, an additional virtual heat source term is applied to the right hand side of energy equation at the solid object region to maintain the solid object at the prescribed temperature all the time. The current method was validated by some benchmark forced and natural convection problems such as a uniform flow past a heated circular cylinder, and a heated circular cylinder inside a square enclosure. We further demonstrated this method by studying a mixed convection problem involving a heated circular cylinder moving inside a square enclosure. Our current method avoids the otherwise requested dynamic grid generation in traditional method and shows great efficiency in the computation of thermal and flow fields caused by fluid-structure interaction.

Multiplicative noise removal is a challenging problem in image restoration. In this paper, by applying Box-Cox transformation, we convert the multiplicative noise removal problem into the additive noise removal problem and the block matching three dimensional (BM3D) method is applied to get the final recovered image. Indeed, BM3D is an effective method to remove additive Gaussian white noise in images. A maximum likelihood method is designed to determine the parameter in the Box-Cox transformation. We also present the unbiased inverse transform for the Box-Cox transformation which is important. Both theoretical analysis and experimental results illustrate clearly that the proposed method can remove multiplicative noise very well especially when multiplicative noise is heavy. The proposed method is superior to the existing methods for multiplicative noise removal in the literature.

In this paper the numerical method for solution of an aeroelastic model describing the interactions of air flow with vocal folds is described. The flow is modelled by the incompressible Navier-Stokes equations spatially discretized with the aid of the stabilized finite element method. The motion of the computational domain is treated with the aid of the Arbitrary Lagrangian Eulerian method. The structure dynamics is replaced by a mechanically equivalent system with the two degrees of freedom governed by a system of ordinary differential equations and discretized in time with the aid of an implicit multistep method and strongly coupled with the flow model. The influence of inlet/outlet boundary conditions is studied and the numerical analysis is performed and compared to the related results from literature.

The high order inverse Lax-Wendroff (ILW) procedure is extended to boundary treatment involving complex geometries on a Cartesian mesh. Our method ensures that the numerical resolution at the vicinity of the boundary and the inner domain keeps the fifth order accuracy for the system of the reactive Euler equations with the two-step reaction model. Shock wave propagation in a tube with an array of rectangular grooves is first numerically simulated by combining a fifth order weighted essentially non-oscillatory (WENO) scheme and the ILW boundary treatment. Compared with the experimental results, the ILW treatment accurately captures the evolution of shock wave during the interactions of the shock waves with the complex obstacles. Excellent agreement between our numerical results and the experimental ones further demonstrates the reliability and accuracy of the ILW treatment. Compared with the immersed boundary method (IBM), it is clear that the influence on pressure peaks in the reflected zone is obviously bigger than that in the diffracted zone. Furthermore, we also simulate the propagation process of detonation wave in a tube with three different widths of wall-mounted rectangular obstacles located on the lower wall. It is shown that the shock pressure along a horizontal line near the rectangular obstacles gradually decreases, and the detonation cellular size become large and irregular with the decrease of the obstacle width.

A conservative modification to the ghost fluid method (GFM) is developed for compressible multiphase flows. The motivation is to eliminate or reduce the conservation error of the GFM without affecting its performance. We track the conservative variables near the material interface and use this information to modify the numerical solution for an interfacing cell when the interface has passed the cell. The modification procedure can be used on the GFM with any base schemes. In this paper we use the fifth order finite difference WENO scheme for the spatial discretization and the third order TVD Runge-Kutta method for the time discretization. The level set method is used to capture the interface. Numerical experiments show that the method is at least mass and momentum conservative and is in general comparable in numerical resolution with the original GFM.

In this paper we show that entropy can be used within a functional for the stress relaxation time of solid materials to parametrise finite viscoplastic strain-hardening deformations. Through doing so the classical empirical recovery of a suitable irreversible scalar measure of work-hardening from the three-dimensional state parameters is avoided. The success of the proposed approach centres on determination of a rate-independent relation between plastic strain and entropy, which is found to be suitably simplistic such to not add any significant complexity to the final model. The result is sufficiently general to be used in combination with existing constitutive models for inelastic deformations parametrised by one-dimensional plastic strain provided the constitutive models are thermodynamically consistent. Here a model for the tangential stress relaxation time based upon established dislocation mechanics theory is calibrated for OFHC copper and subsequently integrated within a two-dimensional moving-mesh scheme. We address some of the numerical challenges that are faced in order to ensure successful implementation of the proposedmodel within a hydrocode. The approach is demonstrated through simulations of flyer-plate and cylinder impacts.

A multigrid method is proposed to compute the ground state solution of Bose-Einstein condensations by the finite element method based on the multilevel correction for eigenvalue problems and the multigrid method for linear boundary value problems. In this scheme, obtaining the optimal approximation for the ground state solution of Bose-Einstein condensates includes a sequence of solutions of the linear boundary value problems by the multigrid method on the multilevel meshes and some solutions of nonlinear eigenvalue problems some very low dimensional finite element space. The total computational work of this scheme can reach almost the same optimal order as solving the corresponding linear boundary value problem. Therefore, this type of multigrid scheme can improve the overall efficiency for the simulation of Bose-Einstein condensations. Some numerical experiments are provided to validate the efficiency of the proposed method.

We couple different flow models, i.e. a finite element solver for the Navier-Stokes equations and a Lattice Boltzmann automaton, using the framework Peano as a common base. The new coupling strategy between the meso- and macroscopic solver is presented and validated in a 2D channel flow scenario. The results are in good agreement with theory and results obtained in similar works by Latt et al. In addition, the test scenarios show an improved stability of the coupled method compared to pure Lattice Boltzmann simulations.

An unstructured nodal spectral-element method for the Navier-Stokes equations is developed in this paper. The method is based on a triangular and tetrahedral rational approximation and an easy-to-implement nodal basis which fully enjoys the tensorial product property. It allows arbitrary triangular and tetrahedral mesh, affording greater flexibility in handling complex domains while maintaining all essential features of the usual spectral-element method. The details of the implementation and some numerical examples are provided to validate the efficiency and flexibility of the proposed method.

A novel method is developed for solving the inverse problem of reconstructing the shape of an interior cavity. The boundary of the cavity is assumed to be a small and smooth perturbation of a circle. The incident field is generated by a point source inside the cavity. The scattering data is taken on a circle centered at the source. The method requires only a single incident wave at one frequency. Using a transformed field expansion, the original boundary value problem is reduced to a successive sequence of two-point boundary value problems and is solved in a closed form. By dropping higher order terms in the power series expansion, the inverse problem is linearized and an explicit relation is established between the Fourier coefficients of the cavity surface function and the total field. A nonlinear correction algorithm is devised to improve the accuracy of the reconstruction. Numerical results are presented to show the effectiveness of the method and its ability to obtain subwavelength resolution.

A method for gyrokinetic simulation of low frequency (lower than the cyclotron frequency) magnetic compressional modes in general geometry is presented. The gyrokinetic-Maxwell system of equations is expressed fully in terms of the compressional component of the magnetic perturbation, δB∥, with finite Larmor radius effects. This introduces a “gyro-surface” averaging of δB∥ in the gyrocenter equations of motion, and similarly in the perpendicular Ampere’s law, which takes the form of the perpendicular force balance equation. The resulting system can be numerically implemented by representing the gyro-surface averaging by a discrete sum in the configuration space. For the typical wavelength of interest (on the order of the gyroradius), the gyro-surface averaging can be reduced to averaging along an effective gyro-orbit. The phase space integration in the force balance equation can be approximated by summing over carefully chosen samples in the magnetic moment coordinate, allowing for an efficient numerical implementation.

Consider the acoustic wave scattering by an impenetrable obstacle in two dimensions, where the wave propagation is governed by the Helmholtz equation. The scattering problem is modeled as a boundary value problem over a bounded domain. Based on the Dirichlet-to-Neumann (DtN) operator, a transparent boundary condition is introduced on an artificial circular boundary enclosing the obstacle. An adaptive finite element based on a posterior error estimate is presented to solve the boundary value problem with a nonlocal DtN boundary condition. Numerical experiments are included to compare with the perfectly matched layer (PML) method to illustrate the competitive behavior of the proposed adaptive method.

We study compact finite difference methods for the Schrodinger-Poisson equation in a bounded domain and establish their optimal error estimates under proper regularity assumptions on wave function and external potential V(x). The Crank-Nicolson compact finite difference method and the semi-implicit compact finite difference method are both of order Ҩ(h4 + τ2) in discrete l2,H1 and l∞ norms with mesh size h and time step t. For the errors ofcompact finite difference approximation to the second derivative andPoisson potential are nonlocal, thus besides the standard energy method and mathematical induction method, the key technique in analysisis to estimate the nonlocal approximation errors in discrete l∞ and H1 norm by discrete maximum principle of elliptic equation and properties of some related matrix. Also some useful inequalities are established in this paper. Finally, extensive numerical re-sults are reported to support our error estimates of the numerical methods.

The evolution of precipitates in stressed solids is modeled by coupling a quasi-steady diffusion equation and a linear elasticity equation with dynamic boundary conditions. The governing equations are solved numerically using a boundary integral method (BIM). A critical step in applying BIM is to develop fast algorithms to reduce the arithmetic operation count of matrix-vector multiplications. In this paper, we develop a fast adaptive treecode algorithm for the diffusion and elasticity problems in two dimensions (2D). We present a novel source dividing strategy to parallelize the treecode. Numerical results show that the speedup factor is nearly perfect up to a moderate number of processors. This approach of parallelization can be readily implemented in other treecodes using either uniform or non-uniform point distribution. We demonstrate the effectiveness of the treecode by computing the long-time evolution of a complicated microstructure in elastic media, which would be extremely difficult with a direct summation method due to CPU time constraint. The treecode speeds up computations dramatically while fulfilling the stringent precision requirement dictated by the spectrally accurate BIM.

This paper deals with boundary optimal control problems for the heat and Navier-Stokes equations and addresses the issue of defining controls in function spaces which are naturally associated to the volume variables by trace restriction. For this reason we reformulate the boundary optimal control problem into a distributed problem through a lifting function approach. The stronger regularity requirements which are imposed by standard boundary control approaches can then be avoided. Furthermore, we propose a new numerical strategy that allows to solve the coupled optimality system in a robust way for a large number of unknowns. The optimality system resulting from a finite element discretization is solved by a local multigrid algorithm with domain decomposition Vanka-type smoothers. The purpose of these smoothers is to solve the optimality system implicitly over subdomains with a small number of degrees of freedom, in order to achieve robustness with respect to the regularization parameters in the cost functional. We present the results of some test cases where temperature is the observed quantity and the control quantity corresponds to the boundary values of the fluid temperature in a portion of the boundary. The control region for the observed quantity is a part of the domain where it is interesting to match a desired temperature value.

A consistent dual-mesh hybrid LES/RANS framework for turbulence modeling has been proposed recently (H. Xiao, P. Jenny, A consistent dual-mesh framework for hybrid LES/RANS modeling, J. Comput. Phys. 231 (4) (2012)). To better enforce componentwise Reynolds stress consistency between the LES and the RANS simulations, in the present work the original hybrid framework is modified to better exploit the advantage of more advanced RANS turbulence models. In the new formulation, the turbulent stresses in the filtered equations in the under-resolved regions are directly corrected based on the Reynolds stresses provided by the RANS simulation. More precisely, the new strategy leads to implicit LES/RANS consistency, where the velocity consistency is achieved indirectly via imposing consistency on the Reynolds stresses. This is in contrast to the explicit consistency enforcement in the original formulation, where forcing terms are added to the filtered momentum equations to achieve directly the desired average velocity and velocity fluctuations. The new formulation keeps the averaging procedure for the filtered quantities and at the same time preserves the ability of the original formulation to conform with the physical differences between LES and RANS quantities. The modified formulation is presented, analyzed, and then evaluated for plane channel flow and flow over periodic hills. Improved predictions are obtained compared with the results obtained using the original formulation.

This paper presents an approach using the method of separation of variables applied to 2D Helmholtz equations in the Cartesian coordinate. The solution is then computed by a series solutions resulted from solving a sequence of 1D problems, in which the 1D solutions are computed using pollution free difference schemes. Moreover, non-polluted numerical integration formulae are constructed to handle the integration due to the forcing term in the inhomogeneous 1D problems. Consequently, the computed solution does not suffer the pollution effect. Another attractive feature of this approach is that a direct method can be effectively applied to solve the tridiagonal matrix resulted from numerical discretization of the 1D Helmholtz equation. The method has been tested to compute 2D Helmholtz solutions simulating electromagnetic scattering from an open large cavity and rectangular waveguide.

This paper presents a fast surface voxelization technique for the mapping of tessellated triangular surface meshes to uniform and structured grids that provide a basis for CFD simulations with the lattice Boltzmann method (LBM). The core algorithm is optimized for massively parallel execution on graphics processing units (GPUs) and is based on a unique dissection of the inner body shell. This unique definition necessitates a topology based neighbor search as a preprocessing step, but also enables parallel implementation. More specifically, normal vectors of adjacent triangular tessellations are used to construct half-angles that clearly separate the per-triangle regions. For each triangle, the grid nodes inside the axis-aligned bounding box (AABB) are tested for their distance to the triangle in question and for certain well-defined relative angles. The performance of the presented grid generation procedure is superior to the performance of the GPU-accelerated flow field computations per time step which allows efficient fluid-structure interaction simulations, without noticeable performance loss due to the dynamic grid update.

The coordinate transformation offers a remarkable way to design cloaks that can steer electromagnetic fields so as to prevent waves from penetrating into the cloaked region (denoted by Ω0, where the objects inside are invisible to observers outside). The ideal circular and elliptic cylindrical cloaked regions are blown up from a point and a line segment, respectively so the transformed material parameters and the corresponding coefficients of the resulted equations are highly singular at the cloaking boundary ∂Ω0. The electric field or magnetic field is not continuous across ∂Ω0. The imposition of appropriate cloaking boundary conditions (CBCs) to achieve perfect concealment is a crucial but challenging issue.

Based upon the principle that a well-behaved electromagnetic field in the original space must be well-behaved in the transformed space as well, we obtain CBCs that intrinsically relate to the essential “pole” conditions of a singular transformation. We also find that for the elliptic cylindrical cloak, the CBCs should be imposed differently for the cosine-elliptic and sine-elliptic components of the decomposed fields. With these at our disposal, we can rigorously show that the governing equation in Ω0 can be decoupled from the exterior region , and the total fields in the cloaked region vanish under mild conditions. We emphasize that our proposal of CBCs is different from any existing ones.

Using the exact circular (resp., elliptic) Dirichlet-to-Neumann (DtN) non-reflecting boundary conditions to reduce the unbounded domain to a bounded domain, we introduce an accurate and efficient Fourier-Legendre spectral-element method (FLSEM) (resp., Mathieu-Legendre spectral-element method (MLSEM)) to simulate the circular cylindrical cloak (resp., elliptic cylindrical cloak). We provide ample numerical results to demonstrate that the perfect concealment of waves can be achieved for the ideal circular/elliptic cylindrical cloaks under our proposed CBCs and accurate numerical solvers.