The numerical simulation of non conservative system is a difficult challenge for two reasons at least. The first one is that it is not possible to derive jump relations directly from conservation principles, so that in general, if the model description is non ambiguous for smooth solutions, this is no longer the case for discontinuous solutions. From the numerical view point, this leads to the following situation: if a scheme is stable, its limit for mesh convergence will depend on its dissipative structure. This is well known since at least [1]. In this paper we are interested in the “dual” problem: given a system in non conservative form and consistent jump relations, how can we construct a numerical scheme that will, for mesh convergence, provide limit solutions that are the exact solution of the problem. In order to investigate this problem, we consider a multiphase flow model for which jump relations are known. Our scheme is an hybridation of Glimm scheme and Roe scheme.

An accurate cartesian method is devised to simulate incompressible viscous flows past an arbitrary moving body. The Navier-Stokes equations are spatially discretized onto a fixed Cartesian mesh. The body is taken into account via the ghost-cell method and the so-called penalty method, resulting in second-order accuracy in velocity. The accuracy and the efficiency of the solver are tested through two-dimensional reference simulations. To show the versatility of this scheme we simulate a three-dimensional self propelled jellyfish prototype.

In this paper we develop a class of Implicit-Explicit Runge-Kutta schemes for solving the multi-scale semiconductor Boltzmann equation. The relevant scale which characterizes this kind of problems is the diffusive scaling. This means that, in the limit of zero mean free path, the system is governed by a drift-diffusion equation. Our aim is to develop a method which accurately works for the different regimes encountered in general semiconductor simulations: the kinetic, the intermediate and the diffusive one. Moreover, we want to overcome the restrictive time step conditions of standard time integration techniques when applied to the solution of this kind of phenomena without any deterioration in the accuracy. As a result, we obtain high order time and space discretization schemes which do not suffer from the usual parabolic stiffness in the diffusive limit. We show different numerical results which permit to appreciate the performances of the proposed schemes.

The carbuncle phenomenon has been regarded as a spurious solution produced by most of contact-preserving methods. The hybrid method of combining high resolution flux with more dissipative solver is an attractive attempt to cure this kind of non-physical phenomenon. In this paper, a matrix-based stability analysis for 2-D Euler equations is performed to explore the cause of instability of numerical schemes. By combining the Roe with HLL flux in different directions and different flux components, we give an interesting explanation to the linear numerical instability. Based on such analysis, some hybrid schemes are compared to illustrate different mechanisms in controlling shock instability. Numerical experiments are presented to verify our analysis results. The conclusion is that the scheme of restricting directly instability source is more stable than other hybrid schemes.

We present a high precision frequency determination method for digitized NMR FID signals. The method employs high precision numerical integration rather than simple summation as in many other techniques. With no independent knowledge of the other parameters of a NMR FID signal (phase ф, amplitude A, and transverse relaxation time T2) this method can determine the signal frequency f0 with a precision of if the observation time T ≫ T2. The method is especially convenient when the detailed shape of the observed FT NMR spectrum is not well defined. When T2 is +∞ and the signal becomes pure sinusoidal, the precision of the method is which is one order more precise than the ±1 count error induced precision of a typical frequency counter. Analysis of this method shows that the integration reduces the noise by bandwidth narrowing as in a lock-in amplifier, and no extra signal filters are needed. For a pure sinusoidal signal we find from numerical simulations that the noise-induced error in this method reaches the Cramer-Rao Lower Band (CRLB) on frequency determination. For the damped sinusoidal case of most interest, the noise-induced error is found to be within a factor of 2 of CRLB when the measurement time T is 2 or 3 times larger than T2. We discuss possible improvements for the precision of this method.

Quantum Monte Carlo data are often afflicted with distributions that resemble lognormal probability distributions and consequently their statistical analysis cannot be based on simple Gaussian assumptions. To this extent a method is introduced to estimate these distributions and thus give better estimates to errors associated with them. This method entails reconstructing the probability distribution of a set of data, with given mean and variance, that has been assumed to be lognormal prior to undergoing a blocking or renormalization transformation. In doing so, we perform a numerical evaluation of the renormalized sum of lognormal random variables. This technique is applied to a simple quantum model utilizing the single-thread Monte Carlo algorithm to estimate the ground state energy or dominant eigenvalue of a Hamiltonian matrix.

We develop the dimension-reduced hyperbolic moment method for the Boltzmann equation, to improve solution efficiency using a numerical regularized moment method for problems with low-dimensional macroscopic variables and high-dimensional microscopic variables. In the present work, we deduce the globally hyperbolic moment equations for the dimension-reduced Boltzmann equation based on the Hermite expansion and a globally hyperbolic regularization. The numbers of Maxwell boundary condition required for well-posedness are studied. The numerical scheme is then developed and an improved projection algorithm between two different Hermite expansion spaces is developed. By solving several benchmark problems, we validate the method developed and demonstrate the significant efficiency improvement by dimension-reduction.

Correlation functions in the (n) models below the critical temperature are considered. Based on Monte Carlo (MC) data, we confirm the fact stated earlier by Engels and Vogt, that the transverse two-plane correlation function of the (4) model for lattice sizes about L = 120 and small external fields h is very well described by a Gaussian approximation. However, we show that fits of not lower quality are provided by certain non-Gaussian approximation. We have also tested larger lattice sizes, up to L = 512. The Fourier-transformed transverse and longitudinal two-point correlation functions have Goldstone mode singularities in the thermodynamic limit at k → 0 and h = +0, i.e., G⊥ (k) ≃ ak–λ⊥ and G‖(k)≃bk–λ‖, respectively. Here a and b are the amplitudes, k = |k| is the magnitude of the wave vector k. The exponents λᚆ, λ‖ and the ratio bM2/a2, where M is the spontaneous magnetization, are universal according to the GFD (grouping of Feynman diagrams) approach. Here we find that the universality follows also from the standard (Gaussian) theory, yielding bM2/a2=(n−1)/16. Our MC estimates of this ratio are 0.06±0.01 for n=2, 0.17±0.01 for n = 4 and 0.498±0.010 for n = 10. According to these and our earlier MC results, the asymptotic behavior and Goldstone mode singularities are not exactly described by the standard theory. This is expected from the GFD theory. We have found appropriate analytic approximations for G⊥(k) and G‖(k), well fitting the simulation data for small k. We have used them to test the Patashinski-Pokrovski relation and have found that it holds approximately.

This paper concerns the electromagnetic scattering by arbitrary shaped three dimensional imperfectly conducting objects modeled with non-constant Leontovitch impedance boundary condition. It has two objectives. Firstly, the intrinsically well-conditioned integral equation (noted GCSIE) proposed in [30] is described focusing on its discretization. Secondly, we highlight the potential of this method by comparison with two other methods, the first being a two currents formulation in which the impedance condition is implicitly imposed and whose the convergence is quasi-optimal for Lipschitz polyhedron, the second being a CFIE-like formulation [14]. In particular, we prove that the new approach is less costly in term of CPU time and gives a more accurate solution than that obtained from the CFIE formulation. Finally, as expected, It is demonstrated that no preconditioner is needed for this formulation.

A weak Galerkin (WG) method is introduced and numerically tested for the Helmholtz equation. This method is flexible by using discontinuous piecewise polynomials and retains the mass conservation property. At the same time, the WG finite element formulation is symmetric and parameter free. Several test scenarios are designed for a numerical investigation on the accuracy, convergence, and robustness of the WG method in both inhomogeneous and homogeneous media over convex and non-convex domains. Challenging problems with high wave numbers are also examined. Our numerical experiments indicate that the weak Galerkin is a finite element technique that is easy to implement, and provides very accurate and robust numerical solutions for the Helmholtz problem with high wave numbers.

We present a multi-phase image segmentation method based on the histogram of the Gabor feature space, which consists of a set of Gabor-filter responses with various orientations, scales and frequencies. Our model replaces the error function term in the original fuzzy region competition model with squared 2-Wasserstein distance function, which is a metric to measure the distance of two histograms. The energy functional is minimized by alternative minimization method and the existence of closed-form solutions is guaranteed when the exponent of the fuzzy membership term being 1 or 2. We test our model on both simple synthetic texture images and complex natural images with two or more phases. Experimental results are shown and compared to other recent results.