NUMERICAL EXPERIMENTS FOR NONLINEAR BURGER’S PROBLEM

The article contains the results of computational experiments for the non-homogeneous Burger’s problem and finding its solution by using the non-classical variational-Cole-Hopf transformation approach. On using exact linearization via Cole-Hopf transformation, as well as the application of the non-classical variational approach, then the non-homogeneous Burger’s problem has been solved. The solution which is obtained by this approach is in a compact form so that the original nonlinear solution is easy to be approximated. The accuracy of this method is dependent on the types of selected basis and its number.

Efficient Localized Nonlinear Solution Strategies for Unconventional-Reservoir Simulation with Complex Fractures

It is very challenging to simulate unconventional reservoirs efficiently and accurately. Transient flow can last for a long time and sharp solution (pressure, saturation, compositions) gradients are induced because of the severe permeability contrast between fracture and matrix. Although high-resolution models for well and fracture are required to achieve adequate resolution, they are computationally too demanding for practical field models with many stages of hydraulic fracture. The paper aims to innovate localization strategies that take advantage of locality on timestep and Newton iteration levels. The strategies readily accommodate to complicated flow mechanisms and multiscale fracture networks in unconventional reservoirs. Large simulation speed-up can be obtained if performing localized computations only for the solution regions that will change. We develop an a-priori method to exploit the locality, based on the diffusive character of the Newton updates of pressure. The method makes adequate estimate of the active computational gridblock for the next iterate. The active gridblock set marks the ones need to be solved, and then the solution to local linear system is accordingly computed. Fully Implicit Scheme is used for time discretization. We study several challenging multi-phase and compositional model cases with explicit fractures. The test results demonstrate that significant solution locality of variables exist on timestep and iteration levels. A nonlinear solution update usually has sparsity, and the nonlinear convergence is restricted by a limited fraction of the simulation model. Through aggressive localization, the proposed methods can prevent overly conservative estimate, and thus achieve significant computational speedup. In comparison to a standard Newton method, the novel solver techniques achieve greatly improved solving efficiency. Furthermore, the Newton convergence exhibits no degradation, and there is no impact on the solution accuracy. Previous works in the literature largely relate to the meshing aspect that accommodates to horizontal wells and hydraulic fractures. We instead develop new nonlinear strategies to perform localization. In particular, the adaptive DD method produces proper domain partitions according to the fluid flow and nonlinear updates. This results in an effective strategy that maintains solution accuracy and convergence behavior.

Computation of Nonlinear Thermoelectric Effects with Adaptive Methods

An implementation for a fully automatic adaptive finite element method (AFEM) for computation of nonlinear thermoelectric problems in three dimensions is presented. Adaptivity of the nonlinear solvers is based on the well-established hp-adaptivity where the mesh refinement and the polynomial order of elements are methodically controlled to reduce the discretization errors of the coupled field variables temperature and electric potential. A single mesh is used for both fields and the nonlinear coupling of temperature and electric potential is accounted in the computation of a posteriori error estimate where the residuals are computed element-wise. Mesh refinements are implemented for tetrahedral mesh such that conformity of elements with neighboring elements is preserved. Multiple nonlinear solution steps are assessed including variations of the fixed-point method with Anderson acceleration algorithms. The Barzilai-Borwein algorithm to optimize the nonlinear solution steps are also assessed. Promising results have been observed where all the nonlinear methods show the same accuracy with the tendency of approaching convergence with more elements refining. Anderson acceleration is the most efficient among the nonlinear solvers studied where its total computing time is less than half of the more conventional fixed-point iteration.

Pointwise Convergence of the Schrödinger Flow

In this paper we address the question of the pointwise almost everywhere limit of nonlinear Schrödinger flows to the initial data, in both the continuous and the periodic settings. Then we show how, in some cases, certain smoothing effects for the non-homogeneous part of the solution can be used to upgrade to a uniform convergence to zero of this part, and we discuss the sharpness of the results obtained. We also use randomization techniques to prove that with much less regularity of the initial data, both in continuous and the periodic settings, almost surely one obtains uniform convergence of the nonlinear solution to the initial data, hence showing how more generic results can be obtained.

Effect of Noise and Nonlinearities on Thermoacoustics of Can-Annular Combustors

Can-annular combustors consist of N distinct cans setup symmetrically around the axis of the gas turbine. Each can is connected to the turbine inlet by means of a transition duct. At the turbine inlet, a small gap between the neighboring transition ducts allows acoustic communication between the cans. Thermoacoustic pulsations in the cans are driven by the respective flames, but also the communication between neighboring cans through the gap plays a significant role. In this study, we focus on the effect of the background noise intensity and of the nonlinear flame saturation. We predict how usually clusters of thermoacoustic modes are unstable in the linear regime and compete with each other in the nonlinear regime, each cluster consisting of axial, azimuthal and push-pull modes. Since linear theory cannot predict the nonlinear solution, stochastic simulations are run to study the nonlinear solution in a probabilistic sense. One outcome of these simulations is the various pulsation patterns, which are in principle different from one can to the next. We recover how not only a stronger flame response in one can gives rise to the phenomenon of mode localization, but also how the nonlinearity of the flame saturation and the competition between modes have an effect on the nonlinear mode shape. We finally predict the coherence and phase between cans on the linearized system subject to noise, and compare the predictions with engine measurements, in terms of spectra of amplitude in each can and coherence and phase, observing a good match.