In this paper, we present an efficient rescaling scheme for computing the long-time dynamics of expanding interfaces. The idea is to design an adaptive time-space mapping such that in the new time scale, the interfaces evolves logarithmically fast at early growth stage and exponentially fast at later times. The new spatial scale guarantees the conservation of the area/volume enclosed by the interface. Compared with the original rescaling method in [J. Comput. Phys. 225(1) (2007) 554–567], this adaptive scheme dramatically improves the slow evolution at early times when the size of the interface is small. Our results show that the original three-week computation in [J. Comput. Phys. 225(1) (2007) 554–567] can be reproduced in about one day using the adaptive scheme. We then present the largest and most complicated Hele-Shaw simulation up to date.

A thermal lattice Boltzmann method (LBM) for two-phase fluid flows in nucleate pool boiling process is proposed. In the present method, a new function for heat transfer is introduced to the isothermal LBM for two-phase immiscible fluids with large density differences. The calculated temperature is substituted into the pressure tensor, which is used for the calculation of an order parameter representing two phases so that bubbles can be formed by nucleate boiling. By using this method, two-dimensional simulations of nucleate pool boiling by a heat source on a solid wall are carried out with the boundary condition for a constant heat flux. The flow characteristics and temperature distribution in the nucleate pool boiling process are obtained. It is seen that a bubble nucleation is formed at first and then the bubble grows and leaves the wall, finally going up with deformation by the buoyant effect. In addition, the effects of the gravity and the surface wettability on the bubble diameter at departure are numerically investigated. The calculated results are in qualitative agreement with other theoretical predictions with available experimental data.

In this paper, a novel implementation of immersed interface method combined with Stokes solver on a MAC staggered grid for solving the steady two-fluid Stokes equations with interfaces. The velocity components along the interface are introduced as two augmented variables and the resulting augmented equation is then solved by the GMRES method. The augmented variables and /or the forces are related to the jumps in pressure and the jumps in the derivatives of both pressure and velocity, and are interpolated using cubic splines and are then applied to the fluid through the jump conditions. The Stokes equations are discretized on a staggered Cartesian grid via a second order finite difference method and solved by the conjugate gradient Uzawa-type method. The numerical results show that the overall scheme is second order accurate. The major advantages of the present IIM-Stokes solver are the efficiency and flexibility in terms of types of fluid flow and different boundary conditions. The proposed method avoids solution of the pressure Poisson equation, and comparisons are made to show the advantages of time savings by the present method. The generalized two-phase Stokes solver with correction terms has also been applied to incompressible two-phase Navier-Stokes flow.

We consider the relativistic Euler equations governing spherically symmetric, perfect fluid flows on the outer domain of communication of Schwarzschild space-time, and we introduce a version of the finite volume method which is formulated from the geometric formulation (and thus takes the geometry into account at the discretization level) and is well-balanced, in the sense that it preserves steady solutions to the Euler equations on the curved geometry under consideration. In order to formulate our method, we first derive a closed formula describing all steady and spherically symmetric solutions to the Euler equations posed on Schwarzschild spacetime. Second, we describe a geometry-preserving, finite volume method which is based from the family of steady solutions to the Euler system. Our scheme is second-order accurate and, as required, preserves the family of steady solutions at the discrete level. Numerical experiments are presented which demonstrate the efficiency and robustness of the proposed method even for solutions containing shock waves and nonlinear interacting wave patterns. As an application, we investigate the late-time asymptotics of perturbed steady solutions and demonstrate its convergence for late time toward another steady solution, taking the overall effect of the perturbation into account.

Flexible discretization techniques for the approximative solution of coupled wave propagation problems are investigated. In particular, the advantages of using non-matching grids are presented, when one subregion has to be resolved by a substantially finer grid than the other subregion. We present the non-matching grid technique for the case of a mechanical-acoustic coupled as well as for acoustic-acoustic coupled systems. For the first case, the problem formulation remains essentially the same as for the matching situation, while for the acoustic-acoustic coupling, the formulation is enhanced with Lagrange multipliers within the framework of Mortar Finite Element Methods. The applications will clearly demonstrate the superiority of the Mortar Finite Element Method over the standard Finite Element Method both concerning the flexibility for the mesh generation as well as the computational time.

The implementation of a turbulent gas-kinetic scheme into a finite-volume RANS solver is put forward, with two turbulent quantities, kinetic energy and dissipation, supplied by an allied turbulence model. This paper shows a number of numerical simulations of flow cases including an interaction between a shock wave and a turbulent boundary layer, where the shock-turbulent boundary layer is captured in a much more convincing way than it normally is by conventional schemes based on the Navier-Stokes equations. In the gas-kinetic scheme, the modeling of turbulence is part of the numerical scheme, which adjusts as a function of the ratio of resolved to unresolved scales of motion. In so doing, the turbulent stress tensor is not constrained into a linear relation with the strain rate. Instead it is modeled on the basis of the analogy between particles and eddies, without any assumptions on the type of turbulence or flow class. Conventional schemes lack multiscale mechanisms: the ratio of unresolved to resolved scales – very much like a degree of rarefaction – is not taken into account even if it may grow to non-negligible values in flow regions such as shocklayers. It is precisely in these flow regions, that the turbulent gas-kinetic scheme seems to provide more accurate predictions than conventional schemes.

Newtonian and non-Newtonian dispersed phase droplet formation in non-Newtonian continuous phase in T-junction and cross junction microchannels are investigated by the immiscible lattice BGK model. The effects of the non-Newtonian fluid power-law exponent, viscosity and interfacial tension on the generation of the droplet are studied. The final droplet size, droplet generation frequency, and detachment point of the droplet change with the power-law exponent. The results reveal that it is necessary to take into account the non-Newtonian rheology instead of simple Newtonian fluid assumption in numerical simulations. The present analysis also demonstrates that the lattice Boltzmann method is of potential to investigate the non-Newtonian droplet generation in multiphase flow.

The modified ghost fluid method (MGFM), due to its reasonable treatment for ghost fluid state, has been shown to be robust and efficient when applied to compressible multi-medium flows. Other feasible definitions of the ghost fluid state, however, have yet to be systematically presented. By analyzing all possible wave structures and relations for a multi-medium Riemann problem, we derive all the conditions to define the ghost fluid state. Under these conditions, the solution in the real fluid region can be obtained exactly, regardless of the wave pattern in the ghost fluid region. According to the analysis herein, a practical ghost fluid method (PGFM) is proposed to simulate compressible multi-medium flows. In contrast with the MGFM where three degrees of freedomat the interface are required to define the ghost fluid state, only one degree of freedomis required in this treatment. However, when these methods proved correct in theory are used in computations for the multi-medium Riemann problem, numerical errors at the material interface may be inevitable. We show that these errors are mainly induced by the single-medium numerical scheme in essence, rather than the ghost fluid method itself. Equipped with some density-correction techniques, the PGFM is found to be able to suppress these unphysical solutions dramatically.

This paper is devoted to the discussion of numerical methods for solving two-dimensional time-fractional advection-diffusion equations. Two different three-point combined compact alternating direction implicit (CC-ADI) schemes are proposed and then, the original schemes for solving the two-dimensional problems are divided into two separate one-dimensional cases. Local truncation errors are analyzed and the unconditional stabilities of the obtained schemes are investigated by Fourier analysis method. Numerical experiments show the effectiveness and the spatial higher-order accuracy of the proposed methods.

We extend classic Sommerfeld and Silver-Müller radiation conditions for bounded scatterers to acoustic and electromagnetic fields propagating over three isotropic homogeneous layers in three dimensions. If X= (x1,x2,x3)ϵℝ3, with x3 denoting the direction orthogonal to the layers, standard conditions only hold for the outer layers in the region ∣x3∣ > ∣∣x∣γ, for γϵ(1/4,1/2) and x large. For ∣x3∣ < ∣∣x∣∣γ and inside the slab, asymptotic behavior depends on the presence of surface or guided modes given by the discrete spectrum of the associated operator.

The important task of evaluating the impact of random parameters on the output of stochastic ordinary differential equations (SODE) can be computationally very demanding, in particular for problems with a high-dimensional parameter space. In this work we consider this problem in some detail and demonstrate that by combining several techniques one can dramatically reduce the overall cost without impacting the predictive accuracy of the output of interests. We discuss how the combination of ANOVA expansions, different sparse grid techniques, and the total sensitivity index (TSI) as a pre-selective mechanism enables the modeling of problems with hundred of parameters. We demonstrate the accuracy and efficiency of this approach on a number of challenging test cases drawn from engineering and science.

We construct and analyze conservative local discontinuous Galerkin (LDG) methods for the Generalized Korteweg-de-Vries equation. LDG methods are designed by writing the equation as a system and performing separate approximations to the spatial derivatives. The main focus is on the development of conservative methods which can preserve discrete versions of the first two invariants of the continuous solution, and a posteriori error estimates for a fully discrete approximation that is based on the idea of dispersive reconstruction. Numerical experiments are provided to verify the theoretical estimates.

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.

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.

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.

In this paper, we investigate the effects of kinematic transports on the nematic liquid crystal system numerically and theoretically. The model we used is a “1+2” elastic continuum model simplified from the Ericksen-Leslie system. The numerical experiments are carried out by using a Legendre-Galerkin spectral method which can preserve the energy law in the discrete form. Based on this highly accurate numerical approach we find some interesting and important relationships between the kinematic transports and the characteristics of the flow. We make some analysis to explain these results. Several significant scaling properties are also verified by our simulations.

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.

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 general analysis of the hydrodynamic limit of multi-relaxation time lattice Boltzmann models is presented. We examine multi-relaxation time BGK collision operators that are constructed similarly to those for the MRT case, however, without explicitly moving into a moment space representation. The corresponding ‘moments’ are derived as left eigenvectors of said collision operator in velocity space. Consequently we can, in a representation independent of the chosen base velocity set, generate the conservation equations. We find a significant degree of freedom in the choice of the collision matrix and the associated basis which leaves the collision operator invariant. We explain why MRT implementations in the literature reproduce identical hydrodynamics despite being based on different orthogonalization relations. More importantly, however, we outline a minimal set of requirements on the moment base necessary to maintain the validity of the hydrodynamic equations. This is particularly useful in the context of position and time-dependent moments such as those used in the context of peculiar velocities and some implementations of fluctuations in a lattice-Boltzmann simulation.

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.