Hybrid Finite Differences - Lattice-Boltzmann Method for Multiphase Flows: Application to a Jet Simulation

We present a discontinuous Galerkin (DG)–finite difference (FD) hybrid scheme that allows high-order shock capturing with the DG method for general relativistic magnetohydrodynamics. The hybrid method is conceptually quite simple. An unlimited DG candidate solution is computed for the next time step. If the candidate solution is inadmissible, the time step is retaken using robust FD methods. Because of its a posteriori nature, the hybrid scheme inherits the best properties of both methods. It is high-order with exponential convergence in smooth regions, while robustly handling discontinuities. We give a detailed description of how we transfer the solution between the DG and FD solvers, and the troubled-cell indicators necessary to robustly handle slow-moving discontinuities and simulate magnetized neutron stars. We demonstrate the efficacy of the proposed method using a suite of standard and very challenging 1D, 2D, and 3D relativistic magnetohydrodynamics test problems. The hybrid scheme is designed from the ground up to efficiently simulate astrophysical problems such as the inspiral, coalescence, and merger of two neutron stars.

One of the major problems in the modeling and inversion of marine controlled-source electromagnetic (CSEM) data is related to the need for accurate representation of very complex geoelectrical models typical for marine environment. At the same time, the corresponding forward-modeling algorithms should be powerful and fast enough to be suitable for repeated use in hundreds of iterations of the inversion and for multiple transmitter/receiver positions. To this end, we have developed a novel 3D modeling and inversion approach, which combines the advantages of the finite-difference (FD) and integral-equation (IE) methods. In the framework of this approach, we have solved Maxwell’s equations for anomalous electric fields using the FD approximation on a staggered grid. Once the unknown electric fields in the computation domain of the FD method are computed, the electric and magnetic fields at the receivers are calculated using the IE method with the corresponding Green’s tensor for the background conductivity model. This approach makes it possible to compute the fields at the receivers accurately without the need of very fine FD discretization in the vicinity of the receivers and sources and without the need for numerical differentiation and interpolation. We have also developed an algorithm for 3D inversion based on the hybrid FD-IE method. In the case of the marine CSEM problem with multiple transmitters and receivers, the forward modeling and the Fréchet derivative calculations are very time consuming and require using large memory to store the intermediate results. To overcome those problems, we have applied the moving sensitivity domain approach to our inversion. A case study for the 3D inversion of towed streamer EM data collected by PGS over the Troll field in the North Sea demonstrated the effectiveness of the developed hybrid method.

This paper presents a mixed finite volume and finite difference solver with results showing the solitary wave interactions and shoaling process by solving a set of conservative forms of Boussinesq equations. A second order accurate finite volume scheme is applied to the conservative terms of the governing equations while up to the second order finite difference formulations are used to discretize the dispersive source terms with higher order derivatives. The limiters and surface gradient method are implemented in the model to remove the unwanted spurious oscillations and preserve the still water condition without introducing errors at the interfaces. The performance of the present numerical solver is tested with results of head on collisions and shoaling of solitary waves compared against those from finite element models that were developed based on fully nonlinear weakly dispersive and weakly nonlinear weakly dispersive forms of the Boussinesq equations as well as analytical solutions and experimental observations.

A Galilean invariant model of the lattice Boltzmann method for multiphase fluid flows using free-energy approach

Dispersion modeling in pore networks: A comparison of common pore-scale models and alternative approaches

This paper aims to propose a hybrid approach of two powerful methods, namely the differential transform and finite difference methods, to obtain the solution of the coupled Whitham-Broer-Kaup-Like equations which arises in shallow-water wave theory. The capability of the method to such problems is verified by taking different parameters and initial conditions. The numerical simulations are depicted in 2D and 3D graphs. It is shown that the used approach returns accurate solutions for this type of problems in comparison with the analytic ones.

A Boussinesq-type wave model is developed in this paper to simulate the interaction of coastal waves with bottom-mounted porous structures. The governing equations are rewritten in the conservative form to facilitate the use of hybrid finite volume (FV) and finite difference (FD) method. Higher-order slope terms are also inserted into the equations to account for rapidly varying bathymetry. The convective flux is approximated using the FV method, while the remaining terms are discretized using the FD method in a uniform rectangle grid system. The time integration is implemented using the third order Runge–Kutta method with an adaptive time step. A single GPU parallel computation is also implemented to save computation costs. The numerical model is validated against a series of experimental datasets, including data acquired in a new laboratory experiment. The predictions are in overall agreement with the measurements, proving that the model is capable of handling wave interaction with porous structures in the coastal region for a wide range of scenarios.

Modeling positive surge propagation in open channels using the Serre-Green-Naghdi equations

Finite volume scheme for the solution of 2D extended Boussinesq equations in the surf zone

Benchmarking of different numerical models simulating groundwater flow and contaminant mass transport is the aim of the present study, in order to determine criteria for the selection of numerical model(s) that could be better tailored to the needs of a specific region. This analysis aims at evaluating the performance of a finite difference-based numerical model (MODFLOW-ΜΤ3DMS), a finite element-based numerical model (FEFLOW), and a hybrid finite element-finite difference coupling numerical model (Princeton Transport Code-PTC), all developed to simulate groundwater flow and nitrate mass transport in an alluvial aquifer. The evaluation of the models’ performance is assessed based on statistical measures and graphical performance analysis of the model point predictions to the observed values. The outcome of the analysis showed that among the three groundwater simulation models, FEFLOW algorithm exhibited the best performance in simulating both groundwater level and nitrate mass distribution. All simulation algorithms were found to offer different advantages, so in principle the selection of the appropriate model(s) should be done in accordance with the problem’s characteristics and/ or in a complementary way in order to achieve accurate representation of the aquifer system and thus optimal groundwater resources management. Even though the selection of the most suitable groundwater simulation algorithm is case-oriented, however, fractional gross error (FGE) was proven to be a reliable indicator that could be used by modelers to select the most suitable groundwater algorithm based on the available groundwater data.

One of the fundamental problems in simulating the motion of sharp interfaces between immiscible fluids is a description of the transition that occurs when the interfaces merge and reconnect. It is well known that classical methods involving sharp interfaces fail to describe this type of phenomena. Following some previous work in this area, we suggest a physically motivated regularization of the Euler equations which allows topological transitions to occur smoothly. In this model, the sharp interface is replaced by a narrow transition layer across which the fluids may mix. The model describes a flow of a binary mixture, and the internal structure of the interface is determined by both diffusion and motion. An advantage of our regularization is that it automatically yields a continuous description of surface tension, which can play an important role in topological transitions. An additional scalar field is introduced to describe the concentration of one of the fluid components and the resulting system of equations couples the Euler (or Navier–Stokes) and the Cahn–Hilliard equations. The model takes into account weak non–locality (dispersion) associated with an internal length scale and localized dissipation due to mixing. The non–locality introduces a dimensional surface energy; dissipation is added to handle the loss of regularity of solutions to the sharp interface equations and to provide a mechanism for topological changes. In particular, we study a non–trivial limit when both components are incompressible, the pressure is kinematic but the velocity field is non–solenoidal (quasi–incompressibility). To demonstrate the effects of quasi–incompressibility, we analyse the linear stage of spinodal decomposition in one dimension. We show that when the densities of the fluids are not perfectly matched, the evolution of the concentration field causes fluid motion even if the fluids are inviscid. In the limit of infinitely thin and well–separated interfacial layers, an appropriately scaled quasi–incompressible Euler–Cahn–Hilliard system converges to the classical sharp interface model. In order to investigate the behaviour of the model outside the range of parameters where the sharp interface approximation is sufficient, we consider a simple example of a change of topology and show that the model permits the transition to occur without an associated singularity.

This is the first year in the proposed three-year effort to develop high-resolution numerical methods for multi-phase flow in hierarchical porous media. The issues being addressed in this research are: Computational efficiency: Field-scale simulation of enhanced oil recovery, whether for energy production or aquifer remediation, is typically highly under-resolved. This is because rock transport properties vary on many scales, and because current numerical methods have low resolution. Effective media properties: Since porous media are formed through complex geologic processes, they involve significant uncertainty and scale-dependence. Given this uncertainty, knowledge of ensemble averages of flow in porous media can be preferable to knowledge of flow in specific realizations of the reservoir. However, current models of effective properties do not represent the observed behavior very well. Relative permeability models present a good example of this problem. In practice, these models seldom provide realistic representations of hysteresis, interfacial tension effects or three-phase flow; there are no models that represent well all three effects simultaneously. Wave propagation: It is common in the petroleum industry to assume that the models have the same well-posedness properties as the physical system. An example of this fallacy is given by the three-phase relative permeability models; they were widely assumed by the petroleum community to produce hyperbolic systems for the Buckley-Leverett equations, but later the mathematics community proved that these models inherently produce local elliptic regions. Since numerical methods must use the models for computations, oscillations that develop could erroneously be attributed to numerical error rather than modeling difficulties. During this year, we have made significant progress on several tasks aimed at addressing these issues.

Three dimensional simulator of the single electron transistor based on ISIS structure

In the present study, the hybrid differential transformation and finite difference method is applied to analyze the dynamic behavior of the nonlinear micro-electromechanical circular plate actuated by combined DC/AC loading schemes. The analysis takes account of the axial residual stress and hydrostatic pressure acting on micro circular plate upper surface. The dynamic response of the plate as a function of the magnitude of the AC driving voltage is explored. Moreover, the effect of the initial gap height on the pull-in voltage of the plate is systematically explored.

Buoyancy-driven natural-convection heat transfer in enclosures has been the subject of considerable research with applications to electronic packaging, solar collectors, and shipping containers for spent nuclear fuel. A numerical study has been carried out to predict combined natural-convection and radiation heat transfer in the annular region between concentric tubes. The inner tube was volumetrically heated. Both tubes were of finite conductance. The surfaces of the annular region were diffuse and gray. The gas in the annulus was assumed to be nonparticipating. A newly developed hybrid finite element finite difference method was used for the study