A cheapest stable nonconforming finite element method is presented for solving the incompressible flow in a square cavity without smoothing the corner singularities. The stable cheapest nonconforming finite element pair based on P1×P0 on rectangularmeshes [29] is employed with a minimal modification of the discontinuous Dirichlet data on the top boundary, where is the finite element space of piecewise constant pressures with the globally one-dimensional checker-board pattern subspace eliminated. The proposed Stokes elements have the least number of degrees of freedom compared to those of known stable Stokes elements. Three accuracy indications for our elements are analyzed and numerically verified. Also, various numerous computational results obtained by using our proposed element show excellent accuracy.

We perform a comparison in terms of accuracy and CPU time between second order BDF semi-Lagrangian and Lagrange-Galerkin schemes in combination with high order finite element method. The numerical results show that for polynomials of degree 2 semi-Lagrangian schemes are faster than Lagrange-Galerkin schemes for the same number of degrees of freedom, however, for the same level of accuracy both methods are about the same in terms of CPU time. For polynomials of degree larger than 2, Lagrange-Galerkin schemes behave better than semi-Lagrangian schemes in terms of both accuracy and CPU time; specially, for polynomials of degree 8 or larger. Also, we have performed tests on the parallelization of these schemes and the speedup obtained is quasi-optimal even with more than 100 processors.

A front tracking method combined with the real ghost fluid method (RGFM) is proposed for simulations of fluid interfaces in two-dimensional compressible flows. In this paper the Riemann problem is constructed along the normal direction of interface and the corresponding Riemann solutions are used to track fluid interfaces. The interface boundary conditions are defined by the RGFM, and the fluid interfaces are explicitly tracked by several connected marker points. The Riemann solutions are also used directly to update the flow states on both sides of the interface in the RGFM. In order to validate the accuracy and capacity of the new method, extensive numerical tests including the bubble advection, the Sod tube, the shock-bubble interaction, the Richtmyer-Meshkov instability and the gas-water interface, are simulated by using the Euler equations. The computational results are also compared with earlier computational studies and it shows good agreements including the compressible gas-water system with large density differences.

How to effectively simulate the interaction between fluid and solid elements of different sizes remains to be challenging. The discrete element method (DEM) has been used to deal with the interactions between solid elements of various shapes and sizes, while the material point method (MPM) has been developed to handle the multiphase (solid-liquid-gas) interactions involving failure evolution. A combined MPM-DEM procedure is proposed to take advantage of both methods so that the interaction between solid elements and fluid particles in a container could be better simulated. In the proposed procedure, large solid elements are discretized by the DEM, while the fluid motion is computed using the MPM. The contact forces between solid elements and rigid walls are calculated using the DEM. The interaction between solid elements and fluid particles are calculated via an interfacial scheme within the MPM framework. With a focus on the boundary condition effect, the proposed procedure is illustrated by representative examples, which demonstrates its potential for a certain type of engineering problems.

This paper studies the magneto-heat coupling model which describes iron loss of conductors and energy exchange between magnetic field and Ohmic heat. The temperature influences Maxwell's equations through the variation of electric conductivity, while electric eddy current density provides the heat equation with Ohmic heat source. It is in this way that Maxwell's equations and the heat equation are coupled together. The system also incorporates the heat exchange between conductors and cooling oil which is poured into and out of the transformer. We propose a weak formulation for the coupling model and establish the well-posedness of the problem. The model is more realistic than the traditional eddy current model in numerical simulations for large power transformers. The theoretical analysis of this paper paves a way for us to design efficient numerical computation of the transformer in the future.

This study aims to develop a numerical scheme in collocated Cartesian grids to solve the level set equation together with the incompressible two-phase flow equations. A seventh-order accurate upwinding combined compact difference (UCCD7) scheme has been developed for the approximation of the first-order spatial derivative terms shown in the level set equation. Developed scheme has a higher accuracy with a three-point grid stencil to minimize phase error. To preserve the mass of each phase all the time, the temporal derivative term in the level set equation is approximated by the sixth-order accurate symplectic Runge-Kutta (SRK6) scheme. All the simulated results for the dam-break, Rayleigh-Taylor instability, bubble rising, two-bubble merging, and milkcrown problems in two and three dimensions agree well with the available numerical or experimental results.

Arterial diseases such as aneurysm and stenosis may result from the mechanical and/or morphological change of an arterial wall structure and correspondingly altered hemodynamics. The development of a 3D computational model of blood flow can be useful to study the hemodynamics in major blood vessels and may provide an insight into the noninvasive technique to detect arterial diseases in early stage. In this paper, we present a three-dimensional model of blood flow in the aorta, which is based on the immersed boundary method to describe the interaction of blood flow with the aortic wall. Our simulation results show that the hysteresis loop is evident in the pressure-diameter relationship of the normal aorta when the arterial wall is considered to be viscoelastic. In addition, it is shown that flow patterns and pressure distributions are altered in response to the change of aortic morphology.

Scale similarity is found in many natural phenomena in the universe, from fluid dynamics to astrophysics. In large eddy simulations of turbulent flows, some sub-grid scale (SGS) models are based on scale similarity. The earliest scale similarity SGS model was developed by Bardina et al., which produced SGS stresses with good correlation to the true stresses. In the present study, we perform a mathematical analysis of scale similarity. The analysis has revealed that the ratio of the resolved stress to the SGS stress is γ2, where γ is the ratio of the second filter width to the first filter width, under the assumption of small filter width. The implications of this analysis are discussed in the context of large eddy simulation.

Incorrect propagation speed of discontinuities may occur by straightforward application of standard dissipative schemes for problems that contain stiff source terms for underresolved grids even for time steps within the CFL condition. By examining the dissipative discretized counterpart of the Euler equations for a detonation problem that consists of a single reaction, detailed analysis on the spurious wave pattern is presented employing the fractional step method, which utilizes the Strang splitting. With the help of physical arguments, a threshold values method (TVM), which can be extended to more complicated stiff problems, is developed to eliminate the wrong shock speed phenomena. Several single reaction detonations as well as multispecies and multi-reaction detonation test cases with strong stiffness are examined to illustrate the performance of the TVM approach.

In two-phase flow simulations, a difficult issue is usually the treatment of surface tension effects. These cause a pressure jump that is proportional to the curvature of the interface separating the two fluids. Since the evaluation of the curvature incorporates second derivatives, it is prone to numerical instabilities. Within this work, the interface is described by a level-set method based on a discontinuous Galerkin discretization. In order to stabilize the evaluation of the curvature, a patch-recovery operation is employed. There are numerous ways in which this filtering operation can be applied in the whole process of curvature computation. Therefore, an extensive numerical study is performed to identify optimal settings for the patch-recovery operations with respect to computational cost and accuracy.

This paper analyzes an abstract two-level algorithm for hybridizable discontinuous Galerkin (HDG) methods in a unified fashion. We use an extended version of the Xu-Zikatanov (X-Z) identity to derive a sharp estimate of the convergence rate of the algorithm, and show that the theoretical results also are applied to weak Galerkin (WG) methods. The main features of our analysis are twofold: one is that we only need the minimal regularity of the model problem; the other is that we do not require the triangulations to be quasi-uniform. Numerical experiments are provided to confirm the theoretical results.

We study the inverse problem of recovering the scatterer shape from the far-field pattern(FFP) in the presence of noise. Furthermore, only a discrete partial aperture is usually known. This problem is ill-posed and is frequently addressed using regularization. Instead, we propose to use a direct approach denoising the FFP using a filtering technique. The effectiveness of the technique is studied on a scatterer with the shape of the ellipse with a tower. The forward scattering problem is solved using the finite element method (FEM). The numerical FFP is additionally corrupted by Gaussian noise. The shape parameters are found based on a least-square error estimator. If ũ∞ is a perturbation of the FFP then we attempt to find Γ, the scatterer shape, which minimizes ∣∣ũ∞ − ũ∞∣∣ using the conjugate gradient method for the denoised FFP

In this article, we discuss the stability of soft quasicrystalline phases in a coupled-mode Swift-Hohenberg model for three-component systems, where the characteristic length scales are governed by the positive-definite gradient terms. Classic two-mode approximation method and direct numerical minimization are applied to the model. In the latter approach, we apply the projection method to deal with the potentially quasiperiodic ground states. A variable cell method of optimizing the shape and size of higher-dimensional periodic cell is developed to minimize the free energy with respect to the order parameters. Based on the developed numerical methods, we rediscover decagonal and dodecagonal quasicrystalline phases, and find diverse periodic phases and complex modulated phases. Furthermore, phase diagrams are obtained in various phase spaces by comparing the free energies of different candidate structures. It does show not only the important roles of system parameters, but also the effect of optimizing computational domain. In particular, the optimization of computational cell allows us to capture the ground states and phase behavior with higher fidelity. We also make some discussions on our results and show the potential of applying our numerical methods to a larger class of mean-field free energy functionals.

Algorithms for adaptive mesh refinement using a residual error estimator are proposed for fluid flow problems in a finite volume framework. The residual error estimator, referred to as the ℜ-parameter is used to derive refinement and coarsening criteria for the adaptive algorithms. An adaptive strategy based on the ℜ-parameter is proposed for continuous flows, while a hybrid adaptive algorithm employing a combination of error indicators and the ℜ-parameter is developed for discontinuous flows. Numerical experiments for inviscid and viscous flows on different grid topologies demonstrate the effectiveness of the proposed algorithms on arbitrary polygonal grids.

Free energy lattice Boltzmann methods are well suited for the simulation of two phase flow problems. The model for the interface is based on well understood physical grounds. In most cases a numerical interface is used instead of the physical one because of lattice resolution limitations. In this paper we present a framework where we can both follow the droplet behavior in a coarse scale and solve the interface in a fine scale simultaneously. We apply the method for the simulation of a droplet using an interface to diameter size ratio of 1 to 280. In a second simulation, a small droplet coalesces with a 42 times larger droplet producing on it only a small capillary wave that propagates and dissipates.

Acupuncture has been in clinical practice in China for thousands of years and its analgesia effect is worldwide accepted. However, the mechanism of acupuncture effect is not well understood. The study focus on signaling pathways induced by acupuncture, analyzes the cooperative action of the acupoints’ structure and the associated chemical mediators during acupuncture, establishes a mathematical model clarifying the roadmap of electroneurographic signal startup and transmission mechanism induced by acupuncture, quantitatively analyzing the response in acupoints to acupuncture. These work contribute to reveal the activation and transmission mechanism of neural signals induced by acupuncture from systems biology perspective, lay the foundation for the integration of acupuncture theory and modern science and further guide the clinical treatment and experimental research of acupuncture.

The importance of eliminating errors in grid-metric evaluation for high-order difference schemes has been widely recognized in recent years, and it is known from the proof by Vinokur and Yee (NASA TM 209598, 2000) that when conservative derivations of grid metric are used by Thomas, Lombard and Neier (AIAA J., 1978, 17(10) and J. Spacecraft and rocket, 1990, 27(2)), errors caused by metric evaluation could be eliminated by linear schemes when flux splitting is not considered. According to the above achievement, central schemes without the use of flux splitting could fulfill the requirement of error elimination. Difficulties will arise for upwind schemes to attain the objective when the splitting is considered. In this study, further investigations are made on three aspects: Firstly, an idea of central scheme decomposition is introduced, and the procedure to derive the central scheme is proposed to evaluate grid metrics only. Secondly, the analysis has been made on the requirement of flux splitting to acquire free-stream preservation, and a Lax-Friedrichs-type splitting scheme is proposed as an example. Discussions about current study with that by Nonomura et al. (Computers and Fluids, 2015, 107) have been made. Thirdly, for half-node- or mixed-type schemes, interpolations should be used to derive variables at half nodes. The requirement to achieve metric identity on this situation is analyzed and an idea of directionally consistent interpolation is proposed, which is manifested to be indispensable to avoid violations of metric identity and to eliminate metric-caused errors thereafter. Two numerical problems are tested, i.e., the free-stream and vortex preservation on wavy, largely randomized and triangular-like grids. Numerical results validate aforementioned theoretical outcomes.

Moving mesh methods provide an efficient way of solving partial differential equations for which large, localised variations in the solution necessitate locally dense spatial meshes. In one-dimension, meshes are typically specified using the arclength mesh density function. This choice is well-justified for piecewise polynomial interpolants, but it is only justified for spectral methods when model solutions include localised steep gradients. In this paper, one-dimensional mesh density functions are presented which are based on a spatially localised measure of the bandwidth of the approximated model solution. In considering bandwidth, these mesh density functions are well-justified for spectral methods, but are not strictly tied to the error properties of any particular spatial interpolant, and are hence widely applicable. The bandwidth mesh density functions are illustrated in two ways. First, by applying them to Chebyshev polynomial approximation of two test functions, and second, through use in periodic spectral and finite-difference moving mesh methods applied to a number of model problems in acoustics. These problems include a heterogeneous advection equation, the viscous Burgers’ equation, and the Korteweg-de Vries equation. Simulation results demonstrate solution convergence rates that are up to an order of magnitude faster using the bandwidth mesh density functions than uniform meshes, and around three times faster than those using the arclength mesh density function.

We study the exciton diffusion in organic semiconductors from a macroscopic viewpoint. In a unified way, we conduct the equivalence analysis between Monte-Carlo method and diffusion equation model for photoluminescence quenching and photocurrent spectrum measurements, in both the presence and the absence of Förster energy transfer effect. Connections of these two models to Stern-Volmer method and exciton-exciton annihilation method are also specified for the photoluminescence quenching measurement.

The fine description of complex fluids can be carried out by describing the evolution of each individual constituent (e.g. each particle, each macromolecule, etc.). This procedure, despite its conceptual simplicity, involves many numerical issues, the most challenging one being that related to the computing time required to update the system configuration by describing all the interactions between the different individuals. Coarse grained approaches allow alleviating the just referred issue: the system is described by a distribution function providing the fraction of entities that at certain time and position have a particular conformation. Thus, mesoscale models involve many different coordinates, standard space and time, and different conformational coordinates whose number and nature depend on the particular system considered. Balance equation describing the evolution of such distribution function consists of an advection-diffusion partial differential equation defined in a high dimensional space. Standard mesh-based discretization techniques fail at solving high-dimensional models because of the curse of dimensionality. Recently the authors proposed an alternative route based on the use of separated representations. However, until now these approaches were unable to address the case of advection dominated models due to stabilization issues. In this paper this issue is revisited and efficient procedures for stabilizing the advection operators involved in the Boltzmann and Fokker-Planck equation within the PGD framework are proposed.